ep625b.m

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

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

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