matlab-线性控制系统的数学模型

第3

章

线性控制系统的数学模型薛定宇著《控制系统计算机辅助设计---MATLAB语言与应用》第二版，清华大学出版社2006CAI课件开发：张望舒哈尔滨工程大学薛定宇东北大学2/4/20231控制系统计算机辅助设计---MATLAB语言与应用系统的数学模型系统数学模型的重要性系统仿真分析必须已知数学模型系统设计必须已知数学模型本课程数学模型是基础系统数学模型的获取建模方法：从已知的物理规律出发，用数学推导的方式建立起系统的数学模型辨识方法：由实验数据拟合系统的数学模型2/4/20232控制系统计算机辅助设计---MATLAB语言与应用系统数学模型的分类系统模型非线性线性连续离散混合单变量多变量定常时变2/4/20233控制系统计算机辅助设计---MATLAB语言与应用主要内容线性连续系统的数学模型与MATLAB表示线性离散时间系统的数学模型方框图描述系统的化简系统模型的相互转换线性系统的模型降阶线性系统的模型辨识本章要点简介2/4/20234控制系统计算机辅助设计---MATLAB语言与应用3.1连续线性系统的数学

模型与MATLAB表示3.1.1线性系统的状态方程模型3.1.2线性系统的传递函数模型3.1.3线性系统的零极点模型3.1.4多变量系统的传递函数矩阵模型2/4/20235控制系统计算机辅助设计---MATLAB语言与应用3.1.1线性连续系统数学模型及MATLAB表示线性系统的传递函数模型为阶次，为常数，物理可实现2/4/20236控制系统计算机辅助设计---MATLAB语言与应用传递函数的引入

Pierre-SimonLaplace

(1749--1827)，法国数学家

Laplace变换Laplace变换的一条重要性质：若则2/4/20237控制系统计算机辅助设计---MATLAB语言与应用传递函数表示数学方式MATLAB输入语句2/4/20238控制系统计算机辅助设计---MATLAB语言与应用传递函数输入举例例3-1输入传递函数模型MATLAB输入语句在MATLAB环境中建立一个变量G2/4/20239控制系统计算机辅助设计---MATLAB语言与应用另外一种传递函数输入方法例3-2如何处理如下的传递函数？定义算子，再输入传递函数2/4/202310控制系统计算机辅助设计---MATLAB语言与应用应该根据给出传递函数形式选择输入方法例3-3输入混合运算的传递函数模型显然用第一种方法麻烦，所以2/4/202311控制系统计算机辅助设计---MATLAB语言与应用MATLAB的传递函数对象2/4/202312控制系统计算机辅助设计---MATLAB语言与应用传递函数属性修改例3-4延迟传递函数，即若假设复域变量为，则2/4/202313控制系统计算机辅助设计---MATLAB语言与应用传递函数参数提取由于使用单元数组，直接用不行有两种方法可以提取参数这样定义的优点：可以直接描述多变量系统第i输入对第j输入的传递函数2/4/202314控制系统计算机辅助设计---MATLAB语言与应用3.1.2线性系统的状态方程模型状态方程模型状态变量，阶次n

，输入和输出非线性函数：一般非线性系统的状态方程描述2/4/202315控制系统计算机辅助设计---MATLAB语言与应用线性状态方程时变模型线性时不变模型(lineartimeinvariant,LTI)2/4/202316控制系统计算机辅助设计---MATLAB语言与应用线性时不变模型的MATLAB描述MATLAB输入方法矩阵是

方阵，为矩阵为矩阵，为矩阵可以直接处理多变量模型给出矩阵即可注意维数的兼容性2/4/202317控制系统计算机辅助设计---MATLAB语言与应用例3-52/4/202318控制系统计算机辅助设计---MATLAB语言与应用带时间延迟的状态方程数学模型MATLAB输入语句其他延迟属性：ioDelay2/4/202319控制系统计算机辅助设计---MATLAB语言与应用3.1.3线性系统的零极点模型零极点模型是因式型传递函数模型零点、极点和增益零极点模型的

MATLAB表示2/4/202320控制系统计算机辅助设计---MATLAB语言与应用例3-5零极点模型MATLAB输入方法另一种输入方法2/4/202321控制系统计算机辅助设计---MATLAB语言与应用3.1.4多变量系统传递函数矩阵模型传递函数矩阵为第i输出对第j输入的传递函数可以先定义子传递函数，再由矩阵定义2/4/202322控制系统计算机辅助设计---MATLAB语言与应用例3-7多变量模型2/4/202323控制系统计算机辅助设计---MATLAB语言与应用3.2线性离散时间系统的数学模型单变量系统：差分方程取代微分方程主要内容离散传递函数离散状态方程2/4/202324控制系统计算机辅助设计---MATLAB语言与应用3.2.1离散传递函数模型数学表示（Z变换代替Laplace变换）MATLAB表示（采样周期）算子输入方法：2/4/202325控制系统计算机辅助设计---MATLAB语言与应用例3-8离散传递函数，采样周期MATLAB输入方法另一种输入方法2/4/202326控制系统计算机辅助设计---MATLAB语言与应用离散延迟系统与输入数学模型延迟为采样周期的整数倍MATLAB输入方法2/4/202327控制系统计算机辅助设计---MATLAB语言与应用滤波器型描述方法滤波器型离散模型分子、分母除以记，则2/4/202328控制系统计算机辅助设计---MATLAB语言与应用MATLAB表示方法例3-92/4/202329控制系统计算机辅助设计---MATLAB语言与应用3.2.2离散状态方程模型数学形式注意兼容性MATLAB表示方法2/4/202330控制系统计算机辅助设计---MATLAB语言与应用离散延迟系统的状态方程数学模型MATLAB表示方法2/4/202331控制系统计算机辅助设计---MATLAB语言与应用3.3方框图描述系统的化简单环节模型前面已经介绍了实际系统为多个环节互连如何解决互连问题，获得等效模型？主要内容控制系统的典型连接结构节点移动时的等效变换复杂系统模型的简化2/4/202332控制系统计算机辅助设计---MATLAB语言与应用3.3.1控制系统的典型连接结构系统串、并联串联传递函数并联传递函数2/4/202333控制系统计算机辅助设计---MATLAB语言与应用串、并联状态方程模型串联系统的状态方程并联系统的状态方程2/4/202334控制系统计算机辅助设计---MATLAB语言与应用串、并联系统的MATLAB求解若一个模型为传递函数、另一个为状态方程，如何处理？将二者变换成同样结构再计算基于MATLAB的计算方法串联注意次序：多变量系统并联优点，无需实现转换2/4/202335控制系统计算机辅助设计---MATLAB语言与应用系统的反馈连接反馈连接正反馈负反馈2/4/202336控制系统计算机辅助设计---MATLAB语言与应用状态方程的反馈等效方法其中若2/4/202337控制系统计算机辅助设计---MATLAB语言与应用反馈连接的MATLAB求解LTI模型符号运算(置于@sym目录)2/4/202338控制系统计算机辅助设计---MATLAB语言与应用例3-102/4/202339控制系统计算机辅助设计---MATLAB语言与应用例3-11控制器为对角矩阵2/4/2023403.3.2节点移动时的等效变换考虑模型难点：A点在回路间，移至输出端2/4/202341控制系统计算机辅助设计---MATLAB语言与应用节点移动2/4/202342控制系统计算机辅助设计---MATLAB语言与应用3.3.3复杂系统模型的简化例3-12原系统可以移动新支路模型2/4/202343控制系统计算机辅助设计---MATLAB语言与应用得出2/4/202344控制系统计算机辅助设计---MATLAB语言与应用例3-13电机拖动模型

2/4/202345控制系统计算机辅助设计---MATLAB语言与应用信号单独输入得出另一个传递函数2/4/202346控制系统计算机辅助设计---MATLAB语言与应用最终得出传递函数矩阵2/4/202347控制系统计算机辅助设计---MATLAB语言与应用3.4系统模型的相互转换前面介绍的各种模型之间的相互等效变换主要内容连续模型和离散模型的相互转换系统传递函数的获取控制系统的状态方程实现状态方程的最小实现传递函数与符号表达式的相互转换2/4/202348控制系统计算机辅助设计---MATLAB语言与应用3.4.1连续模型和离散模型的相互转换连续状态方程的解析阶采样周期选择2/4/202349控制系统计算机辅助设计---MATLAB语言与应用这样可以得出离散模型记则可以得出离散状态方程模型MATLAB函数直接求解2/4/202350控制系统计算机辅助设计---MATLAB语言与应用还可以采用Tustin变换（双线性变换）例3-14双输入模型，2/4/202351控制系统计算机辅助设计---MATLAB语言与应用输入模型、变换模型2/4/202352控制系统计算机辅助设计---MATLAB语言与应用例3-15时间延迟系统的离散化MATLAB求解零阶保持器变换变换结果2/4/202353控制系统计算机辅助设计---MATLAB语言与应用Tustin变换数学表示其他转换方法FOH一阶保持器matched单变量系统零极点不变imp脉冲响应不变准则2/4/202354控制系统计算机辅助设计---MATLAB语言与应用离散模型连续化对前面的变换求逆Tustin反变换MATLAB求解(无需)2/4/202355控制系统计算机辅助设计---MATLAB语言与应用例3-16对前面的连续状态方程模型离散化，对结果再连续化，则

可以基本上还原连续模型2/4/202356控制系统计算机辅助设计---MATLAB语言与应用3.4.2系统传递函数的获取已知状态方程两端Laplace变换则2/4/202357控制系统计算机辅助设计---MATLAB语言与应用因此可以得出传递函数难点基于Fadeev-Fadeeva算法能得出更好结果由零极点模型，直接展开分子分母用MATLAB统一求解2/4/202358控制系统计算机辅助设计---MATLAB语言与应用例3-17多变量模型，求传递函数矩阵2/4/202359控制系统计算机辅助设计---MATLAB语言与应用3.4.3控制系统的状态方程实现由传递函数到状态方程的转换不同状态变量选择，结果不唯一默认变换方式，采用MATLAB函数G可以是传递函数、状态方程和零极点模型适用于有延迟的、离散的或多变量模型2/4/202360控制系统计算机辅助设计---MATLAB语言与应用例3-18连续多变量模型状态方程获取2/4/202361控制系统计算机辅助设计---MATLAB语言与应用得出的状态方程模型ioDelay矩阵2/4/202362控制系统计算机辅助设计---MATLAB语言与应用该模型可以转换回传递函数矩阵得出的转换结果2/4/202363控制系统计算机辅助设计---MATLAB语言与应用均衡实现（banlancedrealization）由一般状态方程输入输出关系显著程度不明显，需要进一步变换均衡实现是一种很有用的方式用MATLAB直接求解得出均衡实现的模型得出排序的Gram矩阵2/4/202364控制系统计算机辅助设计---MATLAB语言与应用例原系统模型引入(内部坐标变换)2/4/202365控制系统计算机辅助设计---MATLAB语言与应用3.4.4状态方程的最小实现例3-19观察传递函数模型未见有何特殊求取零极点模型2/4/202366控制系统计算机辅助设计---MATLAB语言与应用得出结果相同位置的零极点，可以对消问题：状态方程如何处理？MATLAB解决方法2/4/202367控制系统计算机辅助设计---MATLAB语言与应用例3-20多变量模型不能直接看出是否最小实现2/4/202368控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解2/4/202369控制系统计算机辅助设计---MATLAB语言与应用3.4.5传递函数与符号表达式

的相互转换传递函数到符号表达式表达式到传递函数置于@sym目录下2/4/202370控制系统计算机辅助设计---MATLAB语言与应用3.5线性系统模型降阶用低阶模型近似高阶模型和最小实现不同最早由EdwardJ.Davison提出(1966)主要内容与Routh算法时间延迟模型的近似带有延迟的最优降阶算法状态空间的降阶算法2/4/202371控制系统计算机辅助设计---MATLAB语言与应用3.5.1降阶算法

与Routh

降阶算法原始模型寻求降阶模型假设2/4/202372控制系统计算机辅助设计---MATLAB语言与应用展开原模型其中时间矩量可以递推求出若已知状态方程模型2/4/202373控制系统计算机辅助设计---MATLAB语言与应用时间矩量的MATLAB求解降阶思想：保留前时间矩量2/4/202374控制系统计算机辅助设计---MATLAB语言与应用对比系数，则2/4/202375控制系统计算机辅助设计---MATLAB语言与应用这样可以得出2/4/202376控制系统计算机辅助设计---MATLAB语言与应用降阶求解函数2/4/202377控制系统计算机辅助设计---MATLAB语言与应用例3-21原始模型Padé

近似结果2/4/202378控制系统计算机辅助设计---MATLAB语言与应用例3-22反例零极点模型求取稳定模型2/4/202379控制系统计算机辅助设计---MATLAB语言与应用Padé

近似不稳定降阶模型Padé

不能保证降阶模型的稳定性不稳定降阶模型可能得出稳定降阶模型2/4/202380控制系统计算机辅助设计---MATLAB语言与应用Routh

降阶方法与实例Routh算法（较烦琐，从略）2/4/202381控制系统计算机辅助设计---MATLAB语言与应用Routh算法的最大特色：稳定系统降阶后能保证降阶模型稳定性例3-23仍考虑稳定模型2/4/202382控制系统计算机辅助设计---MATLAB语言与应用3.5.3时间延迟模型的Padé

近似纯延迟的Padé近似方法近似函数纯滞后逼近2/4/202383控制系统计算机辅助设计---MATLAB语言与应用编写MATLAB函数其中r/m任意选择可以选择0/m，以避免非最小相位模型2/4/202384控制系统计算机辅助设计---MATLAB语言与应用例3-24纯延迟模型MATLAB求解拟合结果2/4/202385控制系统计算机辅助设计---MATLAB语言与应用例3-25已知带有延迟的线性模型可以得出近似模型2/4/202386控制系统计算机辅助设计---MATLAB语言与应用3.5.4带有时间延迟系统的

次最优降阶算法降阶模型的降阶效果误差定义ISE准则2/4/202387控制系统计算机辅助设计---MATLAB语言与应用原模型降阶模型降阶误差定义2/4/202388控制系统计算机辅助设计---MATLAB语言与应用参数向量误差MATLAB实现（从略）调用格式2/4/202389控制系统计算机辅助设计---MATLAB语言与应用例3-26对给出的传递函数进行降阶研究可以给出下面的语句得出的降阶模型为2/4/202390控制系统计算机辅助设计---MATLAB语言与应用例3-27已知高阶模型可以给出如下命令得出的降阶模型2/4/202391控制系统计算机辅助设计---MATLAB语言与应用3.5.4状态方程模型的降阶算法均衡实现模型的降阶算法2/4/202392控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解函数例3-282/4/202393控制系统计算机辅助设计---MATLAB语言与应用基于Schur

均衡实现模型的降阶算法MATLAB求解函数例3-29高阶传递函数思路：先转换成状态方程，再降阶2/4/202394控制系统计算机辅助设计---MATLAB语言与应用模型输入与降阶Schur

降阶模型2/4/202395控制系统计算机辅助设计---MATLAB语言与应用最优Hankel

范数的降阶模型近似MATLAB求解函数例3-30仍采用前面模型2/4/202396控制系统计算机辅助设计---MATLAB语言与应用降阶算法综述状态方程方法不能任意选择分母分子阶次，而很多传递函数方法可以降阶效果比较，下章给出时域响应比较频域响应比较降阶模型的应用仿真应用（用途越来越小）控制器设计应用2/4/202397控制系统计算机辅助设计---MATLAB语言与应用3.6线性系统的模型辨识模型辨识由已知实测数据获得系统模型的方法实测数据时域响应数据、频率响应数据主要内容离散系统辨识方法辨识信号生成多变量系统辨识离散系统在线辨识2/4/202398控制系统计算机辅助设计---MATLAB语言与应用3.6.1离散系统的模型辨识离散传递函数模型对应的差分方程模型2/4/202399控制系统计算机辅助设计---MATLAB语言与应用已知实测信号输入输出由数据可以得出2/4/2023100控制系统计算机辅助设计---MATLAB语言与应用矩阵形式定义残差最小指标最小二乘解2/4/2023101控制系统计算机辅助设计---MATLAB语言与应用系统辨识工具箱求解T为结构体变量，T.a,T.b,tf(T)当然由前面的公式也能直接求解2/4/2023102控制系统计算机辅助设计---MATLAB语言与应用例3-31实测数据2/4/2023103控制系统计算机辅助设计---MATLAB语言与应用基于MATLAB的求解2/4/2023104控制系统计算机辅助设计---MATLAB语言与应用数学形式辨识模型的提取还可以写成2/4/2023105控制系统计算机辅助设计---MATLAB语言与应用还可以由下面语句求解辨识结果2/4/2023106控制系统计算机辅助设计---MATLAB语言与应用直接辨识方法辨识结果辨识界面：ident2/4/2023107控制系统计算机辅助设计---MATLAB语言与应用3.6.2离散系统辨识信号的生成问题：什么样信号激励系统，辨识效果最好？有丰富频率信息的信号最好，如PRBS伪随机二进制序列pseudo-randombinarysequence频率丰富值为可重复构建MATLAB直接生成2/4/2023108控制系统计算机辅助设计---MATLAB语言与应用例3-32生成63个点的PRBS信号辨识效果残差明显减小2/4/2023109控制系统计算机辅助设计---MATLAB语言与应用连续系统的辨识可以考虑的方法连续频率拟合方法，不惟一离散方法，再转换成连续模型例3-332/4/2023110控制系统计算机辅助设计---MATLAB语言与应用正弦信号激励辨识结果问题原因：输入频率单一2/4/2023111控制系统计算机辅助设计---MATLAB语言与应用3.6.3多变量离散系统的辨识离散传递函数矩阵模型其中例3-342/4/2023112控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解2/4/2023113控制系统计算机辅助设计---MATLAB语言与应用得出的高阶模型应该最小实现辨识结果2/4/2023114控制系统计算机辅助设计---MATLAB语言与应用3.6.4离散系统的递推最小二乘辨识在系统运行过程中实时获取系统参数，而不是象前面介绍的方法那样一次性获得模型，适合于变参数模型的实时控制广泛应用于自适应控制渐近地逼近参数真值这里介绍算法，仿真研究将在后面介绍2/4/2023115控制系统计算机辅助设计---MATLAB语言与应用递推最小二乘辨识传递函数模型差分方程模型待辨识参数2/4/2023116控制系统计算机辅助设计---MATLAB语言与应用递推初值和加权矩阵输入输出数据向量递推公式2/4/2023117控制系统计算机辅助设计---MATLAB语