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

jumboee

clear all
/ [1 R' S% f: F+ k%-------------------------------------------------------------
N=500;
p=0.05;
' ?( U% E. E% a& ^& [9 Cf=1/50;
n=1:N;
6 ^( x/ ]% @7 R; O' ju=rand(1,N);
" |- t3 o  k1 Z) T+ o1 {s=sin(2*pi*f*0:N-1);
0 F1 l5 m$ N: R4 w8 Rx=u(1:N)+s;
$ z4 @' ]) V; G7 E* C/ ~subplot(311)
, T5 N) Y# P7 I& Y% Xplot(x(1:500));grid on;
ylabel('x(n)')
7 m$ D1 \# A5 ?+ M: B, U: [0 ~%---------------------------------------------------------------
0 v1 G5 x# I" z. pfp=200;
fs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
rs=35;  %阻带衰减不大于35dB
5 p: q1 f0 R* W' c/ L4 ^4 v. ]: {wp=2*fp/Fs;
ws=2*fs/Fs;
syms a w real
0 `4 j3 o; s3 c- I% I2 E[nz,wn]=buttord(wp,ws,rp,rs);
  W. V. C: a' C/ c4 X1 _[bz,az]=butter(nz,wp,'low');
: Q8 N2 @; J( m. U7 P: }4 V1 t[h,w]=freqz(bz,az,1024,Fs);
4 }; B- S, n3 i3 ?7 ^ih=ifourier(h);
7 B  \* |+ J4 s7 [- ^subplot(312)
plot(x,h);grid;
# I  H) I- v$ D% u, |4 o! A% vylabel('Low-pass DF')
%--------------------------------------------------------
y=fftfilt(h,x);
  w' {1 c8 P# L8 n2 Y9 e. Dsubplot(313)
plot(y);

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