资源描述：
《现代控制理论2》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《现代控制理论》实验报告实验名称倒立摆系统状态反馈设计与仿真课程名称现代控制理论院系台P:专业班级：学生姓名：学号：同组人：实验台号:指导教师：成绩实验日期：年月曰实验二倒立摆系统状态反馈设计与仿真一、实验目的1、利用MATLAB分析倒立摆系统的能控性和能观性。2、对倒立摆线性化系统进行状态反馈设计并在simulink中进行仿真，比较线性模型与非线性模型在不同初始值卜的响应。3、对倒立摆线性化系统设计基于观测器的状态反馈，迕simulink中进行仿真，并比较线性模型与非线性模型在不同初始值卜*的响应；与
2、原状态反馈控制进行比较。4、运用现代控制理论解决实际问题。二、实验类型设计性三、实验设备计算机四、实验原理当系统状态完全能控时，可以采用极点配置方法设计状态反馈《=-似使得闭环极点等于指定的极点，从而达到期望的闭环性能。当状态不可测时，可以采用状态观测器利用系统过去输出及输出估计当前状态，从而实现基于状态观测器的状态反馈。五、实验内容和要求（一）内容取似=l，m=0.1，g=10，/=1,倒立摆系统的线性化状态空间模型为:II070111ooOoooO1其中，X,=y,x2=y,x3=0,x4=0试利用
3、matUb进行以不分析：(1)分析系统的可控性,可观性;⑵对该系统设计状态反馈控制，对比倒立摆非线性模型与线性化模型在初始摆角小扰动及大扰动两种情况下的响应曲线，并分析结果；(3)对该系统设计基于观测器的状态反馈控制，对比倒立摆非线性模型与线性化模型在初始摆角小扰动及大扰动两种情况下的响应曲线，并分析结果。(二)要求1、预习现代控制理论知识点，编制相应实现的MATLAB程序。2、预习倒立摆系统的原理及结构。3、预习使用simulink工具箱。[实验步骤]:(1)能控性分析：A=[0100;00-10;0
4、001:00110];B=[0;l;0;-l];n=rank(ctrb(A,B))n=4系统满秩，说明系统为完全能控系统，可以对其进行闭环极点的任意配置。能观性分析：•I234&^v&&ooIL~-a)如果系统输出为小车的位移，则或者说C=b°°°]进行能观性分析A=[0100;00-10;0001;00110];c=[i000];n二rank（obsv（A，0）系统满秩，说明系统为完全能观系统，可以设计观测器，且观测器极点可任意配ab)如果系统输出为倒立摆的摆角，则■&4」,或者说c=[001o]进行
5、能观性分析A=[0100;00-10;0001;00110];C=[0010];n=rank(obsv(A,C))n=2系统不满秩，说明系统不是完全能观的，观测器极点不可任意配置。(2)当倒立摆系统所有状态都能测M吋，可以采用状态反馈控制达到期增性能。取闭环期望极点为：51,=-6,s2=-6.5,53=—7,54=—7.5在matlab中输入A=[0100;00-10;0001;00110];P=[-6-6.5-7-7.5];Replace(A,B,P)得到计算结果为：K=-204.7500-122.
6、1750-488.5000一149.1750•线性系统闭环仿真在状态反馈么=一&&下，倒立摆线性化模型在simulink的闭环仿真图如下:a）线性模型在无外部扰动、初始摆角为1时,响应如下图所示:b）线性模型在外部脉冲输入f=lN、初始摆角为吋，响应如卜阁所示:Time(sec.)Time(sec.)•非线性模型闭环仿真在状态反馈&=一夂&*下，倒立摆非线性模型在simulink的闭环仿真图如下:K1a）非线性模型在无外部扰动、初始摆角为1时，响应如下图虚线所示:0.6■02i00.20.40.60.8
7、11.21.41.61.82与线性模型响应（实线）相比，相差不大。b）非线性模型在外部脉冲输入f=lN、初始®角为0时，响!、V:如下图所示:0.02'0-0.02-0.041111111111100.20.40.60.811.21.41.61.82Time（sec.）与线性模型响应（实线）相比，十分接近。C）非线性模型在无外部扰动、初始摆角为2（r时，响应如K阁虚线所示：0.5X0-051111111111100.20.40.60.811.21.41.61.82与线性模型响应（实线）相比，超调更人。d
8、）非线性模型在外部脉冲输入f=10N、初始摆角为&时，响皮如下图所示:与线性模型响应（实线）相比，超调更大。e）当初始角为时，非线性模型不稳定。（3）观测器设计与仿真如果仅小车位移能被测朵，则OOrnu-II或者说00o]由（1）知（A，C）完全能观，从而可以任意配置观测器极点。取闭环期望极点为：5,=—20,59=—21,53=—22,54=—23在matlab中输入A=[0100;00-10;0001;00110];C=[l000];P
