m3.m

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

fp=[480,520];
fs=[450,550];
rp=3;
rs=20;
fo=10000;
wp=2*pi*fp/fo;
ws=2*pi*fs/fo;
[n,wn]=cheb1ord(wp/pi,ws/pi,rp,rs);
[b,a]=cheby1(n,rp,wn);
[h,w]=freqz(b,a,128,10000)
mag=abs(h);
pha=unwrap(angle(h));
subplot(321);
plot(w,mag);
subplot(322);
plot(w,pha);
hi=impz(b,a);
subplot(324);
plot(hi);
n=0:500;
t=n/fo;
x1=2*square(2*pi*500*t);
subplot(323);
plot(t,x1);
yi=conv(hi,x1);
subplot(326);
plot(yi);
y1=filter(b,a,x1);
subplot(325);
stem(y1);