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

jumboee

clear all
% J1 h" l( u7 v5 K%-------------------------------------------------------------
! A! Z% C' Y7 B1 @( g1 i8 @3 B/ ]N=500;
p=0.05;
f=1/50;
& c6 z: K& a0 z: ]9 C0 bn=1:N;
u=rand(1,N);
& |: d( a- u" q& I0 B0 G/ B5 ~6 `s=sin(2*pi*f*0:N-1);
( P8 y" p& m5 ox=u(1:N)+s;
& G/ E6 @' B/ V) Gsubplot(311)
% s% A; g% t' G  H# gplot(x(1:500));grid on;
  Z$ w. `; D4 U5 n* q5 |% I( Rylabel('x(n)')
7 D7 ^, d2 @; Q' S3 k9 o# t9 e7 Y- u%---------------------------------------------------------------
fp=200;
fs=400;Fs=1200;
" ?0 n5 K3 [# z/ {1 O. Nrp=3;  %通带衰减不大于3dB
2 d3 ^  ]" `- B. Q2 Yrs=35;  %阻带衰减不大于35dB
% ^# A! G' ^8 F  D2 [wp=2*fp/Fs;
) _$ D* G0 [- E' h: f2 p, e/ Kws=2*fs/Fs;
! {" p$ _/ `: \  esyms a w real
4 e! [1 @- R6 w6 m) N[nz,wn]=buttord(wp,ws,rp,rs);
+ g8 M2 p4 T1 |8 i1 T[bz,az]=butter(nz,wp,'low');
) H9 J6 U+ j8 m9 q( o. H. o0 |[h,w]=freqz(bz,az,1024,Fs);
ih=ifourier(h);
- i7 l! V$ L9 c$ o2 D; d' C/ l% ]subplot(312)
0 J' g  ]4 I  hplot(x,h);grid;
' o$ S% l) ]. x, G9 Pylabel('Low-pass DF')
3 V4 u- v" y$ y2 d* A: b: c* F- T%--------------------------------------------------------
y=fftfilt(h,x);
subplot(313)
plot(y);
! K5 Z+ H+ o6 m4 b7 O$ }( u
" B, d: u4 c: e9 O5 e画出的图不对，好像是逆傅里叶变换用的不太对似的。