音频降噪Matlab仿真.docx

: .
音频降噪Matlab仿真
旗砌钥秘况孝
数字信号处理大作业班级：1401012
姓名：齐: .
音频降噪Matlab仿真
旗砌钥秘况孝
数字信号处理大作业班级：1401012
姓名：齐翔募学号：**********
输入信号的时域波形及其功率谱密度：
F
QFigure3|口||回|「S3|
FileEditViewInsertToolsDesktopWindowHelp华
□巳M鼠Q乃®曼曰国|■
012-O0
204060冏1QQ1201401601和200
(mp)，pnl-UEEzEn』lga)ds」«Mod
滤波后的皙音信号功率谱密度
+
叠加噪声后的音频信号的时域图形及功率谱密度：
Figure2
FileEditViewInsertToolsDesktopWindowHelp
臼9溢feW策O为受届，|@|曰医］■□加噪声后的语音倍号时域波形
.21Q1?
.
204060SO100120140160160200
ffip)SP5CBS巨En」ID・.」•鼻目d
o*1
加噪声后的语音信号功率谱密度

经过带通滤波器的音频信号的时域和功率谱密度：
程序解读:
clc;clearall;closeall;[wav,fs]=audioread('');t_end=1/fs*length(wav);%计算声音的时间长度Fs=50000;%仿真系统采样率t=1/Fs:1/Fs:t_end;%仿真系统采样时间点%利用插值函数将音频信号的采样率提升为Fs=50KHzwav=interp1([1/fs:1/fs:t_end],wav,t,'spline');%设计300Hz〜3400Hz的带通预滤波器H(z)[fenzi,fenmu]=butter(6,[3003400]/(Fs/2));nt=wgn(1,length(t),);%噪声nt=nt/(max(abs(nt)));%日一化噪声wav_noise=wav+nt;%对音频信号进行滤波wav_after=filter(fenzi,fenmu,wav_noise);figure(1);subplot(2,1,1);plot(wav(53550:53750));title('语音信号时域波形');axis([0200-]);subplot(2,1,2);psd