【正文】 前面程序運(yùn)行的結(jié)果，直接可以在Matlab的“Workspace”里面查找到設(shè)計的FIR濾波器的脈沖相應(yīng)序列h_bs，其值為： 系數(shù)計算好以后，將濾波器結(jié)構(gòu)用visio或其他繪圖軟件畫好，填好系數(shù)即可。 在MATLAB中，經(jīng)過sound(x,fs,bits)函數(shù)，對經(jīng)過Flattopwin窗設(shè)計的FIR濾波器之后的音樂信號進(jìn)行回放，可以聽出濾波之后的音樂信號跟原始音樂信號一樣清晰，完全濾除掉了噪聲的干擾，通過下面語句來進(jìn)行音樂信號回放比較： sound (y_fil,fs,bits)。=20*log10((mag+eps)/max(mag))(eps，mag分別為濾波前后的噪聲頻率)中對濾波前后噪聲頻率的要求。)。wavwrite(y,fs,bits,39。39。 %播放原始語音與單頻噪聲混疊之后的信號sound(y_fil,fs,bits)。sound(x,fs,bits)。)。title(39。幅度譜39。)。xlabel(39。subplot(3,2,6)。濾波后語音信號39。)。ylabel(39。時間(t)39。plot(t, y_fil)。axis([0,5000,0,300])。加干擾后的語音信號幅度譜39。)。ylabel(39。頻率(f)39。plot(f,Y)。)。title(39。幅度39。)。xlabel(39。subplot(3,2,3)。)。title(39。幅度譜39。)。xlabel(39。subplot(3,2,2)。原始語音信號39。)。ylabel(39。時間(t)39。plot(t,x)。)。 % 計算頻譜取前一半figure(4)suptitle(39。% 用設(shè)計好的濾波器對y進(jìn)行濾波Y_fil=fft(y_fil)。 設(shè)計的濾波器圖形 信號濾波處理濾波器設(shè)計完成后，在MATLAB平臺上用函數(shù)filter實(shí)現(xiàn)濾波。)。ylabel(39。n39。)。title(39。subplot(2,2,4)。)。ylabel(39。w/pi39。)。title(39。subplot(2,2,3)。幅度mag39。)。xlabel(39。濾波器幅度響應(yīng)圖39。plot(w,mag)。,1)。,39。,39。,39。,39。line([wsu,wsu],[350,50],39。LineWidth39。39。linestyle39。r39。color39。,1)。,39。,39。,39。,39。line([0,],[Rp,Rp],39。LineWidth39。39。linestyle39。r39。color39。)。ylabel(39。w/pi39。)。title(39。% 調(diào)用自編函數(shù)計算濾波器的頻率特性figure(3)subplot(2,2,1)。.*hd_bs。 % 產(chǎn)生M階的Flattopwin窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)ideal_lp(wcu,M)。 % 計算窗設(shè)計該濾波器時需要的階數(shù)n=0:M1。wsu=fsu/fs*2*pi。dw=df/fs*2*pi。% 計算上下邊帶中心頻率，和頻率間隔 wcd=fcd/fs*2*pi。fcu=(fpu+fsu)/2。As=100。fpu=3800。fsd=1950。在本次的課程設(shè)計中我所采用的就是利用Flattopwin窗函數(shù)來設(shè)計FIR濾波器。)。ylabel(39。f39。)。grid on%title(39。) 。ylabel(39。頻率（單位Hz）39。)。plot(f,Y)title(39。)。ylabel(39。t39。)。grid on%title(39。)。ylabel(39。時間(單位：s)39。)。title(39。subplot(2,2,3)。X39。)。xlabel(39。原始語音信號的頻域圖39。) 。ylabel(39。頻率（單位：Hz）39。)。axis([0,3000,0,300])title(39。subplot(2,2,2)。x39。)。xlabel(39。原始語音信號的時域圖39。axis([0,2,1,1])。幅度39。)。xlabel(39。原始語音信號39。 原始語音信號經(jīng)過加噪處理后，可以通過Matlab軟件的plot函數(shù)分別對原始語音信號和加噪之后的混疊信號進(jìn)行時域和頻域的圖形對比,獲得最為直觀的差異，具體程序如下：subplot(2,2,1)。 % 截取前半部分deltaf=fs/N。 % 對原始信號和加噪信號進(jìn)行fft變換,取幅度譜X=X(1:N/2)。 % 應(yīng)該可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x))。 y=x+*sin(fn*2*pi*t)。 % 單頻噪聲頻率，2000Hz t=0:1/fs:(N1)/fs。 % 按指定的采樣率和每樣本編碼位數(shù)回放N=length(x)。% % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數(shù)。d:\39。 采集完成后在信號中加入一個單頻噪聲，設(shè)計的任務(wù)即為從含噪信號中濾除單頻噪聲，還原原始信號。 是否達(dá)到去噪 效果？ 無效 結(jié)束 有效 錄制語音信號在Windows XP環(huán)境下，打開“開始——程序——附件——錄音機(jī)”，在錄音機(jī)程序的“文件——屬性——立即轉(zhuǎn)換”中將聲音選定屬性選擇為“，8位，單聲道，7kB/s”, XP系統(tǒng)錄音機(jī)設(shè)置錄制語音信號“大家好，我是徐康”，時間在1～2s左右，然后保存為“.wav”格式的音頻文件然后在Matlab軟件平臺下，利