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

jumboee

clear all
0 s; W, q, k  @$ g%-------------------------------------------------------------
. v2 l' Q( W& f% y: j7 L& ?N=500;
) @: v# o- L: g  a, y4 Ap=0.05;
9 E# p% q1 c* T8 P. }, l/ ]f=1/50;
n=1:N;
u=rand(1,N);
s=sin(2*pi*f*0:N-1);
x=u(1:N)+s;
3 K6 Y4 @  @# Hsubplot(311)
plot(x(1:500));grid on;
ylabel('x(n)')
. V7 _9 ]! M8 D& D! j! X" M%---------------------------------------------------------------
; V9 t( C) S3 G) ^- sfp=200;
) @9 ^* D4 T; Ifs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
! q# Y' Q. P' _' Q$ m- N: ^+ Ers=35;  %阻带衰减不大于35dB
wp=2*fp/Fs;
" ]& A- W8 B# U7 {0 Zws=2*fs/Fs;
syms a w real
[nz,wn]=buttord(wp,ws,rp,rs);
[bz,az]=butter(nz,wp,'low');
[h,w]=freqz(bz,az,1024,Fs);
ih=ifourier(h);
subplot(312)
" N' P  E. b  H  `3 _+ wplot(x,h);grid;
ylabel('Low-pass DF')
! R. A9 |( [$ M& N' Q& I3 _; p%--------------------------------------------------------
y=fftfilt(h,x);
subplot(313)
1 j( O% b. ~& i- J( l0 nplot(y);

: A7 `8 K, i9 e8 B) f" T画出的图不对，好像是逆傅里叶变换用的不太对似的。