《MATLAB程序原始语音信号波形与加噪信号波形》由会员分享,可在线阅读,更多相关《MATLAB程序原始语音信号波形与加噪信号波形(2页珍藏版)》请在金锄头文库上搜索。
1、%用 MATLAB对原始语音信号进行分析,画出它的时域波形和频谱fs=22050;x1 =wavread(C:Documents and SettingsAdministrator桌面1.wav ) ;sound(x1,22050); %播放语音信号y1=fft(x1,1024); %对信号做 1024点 FFT变换f=fs*(0:511)/1024;figure(1)plot(x1) %做原始语音信号的时域图形title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1) %绘制原始语音信号的频率响应图title(频率响应图)f
2、igure(3)subplot(2,1,1);plot(abs(y1(1:512) %做原始语音信号的 FFT频谱图title(原始语音信号 FFT频谱)subplot(2,1,2);plot(f,abs(y1(1:512);title(原始语音信号频谱)xlabel(Hz);ylabel(幅值);%给原始的语音信号加上一个高频余弦噪声,频率为(56Khz 任选) 。对加噪后的语音进行分析,并画出其信号时域和频谱图。fs=22050;x1=wavread(C:Documents and SettingsAdministrator桌面1.wav );f=fs*(0:511)/1024;t=0:1
3、/fs:(length(x1)-1)/fs;Au=0.03;d=Au*cos(2*pi*5000*t);x2=x1+d;sound(x2,22050);wavwrite(x2,8000,16, C:Documents and SettingsAdministrator桌面2.wav);y2=fft(x2,1024);figure(4);plot(t,x2)title(加噪后的信号);xlabel(time n);ylabel(fuzhi n);figure(5)subplot(2,1,1);plot(f,abs(y1(1:512);title(原始语音信号频谱);xlabel(Hz);ylabel(fuzhi);subplot(2,1,2);plot(f,abs(y2(1:512);title(加噪后的信号频谱);xlabel(Hz);ylabel(fuzhi);