倒立摆数学模型

1单级倒立摆的数学模型的建立：小车由电机通过同步带驱动在滑杆上来回运动，保持摆杆平衡。电机编码器和角编码器向运动卡反馈小车和摆杆位置（线位移和角位移）。导轨截面成H型，小车在轨道上可以自由滑动，其在轨道上的有效运行长度为1米。轨道两端装有电气限位开关，以防止因意外失控而撞坏机构。图1单级倒立摆系统数学模型倒立摆系统的模型参数如下[]：M小车质量1.096Kg；m摆杆质量0.109Kgb小车摩擦系数0.1N/m/secI摆杆质量0.0034kg*m*ml摆杆转动轴心到杆质心的长度0.25mT采样频率0.005s下面N和P为小车与摆杆相互作用力的水平和垂直方向的分量。分析小车水平方向所受的合力，可得到方程为：NxbFxM(1)由摆杆水平方向的受力进行分析可以得到下面等式：sincossin222mlmlxmNlxdtdmN（2）把这个等式代入（1）式中，得到系统的第一个运动方程：FmlmlxbxmMsincos2（3）为了推出系统的第二个运动方程，对摆杆垂直方向的合力进行分析，得到下面的方程：cos22ldtdmmgPcossin2mlmlmgP（4）力矩平衡方程如下：INlPlcossin（5）方程中力矩的方向，由于，sinsin,coscos，故等式前面有负号。合并这两个方程，约去P和N，得到第二个运动方程：cossin2xmlmglmlI（6）假设与1（单位是弧度）相比很小，即1，则可进行近似处理：0,sin,1cos2dtd用u代表被控对象的输入力，线性化后两个运动方程如下：umlxbxmMxmlmglmlI2（7）对方程（7）进行拉普拉斯变换，得到：)()()()()()()(22222sUssmlssbXssXmMssmlXsmglssmlI（8）（推到时假设初始条件为0）则，摆杆角度和小车位移的传递函数为：mglsmlImlssXs222)()()(将上述参数代入，摆杆角度和小车位移的传递函数为：26705.00102125.002725.0)()(22sssXs摆杆角度和小车加速度之间的传递函数为：mglsmlImlsAs22)()(将上述参数代入，摆杆角度和小车加速度之间的传递函数为：26705.00102125.002725.0)()(22sssAs摆杆角度和小车所受外界作用力的传递函数：22432222()()()()()()mlssqbImlMmmglbmglFsssssqqqqMmImlml将上述参数代入，摆杆角度和小车所受外界作用力的传递函数：32()2.35655()0.088316727.91692.30942ssFssss以外界作用力作为输入的系统状态空间表达式为：222222222201000()00()()()00010()00()()()xxImlbmglImlxxIMmMmlIMmMmlIMmMmlumlbmglMmmlIMmMmlIMmMmlIMmMml1000000100xxxyu将上述参数代入，以外界作用力作为输入的系统状态空间表达式为：0100000.08831670.62931700.8831670001000.23565527.828502.356551000000100xxxxuxxxyu以小车加速度作为输入的系统系统状态空间表达式：'0100000001000103300044xxxxugll'1000000100xxxyu将上述参数代入，以小车加速度作为输入的系统系统状态空间表达式：0100000001000100029.4031000000100xxxxuxxxyu2系统的可控性、可观测性分析对于连续时间系统：BuAXXDuCXy系统状态完全可控的条件为：当且仅当向量组BAABBn1,...,,是线性无关的，或n×n维矩阵BAABBn1的秩为n。系统的输出可控条件为：当且仅当矩阵DBCABCACABCBn12的秩等于输出向量y的维数。应用以上原理对输入为加速度输出为摆杆与竖直方向的角度的夹角时的系统进行可控性分析：04.2900100000000010A3010B01000001C00D带入上式，计算得：2301001000:::00010029.40SBABABAB4rankS23:::VCCACACA带入数值得：4rankV令0SIA得系统的开环特征方程为（0，0，5.42，-5.42）系统状态可控性矩阵的秩=4=系统的状态变量的维数，系统的输出完全可控性矩阵的秩=2=系统输出向量y的维数，所以系统可控。可观测性矩阵的秩=4=矩阵A的维数，所以系统可观测。系统有一个极点位于s又半平面上，有两个极点位于坐标原点。所以系统不稳定。因此可以对系统进行控制器的设计，使系统稳定。2.3系统阶跃响应分析上面已经得到系统的状态方程，先对其进行阶跃性分析，在Matlab中键入以下指令：clear;A=[0100;0000;0001;0029.40];B=[0103]';C=[1000;0100];D=[00]';step(A,B,C,D)得到如下图：可以看出，在单位阶跃响应作用下，小车位置和摆杆角度都是发散的。因此接下来，我要进行控制器的设计。