untitled3.m

来自「我的数字信号处理课程设计」· M 代码 · 共 16 行

[y,fs,bits]=wavread('b.wav');
wp=0.1*pi;
ws=0.3*pi;
Rp=30;
Rs=120;
Fs=fs;
Ts=1/Fs;
wp1=2/Ts*tan(wp/2);                  %将模拟指标转换成数字指标
ws1=2/Ts*tan(ws/2); 
[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');      %选择滤波器的最小阶数
[Z,P,K]=buttap(N);                     %创建butterworth模拟滤波器
[Bap,Aap]=zp2tf(Z,P,K);
[b,a]=lp2lp(Bap,Aap,Wn);   
[bz,az]=bilinear(b,a,Fs);         %用双线性变换法实现模拟滤波器到数字滤波器的转换
[H,W]=freqz(bz,az);                 %绘制频率响应曲线
plot(W*Fs/(2*pi),abs(H))