基于matlab的倒立摆模糊控制课程设计报告

智能控制理论及应用课程设计报告题目基于MATLAB的倒立摆模糊控制院系西北民族大学电气工程学院专业班级10级自动化（3）班学生姓名蔡余敏学号P101813455指导教师刁晨201310基于MATLAB的倒立摆模糊控制作者蔡余敏指导老师刁晨摘要倒立摆的控制问题就是使摆杆尽快地达到一个平衡位置，并且使之没有大的振荡和过大的角度和速度。当摆杆到达期望的位置后，系统能克服随机扰动而保持稳定的位置。本文主要针对较为简单的单级倒立摆控制系统而进行的设计分析。通过建立微分方程模型，求出相关参数，设计出对应的模糊控制器，并运用MATLAB软件进行系统模型的软件仿真，从而达到预定控制效果。目前，一级倒立摆的研究成果应用于火箭发射推进器和控制卫星的飞行状态等航空航天领域。关键词单级倒立摆；微分方程；模糊控制；MATLAB仿真1背景分析倒立摆控制系统是一个复杂的、不稳定的、非线性系统，是进行控制理论教学及开展各种控制实验的理想实验平台。对倒立摆系统的研究能有效的反映控制中的许多典型问题如非线性问题、鲁棒性问题、镇定问题、随动问题以及跟踪问题等。通过对倒立摆的控制，用来检验新的控制方法是否有较强的处理非线性和不稳定性问题的能力。同时，其控制方法在军工、航天、机器人和一般工业过程领域中都有着广泛的用途，如机器人行走过程中的平衡控制、火箭发射中的垂直度控制和卫星飞行中的姿态控制等。正是由于倒立摆系统的特殊性，许多不同领域的专家学者在检验新提出理论的正确性和实际可行性时，都将倒立摆系统作为实验测试平台。再将经过测试后的控制理论和控制方法应用到更为广泛的领域中去。现代控制理论已经在工业生产过程、军事科学、航空航天等许多方面都取得了成功的应用。例如极小值原理可以用来解决某些最优控制问题；利用卡尔曼滤波器可以对具有有色噪声的系统进行状态估计；预测控制理论可以对大滞后过程进行有效的控制。但是它们都有一个基本的要求需要建立被控对象的精确数学模型。随着科学技术的迅猛发展，各个领域对自动控制控制精度、响应速度、系统稳定性与适应能力的要求越来越高，所研究的系统也日益复杂多变。然而由于一系列的原因，诸如被控对象或过程的非线性、时变性、多参数间的强烈耦合、较大的随机干扰、过程机理错综复杂、各种不确定性以及现场测量手段不完善等，难以建立被控对象的精确模型。虽然常规自适应控制技术可以解决一些问题，但范围是有限的。对于像二级倒立摆这样的非线性、多参数、强耦合的被控对象，使用传统控制理论难以达到良好的控制性能。而模糊控制理论能够克服这些困难，达到实际设计要求。所以说，对倒立摆系统控制理论的研究不仅具有理论研究价值，也具有相当的实际工程应用价值。2提出控制问题倒立摆控制问题是展示智能控制方法由于传统控制方法的典型范例。一级倒立摆的背景源于火箭发射助推器；二级倒立摆于双足机器人控制有关。这里只讨论一级倒立摆的控制问题。有一个倒立摆控制系统如图21所示。图21倒立摆控制系统它由小车和倒立摆构成，小车在控制器的作用下，沿滑轨在水平方向运动，使倒立摆在垂直平面内稳定。倒立摆的控制问题就是使摆杆尽快地达到一个平衡位置，并且使之没有大的振荡和过大的角度和速度。当摆杆到达期望的位置后，系统能克服随机扰动而保持稳定的位置。倒立摆系统的输入为小车的位移（即位置）和摆杆的倾斜角度期望值，计算机在每一个采样周期中采集来自传感器的小车与摆杆的实际位置信号，与期望值进行比较后，通过控制算法得到控制量，再经数模转换驱动直流电机实现倒立摆的实时控制。直流电机通过皮带带动小车在固定的轨道上运动，摆杆的一端安装在小车上，能以此点为轴心使摆杆能在垂直的平面上自由地摆动。作用力平行于铁轨的方向作用于小车，使杆绕小车上的轴在竖直平面内旋转，小车沿着水平铁轨运动。当没有作用力时，摆杆处于垂直的稳定的平衡位置（竖直向下）。为了使杆子摆动或者达到竖直向上的稳定，需要给小车一个控制力，使其在轨道上被往前或朝后拉动。3建立控制模型首先假设摆杆为刚体；忽略摆杆与支点之间的摩擦；忽略小车与导轨之间的摩擦。一级倒立摆系统可抽象成小车与匀质杆组成的系统，假设M为小车的质量；M为摆杆质量；为摆杆转动轴心到杆质心的长度；LI为摆杆惯量；U为加在小车上的力；为小车位置；为摆杆与X垂直向上方向的夹角。应用NEWTON第二定律的方法可得到系统方向的运动方程为X21MAFMGXXX13COSSINI2222DTDTMLDTXLTT规定逆时针方向的力矩为正，以摆与小车的连接点为原点，列出摆的力矩方程考虑到摆的惯性力矩，求得系统的运动方程为（未考虑摆旋转的摩擦阻力矩）23COSSIN3434SINCOS222222DTXGDTLMLJGLLDTXT由式（31）和式（32）可得3COSSIN34COS2SIN1SIN22222MMDTLDTLFDTXLMMFDTMLGDT4设计模糊控制器41确定输入输出变量以摆角、摆角角速度、小车位移、速度为状态变量。XX将这些状态变量作为控制器输入量，以作用在小车的力F作为模糊控制器输出量。所以分别在四个输入变量的空间建立相应的隶属度函数。图41四输入变量单级倒立摆以下分别在MATLAB的仿真图中来说明情况。确定摆角的论域03，03，将其划分为两个语言变量“大”和“小”，隶属度函数如图42；摆角角速度的论域1,1，划分两个语言变量为“快”和“慢”；隶属度函数如图43；小车位移的论域3，3，划分两个语言变量为“远”和“近”；X隶属度函数如图44；速度的论域3，3，划分两个语言变量X为“快”和“慢”；隶属度函数如图45；输出变量的论域为10，10，如图46。图42摆角的隶属度函数图43摆角角速度的隶属度函数图44位移的隶属度函数X图45速度的隶属度函数X图46输出变量42设计模糊规则库这里选取TS的控制器，控制器根据这4个输入变量，综合得出作用于小车的控制信号。然后，列出每种输入所对应的输出量的模糊规则，共计设置了16条规则1IFIN1ISIN1MF1ANDIN2ISIN2MF1ANDIN3ISIN3MF1ANDIN4ISIN4MF1THENOUTISMF112IFIN1ISIN1MF1ANDIN2ISIN2MF1ANDIN3ISIN3MF1ANDIN4ISIN4MF2THENOUTISMF213IFIN1ISIN1MF1ANDIN2ISIN2MF1ANDIN3ISIN3MF2ANDIN4ISIN4MF1THENOUTISMF314IFIN1ISIN1MF1ANDIN2ISIN2MF1ANDIN3ISIN3MF2ANDIN4ISIN4MF2THENOUTISMF415IFIN1ISIN1MF1ANDIN2ISIN2MF2ANDIN3ISIN3MF1ANDIN4ISIN4MF1THENOUTISMF516IFIN1ISIN1MF1ANDIN2ISIN2MF2ANDIN3ISIN3MF1ANDIN4ISIN4MF2THENOUTISMF617IFIN1ISIN1MF1ANDIN2ISIN2MF2ANDIN3ISIN3MF2ANDIN4ISIN4MF1THENOUTISMF718IFIN1ISIN1MF1ANDIN2ISIN2MF2ANDIN3ISIN3MF2ANDIN4ISIN4MF2THENOUTISMF819IFIN1ISIN1MF2ANDIN2ISIN2MF1ANDIN3ISIN3MF1ANDIN4ISIN4MF1THENOUTISMF9110IFIN1ISIN1MF2ANDIN2ISIN2MF1ANDIN3ISIN3MF1ANDIN4ISIN4MF2THENOUTISMF10111IFIN1ISIN1MF2ANDIN2ISIN2MF1ANDIN3ISIN3MF2ANDIN4ISIN4MF1THENOUTISMF11112IFIN1ISIN1MF2ANDIN2ISIN2MF1ANDIN3ISIN3MF2ANDIN4ISIN4MF2THENOUTISMF12113IFIN1ISIN1MF2ANDIN2ISIN2MF2ANDIN3ISIN3MF1ANDIN4ISIN4MF1THENOUTISMF13114IFIN1ISIN1MF2ANDIN2ISIN2MF2ANDIN3ISIN3MF1ANDIN4ISIN4MF2THENOUTISMF14115IFIN1ISIN1MF2ANDIN2ISIN2MF2ANDIN3ISIN3MF2ANDIN4ISIN4MF1THENOUTISMF15116IFIN1ISIN1MF2ANDIN2ISIN2MF2ANDIN3ISIN3MF2ANDIN4ISIN4MF2THENOUTISMF161其中，IN1MF1和IN1MF2表示角是“大”和“小”，IN2MF1和IN2MF2表示角速度是“快”和“慢”，IN3MF1和IN3MF2表示位移的X“远”和“近”、IN4MF1和IN4MF2表示速度是“快”和“慢”。X在多输入多输出的模糊逻辑系统中，常可能会有很多规则，并且这些规则的前提部分和结论部分也可以由许多语句组成。前面介绍的IFTHEN规则只是最简单的形式，是最基本的模糊系统单元。对于复杂的模糊逻辑系统，如果满足“输出一输入解耦的，各输入间、输出间分别是独立的”，就可以用连接词“AND”、“OR”和“ALSO”连接的一系列简单的模糊规则组成的模糊规则库来表示。其中，“AND”和“OR”用来连接同一规则的多个输入或多个输出，“ALSO”用来连接多条不同的模糊规则。模糊推理的规则通常来源于专家的知识，对于多输入多输出系统，其规则库可以看成由N个子规则库所组成，每一个子规则库由N个多输入单输出规则所组成，而且每个子规则库是相互独立的，因此通常只要考虑一个多输入单输出规则库的模糊推理问题，分别求出每一个单输出，组合起来就是最终结果。然后将上述实施规则录入到MATLAB仿真结构中，具体如下图47所示。图47模糊规则录入结束后要对所指定的输入输出规则进行保存，以便在具体的仿真过程中实现调用，保存的文件扩展名为单级倒立摆CAYMIFIS。图48模糊推理规则通过VIEW直观显示图此外，在建立起具体的模糊推理规则以后，可以通过VIEW功能来得到规则的直观显示，这些图形可以帮助更好的理解建立的模糊推理规则。如图48所示。设置完后保存为单级倒立摆CAYMIFIS文件。对应的源代码如下SYSTEMNAME单级倒立摆CAYMITYPESUGENOVERSION20NUMINPUTS4NUMOUTPUTS1NUMRULES16ANDMETHODPRODORMETHODPROBORIMPMETHODPRODAGGMETHODSUMDEFUZZMETHODWTAVERINPUT1NAMEIN1RANGE0303NUMMFS2MF1IN1MF1GBELLMF,032503MF2IN1MF2GBELLMF,032503INPUT2NAMEIN2RANGE11NUMMFS2MF1IN2MF1GBELLMF,1251MF2IN2MF2GBELLMF,1251INPUT3NAMEIN3RANGE33NUMMFS2MF1IN3MF1GBELLMF,3253MF2IN3MF2GBELLMF,3253INPUT4NAMEIN4RANGE33NUMMFS2MF1IN4MF1GBELLMF,3253MF2IN4MF2GBELLMF,3253OUTPUT1NAMEOUTRANGE01NUMMFS16MF1MF1CONSTANT,0MF2MF2CONSTANT,0125MF3MF3CONSTANT,025MF4MF4CONSTANT,0375MF5MF5CONSTANT,05MF6MF6CONSTANT,0625MF7MF7CONSTANT,075MF8MF8CONSTANT,0875MF9MF9CONSTANT,1MF10MF10CONSTANT,0MF11MF11CONSTANT,0166666666666667MF12MF12CONSTANT,0333333333333333MF13MF13CONSTANT,05MF14MF14CONSTANT,0666666666666667MF15MF15CONSTANT,0833333333333333MF16MF16CONSTANT,1RULES1111,1111112,2111121,3111122,4111211,5111212,6111221,7111222,8112111,9112112,10112121,11112122,12112211,13112212,14112221,15112222,1611根据控制器和倒立摆模型的分析，通过MATLAB可以得到单级倒立摆控制系统的结构如图所示。图49单级倒立摆控制系统的结构5根据控制模型，分析动稳态特性根据控制器和倒立摆模型得出图49所示的单级倒立摆控制系统的结构后，设置不同的隶属度函数，分别对上述过程进行仿真，得到不同情况下的线速度和角位移的变化情况，说明选取不同的隶属度函数曲线在控制效果上有着明显的差异，其中系统默认的隶属度函数是较好的控制方法。通过不断调试，达到最佳状态，在MATLAB中进行仿真，设定方波为跟踪的期望信号，得到结果如图51所示。图52所示为仿真运行时的图形用户界面。图51小车对方波期望值的跟踪情况图52仿真运行时的图形用户界面从以上的图51可知，绿色规则的方波期望值代表小车位置的期望，红色不规则的近似方波为小车的实际位置。从以上的图52可知，绿色的三角形为小车位置的期望，红色的方块为小车的实际位置，黄色的粗斜线为摆杆实际位置。图53小车在不同时刻跟随情况从以上的图53的四个不同情况可知，当小车的实际位置和期望位置相距最远时，跟踪误差达到最大值；而当小车的实际位置和期望位置相距最近时，跟踪误差达趋于零；尤其是在小车的实际位置和期望位置重合，并且摆杆垂直时，跟踪误差为零。图53所示不同时刻小车和摆杆的运动轨迹图。图53不同时刻小车和摆杆的运动轨迹6分析数据，得出结论通过以上仿真实践的过程可以看出，对于倒立摆一类的包含非线性以及不稳定因素的控制系统模型，采用模糊控制的方法具有很强的适应性以及理想的控制效果，在此基础上，选用恰当的隶属度函数也是实现优化控制的重要条件。倒立摆系统的输入