《现代控制理论 刘豹 .pptx》由会员分享,可在线阅读,更多相关《现代控制理论 刘豹 .pptx(41页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第二章第二章 控制系统状态空间表达式的解控制系统状态空间表达式的解主要内容 线性定常连续时间齐次状态方程的解 状态转移阵及其求解 线性定常连续时间非齐次状态方程的解线性连续时间系统的离散化线性离散系统状态空间表达式的解 第1页/共41页第二章第二章 控制系统状态空间表达式的解控制系统状态空间表达式的解2-1 线性定常齐次状态方程的解系统的唯一解证明:假设x(t)为矢量幂级数形式第2页/共41页2-2-1 1 线性定常齐次状态方程的解线性定常齐次状态方程的解同幂次项系数相等证明:假设x(t)为矢量幂级数形式第3页/共41页2-2-1 1 线性定常齐次状态方程的解线性定常齐次状态方程的解方程的解可
2、写为t=0,b0=x0第4页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩状态转移矩阵阵状态转移矩阵状态转移阵是nxn的状态转移阵的概念、性质、求法第5页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的性质第6页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的性质第7页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的性质性质1分步走跃步走第8页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的性质性质2性质3性质4从定义得
3、证从定义得证第9页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的性质问题1:已知某状态转移阵如下,求其逆阵问题2:已知某状态转移阵如上,求其A阵提示:第10页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵特殊矩阵的状态转移阵1 A为对角型第11页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵第12页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵特殊矩阵的状态转移阵2 可化为对角型第13页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵证
4、明以上结论第14页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的计算方法1方法2第15页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵状态转移阵的计算例题:方法1解第16页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵例题:第17页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵例题:方法2特征值:第18页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵例题第19页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩
5、阵例题第20页/共41页2-2 2-2 矩阵指数函数矩阵指数函数-状态转移矩阵状态转移矩阵凯莱-哈密顿定理:方阵A满足其自身的特征方程,即它会产生什么结论?第21页/共41页2-3 2-3 线性定常非齐次方程的解线性定常非齐次方程的解证明第22页/共41页2-3 2-3 线性定常非齐次方程的解线性定常非齐次方程的解例题2-8第23页/共41页2-3 2-3 线性定常非齐次方程的解线性定常非齐次方程的解例题2-8第24页/共41页2-4 2-4 连续时间状态方程的离散化连续时间状态方程的离散化连续时间系统 K-KT u(K)=constant第25页/共41页2-4 2-4 连续时间状态方程的离
6、散化连续时间状态方程的离散化连续时间系统 近似离散化第26页/共41页2-4 2-4 连续时间状态方程的离散化连续时间状态方程的离散化近似离散化第27页/共41页2-4 2-4 连续时间状态方程的离散化连续时间状态方程的离散化例题:2-13 离散化状态方程分析选择不同的采样周期T.的影响第28页/共41页2-4 2-4 连续时间状态方程的离散化连续时间状态方程的离散化第29页/共41页2-4 2-4 离散时间系统状态方程求解离散时间系统状态方程求解应用matlab函数G,H=c2d(A,B,T)实现从连续时间到离散时间的模型转换第30页/共41页2-52-5离散时间系统状态方程求解离散时间系统
7、状态方程求解线性定常系统用迭代法解矩阵差分方程第31页/共41页2-52-5离散时间系统状态方程求解离散时间系统状态方程求解K时刻的状态只与此时刻以前的输入采样值有关,与该时刻的输入采样值无关第32页/共41页2-6 2-6 应用Matlab的系统运动分析方法1:matlab函数来实现(1)expm(A*t)计算给定时刻t的状态转移阵;函数调用格式:fait=expm(A*t);(2)step(sys)求取阶跃输入时系统状态响应;函数调用格式:y,t,x=step(sys)(3)lsim(sys)求取任意输入时系统状态响应 函数调用格式:y,t,x=lsim(sys,u,t)(4)initia
8、l(sys)求取零输入时系统状态响应 函数调用格式:y,t,x=initial(sys,x0)第33页/共41页2-6 2-6 应用Matlab的系统运动分析 matlab的m文件文本如下:A=0-2;1-3;B=2;0;fait02=expm(A*0.2)fait04=expm(A*0.4)求下列系统在t=0.2,0.4秒时的状态转移阵第34页/共41页2-6 2-6 应用Matlab的系统运动分析 matlab的m文件文本如下:A=-20,19,-20;19,-21,20;40-40-40;B=0 1 2;C=1 0 2;D=0;sys=ss(A,B,C,D);y,t,x=step(sys
9、);plot(t,x),gridfgure(2);plot(t,y),grid求下列系统在u=1(t)时的状态响应和输出响应第35页/共41页2-6 2-6 应用Matlab的系统运动分析 matlab的m文件文本如下:A=-20,19,-20;19,-21,20;40-40-40;B=0 1 2;C=1 0 2;D=0;sys=ss(A,B,C,D);t=0:0.02:4;u=1+exp(-t).*cos(5*t);y,t,x=lsim(sys,u,t);plot(t,x),gridfigure(2);plot(t,y),grid求上例系统在 时的状态响应和输出响应x0=1 0 0;y,t,
10、x=initial(sys,x0);plot(t,x)第36页/共41页2-6 2-6 应用Matlab的系统运动分析方法2:应用matlab的simulink中的模型库建立mdl文件求下列系统的状态响应第37页/共41页2-6 2-6 应用Matlab的系统运动分析方法2:应用matlab的simulink中的模型库建立mdl文件求下列非线性系统的状态响应第38页/共41页第二章第二章 作业(习题)作业(习题)2-3,2-4(2),2-52-6,第39页/共41页大作业大作业1 1(1).推导倒立摆系统的状态空间模型(非线性)。(2).分别求在摆角为零(倒立摆)和1800(正摆)时的线性化状态空间模型。(3).采用matlab 求(1)非线性模型在摆角为00附近和1800附近时状态响应曲线。(4).采用matlab 分别求(2)模型在摆角为00附近时状态响应曲线。提示:仿真时间可以选择很小。(5).(3)和(4)的计算结果有什么关系?目的:强化建模和线性化知识,并能借助计算平目的:强化建模和线性化知识,并能借助计算平台完成实际系统的模型分析。台完成实际系统的模型分析。第40页/共41页感谢您的观看!第41页/共41页