我学过数字信号处理，当年用的是胡广书的教材。现在复习一下。
其中第二章第12题，求幅频响应和相频响应等。
4.jpg
5.jpg
MATLAB写了个脚本，运行以后幅频响应还可以，相频响应就差的很多，解卷绕也不行。
我的代码：
b=[0,0.5,-0.5,0];
a=[1,-(sqrt(2)+0.9),(0.9*sqrt(2)+1),-0.9];
[H,w]=freqz(b,a,64,'whole',1);
Hr=10*log10(abs(H));
Hphase=angle(H)%;Hphase=unwrap(Hphase);
n=1:64;
subplot(221);plot(n,Hr);grid on;
subplot(225);plot(n,Hphase);grid on;
画出的图形
6.jpg
高手给看看，哪里有问题？谢谢！