语音信号的频谱分析与FIR滤波器的设计.docVIP

目?????录 设计 （2）时域信号的FFT分析....................(4) （3）基于窗函数的FIR数字滤波器设计法......(4) 2.实验工具..................................(5) 三．仿真程序及频谱分析 1.语音信号的频谱分析.........................(5) 2.FFT频谱分析...............................（7） 3.FIR数字低通滤波器的设计....................(8) 4.FIR数字高通滤波器的设计...................(10) 5.FIR数字带通滤波器的设计...................(13) 6.保存滤波后的声音文件.......................(15) 五． FFT频谱分析 程序： [x1,fs,bits]=wavread(sy.wav); %读取语音信号的数据，赋给变量x1 y1=fft(x1,1024); %对采样后信号做1024点FFT变换 f=fs*(0:511)/1024; figure(1) subplot(2,1,1); plot(f,abs(y1(1:512))); title(原始语音信号频谱) xlabel(频率/Hz); ylabel(幅值) subplot(2,1,2); plot(abs(y1(1:1024))) %采样后信号的FFT频谱图 title(原始语音信号FFT频谱) xlabel(点数N); ylabel(幅值); 频谱图 3.FIR数字低通滤波器的设计 （1）程序： fs=10000; x1=wavread(sy.wav); wp=2*pi*1000/fs; ws=2*pi*1200/fs; Rp=1; Rs=100; wdelta=ws-wp; N=ceil(8*pi/wdelta); %取整 wn=(wp+ws)/2; [b,a]=fir1(N,wn/pi,hamming(N+1)); %选择窗函数，并归一化截止频率 figure(1) freqz(b,a,512); title(FIR低通滤波器); f2=filter(b,a,x1); figure(2) subplot(2,1,1) plot(x1) title(FIR低通滤波器滤波前的时域波形); subplot(2,1,2) plot(f2); title(FIR低通滤波器滤波后的时域波形); sound(f2,7000); %播放滤波后的语音信号 F0=fft(f2,1024); f=fs*(0:511)/1024; figure(3) y2=fft(x1,1024); subplot(2,1,1); plot(f,abs(y2(1:512))); title(FIR低通滤波器滤波前的频谱) xlabel(频率/Hz); ylabel(幅值); subplot(2,1,2) F2=plot(f,abs(F0(1:512))); title(FIR低通滤波器滤波后的频谱) xlabel(频率/Hz); ylabel(幅值); 频谱图 FIR高通滤波器的设计 （1）程序 fs=22050; x1=wavread(sy.wav); wp=2*pi*5000/fs; ws=2*pi*4800/fs; Rp=1; Rs=100; wdelta=wp-ws; N=ceil(8*pi/wdelta); %取整 wn=(wp+ws)/2; [b,a]=fir1(N,wn/pi,high); figure(1) freqz(b,a,512); title(FIR高通滤波器); f2=filter(b,a,x1); figure(2) subplot(2,1,1) plot(x1) title(FIR高通滤波器滤波前的时域波形); subplot(2,1,2) plot(f2); title(FIR高通滤波器滤波后的时域波形); sound(f2,20000); %播放滤波后的语音信号 F0=fft(f2,1024); f=fs*(0:511)/1024; figure(3) y2=fft(x1,1024); subplot(2,1,1); plot(f,abs(y2(1:512))); title(FIR高通滤波器滤波前的频谱) xlabel(频率/Hz); yla