fig622.m

高西全2009年最新版<数字信号处理>源码程序

%fig622.m
wp=2*pi*2000;ws=2*pi*3000;Rp=0.1;As=60;   %设置指标参数
%设计椭圆滤波器
[N,wpo]=ellipord(wp,ws,Rp,As,'s');      %计算椭圆低通模拟滤波器阶数和通带边界频率
[B,A]=ellip(N,Rp,As,wpo,'s');           %计算低通模拟滤波器系统函数系数
k=0:1023;fk=0:6000/1024:8000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
subplot(2,2,1);
plot(fk,abs(Hk));grid on
xlabel('f(Hz)');ylabel('|H_a(jf)|')
axis([0,8000,0,1.1])
subplot(2,2,2);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('f(Hz)');ylabel('A(f)(dB)')
axis([0,8000,-100,5])