基于matlab對語音信號進(jìn)行頻譜分析及濾波-資料下載頁
2024-11-16 20:22本頁面
【導(dǎo)讀】數(shù)字信號處理是一門以算法為核心，理論和實(shí)踐性較強(qiáng)的學(xué)科。數(shù)字信號處理課程是在學(xué)習(xí)完數(shù)字信號處理的相關(guān)理論后，進(jìn)行的綜合性。下面對各步驟加以具體說明。利用Windows下的錄音機(jī)，錄制一段自己的話音，時(shí)間在1s內(nèi)。于Matlab軟件平臺下，利用函數(shù)[x,fs,bits]=wavread對語音信號進(jìn)行。采樣，記住采樣頻率和采樣點(diǎn)數(shù)。通過wavread函數(shù)的使用，理解采樣頻率、＝3200Hz，As＝100dB，Ap＝1dB。要求用窗函數(shù)法和雙線性變換法設(shè)計(jì)上面要求的3種濾波器。IIR濾波器；利用Matlab中的函數(shù)freqz畫出各濾波器的頻率響應(yīng)。要求在一個(gè)窗口同時(shí)畫出濾波前后的波形及頻譜。在Matlab中，調(diào)用函數(shù)sound對聲音進(jìn)行回放。bits)；可以感覺濾波前后的聲音有變化。為了使編制的程序操作方便，設(shè)計(jì)處理系統(tǒng)的用戶界面。頻率響應(yīng)，選擇音頻信號等。MATLAB是一款功能超強(qiáng)的數(shù)學(xué)軟件，應(yīng)用于各個(gè)行業(yè)。操作起來很人性化。其中所附帶的GUI面向?qū)ο蟮挠押镁幊谭绞礁阌趯W(xué)
【正文】 圖及濾波圖： 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 19 3. FIR 帶 通濾波器頻率相位特性響應(yīng)圖及濾波圖： 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 20 4． IIR 低通濾波器頻率相位特性響應(yīng)圖及濾波圖： 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 21 5． IIR 高通濾波器頻率相位特性響應(yīng)圖及濾波圖： 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 22 6． IIR 帶通濾波器頻率相位特性響應(yīng)圖及濾波圖： 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 23 6． 2 結(jié)果分析： （詳見設(shè)計(jì)流程）。 本次 數(shù)字信號處理課程設(shè)計(jì)， GUI 只是一個(gè)工具，一個(gè)托體，通過 GUI 面向?qū)ο蟮娜诵曰喕诵盘柼幚淼某橄笮?；便于理解。?MATLAB 函數(shù)的調(diào)用設(shè)計(jì)濾波器以及對讀入音頻處理濾波，得到其時(shí) 域，頻域圖。 七、 設(shè)計(jì) 經(jīng)驗(yàn)總結(jié) Malta 無疑是一款功能超強(qiáng)的數(shù)學(xué)軟件，基于它可實(shí)現(xiàn)多領(lǐng)域的科學(xué)研究和開發(fā)。在本次課程設(shè)計(jì)中 ,收獲最多的是如何看待問題以及如何解決問題。一個(gè)很現(xiàn)實(shí)的問題就是搞科研開發(fā)，遇到問題是必不可免，關(guān)鍵是如何看待這樣的問題，持怎樣的心態(tài)，持怎樣的目的。再就是問題是遇到了，要如何解決以及用怎樣最優(yōu)的方案去處理，都得下心思去思考和摸索。再一個(gè)很現(xiàn)實(shí)的問題是，要學(xué)會學(xué)習(xí)，學(xué)會自學(xué)。畢竟社會日新月異，科技蓬勃發(fā)展，新領(lǐng)域新知識不斷涌現(xiàn)， 怎么樣才能在社會保持有效的立足，怎樣才能翻弄潮流，面對社會 變遷做到運(yùn)籌帷幄之中呢？在我看來，關(guān)鍵要保持良好心態(tài)，保持一種積極永不妥協(xié)永不放棄的信念，不斷進(jìn)取不斷開拓。 在學(xué)好基本課內(nèi)的專業(yè)知識之外，應(yīng)盡量多補(bǔ)充課外的和本專業(yè)相協(xié)調(diào)的專業(yè)知識。數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 24 強(qiáng)勢未來就業(yè)利劍。 參考文獻(xiàn)： 【 1】 MATLAB 與科學(xué)計(jì)算 【 2】 精通基于 MATLAB 的 GUI 設(shè)計(jì) 【 3】 MATLAB 在電子類學(xué)科的應(yīng)用 參考網(wǎng)址 ： 【 1】 【 2】 【 3】 附錄 ： 方案主體代碼： global y1。 global fs。 global fe。 global bits。 [fname,dirpath]=uigetfile(39。*.wav39。)。 myfile=[dirpath fname] [y1,fe,bits] = wavread(myfile)。 [H,f]=freqz(b,1,512,fs)。 % axes()。 plot(f,20*log10(abs(H)))。 xlabel(39。頻率 /Hz39。)。 ylabel(39。幅值 39。)。 grid on。 title(39。濾波器頻率響應(yīng)曲線 39。)。 fb=str2num(get(,39。string39。))。 fc=str2num(get(,39。string39。))。 FB=str2num(get(,39。string39。))。 FC=str2num(get(,39。string39。))。 As=str2num(get(,39。string39。))。 Ap=str2num(get(,39。string39。))。 fs=str2num(get(,39。string39。))。 temp=get(,39。value39。)。 num=get(,39。value39。)。 if num==2。 %FIR 數(shù)字信號處理 課程設(shè)計(jì) 報(bào)告 25 switch (temp) case 2 [n,wn,beta]=kaiserord([fb fc],[1 0],[ ],fs)。 b=fir1(n,wn,kaiser(n+1,beta))。 z=fftfilt(b,y1)。 case 3 [N,Wn,beta,typ] = kaiserord( [fc fb], [1 0], [ ],fs )。 b = fir1(N, Wn,39。high39。)。 z=fftfilt(b,y1)。 case 4 [N,Wn,beta,ftype]=kaiserord([fc fb FB FC],[0 1 0],[ ],2*fs)。 W1=fb/(fs/2)。 W2=FB/(fs/2)。 b=fir1(N,[W1 W2],kaiser(N+1,beta))。 z=fftfilt(b,y1)。 end elseif num==3 % switch temp case 2 wc=2*fc/fs。 wb=2*fb/fs。 [n,wn]=ellipord(wb,wc,Ap,As)。 [b,a]=ellip(n,Ap,As,wn)。 z=filter(b,a,y1)。 case 3 Wb=fb/(fs/2)。 Wc=fc/(fs/2)。 [n,Wn]=ellipord(Wb,Wc,Ap,As)。 [b,a]=ellip(n,Ap,As,Wn,39。high39。)。 % z=filter(b,a,y1)。 case 4 wc=[2*fc/fs,2*FC/fs]。 wb=[2*fb/fs,2*FB/fs]。 [n,wn]=buttord(wb,wc,Ap,As)。 [b,a]=butter(n,wn)。 Z=filter(b,a,y1)。 end en
法律信息相關(guān)推薦
鄂ICP備17016276號-1