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

《实验二 系统的时域分析报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、班级：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:50y

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

3、y3')subplot(2,2,4)t=0:0.001:100sys=tf([11],[123]);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]时系统的冲激响应、阶跃响应和零状态响应；程序代码：subplot(3,1,1)n=0:1:30h=impz([1],[1-10.9],n)stem(n,h)%绘

4、制系统的冲激响应title('冲激响应冲激响应')xlabel('n')ylabel('h')subplot(3,1,2)n=0:1:30g=stepz([1],[1-10.9],31)stem(n,g)%绘制阶跃响应title('阶跃响应')xlabel('n')ylabel('g')subplot(3,1,3)n=0:1:30f=n.^2y=filter([1],[1-10.9],f)stem(n,y)%绘制零状态响应title('零状态响应')xlabel('n')ylabel('y')运行结果显示与分析：3．已知x(n)=-1，h(n)＝0，y(n)=x(n)*h

5、(n)，求y(n)，并用分别画出离散序列x(n)、h(n)和y(n)的图形。程序代码：subplot(3,1,1)n=-1:1:3x=[12345]stem(n,x)%绘制离散序列x(n)title('离散序列x(n)')xlabel('n')ylabel('x')subplot(3,1,2)n=0:1:5h=[623642]stem(n,h)%%绘制离散序列h(n)title('离散序列h(n)')xlabel('n')ylabel('h')subplot(3,1,3)x=[12345]h=[623642]k=-1:1:8y=conv(x,h)stem(k,y)%%绘制

6、离散序列yn)title('离散序列y(n)')xlabel('k')ylabel('y')y=6142542635055522810运行结果显示与分析：4．求的部分分式展开。程序代码：b=[1.50.98-2.6081.2-0.144]a=[1-1.40.6-0.072][r,p,k]=residuez(b,a)运行结果显示与分析：r=0.70000.50000.3000p=0.60000.60000.2000k=02三、实验总结及心得：通过MATLAB函数库里的相应函数，来绘制连续和离散的系统冲激响应、阶跃响应及零状态响应的波形，以及卷积序列的画图，既了解了相应函数的

7、应用，又加深了对这些响应的理解。