音頻降噪Matlab仿真

1、-作者xxxx-日期xxxx音頻降噪Matlab仿真【精品文檔】數(shù)字信號(hào)處理大作業(yè)班 級(jí)： 1401012_ 姓 名： 齊翔奡 _ 輸入信號(hào)的時(shí)域波形及其功率譜密度：疊加噪聲后的音頻信號(hào)的時(shí)域圖形及功率譜密度：經(jīng)過(guò)帶通濾波器的音頻信號(hào)的時(shí)域和功率譜密度：程序解讀：clc;clear all;close all; wav,fs=audioread('GDGvoice8000.wav');t_end=1/fs *length(wav); % 計(jì)算聲音的時(shí)間長(zhǎng)度Fs=50000; % 仿真系統(tǒng)采樣率t=1/Fs:1/Fs:t_end; % 仿真系統(tǒng)采樣時(shí)間點(diǎn)% 利用插值函數(shù)將音頻信號(hào)

2、的采樣率提升為Fs=50KHz wav=interp1(1/fs:1/fs:t_end,wav,t,'spline');% 設(shè)計(jì)300Hz3400Hz的帶通預(yù)濾波器H(z) fenzi,fenmu=butter(6,300 3400/(Fs/2);nt = wgn(1,length(t),0.1); % 噪聲nt=nt/(max(abs(nt);%歸一化噪聲wav_noise = wav + nt;% 對(duì)音頻信號(hào)進(jìn)行濾波 wav_after = filter(fenzi,fenmu,wav_noise);figure(1); subplot(2,1,1); plot(wav(5

3、3550:53750); title('語(yǔ)音信號(hào)時(shí)域波形');axis(0 200 -0.3 0.3);subplot(2,1,2); psd(wav, 10000, Fs);title('語(yǔ)音信號(hào)功率譜密度');axis(0 25000 -20 10);figure(2); subplot(2,1,1); plot(wav_noise(53550:53750); title('加噪聲后的語(yǔ)音信號(hào)時(shí)域波形');axis(0 200 -0.3 0.3);subplot(2,1,2); psd(wav_noise, 10000, Fs);title(

4、'加噪聲后的語(yǔ)音信號(hào)功率譜密度');axis(0 25000 -20 10);figure(3); subplot(2,1,1); plot(wav_after(53550:53750); title('濾波后的語(yǔ)音信號(hào)時(shí)域波形');axis(0 200 -0.3 0.3);subplot(2,1,2); psd(wav_after, 10000, Fs);title('濾波后的語(yǔ)音信號(hào)功率譜密度');axis(0 25000 -20 10);sound(wav/max(wav), Fs); % 播放解調(diào)音頻sound(wav_after/max(wav_after), Fs); % 播放解調(diào)音頻w