%Samp4_1
clf;
b=[1 2 1]; a=[1 0.5 0.25];   %傳遞函數分子和分母多項式系數
m=0:length(b)-1;     n=0:length(a)-1;  %序列序號
K=512;   %采用512個點繪制幅頻特性
k=1:K;w=pi*k/K;     %轉換為歸一化頻率，Nyquist頻率對應于pi。
num=b*exp(-j*m'*w);    %分子多項式值序列
den=a*exp(-j*n'*w);   %分母多項式值序列
H=num./den;        %計算復數頻率響應
magH=20*log10(abs(H));   %轉換為分貝
phaH=angle(H);   %提取復數頻率響應的相位
subplot(2,1,1),plot(w/pi,magH);   %繪制幅頻響應
xlabel('w/pi');ylabel('|H(jw)|/dB');title('幅頻圖')
grid on;
subplot(2,1,2)
plot(w/pi,phaH*180/pi);  %繪制相頻響應
xlabel('w/pi');ylabel('相位角/^o');
title('相頻圖');grid on;