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

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

2、， sys 是 LTI 系统模型，用来表示微分方程，差分方程或状态方程。其调用格式sys=tf(b,a)式中 ， b 和 a 分别是微分方程的右端和左端系数向量。例如，对于以下方程:a3 y'''(t )a2 y '' (t )a1 y' (t )a0 y(t )b3 f ''' (t )b2 f '' (t)b1 f ' (t )b0 f (t)可用 aa3 , a2 , a1, a0 ; bb3 , b2 ,b1, b0 ;systf (b, a) 获得其 LTI 模型 。注意 ，如果微分方程的

3、左端或右端表达式中有缺项，则其向量a 或 b 中的对应元素应为零 ，不能省略不写 ，否则出错 。.专业 .专注.例 3-1已知某 LTI 系统的微分方程为y (t)+ 2y (t)+100y(t)=f(t)其中 ， y(0)y' (0)0, f (t )10sin(2t ) ，求系统的输出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(

4、sec)');ylabel('y(t)');2 连续时间系统冲激响应和阶跃响应的求解在 MATLAB 中，对于连续 LTI 系统的冲激响应和阶跃响应 ，可分别用控制系统工具箱提供的函数 impluse 和 step 来求解 。 其调用格式为y=impluse(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

5、(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('Time(sec)');ylabel('g(t)');3. 用 MATLAB 实现连续时间信号的卷积信号的卷积运算有符号算法和数值算法，此处采用数值计算法，需调用MATLAB的conv( ) 函数近似计算信号的卷积积分。连续信号的卷积积分定义是f (t )f 1 (t)

6、f 2 (t )f1 ( ) f 2 (t)d如果对连续信号f1 (t) 和 f 2 (t ) 进行等时间间隔均匀抽样 ，则 f 1 (t) 和 f 2 (t ) 分别变为离散时间信号 f1 (m) 和 f 2 (m) 。 其中 ， m 为整数 。 当足够小时 ， f 1 (m) 和 f2 (m )既为连续时间信号f1 (t ) 和 f 2 (t) 。 因此连续时间信号卷积积分可表示为f (t )f 1 (t) f 2 (t )f1 () f 2 (t)df1 ( m )f 2 (tm )lim0 m.专业 .专注.采用数值计算时 ，只求当 tn时卷积积分f (t ) 的值 f (n) ，其中

7、 ， n 为整数 ，既f (n )f1 ( m )f 2 (nm )mf1 (m )f 2( nm)m其中 ，f 1 (m )f 2( nm) 实际就是离散序列f 1 (m ) 和 f 2 (m ) 的卷积和 。 当m足够小时 ，序列 f (n ) 就是连续信号f (t) 的数值近似 ，既f (t)f ( n ) f 1 (n)f 2 (n)上式表明 ，连续信号f 1 (t) 和 f 2 (t ) 的卷积 ，可用各自抽样后的离散时间序列的卷积再乘以抽样间隔。抽样间隔越小 ，误差越小 。例 3-3 用数值计算法求f 1( t)(t)(t2) 与 f2 (t )e 3t ( t) 的卷积积分 。解

8、：因为 f 2 ( t)e 3 t (t ) 是一个持续时间无限长的信号，而计算机数值计算不可能计算真正的无限长信号，所以在进行f 2 (t ) 的抽样离散化时，所取的时间范围让f 2 (t) 衰减到足够小就可以了 ，本例取 t2.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), grid on;axis(-1,2.5,-0.2,1.2); tit

9、le('f1(t)'); xlabel('t').专业 .专注.subplot(222), plot(t,f4), grid on;axis(-1,2.5,-0.2,1.2); title('f2(t)'); xlabel('t')subplot(223), plot(tt,f), grid on;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).*f2;subplot(224

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

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

12、s=tf(m,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=1 2 2;m=1 0;sys=tf(m,n);h=impulse(sys,t);figure;plot(t,h);xlabel('

13、;Time(sec)');ylabel('h(t)');.专业 .专注.g=step(sys,t);figure;plot(t,g);xlabel('Time(sec)');ylabel('g(t)');.专业 .专注.4 画出信号卷积积分f1 (t )f 2 (t) 的波形 ， f1 (t )f2 (t )(t)(t1)dt=0.01; t=-1:dt:3;f1=(t>=0);f2=(t>=1);f3=f1-f2;f=conv(f3,f3)*dt;n=length(f);.专业 .专注.tt=(0:n-1)*dt-2;sub

14、plot(221);plot(t,f1);grid on;axis(-1,3,-0.2,1.2);title('f1(t)');xlabel('t')subplot(222);plot(t,f2);grid on;axis(-1,3,-0.2,1.2);title('f2(t)');xlabel('t')subplot(223);plot(t,f3);grid on;axis(-1,3,-0.2,1.2);title('f3(t)');xlabel('t');subplot(224);plot(tt,f);grid on;.专业 .专注.title('f(t)=f3(t)*f3(t)');xlabel('t');四、实验要求1. 在 MATLAB 中输入程序 ，验证