音频降噪Matlab仿真

1、音频降噪Matlab仿真数字信号处理大作业班级： 1401012_姓名： 齐翔奡 _学 号入信号的时域波形及其功率谱密度：FQ Figure 3| a | B | E3 |File Edit View Insert Tools Desktop Window Help爼曰Pd|头卩哽貳T(a|日匡l| 口滤波后的语音信号时域派形20 4Q 60 SO 1QQ 1201401602Q0.21 oiOL0 12-O0(HP)空 Pru-UB 如乏一 Enlo#d/3心芬00滤波后的谊音信号功率谱密度Frequency斗x10叠加噪声后的音频信号的时域图形及功率谱密 度:Fi

2、le Edit View Insert Tools Desktop Win dow Helpde加噪声后的语音信号时域波形204060 SO 1001201401601S02002-10 12- o nu 00.ffim PECBE 运 En75wdn加噪声后的语音信号功率谱密度0 51.5Frequency22.5斗o*1x10经过带通滤波器的音频信号的时域和功率谱密度:Q Figure 1File Edit View Insert Tools Desktop Window Helpna afe v - ” 厘亡 n语音佶号时域液形204060 SO 1001201401601S02002-

3、10 12- o nu 00.(tDpj pnlcBE 运 En 二owdtn语音信号功率谱密度0 51.5Frequency22.5x10斗程序解读:clc;clear all;close all;wav,fs=audioread(GDGvoice8000.wav);t_en d=1/fs *le ngth(wav); % 计算声音的时间 长度Fs=50000;%仿真系统采样t=1/Fs:1/Fs:t_e nd;间点%仿真系统采样时%利用插值函数将音频信号的采样率提升为Fs=50KHzwav=i nterp1(1/fs:1/fs:t_e nd,wav,t,spli ne);%设计300Hz3

4、400Hz的带通预滤波器H(z) fen zi,fe nmu=butter(6,300 3400/(Fs/2);nt = wgn(1,length(t),0.1); % 噪声nt=n t/(max(abs( nt);%归一化噪声 wav_ no ise = wav + nt;%对音频信号进行滤波wav_after = filter(fe nzi,fe nmu,wav_ noise);figure(1);subplot(2,1,1); plot(wav(53550:53750); title(语音信号时域波形);axis(0 200 -0.3 0.3);subplot(2,1,2); psd(w

5、av, 10000, Fs);title(语音信号功率谱密度);axis(0 25000 -20 10);figure(2);subplot(2,1,1); plot(wav_noise(53550:53750);title(加噪声后的语音信号时域波形);axis(0 200 -0.3 0.3);subplot(2,1,2);psd(wav_noise, 10000, Fs);title(加噪声后的语音信号功率谱密度);axis(0 25000 -20 10);figure(3);subplot(2,1,1); plot(wav_after(53550:53750);title(滤波后的语音信号时域波形);axis(0 200 -0.3 0.3);subplot(2,1,2);psd(wav_after, 10000, Fs);title(滤波后的语音信号功率谱密度);axis(0 25000 -20 10);sou nd(wav/max(wav), Fs);% 播放解调音频soun