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

jumboee

clear all
1 l* J1 X0 o4 W8 G. ?- C. z%-------------------------------------------------------------
N=500;
* U% y0 s$ z4 Y2 J# H" Jp=0.05;
' |! L. D5 E- J7 r- @f=1/50;
n=1:N;
8 l  M' c" ?- m7 a$ _& C+ ?6 Du=rand(1,N);
s=sin(2*pi*f*0:N-1);
2 |( O9 _: `3 ^9 N2 r8 G& S# {x=u(1:N)+s;
subplot(311)
1 c( T6 ?5 D- ~! z* yplot(x(1:500));grid on;
& n7 e; J" H8 Q1 x6 e+ i& l4 qylabel('x(n)')
%---------------------------------------------------------------
fp=200;
fs=400;Fs=1200;
3 y/ \( Q6 f9 u+ ^. r. d0 orp=3;  %通带衰减不大于3dB
rs=35;  %阻带衰减不大于35dB
wp=2*fp/Fs;
3 t/ v+ D6 [9 _8 V, wws=2*fs/Fs;
( d' F8 i  ]. @syms a w real
8 A3 I; Y, K/ v3 f[nz,wn]=buttord(wp,ws,rp,rs);
[bz,az]=butter(nz,wp,'low');
# l* m$ y( F- q* g, ^[h,w]=freqz(bz,az,1024,Fs);
' k* l  @& {  D  P9 _ih=ifourier(h);
subplot(312)
: R' s2 Q. Y+ e. G5 [, {plot(x,h);grid;
9 K( p  H6 _) r4 q* K9 x  F5 Xylabel('Low-pass DF')
%--------------------------------------------------------
y=fftfilt(h,x);
subplot(313)
plot(y);

" L8 Y1 ^" r# L, \画出的图不对，好像是逆傅里叶变换用的不太对似的。