matlab語音信號(hào)采集課程設(shè)計(jì)-文庫(kù)吧
2025-01-01 05:28 本頁面
【正文】 程序?yàn)榻忉寛?zhí)行，所以速度較慢。 （8）功能強(qiáng)大的工具箱是MATLAB的另一特色。MATLAB包含兩個(gè)部分：核心部分和各種可選的工具箱。核心部分中有數(shù)百個(gè)核心內(nèi)部函數(shù)。其工具箱又分為兩類：功能性工具箱和學(xué)科性工具箱。功能性工具箱主要用來擴(kuò)充其符號(hào)計(jì)算功能，圖示建模仿真功能，文字處理功能以及與硬件實(shí)時(shí)交互功能。功能性工具箱用于多種學(xué)科。而學(xué)科性工具箱是專業(yè)性比較強(qiáng)的。（9）源程序的開放性。開放性也許是MATLAB最受人們歡迎的特點(diǎn)。除內(nèi)部函數(shù)以外，所有MATLAB的核心文件和工具箱文件都是可讀可改的源文件，用戶可通過對(duì)源文件的修改以及加入自己的文件構(gòu)成新的工具箱。三．課程設(shè)計(jì)任務(wù) 利用Windows下的錄音機(jī)，錄制一段自己的話音，時(shí)間在3秒左右，然后在Matlab軟件平臺(tái)下，利用函數(shù)wavread對(duì)語音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)，采樣率8000Hz，單通道，格式*.wav。注：如果沒有辦法錄制自己語音，也可以用現(xiàn)成的語音也可以，文件格式 *.wav 在Matlab中，可以利用函數(shù)fft對(duì)信號(hào)進(jìn)行快速傅立葉變換，得到信號(hào)的頻譜特性，要求學(xué)生首先畫出語音信號(hào)的時(shí)域波形，然后對(duì)語音信號(hào)進(jìn)行頻譜分析。； 給定濾波器的性能指標(biāo)如下： (1)低通濾波器的性能指標(biāo)：fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB. (2)高通濾波器的性能指標(biāo)：fb=2700Hz, fc=3000Hz,As=100dB,Ap=1dB. (3)帶通濾波器的性能指標(biāo)： fb1=1200Hz, fb2=3000Hz, As=100dB,Ap=1dB. fc1=1000Hz, fc2=3200Hz,As=100dB, Ap=1dB. 采用窗函數(shù)法和雙線性變換法設(shè)計(jì)上面要求的3種濾波器，并畫出濾波器的頻率響應(yīng)； 然后用自己設(shè)計(jì)的濾波器對(duì)采集到的信號(hào)進(jìn)行濾波，畫出濾波后信號(hào)的時(shí)域波形及頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比，分析信號(hào)的變化； ，分析濾波前后的語音變化； 四．濾波器設(shè)計(jì)原理及設(shè)計(jì)方案根據(jù)過渡帶寬及阻帶衰減要求，選擇窗函數(shù)的類型并估計(jì)窗口長(zhǎng)度N（或階數(shù)M=N1），窗函數(shù)類型可根據(jù)最小阻帶衰減As獨(dú)立選擇，因?yàn)榇翱陂L(zhǎng)度N對(duì)最小阻帶衰減As沒有影響，在確定窗函數(shù)類型以后，可根據(jù)過渡帶寬小于給定指標(biāo)確定所擬用的窗函數(shù)的窗口長(zhǎng)度N，設(shè)待求濾波器的過渡帶寬為Δw，它與窗口長(zhǎng)度N近似成反比，窗函數(shù)類型確定后，其計(jì)算公式也確定了，不過這些公式是近似的，得出的窗口長(zhǎng)度還要在計(jì)算中逐步修正，原則是在保證阻帶衰減滿足要求的情況下，盡量選擇較小的N，在N和窗函數(shù)類型確定后，即可調(diào)用MATLAB中的窗函數(shù)求出窗函數(shù)wd（n）。根據(jù)待求濾波器的理想頻率響應(yīng)求出理想單位脈沖響應(yīng)hd(n)，如果給出待求濾波器頻率應(yīng)為Hd，則理想的單位脈沖響應(yīng)可以用下面的傅里葉反變換式求出：在一般情況下，hd(n)是不能用封閉公式表示的，需要采用數(shù)值方法表示；從w=0到w=2π采樣N點(diǎn)，采用離散傅里葉反變換（IDFT）即可求出。用窗函數(shù)wd(n)將hd(n)截?cái)?，并進(jìn)行加權(quán)處理，得到如果要求線性相位特性， 則h(n)還必須滿足：根據(jù)上式中的正、 負(fù)號(hào)和長(zhǎng)度N的奇偶性又將線性相位FIR濾波器分成四類。 要根據(jù)所設(shè)計(jì)的濾波特性正確選擇其中一類。 例如， 要設(shè)計(jì)線性相位低通特性可選擇h(n)=h(N1n)一類，而不能選h(n)=h(N1n)一類。 驗(yàn)算技術(shù)指標(biāo)是否滿足要求，為了計(jì)算數(shù)字濾波器在頻域中的特性，可調(diào)用freqz子程序，如果不滿足要求，可根據(jù)具體情況，調(diào)整窗函數(shù)類型或長(zhǎng)度，直到滿足要求為止。脈沖響應(yīng)不變法的主要缺點(diǎn)是產(chǎn)生頻率響應(yīng)的混疊失真。這是因?yàn)閺腟平面到Ｚ平面是多值的映射關(guān)系所造成的。為了克服這一缺點(diǎn)，可以采用非線性頻率壓縮方法，將整個(gè)頻率軸上的頻率范圍壓縮到π/T～π/T之間，再用z=esT轉(zhuǎn)換到Z平面上。也就是說，第一步先將整個(gè)S平面壓縮映射到S1平面的π/T～π/T一條橫帶里；第二步再通過標(biāo)準(zhǔn)變換關(guān)系z(mì)=es1T將此橫帶變換到整個(gè)Z平面上去。這樣就使S平面與Z平面建立了一一對(duì)應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖1所示。圖1雙線性變換的映射關(guān)系為了將S平面的整個(gè)虛軸jΩ壓縮到S1平面jΩ1軸上的π/T到π/T段上，可以通過以下的正切變換實(shí)現(xiàn) （1）式中,T仍是采樣間隔。當(dāng)Ω1由π/T經(jīng)過0變化到π/T時(shí)，Ω由∞經(jīng)過0變化到+∞，也即映射了整個(gè)jΩ軸。將式（1）寫成將此關(guān)系解析延拓到整個(gè)S平面和S1平面，令jΩ=s，jΩ1=s1，則得再將S1平面通過以下標(biāo)準(zhǔn)變換關(guān)系映射到Z平面z=es1T從而得到S平面和Z平面的單值映射關(guān)系為： (2) (3)式（2）與式（3）是S平面與Z平面之間的單值映射關(guān)系，這種變換都是兩個(gè)線性函數(shù)之比，因此稱為雙線性變換式（1）與式（2）的雙線性變換符合映射變換應(yīng)滿足的兩點(diǎn)要求。首先,把z=ejω，可得 (4)即S平面的虛軸映射到Z平面的單位圓。其次，將s=σ+jΩ代入式（4），得因此由此看出，當(dāng)σ0時(shí)，|z|1；當(dāng)σ0時(shí)，|z|1。也就是說，S平面的左半平面映射到Z平面的單位圓內(nèi)，S平面的右半平面映射到Z平面的單位圓外，S平面的虛軸映射到Z平面的單位圓上。因此，穩(wěn)定的模擬濾波器經(jīng)雙線性變換后所得的數(shù)字濾波器也一定是穩(wěn)定的。雙線性變換法優(yōu)缺點(diǎn):雙線性變換法與脈沖響應(yīng)不變法相比，其主要的優(yōu)點(diǎn)是避免了頻率響應(yīng)的混疊現(xiàn)象。這是因?yàn)镾平面與Z平面是單值的一一對(duì)應(yīng)關(guān)系。S平面整個(gè)jΩ軸單值地對(duì)應(yīng)于Z平面單位圓一周，即頻率軸是單值變換關(guān)系。這個(gè)關(guān)系
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1