控制系统模型与转换

控制系统模型与转换2026/8/172系统的数学模型系统数学模型的重要性系统仿真分析必须已知数学模型系统设计必须已知数学模型本课程数学模型是基础系统数学模型的获取建模方法：从已知的物理规律出发，用数学推导的方式建立起系统的数学模型辨识方法：由实验数据拟合系统的数学模型2026/8/173系统数学模型的分类系统模型非线性线性连续离散混合单变量多变量定常时变2026/8/174主要内容线性连续系统的数学模型与MATLAB表示线性离散时间系统的数学模型方框图描述系统的化简系统模型的相互转换线性系统的模型降阶线性系统的模型辨识本章要点简介2026/8/1753.1连续线性系统的数学

模型与MATLAB表示3.1.1线性系统的状态方程模型3.1.2线性系统的传递函数模型3.1.3线性系统的零极点模型3.1.4多变量系统的传递函数矩阵模型2026/8/1763.1.1线性连续系统数学模型及MATLAB表示线性系统的传递函数模型为阶次，为常数，物理可实现大家有疑问的，可以询问和交流可以互相讨论下，但要小声点2026/8/178传递函数的引入Pierre-SimonLaplace(1749--1827)，法国数学家

Laplace变换Laplace变换的一条重要性质：若则2026/8/179传递函数表示数学方式MATLAB输入语句2026/8/1710传递函数输入举例例3-1输入传递函数模型MATLAB输入语句在MATLAB环境中建立一个变量G2026/8/1711当传递函数的分子或分母由若干个多项式乘积表示时，可由MATLAB提供的多项式乘法运算函数conv()来处理，以便获得分子和分母多项式向量，此函数的调用格式为

c=conv(a,b)其中a和b分别为由两个多项式系数构成的向量，而c为a和b多项式的乘积多项式系数向量。conv()函数的调用是允许多级嵌套的。2026/8/1712另外一种传递函数输入方法例3-2若给定系统的传递函数为解则可以将其用下列MATLAB语句表示>>num=4*conv([12],[166])>>den=conv([10],conv([11],conv([11],conv([11],[1325]))))2026/8/1713对具有r个输入和m个输出的多变量系统，可把m×r的传递函数阵G(s)写成和单变量系统传递函数相类似的形式，即

式中B0,B1,…,Bn

均为m×r

实常数矩阵，分母多项式为该传递函数阵的特征多项式。在MATLAB控制系统工具箱中，提供了表示单输入多输出系统的表示方法，即

num=[B0B1…Bn];den=[1a1

a2

…an]其中分子系数包含在矩阵num中，num行数与输出y

的维数一致，每行对应一个输出，den是行向量,为传递函数阵公分母多项式系数。2026/8/1714例3-3

对于单输入多输出系统

解则可将其用下列MATLAB语句表示>>num=[0032;1025];den=[3521];2026/8/1715MATLAB的传递函数对象2026/8/1716传递函数属性修改例3-4延迟传递函数，即2026/8/1717传递函数参数提取由于使用单元数组，直接用不行有两种方法可以提取参数这样定义的优点：可以直接描述多变量系统第i输入对第j输出的传递函数2026/8/17183.1.2线性系统的状态方程模型状态方程模型状态变量，阶次n

，输入和输出非线性函数：一般非线性系统的状态方程描述2026/8/1719线性状态方程时变模型线性时不变模型(lineartimeinvariant,LTI)2026/8/1720线性时不变模型的MATLAB描述MATLAB输入方法矩阵是方阵，为矩阵为矩阵，为矩阵可以直接处理多变量模型给出矩阵即可注意维数的兼容性2026/8/1721例3-52026/8/1722带时间延迟的状态方程数学模型MATLAB输入语句其他延迟属性：ioDelay2026/8/17233.1.3线性系统的零极点模型零极点模型是因式型传递函数模型零点、极点和增益零极点模型的

MATLAB表示2026/8/1724例3-5零极点模型MATLAB输入方法另一种输入方法2026/8/17253.1.4多变量系统传递函数矩阵模型传递函数矩阵为第i输出对第j输入的传递函数可以先定义子传递函数，再由矩阵定义2026/8/1726例3-7多变量模型2026/8/17273.2线性离散时间系统的数学模型单变量系统：差分方程取代微分方程主要内容离散传递函数离散状态方程2026/8/17283.2.1离散传递函数模型数学表示（Z变换代替Laplace变换）MATLAB表示（采样周期）算子输入方法：2026/8/1729例3-8离散传递函数，采样周期MATLAB输入方法另一种输入方法2026/8/1730离散延迟系统的输入数学模型延迟为采样周期的整数倍MATLAB输入方法2026/8/1731滤波器型描述方法滤波器型离散模型分子、分母除以记，则2026/8/1732MATLAB表示方法例3-92026/8/17333.2.2离散状态方程模型数学形式注意兼容性MATLAB表示方法2026/8/1734离散延迟系统的状态方程数学模型MATLAB表示方法2026/8/17353.3方框图描述系统的化简单环节模型前面已经介绍了，实际系统为多个环节互连，如何解决互连问题，获得等效模型？主要内容控制系统的典型连接结构节点移动时的等效变换复杂系统模型的简化2026/8/17363.3.1控制系统的典型连接结构系统串、并联串联传递函数并联传递函数2026/8/1737串、并联状态方程模型串联系统的状态方程并联系统的状态方程2026/8/1738串、并联系统的MATLAB求解若一个模型为传递函数、另一个为状态方程，如何处理？将二者变换成同样结构再计算基于MATLAB的计算方法串联注意次序：多变量系统并联优点，无需实现转换2026/8/1739系统的反馈连接反馈连接正反馈负反馈2026/8/1740状态方程的反馈等效方法其中若2026/8/1741反馈连接的MATLAB求解LTI模型符号运算(置于@sym目录)2026/8/1742例3-102026/8/1743例3-11控制器为对角矩阵2026/8/17442026/8/17453.3.2节点移动时的等效变换考虑模型难点：A点在回路间，移至输出端2026/8/1746节点移动2026/8/17473.3.3复杂系统模型的简化例3-12原系统可以移动新支路模型2026/8/1748得出2026/8/1749例3-13电机拖动模型

2026/8/1750

信号单独输入得出另一个传递函数2026/8/1751最终得出传递函数矩阵2026/8/17523.4系统模型的相互转换前面介绍的各种模型之间的相互等效变换主要内容连续模型和离散模型的相互转换系统传递函数的获取控制系统的状态方程实现状态方程的最小实现传递函数与符号表达式的相互转换2026/8/17533.4.1连续模型和离散模型的相互转换连续状态方程的解析解采样周期选择2026/8/1754根据上式可得系统离散化后的状态方程（设T-τ=t）或

式中令τ=τ-kT，则上式可得2026/8/1755输出变量的差分方程，可由输出方程直接确定，即有

则可以得出离散状态方程模型MATLAB函数直接求解2026/8/1756还可以采用Tustin变换（双线性变换）例3-14双输入模型，2026/8/1757输入模型、变换2026/8/1758模型2026/8/1759例3-15时间延迟系统的离散化MATLAB求解零阶保持器变换变换结果2026/8/1760Tustin变换数学表示其他转换方法FOH一阶保持器matched单变量系统零极点不变imp脉冲响应不变准则2026/8/1761离散模型连续化对前面的变换求逆Tustin反变换MATLAB求解(无需)2026/8/1762例3-16对前面的连续状态方程模型离散化，对结果再连续化，则

可以基本上还原连续模型2026/8/17633.4.2系统传递函数的获取已知状态方程两端Laplace变换则2026/8/1764因此可以得出传递函数难点基于Fadeev-Fadeeva算法能得出更好结果由零极点模型，直接展开分子分母用MATLAB统一求解2026/8/1765例3-17多变量模型，求传递函数矩阵2026/8/17663.4.3控制系统的状态方程实现由传递函数到状态方程的转换不同状态变量选择，结果不唯一默认变换方式，采用MATLAB函数G可以是传递函数、状态方程和零极点模型适用于有延迟的、离散的或多变量模型2026/8/1767例3-18连续多变量模型状态方程获取2026/8/1768得出的状态方程模型ioDelay矩阵2026/8/1769该模型可以转换回传递函数矩阵得出的转换结果2026/8/17703.4.4状态方程的最小实现例3-19观察传递函数模型未见有何特殊求取零极点模型2026/8/1771得出结果相同位置的零极点，可以对消问题：状态方程如何处理？MATLAB解决方法2026/8/1772例3-20多变量模型不能直接看出是否最小实现2026/8/1773MATLAB求解2026/8/17743.5线性系统模型降阶用低阶模型近似高阶模型和最小实现不同最早由EdwardJ.Davison提出(1966)主要内容与Routh算法时间延迟模型的近似带有延迟的最优降阶算法状态空间的降阶算法2026/8/17753.5.1降阶算法

与Routh降阶算法原始模型寻求降阶模型假设2026/8/1776展开原模型其中时间矩量可以递推求出若已知状态方程模型2026/8/1777时间矩量的MATLAB求解降阶思想：保留前时间矩量2026/8/1778对比系数，则2026/8/1779这样可以得出2026/8/1780

降阶求解函数2026/8/1781例3-21原始模型Padé

近似结果2026/8/1782例3-22反例零极点模型求取稳定模型2026/8/1783Padé

近似不稳定降阶模型Padé

不能保证降阶模型的稳定性不稳定降阶模型可能得出稳定降阶模型2026/8/1784Routh降阶方法与实例Routh算法（较烦琐，从略）2026/8/1785Routh算法的最大特色：稳定系统降阶后能保证降阶模型稳定性例3-23仍考虑稳定模型2026/8/17863.5.2时间延迟模型的Padé

近似纯延迟的Padé近似方法近似函数纯滞后逼近2026/8/1787编写MATLAB函数其中r/k任意选择2026/8/1788例3-24纯延迟模型MATLAB求解拟合结果2026/8/17893.5.3带有时间延迟系统的

次最优降阶算法降阶模型的降阶效果原模型降阶模型降阶误差定义2026/8/1790降阶模型的降阶效果误差定义ISE准则2026/8/1791参数向量误差MATLAB实现（从略）调用格式2026/8/1792例3-26对给出的传递函数进行降阶研究可以给出下面的语句得出的降阶模型为2026/8/1793例3-27已知高阶模型可以给出如下命令得出的降阶模型2026/8/1794降阶算法综述降阶效果比较，下章给出时域响应比较频域响应比较降阶模型的应用仿真应用（用途越来越小）控制器设计应用2026/8/17953.6线性系统的模型辨识模型辨识由已知实测数据获得系统模型的方法实测数据时域响应数据、频率响应数据主要内容离散系统辨识方法辨识信号生成多变量系统辨识离散系统在线辨识2026/8/17963.6.1离散系统的模型辨识离散传递函数模型对应的差分方程模型2026/8/1797已知实测信号输入输出由数据可以得出2026/8/1798矩阵形式定义残差最小指标最小二乘解2026/8/1799系统辨识工具箱求解T为结构体变量，T.a,T.b,tf(T)当然由前面的公式也能直接求解2026/8/17100例3-31实测数据2026/8/17101基于MATLAB的求解2026/8/17102数学形式辨识模型的提取还可以写成2026/8/17103还可以由下面语句求解辨识结果2026/8/17104直接辨识方法辨识结果辨识界面：ident2026/8/171053.6.2离散系统辨识信号的生成问题：什么样信号激励系统，辨识效果最好？有丰富频率信息的信号最好，如PRBS伪随机二进制序列pseudo-randombinarysequence频率丰富