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

1、分数:任课教师签字：华北电力大学研究生结课作业学年学期：第一学年第一学期课程名称：线性系统理论学生姓名：学号：提交时目录目录1.1研究背景及意义3.2弹簧-质量-阻尼模型32.1 系统的建立32.1.1 系统传递函数的计算42.2 系统的能控能观性分析62.2.1 系统能控性分析62.2.2 系统能观性分析72.3 系统的稳定性分析7.2.3.1 反响限制理论中的稳定性分析方法72.3.2 利用Matlab分析系统稳定性82.3.3 Simulink仿真结果92.4 系统的极点配置102.4.1 状态反响法102.4.2 输出反响法112.4.3 系统极点配置112.5 系统的状态观测器1.3

2、2.6 利用离散的方法研究系统的特性152.6.1 离散化定义和方法152.6.2 零阶保持器162.6.3 一阶保持器172.6.4 双线性变换法183 .总结184 .参考文献1.9弹簧-质量-阻尼系统的建模与限制系统设计1研究背景及意义弹簧、阻尼器、质量块是组成机械系统的理想元件.由它们组成的弹簧-质量-阻尼系统是最常见的机械振动系统,在生活中具有相当广泛的用途,缓冲器就是其中的一种.缓冲装置是吸收和耗散过程产生能量的主要部件,其吸收耗散能量的水平大小直接关系到系统的平安与稳定.缓冲器在生活中处处可见,例如我们的汽车减震装置和用来消耗碰撞能量的缓冲器,其缓冲系统的性能直接影响着汽车的稳定

3、与驾驶员平安；另外,天宫一号在太空实现交会对接时缓冲系统的稳定与否直接影响着交会对接的成功.因此,对弹簧-质量-阻尼系统的研究有着非常深的现实意义.2弹簧-质量-阻尼模型数学模型是定量地描述系统的动态特性,揭示系统的结构、参数与动态特性之间关系的数学表达式.其中,微分方程是根本的数学模型,不管是机械的、液压的、电气的或热力学的系统等都可以用微分方程来描述.微分方程的解就是系统在输入作用下的输出响应.所以,建立数学模型是研究系统、预测其动态响应的前提.通常情况下,列写机械振动系统的微分方程都是应用力学中的牛顿定律、质量守恒定律等.弹簧-质量-阻尼系统是最常见的机械振动系统.机械系统如图2.1所示

4、,图2-1弹簧-质量-阻尼系统机械结构简图其中口卜m?表示小车的质量,q表示缓冲器的粘滞摩擦系数,4表示弹簧的弹性系数,表示小车所受的外力,是系统的输入即示小车的位移,是系统的输出,即i=1,2o设缓冲器的摩擦力与活塞的速度成正比,其中叫=lkg,1nL2kg,及二%二1.0肌由,卜=300班加,yc广3Ws/cm,y6M/cm.2.1 系统的建立由图2.1,根据牛顿第二定律,分别分析两个小车的受力情况,建立系统的动力学模型如下：那么有:嘲+瑞+q仅-+叫+-切二F®对吨有：叫用谓小2(卜M+匾+磷27卜明联立得到:叫:1nx+(q+与况-遥+011+k2)Xt-k2X2=Ft(t)

5、对m：:m2XE+(Cj+口)%-5i+8+k温-k2Kl=F2(t)得出状态空间表达式:所以,状态空间表达式为:由此可以得出：叫二1,叱二2,&二k100,k?=300,q二勺二3,q二6代入数据得：2.1.1 系统传递函数的计算在Matlab中,函数ss2tf给出了状态空间模型所描述系统的传递函数,其一般形式是num,den=ss2tf(A,B,C,D,iu),其中iu是输入值.用Matlab将状态空间表达式表示为传递函数：在输入1单独作用的情况下A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00

6、;00;num,den=ss2tf(A,B,C,D,1)运行程序,得到：num=0-0.00001.00004.5000200.00000-0.0000-0.00003.0000150.0000den=1.0e+004*0.00010.00140.06230.18003.5000在输入2单独作用的情况下:A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;num,den=ss2tf(A,B,C,D,2)运行程序,得到：num=0-0.0000-0.00003.0000150.00000-0.00000.

7、50004.5000200.0000den=1.0e+004*0.00010.00140.06230.18003.5000由此可知：的传递函数是:位移对外力位的传递函数是:位移X吊对外力Fw的传递函数是:位移Aw,对外力的传递函数是:2.2 系统的能控能观性分析在反响限制理论中只讨论输入量对输出量的限制.而这两个量的关系唯一地由系统的传递函数所确定.一个稳定的系统,一定能控.同时,系统的输出量本身就是我们想要限制的量,对于一个实际的系统来说,输出量当然是可以被观测到的,因此在反响限制理论中没有必要设立能控和能观这两个概念.然而在现代限制理论中,能控和能观是两个重要的根本概念.我们把反映系统内部

8、运动状态的状态向量作为被控量,而且它们不一定是实际上可观测到的物理量,至于输出量那么是状态向量的线性组合,这就产生了从输入量到状态量的能控性问题和从输出量到状态量的能观测性问题.在现代限制中,分析和设计一个限制系统,必须研究这个系统的能控性和能观性.状态方程描述了输入U引起状态X(t)的变化过程；输出方程那么描述了由状态变化引起的输出Y(t)的变化.能控性和能观性正是分别分析U对状态X(t)的限制水平以及Y(t)对X(t)的反响水平.2.2.1 系统能控性分析设线性定常系统的状态方程为式中Anxn矩阵BnXr矩阵Cmxn矩阵Dmxr矩阵系统能控的充分必要条件为：能控判别阵U二B的秩R(Up=n

9、,用Matlab计算能控矩阵的秩,从而对该系统的能控性进行判别,程序为：A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;Qc=ctrb(A,B)R1=rank(Qc)运行程序,得到：R1=4等于矩阵行数,由此可以判断,系统是完全能控的.2.2.2 系统能观性分析设线性定常系统的状态方程为：式中Anxn矩阵BnXr矩阵Cmxn矩阵Dmxr矩阵能观的充分必要条件为：能观判别阵UD=CCA“.CA"T卜勺秩R(U0)=n,下面,用Matlab计算能控矩阵的秩,从而对该系统的能控性进行判断：A=0

10、010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;Qo=obsv(A,C)R2=rank(Qo)运行程序,得到：R2=4满秩,因此可以判断,该系统是完全能观的.综上所述,这是一个既能控又能观的系统.2.3 系统的稳定性分析2.3.1 反响限制理论中的稳定性分析方法稳定性是一个系统可以被采用的最根本的条件,是系统的固有属性.稳定系统的定义如下：设限制系统处于某一起始的平衡状态,在外力的作用下,它离开了平衡状态,当外作用消失后,如果经过足够长的时间它能够恢复到起始的平衡状态,那么称这样的系统为稳定的系统,否那么

11、称为不稳定的系统.由稳定性的定义可见,稳定性是系统去掉外力作用后自身的一种恢复水平,所以是系统的一种固有特性.对于线性定常系统,它取决于系统本身的结构和参数,而与初始条件和外界作用无关.线性定常系统稳定的充分必要条件是：闭环系统特征方程的所有特征根为负实数或具有负实部的共腕复数,即所有特征根位于复平面的左半平面.只要有一个闭环特征根分布在右半平面上,系统就是不稳定的；如果没有右半平面的根,但有纯虚根,那么系统是临界稳定的；在工程上,处于不稳定和临界稳定的线性定常系统是不能采用的1.在古典限制系统中,我们判断系统的稳定性经常用劳斯-赫尔维茨代数判据、时域分析法、根轨迹法、频域分析法等方法,但那只

12、针对低阶系统.实际的工业生产中,经常会遇见一些特别复杂的系统.这时古典限制理论中的方法就有点捉襟见肘了.1892年俄国学者李雅普诺夫提出的稳定性理论是确定系统稳定性的更一般性理论,它采用了状态向量描述,不仅适用于单变量、线性、定常的系统,而且适用于多变量,非线性、时变的系统.李雅普诺夫理论在建立一系列关于稳定性概念的根底上,提出了判断系统稳定性的两种方法：一种方法是利用线性系统微分方程的解来判断系统稳定性,称为李雅普诺夫第一法或间接法；另一种方法是首先利用经验和技巧来构造李雅普诺夫函数,进而利用李雅普诺夫函数来判断系统稳定性,称为李雅普诺夫第二法或直接法.2.3.2 利用Matlab分析系统稳

13、定性随着计算机技术的开展,在现代限制理论中,我们经常采用Matlab判断系统的稳定性.对于线性定常系统,典型的系统输入信号类型有脉冲、阶跃、斜坡、加速度、正弦信号.系统的稳定性是对任何输入信号而言,即假设一个系统是稳定的,那么其在任何输入信号情况下对应的输出曲线是收敛的.然而,阶跃信号包含了另外几种常见输入信号的特性,所以我们常通过观察系统的单位阶跃响应曲线判断判断系统的稳定性.假设系统的单位阶跃响应是收敛的,那么系统一般是收敛的；否那么,是发散的.在Matlab中输入相应系统的状态空间表达式矩阵来求取系统的特征值：A=0010;0001;-400300-96;150-2003-4.5;B=0

14、0;00;10;00.5;C=1000;0100;D=00;00;eig(A)运行程序,得到：ans=-5.7735+22.3859i-5.7735-22.3859i-0.9765+8.0332i-0.9765-8.0332i由此可以知道,经计算得出A阵的所有特征根均在复平面的左半平面,因此得出该系统是稳定的.给系统加起阶跃信号：A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;step(A,B,C,D)结果如下图2-2阶跃响应曲线由图可以看出,在阶跃响应下,系统在一定时间内收敛于某一固定值,因此可以

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

16、致,系统稳定,但是超调量较大,调节时间较长.(2)仅有年作用,系统的输入如下所示：图2-5u2作用时响应曲线图中,绿色为输出1的曲线,蓝色为输出2的曲线.经分析：同样,此曲线与对应的Matlab曲线一致,系统稳定,但是超调量较大,调节时间较长.在口,11共同作用下,系统的输出如下列图所示：图2-6u1、u2同时作用时响应曲线图中绿色为输出1的曲线,蓝色为输出2的曲线.经分析：此曲线与Matlab曲线一致,系统稳定,但是超调量较大,调节时间较长.需要进行极点配置,使系统得到更好的性能.2.4 系统的极点配置限制系统的性能主要取决于系统极点在根平面上的分布.因此,在系统设计中,通常是根据对系统的品

17、质要求,规定闭环极点应有的分布情况.所谓的极点配置就是,就是通过选择反响矩阵K,将闭环系统的极点恰好配置在根平面上所期望的位置,以获得所希望的动态性能.2.4.1 状态反响法极点问题首先解决是否能通过状态反响来实现给定的极点配置,即在什么条件下才有可能根据规定的要求来配置极点.其次是,这样的反响阵K如何确定的问题.图2-7状态反响示意图(1)采用状态反响配置系统极点条件：系统(A,B,C)采用状态反响,任意配置其闭环系统极点的充要条件为：系统,(A,B,C)完全能控.假设系统不是完全能控的,就必须按能控性分解,只能任意配置可控的极点.(2)极点配置的方法：假设原系统A,B可控,那么采用状态反响

18、阵K,有ABKBH控.设原系统的特征方程为sn,an一-aisa0=0o一010其中K=1"*n,那么有:ABK=-ankn配置后的闭环特征方程为：Sn+an+knsni+a1+kis+a0+k0=0.假设闭环系统希望的极点为'25n得到：f九=s-%s-%s-%=sn+rnsn+r1s+r0o为使系统到达希望性能,比照式1和式2中系数,使之相等,即可求得状态反响阵K=k0,k1,kn/L采用状态反响配置系统极点不改变系统可控性,它不能影响系统中不可控局部模块.2.4.2输出反响法图2-8输出反响示意图对于完全能控的单变量系统,不能采用输出线性反响来实现闭环系统极点任意配置.

19、不能任意配置极点,正是输出线性反响的根本弱点.为了克服这个弱点,在经典限制理论中,往往采取引入附加校正网络,通过增加开环零极点的方法改变根轨迹走向,从而使其落在指定的期望位置上.对于完全能控的单变量系统'A,B,C,通过带动态补偿器的输出反响时限极点任意配置的充要条件是：1.系统完全能观测；2.动态补偿器的阶数为n-1.2.4.2系统极点配置在现代限制理论中是用系统内部的状态来描述系统的,所以经常从系统的状态引出信号作为反响量.利用状态反响只能改变系统能控局部的极点,而不能改变不能控局部的极点,因此利用状态反响进行极点配置的充分必要条件是系统必须是完全能控的.对一个可控系统,在采用状态

20、反响后,可以实现闭环极点的任意配置,即通过状态反响的方法,使闭环系统的极点位于任意期望的位置上.对于其中x是状态变量(n维),u是限制信号,这里选取限制信号为因此T/YA-BKhg系统的稳态响应和瞬态响应特性由矩阵a_bk的特征决定虽然理论上系统的闭环极点离S左半平面越远越好,但是在工业生产实践中,系统极点离左半平面越远,系统的运动状态就变化的越快,这就要求执行机构快速运作,即使再妤的执行元件也会短时间内被损坏掉.所以新的极点的绝对值大约是原系统极点绝对值的3至4倍左右.取P1=-15+40i;P2=-15-40i;P3=-3+10i;P4=-3-10i;利用Matlab进行极点配置,希望可以

21、减小超调量,缩短稳定时间以优化系统.Matlab程序如下：A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;p=-15+40i,-15-40i,-3+10i,-3-10i;k=place(A,B,p)step(A-B*k,B,C,D)运行程序,得到：k=-234.6522131.851214.45616.3957643.3762-89.97656.765836.0878StepResponsex10-From:In(1)From:In(2)earpmA5o2525ol).kuo-DT-1.5cL00.

22、511.5200.511.52Time(sec)505°12too-OT图2-9稳态响应曲线由响应曲线可以看出该系统重新配置极点后,具有较快的调节时间,而且也减少了超调量,改善了系统的动态性能与稳态性能.2.5系统的状态观测器图2-10状态观测器示意图通过状态观测器可以任意配置系统的极点,从而使闭环系统具有期望的稳态和动态性能.但在工业生产中,系统的状态变量并非都是物理量,或者是难以测得的量.这样一来,系统的所有状态变量未必都可以直接测量得到,因此,状态反响这种限制方式在许多实际限制问题中往往难以直接应用和实现.状态观测器就是利用系统的外部输入输出信息来确定系统内部的状态,进而,在系

23、统的极点配置状态反响中,用观测器得到的状态估计值代替系统的真实状态.下列图为状态观测器的结构图：图2-11状态观测器示意图使用MATLAB为本系统设置状态观测器,选用极点配置时的极点,程序如下列图所示：A=0010;0001;-400300-96;150-2003-4.5;C=1000；0100；D=00;00;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=

24、ss(A2,B,C,D)sys2=ss(A3,B,C,D)运行上面程序,得到:L1=7.083330.0895-30.579615.4167-41.6552-96.5401168.1877200.0790A2=-7.0833-30.08951.000030.5796-15.4167-358.3448396.5401-9.0000-18.1877-400.07903.000001.00006.0000-4.5000L2=3.7432-7.1200-21.4563-3.7432190.989493.5822115.5037-24.2083A3=21.45633.74321.0000-655.579

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

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

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

28、使其对连续信号中的最高频率分量,在一个周期内a被采样2次以上(上半周与下半周都至少采样一次),那么采样后的脉冲序列中将包含了连续信号的全部信息.但是,在仿真中所遇到的大多数被再现信号是没有频带限的,所以一般取采样频率再现信号主要频带中的最高频率的510倍.在离散限制系统的设计过程中,采样周期确实定依据的是现场检测的被调量信号的频率,对于频率较高的信号,采样周期的设定就小,而对于变化过程较慢的低频信号,采样周期的设定可以大一些.有关概念在工程上的实际应用会有专门的内容介绍.线性连续系统状态方程离散化的实质是将矩阵微分方程化为矩阵差分方程,它是描述多输入多输出离散系统的一种方便的数学模型.在推导离

29、散化系统的方程时,假定系统是周期性采样,并且采样脉冲宽度远小于采样周期,采样周期T的选择满足香农采样定理,还假设系统具有零阶保持特性,即在两个采样瞬间之间,采样值不变,并等于前一个采样时刻的值.通常离散化的方法有很多,伤J如欧拉法,梯形法,龙哥-库塔(Runge-Kutta)法,阿达姆斯(Adams)法等等.下面主要运用三种方法来对系统进行离散化并运用计算机进行模拟系统的特性,分析不同采样周期对系统的影响效果.2.6.2 零阶保持器零阶保持器可以将脉冲序列变成连续的方波信号,即将前一个采样周期的数值保存到下一个采样点到来的时候.在Matlab中输入函数如下：A=0010;0001;-40030

30、0-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;p=-15+40i,-15-40i,-3+10i,-3-10i;k=place(A,B,p);H,I,J,K=c2dm(A-B*k,B,C,D,0.1,'zoh')dstep(H,I,J,K)分别设置采样时间为0.1s,0.05s,0.01s,运行程序,得到下列图:t=0.1st=0.05sStepResponset=0.01s图2-12零阶保持器离散化2.6.3 一阶保持器采用一阶保持器进行离散化,程序如下：A=0010;0001;-400300-96;150-2003-4.5;B=00;00;10;00.5;C=1000;0100;D=00;00;p=-15+40i,-15-40i,-3+10i,-3-10i;k=place