matlab三维图格点,基于MATLAB的语音信号的采集与处理详解

5-3：巴特沃斯滤波器的设计：

wp=0.25*pi; %通带截止频率

ws=0.3*pi; %阻带截止频率

Rp=1; %通带最大衰减(db)

Rs=15; %阻带最大衰减(db)

Fs=22050;

Ts=1/Fs;

wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标

ws1=2/Ts*tan(ws/2); %将模拟指标转换成数字指标

[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s'); %选择滤波器的最小阶数(估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc)

[Z,P,K]=buttap(N); %创建Butterworth低通滤波器原型

[Bap,Aap]=zp2tf(Z,P,K); %将零极点增益转换为普遍分子，分母

[b,a]=lp2lp(Bap,Aap,Wn); %将普遍的分子和分母转换为以Wn为截止频率

[bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换(模拟转换为数字)

[H,W]=freqz(bz,az); %求频率响应

figure(6)

plot(W*Fs/(2*pi),abs(H)) %绘制Butterworth低通滤波器频率响应曲线

grid %添加图格

xlabel('频率／Hz')

ylabel('频率响应幅度')

title('Butterworth')