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

jumboee

clear all
%-------------------------------------------------------------
% r. Y6 T) t( s& f' uN=500;
p=0.05;
f=1/50;
5 H- O" D  i3 G0 ?$ Y$ in=1:N;
5 `) i7 O% |' mu=rand(1,N);
$ w, ]: [! v" j( X  W, Ps=sin(2*pi*f*0:N-1);
x=u(1:N)+s;
subplot(311)
9 S  @  [" Z5 p7 Tplot(x(1:500));grid on;
ylabel('x(n)')
2 {' X' J& I0 }%---------------------------------------------------------------
- v. c$ D- i* H" ]fp=200;
* h: M% [( l& x5 rfs=400;Fs=1200;
% X6 E  o) T  f! b7 brp=3;  %通带衰减不大于3dB
rs=35;  %阻带衰减不大于35dB
$ N8 x& X) L) ?) q# N- Wwp=2*fp/Fs;
) z( l; u5 }  ]: {# }% kws=2*fs/Fs;
syms a w real
[nz,wn]=buttord(wp,ws,rp,rs);
' Z( J, Z/ S# c& g" y1 g0 [# B& j[bz,az]=butter(nz,wp,'low');
# c* w& ?, V' ?. }% _  k2 ~( p[h,w]=freqz(bz,az,1024,Fs);
( c6 ^: z" N, l7 @ih=ifourier(h);
subplot(312)
  x7 v2 o( R$ b4 H6 U! I+ R- Bplot(x,h);grid;
( U) K5 F/ ^9 }4 v3 i% Dylabel('Low-pass DF')
/ G" F; p1 s0 I' i9 G+ C%--------------------------------------------------------
: h# b1 N7 Q" C) |+ ]y=fftfilt(h,x);
subplot(313)
2 t' d/ D( p/ S; Splot(y);
4 n! s! E; u/ a2 A* Z' Y
画出的图不对，好像是逆傅里叶变换用的不太对似的。