ep622.m

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

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

wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;   %设置滤波器参数
[N,wc]=buttord(wp,ws,Rp,As,'s');        %计算滤波器阶数N和3dB截止频率
[B,A]=butter(N,wc,'s');                 %计算滤波器系统函数分子分母多项式系数
k=0:511;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,-40,5])