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

jumboee

clear all
/ L) O" v& p1 Q; I: S7 r% ]% H5 }%-------------------------------------------------------------
N=500;
, T4 ~; P) x0 ?) d3 U. J1 gp=0.05;
f=1/50;
6 D9 o$ u. \. h7 T1 v; N; |0 E- {n=1:N;
u=rand(1,N);
s=sin(2*pi*f*0:N-1);
x=u(1:N)+s;
8 q7 |& \7 [( f# S0 ^. qsubplot(311)
plot(x(1:500));grid on;
* q  ]: X+ {; \' ^3 u! Fylabel('x(n)')
%---------------------------------------------------------------
8 Z1 h# W! ?$ \0 b) Tfp=200;
fs=400;Fs=1200;
rp=3;  %通带衰减不大于3dB
, v# V% F( C5 w; A/ A6 P, @  Crs=35;  %阻带衰减不大于35dB
wp=2*fp/Fs;
ws=2*fs/Fs;
( {( S' X. K/ P' O5 E5 ?0 p8 qsyms a w real
[nz,wn]=buttord(wp,ws,rp,rs);
[bz,az]=butter(nz,wp,'low');
& Q( F( ^" y! D  r; a. N4 g[h,w]=freqz(bz,az,1024,Fs);
ih=ifourier(h);
8 o) M6 R5 d: B7 d2 Nsubplot(312)
plot(x,h);grid;
ylabel('Low-pass DF')
9 |/ \# F0 @  e: b1 Z. v; Q1 l8 I%--------------------------------------------------------
2 ]: N6 Y, x3 Z( R4 Cy=fftfilt(h,x);
% z/ c* a  w& q- L" X. _# C9 \subplot(313)
plot(y);
7 R5 l; u: U: J* U
' _1 D; F; D3 S# |# {$ Q+ q1 I画出的图不对，好像是逆傅里叶变换用的不太对似的。