sssss.m

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

[y,fs,bits]=wavread('b.wav');

m=fft(y,512);  %取y的512点fft变换。
p=(0:255)/256*(fs/2);  

wp=0.2*pi;ws=0.3*pi;Rp=30;Rs=120; %Fs=fs;Ts=1/Fs;
wp1=2*fs*tan(wp/2);                  %将模拟指标转换成数字指标
ws1=2*fs*tan(ws/2); 
[N,Wc]=buttord(wp1,ws1,Rp,Rs,'s');     
[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);               
f=filter(bz,az,y);
F=fft(f,512); %求滤波后信号的fft变换
subplot(221);plot(y); title('采样语音时域图');%绘采样后信号的时域图
subplot(222);plot(f);title('滤波语音时域图'); %绘滤波后信号的频谱图
subplot(223);plot(p,abs([m(1:256)']));title('采样语音频谱图'); %画采样后信号的频谱图
subplot(224);plot(p,abs([F(1:256)']));title('滤波语音频谱图') %画滤波后信号的频谱图