胡广书的教材第二章第12题相频响应的问题

[chen3bing]

我学过数字信号处理，当年用的是胡广书的教材。现在复习一下。

其中第二章第12题，求幅频响应和相频响应等。
3655

3656
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;

画出的图形
3657
高手给看看，哪里有问题？谢谢！