離散系統(tǒng)分析

1、YUNNAN NORMAL UNIVERSITY本科學(xué)生實(shí)驗(yàn)報(bào)告學(xué)號姓名學(xué)院物電學(xué)院專業(yè)、班級1實(shí)驗(yàn)課程名稱教師及職稱開課學(xué)期2014 至 2015 學(xué)年 下 學(xué)期填報(bào)時(shí)間2015 年 5 月 14 日云南師范大學(xué)教務(wù)處編印實(shí)驗(yàn)序號 實(shí)驗(yàn)時(shí)間 實(shí)驗(yàn)預(yù)習(xí)實(shí)驗(yàn)名稱2015,5,13離散系統(tǒng)分析同析3-3121 .實(shí)驗(yàn)?zāi)康纳羁汤斫怆x散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)在分析離散系統(tǒng)的時(shí)域特性、頻域特性以及穩(wěn)定性中 的重要作用及意義，熟練掌握利用MATLAB分析離散系統(tǒng)的時(shí)域響應(yīng)、頻響特性和零極點(diǎn)的方法。掌握利用 DTFT和DFT確定系統(tǒng)特性的原理和方法。2 .實(shí)驗(yàn)原理、實(shí)驗(yàn)流程或裝置示意圖MATLAB提供了許多可用

2、于分析線性時(shí)不變連續(xù)系統(tǒng)的函數(shù)，主要包含有系統(tǒng)函數(shù)、系統(tǒng)時(shí)域響應(yīng)、系統(tǒng)頻域響應(yīng)等分析函數(shù)。1 .離散系統(tǒng)的時(shí)域響應(yīng)離散時(shí)間LTI系統(tǒng)可用如下的線性常系數(shù)差分方程來描述：yk aiyk 1aNiyk (N 1) aNyk Nboxk b1xk 1bM 1xk (M 1) bM xk M已知輸入信號xk以及系統(tǒng)初始狀態(tài) y-1 , y-2，就可以求出系統(tǒng)的響應(yīng)。MATLAB提供了計(jì)算差分方程的數(shù)值解的函數(shù)，可以計(jì)算上述 差分方程描述的離散系統(tǒng)的單位脈沖響應(yīng)、零狀態(tài)響應(yīng)和完全響應(yīng)等。在調(diào)用MATLAB函數(shù)時(shí)，需要利用描述該離散系統(tǒng)的系數(shù)函數(shù)。對差分方程進(jìn)行：變換即可得系統(tǒng)函數(shù)：1Y(z)均型H (Z

3、) 1X(z) 1 a1z.(M 1), MbM 1zbMz b(z)(N 1)NaNz ) aNza(z)在MATLAB中可使用向量a和向量b分別保存分母多項(xiàng)式和分子多項(xiàng)式的系數(shù):a 1, a，aN 1, aN b b0, b1, bM 1, bM 這些系數(shù)均從z0按z的降哥排列。離散系統(tǒng)的單位沖激響應(yīng)hk的計(jì)算h,k = impz(b, a):計(jì)算系統(tǒng)的單位脈沖響應(yīng)hk和相應(yīng)的時(shí)間向量 k；也可簡寫為:h = impz(b, a)。 其中:a 1, a1, aN 1, aN b b0,b1, bM 1,bM h,k = impz(b, a, n):計(jì)算n點(diǎn)單位脈沖響應(yīng) hk; 也可簡寫為：

4、h = impz(b, a, n)。impz(b, a):繪制單位脈沖響應(yīng)hk的圖形。離散系統(tǒng)響應(yīng)yk的計(jì)算y = filter(b, a, x)：計(jì)算系統(tǒng)在輸入 x作用下的零狀態(tài)響應(yīng)yk;a 1, a1, aN 1, aN b 聞“，,bM 1,bMy = filter(b, a, x, zi)yk。計(jì)算系統(tǒng)在輸入x和初始狀態(tài)作用下的完全響應(yīng)zi是由系統(tǒng)的初始狀態(tài)經(jīng)過filtic函數(shù)轉(zhuǎn)換而得到的初始條件：zi= filtic (b, a, Y0) , Y0為系統(tǒng)的初始狀態(tài)，丫0= y-1 , y-2, y-3,。2 .離散系統(tǒng)的系統(tǒng)函數(shù)零極點(diǎn)分析離散LTI系統(tǒng)的系統(tǒng)函數(shù) H(z)可以表示為零

5、極點(diǎn)形式:H(z)組：z1X(z) 1 a1zb z (M 1) b z MbM 1zbM z(N 1)NaN iz3nzk (z z)(z z2).(z zm )(z P1)(z P2)(z Pn)使用MATLAB提供的roots函數(shù)計(jì)算離散系統(tǒng)的零極點(diǎn)；使用zplane函數(shù)繪制離散系統(tǒng)的零極點(diǎn)分布圖。注意：在利用這些函數(shù)時(shí)， 要求H(z)的分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)的個(gè)數(shù)相等，若不等則需要補(bǔ)零。3 .離散系統(tǒng)的頻率響應(yīng)當(dāng)離散因果LTI系統(tǒng)的系統(tǒng)函數(shù)H(z)的極點(diǎn)全部位于z平面單位圓內(nèi)時(shí)，系統(tǒng)的頻率響應(yīng)可 由H求出,即H(ej ) H(z)L|H(ej|ej()H, w=freqz(b,

6、a, n)計(jì)算系統(tǒng)的n點(diǎn)頻率響應(yīng)H , w為頻率點(diǎn)向量。H=freqz(b, a, w):計(jì)算系統(tǒng)在指定頻率點(diǎn)向量w上的頻響；freqz(b,a): 繪制頻率響應(yīng)曲線。其中：b和a分別為系統(tǒng)函數(shù) H(z)的分子分母系數(shù)矩陣；4 .利用DTFT和DFT確定離散系統(tǒng)的特性 _在很多情況下，需要根據(jù)LTI系統(tǒng)的輸入和輸出對系統(tǒng)進(jìn)行辨識，即通過測量系統(tǒng)在已知輸入xk激勵(lì)下的響應(yīng)yk來確定系統(tǒng)的特性。若系統(tǒng)的脈沖響應(yīng)為hk,由于存在yk=xk* hk,因而可在時(shí)域通過解卷積方法求解hk。在實(shí)際應(yīng)用中，進(jìn)行信號解卷積比較困難。因此，通常從頻域來分析系統(tǒng)，這樣就可以將時(shí)域的卷積轉(zhuǎn)變?yōu)轭l域的乘積，從而通過分析

7、系統(tǒng)輸入序列和輸出序列的頻率特性來確定系統(tǒng)的頻率特性H(ej ),再由H(ej )得到系統(tǒng)的脈沖響應(yīng)hk。若該LTI系統(tǒng)輸入xk的DTFT為X(ej )，系統(tǒng)輸出yk的DTFT為丫(ej ),則系統(tǒng)的Y(ej )X(ej )H(ej )頻率特性可表不為:N個(gè)等間隔頻率點(diǎn)上的樣點(diǎn)值:有限長序列的DTFT可以利用FFT計(jì)算出其在區(qū)間內(nèi)的X = fft(x,N);Y = fft(y,N);再利用 H = Y./X和h = ifft(H,N)可以得到系統(tǒng)的單位脈沖響應(yīng)hk3 .實(shí)驗(yàn)設(shè)備及材料計(jì)算機(jī)一臺及MATLAB仿真軟件。4 .實(shí)驗(yàn)方法步驟及注意事項(xiàng)實(shí)驗(yàn)方法步驟：先打開電腦，然后再打開 MATLAB

8、仿真軟件，在Blank M-File中輸入程序，然后再編譯運(yùn) 行程序，直到程序能編譯運(yùn)行為止。(1) .為了省時(shí)間以及編譯的方便性，程序應(yīng)該在 Blank M-File中輸入，而不應(yīng)該在 CommandWindow中直接運(yùn)行；(2) .在使用MATLAB時(shí)應(yīng)注意中英輸入法的切換，在中文輸入法輸入程序時(shí)得到的程序是錯(cuò) 誤的；(3) . MATLAB中兩個(gè)信號相乘表示為 x.*u,中間有個(gè).，同樣兩個(gè)信號相除也是如此，也就是在實(shí)驗(yàn)中要注意乘和點(diǎn)乘的區(qū)別。二.實(shí)驗(yàn)內(nèi)容1 .實(shí)驗(yàn)現(xiàn)象與結(jié)果1.1. 知某LTI系統(tǒng)的差分方程為：yk 1.143yk 1 0.412yk 20.0675xk 0.1349x

9、k 1 0.0675xk 2(1)初始狀態(tài)y 1 1, y 2 2 ,輸入 xk uk計(jì)算系統(tǒng)的完全響應(yīng)。(2)當(dāng)以下三個(gè)信號分別通過系統(tǒng)時(shí)，分別計(jì)算離散系統(tǒng)的零狀態(tài)響應(yīng)：x1k cos(k)uk; x2k cos(k)uk; x3k cos(7 k)uk10510(3)該系統(tǒng)具有什么特性？(1) a=1,-1.143,0.412;b=0.0675,0.1349,0.0675;N=100;x=ones(1,N);zi=filtic(b,a,1,2); y=filter(b,a,x,zi) stem(y);(2) a=1,-1.143,0.412;b=0.0675,0.1349,0.0675;N

10、=100;k=1:N; x1=cos(pi/10*k); y1=filter(b,a,x1) stem(y1);x2=cos(pi/5*k); y2=filter(b,a,x2) stem(y2);x3=cos(7*pi/10*k);y3=filter(b,a,x3)stem(y3);2.已知某因果LTI離散系統(tǒng)的系統(tǒng)函數(shù)為0.03571 0.1428* zA 1 0.2143* zA 2 0.1428* zA 3 0.03571* zA 4H(z)=1 1.035* zA 1 0.8264 *zA 2 0.2605* zA 3 0.04033* zA 4(1)計(jì)算系統(tǒng)的單位脈沖響應(yīng)。(2)當(dāng)

11、信號xk=uk+cos(兀k/4)*uk+ cos(兀k/2)*uk通過系統(tǒng)時(shí)，計(jì)算系統(tǒng)的零狀態(tài)響應(yīng)。(1) clc,clear,closeallN=40;a=1,-1.035,0.826,-0.2605,0.04033;b=0.03571,0.1428,0.2143,0.1428,0.03571;y=impz(b,a,N);(2) clc,clear,closeallN=40;k=0:0,01:100;a=1,-1,035,0.826,-0.2605,0.04033;b=0,03571,0.1428,0.2143,0.1428,0.03571;x1=1+cos(pi*k/4)+ cos(pi

12、*k/2);x=ones(x1,N);y=filter(b,a,x)y =Columns 1 through 80.0357 0.0357 0.03570.0357 0.03570.2155 0.2155 0.21550.2155 0.21550.5863 0.5863 0.58630.5863 0.5863Columns 9 through 160.0357 0.0357 0.03570.0357 0.03570.2155 0.2155 0.21550.2155 0.21550.5863 0.5863 0.58630.5863 0.5863Columns 17 through 240.035

13、7 0.0357 0.03570.0357 0.03570.2155 0.2155 0.21550.2155 0.21550.5863 0.5863 0.58630.5863 0.5863Columns 25 through 320.0357 0.0357 0.03570.0357 0.03570.2155 0.2155 0.21550.2155 0.21550.0357 0.0357 0.03570.2155 0.2155 0.21550.5863 0.5863 0.58630.0357 0.0357 0.03570.2155 0.2155 0.21550.5863 0.5863 0.586

14、30.0357 0.0357 0.03570.2155 0.2155 0.21550.5863 0.5863 0.58630.0357 0.0357 0.03570.2155 0.2155 0.21550.5863 0.58630.5863 0.58630.58630.5863 0.5863 0.5863Columns 33 through 400.0357 0.0357 0.03570.0357 0.03570.2155 0.2155 0.21550.2155 0.21550.5863 0.5863 0.58630.5863 0.58630.0357 0.0357 0.03570.2155

15、0.2155 0.21550.5863 0.5863 0.58632 .對實(shí)驗(yàn)現(xiàn)象、實(shí)驗(yàn)結(jié)果的分析及其結(jié)論思考題：離散時(shí)間系統(tǒng)的穩(wěn)定性與系統(tǒng)零點(diǎn)無關(guān)，與系統(tǒng)的極點(diǎn)有關(guān)，而系統(tǒng)零點(diǎn)則影響系統(tǒng)單位脈沖響應(yīng)的幅度和相位。理解系統(tǒng)的零極點(diǎn)與系統(tǒng)的穩(wěn)定性之間的關(guān)系有利于對系統(tǒng)的理解。如果給定系統(tǒng)函數(shù) H(s),或給定系統(tǒng)微分方程(可以求出系統(tǒng)函數(shù))，通過系統(tǒng)函數(shù)可以零 極點(diǎn)圖判斷系統(tǒng)的穩(wěn)定性。結(jié)論：當(dāng)離散LTI系統(tǒng)的系統(tǒng)函數(shù) H(Z)的ROC包含單位圓時(shí)，系統(tǒng)的頻率響應(yīng)是系統(tǒng)函數(shù)在單位圓上的Z變換，系統(tǒng)函數(shù)與系統(tǒng)的穩(wěn)定性表現(xiàn)在對于離散LTI系統(tǒng)，系統(tǒng)穩(wěn)定的充分必要條件是H(Z)的ROC&含Z平面上的單位圓。若離散 LTI系統(tǒng)是因果系統(tǒng)，由于 hk是因果 序列，所以穩(wěn)定系統(tǒng)的RO3、須包含單位圓及單位圓外Z平面的全部區(qū)域，由于 H(Z)在其ROM不能有極點(diǎn)，故因果的離散LTI系統(tǒng)穩(wěn)定的充要條件是 H(Z)的極點(diǎn)全部位于 Z平面的 單位圓內(nèi)。總結(jié)：通過本次試驗(yàn)可以得出需要注意的是，離散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)可能有兩種形式，一種是分子和分母多項(xiàng)式均按 Z的正募次項(xiàng)排列，另一種是分子和分母多項(xiàng)式均按Z的負(fù)哥次項(xiàng)排列。無論采用哪種排列