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

jumboee

clear all
%-------------------------------------------------------------
N=500;
p=0.05;
3 _' e5 ^) P) u; D, a# Zf=1/50;
n=1:N;
$ _& B1 A4 l% o8 Y9 su=rand(1,N);
+ [8 A: e6 ]1 \' p/ s+ C: _$ r  Ms=sin(2*pi*f*0:N-1);
, P' Q9 x8 Q2 y  ^: x( Gx=u(1:N)+s;
! j: A; p* Q5 {# E$ Gsubplot(311)
/ J9 G# v$ e  |- h3 E0 w  Qplot(x(1:500));grid on;
& h1 l& s3 ^: _+ ?. Oylabel('x(n)')
%---------------------------------------------------------------
$ n2 l: D: \& Q. Dfp=200;
fs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
rs=35;  %阻带衰减不大于35dB
wp=2*fp/Fs;
  Q; g3 q# y: x( b( G7 E. u4 i4 B3 Ews=2*fs/Fs;
syms a w real
[nz,wn]=buttord(wp,ws,rp,rs);
& C7 B) _" X+ O; ~  d/ `- m% l[bz,az]=butter(nz,wp,'low');
, ]4 |1 w" W, J2 M3 q8 x# J: C% s[h,w]=freqz(bz,az,1024,Fs);
8 _4 {7 K& Y3 K( j( g$ eih=ifourier(h);
subplot(312)
6 }# z: y, Y: p' R) p* ~plot(x,h);grid;
ylabel('Low-pass DF')
%--------------------------------------------------------
/ K+ ?1 v, m( h% O( n/ Ry=fftfilt(h,x);
+ S: ?  f( C) d. x8 `% Jsubplot(313)
plot(y);

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