实验二系统的时域分析报告

1、word 班级： 10通信2班 学号： Xb10680218 姓名： 张泽星 实验时间：2012年4月23日 地点：3号楼三楼 实验二 系统的时域分析一、实验目的：1. 掌握利用MATLAB求解连续和离散LTI系统的零状态响应方法。2. 加深对连续和离散LTI系统中冲激响应概念的理解，掌握其求解方法。3. 掌握离散卷积和的求解方法。4. 了解局部分式展开的求解方法。二、实验内容及运行结果：1给定微分方程y(t)+2y(t)+3y(t)=f(t)+f(t)，鼓励f(t)=et，绘制鼓励f(t)、系统冲激响应、阶跃响应及零状态响应的波形。程序代码subplot(2,2,1)t=0:0.001:50

2、y1=exp(t);plot(t,y1) %绘制鼓励f(t)=et的波形title('鼓励f(t)') %添加标题 xlabel('t')ylabel('y1') subplot(2,2,2)sys=tf(1 1,1 2 3);y2=impulse(sys,t)plot(t,y2) %绘制系统冲激响应的波形的波形 title('系统冲激响应') xlabel('t')ylabel('y2')subplot(2,2,3)sys=tf(1 1,1 2 3);y3=step(sys,t)plot(t,y3

3、) %绘制阶跃响应的波形 title('阶跃响应')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) %绘制零状态响应的波形title('零状态响应')xlabel('t')ylabel('y4')运行结果显示与分析2设离散系统可由以下差分方程表示：y(n)-y(n-1)+0.9y(n-2)=f(n)，且f(n)=n2。试分别绘制n=0:30时

4、系统的冲激响应、阶跃响应和零状态响应；程序代码：subplot(3,1,1)n=0:1:30h=impz(1,1 -1 0.9,n)stem(n,h) %绘制系统的冲激响应title('冲激响应冲激响应')xlabel('n')ylabel('h')subplot(3,1,2)n=0:1:30g=stepz(1,1 -1 0.9,31)stem(n,g) %绘制阶跃响应title('阶跃响应')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) %绘制零状态响应title('零状态响应')xlabel('n')ylabel('y')运行结果显示与分析：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运行结果显示与分析：4求的局部分式展开。程序代码：b=1.5 0.98 -2.608 1.2 -0.144a=1 -1.4 0.6 -0.072r,p,k=residuez(b,a)运行结果显示与分析：r = 0.7000 0.5000 0.3000p = 0.6000 0.6000 0.2000