控制系统设计9195914.docx

1、控制系统设计9195914(此文档为word格式，下载后您可任意编辑修改！)第1章 概述控制系统设计是控制科学与工程系学生的一门必修专业基础课，课程中的一些概念相对比较抽象，如系统的稳定性、可控性、收敛速度和抗干扰能力等。倒立摆系统是一个典型的非线性、强耦合、多变量和不稳定系统，作为控制系统的被控对象，它是一个理想的教学实验设备，许多抽象的控制概念都可以通过倒立摆直观地表现出来。本课程设计的目的是让学生以一阶倒立摆为被控对象，了解用古典控制理论设计控制器（如PID控制器）的设计方法和用现代控制理论设计控制器（最优控制）的设计方法，掌握MATLAB仿真软件的使用方法及控制系统的调试方法，加深学生

2、对所学课程的理解，培养学生理论联系实际的能力。本课程设计的被控对象采用固高公司生产的GIP-100-L型一阶倒立摆系统，课程设计包括三方面的内容：（1）建立一阶倒立摆的线性化数学模型；（2）倒立摆系统的PID控制器设计、MATLAB仿真及实物调试；（3）倒立摆系统的最优控制器设计、MATLAB仿真及实物调试。1.1 实验设备简介一阶倒立摆系统的结构示意图如图1所示。图1 一阶倒立摆结构示意图系统组成框图如图2所示。图2 一阶倒立摆系统组成框图系统是由计算机、运动控制卡、伺服机构、倒立摆本体和光电码盘几大部分组成的闭环系统。光电码盘1将小车的位移、速度信号反馈给伺服驱动器和运动控制卡，摆杆的角度

3、、角速度信号由光电码盘2反馈给运动控制卡。计算机从运动控制卡中读取实时数据，确定控制决策（小车运动方向、移动速度、加速度等），并由运动控制卡来实现该控制决策，产生相应的控制量，使电机转动，通过皮带，带动小车运动，保持摆杆平衡。1.2 设计内容1建立一阶倒立摆数学模型在自动控制理论课程中，有一章专门讲述控制系统的数学模型的建立方法，并将非线性数学模型在一定条件下化简成线性数学模型，在此以一阶倒立摆为例，建立其数学模型，并在摆角附近将其非线性数学模型线性化，学生通过实际数学模型的推导，加深对所学内容的理解。2. 控制系统的MATLAB仿真自动控制理论（古典部分）中所讲的控制器的设计方法很多，如根轨

4、迹设计法、频率特性设计法和PID设计法，在实际系统中PID控制器应用最多，在本课程设计中选择PID控制器，PID控制器的特点是只能对单变量（此处为摆杆角度）进行控制。在现代控制理论中，普遍采用的控制方法是最优控制（LQR控制），最优控制可对多变量进行控制（如同时控制摆杆角度和小车位置），由于现代控制理论课时较少，对最优控制的讲解也很少，在这里通过对倒立摆的控制，让学生理解单变量控制和多变量控制的差别。本部分课程设计的目的是学习PID控制器和LQR控制器的设计方法，了解控制器参数对系统性能指标的影响，学会根据控制指标要求和实际响应调整控制器的参数，加深学生对所学内容的理解。系统的MATLAB控制

5、软件已编好，学生在进行控制系统设计时，只需将控制器参数输入，就可观察到控制系统的输出，仿真的目的一方面是让学生得到满足系统性能指标的控制器参数，另一方面是让学生将理论分析与仿真结果进行对比，更直观地理解各参数对控制性能的影响。3. 倒立摆控制系统实物调试倒立摆系统实验装置如图3所示。具体实验步骤如下：（1）将小车推到导轨正中间位置，并且使摆杆处于自由下垂的静止状态。（2）给计算机和电控箱通电。（3）打开计算机，在DOS操作系统下启动程序，并使系统处于准备状态。（4）起摆：由于PID控制只能控制摆杆的摆角，不能控制小车的位置，若用PID控制方式起摆，可能在摆杆竖起前就已出现“撞墙”现象，因此先采

6、用最优控制方式起摆，等摆杆立起来并稳定下来之后，选择控制器菜单中的控制方式，输入参数，观察小车和摆杆的运动。由于PID控制器只对摆杆角度进行控制，所以在PID控制中小车可能向一个方向运动，此时需用手轻轻扶一下摆杆，以避免小车“撞墙”。最优控制方式可同时对摆杆角度和小车位置进行控制，因此不会出现“撞墙”现象。（）观察控制效果：用金属棒碰一下摆杆，观察倒立摆在干扰信号作用下的输出响应。若不能达到指标要求，分析原因，重新设计，直到对实际系统的控制达到满意的结果。具体实验操作见附录1。图3 倒立摆实验装置1.3 实验准备及注意事项1为了安全起见，在进行系统连线、拆卸与安装前，必须关闭系统所有电源。2使

7、用前请仔细检查连线，确保接线正确无误。如果码盘连线接反或断线，将会发生冲撞。3为了避免设备失控时造成人身伤害，操作时人员应该与设备保持安全距离，一定不要站在摆的两端。4如果发生异常，按下空格键，系统会提示“急停按钮被触发，无法继续控制系统，按任意键退出程序。”5超速时，系统会自动关闭伺服电机，并出现“小车失速，系统被终止，按任意键退出程序。”1.4 课程设计任务书哈尔滨工业大学课程设计任务书 姓 名： 院 （系）： 专 业： 班 号： 任务起至日期： 年 月 日至 年 月 日 课程设计题目： 一阶倒立摆控制器设计 已知技术参数和设计要求：本课程设计的被控对象采用固高公司的一阶倒立摆系统GIP-

8、100-L。系统内部各相关参数为：小车质量 0.5 Kg ；摆杆质量0.2 Kg ；小车摩擦系数0.1 N/m/sec ；摆杆转动轴心到杆质心的长度0.3 m ；摆杆惯量0.006 kg*m*m ；采样时间0.005秒。设计要求：1推导出系统的传递函数和状态空间方程。用Matlab进行脉冲输入仿真，验证系统的稳定性。2设计PID控制器，使得当在小车上施加1N的脉冲信号时，闭环系统的响应指标为：（1）稳定时间小于5秒（2）稳态时摆杆与垂直方向的夹角变化小于0.1 弧度3设计最优控制器，使得当在小车上施加0.2m的阶跃信号时，闭环系统的响应指标为：（1）摆杆角度和小车位移的稳定时间小于5秒（2）的

9、上升时间小于1秒（3）的超调量小于20度（0.35弧度）（4）稳态误差小于2%。 工作量：1. 建立一阶倒立摆的线性化数学模型；2. 倒立摆系统的PID控制器设计、MATLAB仿真及实物调试；3. 倒立摆系统的最优控制器设计、MATLAB仿真及实物调试。 工作计划安排： 同组设计者及分工： 各项工作独立完成 指导教师签字_ 年 月 日 教研室主任意见： 教研室主任签字_ 年 月 日*注：此任务书由课程设计指导教师填写。第2章 一阶倒立摆的数学模型设计目的：建立一阶倒立摆系统的数学模型，并在摆角附近将其非线性数学模型线性化，学生通过实际数学模型的推导，加深对系统建模和模型线性化问题的理解。对系统

10、加入输入信号，进行Matlab 仿真，理解不稳定极点对系统稳定性的影响。设计要求：写出系统的动态方程，得出传递函数和状态空间方程。用Matlab进行脉冲输入仿真，验证系统的稳定性。2.1 一阶倒立摆数学模型的推导对系统建立数学模型是系统分析、设计的前提，而一个准确又简练的数学模型将大大简化后期的工作。为了简化系统分析，在实际的模型建立过程中，要忽略空气流动阻力，以及各种次要的摩擦阻力。这样，可将倒立摆系统抽象成小车和匀质刚性杆组成的系统，如下图所示。本系统内部各相关参数定义如下： 小车质量 摆杆质量 小车摩擦系数 摆杆转动轴心到杆质心的长度 摆杆惯量 加在小车上的力 小车位置 摆杆与垂直向上方

11、向的夹角 摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）下图是系统中小车和摆杆的受力分析图。其中，和为小车与摆杆相互作用力的水平和垂直方向的分量。注意：在实际倒立摆系统中检测和执行装置的正负方向已经完全确定，因而矢量方向定义如图，图示方向为矢量正方向。应用Newton方法来建立系统的动力学方程过程如下：分析小车水平方向所受的合力，可以得到以下方程：由摆杆水平方向的受力进行分析可以得到下面等式：即： 把这个等式代入上式中，就得到系统的第一个运动方程： （2-1）为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程：即： 力矩平衡方程如下：注意：此方程中力矩

12、的方向，由于，故等式前面有负号。合并这两个方程，约去和，得到第二个运动方程： （2-2）1微分方程模型设，当摆杆与垂直向上方向之间的夹角与1（单位是弧度）相比很小，即 时，则可以进行近似处理：，。为了与控制理论的表达习惯相统一，即一般表示控制量，用来代表被控对象的输入力，线性化后得到该系统数学模型的微分方程表达式： (2-3) 2传递函数模型对方程组（2-3）进行拉普拉斯变换，得到 （2-4）注意：推导传递函数时假设初始条件为0。由于输出为角度，求解方程组（2-4）的第一个方程，可以得到把上式代入方程组（2-4）的第二个方程，得到整理后得到以输入力为输入量，以摆杆摆角为输出量的传递函数：其中

13、若取小车位移为输出量，可得传递函数：3状态空间数学模型由现代控制理论原理可知，控制系统的状态空间方程可写成如下形式：方程组（2-3）对解代数方程，得到如下解：整理后得到系统状态空间方程：以上就是一阶倒立摆小车系统的状态空间表达式。2.2 系统MATLAB仿真和开环响应1传递函数在Matlab中，拉普拉斯变换后得到的传递函数可以通过计算并输入分子和分母矩阵来实现。假设系统内部各相关参数为： 小车质量 0.5 Kg 摆杆质量 0.2 Kg 小车摩擦系数 0.1 N/m/sec 摆杆转动轴心到杆质心的长度 0.3 m 摆杆惯量 0.006 kg*m*m 采样时间 0.005秒文件trans.m用于求

14、系统传递函数、传递函数的极点以及开环脉冲响应。% trans.m % 倒立摆传递函数、开环极点及开环脉冲响应% 输入倒立摆传递函数 G(S)=num/denM = 0.5;m = 0.2;b = 0.1;I = 0.006;g = 9.8;l = 0.3;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 % 求传递函数的脉冲响应并显示t=0:0.005:5;impulse(num,den,t) % 显示范围：横坐标0-1，纵坐标0-60，此条语句参数可根据仿真输出曲线调整axis(0 1 0 60)grid% end 执行上面的文件，得到系统传递函数的分子（num）与分母（den）多项式的Matlab 表示及系统的开环极点（s），显示结果如下所示： transnum = 4.5455 0 0den = 1.0000 0.1818 -31.1818 -4.4545 0s = -5.6041 5.5651 -0.1428 0 由此可知，系统传递函数的多项式表达式为：系统的开环极点为、，由于