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

jumboee

clear all
%-------------------------------------------------------------
, `- o+ I' I5 V  n% {/ D8 MN=500;
3 Q5 H. E# N4 ^9 f8 C3 y; Xp=0.05;
f=1/50;
" F6 `9 f& ?8 P1 @- B6 P$ u4 g; Zn=1:N;
& ]$ d. _5 j& c3 I7 x6 ]u=rand(1,N);
s=sin(2*pi*f*0:N-1);
* W2 i% f# t/ B$ ^- U2 Ux=u(1:N)+s;
2 y0 Y3 t9 ?7 l7 o4 [subplot(311)
plot(x(1:500));grid on;
ylabel('x(n)')
%---------------------------------------------------------------
% F5 k  i& H3 G% O6 G- [4 B4 n0 K: {fp=200;
: r! ~4 A( x0 V$ B; o  T7 ]$ cfs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
1 u8 @" m1 _! K& y3 Grs=35;  %阻带衰减不大于35dB
. Q5 x0 w6 ]' V9 |wp=2*fp/Fs;
ws=2*fs/Fs;
, j2 O! J) u* K0 s  q. @0 \syms a w real
* V- }# T6 p" f/ }[nz,wn]=buttord(wp,ws,rp,rs);
* d8 {, {6 [$ V% e! i[bz,az]=butter(nz,wp,'low');
2 c5 G* L$ r7 m- P) q- s) _9 }[h,w]=freqz(bz,az,1024,Fs);
ih=ifourier(h);
subplot(312)
& ~. q, y' T6 D6 W6 w# [3 fplot(x,h);grid;
ylabel('Low-pass DF')
%--------------------------------------------------------
4 Z) T5 e+ x' i2 X& Q; N1 O" {y=fftfilt(h,x);
subplot(313)
plot(y);

: \8 f  C$ q# t8 ~1 {画出的图不对，好像是逆傅里叶变换用的不太对似的。