现代控制工程大作业
现代控制大作业
现代控制最后大作业一、控制系统数学模型控制系统的闭环数学模型为150G= -----------------------s^3 + 11s^2 + 30 s+150用matlab进行系统仿真,确定系统时域性能指标;G=tf([150],[1,11,30,150])Transfer function:150---------------------------s^3 + 11 s^2 + 30s+150step(G)由上图可知:超调量Mp=50.1% 上升时间tr=0.334s峰值时间tp=0.911s 调节时间ts=5.1s下面采用两种控制器分别对上述系统进行控制,使得校正后的系统满足M 。
20%P二模型校正1、采用经验试凑法矫正首先,利用Simulink构建原系统,并观察输出情况。
第二,在前向通路加入PID矫正单元,并利用试凑法调节PID参数,直至系统超调量满足要求。
当Kp=1.08,Ki=0.1,Kd=1时,系统超调量Mp=19.8%满足矫正要求2、采用根轨迹串联校正控制器用根据轨迹法进行超前校正的一般步骤为:1)根据对系统静态性能指标和动态性能指标的要求,分析确定希望的开环增益和闭环主导极点的位置。
2)画出校正前系统的根轨迹,判断希望的主导极点位于原系统的根轨迹左侧,以确定是否应加超前校正装置。
3)根据题目要求解出超前校正网络在闭环主导极点处应提供的相位超前角。
4)根据图解法求得Gc(s)的零点和极点,进而求出校正装置的参数。
5)画出校正后系统的根轨迹,校核闭环主导极点是否符合设计要求。
M ,校正装置传递函数为Gc(s)=Kc。
采用串联校正,使其满足条件20%P经以上分析可知:在原有的系统中串联一个合适的具有比例性质的控制器,可使得系统的根轨迹向左平移,使得系统的稳定性变得更好,可以实现在原有系统的基础上使得a在更大范围内都满足题意。
绘制根轨迹程序为:150KcG= ------------------------------s^3 + 11 s^2 + 30 s+150G=zpk([],[0 -5 -6],[150]);figure(1);rlocus(G);figure(2);step(feedback(0.4*G,1));原系统传函的根轨迹图矫正后的阶跃响应图通过以上图三可以知道:超调量Mp=14.5%上升时间tr=0.612峰值时间tp=1.5s 调节时间ts=3.04s通过以上的校正后可以看出超调量由原来的50%变为14.5%可以满足条件,实现校正,此时Kc=0.4。
现代控制理论大作业
现代控制理论大作业一、位置控制系统----双电位器位置控制系统由系统分析可知,系统的开环传递函数:2233.3s =s s 2*0.07s*s 205353G()(+1)*(++1)另:该系统改进后的传递函数:223.331s =s s 2*0.07s*s 3455353G ()(+1)*(++1)1、时域数学模型<1>稳定性>> s=tf('s');>> G=33.3/(s*(s/20+1)*(s^2/53^2+2*0.07*s/53+1)); >>sys=feedback(G,1); >> sysTransfer function:9.915e007 -----------------------------------------------------------53 s^4 + 1453 s^3 + 1.567e005 s^2 + 2.978e006 s + 9.915e007>> pzmap(sys)由零极点图可知,该系统有四个极点,没有零点,其中两个在左半s 开平面上,两个在s 平面的虚轴处,则,四个极点的坐标分别是:>> p=pole(sys)p =0.0453 +45.2232i0.0453 -45.2232i-13.7553 +26.9359i-13.7553 -26.9359i系统的特征方程有的根中有两个处于s的右半平面,系统处于不稳定状态<2>稳态误差分析稳态误差分析只对稳定的系统有意义,系统(G)处于不稳定状态,所以不做分析。
改进后系统(G1)如下,求其特征方程的极点:>> s=tf('s');>> G1=3.33/(s*(s/345+1)*(s^2/53^2+2*0.07*s/53+1));>> sys2=feedback(G1,1);>>p=pole(sys2);p =1.0e+002 *-3.4492-0.0206 + 0.5258i-0.0206 - 0.5258i-0.0338可以看出,改进后的传递函数G1的四个极点都在s平面的右半开平面上,则系统G1是稳定的,故对此系统做稳态误差分析:由系统G1的开环传递函数在原点处有一个极点,故属于1型系统。
现代控制理论大作业
分类号:TH89 单位代码:10110学号:中北大学综合调研报告题目: 磁盘驱动器读写磁头的定位控制系别: 计算机科学与控制工程学院专业年级: 电气工程与智能控制2014级姓名: 何雨贾晨凌朱雨薇贾凯张钊中袁航学号: 14070541 39/03/04/16/33/47指导教师: 靳鸿教授崔建峰讲师2017年5月7日摘要硬盘驱动器作为当今信息时代不可缺少的存储设备,在人们日常生活中正扮演着越来越重要的角色,同时它也成为信息时代科学技术飞速发展的助推器。
然而,随着信息量的日益增长,人们对硬盘驱动器存储容量的要求越来越高。
但另一方面由于传统硬盘驱动器的低带宽、低定位精度,导致磁头很难准确地定位在目标磁道中心位置,从而限制了存储容量的持续增加。
自IBM公司于1956年向全球展示第一台磁盘存储系统R.AMAC以来,随着存储介质、磁头、电机及半导体芯片等相关技术的不断发展,硬盘的存储容量成倍增长、读写速度不断提高。
要保证可靠的读写性能,盘片的转速控制和磁头的定位控制问题具有重要意义。
其中磁头的定位控制主要包括寻道控制与定位跟踪控制两个问题,如PID控制、自适应控制、模态切换控制等,这些控制方法大大提高了硬盘磁头伺服系统的性能。
为达到更高的精度,磁头双级驱动模型成近年的研究热点,多种控制策略已有相关报道,但目前仍处于实验水平。
关键词: 磁盘驱动器;磁头;定位;控制AbstractHard disk drive (HDD), acted as requisite storage equipment in current information age,plays a more and more vital role in people’s daily life, and it becomes a roll booster in rapid development of science and technology. However, with the increase of information capacity, we put forward a severe request for HDD data storage capacity. Unfortunately, due to the low bandwidth, low positioning accuracy in conventional HDD, magnetic head is hard to be positioned onto the destination track center, thus it limits the continuing increase in storage capacity.Since IBM brought the first disk-the random access memory accounting machine(RAMAC) to market in 1956, the storage capacity and read/write speed have continuously increased along with the development of the techniques of media,read/write head, actuators and semiconducting chips. The problems of R/W head's settling control is definitely important in order to ensure the reliability of read and write performance. Track seeking and track following are two main stages of the hard disk servo system. Researchers have developed kinds of control strategies to implement the servo control from PID control to advanced control methods.Dual-stage actuator has attracted many researchers and engineers for its broaderbandwidth compared with single-stage actuator.Key Words:Hard Disk Drive;Heads; Location; Control第1章磁盘驱动器的介绍自上世纪50年代计算机发明以来,随着科技的进步,软硬件技术都获得了相当大的发展。
现代控制理论大作业课件
自适应控制设计需要建立合适的自适应算法和控制器结构。
详细描述
自适应控制设计需要建立合适的自适应算法和控制器结构, 以确保控制器能够实时地调整自身参数并保证系统的最优 性能。此外,还需要对算法和控制器结构进行充分的验证 和测试,以确保其在实际应用中的可靠性和有效性。
06
非线性系统理论
非线性系统的基本性质
齐次性
线性系统的输出与输入成正比,即当输入信 号乘以一个常数时,输出信号也乘以相同的 常数。
叠加性
对于线性系统的多个输入信号,其输出信号等于各 个输入信号单独作用于系统所产生的输出信号之和。
时不变性
线性系统的输出信号与输入信号在不同时刻 的比值保持恒定,即系统特性不随时间变化。
线性系统的状态空间表示
平衡点稳定性
分析非线性系统平衡点的稳定性,可以通过求 解系统的线性化方程来获得。
动态稳定性
动态稳定性是指非线性系统在受到扰动后恢复稳定的能力。
非线性系统的控制设计
状态反馈控制
通过测量系统的状态变量, 并使用状态反馈控制器来 控制非线性系统的输出。
滑模控制
滑模控制是一种变结构控 制方法,通过设计滑模面 和滑模控制器来实现对非
最优控制问题的描述
01
最优控制问题是在给定初始和终端状态约束下,寻找一个控制输 入,使得系统状态在满足约束条件下,某个性能指标达到最优。
02
性能指标通常包括系统状态和控制输入的能量、时间和成 本等。
03
约束条件可以是系统状态、控制输入或性能指标的限制。
极小值原理
极小值原理是求解最优控制问题的一种方法,它基于动态规划的思想,通 过求解一系列的子问题来找到最优解。
状态方程
描述系统内部状态变量随时间变化的数学方程,通常 表示为矩阵形式。
现代控制理论大作业7页
现代控制理论大作业老师:周晓敏姓名:李维奇班级:机研141班学号:s2*******2019年1月一.系统的工程背景及物理描述超精密机床是实现超精密加工的关键设备,而环境振动又是影响超精密加工精度的重要因素。
为了充分隔离基础振动对超精密机床的影响,目前国内外均采用空气弹簧作为隔振元件,并取得了一定的效果,但是这属于被动隔振,这类隔振系统的固有频率一般在2Hz 左右。
上图表示了亚微米超精密车床隔振控制系统的结构原理,其中被动隔振元件为空气弹簧,主动隔振元件为采用状态反馈控制策略的电磁作动器。
上图表示一个单自由度振动系统,空气弹簧具有一般弹性支承的低通滤波特性,其主要作用是隔离较高频率的基础振动,并支承机床系统;主动隔振系统具有高通滤波特性,其主要作用是有效地隔离较低频率的基础振动。
主、被动隔振系统相结合可有效地隔离整个频率范围内的振动。
床身质量的运动方程为:p F ——空气弹簧所产生的被动控制力a F ——作动器所产生的主动控制力假设空气弹簧内为绝热过程,则被动控制力可以表示为:电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及永久磁铁和电磁铁之间的间隙面积有关,这一关系具有强非线性。
由于系统工作在微振动状况,且在低于作动器截止频率的低频范围内,因此主动控制力可近似线性化地表示为:其中,电枢电流Ia 满足微分方程: 1.性能指标:闭环系统单位阶跃响应的:超调量不大于5%;过渡过程时间不大于0.5秒(∆=0.02)2.实际给定参数:某一车床的已知参数3.开环系统状态空间数学模型的推导过程:对式0y s s =-两边求二次导,.....011()({1[/()]})n p a r r r e e e a y s F F c y k y p V V A y A k I m m ==-+=-++-++对上式再求一次导,其中1/()r r r e e p V V A y A η⎧⎫''⎡⎤=-+⎨⎬⎣⎦⎩⎭则,又由,代入 00(,)()a e emy cy k y my cy k y L R E I y u t k k ηη++++++--+=,即 令状态变量为 , 得系统开环的状态方程为:1223003123e x x x x Rk Lk Rc k Lc Rm x x x x u Lm Lm Lm Lm ⎧⎪=⎪=⎨⎪++⎪=----⎩于是状态空间表达式为:[]1122003312301000010100e x x x x u Rk Lk Rc x k Lc Rm x Lm Lm Lm Lm x y x x ⎧⎡⎤⎡⎤⎪⎢⎥⎢⎥⎡⎤⎡⎤⎪⎢⎥⎢⎥⎢⎥⎢⎥=+⎪⎢⎥⎢⎥⎢⎥⎢⎥⎪⎢⎥⎢⎥⎢⎥⎢⎥++⎪⎣⎦⎣⎦----⎢⎥⎢⎥⎨⎣⎦⎣⎦⎪⎪⎡⎤⎪⎢⎥=⎪⎢⎥⎪⎢⎥⎣⎦⎩代入系统参数,二、系统的定性分析系统的能控能观性根据其能控性矩阵和能观性矩阵是否满秩来判断。
现代控制理论大作业
现代控制理论直流电动机模型的分析姓名:李志鑫班级:测控1003学号:20100203030921直流电动机的介绍1.1研究的意义直流电机是现今工业上应用最广的电机之一,直流电机具有良好的调速特性、较大的启动转矩、功率大及响应快等优点。
在伺服系统中应用的直流电机称为直流伺服电机,小功率的直流伺服电机往往应用在磁盘驱动器的驱动及打印机等计算机相关的设备中,大功率的伺服电机则往往应用在工业机器人系统和CNC铣床等大型工具上。
[1]1.2直流电动机的基本结构直流电动机具有良好的启动、制动和调速特性,可以方便地在宽范围内实现无级调速,故多采用在对电动机的调速性能要求较高的生产设备中。
直流伺服电机的电枢控制:直流伺服电机一般包含3个组成部分:-图1.1①磁极:电机的定子部分,由磁极N—S级组成,可以是永久磁铁(此类称为永磁式直流伺服电机),也可以是绕在磁极上的激励线圈构成。
②电枢:电机的转子部分,为表面上绕有线圈的圆形铁芯,线圈与换向片焊接在一起。
③电刷:电机定子的一部分,当电枢转动时,电刷交替地与换向片接触在一起。
直流电动机的启动电动机从静止状态过渡到稳速的过程叫启动过程。
电机的启动性能有以下几点要求:1)启动时电磁转矩要大,以利于克服启动时的阻转矩。
2)启动时电枢电流要尽可能的小。
3)电动机有较小的转动惯量和在加速过程中保持足够大的电磁转矩,以利于缩短启动时间。
直流电动机调速可以有:(1)改变电枢电源电压;(2)在电枢回路中串调节电阻;(3)改变磁通,即改变励磁回路的调节电阻Rf以改变励磁电流。
本文章所介绍的直流伺服电机,其中励磁电流保持常数,而有电枢电流进行控制。
这种利用电枢电流对直流伺服电机的输出速度的控制称为直流伺服电机的电枢控制。
如图1.2Bm电枢线路图1.2——定义为电枢电压(伏特)。
——定义为电枢电流(安培)。
——定义为电枢电阻(欧姆)。
——定义为电枢电感(亨利)。
——定义为反电动势(伏特)。
现代控制理论大作业
2019/11/21
2
课题背景
系统工作原理
Accelerometer gyroscope 遥控器
电源
左轮 软件编码
AD 采样
滤波
MCU
驱动器1 驱动器2
无线模块
软件编码 右轮
减速机构 左电机
右电机 减速机构
2019/11/21
3
系统工作原理 前进(后仰)
后退(前倾)
2019/11/21
课题背景
后退(纠正后仰) 前进(纠正前倾)
Matlab计算程序:
pole=[-1,-2,-3,-4]; K=place(A,B,pole) A1=A-B*K; pole=[-3,-4,-5,-6]; K=place(A,B,pole) A2=A-B*K; pole=[-8,-9,-10,-11]; K=place(A,B,pole) A3=A-B*K;
figure; hold on; plot(curve1(:,2),'color','blue'); plot(curve2(:,2),'color','red'); plot(curve3(:,2),'color','green'); xlabel(‘时间(s)'); ylabel(‘摆动角度(rad)'); hold off
2019/11/21
N C,CA,CA2 ,CA3
计算结果:
rankc = 4 ranko = 4
结论: k(M)=4 系统完全能控 Rank(N)=4 系统完全能观
稳定性分析
运用Matlab解出矩阵A的特征值如下:
(完整)现代控制理论-大作业-倒立摆
摘要倒立摆系统是一个复杂的、高度非线性的、不稳定的高阶系统,是学习和研究现代控制理论最合适的实验装置。
倒立摆的控制是控制理论应用的一个典型范例,一个稳定的倒立摆系统对于证实状态空间理论的实用性是非常有用的.本文主要研究的是二级倒立摆的极点配置方法,首先用Lagrange方程建立了二级倒立摆的数学模型,然后对二级倒立摆系统的稳定性进行了分析和研究,并给出了系统能控能观性的判别。
基于现代控制理论中的极点配置理论,根据超调量和调整时间来配置极点,求出反馈矩阵并利用Simulink对其进行仿真,得到二级倒立摆的变化曲线,实现了对闭环系统的稳定控制。
关键词:二级倒立摆;极点配置;Simulink目录1.绪论 (1)2 数学模型的建立和分析 (1)2。
1 数学建模的方法 (1)2。
2 二级倒立摆的结构和工作原理 (2)2。
3 拉格朗日运动方程 (3)2。
4推导建立数学模型 (3)3 二级倒立摆系统性能分析 (9)3.1 稳定性分析 (9)3。
2 能控性能观性分析 (10)4 状态反馈极点配置 (11)4。
1 二级倒立摆的最优极点配置1 (11)4.2 二级倒立摆最优极点配置2 (12)5。
二级倒立摆matlab仿真 (14)5。
1 Simulink搭建开环系统 (14)5.2 开环系统Simulink仿真结果 (14)5.3 Simulink搭建极点配置后的闭环系统 (15)5.4极点配置Simulink仿真结果 (16)5.4。
2 第二组极点配置仿真结果 (18)6。
结论 (19)7.参考文献 (20)附录一 (21)1.绪论倒立摆最初诞生于麻省理工学院,仅有一级摆杆,另一端铰接于可以在直线导轨上自由滑动的小车上.后来在此基础上,人们又进行拓展,设计出了直线二级倒立摆、环型倒立摆、平面倒立摆、柔性连接倒立摆、多级倒立摆等实验设备。
在控制理论的发展过程中,为验证某一理论在实际应用中的可行性需要按其理论设计的控制器去控制一个典型对象来验证。
控制工程基础大作业
控制工程基础大作业第一篇:控制工程基础大作业控制工程基础(B)大作业一、查阅相关资料,撰写读书报告通过广泛的资料查阅,以论文形式就控制领域发展状况提交相关读书报告,以培养学生文献检索和阅读能力。
要求:至少阅读5篇2008年以后的文献资料,并在上交时附上相应的参考文献或网址,阅读报告的字数不得少于3000汉字。
二、应用Matlab软件绘制Nyquist图及Bode图1、自己从教材上的例题或者课后作业中选择一个系统,绘制其Bode图。
2、已知系统的开环传递函数为G(s)=100k,用Matlab分s(s+5)(s+10)别绘制k=1,8,20时系统Nyquist 图,并判断系统的稳定性。
试分析k对系统稳定性的影响。
三、控制系统典型环节性能分析熟悉Matlab软件Simulink的基本使用方法,利用Simulink建立各典型环节的仿真模型,并通过仿真得到各典型环节的单位阶跃响应曲线,给出各典型环节相关参数变化对典型环节动态性能的影响。
要求:每个学生的文件名包括姓名、班级和学号(1~2位数字)。
该部分主要是考察大家的自学能力,网上有很多资料可下载、查阅。
大作业以报告形式撰写,由各班长(或者学习委员)负责收齐电子文档,并于6月12日上交,其中电子文档以班为单位打包后发送至邮箱:*****************附:本课程成绩评定方法考勤、课堂表现、平时作业、实验报告、大作业及期末考试相结合,其中:考勤、课堂表现、平时作业与实验报告占20%,大作业占20%,期末考试成绩占60%。
任课教师:江老师2013年5月第二篇:控制工程基础试卷序号一、试写出图示机械装置(物体与地面无摩擦)在外力F作用下的运动(微分)方程,并求该系统的传递函数Y0(s)/F(s)。
(10分)二、试求图示框图的传递函数C(s)/R(s)(方法不限,10分)。
序号三、某单位负反馈控制系统如图。
已知K=4,设输入信号为单位阶跃函数,求:(15分)(1)闭环传递函数C(s)/R(s);(2)单位阶跃响应c(t)(3)调整时间ts;(4)最大超调量。
控制工程大作业
2014年控制工程大作业要求:a)每位同学对每个题目按照要求设计控制方案,画出带控制点的工艺流程图。
b)详细说明设计的控制方案是如何运作的,如何能满足工艺要求。
c)报告采用PDF文档,正文为五号宋体,标题可适当加大字号或用粗体。
纸张采用A4。
d)所有报告均在首页标明姓名、学号、班级,文件以“学号.pdf”的格式保存。
1.考虑如题图1所示的系统,希望将质量分数为50%的NaOH溶液稀释为30%。
NaOH控制阀受其它控制器(未在图中展示)控制。
由于50% NaOH溶液的流量变化频繁,希望设计一个比值控制方案,通过控制稀释水的流量以达到稀释要求。
50%NaOH溶液的标称流量(即设计值)为100 kg/hr。
假设选用的两个流量变送器为质量流量计,其输出信号与质量流量成线性比例关系。
50%NaOH溶液的变送器量程为0到200 kg/hr,水的变送器量程为0到160 kg/hr。
请详细描述其比值控制方案,并计算控制系统所涉及的比值系数。
2题图1 混合过程2. 考虑如题图2所示的制氢反应器。
同锅炉一样,该反应器需要通过燃料和空气进行燃烧才能提供反应所需要的能量(为简单起见,温度控制系统并未全部在图中画出)。
甲烷和水蒸气反应产生氢气,反应方程为4222CH +2H O CO +4H该反应发生在反应器的内部管道中,管道内填充着反应所需的催化剂。
为避免催化剂结焦,保证反应混合物中水蒸气富余十分重要。
如果碳沉积在催化剂表面,将会使催化剂中毒。
通过保证输入混合物中水蒸气富余,就可以避免此种现象。
然而,过多的水蒸气又会引起能量消耗过大。
工程部门估算,应当保证的最佳甲烷-水比为R 1(质量比)。
请设计一套控制方案,可使实际比例维持在最佳比例;并且在生产过程中,无论产量增加或减少,都要保证反应混合物中水蒸气富余(即高于最佳比例)。
注意:甲烷需求量由另外专门的信号控制。
题图2 制氢反应器4.考虑如题图4所示的混合池。
三种原料以一个确定的比例注入水箱之中形成一种混合液,供后续过程使用。
现代控制理论作业
现代控制理论大作业要求:(1)自选一实际物理对象进行研究,建立实际物理系统的状态空间模型;(2)进行原系统的定性分析,包括稳定性、能控性、能观性分析;(3)根据系统提出的性能指标要求(如超调量、超调时间、调节时间等动态 性能指标以及稳态误差等稳态性能指标),进行原系统的仿真分析,和要求的性能指标做对比;(4)对不稳定系统且能镇定的系统,进行镇定控制;(5)对未达到性能指标要求的系统进行状态反馈控制设计,满足系统性能指 标要求;(6)设计状态观测器观测所有状态;(7)设计降阶状态观测器;(可选)(8)最优控制;(9)体会及对课程建议。
1实际物理模型:如图1所示,为一交接车前后连接振动简化模型。
设计一个调节器系统使得在无扰动的情况下,系统保持在零位置上(y1=0)。
其中m1=1,m2=2,k=36,b=0.62系统的描述方程:)()(m )()(m 212122121211y y b y y k yu y y b y y k y-+-=+-+-= 其空间状态模型为:设:。
,,,24132211y x yx y x y x ====[]⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡432121432143210001u 01003.03.018186.06.0-3636-10000100x x x x y y x x x x x x x x 3分析与求解过程:由根轨迹和特征根(a = -0.4500 + 7.3347i -0.4500 - 7.3347i -0.000 0 )知虽实根都为负数但都靠近零轴,是李雅普诺夫定义下的稳定,但存在震荡,所以把希望闭环极点配置在10-s ,10-s ,32-2-s ,322-s ===+=和把最小阶观测器希望极点配置在16-s ,15-s ==来改善系统的性能。
北航-现代控制理论结课大作业
1. 控制系统任务的物理描述为了满足飞机品质的要求,飞机的纵向运动和横侧向运动都需要有能够连续工作的阻尼器,以用来调整飞机的飞行姿态,避免其出现不必要的俯仰和倾斜。
维持飞机纵向运动的阻尼器称为俯仰阻尼器,维持飞机横侧向运动的阻尼器称为偏航阻尼器。
本次课程大作业旨在通过运用Matlab 的经典控制系统设计工具对某型飞机偏航阻尼器进行控制系统的设计。
2. 控制系统对象的数学模型巡航状态下,某型飞机侧向运动的状态空间模型为:111121314111222122232421221313233343132234142434441424()1()()()()2()()()3()()4t x t a a a a b b t x t a a a a b b u t a a a a b b u t x t t a a a a b b x t t x x x x ••••⎡⎤⎢⎥⎢⎥⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦ 1111213141221222324234()()()()()()x t c c c c y t x t c c c c y t x t x t ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦式中: 1()x t :侧滑角(单位为rad )2()x t :偏航角速度(单位为/rad s )3()x t :滚转角速度(单位为/rad s )4()x t :倾斜角(单位为rad )输入向量及输出向量分别为:1()u t :方向舵偏角(单位为rad )2()u t :副翼偏角(单位为rad )1()y t :偏航角速度(单位为/rad s )2()y t :倾斜角(单位为rad )设飞机巡航飞行时的速度为0.8马赫,高度为40000英尺,此时模型的参数为:0.055800.99680.08020.04150.598000.11500.031803.050000.38800.4650000.0805 1.00000A --⎡⎤⎢⎥--⎢⎥=⎢⎥--⎢⎥⎣⎦ 0.007290.00000.475000.007750.15300.1430000B ⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦0 1.000000 1.0C ⎡⎤=⎢⎥⎣⎦0000D ⎡⎤=⎢⎥⎣⎦首先输入飞机状态空间模型参数。
现代控制理论课程设计(大作业)
现代控制理论课程设计报告题目打印机皮带驱动系统能控能观和稳定性分析项目成员史旭东童振梁沈晓楠专业班级自动化112指导教师何小其分院信息分院完成日期 2014-5-28目录1. 课程设计目的 (4)2.课程设计题目描述和要求 (4)3.课程设计报告内容 (4)3.1 原理图 (4)3.2 系统参数取值情况 (5)3.3 打印机皮带驱动系统的状态空间方程 (5)4. 系统分析 (8)4.1 能控性分析 (8)4.2 能观性分析 (8)4.3 稳定性分析 (9)5. 总结 (11)项目组成员具体分工打印机皮带驱动系统能控能观和稳定性分析课程设计的内容如下:1.课程设计目的综合运用自控现代理论分析皮带驱动系统的能控性、能观性以及稳定性,融会贯通并扩展有关方面的知识。
加强大家对专业理论知识的理解和实际运用。
培养学生熟练运用有关的仿真软件及分析,解决实际问题的能力,学会应用标准、手册、查阅有关技术资料。
加强了大家的自学能力,为大家以后做毕业设计做很好的铺垫。
2.课程设计题目描述和要求(1)环节项目名称:能控能观判据及稳定性判据(2)环节目的:①利用MATLAB分析线性定常系统的可控性和客观性。
②利用MATLAB进行线性定常系统的李雅普诺夫稳定性判据。
(3)环节形式:课后上机仿真(4)环节考核方式:根据提交的仿真结果及分析报告确定成绩。
(5)环节内容、方法:①给定系统状态空间方程,对系统进行可控性、可观性分析。
②已知系统状态空间方程,判断其稳定性,并绘制出时间响应曲线验证上述判断。
3.课程设计报告内容3.1 原理图在计算机外围设备中,常用的低价位喷墨式或针式打印机都配有皮带驱动器。
它用于驱动打印头沿打印页面横向移动。
图1给出了一个装有直流电机的皮带驱动式打印机的例子。
其光传感器用来测定打印头的位置,皮带张力的变化用于调节皮带的实际弹性状态。
图1 打印机皮带驱动系统3.2 系统参数取值情况表1打印装置的参数3.3 打印机皮带驱动系统的状态空间方程图2 打印机皮带驱动模型状态空间建模及系统参数选择。
现代控制理论结课大作业
现代控制理论结课大作业一、引言现代控制理论是现代科学技术的重要组成部分,广泛应用于工程控制系统中。
在控制理论课程的学习过程中,结课大作业是一项重要的任务。
本文将介绍现代控制理论结课大作业的相关要求和设计思路。
二、研究背景现代控制理论是控制理论的一个重要分支,它主要研究控制系统的建模、分析和设计方法。
通过运用数学和工程技术知识,利用现代控制理论可以对各种系统进行精确的描述和控制。
因此,现代控制理论在自动控制领域具有广泛的应用。
三、大作业要求现代控制理论结课大作业要求学生能够独立选择一个控制系统并进行详细的研究和设计。
具体要求如下: 1. 选择一个真实的控制系统作为研究对象;2. 系统建模:根据实际情况,选择合适的建模方法,将系统转化为数学模型;3. 系统分析:通过分析系统模型,对系统的稳定性、鲁棒性等进行评估; 4. 系统设计:基于现代控制理论的设计思想,设计适合该系统的控制器; 5. 系统仿真:利用仿真软件对设计的控制系统进行验证和优化; 6. 结果分析和总结:对仿真结果进行分析,总结设计过程和经验教训。
四、设计思路在完成现代控制理论结课大作业时,需要有清晰的设计思路和步骤。
以下是一个可能的设计思路供参考: 1. 选择合适的控制系统:可以选择一个典型的工业控制系统,或者选择一个与个人兴趣相关的系统; 2. 进行系统建模:根据系统的实际情况,选择适合的建模方法,如状态空间法、传递函数法等;3. 系统分析:利用控制理论的知识和工具,分析系统的稳定性、鲁棒性,确定系统的可控性和可观性等性能指标;4. 系统设计:基于现代控制理论,设计一个合适的控制器结构,并选择适当的控制参数;5. 系统仿真:利用仿真软件,对设计的控制系统进行仿真验证,观察系统的响应特性和控制性能; 6.结果分析和总结:根据仿真结果,分析系统的优点和不足之处,并总结设计过程中的经验教训。
五、实例分析下面以一个简单的倒立摆系统为例,介绍如何完成现代控制理论结课大作业。
现代控制理论大作业
现代控制理论大作业现代控制理论大作业1.解:(1).选取状态变量为:x1=y, x2=y’,x3=y’’由题可得:a2=1 , a1=4, a0=5所以x3’=-5x1-4x2-x3+3u系统的状态方程为:x1’=x2x2’=x3x3’=-5x1-4x2-x3+3u输出方程为:y=x1将微分方程表达为矩阵形式即得其状态空间表达式:[x1’; x2’; x3’]=[0,1,0;0,0,1;-5, -4, -1][x1;x2;x3]+[0;0;3]uy=[1, 0, 0][x1;x2;x3](2).选取系统的状态变量为:x1=y-h0ux2=x1’-h1u=y’-h0u’-h1ux3=x2’-h2u=y’’-h0u’’-h1u’-h2u 由题可得:a0=0, a1=3/2, a2=0b0=-1/2, b1=0, b2=1/2, b3=0所以:[h0;h1;h2;h3]=[1 0 0 0;0 1 0 0;3/2 0 1 0;0 3/2 0 1]^-1*[0 ;1/2;0 ;-1/2]=[0;1/2;0;-5/4]取状态变量为:x1=y-h0u=yx2=x1’-h1u=x1’-1/2ux3=x2’-h2u=x2’所以该系统的状态空间表达式为:[x1’;x2’;x3’]=[0 1 0;0 0 1;0 -3/2 0][x1;x2;x3]+[0;1/2;0;-5/4]uy=[1 ,0, 0][x1;x2;x3](3)由题可得:a2=2, a1=3, a0=5;b3=5, b2=0, b1=0, b0=7所以[h0;h1;h2;h3]=[1 0 0 0;2 1 0 0;3 2 1 0;5 3 2 1]^-1*[5;0;0;7] =[5;-10;5;2]取状态变量为:x1=y-h0u=y-5ux2=x1’-h1u=x1’x3=x2’-h2u=x2’所以该系统的状态空间表达式为:[x1’;x2’;x3’;]=[0 1 0;0 0 1;-5 -3 -2][x1;x2;x3]+[5;-10;5;2]u2.经典控制理论是建立在常微分方程稳定性理论和以拉普拉斯变换为基础的根轨迹和奈奎斯特判断理论之上。
现代控制理论大作业 北科
现代控制理论大作业分析对象:汽车悬架系统指导老师:周晓敏专业:机械工程姓名:白国星学号:S2*******1.建模悬架是车轮或车桥与汽车承载部分之间具有弹性的连接装置的总称,具有传递载荷、缓和冲击、衰减振动以及调节汽车行驶中的车身位置等作用。
传统汽车悬驾系统是被动悬驾,其参数不能改变,无法控制其对不同路面激励的响应,因此对不同路面的适应性较差。
为提高汽车的行驶平顺性、操纵稳定性和制动性等性能,人们开始用主动悬架系统来代替传统的被动悬架系统。
主动悬架系统能根据路面的情况通过一个动力装置改变悬挂架的参数,改善汽车的各方面性能。
对悬驾系统进行仿真计算首先要建立悬驾系统动力学模型,随后对所建立的模型进行仿真分析。
为了简化模型,取汽车的一个车轮的悬驾系统进行研究,该模型可简化为一维二自由度的弹簧阻尼质量系统,图1所示为该模型的模拟图。
图1 悬架系统模型的模拟图其中u 为动力装置的作用力,w 为路面位移,x1为车身位移,x2为悬驾位移,用车身位移来度量车身的振动情况,并视为系统的输出。
路面状况以w 为尺度,并视为系统的一个干扰输入。
当汽车从平面落入坑时,w 可用一个阶跃信号来模拟。
u 为主动悬架的作用力,它是系统的控制量。
进行受力分析,由牛顿第二规律可得车身悬架系统的动力学方程为:()()()()()1121212212122s s t m x K x x b x x um x K x x b x x u K w x ⎧=-+-+⎪⎨=-+--+-⎪⎩ 设系统状态变量为:[]1212x x x x x =则上面系统动力学方程可改写为状态空间表达式:x Ax Buy Cx Du=+⎧⎨=+⎩ 其中:()1111222200100001sss t s K K b b A m m m m K K K b b m m m m ⎡⎤⎢⎥⎢⎥⎢⎥--=⎢⎥⎢⎥-+⎢⎥-⎢⎥⎣⎦1220000101t B m K m m ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦[]1000C = []00D = u u w ⎡⎤=⎢⎥⎣⎦Matlab 系统模型程序代码: m1=800;m2=320;ks=10000;b=30000; kt=10*ks;a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1; ks/m2 -(ks+kt)/m2 b/m2 -b/m2]; b=[0 0;0 0;1/m1 0;-1/m2 kt/m2]; c=[1 0 0 0]; d=[0 0]; G=ss(a,b,c,d)[num,den]=ss2tf(a,b,c,d,1);tfsys=tf(num,den)eig(a)figure(1)impulse(G)figure(2)step(G)系统状态空间模型::系统传递函数:系统矩阵特征值:系统矩阵特征值都在左半平面,故系统稳定。
现代控制理论大作业资料
现代控制理论(主汽温对象模型)班级:学号:姓名:目录一. 背景及模型建立1.火电厂主汽温研究背景及意义2.主汽温对象的特性3.主汽温对象的数学模型二.分析1.状态空间表达2.化为约当标准型状态空间表达式并进行分析3.系统状态空间表达式的求解4.系统的能控性和能观性5.系统的输入输出传递函数6.分析系统的开环稳定性7.闭环系统的极点配置8.全维状态观测器的设计9.带状态观测器的状态反馈控制系统的状态变量图10.带状态观测器的闭环状态反馈控制系统的分析三.结束语1.主要内容2.问题及分析3.评价一.背景及模型建立1.火电厂主汽温研究背景及意义火电厂锅炉主汽温控制决定着机组生产的经济性和安全性。
由于锅炉的蒸汽容量非常大、过热汽管道很长,主汽温调节对象往往具有大惯性和大延迟,导致锅炉主汽温控制存在很多方面的问题,影响机组的整个工作效率。
主汽温系统是表征锅炉特性的重要指标之一,主汽温的稳定对于机组的安全运行至关重要。
其重要性主要表现在以下几个方面:(1) 汽温过高会加速锅炉受热面以及蒸汽管道金属的蠕变,缩短其使用寿命。
例如,12CrMoV 钢在585℃环境下可保证其应用强度的时间约为10万小时,而在 595℃时,其保证应用强度的时间可能仅仅是 3 万小时。
而且一旦受热面严重超温,管道材料的强度将会急剧下降,最终可能会导致爆管。
再者,汽温过高也会严重影响汽轮机的汽缸、汽门、前几级喷嘴和叶片、高压缸前轴承等部件的机械强度,从而导致设备损坏或者使用年限缩短。
(2) 汽温过低,会使得机组循环热效率降低,增大煤耗。
根据理论估计可知:过热汽温每降低10℃,会使得煤耗平均增加0.2%。
同时,汽温降低还会造成汽轮机尾部的蒸汽湿度增大,其后果是,不仅汽轮机内部热效率降低,而且会加速汽轮机末几级叶片的侵蚀。
此外,汽温过低会增大汽轮机所受的轴向推力,不利于汽轮机的安全运行。
(3) 汽温变化过大会使得管材及有关部件产生疲劳,此外还将引起汽轮机汽缸的转子与汽缸的胀差变化,甚至产生剧烈振动,危及机组安全运行。
现代控制理论大作业悬挂
现代控制理论大作业-悬挂1. 引言悬挂系统是一种常见的现代控制系统应用。
它可以在车辆、建筑物等领域中发挥重要的作用。
悬挂系统的目标是提供良好的车辆行驶舒适性和稳定性。
本文将通过运用现代控制理论来设计一个悬挂系统,并讨论其性能和优化。
2. 悬挂系统的基本原理悬挂系统的基本原理是通过悬挂系统的控制,使得车辆的驾驶舒适性和稳定性得到提升。
在悬挂系统中,主要包括悬挂弹簧、减振器和控制系统等组成部分。
2.1 悬挂弹簧悬挂弹簧是悬挂系统的核心组件之一。
它可以吸收车辆行驶中的冲击和震动,减少对车辆和乘客的影响。
悬挂弹簧的刚度决定了悬挂系统的刚度和响应速度。
2.2 减振器减振器是悬挂系统的另一个重要组成部分。
它可以通过控制阻尼力来减少车辆的摇晃和倾斜,提高车辆的稳定性和舒适性。
减振器的阻尼系数决定了减振器的性能。
2.3 控制系统控制系统是悬挂系统的关键组成部分。
它可以通过传感器检测车辆的状态和路况,并通过控制器来控制悬挂弹簧和减振器的工作方式。
控制系统的设计和算法决定了悬挂系统的性能和效果。
3. 现代控制理论在悬挂系统中的应用现代控制理论在悬挂系统中有广泛的应用。
其中,常用的现代控制技术包括PID控制器、状态反馈控制、模糊控制和自适应控制等。
3.1 PID控制器PID控制器是一种经典的控制算法,它通过调节比例、积分和微分三个参数,来控制悬挂系统的动态响应。
PID控制器具有简单、稳定和易于理解的特点,被广泛应用于悬挂系统的控制。
3.2 状态反馈控制状态反馈控制是一种基于系统状态的控制方法。
它通过测量系统状态,并通过反馈控制器来调节悬挂系统的输出。
状态反馈控制可以提供更好的系统性能和稳定性,并且可以根据需要进行调整。
3.3 模糊控制模糊控制是一种基于模糊逻辑的控制方法。
它可以根据模糊规则来调节悬挂系统的输出,并考虑到多个因素的影响。
模糊控制可以提供更好的鲁棒性和适应性。
3.4 自适应控制自适应控制是一种根据系统的动态特性来调整控制参数的方法。
现代控制工程大作业
目 录题目一: (1)一、系统设计及仿真分析 ............................................... 错误!未定义书签。
1.确定被控系统状态空间体现式 (1)2.系统能控、能观性鉴别 (1)3.系统极点配置 (1)4.确定状态反馈增益矩阵 (2)5.确定输入变换线性放大器K (2)6.验证跟踪单位斜坡输入信号旳稳态误差 (2)7、运用SIMULINK 建立控制系统旳动态仿真模型 ....................错误!未定义书签。
(1)跟踪单位阶跃信号旳动态仿真分析 ..................................错误!未定义书签。
(2)跟踪单位斜坡信号旳动态仿真分析 ..................................错误!未定义书签。
二、采用全维状态观测器旳状态反馈系统....................... 错误!未定义书签。
1.配置闭环系统状态观测器极点 (4)2.确定观测器偏差反馈增益矩阵G (4)3.运用SIMULINK 建立系统动态仿真模型 (5)(1)对单位阶跃输入信号仿真分析 ..........................................错误!未定义书签。
(2)对单位斜坡输入信号仿真分析 ..........................................错误!未定义书签。
4.与采用直接状态反馈旳控制系统仿真成果比较 (7)(1)0)0(=x 且0)0(ˆ=x时旳仿真成果比较 ...............................错误!未定义书签。
(2)0)0(=x 且[]222)0(ˆ=x时旳仿真成果比较 ...................错误!未定义书签。
5.单位阶跃输入作用下状态估值误差收敛性分析 (8)(1)收敛过程分析......................................................................错误!未定义书签。
现代控制理论大作业—弹簧系统
已知由弹簧-小车组成的系统。
如下图:其中两小车质量M 1=1kg M 2=1kg ,弹簧系数K=1,u 是作用在小车M 1的外力,y 1和y 2分别是小车位移,不计任何摩擦。
由以上条件可以进行受力分析得: u + M 11y ••= K(y 1-y 2)1y ••= y 1-y 2 -uK y1+ M 22y ••=K y 2 2y ••=-2y 1+2y 2因此可以得到该系统的状态空间模型设1x = y 1, 2x =1y •, 3x = y 2,4x = 2y •,则有x •=A x +B u ,y=C x ,其中:A = --⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦0100101000012020 B=100⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥⎣⎦0 C = []1000 D= []01. 先对以上模型稳定性分析如下:M 2M 1KKy 2y 1u从图中可知矩阵A 的特征值为-1.732 ,1.732, 0,0,有正值,因此开环系统是不稳定的。
2. 判断系统可控性与可观性可见系统是可控和可观的。
3. 状态反馈控制器设计有步骤1可知系统是能控的。
因此,可以通过状态反馈来任意配置闭环系统的极点。
特别是将闭环极点配置在左半开复平面可以保证闭环系统是渐进稳定的。
因此可以设计一个状态反馈控制器u= -K x ,使得闭环系统极点是-2,-1,-1+j ,-1-j ,期望的闭环特征多项式为(λ+2)(λ+1)(λ+1+j )(λ+1-j)= 432510104λλλλ++++而状态反馈控制器所到出的闭环系统特征多项式为det{()I A BK λ--}= det{I λ- 1⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦234 0 1 0 01+k k k k 0 0 0 1-2 2 0 0} =2431k 2k k λλλλ-+432-k-2 其中K=[]2134k k k k。
由以上两个多项式的相等,可得 1k=4 2k = -5 3k= -5 4k= -5。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
目录题目一:设被控系统的传递函数为)12)(6(1)(++=s s s s G o ,现在采用带有输入变换的状态反馈控制(如教材图5-15)改善系统的性能,要求控制系统满足如下性能指标:超调量%5≤σ;调节时间5.0≤s t 秒;跟踪阶跃输入信号的稳态误差0=p e ;跟踪单位斜坡输入信号的稳态误差2.0≤v e 。
(1)根据控制系统综合指标设计状态反馈增益矩阵F 及输入变换线性放大器K ,并在SIMULINK 中建立控制系统的动态仿真模型,对系统的动态和静态性能进行仿真分析。
(2)根据控制系统的响应速度,合理确定状态观测器系统矩阵期望特征值的分布,设计对状态向量x 进行估计的全维状态观测器以实现状态反馈,并保证控制系统的性能仍满足要求;在SIMULINK 中建立采用状态观测器实现状态反馈的控制系统的动态仿真模型,对系统的动态和静态性能进行仿真分析,并与采用直接状态反馈的控制系统仿真结果比较;观察单位阶跃输入作用下状态估计值xˆ与状态真实值x 的误差的收敛过程,并研究状态估计误差收敛速度与状态观测器极点的关系及其对系统性能的影响。
(3)设系统输出量可准确测量,试设计实现状态反馈的降维观测器,在SIMULINK 中建立采用降维状态观测器实现状态反馈的控制系统的动态仿真模型,对系统的动态和静态性能进行仿真分析。
一、系统设计及仿真分析1、确定被控系统状态空间表达式被控系统传递函数为)12)(6(1)(++=s s s s G o则其能观型状态空间表达式为[]⎪⎪⎩⎪⎪⎨⎧=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=xy u x x 10000118107201000 2、系统能控、能观性判别系统能控性判别矩阵[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==1000100012B A ABB Q c 则n rankQ c ==3,所以系统完全能控;又因为其状态空间表达式为能观标准型,所以系统完全能观。
综上所述,此系统既能控又能观。
3、系统极点配置根据系统性能指标要求:%5≤σ,5.0≤s t由%521≤=--ξξπσe,可推出690107.0≥ξ,取7.0=ξ;由5.05.3≤=ns t ξω,7.0=ξ,得10≥n ω,取11=n ω。
根据经典控制理论,二阶系统闭环传递函数标准形式为2222)(n n n s s s ωξωωφ++=,将7.0=ξ,11=n ω代入式0222=++n n s s ωξω,得期望的闭环主导极点为8556.7j 7.72,1±-=*λ;选择一个期望的闭环非主导极点离虚轴为主导极点的5倍以上,即2.466.773-=⨯-=*λ。
则期望的闭环特征多项式为2.559048.8326.61))()(()(23*3*2*1*-++=---=s s s s s s s p λλλ4、确定状态反馈增益矩阵F由上已知,系统完全能控,所以可以通过状态反馈任意配置系统闭环极点。
应用MATLAB 极点配置函数求解状态反馈增益矩阵,代码如下:A=[0 0 0;1 0 -72;0 1 -18];B=[1;0;0]; P=[+j*;*;]; F=place(A,B,P)求得状态反馈增益矩阵[][]8.28883.246.43321-==f f f F5、确定输入变换线性放大器K根据题目要求,跟踪单位阶跃参考输入信号的稳态误差0)(11=+--=-BK BF A C e p ,求得2.5590=K ,则带有输入变换的状态反馈系统传递函数2.559048.8326.612.5590)()(23*+++==s s s s p K s G F6、验证跟踪单位斜坡输入信号的稳态误差v e将带有输入变换的状态反馈系统等效为具有单位负反馈的闭环系统,系统方框图如下所示:图1-1 等效单位负反馈闭环系统方框图则系统传递函数2.559048.8326.612.5590)()(1)()(23+++==+=s s s s G s G s G s F φ 求得,ss s s G 48.8326.612.5590)(23++=对于单位斜坡输入信号,2/1)(s s R =;则2.0149.0)(1)(lim)(lim 00<≈+==→→s G s sR s sE e s s v所以,所设计系统能满足要求系统性能指标。
7、利用SIMULINK 建立控制系统的动态仿真模型图1-2 带有输入变换的状态反馈系统(1)跟踪单位阶跃信号的动态仿真分析为便于输出图形数据的观察比较,在t=处输入单位阶跃响应,仿真结果如下图所示:图1-3 单位阶跃响应图1-4 单位阶跃响应局部放大图静态性能分析:由图1-3可知,系统满足跟踪阶跃输入信号的稳态误差0=p e 。
动态性能分析:从局部放大图1-4中易知,超调量%51105.1=-<σ,调节时间5.0≤s t 秒(注意:为便于观察,本实验中单位阶跃信号在t=时刻输入),系统满足设计要求。
(2)跟踪单位斜坡信号的动态仿真分析在t=0处输入单位阶跃响应,仿真结果如下图所示:图1-5 单位斜坡响应图1-6 单位斜坡响应局部放大图由图t=1s 的时刻知,跟踪单位斜坡信号的稳态误差()2.08.01=-<v e ,故系统满足设计要求。
二、采用全维状态观测器的状态反馈系统基于复合系统特征值的分离特性,只要被控系统能控能观,则用状态观测器估值形成状态反馈时,可对系统的状态反馈控制器及状态观测器分别按各自的要求进行独立设计,由(一)中已求得[]8.28883.246.43-=F 。
1、配置闭环系统状态观测器极点系统完全能观,所以闭环系统状态观测器极点可任意配置。
按通常选择观测器的响应速度比所考虑的状态反馈闭环系统快2~5倍这一选择原则,取观测器期望极点为19)7.7(5.22,1-≈-⨯=*λ,116)2.46(5.23-≈-⨯=*λ2、确定观测器偏差反馈增益矩阵G应用MATLAB 极点配置函数求解观测器偏差反馈增益矩阵G ,代码如下: A=[0 0 0;1 0 -72;0 1 -18],C=[0 0 1]; P=[-19;-19;-116]; Gt=acker(A',C',P); G=Gt'求得偏差反馈增益矩阵[][]T ==136469741876321g g g G3、利用SIMULINK 建立系统动态仿真模型观测器的状态方程为:u y x Bu Gy x GC A x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=++-=001136469741876ˆ154104769014187600ˆ)(ˆ利用SIMULINK 建立系统动态仿真模型如下:图1-7 采用全维状态观测器的状态反馈系统(注:Scope1可观测误差的收敛过程)(1)对单位阶跃输入信号仿真分析在t=时刻输入单位阶跃信号,对于0)0(=x 和0)0(ˆ=x的情况下仿真结果如下图所示:图1-8 0)0(=x 且0)0(ˆ=x时单位阶跃响应图1-9 0)0(=x 且0)0(ˆ=x时单位阶跃响应局部放大图 静态性能分析:由图1-8可知,系统满足跟踪阶跃输入信号的稳态误差0=p e 。
动态性能分析:从局部放大图1-9中易知,超调量%51105.1=-<σ,调节时间5.0≤s t 秒(注意:为便于观察,本实验中单位阶跃信号在t=时刻输入),系统满足设计要求。
(2)对单位斜坡输入信号仿真分析在t=0时刻输入单位斜坡信号,对于各状态变量初值0)0(=x 且0)0(ˆ=x 情况下的仿真结果如下图所示:图1-100)0(=x 且0)0(ˆ=x时单位斜坡信号响应图1-110)0(=x 且0)0(ˆ=x时单位斜坡信号响应局部放大图 由图t=1s 的时刻知,跟踪单位斜坡信号的稳态误差()2.08.01=-<v e ,故系统满足设计要求。
4、与采用直接状态反馈的控制系统仿真结果比较(1)0)0(=x 且0)0(ˆ=x 时的仿真结果比较(a )直接状态反馈的单位阶跃响应 (b )采用状态观测器的单位阶跃响应(c )直接状态反馈的单位斜坡响应 (d )采用状态观测器的单位斜坡响应图1-12 0)0(=x 且0)0(ˆ=x的仿真结果比较 从以上仿真结果图可以看出,在各状态变量0)0(=x 且0)0(ˆ=x的情况下,直接状态反馈与采用状态观测的反馈系统仿真结果几乎一样。
(2) 0)0(=x 且 []222)0(ˆ=x时的仿真结果比较(a) 直接状态反馈的单位阶跃响应 (b) 采用状态观测器的单位阶跃响应(c) 直接状态反馈的单位斜坡响应 (d) 采用状态观测器的单位斜坡响应图1-13 0)0(=x 且[]222)0(ˆ=x的仿真结果比较 由图1-13可知,尽管初始时刻0)0(=x 和[]222)0(ˆ=x 存在差异,但随着时间的推移,初始值不同情况下的响应也将趋于同样的稳定状态。
这其实与后面讨论的估值误差的收敛性质有关。
对于复杂的实际工程应用,一般不可能采用直接状态反馈,而是采用观测器进行状态反馈以获得控制精度并且降低成本。
5、单位阶跃输入作用下状态估值误差收敛性分析(1)收敛过程分析对0)0(ˆ)0(==xx 和0)0(=x ,[]222)0(ˆ=x 两种情况下的系统进行收敛性仿真分析,结果如下所示:图1-14 0)0(ˆ)0(==xx 时误差收敛情况图1-15 0)0(=x ,[]222)0(ˆ=x时误差收敛情况 由图可知,当初始值相同,即0)0(ˆ)0(==xx 时,状态估值没有误差;当初始值不同,即0)0(=x ,[]222)0(ˆ=x时,开始一段时期存在估值误差,但是,随着时间的推移)(ˆ)(t x t x =,即估值误差收敛到0。
从理论上讲,状态观测器之所以出现估值误差是因为观测器与系统初始状态在实际工程中难以做到完全一致所造成的。
估值误差为())]0()0(ˆ[)(x xe t t GC A x -=∆- 根据上式,当初始时刻00=t ,)0(ˆ)0(xx =时,显然根据上式0)(=∆t x ;当0)0(=x ,[]222)0(ˆ=x时,由于我们通过选择输出偏差反馈增益矩阵G 使GC A -的所有特征值均位于复平面的左半开平面,尽管初始时刻00=t 时,0)0(=x ,[]222)0(ˆ=x存在差异,但观测器的状态)(ˆt x仍将以一定精度和速度渐进逼近系统的实际状态)(t x 。
因此,仿真结果与理论分析相符。
(2)收敛速度与状态观测器极点的关系及对系统性能的影响 现重新选取全维观测器期望极点为35)7.7(5.42,1-≈-⨯=*λ,208)2.46(5.43-≈-⨯=*λ应用MATLAB 极点配置函数求解观测器偏差反馈增益矩阵G ,代码如下: A=[0 0 0;1 0 -72;0 1 -18],C=[0 0 1]; P=[-35;-35;-208]; Gt=acker(A',C',P); G=Gt'求得偏差反馈增益矩阵[][]T==260 15713 2548003211g g g G将1G 代入采用全维状态观测器的状态反馈系统中,并在SIMULINK 中进行动态仿真,当0)0(=x , []222)0(ˆ=x时,各状态变量的收敛情况如下图所示:图1-16 208,3532,1-=-=**λλ时的各状态变量的收敛情况由(1)中已知116,1932,1-=-=**λλ情况下的误差收敛情况如下图所示:图1-17 116,1932,1-=-=**λλ时的各状态变量的收敛情况从图1-16和图1-17中各状态变量的收敛情况可以看出观测器极点取不同值时,)(ˆt x的误差收敛情况。