国产 无码 综合区,色欲AV无码国产永久播放,无码天堂亚洲国产AV,国产日韩欧美女同一区二区

【數(shù)字信號(hào)處理課程設(shè)計(jì)】基于MATLAB實(shí)現(xiàn)語(yǔ)音信號(hào)的采集與處理(偏重濾波)

這篇具有很好參考價(jià)值的文章主要介紹了【數(shù)字信號(hào)處理課程設(shè)計(jì)】基于MATLAB實(shí)現(xiàn)語(yǔ)音信號(hào)的采集與處理(偏重濾波)。希望對(duì)大家有所幫助。如果存在錯(cuò)誤或未考慮完全的地方,請(qǐng)大家不吝賜教,您也可以點(diǎn)擊"舉報(bào)違法"按鈕提交疑問(wèn)。

目錄

一、目標(biāo)與任務(wù)

二、原理介紹

2.1 錄音原理

2.2 濾波器的設(shè)計(jì)原理及設(shè)計(jì)方法

2.3 IIR 數(shù)字濾波器設(shè)計(jì)原理

2.4 雙線性變換法

三、GUI界面設(shè)計(jì)與實(shí)現(xiàn)

四、基于MATLAB仿真

4.1實(shí)驗(yàn)過(guò)程

4.2 結(jié)果分析

五、總結(jié)

5.1 函數(shù)用法總結(jié)

5.2?心得體會(huì)

六、參考文獻(xiàn)


這個(gè)項(xiàng)目在我的B站上有專(zhuān)門(mén)的視頻演示:

【數(shù)字信號(hào)處理課程設(shè)計(jì)】基于MATLAB實(shí)現(xiàn)語(yǔ)音信號(hào)的采集與處理(側(cè)重濾波)+【通信原理課程設(shè)計(jì)】基于MATLAB實(shí)現(xiàn)的信源編碼系統(tǒng)仿真 +【GUI設(shè)計(jì)】_嗶哩嗶哩_bilibili

一、目標(biāo)與任務(wù)

1、 語(yǔ)音信號(hào)的采集。利用 Windows 下的錄音機(jī),錄制一段自己的話音,時(shí)間在 1s 內(nèi),然后在 Matlab 軟件平臺(tái)下,利用函數(shù) wavread 對(duì)語(yǔ)音信號(hào)進(jìn)行采樣,記住采樣頻率和采樣點(diǎn)數(shù)。

2、語(yǔ)音信號(hào)的頻譜分析。在 Matlab 中,可以利用函數(shù) fft 對(duì)信號(hào)進(jìn)行快速傅立葉變換,得到信號(hào)的 頻譜特性,要求學(xué)生首先畫(huà)出語(yǔ)音信號(hào)的時(shí)域波形,然后對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。

3、設(shè)計(jì)數(shù)字濾波器和畫(huà)出其頻率響應(yīng)給出各濾波器的性能指標(biāo); 給定濾波器的性能指標(biāo)如下:

(1)低通濾波器的性能指標(biāo):fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB.

(2)高通濾波器的性能指標(biāo):fc=4800Hz,fb=5000Hz,As=100dB,Ap=1dB.

(3)帶通濾波器的性能指標(biāo): fb1=1200Hz, fb2=3000Hz,fc1=1000Hz, fc2=3200Hz,As=100dB,Ap=1dB. 采用窗函數(shù)法和雙線性變換法設(shè)計(jì)上面要求的 3 種濾波器,并畫(huà)出濾波 器的頻率響應(yīng);

4、用濾波器對(duì)信號(hào)進(jìn)行濾波 ,然后用自己設(shè)計(jì)的濾波器對(duì)采集到的信號(hào)進(jìn)行濾波,畫(huà)出濾波后信號(hào)的 時(shí)域波形及頻譜,并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比,分析信號(hào)的變化;

5、回放語(yǔ)音信號(hào),分析濾波前后的語(yǔ)音變化;設(shè)計(jì)系統(tǒng)界面 為了使編制的程序操作方便,設(shè)計(jì)處理系統(tǒng)的用戶(hù)界面,在所設(shè)計(jì)的系統(tǒng)界面上可以實(shí)現(xiàn)上述要求中的包括采集、分析、濾波等全部?jī)?nèi)容,并能夠選 擇濾波器的類(lèi)型,輸入濾波器的參數(shù)、顯示濾波器的頻率響應(yīng)等。

二、原理介紹

2.1 錄音原理

????????Matlab 中 有 兩 個(gè) 函 數(shù) 可 以 用 于 實(shí) 現(xiàn) 錄 音 功 能 , 一 個(gè) 是waverecord 函數(shù)調(diào)用 windows 音頻設(shè)備進(jìn)行錄音,主要用于實(shí)時(shí)語(yǔ)音采集。第二個(gè)是使用Matlab 中的 wavread 函數(shù)也可以進(jìn)行語(yǔ)音采集。不過(guò) wavread 函數(shù)只能讀取計(jì)算機(jī)中已經(jīng)有的“.wav” 文件,也就是說(shuō)必須事先使用 windows 的錄音機(jī)等設(shè)備把有關(guān)的語(yǔ)音資料錄制并存儲(chǔ)為“.wav” 文件。顯然 wavread 函數(shù)不能用于實(shí)時(shí)語(yǔ)音信號(hào)采集。另外語(yǔ)音采集的頻率和精度不能改變。 回放時(shí)也有兩個(gè)函數(shù)可以選用。 Sound 函數(shù)和 wavplay,只不過(guò)后者的功能沒(méi)有前者強(qiáng)大。不過(guò)在本設(shè)計(jì)中主要是對(duì)已有的信號(hào)進(jìn)行處理,所以選用wavread 錄放。同時(shí),由于信號(hào)的時(shí)域描述只能反映信號(hào)的幅值隨時(shí)間的變化情況, 除只有一個(gè)頻率分量的簡(jiǎn)諧波外般很難明確揭示信號(hào)的頻率組成和各頻率分量的大小, 因此, 我們要把它轉(zhuǎn)換成頻域來(lái)分析,它能夠提供比時(shí)域信號(hào)波形更直觀,豐富的信息。

2.2 濾波器的設(shè)計(jì)原理及設(shè)計(jì)方法

FIR 濾波器的設(shè)計(jì)問(wèn)題在于尋求一系統(tǒng)函數(shù)H(z) ,使其頻率響應(yīng)H(ejω) 逼近濾波器要求的理想頻率響應(yīng)Hd(ejω) 其對(duì)應(yīng)的單位脈沖響應(yīng)hd(n) 。

1、用窗函數(shù)設(shè)計(jì)FIR 濾波器的基本方法

設(shè)計(jì)思想:從時(shí)域從發(fā),設(shè)計(jì)h(n) ?逼近理想hd(n) 。設(shè)理想濾波器Hd(ejω) 的單位脈沖響應(yīng)為hd(n) 。以低通線性相位 FIR 數(shù)字濾波器為例。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

hd(n) 一般是無(wú)限長(zhǎng)的,且是非因果的,不能直接作為FIR 濾波器的單位脈沖響應(yīng)。要想得到一個(gè)因果的有限長(zhǎng)的濾波器 h(n) ,最直接的方法是截?cái)?em>h(n)=hd(n)w(n) ,即截取為有限長(zhǎng)因果序列,并用合適的窗函數(shù)進(jìn)行加權(quán)作為FIR 濾波器的單位脈沖響應(yīng)。按照線性相位濾波器的要求,h(n) 必須是偶對(duì)稱(chēng)的。對(duì)稱(chēng)中心必須等于濾波器的延時(shí)常數(shù),即

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

用矩形窗設(shè)計(jì)的FIR 低通濾波器,所設(shè)計(jì)濾波器的幅度函數(shù)在通帶和阻帶都呈現(xiàn)出振蕩現(xiàn)象,且最大波紋大約為幅度的 9%,這個(gè)現(xiàn)象稱(chēng)為吉布斯(Gibbs)效應(yīng)。為了消除吉布斯效應(yīng),一般采用其他類(lèi)型的窗函數(shù)。

2、典型的窗函數(shù)

(1)矩形窗(Rectangle Window)

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其頻率響應(yīng)和幅度響應(yīng)分別為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

(2)三角形窗(Bartlett Window)

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其頻率響應(yīng)為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

(3)漢寧(Hanning)窗,又稱(chēng)升余弦窗

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其頻率響應(yīng)和幅度響應(yīng)分別為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

(4)漢明(Hamming)窗,又稱(chēng)改進(jìn)的升余弦窗

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其幅度響應(yīng)為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

(5)布萊克曼(Blankman)窗,又稱(chēng)二階升余弦窗

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其幅度響應(yīng)為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

(6)凱澤(Kaiser)窗

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

其中:β是一個(gè)可選參數(shù),用來(lái)選擇主瓣寬度和旁瓣衰減之間的交換關(guān)系,一般說(shuō)來(lái),β 越大,過(guò)渡帶越寬,阻帶越小衰減也越大。數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab 是第一類(lèi)修正零階貝塞爾函數(shù)。

若阻帶最小衰減表示為As=-20log10?δ ,β 的確定可采用下述經(jīng)驗(yàn)公式:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

若濾波器通帶和阻帶波紋相等即δp=δs 時(shí),濾波器節(jié)數(shù)可通過(guò)

下式確定:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

式中:數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

3.利用窗函數(shù)設(shè)計(jì)FIR 濾波器的具體步驟如下:

(1)按允許的過(guò)渡帶寬度△ω 及阻帶衰減 AS 選擇合適的窗函數(shù),并估計(jì)節(jié)數(shù)N:其中A 由窗函數(shù)的類(lèi)型決定。

(2)由給定的濾波器的幅頻響應(yīng)參數(shù)求出理想的單位脈沖響應(yīng)hd(n)

(3)確定延時(shí)值。

(4)計(jì)算濾波器的單位取樣響應(yīng)h(n)

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

??????(5)驗(yàn)算技術(shù)指標(biāo)是否滿(mǎn)足要求。

2.3 IIR 數(shù)字濾波器設(shè)計(jì)原理

IIR 數(shù)字濾波器是一種離散時(shí)間系統(tǒng),其系統(tǒng)函數(shù)為

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

??????假設(shè) M≤ N, 當(dāng)M> N時(shí), 系統(tǒng)函數(shù)可以看作一個(gè) IIR 的子系統(tǒng)和一個(gè) (M-N)的FIR子系統(tǒng)的級(jí)聯(lián)。 IIR 數(shù)字濾波器的設(shè)計(jì)實(shí)際上是求解濾波器的系數(shù) ak ?和 bk ,它是數(shù)學(xué)上的一種逼近問(wèn)題, 即在規(guī)定意義上 (通常采用最小均方誤差準(zhǔn)則) 去逼近系統(tǒng)的特性。 如果在 S平面上去逼近, 就得到模擬濾波器; 如果在 z平面上去逼近,就得到數(shù)字濾波器。

設(shè)計(jì)通數(shù)字濾波器通常可以歸納為如圖所示的兩種常用方法。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

?圖 2-1 數(shù)字帶通濾波器設(shè)計(jì)的兩種方法

方法 1: 首先設(shè)計(jì)一個(gè)模擬原型低通濾波器,然后通過(guò)頻率變換成所需要的模擬高通濾波器,最后再使用沖激不變法或雙線性變換成相應(yīng)的數(shù)字高通濾波器。

方法 2: 先設(shè)計(jì)一個(gè)模擬原型低通濾波器,然后采用沖激響應(yīng)不變法或雙線性變換法將它轉(zhuǎn)換成數(shù)字原型低通濾波器, 最后通過(guò)頻率變換把數(shù)字原型低通濾波器變換成所需要的數(shù)字高通濾波器。

本課程設(shè)計(jì)采用第一種設(shè)計(jì)方法, 先構(gòu)造一個(gè)巴特沃斯模擬低通濾波器, 然后將模擬低通濾波器轉(zhuǎn)換成模擬高通濾波器, 最后利用雙線性變換將模擬高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器。

2.4 雙線性變換法

為了克服沖激響應(yīng)法可能產(chǎn)生的頻率響應(yīng)的混疊失真, 這是因?yàn)閺?S 平面到Z平面是多值的映射關(guān)系所造成的。 為了克服這一缺點(diǎn), 可以采用非線性頻率壓縮方法,將整個(gè)頻率軸上的頻率范圍壓縮到 - π / T~π / T 之間,再用 z=esT 轉(zhuǎn)換到 Z 平面上。也就是說(shuō),第一步先將整個(gè) S平面壓縮映射到 S1 平面的 - π / T~π/ T 一條橫帶里;第二步再通過(guò)標(biāo)準(zhǔn)變換關(guān)系 z=esT 將此橫帶變換到整個(gè) Z 平面上去。 這樣就使 S平面與 Z 平面建立了一一對(duì)應(yīng)的單值關(guān)系, 消除了多值變換性,

也就消除了頻譜混疊現(xiàn)象,映射關(guān)系如圖 2-2

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 2-2 雙線性變換的映射關(guān)系

為了將 S平面的整個(gè)虛軸 j Ω 壓縮到 S1平面 j Ω 1 軸上的 - π / T 到π / T段上,可以通過(guò)以下的正切變換實(shí)現(xiàn)

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

式中 , T 仍是采樣間隔。

當(dāng)Ω 1 由- π / T 經(jīng)過(guò) 0 變化到 π / T 時(shí), Ω 由- ∞經(jīng)過(guò) 0 變化到 +∞,也即映射了整個(gè) j Ω 軸。將上式寫(xiě)成

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

將此關(guān)系解析延拓到整個(gè) S 平面和 S1平面,令 j Ω =s, j Ω 1=s1,則得

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

再將 S1平面通過(guò)以下標(biāo)準(zhǔn)變換關(guān)系映射到 Z 平面z=esT 從而得到 S 平面和 Z 平面的單值映射關(guān)系為:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

這兩個(gè)關(guān)系式是 S 平面與 Z 平面之間的單值映射關(guān)系, 這種變換都是兩個(gè)線性函數(shù)之比,因此稱(chēng)為雙線性變換

首先 , 把z=eiw 可得

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

即 S 平面的虛軸映射到 Z 平面的單位圓。

其次,將 s=σ+jw ?代入,得

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

因此

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

由此看出,當(dāng) <0 時(shí), | z|<1 ;當(dāng) >0 時(shí), | z|>1 。也就是說(shuō), 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)?S平面與 Z平面是單值的一一對(duì)應(yīng)關(guān)系。S 平面整個(gè) j Ω 軸單值地對(duì)應(yīng)于 Z 平面單位圓一周,即頻率軸是單值變換關(guān)系。

這個(gè)關(guān)系重寫(xiě)如下:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

上式表明, S平面上 Ω 與 Z 平面的 ω 成非線性的正切關(guān)系,如圖 2-3 所示。

由圖 2-3 看出, 在零頻率附近, 模擬角頻率 Ω 與數(shù)字頻率 ω 之間的變換關(guān)系接近于線性關(guān)系;但當(dāng) Ω 進(jìn)一步增加時(shí), ω 增長(zhǎng)得越來(lái)越慢,最后當(dāng) Ω →∞時(shí),ω 終止在折疊頻率 ω =π 處,因而雙線性變換就不會(huì)出現(xiàn)由于高頻部分超過(guò)折疊頻率而混淆到低頻部分去的現(xiàn)象,從而消除了頻率混疊現(xiàn)象。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 2-3 雙線性變換法的頻率變換關(guān)系

三、GUI界面設(shè)計(jì)與實(shí)現(xiàn)

頁(yè)面設(shè)計(jì)如圖3-1所示:

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 3-1 GUI界面

四、基于MATLAB仿真

4.1實(shí)驗(yàn)過(guò)程

(1)加入噪聲。

????????向原始信號(hào)中加入正弦干擾噪聲或高斯噪聲,加入正弦干擾噪聲后,如圖4-1所示。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-1 加入正弦干擾噪聲前后

????????加入高斯干擾噪聲后,如圖4-2所示。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-2 加入高斯噪聲前后

(2)濾波器的頻譜特性

  • IIR線性濾波器

????????巴特沃斯濾波器低通幅頻相頻特性如圖4-3所示

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-3 巴特沃斯低通濾波器幅頻相頻特性

????????車(chē)比雪夫?yàn)V波器高通幅頻相頻特性如圖4-4所示

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-4 車(chē)比雪夫高通濾波器幅頻相頻特性

????????橢圓濾波器帶通幅頻相頻特性如圖4-5所示

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-5 橢圓型帶通濾波器幅頻相頻特性

  • ?FIR窗函數(shù)濾波器

?????????Blackman窗低通抽樣響應(yīng)和頻譜特性如圖4-6所示

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-6 Blackman窗低通抽樣響應(yīng)和頻譜

?????????Kaiser窗高通抽樣響應(yīng)和頻譜特性如圖4-7所示

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-7 kaiser窗低通抽樣響應(yīng)和頻譜

(3)濾波前后波形

????????低通濾波后,如圖4-8所示。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-8 低通濾波

????????高通濾波后,如圖4-9所示。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-9 高通濾波

????????帶通濾波后,如圖4-10所示。

數(shù)字信號(hào)處理課程設(shè)計(jì),課程設(shè)計(jì),數(shù)字信號(hào)處理,matlab

圖 4-10 帶通濾波

4.2 結(jié)果分析

1.通過(guò)本次設(shè)計(jì)發(fā)現(xiàn),當(dāng)抽樣率Fs大于15000Hz時(shí),濾波器的設(shè)計(jì)除低通外,高通和帶通都會(huì)出現(xiàn)異常,主要異常為矩陣奇異值,或者NAN、INF;而Fs小于10000Hz時(shí),低通和帶通設(shè)計(jì)正常,但高通設(shè)計(jì)會(huì)報(bào)錯(cuò)異常,而Fs取11025Hz時(shí),可以實(shí)現(xiàn)設(shè)計(jì)需求。

2.使用巴特沃斯濾波器,低通、高通、帶通濾波效果都不能實(shí)現(xiàn),而切比雪夫1型濾波器和橢圓濾波器除高通外,低通和帶通效果都很明顯。

3.布萊克曼窗和凱塞窗在阻帶衰落為100dB時(shí)設(shè)計(jì)出的濾波器性能幾乎一致,且除了高通外,低通和帶通對(duì)加噪后的音頻濾波效果都很明顯。

4.加正弦噪聲時(shí),由于本設(shè)計(jì)所加單頻主要位于100Hz處,所以通過(guò)修改適當(dāng)?shù)臑V波器參數(shù)指標(biāo),信噪比較大或較小時(shí),高通和帶通都能很好的濾除噪聲。

5.加高斯噪聲時(shí),無(wú)論采用哪種濾波器都不能完美的濾除噪聲,但由于高斯噪聲的隨機(jī)特性,會(huì)發(fā)現(xiàn)加高斯噪聲后頻譜會(huì)出現(xiàn)很多尖峰,所以采用低通濾波器,在信噪比較小,即噪聲干擾很大的時(shí)候,濾波效果會(huì)很明顯。

五、總結(jié)

5.1 函數(shù)用法總結(jié)

1.[n, Wnj= buttord(Wp, Ws,Rp,Rs)

其中Wp和ws分別是通帶和阻帶的拐角頻率(截止頻率),其取值范圍為0至1之間。當(dāng)其值為1時(shí)代表采樣頻率的一半。Rp和Rs分別是通帶和阻帶區(qū)的波紋系數(shù)。

2.chebyl函數(shù)的用法為: [b,a]-chebyl1(n,Rp, Wn,/ftype/)

在使用chebyl函數(shù)設(shè)計(jì)IR濾波器之前,可使用cheblord函數(shù)求出濾波器階數(shù)n和截止頻率Wn, cheblord函數(shù)可在給定濾波器性能的情況下,選擇契比雪夫1型濾波器的最小階和截止頻率Wn.

3.cheblord函數(shù)的用法為:[n, Wnj-cheblord(Wp, Ws,Rp,Rs)

其中Wp和Ws分別是通帶和阻帶的拐角頻率(截止頻率),其取值范圍為0至1之間。當(dāng)其值為1時(shí)代表采樣頻率的一半。Rp和Rs分別是通帶和阻帶區(qū)的波紋系數(shù)。

4.butter函數(shù)的用法為:[b,a]=butter(n,Wn,/ftype/)

其中n代表濾波器階數(shù), Wn代表濾波器的截止頻率,這兩個(gè)參數(shù)可使用buttord函數(shù)來(lái)確定。buttord函數(shù)可在給定濾波器性能的情況下,求出巴特沃斯濾波器的最小階數(shù)n,同時(shí)給出對(duì)應(yīng)的截止頻率Wn.

5,Butterworth模擬低通濾波器函數(shù)語(yǔ)句: [z,p.kj—buttap(n)

其中n為濾波器的階數(shù), z. p、k分別為零點(diǎn)、極點(diǎn)和增益

6. [bz, az]=impinvar (b, atC, FS t, Fp)

式中, b, a為模擬濾波器分子和分母多項(xiàng)式系數(shù)向量: Fs為采樣頻率(所濾波數(shù)據(jù)),單位Hz,缺省時(shí)為1Hz.Fp為預(yù)畸變頻率(Prewarped frequency) ,是一個(gè)“匹配”頻率,在該頻率上,頻率響應(yīng)在變換前后和模擬頻率可精確匹配

在MATLAB中,函數(shù)bilinear采用雙線性變換法實(shí)現(xiàn)模擬s域至數(shù)字z域的映射,直接用于模擬濾波器變換為數(shù)字濾波器。其調(diào)用方式為:

[2zd, pd, kd]-bil inear (a, p, k,Fs)

[numd, dend]-bilinear (num, den, Fs)

式中, z,p分別為模擬濾波器零點(diǎn)、極點(diǎn)列向量: k為模擬濾波器的增益:Fs為采樣頻率,單位Hz. zd, pd, kd為數(shù)字濾波器的零極點(diǎn)和增益.num, den分別為模擬濾波器傳遞函數(shù)分子和分母多項(xiàng)式系數(shù)向量,模擬濾波器傳遞函數(shù)具有下面的形式:

H(s)-mum(5) mum(1)s+ +mum (mn)s+ num (nn+1)

den (s) den()s"..+den(nd )s + den (nd +1)

nund和dend分別為數(shù)字濾波器傳遞函數(shù)分子和分母多項(xiàng)式系數(shù)向量

7,函數(shù)trez用于求數(shù)字逃波器的頻率響應(yīng),其調(diào)用格式為:on[h, w]-]fregz (b, a, nL,'whole']) ;或[h, f]=freqz (b, a, nt,' whole'],Fs); 式中, b,a為數(shù)字濾波器分子和分母多項(xiàng)式的系數(shù), n為復(fù)數(shù)頻率的響應(yīng)點(diǎn)數(shù),為整數(shù),最好為2的冪,缺省時(shí)為512: Fs為采樣頻率,單位Hz,如果給定該值,則f位置輸出為頻率Hz,若沒(méi)有給定,則按角頻率(Angular frequency)給定f的頻率矢量:, whole'表示返回的頻率 或 值包含2平面整個(gè)單位圓頻率矢量,即0—27:缺省時(shí),頻率f或 值包含z平面上半單位圓(0—x)之間等間距n個(gè)點(diǎn)頻率矢量。h為復(fù)頻率響應(yīng): w為n點(diǎn)頻率向量(單位rad); f為n點(diǎn)頻率向量(Hz),函數(shù)返回值缺省時(shí),繪制幅頻響應(yīng)和相頻響應(yīng)圖。

5.2?心得體會(huì)

通過(guò)這次課設(shè), 對(duì)設(shè)計(jì)數(shù)字濾波器的整個(gè)設(shè)計(jì)過(guò)程有了很好的掌握, 懂得了設(shè)計(jì)濾波器的基本方法。 對(duì)IIR濾波器和FIR濾波器的應(yīng)用場(chǎng)景和區(qū)別有了一定了解,熟練了對(duì)巴特沃斯、車(chē)比雪夫、橢圓形濾波器和幾種常見(jiàn)窗函數(shù)的使用。同時(shí)也熟悉了 MATLAB的環(huán)境, 鞏固了相關(guān)知識(shí)。 初步掌握了 MATLAB語(yǔ)言在數(shù)字信號(hào)處理中一些基本庫(kù)函數(shù)的調(diào)用和編寫(xiě)基本程序等應(yīng)用; 熟悉了濾波器設(shè)計(jì)的一般原理, 對(duì)濾波器有了一個(gè)感性的認(rèn)識(shí); 學(xué)會(huì)了數(shù)字高通濾波器設(shè)計(jì)的一般步驟;加深了對(duì)濾波器設(shè)計(jì)中產(chǎn)生誤差的原因以及雙線性變換法優(yōu)缺點(diǎn)的理解和認(rèn)識(shí)??傊?, 使理論聯(lián)系了實(shí)際, 鞏固并深化了對(duì)課本基本知識(shí)的認(rèn)識(shí)和理解, 使理論得以升華。

在做本次課程設(shè)計(jì)的過(guò)程中,我和隊(duì)友曾經(jīng)遇到過(guò)問(wèn)題,在解決問(wèn)題的過(guò)程中鍛煉了自己處理問(wèn)題的能力,彼此相互溝通,通過(guò)查閱書(shū)籍、搜索資料等方法使問(wèn)題得以解決。 培養(yǎng)了我們的耐心,也深深感受到自身知識(shí)的有限。 通過(guò)使所學(xué)應(yīng)用于實(shí)踐,提高了動(dòng)手能力。

六、參考文獻(xiàn)

[1]高西全.數(shù)字信號(hào)處理(第三版).西安電子科技大學(xué)出版社,2018

[2]?王一世,數(shù)字信號(hào)處理,北京理工大學(xué)出版社.?

[3]?吳大正.?信號(hào)與線性系統(tǒng)分析(第四版).?高等教育出版社,2005.8?

[4]?陳垚光.?精通MATLAB?GUI設(shè)計(jì)(第二版)電子工業(yè)出版社,?2011.1?

[5]?陳后金.?信號(hào)分析與處理實(shí)驗(yàn).?高等教育出版社,2006.8?

[6]?趙力.?語(yǔ)音信號(hào)處理.?機(jī)械工業(yè)出版社,2009.7文章來(lái)源地址http://www.zghlxwxcb.cn/news/detail-800948.html

到了這里,關(guān)于【數(shù)字信號(hào)處理課程設(shè)計(jì)】基于MATLAB實(shí)現(xiàn)語(yǔ)音信號(hào)的采集與處理(偏重濾波)的文章就介紹完了。如果您還想了解更多內(nèi)容,請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章,希望大家以后多多支持TOY模板網(wǎng)!

本文來(lái)自互聯(lián)網(wǎng)用戶(hù)投稿,該文觀點(diǎn)僅代表作者本人,不代表本站立場(chǎng)。本站僅提供信息存儲(chǔ)空間服務(wù),不擁有所有權(quán),不承擔(dān)相關(guān)法律責(zé)任。如若轉(zhuǎn)載,請(qǐng)注明出處: 如若內(nèi)容造成侵權(quán)/違法違規(guī)/事實(shí)不符,請(qǐng)點(diǎn)擊違法舉報(bào)進(jìn)行投訴反饋,一經(jīng)查實(shí),立即刪除!

領(lǐng)支付寶紅包贊助服務(wù)器費(fèi)用

相關(guān)文章

  • 畢業(yè)設(shè)計(jì)-基于MATLAB的含噪語(yǔ)音信號(hào)降噪處理系統(tǒng)

    畢業(yè)設(shè)計(jì)-基于MATLAB的含噪語(yǔ)音信號(hào)降噪處理系統(tǒng)

    目錄 前言 課題背景和意義 實(shí)現(xiàn)技術(shù)思路 一、設(shè)計(jì)思路 二、IIR 數(shù)字濾波器的設(shè)計(jì)原理 三、語(yǔ)音信號(hào)的采集、 含噪語(yǔ)音信號(hào)的分析處理 四、總結(jié) 實(shí)現(xiàn)效果圖樣例 最后 ? ? ??大四是整個(gè)大學(xué)期間最忙碌的時(shí)光,一邊要忙著備考或?qū)嵙?xí)為畢業(yè)后面臨的就業(yè)升學(xué)做準(zhǔn)備,一邊

    2024年02月07日
    瀏覽(726)
  • 《FPGA數(shù)字信號(hào)處理》基于FPGA的32點(diǎn)并行FFT/IFFT設(shè)計(jì)

    《FPGA數(shù)字信號(hào)處理》基于FPGA的32點(diǎn)并行FFT/IFFT設(shè)計(jì)

    本人FPGA小白,只關(guān)注FPGA具體功能實(shí)現(xiàn)。如以下描述有誤,望評(píng)論區(qū)指正!以下內(nèi)容,純手打,嚴(yán)禁未經(jīng)過(guò)同意擅自轉(zhuǎn)載,謝謝! 相比于基2算法,基4算法對(duì)本設(shè)計(jì)(32點(diǎn)FFT)運(yùn)算的復(fù)雜度并沒(méi)有顯著的降低。 基2按時(shí)間抽取的FFT蝶形運(yùn)算如下圖所示: 可以發(fā)現(xiàn)輸入端與輸出端

    2024年04月23日
    瀏覽(37)
  • matlab函數(shù)轉(zhuǎn)C++(數(shù)字信號(hào)處理)

    近期主要利用QT完成一個(gè)本科的通信教學(xué)軟件,其中涉及大量matlab轉(zhuǎn)C++的工作,本來(lái)是想利用matlab的Coder模塊進(jìn)行轉(zhuǎn)換的,本人小白不太會(huì)用,還是自己按著matlab內(nèi)置函數(shù)的代碼進(jìn)行轉(zhuǎn)換,函數(shù)寫(xiě)的比較笨,希望大家能夠多多指導(dǎo). 使用的是C++的armadillo矩陣庫(kù)進(jìn)行矩陣的運(yùn)算,

    2024年02月06日
    瀏覽(28)
  • 使用matlab進(jìn)行回聲處理(三重回聲)-數(shù)字信號(hào)處理課設(shè)

    使用matlab進(jìn)行回聲處理(三重回聲)-數(shù)字信號(hào)處理課設(shè)

    1.錄制一段聲音信號(hào),作為原音頻信號(hào),生成頻域與時(shí)域圖形,觀察與分析其時(shí)域 與頻域圖形。 2.對(duì)該音頻信號(hào)進(jìn)行時(shí)域處理,并且實(shí)現(xiàn)對(duì)該聲音信號(hào)添加第一層回聲。 3.對(duì)該音頻信號(hào)再繼續(xù)一次進(jìn)行時(shí)域處理,實(shí)現(xiàn)對(duì)該聲音信號(hào)添加第二層回聲。 4.將三段音頻信號(hào)進(jìn)行合成

    2024年02月03日
    瀏覽(33)
  • 【數(shù)字信號(hào)處理】帶通采樣定理及其MATLAB仿真

    【數(shù)字信號(hào)處理】帶通采樣定理及其MATLAB仿真

    按照奈奎斯特采樣定理(低通采樣),采樣頻率 f s f_{s} f s ? 要大于等于信號(hào)中最高頻率 f m a x f_{max} f ma x ? 的2倍,才可以保證采樣后的數(shù)字信號(hào)通過(guò)DAC轉(zhuǎn)換后,可以無(wú)失真的恢復(fù)為原信號(hào)。然而,如果信號(hào)的頻率分布在某一有限頻帶上,并且信號(hào)的最高頻率 f m a x f_{max} f

    2024年02月16日
    瀏覽(79)
  • 信號(hào)處理之FIR數(shù)字濾波器(Matlab仿真)

    信號(hào)處理之FIR數(shù)字濾波器(Matlab仿真)

    ? ? ? ? 數(shù)字濾波器的作用是濾除不感興趣的信號(hào),留下想要的信號(hào)。數(shù)字濾波器可分為無(wú)限脈沖響應(yīng)(IIR)數(shù)字濾波器、有限脈沖響應(yīng)(FIR)數(shù)字濾波器兩種,兩者各有優(yōu)缺點(diǎn),其中FIR數(shù)字濾波器因其具有良好的線性相位特性受到廣泛應(yīng)用,線性相位是指信號(hào)中各頻率成分的相對(duì)

    2024年02月03日
    瀏覽(29)
  • 數(shù)字信號(hào)處理實(shí)驗(yàn)---LSI系統(tǒng)的分析 Matlab代碼

    數(shù)字信號(hào)處理實(shí)驗(yàn)---LSI系統(tǒng)的分析 Matlab代碼

    1.試用Matlab計(jì)算其幅頻特性和相頻特性,并繪圖。 代碼: n = 0:10; %定義采樣點(diǎn)n w = [0:1:500]*2*pi/500; % [0,pi]軸被分成1002個(gè)點(diǎn) x1 = power(0.9*exp(1i*pi/3),n); %定義輸入序列 x2 = exp(-1i*n); %定義一個(gè)系統(tǒng)的沖激響應(yīng) x = zeros(1,length(w)); %定義空數(shù)組存儲(chǔ)系統(tǒng)的頻域響應(yīng) for i=1:length(x1) ??? x=x

    2024年01月15日
    瀏覽(32)
  • matlab數(shù)字信號(hào)處理實(shí)驗(yàn)(5)時(shí)域采樣與頻域采樣

    一、實(shí)驗(yàn)?zāi)康?1、理解時(shí)域采樣理論與頻域采樣理論; 2、掌握模擬信號(hào)采樣前后頻譜的變化,以及如何選擇采樣頻率才能使采樣后的信號(hào) 不丟失信息; 3、掌握頻率域采樣會(huì)引起時(shí)域周期化的原因,頻率域采樣定理及其對(duì)頻域采樣點(diǎn)數(shù) 選擇的指導(dǎo)作用; 4、對(duì)信號(hào)在某個(gè)表示

    2024年02月07日
    瀏覽(35)
  • 【數(shù)字化處理】仿生假體控制中肌電信號(hào)的數(shù)字化處理研究(Matlab代碼實(shí)現(xiàn))

    【數(shù)字化處理】仿生假體控制中肌電信號(hào)的數(shù)字化處理研究(Matlab代碼實(shí)現(xiàn))

    ????????? 歡迎來(lái)到本博客 ???????? ??博主優(yōu)勢(shì): ?????? 博客內(nèi)容盡量做到思維縝密,邏輯清晰,為了方便讀者。 ?? 座右銘: 行百里者,半于九十。 ?????? 本文目錄如下: ?????? 目錄 ??1 概述 ??2 運(yùn)行結(jié)果 ??3?參考文獻(xiàn) ??4 Matlab代碼及文獻(xiàn) 文獻(xiàn)

    2024年02月12日
    瀏覽(38)
  • 【老生談算法】基于matlab時(shí)域頻域處理的語(yǔ)音信號(hào)變聲處理系統(tǒng)設(shè)計(jì)與算法原理(論文+程序源碼+GUI圖形用戶(hù)界面)——變聲算法

    【老生談算法】基于matlab時(shí)域頻域處理的語(yǔ)音信號(hào)變聲處理系統(tǒng)設(shè)計(jì)與算法原理(論文+程序源碼+GUI圖形用戶(hù)界面)——變聲算法

    大家好,今天給大家介紹基于matlab的語(yǔ)音信號(hào)變聲處理系統(tǒng)設(shè)計(jì)與算法原理(論文+程序源碼)。 運(yùn)用matlab軟件實(shí)現(xiàn)對(duì)聲音的變聲處理,利用離散付里葉變換進(jìn)行頻譜分析;設(shè)計(jì)數(shù)字濾波器組;通過(guò)時(shí)域和頻域方法做出各種音效效果,實(shí)現(xiàn)變速(慢放、快放),變調(diào)(頻譜左

    2024年02月04日
    瀏覽(43)

覺(jué)得文章有用就打賞一下文章作者

支付寶掃一掃打賞

博客贊助

微信掃一掃打賞

請(qǐng)作者喝杯咖啡吧~博客贊助

支付寶掃一掃領(lǐng)取紅包,優(yōu)惠每天領(lǐng)

二維碼1

領(lǐng)取紅包

二維碼2

領(lǐng)紅包