课程设计报告——基于双闭环PID控制的一阶倒立摆控制系统设计.doc

基于双闭环pid控制的一阶倒立摆控制系统设计 目 录摘要i abstract第一章 绪论1.1引言1.2基于双闭环pid控制的一阶倒立摆控制系统的简介第二章 系统模型建立与验证2.1对象模型22 模型验证第三章 基于双闭环pid控制的一阶倒立摆控制系统设计3.1 双闭环pid控制器设计3.2 .控制器的设计3.3 仿真实验 摘要: 对单级倒立摆系统的平衡控制问题进行了研究.。首先建立了系统的数学模型 ,然后采用 pid 方法设计控制器 ,最后对控制系统进行了仿真实验研究. 仿真实验结果不仅证明了pid方案对系统平衡控制的有效性，同时也展示了它们的控制品质和特性.。关键词: 倒立摆;线性控制;仿真abstract: the balance control of a single inverted pendulum system was focused. first , the mathematical model of the system was built , then the pid method was adopted respectively to design the controllers, at last , the simulation experiments were finished, by the simulation experiment results, not only the validity of the three methods were confirmed , but their qualities and characteristics were shown fully also.key words: inverted pendulum , linear control , simulation 第一章 绪论1.1引言倒立摆是典型的快速、多变量、非线性、绝对不稳定系统。早在20世纪50年代,麻省理工学院 (mit)的控制论专家就根据火箭发射助推器原理设计出一阶倒立摆实验设备,此后其控制方法和思路在军工、航天、机器人领域和一般工业过程中都有着广泛的用途,如机器人行走过程中的平衡控制、火箭发射中的垂直度控制、卫星发射架的稳定控制、飞机安全着陆、化工过程控制以及日常生活中所见的任何重心在上、支点在下的控制问题等,均涉及到“立摆问题”。事实上,人们一直在试图寻找不同的控制方法来实现对倒立摆的控制,以便检查或说明该方法对严重非线性和绝对不稳定系统的控制能力。matlab是美国mathworks软件公司于1984年推出的一种用于科学计算的高性能语言。它集数值计算、图形图像显示以及编程于一体,是常用的控制系统分析与设计工具。1990年,math works软件公司为matlab提供了新的控制系统图形化模型输入与仿真工具simulink。这是mat-lab的一个扩展软件模块。该模块提供了一个建模、分析与仿真等多种物理与数学问题软件环境,并为图形用户界面提供了动态系统的结构方块图模型,从而使用户可以既快又方便地对系统进行建模、仿真 而不必写任何代码程序。因此,该工具很快就在控制工程界获得了广泛的认可,并使仿真软件进入了系统模型的图形组态阶段。本文讨论了基于matlab的一阶倒立摆双闭环pid控制系统的建模与仿真,并整定出最优控制参数,同时分析了其鲁棒性。1.2基于双闭环pid控制的一阶倒立摆控制系统的简介1 一阶倒立摆简介本实验的被控对象为固高公司的 gip21002l 型一阶倒立摆系统 ,一阶倒立摆的结构原理图如图所示 ,一阶倒立摆系统的组成框图如图所示。系统包括计算机、运动控制卡、伺服机构、倒立摆本体和光电码盘几大部分 ,组成了一个闭环系统 。光电码盘1 将小车的位移、速度信号反馈给伺服驱动器和运动控制卡 ,摆杆的位置、速度信号由光电码盘 反馈给控制卡。计算机从运动控制卡中读取实时数据 ,确定控制决策(小车向哪个方向移动、移动速度、加速度等) ,并由运动控制卡来实现该控制决策 ,产生相应的控制量 ,使电机转动 ,通过皮带带动小车运动 ,保持摆杆平衡。第二章 系统模型建立2.1对象模型 2.1.1 一阶倒立摆精确模型 当小车的质量m=1kg；倒摆振子的质量m=1kg；倒摆长度2l=0.6m；g=10m/时得2.1.2 一阶倒立摆简化模型 若只考虑在其工作点=0附近（）的细微变化，则可近似为 得到简化的近似模型为 2.1.3 电动机、驱动器及机械传动装置的模型假设：选用日本松下电工msma021型小惯量交流伺服电动机，其有关参数如下：驱动电压：u=0100v 额定功率：pn=200w额定转速：n=3000r/min 转动惯量：j=310-6kg.m2额定转矩：tn=0.64nm 最大转矩：tm=1.91nm电磁时间常数：tl=0.001s 电机时间常数：tm=0.003s 经传动机构变速后输出的拖动力为：f=016n；与其配套的驱动器为：msda021a1a，控制电压：uda=010v。 若忽略电动机的空载转矩和系统摩擦，就可以认为驱动器和机械传动装置均为纯比例环节，并假设这两个环节的增益分别为kd和km。 对于交流电机，传函近似为：由于是小惯性电机，时间常数t1、t2都很小，可以近似为一比例环节kv。电动机、驱动器、机械传动装置三个环节可以合成一个比例环节。 22 模型验证2.2.1模型封装采用仿真实验的方法在matlab的simulink图形仿真环境下进行模型验证试验，其原理图3，其中，上半部分是精确模型仿真图，下半部分为简化模型仿真图。 图3 模型验证原理图利用simulnk压缩子系统功能验证原理图更加简捷地表示为图4的形式。图4 系统封装后的框图其中：fcn:(0.12*u1+0.036*sin(u3)*power(u2,2)-0.9*sin(u3)*cos(u3)/(0.24-0.09*power(cos(u3),2)fcn1:(0.3*cos(u3)u1+0.09*sin(u3)cos(u3)*power(u2,2)-6*sin(u3)/(0.09*power(u3,2)-0.24)fcn2:0.8*u1-0.6*u3fcn3:40*u3-2.0*u12.2.2 实验设计假定在倒立摆初始状态下（=0，x=0），施加微小冲击力作用，有经验知小车向前移动，摆杆倒下。用仿真实验验证正确的数学模型的必要性。2.2.3 绘图子程序% inverted pendulum% test of the discrete model in open loop% singnals recuperation %将导入到xy.mat中的仿真实验数据读出load pid.matt=signals(1,:);f=signals(2,:);x=signals(3,:);q=signals(4,:);xx=signals(5,:);qq=signals(6,:);figure(1)hf=line(t,f(:);xlabel(time)ylabel(force in n)axis(0 1 0 0.12)axet=axes(position,get(gca,position),. xaxislacation,bottom,. yaxislocation,ringt,color,none,. xcolor,k,ycolor,k); ht=line(t,x,color,r,parent,axet); ht=line(t,xx,color,r,parent,axet); ylabel(evolution of the x position) axis(0 1 0 0.1) title(x and xresponse in meters to a f(t) step of 0.1 n) gtext(leftarrow f(t),gtext(x(t)rightarrow), gtext(leftarrow x(t) figure(2) hf=line(t,f(:); xlabel(time) ylabel(force in n) axis(0 1 0 0.12) axet=axes(position,get(gca,position),. xaxislacation,bottom,. yaxislocation,ringt,color,none,. xcolor,k,ycolor,k); ht=line(t,q,color,r,parent,axet); ht=line(t,qq,color,r,parent,axet); ylabel(angle evolution) axis(0 1 -0.3 0) title(response theta(t)andtheta(t) in rd a f(t) step of 0.1 n) gtext(leftarrow f(t),gtext(theta(t)rightarrow), gtext(leftarrow theta(t)仿真实验执行该程序的结果如图所示。可知，在微小冲击力作用下，摆杆倒下，小车位置增加，符合实验设计。同时，可以看出精确模型与简捷模型非常接近，可以用近似模型代替精确模型。 第二章 基于双闭环pid控制的一阶倒立摆控制系统设计 3.1.1 双闭环pid控制器设计从一阶倒立摆的结构图中不难看出，对象传递函数中含有不稳定零极点，该系统为一个“自不稳定的非最小相位系统”。由于一阶倒立摆系统位置伺服控制的核心是“在保证摆杆不到的条件下，使小车位置可控”，因此，依据负反馈闭环原理，将系统小车位置作为“外环”，而降摆杆摆角作为“内环”，则摆角作为外环的一个扰动，能够得到闭环系统的有效抑制。所以，设计一阶倒立摆位置伺服控制如下图所示，剩下的问题就是去解决控制器（校正装置）的d1(s)/d1（s）、d2(s)/d2(s)确定。3.1.2 内环控制器的设计 一）.控制器结构选择考虑到对象为非线性的自不稳定系统，顾采用反馈校正进行控制。控制框图如图其中，ks=1.6为伺服电动机与减速机构的等效模型 ，d2（s）=k为比例环节。下图为不同控制器内环系统的根轨迹。可以看出，采用pd控制的反馈控制器使系统简单、稳定。所以，采用pd结构控制器。 a)为pid控制 b)为pi控制 c)为pd控制。（一）二）.控制器参数整定首先暂定k=-20。这样可以求出内环的传递函数为： 解得： 系统内环传递函数为：三）.系统内环的动态跟随性能仿真实验 根据的道德内环系统的闭环传函，得到仿真模型。如图所示。编写绘图子程序如下：load simu.matt=signals(1,:);x=signals(2,:);hf=line(t,x(:);figure(1)axis(0 2 0 1.2)grid onxlabel(time(s)ylabel(the response of the step siggals(100%)title(respones)得到仿真图如下：（ 3.1.3 外环控制器的设计 可见，系统开环传递函数可视为一个高阶（4阶）且带有不稳定零点的“非最小相位系统”，为了便于设计，需要首先对系统进行一些简化处理（否则，不便利用经典控制理论与方法对它进行设计）。(一）.系统外环模型的降阶 （1）对内环等效闭环传递函数的近似处理 由（1）得：由（2）得：（2）对象模型g1(s)的近似处理 （4）有（3）得 有（4）得经上述处理后，系统被简化为： 近似条件为你：3.2 .控制器的设计图5给出了系统外环前向通道上传函的等效过程，从最终的简化模型上不难看出，这是一个二阶系统；鉴于一阶倒立摆位置伺服控制系统对抗干扰性能与跟随性能的要求（对摆杆长度、质量的变化应具有一定的抑制能力，同时可使小车有效定位），可以将外环系统设计成典型二型的结构形式。同时，系统还应满足前面各环节的近似条件，即系统外环的截至角频率wc1.58.图5 模型简化图为了满足以上对系统的设计要求，不难发现所需要加入的调节器d1(s)也应为pd的形式。设加入的调节器为d1(s)=kp(ts+1),同时，为了系统有较好的跟随性能，采用单位反馈（d1(s)=k=1）来构成外环反馈通道，如下图所示。此时系统开环传函为： 图6 闭环系统结构图取 。则系统开环传函为：再由“典型型”系统bode图特性（）知：则： 有以上可求出外环调节器的两个参数：，这样得到完整的系统仿真 结构图。如图6。 图6 系统仿真结构图3.3 仿真实验综合上述内容，有图5-35所示的simulink仿真系统结构图。需要强调的是，其中的对象模型为精确模型的封装子系统形式。系统仿真绘制子程序及仿真结果如下：3.3.1 画子程序 inverted pidsingnals recuperation将导入到pid.mat中的仿真实验数据读出load pid.matt=signals(1,:);q=signals(2,:);x= signals(3,:);drawing x(t)and theat(t)response signals画小车位置和摆杆角度的响应曲线figure(1)hf=line(t,q(:);grid onxlabel(time)ylabel(angle evolution)anis(0 10 -0.3 1.2)annet=anes(position,get(gca,position),. xaxislocation,bottom,. yaxislocation,right,color,none, xcolor,k,ycolor,k);ht=line(t,x,color,r,parent,axet);ylabel(evolution of the x position)axis(0 10 -0.3 1.2)title(theta(t)and x(t) response to a step input)gtext(leftarrow x(t),gtext(theta(t)uparrow)其中：fcn:0.12u1+0.036sin(u3)power(u2,2) -0.9sin(u3)cos(u3)(0.24-0.09po