控制系统仿真与CAD：第3章 线性控制系统的数学模型

1、2022/7/311控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院第3章 线性控制系统的数学模型2022/7/312控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院系统的数学模型系统数学模型的重要性系统仿真分析必须已知数学模型系统设计必须已知数学模型本课程数学模型是基础系统数学模型的获取建模方法：从已知的物理规律出发，用数学推导的方式建立起系统的数学模型辨识方法：由实验数据拟合系统的数学模型2022/7/313控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院系统数学模型的分类系统模型非线性线性连续离散混合单变量多变量定常时变2022/7/314控制系

2、统计算机辅助设计-MATLAB语言与应用东北大学信息学院主要内容线性连续系统的数学模型与MATLAB表示线性离散时间系统的数学模型方框图描述系统的化简系统模型的相互转换线性系统的模型降阶（自学为主）线性系统的模型辨识（自学为主）飞机运动的动力学建模（补充内容）本章要点简介2022/7/315控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.1 连续线性系统的数学 模型与MATLAB表示3.1.1线性系统的状态方程模型3.1.2 线性系统的传递函数模型3.1.3 线性系统的零极点模型3.1.4 多变量系统的传递函数矩阵模型2022/7/316控制系统计算机辅助设计-MATLAB语

3、言与应用东北大学信息学院3.1.1 线性连续系统数学模型及MATLAB 表示线性系统的传递函数模型 为阶次， 为常数， 物理可实现2022/7/317控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院传递函数的引入 Pierre-Simon Laplace (1749-1827)，法国数学家 Laplace变换 Laplace变换的一条重要性质： 若 则2022/7/318控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院传递函数表示数学方式MATLAB输入语句2022/7/319控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院传递函数输入举例例3-1 输

4、入传递函数模型MATLAB输入语句 在MATLAB环境中建立一个变量 G2022/7/3110控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院另外一种传递函数输入方法例3-2 如何处理如下的传递函数？定义算子 ，再输入传递函数2022/7/3111控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院应该根据给出传递函数形式选择输入方法例3-3 输入混合运算的传递函数模型 显然用第一种方法麻烦，所以2022/7/3112控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院MATLAB的传递函数对象2022/7/3113控制系统计算机辅助设计-MATLAB语言与应

5、用东北大学信息学院传递函数属性修改例3-4 延迟传递函数 ，即若假设复域变量为 ，则2022/7/3114控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院传递函数参数提取由于使用单元数组，直接用 不行有两种方法可以提取参数这样定义的优点：可以直接描述多变量系统第 i 输出对第 j 输入的传递函数2022/7/3115控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.1.2 线性系统的状态方程模型状态方程模型状态变量 ， 阶次 n ，输入和输出非线性函数： 一般非线性系统的状态方程描述2022/7/3116控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学

6、院线性状态方程时变模型线性时不变模型 (linear time invariant, LTI)2022/7/3117控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院线性时不变模型的MATLAB描述MATLAB 输入方法 矩阵是 方阵， 为 矩阵 为 矩阵， 为 矩阵可以直接处理多变量模型给出 矩阵即可注意维数的兼容性2022/7/3118控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-52022/7/3119控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院带时间延迟的状态方程数学模型MATLAB输入语句其他延迟属性：ioDelay2022/7/3

7、120控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.1.3 线性系统的零极点模型零极点模型是因式型传递函数模型零点 、极点 和增益零极点模型的 MATLAB表示2022/7/3121控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-5 零极点模型MATLAB输入方法另一种输入方法2022/7/3122控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.1.4 多变量系统传递函数矩阵模型传递函数矩阵 为第 i 输出对第 j 输入的传递函数可以先定义子传递函数，再由矩阵定义2022/7/3123控制系统计算机辅助设计-MATLAB语言与应用东北大

8、学信息学院例3-7 多变量模型2022/7/3124控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.2 线性离散时间系统的数学模型单变量系统：差分方程取代微分方程主要内容离散传递函数离散状态方程2022/7/3125控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.2.1 离散传递函数模型数学表示 （Z变换代替Laplace变换）MATLAB表示（采样周期 ）算子输入方法：2022/7/3126控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-8 离散传递函数，采样周期MATLAB输入方法另一种输入方法2022/7/3127控制系统计算机辅助

9、设计-MATLAB语言与应用东北大学信息学院离散延迟系统与输入数学模型延迟为采样周期的整数倍MATLAB输入方法2022/7/3128控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院滤波器型描述方法滤波器型离散模型分子、分母除以记 ，则2022/7/3129控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院MATLAB表示方法例3-92022/7/3130控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.2.2 离散状态方程模型数学形式注意兼容性MATLAB表示方法2022/7/3131控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院离散

10、延迟系统的状态方程数学模型MATLAB表示方法2022/7/3132控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.3 方框图描述系统的化简单环节模型前面已经介绍了实际系统为多个环节互连如何解决互连问题，获得等效模型？主要内容控制系统的典型连接结构节点移动时的等效变换复杂系统模型的简化2022/7/3133控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.3.1 控制系统的典型连接结构系统串、并联串联传递函数 并联传递函数2022/7/3134控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院串、并联状态方程模型串联系统的状态方程并联系统的状态方程

11、2022/7/3135控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院串、并联系统的MATLAB求解若一个模型为传递函数、另一个为状态方程，如何处理？将二者变换成同样结构再计算基于MATLAB的计算方法串联 注意次序：多变量系统并联优点，无需实现转换2022/7/3136控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院系统的反馈连接反馈连接正反馈负反馈2022/7/3137控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院状态方程的反馈等效方法其中若2022/7/3138控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院反馈连接的MATLAB

12、求解LTI 模型符号运算 (置于sym目录)2022/7/3139控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-10 2022/7/3140例3-11控制器为对角矩阵2022/7/3141控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.3.2 节点移动时的等效变换考虑模型难点：A点在回路间，移至输出端2022/7/3142控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院节点移动2022/7/3143控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.3.3 复杂系统模型的简化例3-12 原系统可以移动新支路模型2022/7/3

13、144控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院得出2022/7/3145控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-13 电机拖动模型 2022/7/3146控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院 信号单独输入得出另一个传递函数2022/7/3147控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院最终得出传递函数矩阵2022/7/3148控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4 系统模型的相互转换前面介绍的各种模型之间的相互等效变换主要内容连续模型和离散模型的相互转换系统传递函数的获

14、取控制系统的状态方程实现状态方程的最小实现传递函数与符号表达式的相互转换2022/7/3149控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4.1 连续模型和离散模型的相互转换连续状态方程的解析阶采样周期选择2022/7/3150控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院这样可以得出离散模型记则可以得出离散状态方程模型MATLAB函数直接求解2022/7/3151控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院还可以采用Tustin变换（双线性变换）例3-14 双输入模型，2022/7/3152控制系统计算机辅助设计-MATLAB语言与应用

15、东北大学信息学院输入模型、变换模型2022/7/3153控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-15 时间延迟系统的离散化MATLAB求解零阶保持器变换变换结果2022/7/3154控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院Tustin变换数学表示其他转换方法FOH 一阶保持器matched 单变量系统零极点不变imp 脉冲响应不变准则2022/7/3155控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院离散模型连续化对前面的变换求逆Tustin反变换MATLAB求解 (无需 )2022/7/3156控制系统计算机辅助设计-MATL

16、AB语言与应用东北大学信息学院例3-16 对前面的连续状态方程模型离散化，对结果再连续化，则 可以基本上还原连续模型2022/7/3157控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4.2 系统传递函数的获取已知状态方程两端Laplace变换则2022/7/3158控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院因此可以得出传递函数难点基于Fadeev-Fadeeva算法能得出更好结果由零极点模型，直接展开分子分母用MATLAB统一求解2022/7/3159控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-17 多变量模型，求传递函数矩阵2

17、022/7/3160控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4.3 控制系统的状态方程实现由传递函数到状态方程的转换不同状态变量选择，结果不唯一默认变换方式，采用MATLAB函数G可以是传递函数、状态方程和零极点模型适用于有延迟的、离散的或多变量模型2022/7/3161控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-18 连续多变量模型状态方程获取2022/7/3162控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院得出的状态方程模型ioDelay矩阵2022/7/3163控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学

18、院该模型可以转换回传递函数矩阵得出的转换结果2022/7/3164控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院均衡实现 （banlanced realization）由一般状态方程输入输出关系显著程度不明显，需要进一步变换均衡实现是一种很有用的方式用MATLAB直接求解得出均衡实现的模型得出排序的 Gram 矩阵2022/7/3165控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例原系统模型引入 (内部坐标变换)2022/7/3166控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4.4 状态方程的最小实现例3-19 观察传递函数模型未见有何

19、特殊求取零极点模型2022/7/3167控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院得出结果相同位置的零极点，可以对消问题：状态方程如何处理？MATLAB解决方法2022/7/3168控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-20 多变量模型不能直接看出是否最小实现2022/7/3169控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院MATLAB求解2022/7/3170控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.4.5 传递函数与符号表达式的相互转换传递函数到符号表达式表达式到传递函数置于sym目录下2022/7

20、/3171控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.5 线性系统模型降阶用低阶模型近似高阶模型和最小实现不同最早由Edward J. Davison提出(1966)主要内容 与Routh算法时间延迟模型的 近似带有延迟的最优降阶算法状态空间的降阶算法2022/7/3172控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.5.1 降阶算法 与 Routh 降阶算法原始模型寻求降阶模型假设2022/7/3173控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院展开原模型其中时间矩量 可以递推求出若已知状态方程模型2022/7/3174控制系统计算

21、机辅助设计-MATLAB语言与应用东北大学信息学院时间矩量的MATLAB求解 降阶思想：保留前 时间矩量2022/7/3175控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院对比系数，则2022/7/3176控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院这样可以得出2022/7/3177控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院 降阶求解函数2022/7/3178控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-21 原始模型Pad 近似结果2022/7/3179控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例

22、3-22 反例零极点模型求取稳定模型2022/7/3180控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院Pad 近似不稳定降阶模型Pad 不能保证降阶模型的稳定性不稳定降阶模型可能得出稳定降阶模型2022/7/3181控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院Routh 降阶方法与实例Routh算法（较烦琐，从略）2022/7/3182控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院Routh算法的最大特色：稳定系统降阶后能保证降阶模型稳定性例3-23 仍考虑稳定模型2022/7/3183控制系统计算机辅助设计-MATLAB语言与应用东北大学信

23、息学院3.5.3 时间延迟模型的 Pad 近似纯延迟的Pad近似方法近似函数纯滞后逼近2022/7/3184控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院编写 MATLAB 函数其中 r/m 任意选择可以选择 0/m ，以避免非最小相位模型2022/7/3185控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-24 纯延迟模型MATLAB求解拟合结果2022/7/3186控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-25 已知带有延迟的线性模型可以得出近似模型2022/7/3187控制系统计算机辅助设计-MATLAB语言与应用东北大学信息

24、学院3.5.4 带有时间延迟系统的 次最优降阶算法降阶模型的降阶效果误差定义ISE准则2022/7/3188控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院原模型降阶模型降阶误差定义2022/7/3189控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院参数向量误差MATLAB实现（从略）调用格式2022/7/3190控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-26 对给出的传递函数进行降阶研究可以给出下面的语句得出的降阶模型为2022/7/3191控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-27 已知高阶模型可以给出如

25、下命令得出的降阶模型2022/7/3192控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.5.4 状态方程模型的降阶算法均衡实现模型的降阶算法2022/7/3193控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院MATLAB求解函数例3-28 2022/7/3194控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院基于 Schur 均衡实现模型的降阶算法MATLAB求解函数例3-29 高阶传递函数思路：先转换成状态方程，再降阶2022/7/3195控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院模型输入与降阶Schur 降阶模型2022

26、/7/3196控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院最优 Hankel 范数的降阶模型近似MATLAB求解函数例3-30 仍采用前面模型2022/7/3197控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院降阶算法综述状态方程方法不能任意选择分母分子阶次，而很多传递函数方法可以降阶效果比较，下章给出时域响应比较频域响应比较降阶模型的应用仿真应用（用途越来越小）控制器设计应用2022/7/3198控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.6 线性系统的模型辨识模型辨识由已知实测数据获得系统模型的方法实测数据时域响应数据、频率响应数据主

27、要内容离散系统辨识方法辨识信号生成多变量系统辨识离散系统在线辨识2022/7/3199控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.6.1 离散系统的模型辨识离散传递函数模型对应的差分方程模型2022/7/31100控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院已知实测信号输入输出由数据可以得出2022/7/31101控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院矩阵形式定义残差最小指标最小二乘解2022/7/31102控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院系统辨识工具箱求解T 为结构体变量，T.a, T.b, tf(

28、T)当然由前面的公式也能直接求解2022/7/31103控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-31 实测数据2022/7/31104控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院基于MATLAB的求解2022/7/31105控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院数学形式辨识模型的提取还可以写成2022/7/31106控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院还可以由下面语句求解辨识结果2022/7/31107控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院直接辨识方法辨识结果辨识界面：ide

29、nt2022/7/31108控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.6.2 离散系统辨识信号的生成问题：什么样信号激励系统，辨识效果最好？有丰富频率信息的信号最好，如 PRBS伪随机二进制序列 pseudo-random binary sequence频率丰富值为可重复构建MATLAB直接生成2022/7/31109控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院例3-32 生成63个点的PRBS信号辨识效果残差明显减小2022/7/31110控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院连续系统的辨识可以考虑的方法 连续频率拟合方法，不

30、惟一 离散方法，再转换成连续模型例3-33 2022/7/31111控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院正弦信号激励辨识结果问题原因：输入频率单一2022/7/31112控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院3.6.3 多变量离散系统的辨识离散传递函数矩阵模型其中例3-342022/7/31113控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院MATLAB求解2022/7/31114控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院得出的高阶模型应该最小实现辨识结果2022/7/31115控制系统计算机辅助设计-MAT

31、LAB语言与应用东北大学信息学院3.6.4 离散系统的递推最小二乘辨识在系统运行过程中实时获取系统参数，而不是象前面介绍的方法那样一次性获得模型，适合于变参数模型的实时控制广泛应用于自适应控制渐近地逼近参数真值这里介绍算法，仿真研究将在后面介绍2022/7/31116控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院递推最小二乘辨识传递函数模型差分方程模型待辨识参数2022/7/31117控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院递推初值 和加权矩阵输入输出数据向量递推公式2022/7/31118控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院补充

32、内容 飞机的动力学模型坐标系地心惯性坐标系地心固定坐标系航迹坐标系机体坐标系2022/7/31119控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院坐标系之间的转换地心惯性坐标系到地心固定坐标系以上公式使用mathematics软件推导，matlab也具备这个功能2022/7/31120控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院坐标系之间的转换地心固定坐标系到航迹坐标系（L维度，经度）2022/7/31121控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院坐标系之间的转换航迹坐标系到机体坐标系（三个欧拉角）2022/7/31122控制系统计算机辅

33、助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程基本理论牛顿第二定理 动轴定理2022/7/31123控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程xyzP飞机的质心在地心惯性系中2022/7/31124控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程相对航迹坐标系的加速度随地球自转产生的向心加速度科氏加速度2022/7/31125控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程2022/7/31126控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心

34、动力学方程对于飞机2022/7/31127控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程忽略了地球自转，球形地球2022/7/31128控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机的质心动力学方程质心动力学方程与导航方程是不耦合的质心动力学方程质心运动方程（导航方程）2022/7/31129控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机绕质心动力学方程基本理论动量矩定理惯性系中很不方便体坐标系中注意：忽略了地球自转2022/7/31130控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院飞机绕质心运动学方程航迹坐标系到机体坐标系（绕三个欧拉角旋转）2022/7/31131控制系统计算机辅助设计-MATLAB