實(shí)驗(yàn)二_散系統(tǒng)頻率響應(yīng)和零極點(diǎn)分布

1、電信類課程試驗(yàn)報(bào)告電信類課程試驗(yàn)報(bào)告學(xué)院：院：基礎(chǔ)信息工程系別：別：電子信息工程課程名稱：課程名稱：數(shù)字信號(hào)處理姓名：名：學(xué)號(hào)：號(hào)：日期：期：實(shí)驗(yàn)二實(shí)驗(yàn)二實(shí)驗(yàn)名稱：散系統(tǒng)頻率響應(yīng)和零極點(diǎn)分布實(shí)驗(yàn)名稱：散系統(tǒng)頻率響應(yīng)和零極點(diǎn)分布一、實(shí)驗(yàn)?zāi)康囊?、?shí)驗(yàn)?zāi)康耐ㄟ^(guò)MATLAB仿真簡(jiǎn)單的離散時(shí)間系統(tǒng)，研究其時(shí)域特性，加深對(duì)離散系通的脈沖響應(yīng)，頻率響應(yīng)分析和零，極點(diǎn)分布概念的理解。二、主要函數(shù)簡(jiǎn)介二、主要函數(shù)簡(jiǎn)介1.1.計(jì)算離散時(shí)間系統(tǒng)的頻率響應(yīng)函數(shù)

2、計(jì)算離散時(shí)間系統(tǒng)的頻率響應(yīng)函數(shù)[hf]=freqz(banfs)輸入?yún)?shù)：矢量a和b的含義impz函數(shù)中的定義相同，為式（24）中分母，分子的系數(shù)；fs是取樣頻率，以Hz為單位；n為0到fs2的頻率范圍內(nèi)選取的頻率點(diǎn)數(shù)。輸出參數(shù)：h是計(jì)算所得的頻率響應(yīng)值；f是在0到fs2頻率范圍內(nèi)的頻率值。2.2.系統(tǒng)的有理分式形式轉(zhuǎn)化成零極點(diǎn)增益形式的函數(shù)系統(tǒng)的有理分式形式轉(zhuǎn)化成零極點(diǎn)增益形式的函數(shù)[zpk]=tf2zp(ba)輸入?yún)?shù)：b和a與fr

3、eqz相同。輸出參數(shù)：z表示由系統(tǒng)的零點(diǎn)構(gòu)成的矢量，p表示由系統(tǒng)的極點(diǎn)構(gòu)成的矢，k表示系統(tǒng)的增益。3.3.計(jì)算兩個(gè)有限長(zhǎng)序列線性卷積的函數(shù)計(jì)算兩個(gè)有限長(zhǎng)序列線性卷積的函數(shù)y=conv(xh)輸入?yún)?shù)：x和h分別代表需要進(jìn)行線性卷積運(yùn)算的兩個(gè)序列。輸出參數(shù)：y表示x和h的線性卷積結(jié)果。對(duì)于一個(gè)離散時(shí)間系統(tǒng)，如果輸入信號(hào)為x(n)輸出信號(hào)為y(n)在已知此系統(tǒng)的差方方程的情況下，求解y(n)的方法有以下兩種：（1）采用filter函數(shù)（函數(shù)

4、介紹見實(shí)驗(yàn)1）：y=filter(bax)filter函數(shù)的輸出y的長(zhǎng)度與輸入x的長(zhǎng)度相同，假設(shè)此離散時(shí)間系統(tǒng)是一個(gè)無(wú)限長(zhǎng)單位脈沖響應(yīng)系統(tǒng)，但輸入信號(hào)x(n)的長(zhǎng)度y(n)也比較短，這時(shí)就不能較完整地反映出信號(hào)的變化情況，也不能很好的反映系統(tǒng)的特性。在這種情況下，可通過(guò)補(bǔ)零的方法來(lái)解決，即給輸入信號(hào)x(n)補(bǔ)較多的零值點(diǎn)，再利用filter函數(shù)求輸出y(n).（2）采用impz函數(shù)（函數(shù)介紹見實(shí)驗(yàn)1）和conv函數(shù)。首先利用impz函數(shù)

5、求出離散時(shí)間系統(tǒng)的單位脈沖響應(yīng)，然后利用y(n)=x(n)h(n)的關(guān)系調(diào)用conv函數(shù)求（2）若輸入序列）若輸入序列x(n)=δ(n)2δ(n1)3δ(n2)4δ(n3)5δ(n4)，編程求此系統(tǒng)輸出序列，編程求此系統(tǒng)輸出序列y(n)并畫出其波形并畫出其波形程序程序x=[12345000000000000000000000000]num=[0.50.1]den=[11.61.28]y=filter(numdenx)stem(y)xla