【導(dǎo)讀】畫出濾波后信號的時(shí)域波形和頻譜，并對濾波前后的。計(jì)算機(jī)科學(xué)和信息科學(xué)發(fā)展，數(shù)字處理技術(shù)應(yīng)運(yùn)而生并得以快速發(fā)展。人們進(jìn)行信息溝通的主要方式之一，它具有直接、自然、方便等優(yōu)點(diǎn)。語音處理主要是對語音進(jìn)行機(jī)器處理，以達(dá)到傳輸、自。動(dòng)識別、機(jī)器理解等目的。包絡(luò)提取的算法實(shí)現(xiàn)研究，討論了在算法的DSP實(shí)現(xiàn)方法，Matlab語言是一種廣泛應(yīng)用于工程計(jì)算及數(shù)值分析領(lǐng)域的新型高級語言,Matlab功能強(qiáng)大、簡單易學(xué)、編程效率高。特別是Matlab還具有信號分析工具

　　

【正文】 g_Callback(hObject, eventdata, handles) % hObject handle to caiyang (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x1,fs,bits]=wavread(39。D:\39。)。 %讀取信號 x1=x1(:,1)。 %只取單聲道 sound(x1,fs)。 %對聲音的回放 axes()。 plot(x1)。 %做原始語音信號的時(shí)域圖形 title(39。原始語音信號的時(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 y1=fft(x1)。 %做 length(x1)點(diǎn)的 FFT y1=fftshift(y1)。 %平移，是頻率中心為 0 derta_fs = fs/length(x1)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2derta_fs],abs(y1))。%畫出原始語音信號的頻譜圖 title(39。原始語音信號的頻譜 39。,39。fontname39。,39。宋 體 39。)。 % Executes on button press in yanshi. function yanshi_Callback(hObject, eventdata, handles) % hObject handle to yanshi (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x1,fs,bits]=wavread(39。D:\39。)。 x1=x1(:,1)。 %只取單聲道 z=[zeros(2020,1)。x1]。 %對聲音延時(shí) 2020 sound(z,fs)。 武漢理工大學(xué)《 信息群處理課程設(shè)計(jì) 》課程設(shè)計(jì)說明書 20 axes()。 plot(z)。 title(39。延時(shí)后的信號時(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 %畫出延時(shí)后的信號時(shí)域圖 z=fft(x1)。 %做 length(x1)點(diǎn)的 FFT z=fftshift(z)。 %平移，是頻率中心為 0 derta_fs = fs/length(z)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2derta_fs],abs(z))。%畫出延時(shí)后的信號頻譜圖 title(39。延時(shí)后的信號頻譜圖 39。,39。fontname39。,39。宋體 39。)。 % Executes on button press in hunxiang. function hunxiang_Callback(hObject, eventdata, handles) % hObject handle to hunxiang (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 z=[zeros(2020,1)。x]。 %對語音信號進(jìn)行延時(shí) x=[x。zeros(2020,1)]。 %使語音信號與延時(shí)后信號同等長度 y1=x+z。 %信號的混響 sound(y1,fs)。 %回放混響后的信號 axes()。 plot(y1)。 title(39?；祉懙臅r(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 %混響時(shí)域圖 Y1=fft(y1)。 %做 length(y1)點(diǎn)的 FFT Y1=fftshift(Y1)。 %平移，是頻率中心為 0 derta_fs = fs/length(y1)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y1))。%畫出混響頻譜圖 title(39?；祉懙念l譜圖 39。,39。fontname39。,39。宋體 39。)。 武漢理工大學(xué)《 信息群處理課程設(shè)計(jì) 》課程設(shè)計(jì)說明書 20 % Executes on button press in dan. function dan_Callback(hObject, eventdata, handles) % hObject handle to dan (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 a=。 %a取小于等于 1 z=[zeros(2020,1)。x]。 %對 語音信號進(jìn)行延時(shí) x=[x。zeros(2020,1)]。 %使語音信號與延時(shí)后信號同等長度 y2=x+z*a。 %信號經(jīng)單回聲濾波 sound(y2,fs)。 %回放經(jīng)單回聲濾波器后的語音信號 axes()。 plot(y2)。title(39。單回聲濾波的信號時(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 Y2=fft(y2)。 %做 length(y2)點(diǎn)的 FFT Y2=fftshift(Y2)。 %平移，是頻率中心為 0 derta_fs = fs/length(y2)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y2))。%畫出 單回聲信號頻譜圖 title(39。單回聲濾波的信號頻譜圖 39。,39。fontname39。,39。宋體 39。)。 % Executes on button press in duo. function duo_Callback(hObject, eventdata, handles) % hObject handle to duo (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x,fs,bits]=wavread(39。D:\39。)。 x1=x(:,1)。 %變換為單聲道語音信號 delay_t=。 %回聲延遲時(shí)間 武漢理工大學(xué)《 信息群處理課程設(shè)計(jì) 》課程設(shè)計(jì)說明書 20 a=。 %衰減系數(shù) N=3。 %N階延時(shí) Bz=[1。zeros(fs*delay_t,1)*N。a^N]。 Az=[1。zeros(fs*delay_t,1)。a]。 y=filter(Bz, Az ,x1 )。 sound(y,fs,bits)。 %回放經(jīng)多回聲濾波器后的語音信號 axes()。 plot(y)。title(39。多回聲濾波 的信號時(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 Y2=fft(y)。 %做 length(y)點(diǎn)的 FFT Y2=fftshift(Y2)。 %平移，是頻率中心為 0 derta_fs = fs/length(y)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y2))。%畫出多回聲信號頻譜圖 title(39。多回聲濾波的信號頻譜圖 39。,39。fontname39。,39。宋體 39。)。 % Executes on button press in wu. function wu_Callback(hObject, eventdata, handles) % hObject handle to wu (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 a=。 %a取小于等于 1 Bz=[0,0,0,0,0,0,0,0,0,0,1]。 %分子的系數(shù) Az=[1,0,0,0,0,0,0,0,0,0,a]。 %分母的系數(shù) yy1=filter(Bz,Az,x)。 %FIR濾波器進(jìn)行濾波 sound(yy1,fs,bits)。 %回放經(jīng)無限回聲濾波器后的語音信號 axes()。 plot(yy1)。 %無限回聲信號時(shí)域波形 武漢理工大學(xué)《 信息群處理課程設(shè)計(jì) 》課程設(shè)計(jì)說明書 20 title(39。無限個(gè)回聲濾波的信號時(shí)域圖 39。,39。fontname39。,39。宋體 39。)。 YY1=fft(yy1)。 %經(jīng)無限回聲濾波器后的信號做 length(yy1)點(diǎn)的 FFT YY1=fftshift(YY1)。 %平移，是頻率中心為 0 derta_fs = fs/length(yy1)。 %設(shè)置頻譜的間隔，分辨率 axes()。 plot([fs/2:derta_fs:fs/2der