弹簧-质量-阻尼系统的建模与控制系统设计

1、华北电力大学研究生结课作业分 数 : 任课教师签字：学年学期：第一学年第一学期课程名称：线性系统理论学生姓名：学号：提交时间： 2014.11.27华北电力大学（保定）线性系统理论结课报告目录目录2.1研究背景及意义3.2弹簧-质量-阻尼模型 31.1 系统的建立4.1.1.1 系统传递函数的计算51.2 系统的能控能观性分析71.2.1 系统能控性分析 81.2.2 系统能观性分析91.3 系统的稳定性分析 1.01.3.1 反馈控制理论中的稳定性分析方法 1 01.3.2 利用Matlab分析系统稳定性 101.3.3 Simulink 仿真结果1.21.4 系统的极点配置151.4.1

2、状态反馈法 151.4.2 输出反馈法 161.4.3 系统极点配置 1.61.5 系统的状态观测器 1.81.6 利用离散的方法研究系统的特性2 01.6.1 离散化定义和方法201.6.2 零阶保持器221.6.3 一阶保持器 241.6.4 双线性变换法263 .总结284 .参考文献287弹簧-质量-阻尼系统的建模与控制系统设计1研究背景及意义弹簧、阻尼器、质量块是组成机械系统的理想元件。由它们组成的弹簧-质量-阻尼系统是最常见的机械振动系统，在生活中具有相当广泛的用途，缓冲器 就是其中的一种。缓冲装置是吸收和耗散过程产生能量的主要部件，其吸收耗散 能量的能力大小直接关系到系统的安全与

3、稳定。 缓冲器在生活中处处可见，例如 我们的汽车减震装置和用来消耗碰撞能量的缓冲器，其缓冲系统的性能直接影响 着汽车的稳定与驾驶员安全；另外，天宫一号在太空实现交会对接时缓冲系统的 稳定与否直接影响着交会对接的成功。因此，对弹簧 -质量-阻尼系统的研究有着 非常深的现实意义。2弹簧-质量-阻尼模型数学模型是定量地描述系统的动态特性，揭示系统的结构、参数与动态特性 之间关系的数学表达式。其中，微分方程是基本的数学模型，不论是机械的、液 压的、电气的或热力学的系统等都可以用微分方程来描述。微分方程的解就是系 统在输入作用下的输出响应。所以，建立数学模型是研究系统、预测其动态响应 的前提。通常情况下

4、，列写机械振动系统的微分方程都是应用力学中的牛顿定 律、质量守恒定律等。弹簧-质量-阻尼系统是最常见的机械振动系统。机械系统如图2.1所示，图2-1弹簧-质量-阻尼系统机械结构简图其中m工、表示小车的质量，口表示缓冲器的粘滞摩擦系数，用表示弹簧的弹性系数，可 表示小车所受的外力，是系统的输入即%g 三耳5 ,% :匚表示小车的位移，是系统的输出，即 巧的二七地，i=1,2。设缓冲器的摩擦力与活塞的 速度成 正比，其中 = 11kg ,= 2kg , k，1 = ka = ICO N/cm ,k2 = 300N/cm, J = 5 = 3N - s/cm,匚& = 6 N，"c

5、m。2.1系统的建立由图2.1 ,根据牛顿第二定律，分别分析两个小车的受力情况，建立系统的动 力学模型如下：对有：+ cX + 必（匕一先）FkiXi + k2（X± 一） = P式力对口有：+c邑+ c2（X2 一冗+ /&一引=阻式。联立得到：对m工:叫5fl + （j + cj均一 q- + kJ% - 电/二 F<f对%:诏 + （% + 洱一勺氏+但+ %）与 - M 二 - 令勺= X1? xa = Xa ,/二 K【，七二先，Ui = Fl, U2 = Fs;71 =X1,%=沟得出状态空间表达式：兀=小所以，状态空间表达式为:由此可以得出00匕+电b1o

6、% 十 %吗4m-m2已知:代入数据得：0010,A= 0 OOI 40030096150-2003 -4.500 R= 0010- 00；Ip Q 0 0-lo 1 0 0-2.1,1系统传递函数的计算在Matlab中，函数ss2tf给出了状态空间模型所描述系统的传递函数，其一般形式是num,den=ss2tf(A,B,C,D,iu),其中 iu 是输入值。用Matlab将状态空间表达式表示为传递函数：在输入1单独作用的情况下A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0

7、 0;D=0 0;0 0;num,den=ss2tf(A,B,C,D,1)运行程序，得到：num =-0.00001.00004.5000200.0000-0.0000-0.00003.0000150.0000den =1.0e+004 *0.00010.00140.06230.18003.5000在输入2单独作用的情况下:A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;num,den=ss2tf(A,B,C,D,2)运行程序，得到:num =0-0

8、.0000-0.00003.0000 150.0000华北电力大学(保定)线性系统理论结课报告0-0.00000.50004.5000 200.0000den =1.0e+004 *0.00010.00140.06230.18003.5000由此可知：位移X1(t对外力Fg的传递函数是：s2 + 4.5s + 200F 就已s4+14s3 + 623s2 + 180(Js + 3S000位移X2也对外力F工俎的传递函数是：3s + 150Q s4414sa + 623s2 + 1800s -+ 3S000位移X1国对外力F侬的传递函数是：3s + ISO艰 s4+14sa + 623s2 +

9、ISOOs -+ 35000位移Xa(t：对外力的传递函数是：0.5s2 + 4.5s 4- 200s414s3 + 623s2 + IS 00s -+ 350002.2 系统的能控能观性分析在反馈控制理论中只讨论输入量对输出量的控制。而这两个量的关系唯一地由系统的传递函数所确定。一个稳定的系统，一定能控。同时，系统的输出量本 身就是我们想要控制的量，对于一个实际的系统来说，输出量当然是可以被观测 到的，因此在反馈控制理论中没有必要设立能控和能观这两个概念。然而在现代控制理论中，能控和能观是两个重要的基本概念。 我们把反映系 统内部运动状态的状态向量作为被控量， 而且它们不一定是实际上可观测到

10、的物 理量，至于输出量则是状态向量的线性组合， 这就产生了从输入量到状态量的能 控性问题和从输出量到状态量的能观测性问题。在现代控制中，分析和设计一个控制系统，必须研究这个系统的能控性和能 观性。状态方程描述了输入 U引起状态X (t)的变化过程；输出方程则描述 了由状态变化引起的输出 Y (t)的变化。能控性和能观性正是分别分析U对状态X (t)的控制能力以及Y (t)对X (t)的反应能力。2.2.1 系统能控性分析设线性定常系统的状态方程为大 屋二依 5 +BU .y 但 + du M式中Anx n矩阵Bnx r矩阵Cm x n矩阵Dm x r矩阵系统能控的充分必要条件为：能控判别阵口二

11、=BAB胪一】图的秩R(LJc)=n, 用Matlab计算能控矩阵的秩，从而对该系统的能控性进行判别，程序为：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;Qc=ctrb(A,B)R1=rank(Qc)运行程序，得到：R1 =4等于矩阵行数，由此可以判断，系统是完全能控的2.2.2 系统能观性分析设线性定常系统的状态方程为：发k二四m + Bu爆九的三口旧十口口加式中Anx n矩阵Bnx r矩阵Cm x n矩阵Dm x r矩阵能观的充分必要条件为：

12、能观判别阵 =C CACA111的秩R(LJC )=n , 下面，用Matlab计算能控矩阵的秩，从而对该系统的能控性进行判断： A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;Qo=obsv(A,C)R2=rank(Qo)运行程序，得到：R2 =4满秩，因此可以判断，该系统是完全能观的。综上所述，这是一个既能控又能观的系统。9华北电力大学(保定)线性系统理论结课报告2.3 系统的稳定性分析2.3.1 反馈控制理论中的稳定性分析方法稳定性是一个系统可

13、以被采用的最基本的条件，是系统的固有属性。稳定系统的定义如下：设控制系统处于某一起始的平衡状态， 在外力的作用 下，它离开了平衡状态，当外作用消失后，如果经过足够长的时间它能够恢复到 起始的平衡状态，则称这样的系统为稳定的系统，否则称为不稳定的系统。由稳 定性的定义可见，稳定性是系统去掉外力作用后自身的一种恢复能力，所以是系统的一种固有特性。对于线性定常系统，它取决于系统本身的结构和参数，而与 初始条件和外界作用无关。线性定常系统稳定的充分必要条件是：闭环系统特征方程的所有特征根为负 实数或具有负实部的共腕复数，即所有特征根位于复平面的左半平面。 只要有一 个闭环特征根分布在右半平面上，系统就

14、是不稳定的；如果没有右半平面的根， 但有纯虚根，则系统是临界稳定的；在工程上，处于不稳定和临界稳定的线性定 常系统是不能采用的1。在古典控制系统中，我们判断系统的稳定性经常用劳斯-赫尔维茨代数判据、 时域分析法、根轨迹法、频域分析法等方法，但那只针对低阶系统。实际的工业 生产中，经常会遇见一些特别复杂的系统。这时古典控制理论中的方法就有点捉 襟见肘了。1892年俄国学者李雅普诺夫提出的稳定性理论是确定系统稳定性的更一般 性理论，它采用了状态向量描述，不仅适用于单变量、线性、定常的系统，而且 适用于多变量，非线性、时变的系统。李雅普诺夫理论在建立一系列关于稳定性 概念的基础上，提出了判断系统稳定

15、性的两种方法： 一种方法是利用线性系统微 分方程的解来判断系统稳定性，称为李雅普诺夫第一法或间接法；另一种方法是 首先利用经验和技巧来构造李雅普诺夫函数，进而利用李雅普诺夫函数来判断系 统稳定性，称为李雅普诺夫第二法或直接法。2.3.2 利用Matlab分析系统稳定性随着计算机技术的发展，在现代控制理论中，我们经常采用 Matlab判断系 统的稳定性。对于线性定常系统，典型的系统输入信号类型有脉冲、阶跃、斜坡、 加速度、正弦信号。系统的稳定性是对任何输入信号而言，即若一个系统是稳定的，则其在任何输入信号情况下对应的输出曲线是收敛的。然而，阶跃信号包含了另外几种常见输入信号的特性，所以我们常通过

16、观察系统的单位阶跃响应曲线 判断判断系统的稳定性。若系统的单位阶跃响应是收敛的，则系统一般是收敛的； 否则，是发散的。在Matlab中输入相应系统的状态空间表达式矩阵来求取系统的特征值：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;eig(A)运行程序，得到：ans =-5.7735 +22.3859i-5.7735 -22.3859i-0.9765 + 8.0332i-0.9765 - 8.0332i由此可以知道，经计算得出 A阵的所有特征根均在

17、复平面的左半平面，因 此得出该系统是稳定的。给系统加起阶跃信号：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;step(A,B,C,D)结果如下0.01eaMTPmA0.0112tHU o- .OTStep ResponseFrom: In(1)0.0080.0060.0040.0020.008I0.0060.0040.002Time (sec)0From: In(2)图2-2阶跃响应曲线由图可以看出，在阶跃响应下，系统在一定时间内收敛于某一固定值

18、， 因此 可以判断系统是稳定的，但同时我们也可以看出，系统的调节时间比较长，如果 想要减少调节时间，那么需要重新配置极点，对系统进行改进。下面的章节将对 系统进行极点的配置。2.3.3 Simulink 仿真结果根据上述原理，用 Matlab中的Simulink组件进行仿真。根据状态空间表达式，搭建系统模型如下图所示：我们分别对只有输入1作用下和只有输入2作用下的系统使用Simulink进 行仿真，让其与Matlab图像进行对比15图2-3 Simulink模型图(1)仅有力作用时，系统的输出如下图所示图2-4 u1作用时响应曲线图中，绿色为输出1的曲线，蓝色为输出2的曲线。经分析：此曲线与对

19、应Matlab曲线一致，系统稳定，但是超调量较大，调节时间较长。(2)仅有叫作用，系统的输入如下所示:图中，绿色为输出1的曲线，蓝色为输出2的曲线。经分析：同样，此曲线 与对应的Matlab曲线一致，系统稳定，但是超调量较大，调节时间较长。在u工共同作用下，系统的输出如下图所示:图2-6 u1、u2同时作用时响应曲线图中绿色为输出1的曲线，蓝色为输出2的曲线。经分析：此曲线与Matlab 华北电力大学（保定）线性系统理论结课报告曲线一致，系统稳定，但是超调量较大，调节时间较长。需要进行极点配置，使 系统得到更好的性能。2.4 系统的极点配置控制系统的性能主要取决于系统极点在根平面上的分布。因此

20、，在系统设计 中，通常是根据对系统的品质要求，规定闭环极点应有的分布情况。 所谓的极点 配置就是，就是通过选择反馈矩阵 K,将闭环系统的极点恰好配置在根平面上所期望的位置，以获得所希望的动态性能。2.4.1 状态反馈法极点问题首先解决是否能通过状态反馈来实现给定的极点配置，即在什么条件下才有可能按照规定的要求来配置极点。其次是，这样的反馈阵K如何确定的问题。警号电入就3房曾示吉田图2-7状态反馈示意图17（1）采用状态反馈配置系统极点条件:系统（A，B，C） 采用状态反馈，任意配置其闭环系统极点的充要条件为：系 统（A，B，C）完全能控。 若系统不是完全能控的，就必须按能控性分解，只能任 意配

21、置可控的极点。（2）极点配置的方法：若原系统（A，B）可控，则采用状态反馈阵K,有A BK，B可控。nn 1设原系统的特征方程为san 1saisao 0o其中Kk0, k1,，kn 1，则有:010000kiA BKBankn配置后的闭环特征方程为:sn (ani kn1)sn1(a ki)s (a。 kg) 0 .假设闭环系统希望的极点为f ( ) s 1s 2 s1, 2, , n ,得到：nn 1s rn 1sr)s 0o为使系统达到希望性能，对比式（1）和式（2）中系数，使之相等，即可求 得状态反馈阵K k0,k1, ,kn1 。采用状态反馈配置系统极点不改变系统可控 性，它不能影响

22、系统中不可控部分模块。2.4.2 输出反馈法X =小 + Bhx = （A-RK）x +y = Cx= y = Cx图2-8输出反馈示意图对于完全能控的单变量系统，不能采用输出线性反馈来实现闭环系统极点任 意配置。不能任意配置极点，正是输出线性反馈的基本弱点。为了克服这个弱点， 在经典控制理论中，往往采取引入附加校正网络，通过增加开环零极点的方法改 变根轨迹走向，从而使其落在指定的期望位置上。对于完全能控的单变量系统（A，B，C） ,通过带动态补偿器的输出反馈时限极点任意配置的充要条件是：1.系统完全能观测；2.动态补偿器的阶数为n-1。 2.4.2系统极点配置在现代控制理论中是用系统内部的状

23、态来描述系统的， 所以经常从系统的状 态引出信号作为反馈量。利用状态反馈只能改变系统能控部分的极点， 而不能改 变不能控部分的极点，因此利用状态反馈进行极点配置的充分必要条件是系统必 须是完全能控的。对一个可控系统，在采用状态反馈后，可以实现闭环极点的任意配置，即通华北电力大学(保定)线性系统理论结课报告过状态反馈的方法，使闭环系统的极点位于任意期望的位置上。对于其中x是状态变量(n维)，u是控制信号，这里选取控制信号为因此，.二系统的稳态响应和瞬态响应特性由矩阵 A - BK的特征决定虽然理论上系统的闭环极点离S左半平面越远越好，但是在工业生产实践中，系统极点离左半平面越远，系统的运动状态就

24、变化的越快，这就要求执行机构快速运作，即使再妤的执行元件也会短时间内被损坏掉。所以新的极点的绝对值大约是原系统极点绝对值的 3至4倍左右。取P1= -15+40i; P2= -15-40i; P3= -3+10i; P4= -3-10i;利用Matlab进行极点配置，希望可以减小超调量，缩短稳定时间以优化系 统。Matlab程序如下：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;p=-15+40i,-15-40i,-3+10i,-3-10i;k=p

25、lace(A,B,p)step(A-B*k,B,C,D)运行程序，得到：k =-234.6522131.851214.45616.3957643.3762-89.97656.765836.0878#华北电力大学（保定）线性系统理论结课报告21Step Response-1.5 - 00.511.52 0-3X 10From: In(1)ouo- nT5 °12tHU o- .OTeaM«pmAFrom: In(2)0.511.52Time (sec)图2-9稳态响应曲线由响应曲线可以看出该系统重新配置极点后， 具有较快的调节时间，而且也 减少了超调量，改善了系统的动态性能与

26、稳态性能。2.5系统的状态观测器图2-10状态观测器示意图通过状态观测器可以任意配置系统的极点，从而使闭环系统具有期望的稳态和动态性能。但在工业生产中，系统的状态变量并非都是物理量， 或者是难以测 得的量。这样一来，系统的所有状态变量未必都可以直接测量得到，因此，状态反馈这种控制方式在许多实际控制问题中往往难以直接应用和实现。 状态观测器 就是利用系统的外部输入输出信息来确定系统内部的状态， 进而，在系统的极点 配置状态反馈中，用观测器得到的状态估计值代替系统的真实状态。 下图为状态 观测器的结构图：图2-11状态观测器示意图使用MATLAB为本系统设置状态观测器，选用极点配置时的极点，程序如

27、 下图所示：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;p=-15+40i,-15-40i,-3+10i,-3-10i;K1=place(A,B,p)A1=A-B*K1L1=(place(A',C',p)'A2=A-L1*CL2=(place(A1',C',p)'A3=A1-L2*Csys2=ss(A2,B,C,D)sys2=ss(A3,B,C,D)运行上面程序，得到：L1 =7.083330.08

28、95-30.579615.4167-41.6552-96.5401168.1877 200.0790A2 =-7.0833 -30.08951.000030.5796-15.416701.0000-358.3448396.5401-9.00006.0000-18.1877 -400.07903.0000-4.5000L2 =3.7432-7.1200-21.4563-3.7432190.989493.5822115.5037-24.2083A3 =-3.74327.12001.0000021.45633.743201.0000-655.5795 -119.9176-18.285630.9515

29、-81.9216 -402.0612-29.0527 -17.7144其中L1代表没进行状态反馈时的状态观测反馈矩阵，L2代表进行了状态反 馈的状态观测矩阵。2.6利用离散的方法研究系统的特性2.6.1 离散化定义和方法利用数字计算机对线性定常连续系统求数值解是现代科学技术研究中常用 的一种方法，它不但方便，而且精确。由于实际工业生产中线性定常连续系统被 控对象需要在线控制等，必须将连续时间系统的状态方程转化为离散系统的状态 方程，即将矩阵微分方程化成矩阵差分方程，这就是连续系统的离散化。根据离散系统的构成设备不同可以将离散系统分为采样控制系统和数字控制系统：a.采样控制系统：控制系统的构成中

30、选择了采样开关（或者含有开关特性的 设备）0b.数字控制系统：控制系统的控制器选择了专用数字计算机。通常，把系统中的离散信号是脉冲序列形式的离散系统， 称为采样控制系统 或脉冲控制系统；而把数字序列形式的离散系统，称为数字控制系统或计算机控 制系统。采样控制系统：采样控制系统是对来自传感器的连续信息在某些规定上的时 间瞬时值上取值。例如，控制器系统中的误差信号可以是断续连续的脉冲信号， 而相邻两个脉冲之间的误差信息，系统并没有收到。如果在有规律的间隔上，系 统取得了离散信息，则这种采样称为周期采样；反之，如果信息之间的间隔是时 变的，或随机的，则称为非周期采样，或随机采样。在采样控制系统中，把

31、连续信号转变为脉冲序列的过程称为采样过程，简称采样。实现采样的装置称为采样器，或采样开关。用T表示采样周期，单位为s; ffi = 1/T,表示采样频率，单位为1/E ； 5 =工叽=2TT/7表示采样角频率，单 位为rad/s。在采样控制系统中，把脉冲序列转变为连续信号的过程称为信号复 现过程。实现复现过程的装置称为保持器。采样周期的选择满足香农采样定理。 采样周期太大会使信号失真，采样周期 太小则容易造成计算过程的累积偏差或失去采样系统的特性。香农采样定理是在 设计离散系统时必须要遵循的准则，它给出了自采样的离散信号不失真地恢复原 连续信号所必需的理论上的最低采样频率。采样频率应该满足9自

32、2 叫J1HJC即是采样角频率应使其对连续信号中的最高频率分量，在一个周期内 被采样2次以上（上半周与下半周都至少采样一次），则采样后的脉冲序列中将 包含了连续信号的全部信息。但是，在仿真中所遇到的大多数被再现信号是没有 频带限的，所以一般取采样频率再现信号主要频带中的最高频率的510倍。在离散控制系统的设计过程中，采样周期的确定依据的是现场检测的被调量信号的 频率，对于频率较高的信号，采样周期的设定就小，而对于变化过程较慢的低频 华北电力大学(保定)线性系统理论结课报告信号，采样周期的设定可以大一些。有关概念在工程上的实际应用会有专门的内 容介绍。线性连续系统状态方程离散化的实质是将矩阵微分

33、方程化为矩阵差分方程， 它是描述多输入多输出离散系统的一种方便的数学模型。在推导离散化系统的方程时，假定系统是周期性采样，并且采样脉冲宽度远 小于采样周期，采样周期T的选择满足香农采样定理，还假设系统具有零阶保持特性，即在两个采样瞬间之间，采样值不变，并等于前一个采样时刻的值。通常离散化的方法有很多，伤J如欧拉法，梯形法，龙哥-库塔(Runge-Kutta) 法，阿达姆斯(Adams)法等等。下面主要运用三种方法来对系统进行离散化并 运用计算机进行模拟系统的特性，分析不同采样周期对系统的影响效果。2.6.2 零阶保持器零阶保持器可以将脉冲序列变成连续的方波信号，即将前一个采样周期的数值保留到下

34、一个采样点到来的时候。在 Matlab中输入函数如下：A=0 0 1 0;0 0 0 1; -400 300 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;p=-15+40i,-15-40i,-3+10i,-3-10i;k=place(A,B,p);H,I,J,K=c2dm(A-B*kBC,D,0.1,'zoh')dstep(H,I,J,K)分别设置采样时间为0.1s, 0.05s, 0.01s,运行程序，得到下图：23华北电力大学（保定）线性系统理论结课报告35Step Respons

35、el),ku o- .OT12tHU o- .OT eaMTPrnAt=0.1sStep Responsex 10-3From: In(1)From: In(2)12to o- OT eaMTP mAt=0.05sStep Responsex io -3From: In(1)From: In(2)32.5o1X21o 21o5-o- bT XZLO. OTeaMTPmA-1.5 050100150200 0Time (sec)50100150200t=0.01s图2-12零阶保持器离散化2.6.3 一阶保持器采用一阶保持器进行离散化，程序如下：A=0 0 1 0;0 0 0 1; -400 3

36、00 -9 6;150 -200 3 -4.5;B=0 0;0 0;1 0;0 0.5;C=1 0 0 0;0 1 0 0;D=0 0;0 0;p=-15+40i,-15-40i,-3+10i,-3-10i;k=place(A,B,p);H,I,J,K=c2dm(A-B*kBC,D,0.1,'foh')dstep(H,I,J,K)同样，分别设置采样时间为0.1s, 0.05s, 0.01s,运行程序，得到下图:Step Responset=0.1seaMTPmAeaMTP mA32.521.510.5021.510.50-0.5-1-1.501020304050 01020304050Time (sec)t=0.05sStep ResponseFrom: In(1)From: In(2)100150200 050100150200Time (sec)O1X3 5 2 5 1 5 O 2 1l),ku o- .OT50o2 5 1 5 0 5 1 5I , - -1 o o 1- - 12tHU o- .OTeaMTPmAt=0.01