现代控制实验报告
现代控制理论实验(1)
《现代控制理论实验》一．实验目的1、掌握组态王软件开发工程的步骤2、学会设备通讯的连接方法3、开发一个计算机-智能仪表控制系统的单回路控制系统4、绘制一条温度控制曲线二．实验设备计算机，仪表结构实验控制箱，电加热炉，组态王软件三．实验内容与实验步骤第一步：建立一个新工程在工程浏览器中双击“新建工程”图标，建立一个新工程，进入“组态王”工程浏览器。
第二步：设备连接本次实验硬件设备选用理化公司生产的CD901数字温度控制器。
点击工程浏览器树形菜单下面的----设备----COM1---右面的新建图标，选择“智能仪表-理化-CD901-串口”，定义地址为10或11（右侧回路地址为0010，左侧回路地址为0011）。
双击工程浏览器左侧树形菜单中的“COM1”，对智能仪表CD901进行通讯配置：波特率：9600 数据位：8奇偶校验：无停止位： 1通讯超时：3000 豪秒通讯方式：RS232（将通讯配置截图粘贴于此处）第三步：图形界面的设计双击我们设计的画面或新建一个画面，进入“组态王”开发系统，这时开始设计界面。
界面应包含一条温度曲线用以采集并绘出实时温度，两个按钮用以启动和停止程序，两处字符显示用以表示给定温度值和实时温度值，另外需在界面显示编程人的班级、学号、姓名等信息。
用到的控件应包括：一个“X-Y轴曲线”控件，2个矩形、2个按钮，以及必要的文本。
如图所示：（将设计界面截图粘贴于此处）第四步：构造数据库在这个工程中需要定义6个变量，变量属性如下表：其中PV 表示温度的采样值，SV表示温度的给定值，AUTO_TUNING表示自动切换值，H_P表示加热比例带，on 用来设置开始标志，runtime表示当前时间。
第五步：动画连接文本：双击SV对应的文本“sv”，在“模拟量输出”和“模拟量输入”选项卡的“变量名”处选择\\本站点\SV，“模拟量输出”选项的整数位、小数位分别选择3位、0位，“模拟量输入”选项的最大值、最小值处分别选择300、0；双击PV对应的文本“pv”，在“模拟值输出”选项卡的“变量名”处选择“\\本站点\PV”，“模拟量输出”选项的整数位、小数位分别选择3位、2位。
现代控制理论基础实验
现代控制理论基础实验一、 实验目的1. 熟悉MATLAB 的编程以及SIMULINK 仿真工具的使用。
2. 通过实验掌握极点配置及设计状态反馈控制器K 的方法。
3. 深入了解电动机速度控制系统的综合控制方法。
二、 实验内容电动机速度控制系统，设计状态反馈控制器K ，使得系统跟踪单位阶跃指令时无静态误差，超调量s t s 1%,5%<≤σ。
要求写出详细的设计步骤，给出仿真设计系统原理框图，给出仿真的输出波形图和误差波形图。
三、 实验原理控制系统最基本的结构形式是由受控系统和实现反馈控制规律的反馈环节所构成的反馈控制系统。
现代控制理论中，存在两种基本的反馈形式，即状态反馈和输出反馈。
实际情况中，状态反馈具有更好的特性和适应性。
系统动力学的各种特性或各种品质指标，在很大程度上是由系统的极点决定的。
所谓极点配置问题，就是通过状态反馈矩阵K 的选择，使闭环系统的极点，恰好处于所希望的位置。
从线性定常系统运动分析可知，如时域中超调量、过渡过程时间及频域中增益稳定裕度、相位稳定裕度，都被认为等价于系统极点位置，相应综合问题可视为极点配置问题。
四、系统设计1、根据图1计算出电机控制系统的传递函数，并化为状态空间模型图一 受控系统方块图（简化）)(1)10s(0.4s )(5.0]10.05s )(U 3.0)(U 4.0[s Y s Y s s =+-+-可求得受控系统的传递函数：5.2125.502^5.223^5.01.0)()()(++++==s s s s s U s Y s G 系统有一个零点z 1 = -5；用求根函数roots()计算函数极点 >> C=[1 22.5 50.125 2.5];>> roots(C) ans =-20.0000 -2.4490 -0.0510由题意设状态分别为：系统simulink 仿结构如下图二 受控系统simulink 仿真结构图⎪⎪⎪⎩⎪⎪⎪⎨⎧=+--=+=233121*10114.01*]*5.0)4.0[(1005.03.0x s x s x x u x u s x 化为标准形式可得：⎪⎪⎪⎭⎫ ⎝⎛'3'2'1x x x =⎪⎪⎪⎭⎫ ⎝⎛----01.0025.15.25.20020⎪⎪⎪⎭⎫ ⎝⎛321x x x +⎪⎪⎪⎭⎫ ⎝⎛016y=()100 x系统的性能指标：调节时间t s = 76.6s ，上升时间t r = 42.8s ，超调量0%=σ2、确定希望的极点希望的极点数为3，由系统要求超调量低于5%，ts 小于1秒选其中一对为主导极点1s 和2s ，另一个为远极点，并且认为系统的性能主要是由主导极点决定的，远极点所产生的影响很小，可以忽略不计。
现代控制理论实训报告
一、前言随着科技的飞速发展，自动化、智能化已成为现代工业生产的重要特征。
为了更好地掌握现代控制理论，提高自己的实践能力，我参加了现代控制理论实训课程。
本次实训以状态空间法为基础，研究多输入-多输出、时变、非线性一类控制系统的分析与设计问题。
通过本次实训，我对现代控制理论有了更深入的了解，以下是对本次实训的总结。
二、实训目的1. 巩固现代控制理论基础知识，提高对控制系统的分析、设计和调试能力。
2. 熟悉现代控制理论在工程中的应用，培养解决实际问题的能力。
3. 提高团队合作意识，锻炼动手能力和沟通能力。
三、实训内容1. 状态空间法的基本概念：状态空间法是现代控制理论的核心内容，通过建立状态方程和输出方程，描述系统的动态特性。
2. 状态空间法的基本方法：包括状态空间方程的建立、状态转移矩阵的求解、可控性和可观测性分析、状态反馈和观测器设计等。
3. 控制系统的仿真与实现：利用MATLAB等仿真软件，对所设计的控制系统进行仿真，验证其性能。
4. 实际控制系统的分析：分析实际控制系统中的控制对象、控制器和被控量，设计合适的控制策略。
四、实训过程1. 理论学习：首先，我对现代控制理论的相关知识进行了复习，包括状态空间法、线性系统、非线性系统等。
2. 实验准备：根据实训要求，我选择了合适的实验设备和软件，包括MATLAB、控制系统实验箱等。
3. 实验操作：在实验过程中，我按照以下步骤进行操作：（1）根据实验要求，建立控制系统的状态空间方程。
（2）求解状态转移矩阵，并进行可控性和可观测性分析。
（3）设计状态反馈和观测器，优化控制系统性能。
（4）利用MATLAB进行仿真，观察控制系统动态特性。
（5）根据仿真结果，调整控制器参数，提高控制系统性能。
4. 结果分析：通过对仿真结果的分析，我对所设计的控制系统进行了评估，并总结经验教训。
五、实训成果1. 掌握了现代控制理论的基本概念和方法。
2. 提高了控制系统分析与设计能力，能够独立完成实际控制系统的设计。
现代控制理论实验报告
现代控制理论实验报告实验一系统能控性与能观性分析一、实验目的1.理解系统的能控和可观性。
二、实验设备1.THBCC-1型信号与系统·控制理论及计算机控制技术实验平台；三、实验容二阶系统能控性和能观性的分析四、实验原理系统的能控性是指输入信号u对各状态变量x的控制能力，如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间把系统所有的状态引向状态空间的坐标原点，则称系统是能控的。
对于图21-1所示的电路系统，设iL和uc分别为系统的两个状态变量，如果电桥中则输入电压ur能控制iL和uc状态变量的变化，此时，状态是能控的。
反之，当时，电桥中的A点和B点的电位始终相等，因而uc不受输入ur的控制，ur只能改变iL的大小，故系统不能控。
系统的能观性是指由系统的输出量确定所有初始状态的能力，如果在有限的时间根据系统的输出能唯一地确定系统的初始状态，则称系统能观。
为了说明图21-1所示电路的能观性，分别列出电桥不平衡和平衡时的状态空间表达式：平衡时：由式（2）可知，状态变量iL和uc没有耦合关系，外施信号u只能控制iL的变化，不会改变uc的大小，所以uc不能控。
基于输出是uc，而uc与iL无关连，即输出uc中不含有iL的信息，因此对uc的检测不能确定iL。
反之式（1）中iL与uc有耦合关系，即ur的改变将同时控制iL和uc的大小。
由于iL与uc的耦合关系，因而输出uc的检测，能得到iL 的信息，即根据uc的观测能确定iL（ω）五、实验步骤1．用2号导线将该单元中的一端接到阶跃信号发生器中输出2上，另一端接到地上。
将阶跃信号发生器选择负输出。
2．将短路帽接到2K处，调节RP2，将Uab和Ucd的数据填在下面的表格中。
然后将阶跃信号发生器选择正输出使调节RP1，记录Uab和Ucd。
此时为非能控系统，Uab和Ucd没有关系（Ucd始终为0）。
3．将短路帽分别接到1K、3K处，重复上面的实验。
现代控制理论实验报告
现代控制理论实验指导书实验一：线性系统状态空间分析1、模型转换图1、模型转换示意图及所用命令传递函数一般形式：)()(11101110n m a s a s a s a b s b s b s b s G n n n n m m m m ≤++++++++=----K KMATLAB 表示为：G=tf(num,den)，其中num,den 分别是上式中分子，分母系数矩阵。
零极点形式：∏∏==--=n i j mi i ps z s K s G 11)()()( MATLAB 表示为：G=zpk(Z,P ,K)，其中 Z ，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。
传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)；状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu表示对系统的第iu个输入量求传递函数；对单输入iu为1;验证教材P438页的例9-6。
求P512的9-6题的状态空间描述。
>> A=[0 1;0 -2];>> B=[1 0;0 1];>> C=[1 0;0 1];>> D=[0 0;0 0];>> [NUM,DEN] = ss2tf(A,B,C,D,1)NUM =0 1 20 0 0DEN =1 2 0>> [NUM,DEN] = ss2tf(A,B,C,D,2)NUM =0 0 10 1 0DEN =1 2 0给出的结果是正确的，是没有约分过的形式P512 9-6>> [A,B,C,D]=tf2ss([1 6 8],[1 4 3])A =-4 -31 0B =1C =2 5D =12、状态方程求解单位阶跃输入作用下的状态响应：G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x). 零输入响应[y,t,x]=initial(G,x0)其中，x0为状态初值。
现代控制理论实验
现代控制理论实验引言现代控制理论是在工程控制领域中发展起来的一种理论体系，其应用范围非常广泛。
为了帮助学生更好地理解和掌握现代控制理论，学校开设了现代控制理论实验课程。
该实验课程旨在通过具体的实验操作，帮助学生巩固理论知识，培养实际操作能力，并能应用现代控制理论解决实际问题。
本文将介绍现代控制理论实验的内容、目的、实验装置和实验步骤。
实验内容现代控制理论实验主要包括以下内容： 1. PID控制器的设计与实现：通过调节比例、积分和微分参数，设计一个PID 控制器，并将其实现在实验装置上，观察控制效果。
2. 状态反馈控制器的设计与实现：利用状态观测器和状态反馈器，设计一个状态反馈控制器，并将其实现在实验装置上，观察控制效果。
3. 频域方法的应用：通过频域分析方法，设计一个控制器，使得实验装置的频率响应满足特定要求。
4. 鲁棒控制方法的应用：利用鲁棒控制方法设计一个控制器，能够在系统参数变化时保持系统的稳定性和性能。
实验目的现代控制理论实验的主要目的是培养学生的实践能力和问题解决能力。
具体目标包括： 1. 理解现代控制理论的基本原理与方法； 2. 掌握现代控制理论的实验操作技巧； 3. 理解研究现代控制理论的方法和途径； 4. 能够设计、实现和调试现代控制器，并分析控制效果； 5. 学会通过实验结果验证和改进控制算法。
实验装置现代控制理论实验装置主要包括：电机系统、传感器、数据采集卡、计算机控制软件和控制器实现装置。
电机系统是实验装置的核心部件，它模拟了真实的控制对象。
传感器用于感知电机系统的转速、位置或其他关键参数。
数据采集卡负责将传感器采集到的数据传输给计算机进行处理。
计算机控制软件包括了实验的开发工具和界面，可以实时控制电机系统并显示实验结果。
控制器实现装置是通过软件或硬件方式实现控制器，在实验中使用。
实验步骤本节将介绍现代控制理论实验的基本步骤。
具体步骤如下：步骤一：系统建模与参数辨识首先需要对实验装置进行数学建模，并通过实验数据对模型参数进行辨识。
现代控制理论实验
现代控制理论实验华北电力大学实验报告||实验名称状态空间模型分析课程名称现代控制理论基础||专业班级：自动化1203 学生姓名：孟令虎学号：201209020216 成绩：指导教师：刘鑫屏老师实验日期： 2015.4.24一、实验目的l.加强对现代控制理论相关知识的理解；2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验仪器与软件 1. MATLAB7.6 环境三、实验内容1、 模型转换例 1.把传递函数模型转化为状态空间模型3248G s =81912s s s s ++++（）。
解：程序如下num=[4 8]; den=[1 8 19 12];[A,B,C,D]=tf2ss(num,den); G=ss(A,B,C,D) 运行结果： A =-8 -19 -12 1 0 0 0 1 0 B = 1 0 0 C =0 4 8 D =0 结果为112233-8 -19 -1211 0 010 1 00x x x x u x x ∙∙∙⎡⎤⎢⎥⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦，[]1230 4 8x y x x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦例2.把状态空间模型转化为传递函数模型A=0 1 00 0 1-6 -11 -6⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B=001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C= []2 3 0 D=0。
解：程序如下：clearA=[0 1 0;0 0 1;-6 -11 -6]; B=[0;0;1]; C=[3 2 0]; D=0; iu=1;[num,den] = ss2tf(A,B,C,D,iu); sys=tf(num,den) 运行结果为：Transfer function: 2 s + 3---------------------- s^3 + 6 s^2 + 11 s + 62、 状态方程状态解和输出解例1.单位阶跃输入作用下的状态响应A=0 1 00 0 1-6 -11 -6⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B=001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C=[]2 3 0 D=0。
现代电气控制技术实训报告范文（通用7篇）
现代电气控制技术实训报告现代电气控制技术实训报告范文（通用7篇）随着个人的文明素养不断提升，大家逐渐认识到报告的重要性，报告中提到的所有信息应该是准确无误的。
你知道怎样写报告才能写的好吗？以下是小编帮大家整理的现代电气控制技术实训报告范文（通用7篇），仅供参考，欢迎大家阅读。
现代电气控制技术实训报告1一、实训目的本实训环节是集知识、素质和技能训练于一体的应用型课程。
它在运用相关电气理论基础在上的同时又对电工基本素质和技能进行了综合运用培养，通过实训，使学生具备初级电工的基本能力，能满足企业对初级电工的综合需要。
主要内容包括：安全用电常识、常用电工材料的认识和选用、电工基本操作工艺、电气照明与内线工程、常用电工仪表的使用等。
二、实训时间第x周——第x周。
三、实训地点实验楼电工实训室四、实训设备及器材网孔板，电工工具，电度表，功率表，万用表，电压表，兆欧表，电流表，空气开关，熔断器，各种导线，硬塑料管，开关，插座，灯若干。
五、实训内容及要求电工实训一：安全用电常识电工安全操作的各项规定是每一名电工必须遵守的规章制度，它规定对电工的最基本的要求。
电工生产岗位责任制规范了电工的工作范围，是确保电工工艺得以贯彻执行的重要条件。
主要进行电工安全操作规程、预防触电及触电急救基本常识、防雷保护以及电气火灾的扑救等内容。
触电的原因：（1）电气设备的安装过于简陋，不符合安全要求。
（2）电气设备老化，破损严重，维修维护不及时。
（3）作业时不严谨，不遵守电工安全操作规程或粗心大意。
（4）缺乏安全知识电流对人体的伤害：触电对人体的伤害只要是电击和电伤。
点击是触电者直接触了设备的带电部分，电流通过了人的身体，当电流达到一定的数值后，就会将人体击倒；电伤是指触电后皮肤的局部创伤，由于电流的热效应，化学反应，机械效应以及电流的作用下，使熔化和蒸发的金属微粒袭击人体皮肤而遭受灼伤。
影响触电后果的因素：（1）电流强度：一般50ml持续1秒以上，致命。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验二状态空间标准形与控制系统的运动分析一、实验目的1.掌握线性系统的对角线标准形、约当标准形、能控标准形和能观测标准型的表示及相应变换阵的求解。
深入理解状态空间模型的相关理论。
2.掌握利用Matlab进行矩阵指数函数的数值计算和符号计算方法；对定常连续系统和定常离散系统的状态空间模型进行求解，分析其运动规律；对连续系统进行离散化。
二、实验内容及步骤1.将实验一的第2题用对角标准型实现a1 = 5; a0 = 6;b2 = 1; b1 = 2; b0 = 1;c0 = b2;c1 = b1 - a1*c0;c2 = b0 - a0*c0 - a1*c1;G = [0 1;-a0 -a1];H = [c1;c2];C = [1 0];D = [c0];[Q,d] = eig(G);P = inv(Q);Gb = P*G*QHb = P*HCb = C*QDb = D运行结果Gb = -2.0000 0.0000-0.0000 -3.0000Hb =2.236112.6491Cb =0.4472 -0.3162Db =12.已知系统的动态方程如下：0101001061160X X∙⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦，=y[]110X1）求对角标准型实现，并写出实现变换的非奇异阵和变换关系。
A=[0 1 0;0 0 1;-6 -11 -6];B=[1;0;0];C=[1 1 0];D=0;[Q,D]=eig(A)P=inv(Q)%[con_ss,T]=canon(A,B,C,D,'modal')%[con_ss1,T1]=canon(A,B,C,D,'companio n')A1=P*A*QB1=P*B运行结果Q =-0.5774 0.2182 -0.10480.5774 -0.4364 0.3145-0.5774 0.8729 -0.9435D =-1.0000 0 00 -2.0000 00 0 -3.0000P =-5.1962 -4.3301 -0.8660-13.7477 -18.3303 -4.5826-9.5394 -14.3091 -4.7697A1 =-1.0000 -0.0000 0.0000-0.0000 -2.0000 0-0.0000 0.0000 -3.0000B1 =-5.1962-13.7477-9.53942）求可控准型实现（能控规范I形），并写出实现变换的非奇异阵和变换关系。
A=[0 1 0;0 0 1;-6 -11 -6];B=[1;0;0];C=[1 1 0];D=0;Qc=ctrb(A,B)M=eig(A)Q=Qc*[M(2) M(3) 1;M(3) 1 0;1 0 0]P=inv(Q)Ab=P*A*Q,Bb=P*B,Cb=C*Q运行结果Qc = 1 0 00 0 -60 -6 36M =-1.0000-2.0000-3.0000Q =-2.0000 -3.0000 1.0000-6.0000 0 054.0000 -6.0000 0P =0.0000 -0.1667 -0.00000 -1.5000 -0.16671.0000 -4.8333 -0.5000Ab =-9.0000 1.0000 0.0000-40.0000 0 1.0000-144.0000 2.0000 3.0000Bb =0.00001.0000Cb =-8.0000 -3.0000 1.00003.试在Matlab中计算如下矩阵的特征值和广义特征向量。
436102111A---⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦A=[-4 -3 -6;1 0 2;1 1 1]; M=eig(A) [V ,J]=jordan(A) 运行结果 M = -1.0000 -1.0000 -1.0000 V =-3 1 -2 1 0 0 1 0 1 J =-1 1 0 0 -1 0 0 0 -14. 试在Matlab 中将如下状态空间模型变换为约旦规范形。
010*********[100]⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦=x x u y xA = [0 1 0;0 0 1;-6 -11 -6];B = [1;0;0];C = [1 0 0]; [Q,J] = jordan(A); P = inv(Q); Ab = J Bb = P*B Cb = C*P运行结果 Ab =-3 0 0 0 -2 0 0 0 -1 Bb = 9.0000 -12.0000 3.0000 Cb =9.0000 13.5000 4.50005. 将下列状态方程化为约当型，并写出实现变换的非奇异阵和变换关系。
A = [3 -1 1 1 0 0;1 1 -1 -1 0 0;0 0 2 01 1;0 0 02 -1 -1;0 0 0 0 1 1; 0 0 0 0 1 1];B = [1 0;-1 1;2 1;0 -1;0 2;1 0]; [Q,J] = jordan(A) P = inv(Q) Ab = J Bb = P*B 运行结果 Q =2.0000 2.0000 1.0000 0 0 02.0000 0 0A =3 -1 1 1 0 0 1 1 -1 -1 0 0 0 0 2 0 1 1 0 0 0 2 -1 -1 0 0 0 0 1 1 0 0 0 0 1 10 0 0 0 0 1.00000 1.0000 00 0 0 0 -1.0000 0 0 0 0 0.5000 0 0.5000 0 0 0 -0.5000 0 0.5000J =2 1 0 0 0 0 0 2 1 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 0 0 0 0 0 2 P = 0 0.5000 0 0 0 0 0.5000 -0.5000 -0.5000 -0.5000 0 0 0 0 1.0000 1.0000 0 0 0 0 0 0 1.0000 -1.0000 0 0 0 -1.0000 0 0 0 0 0 0 1.0000 1.0000 Ab = 2 1 0 0 0 0 0 2 1 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 0 0 0 0 0 2 Bb = -0.5000 0.5000 0 -0.5000 2.0000 0 -1.0000 2.0000 0 1.00001.00002.0000 三、思考题 将模型实现为某种标准型（可控标准型，可观测标准型和Jondan 标准型）的条件是什么？ 答：将模型实现为某种标准型时，一定要判断模型是否符合某种标准型的要求。
四、实验总结 本次实验的主要目的是：掌握线性系统的对角线标准形、约当标准形、能控标准形和能观测标准型的表示及相应变换阵的求解；掌握利用Matlab 进行矩阵指数函数的数值计算和符号计算方法；对定常连续系统和定常离散系统的状态空间模型进行求解，分析其运动规律；对连续系统进行离散化。
通过实验编写程序实现了预期结果，并进一步巩固了上课时所学习的理论知识及变换方法。
实验三、四 系统的可控性和可观性一、实验目的系统的能控性和能观性关系到系统的极点配置法设计，和最优控制。
通过本实验，掌握判断系统能控性能观法。
二、实验内容和步骤1检验系统1)0,1001,1111,2411=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡-=D C B A2)132210201101311120101y ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎡⎤=⎢⎥-⎣⎦x x u 的能控性和能观性。
%%1.1 A=[1 1;4 -2]; B=[1 -1;1 -1]; C=[1 0; 0 1]; D=0; Co=ctrb(A,B)if(rank(Co)==length(A))disp('系统能控'); elsedisp('系统不能控'); end N=obsv(A,C)if(rank(N)==length(A)) disp('系统能观'); elsedisp('系统不能观');end 运行结果 Co =1 -12 -21 -12 -2 系统不能控 N =1 0 0 1 1 1 4 -2 系统能观 %%1.2A=[1 3 2;0 2 0; 0 1 3]; B=[2 1;1 1;-1 -1]; C=[1 2 0;-1 0 1]; D=0; Co=ctrb(A,B)if(rank(Co)==length(A)) disp('系统能控'); elsedisp('系统不能控'); end N=obsv(A,C)if(rank(N)==length(A)) disp('系统能观'); elsedisp('系统不能观'); end 运行结果 Co =2 13 2 54 1 1 2 2 4 4 -1 -1 -2 -2 -4 -4 系统不能控 N =1 2 0 -1 0 1 1 7 2 -1 -2 1 1 19 8 -1 -6 1 系统能观2试在Matlab 中判定如下离散系统的状态能控性。
1001(1)022()2()1101k k k ⎡⎤⎡⎤⎢⎥⎢⎥+=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦x x u 011(1)()()000k k k ⎡⎤⎡⎤+=+⎢⎥⎢⎥⎣⎦⎣⎦x x u%%2.1G=[1 0 0;0 2 -2;-1 1 0];H=[1;2;1]; Co=ctrb(G,H)if(rank(Co)==length(G)) disp('系统能控'); elsedisp('系统不能控'); End 运行结果 Co =1 1 12 2 2 1 1 1 系统不能控3用格拉姆矩阵判据判断下面系统的可控性与可观性。
12320110()1011[110]k y x---⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦=x x u A = [-1 -2 -3;0 -1 1;1 0 -1]; B = [2;0;1]; C = [1 1 0]; D = 0;sys=ss(A,B,C,D); Wc=gram(sys,'c'); if rank(Wc) == length(A) disp('系统能控') elsedisp('系统不能控') endWo = gram(sys,'o'); if rank(Wo) == length(A) disp('系统能观') elsedisp('系统不能观') end 系统能控 系统能观4 将下面系统化为可控与可观标准型。