exp3_16_1.m

来自「是Digital signal processing--using Matlab」· M 代码 · 共 16 行

%调用函数freqz，其参数分别为分子、分母系数，采样点个数
a=[1 -1.76 1.1829 -0.2781];
b=[0.0181 0.0543 0.0543 0.0181];
m=0:length(a)-1;
n=0:length(b)-1;
k=0:1:500;
w=pi*k/500;
[h,w]=freqz(b,a,500);
realh=real(h);
imagh=imag(h);
magh=abs(h);
angh=angle(h);
subplot(2,2,1),plot(w/pi,realh),title('realh');grid on;
subplot(2,2,2),plot(w/pi,imagh),title('imagh');grid on;
subplot(2,2,3),plot(w/pi,magh),title('magh');grid on;
subplot(2,2,4),plot(w/pi,angh),title('angh');grid on;