1、word 班級： 10通信2班 學(xué)號： Xb10680218 姓名： 張澤星 實驗時間：2012年4月23日 地點：3號樓三樓 實驗二 系統(tǒng)的時域分析一、實驗?zāi)康模?. 掌握利用MATLAB求解連續(xù)和離散LTI系統(tǒng)的零狀態(tài)響應(yīng)方法。2. 加深對連續(xù)和離散LTI系統(tǒng)中沖激響應(yīng)概念的理解，掌握其求解方法。3. 掌握離散卷積和的求解方法。4. 了解局部分式展開的求解方法。二、實驗內(nèi)容及運行結(jié)果：1給定微分方程y(t)+2y(t)+3y(t)=f(t)+f(t)，鼓勵f(t)=et，繪制鼓勵f(t)、系統(tǒng)沖激響應(yīng)、階躍響應(yīng)及零狀態(tài)響應(yīng)的波形。程序代碼subplot(2,2,1)t=0:0.001:50

2、y1=exp(t);plot(t,y1) %繪制鼓勵f(t)=et的波形title('鼓勵f(t)') %添加標(biāo)題 xlabel('t')ylabel('y1') subplot(2,2,2)sys=tf(1 1,1 2 3);y2=impulse(sys,t)plot(t,y2) %繪制系統(tǒng)沖激響應(yīng)的波形的波形 title('系統(tǒng)沖激響應(yīng)') xlabel('t')ylabel('y2')subplot(2,2,3)sys=tf(1 1,1 2 3);y3=step(sys,t)plot(t,y3

3、) %繪制階躍響應(yīng)的波形 title('階躍響應(yīng)')xlabel('t')ylabel('y3')subplot(2,2,4)t=0:0.001:100sys=tf(1 1,1 2 3);f=exp(t);y4=lsim(sys,f,t)plot(t,y4) %繪制零狀態(tài)響應(yīng)的波形title('零狀態(tài)響應(yīng)')xlabel('t')ylabel('y4')運行結(jié)果顯示與分析2設(shè)離散系統(tǒng)可由以下差分方程表示：y(n)-y(n-1)+0.9y(n-2)=f(n)，且f(n)=n2。試分別繪制n=0:30時

4、系統(tǒng)的沖激響應(yīng)、階躍響應(yīng)和零狀態(tài)響應(yīng)；程序代碼：subplot(3,1,1)n=0:1:30h=impz(1,1 -1 0.9,n)stem(n,h) %繪制系統(tǒng)的沖激響應(yīng)title('沖激響應(yīng)沖激響應(yīng)')xlabel('n')ylabel('h')subplot(3,1,2)n=0:1:30g=stepz(1,1 -1 0.9,31)stem(n,g) %繪制階躍響應(yīng)title('階躍響應(yīng)')xlabel('n')ylabel('g')subplot(3,1,3)n=0:1:30f=n.2y=fi

5、lter(1,1 -1 0.9,f)stem(n,y) %繪制零狀態(tài)響應(yīng)title('零狀態(tài)響應(yīng)')xlabel('n')ylabel('y')運行結(jié)果顯示與分析：3x(n)= -1，h(n)0，y(n)=x(n)*h(n)，求y(n)，并用分別畫出離散序列x(n)、h(n)和y(n)的圖形。程序代碼：subplot(3,1,1)n=-1:1:3x=1 2 3 4 5stem(n,x) %繪制離散序列x(n)title('離散序列x(n)')xlabel('n')ylabel('x') subplo

6、t(3,1,2)n=0:1:5h=6 2 3 6 4 2stem(n,h) % %繪制離散序列h(n)title('離散序列h(n)')xlabel('n')ylabel('h') subplot(3,1,3)x=1 2 3 4 5h=6 2 3 6 4 2k=-1:1:8y = conv(x,h)stem(k,y) % %繪制離散序列yn)title('離散序列y(n)')xlabel('k')ylabel('y') y = 6 14 25 42 63 50 55 52 28 10運行結(jié)果顯示與分析：4求的局部分式展開。程序代碼：b=1.5 0.98 -2.608 1.2 -0.144a=1 -1.4 0.6 -0.072r,p,k=residuez(b,a)運行結(jié)果顯示與分析：r = 0.7000 0.5000 0.3000p = 0.6000 0.6000 0.2000