实验三连续时间LTI系统的时域分析报告

1、实验三连续时间LTI系统的时域分析一、实验目的1.学会用MATLAB求解连续系统的零状态响应；2. 学会用MATLAB求解冲激响应及阶跃响应；学会用MATLAB实现连续信号卷积的方法；二、实验原理1.连续时间系统零状态响应的数值计算我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，NM'、ay(t)八bjf(j)(t)i=Sj=9在MATLAB中，控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim。其调用格式y=lsim(sys,f,t)式中，t表示计算系统响应的抽样点向量，f是系统输入信号向量，sys是LTI系统模型，用来表示微分方程，差分方程或状态

2、方程。其调用格式sys=tf(b,a)式中，b和a分别是微分方程的右端和左端系数向量。例如，对于以下方程："'一一一'一a3y(t)a2y(t)ay(t)a°y(t)小3f(t)b?f(t)bif(t)bof(t)可用a=a3,a2,a，a。;b=b3,b2,bi,bo;sys=tf(b,a)获得其LTI模型。注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b中的对应元素应为零，不能省略不写，否则出错。例3-1已知某LTI系统的微分方程为y''(t)+2y'(t)+100y(t)=f(t)其中,y(0)=y'(0)=

3、0,f(t)=10sin(2兀t),求系统的输出y(t).解：显然，这是一个求系统零状态响应的问题。其MATLAB计算程序如下：ts=0;te=5;dt=0.01;sys=tf(1,1,2,100);t=ts:dt:te;f=10*sin(2*pi*t);y=lsim(sys,f,t);plot(t,y);xlabel('Time(sec)');ylabel('y(t)');2.连续时间系统冲激响应和阶跃响应的求解在MATLAB中，对于连续LTI系统的冲激响应和阶跃响应，可分别用控制系统工具箱提供的函数impluse和step来求解。其调用格式为y=implus

4、e(sys,t)y=step(sys,t)式中，t表示计算系统响应的抽样点向量，sys是LTI系统模型。例3-2已知某LTI系统的微分方程为y''(t)+2y'(t)+100y(t)=10f(t)求系统的冲激响应和阶跃响应的波形.解：ts=0;te=5;dt=0.01;sys=tf(10,1,2,100);t=ts:dt:te;h=impulse(sys,t);figure;plot(t,h);xlabel('Time(sec)');ylabel('h(t)');g=step(sys,t);figure;plot(t,g);xlabel(

5、'Time(sec)');ylabel('g(t)');3.用MATLAB实现连续时间信号的卷积信号的卷积运算有符号算法和数值算法，此处采用数值计算法，需调用MATLAB的conv()函数近似计算信号的卷积积分。连续信号的卷积积分定义是f(t)=fjt)f2(t)=f1()f2(t)d如果对连续信号f1(t)和f2(t)进行等时间间隔A均匀抽样，则f(t)和f2(t)分别变为离散时间信号f1(mA)和f2(mA)。其中，m为整数。当足够小时，f(mA)和f2(部)既为连续时间信号f(t)和f2(t)。因此连续时间信号卷积积分可表示为f(t)=f1(t)项)=Ji

6、(性。-)dQO=limf(m)f2(t-m：):采用数值计算时，只求当t=nA时卷积积分f(t)的值f(nA),其中，n为整数，既f(n.：)=Km.：)f2(n.：m：).:m=：oO-"：、'f1(m.：)f2(nm).：m-_:0其中，£f1(mA)，f2(nm)A实际就是离散序列f1(mA)和f2(mA)的卷积和。当m二.：A足够小时，序列f(nA)就是连续信号f(t)的数值近似，既f(t)：f(n')=：fi(n)f2(n)上式表明，连续信号fi(t)和f2(t)的卷积，可用各自抽样后的离散时间序列的卷积再乘以抽样间隔A。抽样间隔越小，误差越小。

7、例3-3用数值计算法求f1(t)=(t)&(t一2)与f2(t)=e项&(t)的卷积积分。解：因为f2(t)=e3&(t)是一个持续时间无限长的信号，而计算机数值计算不可能计算真正的无限长信号，所以在进行f2(t)的抽样离散化时，所取的时间范围让f2(t)衰减到足够小就可以了，本例取t=2.5。程序是dt=0.01;t=-1:dt:2.5;f1=(t>=0);f2=(t>=2);f3=f1-f2f4=exp(-3*t).*f1;f=conv(f3,f4)*dtn=length(f);tt=(0:n-1)*dt-2;subplot(221),plot(t,f3

8、),gridon;axis(-1,2.5,-0.2,1.2);title('f1(t)');xlabel('t')subplot(222),plot(t,f4),gridon;axis(-1,2.5,-0.2,1.2);title('f2(t)');xlabel('t')subplot(223),plot(tt,f),gridon;title('f(t)=f1(t)*f2(t)');xlabel('t')y=(1/3-1./3*exp(-3*t).*f1-(1/3-1./3*exp(-3*(t-2)

9、.*f2;subplot(224),plot(t,y),gridontitle('f(t)=f1(t)*f2(t)');xlabel('t')由于f(t)和f2(t)的时间范围都是从t=-1开始，所以卷积结果的时间范围从t=-2开始，增量还是取样间隔A,这就是语句tt=(0:n-1)*dt-2的由来。三、上机实验内容验证实验原理中所述的相关程序已知描述系统的微分方程和激励信号f(t)如下，试用解析法求系统的零状态响应y(t),并用MATLAB绘出系统零状态响应的时域仿真波形，验证结果是否相同y''(t)+4y'(t)+4y(t)=ff(t

10、)=te+W)£(t)ts=0;te=5;dt=0.01;f1=(t>=0);sys=tf(13,144);t=ts:dt:te;f=exp(-t).*f1;y=lsim(sys,f,t);plot(t,y);xlabel('Time(sec)');ylabel('y(t)');已知描述系统的微分方程如下，试用MATLAB求系统在010秒范围内冲激响应和阶跃响应的数值解，并用绘出系统冲激响应和阶跃响应的时域波形y''(t)+3y'(t)+2y(t)=f(t)(1):t=0:0.01:5;n=132;m=1;sys=tf(m

11、,n);h=impulse(sys,t);figure;plot(t,h);xlabel('Time(sec)');ylabel('h(t)');g=step(sys,t);figure;plot(t,g);xlabel('Time(sec)');ylabel('g(t)');y''(t)+2y'(t)+2y(t)=f'(t)t=0:0.01:5;n=122;m=10;sys=tf(m,n);h=impulse(sys,t);figure;plot(t,h);xlabel('Time(sec

12、)');ylabel('h(t)');g=step(sys,t);figure;plot(t,g);xlabel('Time(sec)');ylabel('g(t)');0.35.0Q511111111100.511.522533.544.55Time(sec)4.画出信号卷积积分f1(t)*f2(t)的波形，f(t)=f2(t)=&(t)&(t1)f2=(t>=1);f3=f1-f2;f=conv(f3,f3)*dt;n=length(f);subplot(221);plot(t,f1);gridon;axis(-

13、1,3,-0.2,1.2);title('f1(t)');xlabel('t')subplot(222);plot(t,f2);gridon;axis(-1,3,-0.2,1.2);title('f2(t)');xlabel('t')subplot(223);plot(t,f3);gridon;axis(-1,3,-0.2,1.2);title('f3(t)');xlabel('t');subplot(224);plot(tt,f);gridon;xlabel('t');Q.5f1(t)1003f2(t)10.5002Q.Ef3(t)1031t四、实验