19t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 24 行

%最高f=30Hz，可取fs=100Hz，即t=0.01n
%s(n)=sin(0.1*pi*n)+sin(0.3*pi*n)+sin(0.6*pi*n);
%s(n)的样本取301点
%注意：这不是双线性变换法，是完全设计法，不过，效果一样。参考P
Rp=0.1;
Rs=40;
wp1=0.2*pi;
wp2=0.4*pi;
wn=[wp1,wp2]/pi;
n=4;
[b,a]=ellip(n,Rp,Rs,wn);%默认时表带通
[db,mag,pha,grd,w]=freqz_m(b,a);
subplot(3,1,1);
plot(w/pi,db);
axis([0,1,-100,5]);
n=0:300;
s=sin(0.1*pi*n)+sin(0.3*pi*n)+sin(0.6*pi*n);
subplot(312);
plot(n,s)
hold on
y=filter(b,a,s);
subplot(313);
plot(n,y)