1、 信號與系統(tǒng) 分析實驗報告 實驗項目名稱：離散線性時不變系統(tǒng)分析； 連續(xù)時間系統(tǒng)分析所屬課程名稱：信號與系統(tǒng)實驗教程實 驗 類 型 ：驗證型指 導(dǎo) 教 師 ：實 驗 日 期 ：2013.06.04班 級 ：學(xué) 號 ：姓 名 ：離散線性時不變系統(tǒng)分析一、實驗?zāi)康?. 掌握離散線性時不變系統(tǒng)的單位序列響應(yīng)、單位階躍響應(yīng)和任意激勵下響應(yīng)的MATLAB求解方法。2.掌握離散線性時不變系統(tǒng)的頻域分析方法；3.掌握離散線性時不變系統(tǒng)的復(fù)頻域分析方法；4. 掌握離散線性時不變系統(tǒng)的零極點分布與系統(tǒng)特性的關(guān)系。二、實驗原理及方法1. 離散線性時不變系統(tǒng)的時域分析描述一個N階線性時不變離散時間系統(tǒng)的數(shù)學(xué)模型是線

2、性常系統(tǒng)差分方程，N階線性時不變離散系統(tǒng)的差分方程一般形式為 （2.1）也可用系統(tǒng)函數(shù)來表示 （2.2）系統(tǒng)函數(shù)反映了系統(tǒng)響應(yīng)和激勵間的關(guān)系。一旦上式中，的數(shù)據(jù)確定了，系統(tǒng)的性質(zhì)也就確定了。特別注意必須進(jìn)行歸一化處理，即。對于復(fù)雜信號激勵下的線性系統(tǒng)，可以將激勵信號在時域中分解為單位序列或單位階躍序列的線性疊加，把這些單元激勵信號分別加于系統(tǒng)求其響應(yīng)，然后把這些響應(yīng)疊加，即可得到復(fù)雜信號作用于系統(tǒng)的零狀態(tài)響應(yīng)。因此，求解系統(tǒng)的單位序列響應(yīng)和單位階躍響應(yīng)尤為重要。由圖2-1可以看出一個離散LSI系統(tǒng)響應(yīng)與激勵的關(guān)系。圖2-1 離散LSI系統(tǒng)響應(yīng)與激勵的關(guān)系(1) 單位序列響應(yīng)（單位響應(yīng)）單位響應(yīng)

3、是指離散線性時不變系統(tǒng)在單位序列激勵下的零狀態(tài)響應(yīng)，因此滿足線性常系數(shù)差分方程(2.1)及零初始狀態(tài)，即 ， (2.3)按照定義，它也可表示為 (2.4)對于離散線性時不變系統(tǒng)，若其輸入信號為，單位響應(yīng)為，則其零狀態(tài)響應(yīng)為 (2.5)可見，能夠刻畫和表征系統(tǒng)的固有特性，與何種激勵無關(guān)。一旦知道了系統(tǒng)的單位響應(yīng)，就可求得系統(tǒng)對任何輸入信號所產(chǎn)生的零狀態(tài)響應(yīng)。MATLAB提供了專門用于求離散系統(tǒng)沖激響應(yīng)的函數(shù)impz()，其調(diào)用格式有h,n=impz(b,a)求解離散系統(tǒng)的單位響應(yīng)，其中，；h,n=impz(b,a,N) 求解離散系統(tǒng)的單位響應(yīng)，采樣點數(shù)由N確定，；impz(b,a) ：在當(dāng)前窗口

4、，用stem(n,h)繪出圖形。（2）單位階躍響應(yīng) 單位階躍響應(yīng)是指離散離散時不變系統(tǒng)在單位階躍序列激勵下的零狀態(tài)響應(yīng)，它可以表示為 (2.6)上式表明，離散線性時不變系統(tǒng)的單位階躍響應(yīng)是單位響應(yīng)的累加和，系統(tǒng)的單位階躍響應(yīng)和系統(tǒng)的單位響應(yīng)之間有著確定的關(guān)系，因此，單位階躍響應(yīng)也能完全刻畫和表征一個線性時不變系統(tǒng)。MATLAB提供了專門用于求離散系統(tǒng)單位階躍響應(yīng)的函數(shù)stepz( )，其調(diào)用格式有s,n=stepz(b,a) ：求解離散系統(tǒng)的單位階躍響應(yīng)，其中，；s,n=stepz(b,a,N) ：求解離散系統(tǒng)的單位階躍響應(yīng)，采樣點數(shù)由N確定，；stepz(b,a) ：在當(dāng)前窗口，用stem(

5、n,s)繪出圖形。（3）任意激勵下的零狀態(tài)響應(yīng)已經(jīng)知道，離散線性時不變系統(tǒng)可用常系數(shù)線性差分方程(2.1)式來描述，Matlab提供的函數(shù)dlsim( )能對上述差分方程描述的離散線性時不變系統(tǒng)的響應(yīng)進(jìn)行仿真，該函數(shù)不僅能繪制指定時間范圍內(nèi)的系統(tǒng)響應(yīng)波形圖，而且還能求出系統(tǒng)響應(yīng)的數(shù)值解。其調(diào)用格式有dlsim(b,a, x) ：求解輸入序列為x的零狀態(tài)響應(yīng)需要特別強(qiáng)調(diào)的是，Matlab總是把由分子和分母多項式表示的任何系統(tǒng)都當(dāng)作是因果系統(tǒng)。所以，利用impz (b,a)，stepz(b,a)，dlsim(b,a,x)函數(shù)求得的響應(yīng)總是因果信號。同時，卷積和也是線性時不變系統(tǒng)求解零狀態(tài)響應(yīng)的重要

6、工具之一。假設(shè)系統(tǒng)的輸入信號為，單位響應(yīng)為，則系統(tǒng)的零狀態(tài)響應(yīng)可由（2.5）式求解。Matlab提供了專門用于求離散系統(tǒng)卷積和的函數(shù)conv( )，其調(diào)用格式有y=conv(x,h) ：求解序列x，h的卷積和，若序列x的長度為n1，序列h的長度為n2，卷積和y的長度為n1+n2-1。這一點需要特別注意，否則，作圖時容易造成橫縱坐標(biāo)長度不匹配。（4）帶初始狀態(tài)的任意激勵下的全響應(yīng)任意激勵下的離散線性時不變系統(tǒng)的全響應(yīng)為零輸入響應(yīng)和零狀態(tài)響應(yīng)之和，表示為 （2.7）Matlab提供了用于求離散系統(tǒng)全響應(yīng)的函數(shù)filter( )，其調(diào)用格式有y=filter( b,a,x) ：求解零狀態(tài)響應(yīng)；y=f

7、ilter( b,a,x,zi) ：求解初始條件為zi的系統(tǒng)的全響應(yīng)，zi向量的長度為max(length(a),length(b)-1，返回值為系統(tǒng)的全響應(yīng)。z = filtic(b,a,y,x)：將初始狀態(tài)轉(zhuǎn)換為初始條件，其中，；z = filtic(b,a,):將初始狀態(tài)轉(zhuǎn)換為初始條件,2 離散線性時不變系統(tǒng)的復(fù)頻域（Z域）分析（1）利用Z變換解差分方程在前面圖2-1中表示了離散系統(tǒng)的響應(yīng)與激勵的關(guān)系，由圖可知，系統(tǒng)的響應(yīng)既可以用時域的方法求解，也可以用Z域的方法求解。當(dāng)已知系統(tǒng)輸入序列的Z變換，系統(tǒng)函數(shù)時，系統(tǒng)響應(yīng)序列的Z變換可由求出。Matlab提供了用于求序列Z變換和Z反變換的函數(shù)

8、，其調(diào)用格式有X=ztrans(x)：求無限長序列x的Z變換，返回Z變換的表達(dá)式，注意這里x，X都是符號表達(dá)式；x=iztrans(X)：求X（z）的Z反變換x(n)，返回Z反變換的表達(dá)式，注意這里x，X都是符號表達(dá)式；r,p,c=residuez(b,a)：把b(z)/a(z)展開成部分分式；b,a=residuez(r,p,c):根據(jù)部分分式的r、p、c數(shù)組，返回有理多項式。（2）系統(tǒng)的零極點分布與系統(tǒng)因果性和穩(wěn)定性的關(guān)系因果系統(tǒng)的單位響應(yīng)一定滿足當(dāng)時，=0，那么其系統(tǒng)函數(shù)的收斂域一定包含點，即點不是極點，極點分布在某個圓的圓內(nèi)，收斂域是圓外區(qū)域。系統(tǒng)穩(wěn)定要求 ，對照z變換定義，系統(tǒng)穩(wěn)定要

9、求收斂域包含單位圓。如果系統(tǒng)因果且穩(wěn)定，收斂域包含點和單位圓，那么收斂域可表示為： (2.8)MATLAB提供了用于求系統(tǒng)零極點的函數(shù)，其調(diào)用格式有roots()：利用多項式求根函數(shù)來確定系統(tǒng)函數(shù)的零極點位置；roots(a)：求極點位置，a為系統(tǒng)函數(shù)分母多項式所構(gòu)成的系數(shù)向量；roots(b)：求零點位置，b為系統(tǒng)函數(shù)分子多項式所構(gòu)成的系數(shù)向量；zplane(b,a)：繪制由行向量b和a構(gòu)成的系統(tǒng)函數(shù)的零極點分布圖；zplane(z,p)：繪制由列向量z確定的零點、列向量p確定的極點構(gòu)成的零極點分布圖。（3）系統(tǒng)的零極點分布與系統(tǒng)頻率響應(yīng)的關(guān)系將式(2.2)因式分解，得到 (2.9)式中，是

10、的零點，是其極點。參數(shù)影響頻率響應(yīng)的幅度大小，影響系統(tǒng)特性的是零點和極點的分布。下面采用幾何方法研究系統(tǒng)零極點分布對系統(tǒng)頻率特性的影響。將式(2.9)的分子、分母同乘以，得到： (2.10)假設(shè)系統(tǒng)穩(wěn)定，將代入上式，得到頻率響應(yīng) (2.11)設(shè)，由式(2.11)得到 (2.12)在平面上，用一根由零點指向單位圓（）上任一點B的向量表示，同樣用一根由極點指向單位圓（）上任一點B的向量來表示，和分別稱為零點矢量和極點矢量，用極坐標(biāo)表示為： ，。將和表示式代入式(2.12)，得到 (2.13) (2.14)系統(tǒng)或者信號的頻率特性由式(2.13)和式(2.14)確定。按照式(2.13)，知道零極點的分

11、布后，可以很容易地確定零極點位置對系統(tǒng)特性的影響。當(dāng)點轉(zhuǎn)到極點附近時，極點矢量長度最短，因而幅度特性可能出現(xiàn)峰值，且極點愈靠近單位圓，極點矢量長度愈短，峰值愈高愈尖銳。如果極點在單位圓上，則幅度特性為，系統(tǒng)不穩(wěn)定。對于零點，情況相反，當(dāng)點轉(zhuǎn)到零點附近時，零點矢量長度變短，幅度特性將出現(xiàn)谷值，且零點愈靠近單位圓，谷值愈接近零。當(dāng)零點在單位圓上時，谷值為零。綜上所述，極點位置主要影響頻響的峰值位置及尖銳程度，零點位置主要影響頻響的谷值位置及形狀。Matlab提供了專門用于求離散系統(tǒng)頻率響應(yīng)的函數(shù)freqz(),其調(diào)用格式如下：H,w = freqz(b,a,n)：返回數(shù)字系統(tǒng)的n點頻率值（復(fù)數(shù)），

12、這n個點均勻地分布在0，上，系統(tǒng)默認(rèn)的采樣點數(shù)目為512點；H,f = freqz(b,a,n,Fs)：用于對在0，F(xiàn)s/2上等間隔采樣n點，采樣點頻率及相應(yīng)的頻響值分別存放在f和H中。H = freqz(b,a,w)：用于對在0，2上進(jìn)行采樣，采樣頻率點由w指定。H = freqz(b,a,f,Fs)：用于對在0，F(xiàn)s上進(jìn)行采樣，采樣頻率點由f指定。freqz(b,a,n)：用于在當(dāng)前圖形窗口中繪制幅頻和相頻特性曲線。下面介紹幾個實用的函數(shù)：mag=abs(H):求解系統(tǒng)的絕對幅頻響應(yīng)；db=20*log10(mag+eps)/max(mag)：求解系統(tǒng)的相對幅頻響應(yīng)；pha=angle(H

13、): 求解系統(tǒng)的相頻響應(yīng)；grd=grpdelay(b,a,w): 求解系統(tǒng)的群延遲。三、實驗內(nèi)容1.設(shè)某LTI系統(tǒng)的，輸入，求系統(tǒng)輸出:主程序如下：clear;clf;clc;n=-5:20; %設(shè)定一個n的觀察范圍h=delta(n-5);x=stepseq(0,-5,20)-stepseq(5,-5,20);y,ny=conv_m(x,n,h,n)subplot(2,2,1);stem(n,x);title('xn');xlabel('n');axis(-5,20,0,1.2);grid on;subplot(2,2,2);stem(n,h);title(

14、'hn');xlabel('n');axis(-5,20,0,1.2);grid on;subplot(2,2,3);stem(ny,y);title('yn');xlabel('n');axis(-5,20,0,1.2);grid on;delta.mfunction y=delta(x)y=(x=0);stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2)if nargin=3%nargin(number of arguments input)；“ =”表示不相等 disp('Us

15、age: Y=stepseq(n0,n1,n2)'); return;else if(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=n2')%n0<n1或者n0>n2或者n1>n2endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程序如下：functiony,ny=conv_m(x,nx,h,nh)%Modified convolution routine for signal processing%y,ny=conv_m(x

16、,nx,h,nh)%y=convolution result%ny=support of y%x=first signal on support nx%nx=support of x%h=second signal on support nh%nh=support of hif nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1); %nys beginingnye=nx(length(x)+nh(length(h); %nys endny=nyb:nye; %ny僅僅為了計算一下結(jié)果y對應(yīng)

17、的橫坐標(biāo)范圍y=conv(x,h); %Convolution and polynomial multiplication;conv為MATLAB的固有函數(shù)%conv_m函數(shù)比conv函數(shù)多用了用于表示橫坐標(biāo)范圍的ny、nx、nh,因為這里假定二維坐標(biāo)范圍nx和nh可能從負(fù)數(shù)開始，所以要重新計算y的橫坐標(biāo)范圍ny；如果nx和nh均為0:N，則可以直接得出ny為0:2N(正如conv函數(shù)中那樣）；在conv_m和conv函數(shù)中，x和h的橫坐標(biāo)范圍都可以實不相等的思考題：若,試寫出與的關(guān)系，并對MATLAB的仿真結(jié)果。解：因為函數(shù)是延時器，所以yn=xn* =x-。2.設(shè)某線性時不變系統(tǒng)的hn=

18、輸入信號為xn=求輸出：a ： b：（a）主程序如下：n=-5:20;u1=stepseq(0,-5,20);u2=stepseq(6,-5,20);%u1=un；u2=un-6%input xnx=u1-u2;%impulse response hnh=n.*x;subplot(3,1,1);stem(n,x);axis(-5 20 0 2);title('Input Sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-5 20 0 6);title('Inpulse Response'

19、);ylabel('hn'); %output responsey1,ny=conv_m(x,n,h,n);%conv_m為自定義求卷積的函數(shù)subplot(3,1,3);stem(ny,y1);title('Output Sequence');xlabel('n');ylabel('y_1n');（b）主程序如下：n=-10:20;u1=stepseq(0,-10,20);u2=stepseq(6,-10,20);%u1=un；u2=un-6%input xnx=u1-u2;u3=stepseq(-5,-10,20);u4=st

20、epseq(1,-10,20);%u3=un+5；u4=un-1x1=u3-u4;%impulse rsponse hn+5h=(n+5).*x1;subplot(3,1,1);stem(n,x);axis(-10 20 0 2);title('Input sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-10 20 0 6);title('Inpulse Response');ylabel('hn+5');%output responsey2,ny=conv_m(x,n

21、,h,n);subplot(3,1,3);stem(ny,y2);title('Output Sequence');xlabel('n');ylabel('y_2n');stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2&

22、#39;)endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程序如下：functiony,cy=conv_m(x,nx,h,nh)if nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1);nye=nx(length(x)+nh(length(x);ny=nyb:nye;y=conv(x,h);3.設(shè),輸入，求系統(tǒng)輸出:n=-5:50;u1=stepseq(0,-5,50);u2=stepseq(10,-5,50);%u1=un;u2=un-10%input xnx=u

23、1-u2;%impluse response hnh=(0.9).n).*u1;subplot(3,1,1);stem(n,x);axis(-5 20 0 2);title('Input Sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-5 20 0 6);title('Inpulse Response');ylabel('hn');%output responsey,ny=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y);title(&

24、#39;Output Sequence');xlabel('n');ylabel('yn');stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2')endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程

25、序如下：functiony,cy=conv_m(x,nx,h,nh)if nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1);nye=nx(length(x)+nh(length(x);ny=nyb:nye;y=conv(x,h);思考題：在此問題中為無限長序列，在計算卷積和時應(yīng)如何處理？解：取盡可能多的包含能量分布集中的具有代表性的點4某離散時間線性時不變系統(tǒng)的差分方程為：a.計算時的系統(tǒng)沖激響應(yīng) b.計算時的系統(tǒng)階躍響應(yīng)a=1,-0.8;b=1;n=-10:50;x=impseq(0,

26、-10,50);h=filter(b,a,x);%MATLAB固有函數(shù)filter表示系統(tǒng)方程ay=bx；x是沖擊信號，則h就是沖擊響應(yīng)subplot(2,1,1);stem(n,h);title('Inpulse Response');xlabel('n');ylabel('hn');x=stepseq(0,-10,50);s=filter(b,a,x); %x是階躍信號，則s就是階躍響應(yīng)subplot(2,1,2);stem(n,s);title('Step Response');xlabel('n');yla

27、bel('sn');impseq.m的源程序如下：functionx,n=impseq(n0,n1,n2)%Generates xn=delta(n-n0)；n1<=n<=n2,n1<=n0<=n2%x,n=impseq(n0,n1,n2)；if nargin=3 disp('Usage:Y=impseq(n0,n1,n2)');return;elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=n2')endn

28、=n1:n2;x=(n-n0)=0;stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2')endn=n1:n2;x=(n-n0)>=0; 連續(xù)時間系統(tǒng)分析1.設(shè)有兩個穩(wěn)定的LTI系統(tǒng)，分別可由下列微分方程來描述：；。請分別畫出它們的系統(tǒng)頻率響應(yīng)的幅值和相位特性曲線。lab41a.m的源程序如下：a=1 3; % y's coeff. in the system equationb=3; % x's coeff.in the system equationfreqs(b,a); % freqs為MATLAB中求連續(xù)時間系統(tǒng)拉普拉斯變換域（s域）的頻率響應(yīng)的固有函數(shù)Lab41b.m的源程序如下：a=3 4 1; % 系數(shù)矩陣b=1 0 5; % 系數(shù)矩陣freqs(b,a); % freqs為MATLAB中求連續(xù)時間系統(tǒng)拉普拉斯變換域（s域）的頻率響應(yīng)的固有函數(shù)思考題：求出系統(tǒng)a的頻率響應(yīng)表達(dá)式；根