倒立摆课程设计(共16页).doc

《倒立摆课程设计(共16页).doc》由会员分享，可在线阅读，更多相关《倒立摆课程设计(共16页).doc（16页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上 控制系统设计课程设计报告题目：倒立摆控制系统设计 班级： 姓名: 学号: 2011年03月09日课程设计（论文）任务书专 业自动化班 级学 生指导教师题 目倒立摆系统的PID控制设计时间2011年 3 月 2 日 至 2011年3月13日 共 2 周设计要求设计（论文）的任务和基本要求，包括设计任务、查阅文献、方案设计、说明书（计算、图纸、撰写内容及规范等）、工作量等内容。1建立一阶倒立摆数学模型2做模型仿真试验（1）给出Matlab仿真程序。（2）给出仿真结果和响应曲线。3倒立摆系统的PID控制算法设计设计PID控制器，使得当在小车上施加1N的脉冲信号时，闭环系统

2、的响应指标为：（1）稳定时间小于5秒（2）稳态时摆杆与垂直方向的夹角变化小于0.1 弧度并作PID控制算法的MATLAB仿真4倒立摆系统的最优控制算法设计用状态空间法设计控制器，使得当在小车上施加0.2m的阶跃信号时，闭环系统的响应指标为：（1）摆杆角度和小车位移的稳定时间小于5秒（2）的上升时间小于1秒（3）的超调量小于20度（0.35弧度）（4）稳态误差小于2%。指导教师签字： 系（教研室）主任签字：2011年 3月 10 日1. 实验设备简介本课程设计的被控对象采用固高公司生产的GIP-100-L型一阶倒立摆系统，课程设计包括三方面的内容：（1）建立一阶倒立摆的线性化数学模型；（2）倒立

3、摆系统的PID控制器设计、MATLAB仿真及实物调试；（3）倒立摆系统的最优控制器设计、MATLAB仿真及实物调试。一阶倒立摆系统的结构示意图如图1所示。 图1 一阶倒立摆结构示意图系统组成框图如图2所示。 图2 一阶倒立摆系统组成框图系统是由计算机、运动控制卡、伺服机构、倒立摆本体和光电码盘几大部分组成的闭环系统。光电码盘1将小车的位移、速度信号反馈给伺服驱动器和运动控制卡，摆杆的角度、角速度信号由光电码盘2反馈给运动控制卡。计算机从运动控制卡中读取实时数据，确定控制决策（小车运动方向、移动速度、加速度等），并由运动控制卡来实现该控制决策，产生相应的控制量，使电机转动，通过皮带，带动小车运动

4、，保持摆杆平衡。2.建立一阶倒立摆数学模型对系统建立数学模型是系统分析、设计的前提，而一个准确又简练的数学模型将大大简化后期的工作。为了简化系统分析，在实际的模型建立过程中，要忽略空气流动阻力，以及各种次要的摩擦阻力。这样，可将倒立摆系统抽象成小车和匀质刚性杆组成的系统，如下图所示。本系统内部各相关参数定义如下： 小车质量 摆杆质量 小车摩擦系数 摆杆转动轴心到杆质心的长度 摆杆惯量 加在小车上的力 小车位置 摆杆与垂直向上方向的夹角 摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）下图是系统中小车和摆杆的受力分析图。其中，和为小车与摆杆相互作用力的水平和垂直方向的分量。注意：在实际倒

5、立摆系统中检测和执行装置的正负方向已经完全确定，因而矢量方向定义如图，图示方向为矢量正方向。应用Newton方法来建立系统的动力学方程过程如下：分析小车水平方向所受的合力，可以得到以下方程：由摆杆水平方向的受力进行分析可以得到下面等式：即：把这个等式代入上式中，就得到系统的第一个运动方程： （2-1）为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程：即：力矩平衡方程如下：注意：此方程中力矩的方向，由于，故等式前面有负号。合并这两个方程，约去和，得到第二个运动方程： （2-2）1微分方程模型设，当摆杆与垂直向上方向之间的夹角与1（单位是弧度）相比很小，即 时，则

6、可以进行近似处理：，。为了与控制理论的表达习惯相统一，即一般表示控制量，用来代表被控对象的输入力，线性化后得到该系统数学模型的微分方程表达式： (2-3) 2传递函数模型对方程组（2-3）进行拉普拉斯变换，得到 （2-4）注意：推导传递函数时假设初始条件为0。由于输出为角度，求解方程组（2-4）的第一个方程，可以得到把上式代入方程组（2-4）的第二个方程，得到整理后得到以输入力为输入量，以摆杆摆角为输出量的传递函数：其中 若取小车位移为输出量，可得传递函数：其中 q=(M+m)(I+ml2)-(ml)2实际系统参数为M=1.096kgm=0.109kgb=0.1N/m/secl=0.25mI=

7、0.0034kg*m23.PID控制原理在模拟控制系统中，控制器最常用的控制规律是PID控制。常规PID控制系统原理框图如下图所示。系统由模拟PID控制器KD(S)和被控对象G(S)组成。PID控制器是一种线性控制器，它根据给定值与实际输出值构成控制偏差将偏差的比例(P)、积分(I)和微分(D)通过线性组合构成控制量，对被控对象进行控制，故称PID控制器。其控制规律为或写成传递函数的形式式中：比例系数；积分时间常数；微分时间常数。在控制系统设计和仿真中，也将传递函数写成式中：比例系数；积分系数；微分系数。简单说来，PID控制器各校正环节的作用如下：（1）比例环节：成比例地反映控制系统的偏差信号

8、，偏差一旦产生，控制器立即产生控制作用，以减少偏差。（2）积分环节：主要用于消除稳态误差，提高系统的型别。积分作用的强弱取决于积分时间常数，越大，积分作用越弱，反之则越强。（3）微分环节：反映偏差信号的变化趋势（变化速率），并能在偏差信号值变得太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减小调节时间。2摆杆角度控制这个控制问题和我们以前遇到的标准控制问题有些不同，在这里输出量为摆杆的位置，它的初始位置为垂直向上，我们给系统施加一个扰动，观察摆杆的响应。系统框图如下：图中是控制器传递函数，是被控对象传递函数。考虑到输入，结构图可以很容易地变换成该系统的输出为其中： 被控

9、对象传递函数的分子项被控对象传递函数的分母项PID控制器传递函数的分子项PID控制器传递函数的分母项被控对象的传递函数是其中 PID控制器的传递函数为只需调节PID控制器的参数，就可以得到满意的控制效果。3小车位置控制小车位置作为输出时，系统框图如下：其中，是摆杆传递函数，是小车传递函数。由于输入信号，所以可以把结构图转换成：其中，反馈环代表我们前面设计的摆杆的控制器。从此框图我们可以看出此处只对摆杆角度进行了控制，并没有对小车位置进行控制。小车位置输出为：其中，分别代表被控对象1和被控对象2传递函数的分子和分母。和代表PID控制器传递函数的分子和分母。下面我们来求，根据前面实验二的推导，有可

10、以推出小车位置的传递函数为其中 可以看出， =，小车的算式可以简化成： 4模型仿真实验MATLAB仿真程序如下:%-trans.m-%-倒立摆传递函数，开环极点及开环脉冲响应-%-输入倒立摆传递函数G（S）=num/den-M=1.096;m=0.109;b=0.1;I=0.0034;g=9.8;l=0.25;q=(M+m)*(I+m*l2)-(m*l)2;%-计算并显示多项式形式的传递函数-num=m*l/q 0 0den=1 b*(I+m*l2)/q -(M+m)*m*g*l/q -b*m*g*l/q 0%计算并显示传递函数的极点pr,p,k=residue(num,den);s=p%-求

11、传递函数的脉冲响应并显示-t=0:0.005:5;impulse(num,den,t)%-显示范围：横坐标0-1，纵坐标0-60，此条语句参数可以根据仿真输出曲线调整-axis(0 1 0 60)%-end- 执行上面文件，得到系统传递函数的分子（num）与分母（den）多项式的Matlab表示及系统的开环极点（s）,显示结果如下所示：num = 2.351 0 0den = 1.0000 0.0883 -27.8285 -2.3094 0s = -5.2780 5.2727 -0.08300 由此可知，系统传递函数的多项式表达式为： 系统的开环极点为S1=-5.2780 ,S2=5.2727

12、, S3=-0.0830, S4=0,由于有一个开环极点位于S平面的右半部，开环系统不稳定。Matlab仿真的开环脉冲响应曲线如下所示，系统不稳定。状态空间法状态空间法可以进行单输入多输出系统设计，因此在这个实验中，我们将尝试同时对摆杆角度和小车位置进行控制。为了更具挑战性，给小车加一个阶跃输入信号。我们用 Matlab 求出系统的状态空间方程各矩阵，并仿真系统的开环阶跃响应。在这里给出一个state.m文件，执行这个文件，Matlab将会给出系统状态空间方程的A，B，C和D矩阵，并绘出在给定输入为一个0.2 m的阶跃信号时系统的响应曲线。state.m：Matlab仿真程序如下:%-stat

13、e.m-%-倒立摆状态方程及开环阶跃响应-%-输入倒立摆相关参数-M=1.096;m=0.109;b=0.1;I=0.0034;g=9.8;l=0.25;%-p用于状态方程计算-p=I*(M+m)+M*m*l2;%-输入倒立摆状态方程并显示-A=0 1 0 0; 0 -(I+m*l2)*b/p (m2*g*l2)/p 0; 0 0 0 1; 0 -(m*l*b)/p m*g*l*(M+m)/p 0 B=0; (I+m*l2)/p; 0; m*l/pC=1 0 0 0; 0 0 1 0D=0; 0%-求开环系统的阶跃响应并显示-T=0:0.05:5;U=0.2*ones(size(T);Y,X=

14、lsim(A,B,C,D,U,T);plot(T,Y)%-显示范围，横坐标0-2，纵坐标0-100，此条语句参数可根据仿真输出曲线调整-axis(0 2 0 100)%-end-A = 0 1.0000 0 0 0 -0.0883 0.6293 0 0 0 0 1.0000 0 -0.2357 27.8285 0B = 0 0.8832 0 2.3566C = 1 0 0 0 0 0 1 0D = 0 0仿真曲线如下图：2. 倒立摆系统的PID控制算法设计摆杆角度控制程序如下:M=1.096;m=0.109;b=0.1;I=0.0034;g=9.8;l=0.25;q=(M+m)*(I+m*l2

15、)-(m*l)2;num1=m*l/q 0 0;den1=1 b*(I+m*l2)/q -(M+m)*m*g*l/q -b*m*g*l/q 0;G1=tf(num1,den1)Kp=100;Ki=1;Kd=35;numPID=Kd Kp Ki;denPID=1 0;num=conv(num1,denPID);den=polyadd(conv(denPID,den1),conv(numPID,num1);r,p,k=residue(num,den);s=pt=0:0.005:5;impulse(num,den,t)axis(0 5 -0.005 0.025)MATLAB仿真分析：本实验是对摆杆角

16、度进行控制，首先将控制器各参数取值为Kp=1、Ki=1、kd=1，可得到图1所示的响应曲线，此时闭环系统不稳定。要想得到稳定系统，可以适当加大比例反馈系数Kp.，取Kp=100,Ki=1,Kd=1仿真结果见图2，由响应曲线看出此时闭环系统处于稳定状态、但响应速度不够快,且振荡次数太多，超调量也大:取Kp=100,Ki=1,Kd=35，系统响应曲线如图3所示、这时的超调量和响应速度已经比较合适，但是还可以继续改善:再增加Kp，取Kp=1000,Ki=1，Kd=35,便可得到响应曲线图4，由图可以看出此时系统的稳态及暂态性能都已经非常理想,其中稳定时间约为0.13s超调量约为2.0%. 图1：Kp

17、=1、Ki=1、Kd=1图2：Kp=100,Ki=1,Kd=1图3：Kp=100,Ki=1,Kd=37 图4：Kp=1000,Ki=1，Kd=35车的位置控制算法如下:M=1.096;m=0.109;b=0.1;I=0.0034;g=9.8;l=0.25;q=(M+m)*(I+m*l2)-(m*l)2;num1=m*l/q 0 0;den1=1 b*(I+m*l2)/q -(M+m)*m*g*l/q -b*m*g*l/q 0;num2=-(I+m*l2)/q 0 m*g*l/q;den2=den1;Kp=95;Ki=1;Kd=30;numPID=Kd Kp Ki;denPID=1 0;num=

18、conv(num2,denPID);den=polyadd(conv(denPID,den2),conv(numPID,num1);r,p,k=residue(num,den);s=pt=0:0.005:5;impulse(num,den,t)axis(0 5 -0.1 0.58)仿真曲线如下：从以上的仿真结果可以看出，当摆角处于很好的闭环控制下时，小车位置处于失控状态，会沿着某一方向运动下去。5.参考资料1 曲延滨. 倒立摆课程设计指导书， 20092胡寿松.自动控制原理（第五版）.科学出版社，2007 3 刘卫国MATLAB程序设计教程北京：中国水利水电出版社，20054 固高科技有限公司.直线倒立摆安装与使用手册R1.0,20055 固高科技有限公司. 倒立摆与自动控制原理实验,2005专心-专注-专业