ep625.m

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

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

wp=1;ws=4;Rp=0.1;As=40;   %设置滤波器指标参数
[N,wc]=buttord(wp,ws,Rp,As,'s');        %计算滤波器G(p)阶数N和3dB截止频率
[B,A]=butter(N,wc,'s'); %计算低通滤波器G(p)系统函数分子分母多项式系数
wph=2*pi*4000;          %高通模拟滤波器通带边界频率
[BH,AH]=lp2hp(B,A,wph); %低通到高通转换
%绘制低通滤波器G(p)的损耗函数曲线
wk=0:0.01:10;
Hk=freqs(B,A,wk);
subplot(2,2,1);
plot(wk,20*log10(abs(Hk)));grid on
xlabel('归一化频率');ylabel('幅度(dB)')
axis([0,10,-80,5]);title('(a) 归一化低通损耗函数')
%绘制高通滤波器的损耗函数曲线
k=0:511;fk=0:6000/512:6000;wk=2*pi*fk;
Hk=freqs(BH,AH,wk);
subplot(2,2,2);
plot(fk,20*log10(abs(Hk)));grid on
xlabel('频率(Hz)');ylabel('幅度(dB)')
axis([0,6000,-80,5]);title('(b) 高通滤波器损耗函数')