m2.m

FIR、IIR数字滤波器的Matlab仿真程序

fp=[480,520];
fs=[450,550];
wp=[480,520]*pi*2;
ws=[450,550]*pi*2;
rp=3;
rs=20;
[n,wn]=buttord(wp,ws,rp,rs,'s');
[b,a]=butter(n,wn,'s');
[ha,w]=freqs(b,a);
ma=abs(ha);
pha=unwrap(angle(ha));
subplot(421);
plot(w/(2*pi),ma);
subplot(423);
plot(w/(2*pi),pha);
fo=5000;
[bn,an]=impinvar(b,a,5000);
[hz,w]=freqs(bn,an);
mz=abs(hz);
phz=unwrap(angle(hz));
subplot(422);
plot(w,mz);
subplot(424);
plot(w,phz);
hi=impz(bn,an);
subplot(425);
plot(hi);
n=0:300;
t=n/fo;
x1=2*square(2*pi*500*t);
subplot(426);
plot(x1);
yi=conv(hi,x1);
subplot(427);
plot(yi);
y1=filter(bn,an,x1);
subplot(428);
plot(y1);