加噪语音信号的滤波matlab.docVIP

《信号与系统》课程设计——加噪语音信号的滤波 组员： 1 读入语音信号Q2，并对其进行时域分析 y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); k=linspace(0,N/44100,N); plot(k,y); title(Q2信号的时域分析图像); xlabel(时间 S); ylabel(幅度); 图像 2.对Q2做频域分析： y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); ff = fft(y); index=linspace(-22050,22050,N); figure(1) plot(index , abs(fftshift(ff))); title(Q2信号的频域波形); xlabel(频率 HZ); ylabel(幅度); 从频谱图像上可以看出，信号的频谱大致在8Khz以下，而且具有较大的直流偏量，信号中低频的成分比较多，越到高频的地方幅度越小。 二． 播放Q2音频，然后加上不同幅度的噪声，再播放加上噪声的音频。 y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); wavplay(y,44100); noise=0.1*randn(size(y)); signal=noise+y; wavplay(signal,44100); noise2=0.2*randn(size(y)); signal2=noise2+y; wavplay(signal2,44100); 我们可以清晰的听出，噪声的幅度越大时，语音信号的质量越差。 对于noise信号进行时域和频域分析 y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); noise=0.1*randn(size(y)); N=length(y); k=linspace(0,N/44100,N); plot(k,noise) title(noise时域图); xlabel(‘时间); ylabel(幅度); y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); noise=0.1*randn(size(y)); ff = fft(noise); index=linspace(-22050,22050,N); figure(3) plot(index , abs(fftshift(ff))); title(noise信号的频域波形); xlabel(频率 HZ); ylabel(幅度); noise信号的频谱图和Q2信号的频谱图相比较，噪声信号从高频到低频都存在，而且幅度大小基本相同，二语音信号基本只有低频信号存在，幅度直接又有明显差别 对加噪语音信号signal进行时域和频域分析： y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); noise=0.1*randn(size(y)); signal=noise+y; k=linspace(0,N/44100,N); figure（1）,plot(k,signal) title(signal的时域图); xlabel(‘时间’); ylabel(幅度); noise2=0.2*randn(size(y)); signal2=noise2+y; figure（2）,plot(k,signal2) title(signal2的时域图); xlabel(‘时间); ylabel(幅度); 从signal和signal2的时域图中可以看出，所加噪声的幅度越大，加噪后的语音信号时域图的幅度越偏离Q2信号时域图对应时间点的幅度 y=wavread(C:\Users\Administrator\Desktop\kecheng\Q2.wav); N=length(y); noise=0.1*randn(size(y)); signal=noise+y; noise2=0.2*randn(size(y)); signal2=noise2+y; ff = fft(signal); index=linspace(-22050,22050,N); figure(1) plot(index , abs(fftshift(ff)));