数字滤波器程序问题，麻烦大家帮看看

jumboee

clear all
%-------------------------------------------------------------
N=500;
8 `! K$ ^3 \  Y3 F  w' l9 gp=0.05;
6 B$ g% q6 u) ?' D" ~* Tf=1/50;
  u9 ~  D' R3 tn=1:N;
u=rand(1,N);
6 Q, F  n2 J" _& Ks=sin(2*pi*f*0:N-1);
x=u(1:N)+s;
8 n* ~6 V& Y- O9 }) d/ t4 L& |subplot(311)
plot(x(1:500));grid on;
$ I! x4 Z" O6 y' ?% D7 `3 ]8 ~ylabel('x(n)')
%---------------------------------------------------------------
& I; w6 _2 P0 l- ifp=200;
fs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
rs=35;  %阻带衰减不大于35dB
wp=2*fp/Fs;
  x  l- P4 p6 ^0 x& n0 u1 Gws=2*fs/Fs;
syms a w real
0 |! J/ U1 N8 T) s( o  V3 O[nz,wn]=buttord(wp,ws,rp,rs);
. d2 k- A# f, H2 T* a[bz,az]=butter(nz,wp,'low');
[h,w]=freqz(bz,az,1024,Fs);
. }/ g9 n+ H/ Z  ~  B1 Jih=ifourier(h);
subplot(312)
2 k" w' ]  i2 O8 c6 _plot(x,h);grid;
. B+ X+ q  F5 r: J$ M+ Vylabel('Low-pass DF')
%--------------------------------------------------------
y=fftfilt(h,x);
subplot(313)
plot(y);

画出的图不对，好像是逆傅里叶变换用的不太对似的。