ep624.m

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

%《数字信号处理--原理、实现及应用》第六章 例6.1.4 设计程序：ep614.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 椭圆滤波器设计程序
wp=2*pi*3000;ws=2*pi*12000;Rp=0.1;As=60;   %设置指标参数
[N,wpo]=ellipord(wp,ws,Rp,As,'s');      %计算椭圆低通模拟滤波器阶数和通带边界频率
[B,A]=ellip(N,Rp,As,wpo,'s');           %计算低通模拟滤波器系统函数系数
fk=0:14000/512:14000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
subplot(2,2,1);
plot(fk/1000,20*log10(abs(Hk)));grid on
xlabel('频率(kHz)');ylabel('幅度(dB)')
axis([0,14,-80,5])