连续时间LTI系统分析
- 格式:doc
- 大小:376.00 KB
- 文档页数:9
连续时间L T I系统分析(总8页) -本页仅作为预览文档封面,使用时请删除本页-实验三 连续时间LTI 系统分析一、实验目的(一)掌握使用Matlab 进行连续系统时域分析的方法1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应(二)掌握使用Matlab 进行连续时间LTI 系统的频率特性及频域分析方法1、学会运用MATLAB 分析连续系统的频率特性2、学会运用MATLAB 进行连续系统的频域分析(三)掌握使用Matlab 进行连续时间LTI 系统s 域分析的方法1、学会运用MATLAB 求拉普拉斯变换(LT )2、学会运用MATLAB 求拉普拉斯反变换(ILT )3、学会在MATLAB 环境下进行连续时间LTI 系统s 域分析二、实验条件装有MATLAB 的电脑三、实验内容(一)熟悉三部分相关内容原理(二)完成作业1、已知某系统的微分方程如下:)(3)()(2)(3)(t e t e t r t r t r +'=+'+''其中,)(t e 为激励,)(t r 为响应。
(1) 用MATLAB 命令求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应(零状态响应分别使用符号法和数值法求解,零输入响应只使用符号法求解);符号法求解零输入响应: >> eq='D2y+3*Dy+2*y=0';>> cond='y(0)=1,Dy(0)=2';>> yzi=dsolve(eq,cond);>> yzi=simplify(yzi)yzi =符号法求解零状态响应:exp(-2*t)*(4*exp(t) - 3)eq1='D2y+3*Dy+2*y=Dx+3*x';eq2='x=exp(-3*t)*heaviside(t)';cond='y=0,Dy=0';yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs)yzs =(exp(-2*t)*(exp(t) - 1)*(sign(t) + 1))/2图像如下:代码:subplot(211)ezplot(yzi,[0,8]);grid ontitle('ÁãÊäÈëÏìÓ¦')subplot(212)ezplot(yzs,[0,8]);grid ontitle('Áã״̬ÏìÓ¦')数值计算法:t=0::10;sys=tf([1,3],[1,3,2]);f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on ;axis([0 10 ]);title('ÊýÖµ¼ÆËã·¨µÄÁã״̬ÏìÓ¦')(2)使用MATLAB命令求出并画出系统的冲激响应和阶跃响应(数值法);用卷积积分法求系统的零状态响应并与(1)中结果进行比较;系统的冲激响应和阶跃响应(数值法):代码:t=0::10;sys=tf([1,3],[1,3,2]);h=impulse(sys,t);g=step(sys,t);subplot(211)plot(t,h),grid on;axis([0 10 ]);title('³å¼¤ÏìÓ¦')subplot(212)plot(t,g),grid on;axis([0 10 ]);title('½×Ô¾ÏìÓ¦'卷积积分法求系统的零状态响应:Ctsconv函数的定义:function[f,t]=ctsconv(f1,f2,t1,t2,dt)f=conv(f1,f2);f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(221)plot(t1,f1);grid onaxis([min(t1),max(t1),min(f1)-abs(min(f1)*,max(f1)+abs(max(f1)*])title('f1(t)');xlabel('t')subplot(222)plot(t2,f2);grid onaxis([min(t2),max(t2),min(f2)-abs(min(f2)*,max(f2)+abs(max(f2)*])title('f2(t)');xlabel('t')subplot(212)plot(t,f);grid onaxis([min(t),max(t),min(f)-abs(min(f)*,max(f)+abs(max(f)*])title('f(t)=f1(t)*f2(t)');xlabel('t')求系统的零状态响应代码:dt=;t1=0:dt:10;f1=exp(-3*t1).*uCT(t1);t2=t1;sys=tf([1,3],[1,3,2]);f2=impulse(sys,t2);[t,f]=ctsconv(f1,f2,t1,t2,dt)如图,根据两图相比较,两种方法做出的零状态响应大体相同。
MATLAB与信号实验——连续LTI系统的时域分析连续LTI系统的时域分析是信号与系统学中的重要课题。
MATLAB作为一种强大的科学计算软件,提供了丰富的工具和函数来进行信号与系统的分析。
下面将介绍MATLAB在连续LTI系统时域分析中的应用。
首先,我们需要了解连续LTI系统的基本概念。
一个连续域线性时不变系统(LTI系统)可以由它的冲激响应完全描述。
冲激响应是系统对单位冲激信号的响应。
在MATLAB中,可以使用impulse函数来生成单位冲激信号。
假设我们有一个连续LTI系统的冲激响应h(t),我们可以使用conv 函数来计算系统对任意输入信号x(t)的响应y(t)。
conv函数实现了卷积运算,可以将输入信号与冲激响应进行卷积运算得到输出信号。
例如,我们假设一个连续LTI系统的冲激响应为h(t) = exp(-t)u(t),其中u(t)是单位阶跃函数。
我们可以使用以下代码生成输入信号x(t)和计算输出信号y(t):```matlabt=-10:0.1:10;%时间范围x = sin(t); % 输入信号h = exp(-t).*heaviside(t); % 冲激响应y = conv(x, h, 'same'); % 计算输出信号```这段代码首先定义了时间范围t,然后定义了输入信号x(t)和冲激响应h(t)。
接下来,使用conv函数计算输入信号和冲激响应的卷积,设置参数’same’表示输出信号与输入信号长度相同。
最后,得到了输出信号y(t)。
在得到输出信号后,我们可以使用MATLAB的绘图功能来可视化结果。
例如,使用以下代码可以绘制输入信号和输出信号的图像:```matlabfigure;plot(t, x, 'b', 'LineWidth', 2); % 绘制输入信号hold on;plot(t, y, 'r', 'LineWidth', 2); % 绘制输出信号xlabel('时间');ylabel('幅度');legend('输入信号', '输出信号');```除了卷积运算外,MATLAB还提供了许多其他函数来进行连续LTI系统的时域分析。
[精品]连续时间LTI系统的频率特性及频域分析连续时间LTI系统(Linear Time-Invariant System)是指可用于描述各种物理和工程系统运动规律的动态系统。
它们由一对连续时变系统(如模型、结构和控制)和一对线性运算符构成,其具有因变量(响应)和自变量(输入)之间的线性关联性、时间不变性、结构连续的性质,并且在响应上呈现出定义的平稳性,因而它们在描述众多系统运动规律中被广泛应用。
对于连续时间LTI系统的频域特性的研究,则涉及这些系统的相位特性、幅频特性、切趾特性等。
同时,也要探讨系统中不同频率分量的传输特性,因为有不同频率分量的信号既可以幅频分析也可以相位分析,可以衡量系统不同频率下的相应响应。
由于连续时间LTI系统在有限频率通道内传播信号时发生了部分信号丢失,因此我们引入了频域分析得到系统频响阻抗。
这样一来,它就可以用来测量系统频带上的增益,系统的模态表现,以及系统的传播属性和可控特性。
在频域分析过程中,由于信号可以被分解为离散频率分量,所以对于单个频率分量来说,有关连续时间LTI系统的分析可以比较容易地完成。
一般情况下,每一个频率分量的传播特性由一个线性系数连接,称之为频响函数,可以衡量一个系统的频率响应情况。
总的来说,对于连续时间LTI系统,研究其频率特性及频域分析具有重要的意义。
他可以提供一个系统的相位特性、幅频特性、切趾特性等详细的分析,而且由于信号可以分解为离散频率分量,因此可以很容易地实现频域分析,并衡量一个系统的频率响应情况。
此外,还可以利用频域分析来测量系统的增益,模态表现,以及系统的传播属性和可控特性,进而提高系统的性能,实现性能的优化。
实验三连续时间LTI系统的时域分析实验报告一、实验目的通过实验三的设计和实现,达到如下目的:1、了解连续时间LTI(线性时不变)系统的性质和概念;2、在时域内对连续时间LTI系统进行分析和研究;3、通过实验的设计和实现,了解连续时间LTI系统的传递函数、共轭-对称性质、单位冲激响应等重要性质。
二、实验原理在常见的线性连续时间系统中,我们知道采用差分方程的形式可以很好地表示出该系统的性质和特点。
但是,在本实验中,我们可以采用微分方程的形式来进行相关的研究。
设系统的输入为 x(t),输出为 y(t),系统的微分方程为:其中,a0、a1、…、an、b0、b1、…、bm为系统的系数,diff^n(x(t))和diff^m(y(t))分别是输入信号和输出信号对时间t的n阶和m阶导数,也可以记为x^(n)(t)和y^(m)(t)。
系统的单位冲激响应函数 h(t)=dy/dx| x(t)=δ(t),则有:其中,h^(i)(t)表示h(t)的第i阶导数定义系统的传递函数为:H(s)=Y(s)/X(s)在时域内,系统的输出y(t)可以表示为:其中,Laplace^-1[·]函数表示Laplace逆变换,即进行s域到t域的转化。
三、实验步骤1、在Simulink中,构建连续时间LTI系统模型,其中系统的微分方程为:y(t)=0.1*x(t)-y(t)+10*dx/dt2、对系统进行单位冲激响应测试,绘制出系统的单位冲激响应函数h(t);4、在S函数中实现系统单位冲激响应函数h(t)的微分方程,并使用ODE45框图绘制出系统单位冲激响应函数h(t)在t=0~10s之间的图像;6、利用数据记录栏,记录系统在不同的参数下的变化曲线、阶跃响应函数u(t)和单位冲激响应函数h(t)的变化规律。
四、实验数据分析1、单位冲激响应测试那么,当输入信号为单位冲激函数δ(t)时,根据系统的微分方程,可以得知输出信号的形式为:即单位冲激响应函数h(t)为一个包含了单位冲激函数δ(t)在内的导数项序列。
MATLAB与信号实验-——-连续LTI系统的时域分析在信号处理中,MATLAB是一个强大的工具,它提供了许多功能,使我们能够模拟和分析各种信号系统。
对于连续LTI系统,时域分析是一个重要的方法,它允许我们直接观察系统的输入和输出信号之间的关系。
下面是一个关于连续LTI系统的时域分析的实验。
一、实验目的本实验的目的是验证连续LTI系统的时域响应,通过使用MATLAB模拟系统,我们可以观察到不同的输入信号产生的输出信号,从而了解系统的特性。
二、实验步骤1.定义系统:首先,我们需要定义我们的连续LTI系统。
这可以通过使用MATLAB中的lti函数来完成。
我们需要提供系统的传递函数,它描述了系统的输入和输出之间的关系。
2.设置输入信号:为了观察系统的行为,我们需要设置一个合适的输入信号。
在MATLAB中,我们可以使用square函数来生成一个方波信号,该信号具有固定的频率和幅度。
3.模拟系统:使用MATLAB的lsim函数,我们可以模拟我们的连续LTI系统。
这个函数将输入信号和系统的传递函数作为参数,然后计算出系统的输出信号。
4.分析结果:我们可以使用MATLAB的图形功能来观察输入和输出信号。
这可以帮助我们理解系统的行为,并验证我们的模型是否正确。
三、实验结果与分析在实验中,我们使用了不同的输入信号(如方波、正弦波等)来测试我们的连续LTI系统。
对于每种输入信号,我们都观察了系统的输出信号,并记录了结果。
通过对比不同的输入和输出信号,我们可以得出以下结论:1.对于方波输入,系统的输出信号是带有延迟的方波,这表明系统对突变信号的响应是瞬时的。
2.对于正弦波输入,系统的输出信号是与输入信号同频同相位的正弦波,这表明系统对正弦波的响应是具有稳定性的。
这些结果验证了连续LTI系统的基本特性:即对于单位阶跃函数(突变信号)的输入,系统的响应是瞬时的;而对于周期性输入(如正弦波),系统的响应具有稳定性。
这些结果与我们在理论上学到的知识相符,从而验证了我们的模型是正确的。
实验三连续时间LTI系统分析一、实验目的(一)掌握使用Matlab进行连续系统时域分析的方法1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应(二)掌握使用Matlab进行连续时间LTI系统的频率特性及频域分析方法1、学会运用MATLAB分析连续系统的频率特性2、学会运用MATLAB进行连续系统的频域分析(三)掌握使用Matlab进行连续时间LTI系统s域分析的方法1、学会运用MATLAB求拉普拉斯变换(LT)2、学会运用MATLAB求拉普拉斯反变换(ILT)3、学会在MATLAB环境下进行连续时间LTI系统s域分析二、实验条件装有MATLAB的电脑三、实验内容(一)熟悉三部分相关内容原理(二)完成作业1、已知某系统的微分方程如下:)(3)()(2)(3)(t e t e t r t r t r +'=+'+''其中,)(t e 为激励,)(t r 为响应。
(1) 用MATLAB 命令求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应(零状态响应分别使用符号法和数值法求解,零输入响应只使用符号法求解);符号法求解零输入响应: >> eq='D2y+3*Dy+2*y=0';>> cond='y(0)=1,Dy(0)=2';>> yzi=dsolve(eq,cond);>> yzi=simplify(yzi)yzi =符号法求解零状态响应:exp(-2*t)*(4*exp(t) - 3)eq1='D2y+3*Dy+2*y=Dx+3*x';eq2='x=exp(-3*t)*heaviside(t)';cond='y(-0.001)=0,Dy(-0.001)=0';yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs)yzs =(exp(-2*t)*(exp(t) - 1)*(sign(t) + 1))/2图像如下:代码:subplot(211)ezplot(yzi,[0,8]);grid ontitle('ÁãÊäÈëÏìÓ¦')subplot(212)ezplot(yzs,[0,8]);grid ontitle('Áã״̬ÏìÓ¦')数值计算法:t=0:0.01:10;sys=tf([1,3],[1,3,2]);f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on;axis([0 10 -0.001 0.3]);title('ÊýÖµ¼ÆËã·¨µÄÁã״̬ÏìÓ¦')(2)使用MATLAB命令求出并画出系统的冲激响应和阶跃响应(数值法);用卷积积分法求系统的零状态响应并与(1)中结果进行比较;系统的冲激响应和阶跃响应(数值法):代码:t=0:0.01:10;sys=tf([1,3],[1,3,2]);h=impulse(sys,t);g=step(sys,t);subplot(211)plot(t,h),grid on;axis([0 10 -0.01 1.1]);title('³å¼¤ÏìÓ¦')subplot(212)plot(t,g),grid on;axis([0 10 -0.01 1.6]);title('½×Ô¾ÏìÓ¦'卷积积分法求系统的零状态响应:Ctsconv函数的定义:function[f,t]=ctsconv(f1,f2,t1,t2,dt)f=conv(f1,f2);f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(221)plot(t1,f1);grid onaxis([min(t1),max(t1),min(f1)-abs(min(f1)*0.2),max(f1)+abs(max(f1)*0.2)])title('f1(t)');xlabel('t')subplot(222)plot(t2,f2);grid onaxis([min(t2),max(t2),min(f2)-abs(min(f2)*0.2),max(f2)+abs(max(f2)*0.2)])title('f2(t)');xlabel('t')subplot(212)plot(t,f);grid onaxis([min(t),max(t),min(f)-abs(min(f)*0.2),max(f)+abs(max(f)*0.2)])title('f(t)=f1(t)*f2(t)');xlabel('t')求系统的零状态响应代码:dt=0.01;t1=0:dt:10;f1=exp(-3*t1).*uCT(t1);t2=t1;sys=tf([1,3],[1,3,2]);f2=impulse(sys,t2);[t,f]=ctsconv(f1,f2,t1,t2,dt)如图,根据两图相比较,两种方法做出的零状态响应大体相同。
(3) 若已知条件同(1),借助MATLAB 符号数学工具箱实现拉普拉斯正反变换的方法求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应,并与(1)的结果进行比较。
普拉斯正反变换的方法求出系统的零状态响应和零输入响应:代码:syms t sRzis=(s+5)/(s^2+3*s+2);rzi=ilaplace(Rzis)rzi =4*exp(-t) - 3*exp(-2*t)et=exp(-3*t)*heaviside(t);es=laplace(et);Rzss=((3+s)*es)/(s^2+3*s+2);rzs=ilaplace(Rzss)rzs =exp(-t) - exp(-2*t)根据图像,同样也能看出拉普拉斯变换法得出的结果相同。
2、已知某RC 网络如下,()r t ()e t +-+-RC(1) 求出该网络的频域系统函数()H j ω;H (jw )=a/(a+jw) 其中a=1/RC(2) 使用MATLAB 命令画出1RC=时系统的幅频特性和相频特性;代码:w=-3*pi:0.01:3*pi;b=[0,1];a=[1,1];h=freqs(b,a,w);subplot(211)plot(w,abs(h)),grid onaxis([-10 10 0 1.1]);title('H(w)µÄ·ùÆµÌØÐÔ')subplot(212)plot(w,angle(h)),grid ontitle('H(W)µÄÏàÆµÌØÐÔ')(3) 若1RC =,且激励信号()sin sin(3)e t t t =+,使用频域分析法求解()r t ,分别画出()e t 和()r t 波形,讨论经传输是否引起失真。
代码:t=0:0.1:20;w1=1;w2=3;H1=1/(1+1i*w1);H2=1/(1+1i*w2);f=sin(t)+sin(3*t);y=abs(H1)*sin(w1*t+angle(H1))+abs(H2)*sin(w2*t+angle(H2));subplot(2,1,1);plot(t,f);grid onylabel('f(t)'),xlabel('Time(s)')title('ÊäÈëÐźŵIJ¨ÐÎ')subplot(2,1,2);plot(t,y);grid onylabel('y(t)'),xlabel('Time(sec)')title('ÎÈ̬ÏìÓ¦µÄ²¨ÐÎ')如图,两组波形进行比较可以明显看出,二者不成线性关系,所以此传输系统失真。
3、已知某系统框图如下,∑+-))((121-+s s K)(1s V )2s V(1) 写出下图所示系统的s 域系统函数()H s ;H (s )=1/(S 2+S-2+K)(2) 使用MATLAB 命令分别用两种方式画出90,1,2,,34K =时该系统的零极点分布图,并由图讨论K 从0增长时,该系统的稳定性变化情况。
代码:b1=[0 1];a1=[1 1 -2];sys1=tf(b1,a1);subplot(321)pzmap(sys1)axis([-2 2 -2 2])b2=[0 1];a2=[1 1 -1];sys1=tf(b2,a2);subplot(322)pzmap(sys1)axis([-2 2 -2 2])b3=[0 1];a3=[1 1 0];sys1=tf(b3,a3);subplot(323)pzmap(sys1)axis([-2 2 -2 2])b4=[0 1];a4=[1 1 0.25];sys1=tf(b4,a4);subplot(324)pzmap(sys1)axis([-2 2 -2 2])b4=[0 1];a4=[1 1 1];sys1=tf(b4,a4);subplot(325)pzmap(sys1)axis([-2 2 -2 2])根据图像,很明显的可以看出,随着K的逐渐增大,系统逐渐稳定。