资源描述：
《现代控制工程课程设计》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、.现代控制工程课程设计目录一、设计目的及要求1设计目的1设计要求1二、状态空间方程建立1三、倒置摆的状态空间模型2四、模型特性分析31.可控性分析32.可观性分析33.稳定性分析44.MATLAB仿真求解4五、倒置摆系统的综合及仿真分析5五、全维状态观测器的设计及仿真分析8六、参考文献13一、设计目的及要求设计目的1、通过课程报告加深理解现代控制理论中的一些基本概念；2、掌握用状态方程描述的线性系统的稳定性、能控性、能观性的分析计算方法；3、掌握状态反馈系统的综合和全维观测器设计方法4、重视理论计算，MATLAB编程计算及SIM
2、ULINK仿真能力，提高计算机编程计算能力。设计要求如图1所示，为单倒置摆系统的原理图。设摆的长度为、质量为m，用铰链安装在质量为M的小车上。小车有一台直流电动机拖动，在水平方向对小车施加控制力u,相对参考系产生位移z。若不给小车施加控制力，则倒置摆会向左或向右倾倒，因此，它是一个不稳定系统。控制的目的是，当倒置摆无论出现向左或向右倾倒时，通过控制直流电动机，使小车在水平方向运动，将倒置摆保持在垂直位置上。二、状态空间方程建立为简化问题，工程上可以忽略一些次要因素。在本例中，我们为了简化问题，方便研究系统空间的设计问题，忽略了摆
3、杆质量、执行电动机惯性以及摆轴、轮轴、轮与接触面之间的摩擦及风力。设小车的瞬时位置为z，倒置摆出现的偏角为θ，则摆心瞬时位置为13。在控制力u的作用下，小车及摆均产生加速运动，根据牛顿第二定律，在水平直线运动方向的惯性力应与控制力u平衡，则有即(1)由于绕摆旋转运动的惯性力矩应与重力矩平衡，因而有即（2）式（1）、式（2）两个方程都是非线性方程，需作线性化处理。由于控制的目的是保持倒置摆直立，因此，在施加合适u的条件下，可以认为、均接近于零，此时≈，,且可以忽略项，于是有（3）+=（4）连联立求解式（3）、式（4），可得（5）（
4、6）消去中间变量θ，可得输入量为u、输出量为z的微分方程为（7）综合上述的分析，可抽象出系统的研究对象为：位移z、小车的速度、摆的角速度θ及其角速度的。系统的研究对象抽象成这四个变量后，接下来就可以根据前面的方程为这四个变量建立空间状态方程，并分析被控对象的特性。三、倒置摆的状态空间模型在上一章节中，我们已经选取了四个研究对象作为状态变量，它们分别为：位移z、小车的速度、摆的角速度θ及其角速度的。Z为输出变量，在考虑，以及式（5）、（6）、（7），可列出倒置摆的状态空间模型表达式为：（8a）(8b)13式中为方便研究，假定系统的
5、参数M=1kg,m=0.1kg,l=1m,,则系统状态方程中参数矩阵为：,，（9）此时倒置摆的状态空间模型表达式为：四、模型特性分析在建立完模型后我们需要对模型进行分析。作为被控制的倒置摆，当它向左或向右倾倒时，能否通过控制作用使它回复到原直立位置，这取决于其能控性。因此我们首先分析它的能控性。1.可控性分析能控性判据：线性定常系统，完全能控的充分必要的条件是能控性矩阵的秩为，即，并将式（9）的有关数据带入该判据，可得（10）因此，单倒置摆的运动状态是可控的。换句话说，这意味着总存在一控制作用u,将非零状态转移到零。2.可观性分
6、析能观测性判据：对于线性定常系统，系统完全能观测的充分必要条件是能观测性矩阵的秩为，即133.稳定性分析由单倒置摆系统的状态方程，可求的其特征方程为：（11）解得特征值为0,0，，-。四个特征值中存在一个正根，两个零根，这说明单倒置摆系统，即被控系统不稳定的。由线性系统用李雅普诺夫稳定性判据时，系统矩阵A必须是非奇异的，系统仅存在唯一的平衡状态，由于该系统矩阵A是奇异的，故也可得出系统是不稳定的。4.MATLAB仿真求解经过上述理论求解后，我们可以得出的MATLAB代码如下：A=[0,1,0,0;0,0,-1,0;0,0,0,1
7、;0,0,11,0];B=[0;1;0;-1];C=[1,0,0,0];D=0;Qc=ctrb(A,B)nc=rank(Qc);ifnc<3disp('系统不可控')elsedisp('系统可控')endQo=obsv(A,C)no=rank(Qo);ifno<3disp('系统不可观')elsedisp('系统可观')end运行结果如下：Qc=010110100-10-11-10-110系统可控Qo=1000010000-10000-1系统可观由上面对系统模型进行分析的结果可知被控系统是具有可控性和可观性，但是被控系统是不稳定
8、的，需对被控系统进行反馈综合，使四个特征值全部位于根平面S左半平面的适当位置，以满足系统的稳定工作已达到良好、静态性能的要求。因此我们需要设计全维状态观测器来使系统到达控制的目的。13五、倒置摆系统的综合及仿真分析采用全状态反馈。取状态变量z、、θ、为反馈信号，
