untitled.m

我的数字信号处理课程设计

[y,fs,bits]=wavread('b.wav');
wp1=8000;ws1=16000;rp=0.5;rs=40;
%wp1=2*fs*tan(wp/2);
%ws1=2*fs*tan(ws/2);
[N,Wc]=buttord(wp1,ws1,rp,rs);
[Z,P,K]=buttap(N);                  %创建butterworth模拟滤波器
[Bap,Aap]=zp2tf(Z,P,K);
[b,a]=lp2lp(Bap,Aap,Wc);   
[bz,az]=bilinear(b,a,fs);           %用双线性变换法实现模拟滤波器到数字滤波器的转换
[H,W]=freqz(bz,az);                 %绘制频率响应曲线
plot(W,abs(H));%plot(W*fs/(2*pi),abs(H))