线性系统实验报告2015.doc

《线性系统实验报告2015.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、线性系统实验报告姓名：院系：学号：导师：2015年12月一、概述日常生活中所见到的任何重心在上、支点在下的杆都可以称之为倒立摆。从工程背景来讲，一级倒立摆的背景源于火箭发射助推器，二级倒立摆与双足机器人控制有关，三级倒立摆应当说由一、二级倒立摆演绎而来，背景相当复杂。倒立摆仿真或实物控制实验是控制领域中用来检验某种控制理论或方法的典型方案。倒立摆研究的意义是，作为一个实验装置，它形象直观，简单，而且参数和形状易于改变；但它又是一个高阶次、多变量、非线性、强耦合、不确定的绝对不稳定系统的被控系统，必须采用十分有效的控制手段才能使之稳定。因此，许多新的控制理论，都通过倒立摆试验对理论加

2、以实物验证，然后在应用到实际工程中去。因此，倒立摆成为控制理论中经久不衰的研究课题，是验证各种控制算法的一个优秀平台，故通过设计倒立摆的控制器，可以对控制学科中的控制理论有一个学习和实践机会。二、实验目的1.建立二级倒立摆系统的数学模型对于倒立摆系统，由于其本身是自不稳定的系统，实验建立模型存在一定的困难，但是经过小心的假设忽略掉一些次要的因素后，倒立摆系统就是一个典型的运动的刚体系统，可以在惯性坐标系内应用经典力学理论建立系统的动力学方程。对于直线二级倒立摆，由于其复杂程度，在这里利用拉格朗日方程推导运动学方程。2.模型线性化由于模型的动力学方程中存在三角函数，因此方程是非线性的

3、，通过小角度线性化处理，将动力学非线性方程变成线性方程，便于后续的工作的进行。3.控制器设计利用Matlab对已经建立好模型的系统进行分析。分析内容主要包括得出原未加控制器时系统的极点分布，系统的能观性，能控性；然后根据上步的结论采用相应的控制策略，设计具体的控制器，加入到系统中，然后通过仿真调节具体参数，得到最优的控制器参数。4.实物仿真通过实验室提供的倒立摆实物以及控制计算机，将设计的控制算法加入到控制器中，通过倒立摆实物来验证设计的控制算法的正确性，并在控制器能够有效控制系统的情况下，改变控制器参数，分析系统控制特性的变化，并且得到系统实时的数据采集，画出图形。三、实验步骤1

4、.根据建模结果仔细计算并寻找合适的理论控制器参数。2.进入matlabcommand窗口，键入仿真文件，进行仿真实验。通过调节参数请仔细观察思考控制器参数对系统瞬态响应和稳态响应的影响。找到几组合适的控制器参数作为实际控制的参数。3.倒立摆系统实物调试具体步骤如下：（1）将小车推到导轨的中央，并且使摆杆处于自由下垂的静止状态，（2）给小车的驱动器的电源接通电源，进入MATLAB的SIUMLINK中，（3）调出实物仿真文件，填入设计的控制器参数，编译，连接实物，双手将倒立摆竖直的缓慢的立起，倒立摆进入到与竖直方向夹角小于10度的范围内，控制器启动，观察控制器的控制效果。（4）将保存的

5、实验数据画出图形，分析实际控制效果与仿真曲线的区别。四、实验报告1.倒立摆系统模型在忽略了空气流动，各种摩擦之后，可将倒立摆系统抽象成小车、匀质杆和质量块组成的系统，如图所示。直线两级倒立摆系统参数：M小车质量1.32Kg摆杆1质量0.04Kg摆杆2质量0.132Kg质量块的质量0.208Kg摆杆1与垂直向上方向的夹角摆杆2与垂直向上方向的夹角摆杆1转动中心到杆质心的距离0.09m摆杆2转动中心到杆质心的距离0.27m作用在系统上的外力利用拉格朗日方程推导运动学方程为：2.写出系统在（0，0，0，0，0，0，0）和（0，pi,0,0,0,0,0）两种平衡状态下的进行泰勒级数展开并线

6、性化得到的状态矩阵。将（3），（4）式对求解代数方程，得到以下两式（5）(6)表示成以下形式：(7)(8)取平衡位置时各变量的初值为零，(9)将（7）式在平衡位置进行泰勒级数展开，并线性化，令:(10)(11)(12)(13)(14)(15)(16)带入（8）式，得到线性化之后的公式(17)将（8）式在平衡位置进行泰勒级数展开，并线性化，令(18)(19)(20)(21)(22)(23)(24)带入（20）式，得到(25)即：(26)(27)现在得到了两个线性微分方程，由于我们采用加速度作为输入，因此还需加上一个方程(28)取状态变量如下：(29)由（26），（27），（29）式得

7、到状态空间方程如下：(30)3.分析这两种状态下系统的开环特性并设计两种状态下的控制器。由以上方程，将以下参数代入(31)求出各个值：(32)得到状态矩阵为：A=0001.00000000001.00000000001.0000000000086.6907-21.61720000-40.311239.4500000B=0001.00006.6402-0.0877C=100000010000001000000100000010000001D=[000000]'4.给出理