離散系統(tǒng)的時域分析實驗報告

1、實驗2離散系統(tǒng)的時域分析一、實驗?zāi)康?、熟悉并掌握離散系統(tǒng)的差分方程表示法2、加深對沖激響應(yīng)和卷積分析方法的理解二、實驗原理在時域中，離散時間系統(tǒng)對輸入信號或者延遲信號進(jìn)行運(yùn)算處理，生成具有所需特性的輸出信號，具體框圖如下:其輸入、輸出關(guān)系可用以下差分方程描述：NM=SPit工和一氐=0;=o輸入信號分解為沖激信號,x?i=Vjcm現(xiàn)丹jkM=3-記系統(tǒng)單位沖激響應(yīng)二二一：-，貝孫統(tǒng)響應(yīng)為如下的卷積計算式:vpi=xt?琳燦丹=22兀加曲円m從=ODdj.=0.Ic=12-3_lV川：Q.A/l當(dāng)時，是有限長度的（），稱系統(tǒng)為系統(tǒng)；反之,稱系統(tǒng)為系統(tǒng)。三、實驗內(nèi)容1、用MATLAB求系統(tǒng)響應(yīng)1

2、)卷積的實現(xiàn)線性移不變系統(tǒng)可由它的單位脈沖響應(yīng)來表征。若已知了單位脈沖響應(yīng)和系統(tǒng)激勵就可通過卷積運(yùn)算來求取系統(tǒng)響應(yīng)，即y(n)x(n)*h(n)程序：x=input(Typeintheinputsequence=);%輸入xh=input(Typeintheimpulseresponsesequence=);%輸入hy=conv(x,h);%對x，h進(jìn)行卷積N=length(y)-1;%求出N的值n=0:1:N;%n從0開始，間隔為1的取值取到N為止disp(outputsequence=);disp(y);%輸出ystem(n,y);%畫出n為橫軸，y為縱軸的離散圖xlabel(Timein

3、dexn);ylable(Amplitude);%規(guī)定x軸y軸的標(biāo)簽輸入為：x=-201-13h=120-1圖形:21CDlllPE-CIE或2)單位脈沖響應(yīng)的線性時不變因果系統(tǒng)y=filter(b,a,x);其中，x和y是長度木濾波器的分子和分母加-取用maTl等的兩個矢量。矢量x表示激勵，矢量4a,b表示系統(tǒng)函數(shù)形式【數(shù)，得到的響應(yīng)為矢量y。例如計算以下系統(tǒng)的單位脈沖響應(yīng)y(n)+07y(n-1)-0.rj)厶廠v/.vf；)程序:N=input(Desiredimpuseresponselength=);b=input(Typeinthevectorb=);a=input(Typeint

4、hevectora=);x=1zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(Timeindexn);ylable(Amplitude);輸入：N=41b=0.8-0.440.360.02a=10.7-0.45-0.6圖形:、以下程序中分別1是否有差別，為個值，使用-05么要使用完成卷積并分卜值，D152025刃35節(jié)inrsince:-:ri5OILIPH4一-dEV程序:圖形OiipitObtaircyCor9lutijrpum-ldLLm因為在y=矢量中，丄若有3編制程丿Cu:pjtW已ieratedayzl:enc和PU3

5、1一遼匚式Tineiceci茅求解下列兩個系統(tǒng)的單位沖激響應(yīng)，分別用實現(xiàn)，并繪出其圖形。vk+0.75v?71H-0.125t2=jxfzs珂川1tk=0.25f城川14-3t?72+坯用一3+jcfx4給出理論計算結(jié)果和程序計算結(jié)果并討論。第一題實現(xiàn)程序N=input(Desiredimpuseresponselength=);b=input(Typeinthevectorb=);a=input(Typeinthevectora=);x=1zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(Timeindexn);ylabel(Am

6、plitude);圖形實現(xiàn)程序b=input(Typeinthevectorb=);a=input(Typeinthevectora=);N=25y=impz(b,a,N);k=0：1：N-1;stem(k,y);xlabel(Timeindexn);ylabel(Amplitude);圖形程序N=input(Desiredb=input(Typen-1第二題：-15實現(xiàn)-J5-JI-Ori_il-Q.LUV510152025responselgth=);a=input(Typeinthevectora=);x=1zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(Timeindexn);ylabel(Amplitude);輸入：Typeinthevectorb=0.250.250.250.25Typeinthevectora=1N=25圖形程序b=input(a=input(N=30n-山！3nl-clEV實現(xiàn)1120EO?ype辛?ypein)theveCtora=);5Tinr己iriceni丄Jy=impz(b,a,N);k=0:1:N-1;stem(k,y);xlabel(Timeindexn);ylabel(Am