離散LSI系統(tǒng)地時域分析報告

1、文檔 ：L張阿讓窗野釈頁麗蛋實驗實用報告課程名稱數(shù)字信號實驗成績_指導(dǎo)教師 _實驗二：離散 LSI系統(tǒng)的時域分析一、實驗內(nèi)容1.知描述某離散LSI系統(tǒng)的差分方程為2y(n)-3y(n-1)+y(n-2)=x(n-1)，分別用impz和dstep函數(shù)、filtic 和filter函數(shù)兩種方法求解系統(tǒng)的單位序列響應(yīng)和單位階躍響應(yīng)。用impz和dstep函數(shù)求解系統(tǒng)的單位序列響應(yīng)和單位階躍響應(yīng)如下a=1,-3/2,1/2;b=0,1/2,0;N=32;n=0:N-1;hn=impz(b,a, n);gn=dstep(b,a ,n);subplot(1,2,1);stem( n,h n,k);titl

2、e(系統(tǒng)的單位序列響應(yīng));ylabel(h( n);xlabel( n);axis(0,N,1.1*mi n(h n),1.1*max(h n);subplot(1,2,2);stem( n,g n,k);title(系統(tǒng)的單位階躍響應(yīng));ylabel(g( n);xlabel( n);axis(0,N,1.1*mi n(g n),1.1*max(g n);實用文檔解：xO仁0;y0仁0;a=1,-3/2,1/2;b=1/2,0,0;N=32 ;n=0:N-1;xi=filtic(b,a,O);x1= n=0;hn=filter(b,a,x1,xi);x2=n=0;gn=filter(b,a,

3、x2,xi);subplot(1,2,1);stem( n,h n,k);title( 系統(tǒng)的單位序列響應(yīng));ylabel(h( n);xlabel( n);axis(0,N,1.1*mi n(h n),1.1*max(h n);系統(tǒng)的單位序列響應(yīng)系統(tǒng)的單位階躍響應(yīng)用函數(shù)filtic 和filter求解離散系統(tǒng)的單位序列響應(yīng)和單位階躍實用文檔subplot(1,2,2);stem( n,g n,k);title(系統(tǒng)的單位階躍響應(yīng));ylabel(g( n);xlabel( n);axis(0,N,1.1*mi n(g n),1.1*max(g n);2. 寫程序描繪下列序列的卷積波形：(1)

4、 f 1(n)=u(n),f 2(n)=u(n-2), (0 n10)n1=0:10;n t=le ngth( n1);f1=on es(1, nt);n2=2:12;n t=le ngth( n2);f2=o nes(1, nt);y, ny=co nvu(f1, n1,f2, n2);系統(tǒng)的單位序列響應(yīng)系統(tǒng)的單位階躍響應(yīng)實用文檔(2) x(n)=sin(n/2),h(n)=(0.5) n1=-3:4*pi;f1=0.84 n1;f2=si n(n 2/2);y, ny=co nvu(f1, n1,f2, n2); subplot(2,2,1);stem( n1,f1); subplot(2

5、,2,2);stem( n2,f2); subplot(2,1,2);stem( ny ,y);subplot(2,2,1);stem( n1,f1);subplot(2,2,2);stem( n2,f2);subplot(2,1,2);stem( ny,y);定義函數(shù)文件調(diào)用部分：fun ctio ny, ny= con vu(f1, n1,f2, n2)nys=n 1(1)+n2(1); nyf=n 1(e nd)+n2(e nd);y=c on v(f1,f2); ny=nys:nyf;(-3實用文檔定義函數(shù)文件調(diào)用部分:fun ctio ny, ny= co nvu(f1, n1,f2

6、, n2)nys=n 1(1)+n2(1); nyf=n 1(e nd)+n2(e nd);y=c on v(f1,f2); ny=nys:nyf;3.知某離散LSI系統(tǒng)的單位序列響應(yīng)為h(n)=3 S (n-3)+0.5 S (n-4)+0.2 S (n-5)+0.7 S (n-6)-0.8S (n-7)求輸入為x(n)=e -0.5n u(n)時的系統(tǒng)響應(yīng)。N=16;n=0:N-1;x=exp(-0.5* n);subplot(2,2,1);stem( n,x);a=1;a實用文檔b=0,0,0,3,0.5,0.2,0.7,0.8;hn=impz(b,a, n); subplot(2,2,

7、2);stem( n,hn) y=c on v(x,h n);subplot(2,1,2);stem(y);4.描述某離散LSI系統(tǒng)的差分方程為y(n)=0.7y(n-1)+2x(n)-x(n-2)，求輸入為x(n)=u(n-3)時的系統(tǒng)響應(yīng)。N=16 ;n 1=3:N+2;f1=zeros(1,3),o nes(1,(N-3); subplot(2,2,1);stem( n1,f1);a=1,0,-0.7,0;b=2,0,-1,0;f2=impz(b,a, n1);subplot(2,2,2);stem( n1,f2);y, ny=co nvu(f1, n1,f2, n1);subplot(

8、2,1,2);stem( ny ,y);10.80.60.40.20051015實用文檔定義函數(shù)文件調(diào)用部分：fun ctio ny, ny= con vu(f1, n1,f2, n2)nys=n 1(1)+n2(1); nyf=n 1(e nd)+n2(e nd);y=c on v(f1,f2); ny=nys:nyf;、思考題答案思考本實驗提出的有關(guān)MATLAB函數(shù)在調(diào)用時應(yīng)注意哪些問題。調(diào)用自定意函數(shù)時按照以下代碼在自定義函數(shù)的取值不是從0開始時就引用一下函數(shù)：fun ctio ny, ny=c on vu(h, nh,x, nx)nys=n h(1)+nx(1); nyf=n h(e

9、nd)+nx(e nd);y=c on v(h,x); ny=nys:nyf;實用文檔里面的變量要 -對應(yīng)代碼寫到function 的文件里面調(diào)用時convu與文件命名一致！、實驗小結(jié)(1)遇到問題及解決方法注意函數(shù)的定義出現(xiàn)問題及時問老師(2)學(xué)到的新的函數(shù)及其用法由離散時間系統(tǒng)的時域分析方法可知，一個離散LSI系統(tǒng)的響應(yīng)與激勵可以用如下框圖表示:xnDiscrete-timeynsystme其輸入2、用函數(shù)impz和dstep求解離散系統(tǒng)的單位脈沖響應(yīng)和單位階躍響應(yīng)。NM、輸出關(guān)系可用以下差分方程描述：、二akyn -k = bkxn -mk =0k =03、 用函數(shù)filtic 和fil

10、ter求解離散系統(tǒng)的單位序列響應(yīng)和單位階躍響應(yīng)。4、用MATLAB實現(xiàn)線性卷積1) 用函數(shù)conv進(jìn)行卷積運算：求解兩個序列的卷積和，關(guān)鍵在于如何確定卷積結(jié)果的時寬區(qū)間。MATLAB提供的求卷積函數(shù)conv默認(rèn)兩個序列的序號均從 n=0開始，卷積結(jié)果y對應(yīng)的序列的序號也從 n=0開 始。2) 非零起始序列的卷積運算：當(dāng)兩個序列不是從 0開始時，必須對conv函數(shù)稍加擴(kuò)展。由卷積原理可知，若待卷積的兩個序列序號分別為x(n) ; nx=nxs:nxf ，h(n)；nh=nhs:nhf，則卷積和 y(n)的序號起點和終點分別為：nys=nxs+nhs， nyf=nxf+nhf。據(jù)此可定義通用卷積函數(shù) convu :fun ctio ny, ny=con vu(h, nh,x, nx)nys=n h(1)+nx(1); nyf=n h(e nd)+nx(e nd);y=c on v