《数字的信号处理实验的.doc》由会员分享,可在线阅读,更多相关《数字的信号处理实验的.doc(4页珍藏版)》请在课桌文档上搜索。
1、wordx,fs=audioread(000.wav);sound(x,fs)n=length(x);x_p=fft(x,n);f=fs*(0:n/2-1)/n;figure(1)subplot(2,1,1); plot(x); title(加噪语音信号采样后的时域波形); xlabel(时间轴)ylabel(幅值A) subplot(2,1,2); plot(f,abs(x_p(1:n/2); title(加噪语音信号采样后的频谱图); xlabel(频率Hz); ylabel(频率幅值);L=length(x);noise=0.1*randn(L,2);x_z=x+noise;sound(
2、x_z,fs)n=length(x);x_zp=fft(x_z,n);f=fs*(0:n/2-1)/n;figure(2)subplot(2,1,1);plot(x_z); title(加噪语音信号时域波形);xlabel(时间轴)ylabel(幅值A)subplot(2,1,2); plot(f,abs(x_zp(1:n/2);title(加噪语音信号频谱图); xlabel(频率Hz); ylabel(频率幅值);L=length(x) ;noise=0.1*randn(L,2)x_z=x+noise;sound(x_z,fs)t=0:0.0001:10;f1=3800;%噪声信号频率y1
3、=0.5*sin(2*pi*f1*t);%噪声信号y2=y1(1,1:5000);%取噪声信号的5000点长度1,1:5000因为y1为列向量y=x+y2;%先对噪声信号转置再把噪声信号参加原始信号Y=fft(y2,4096);%求噪声信号频谱magY=abs(Y);%把傅里叶变换后的复数值取模subplot(2,2,1);plot(y2);title(噪声信号波形);%绘制波形subplot(2,2,2);plot(f,magY(1:2048);title(噪声信号频谱);%绘制波形subplot(2,2,3);plot(y);title(加噪后信号波形);%绘制波形Y1=fft(y,4096);%对加噪后的信号做4096点FFT变换magY1=abs(Y1);%把傅里叶变换后的复数值取模subplot(2,2,4);plot(f,magY1(1:2048);title(加噪后信号频谱);%绘制波形4 / 4