有一个这样的滤波器，b=[1 -2*cos(0.1*pi) 1] ; a=[1 -2*0.99*cos(0.1*pi) 0.99*0.99];
请问，怎样得到它的频率响应特性曲线？
我是这样做的，
fs=1000;
b=[1 -2*cos(0.1*pi) 1] ;
a=[1 -2*0.99*cos(0.1*pi) 0.99*0.99];
[h,f]=freqz(b,a,100,fs);
mag=ans(h);
ph=angle(h);
ph=ph*180/pi;
subplot(211),plot(f,mag);grid%幅度特性曲线
subplot(212),plot(f,ph);grid%相位特性曲线

可是我的结果和别人文章里的结果，差别非常大，尤其是相位特性曲线根本不是一回事。请问这是怎么回事啊？
应该怎么显示出滤波器的频率响应曲线呢？