ep626.m

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

%《数字信号处理--原理、实现及应用》第六章 例6.1.7 设计程序：ep617.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 设计巴特沃斯带模拟滤波器

wp=2*pi*[4000,7000];ws=2*pi*[2000,9000];Rp=1;As=20;   %设置带通滤波器指标参数
[N,wc]=buttord(wp,ws,Rp,As,'s');        %计算带通滤波器阶数N和3dB截止频率
[BB,AB]=butter(N,wc,'s');    %计算带通滤波器系统函数分子分母多项式系数
%绘制带通滤波器的损耗函数曲线
k=0:511;fk=0:15000/512:15000;wk=2*pi*fk;
Hk=freqs(BB,AB,wk);
subplot(2,2,1);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('频率(Hz)');ylabel('幅度(dB)')
axis([0,15000,-80,5]);title('带通滤波器损耗函数')