控制系统仿真与CAD-chapter-3

2020/6/6,1,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,第3章线性控制系统的数学模型,薛定宇著控制系统计算机辅助设计-MATLAB语言与应用第二版，清华大学出版社2006CAI课件开发：张望舒哈尔滨工程大学薛定宇东北大学,2020/6/6,2,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统的数学模型,系统数学模型的重要性系统仿真分析必须已知数学模型系统设计必须已知数学模型本课程数学模型是基础系统数学模型的获取建模方法：从已知的物理规律出发，用数学推导的方式建立起系统的数学模型辨识方法：由实验数据拟合系统的数学模型,2020/6/6,3,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统数学模型的分类,系统模型,非线性,线性,连续,离散,混合,单变量,多变量,定常,时变,2020/6/6,4,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,主要内容,线性连续系统的数学模型与MATLAB表示线性离散时间系统的数学模型方框图描述系统的化简系统模型的相互转换线性系统的模型降阶线性系统的模型辨识本章要点简介,2020/6/6,5,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.1连续线性系统的数学模型与MATLAB表示,3.1.1线性系统的状态方程模型3.1.2线性系统的传递函数模型3.1.3线性系统的零极点模型3.1.4多变量系统的传递函数矩阵模型,2020/6/6,6,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.1.1线性连续系统数学模型及MATLAB表示,线性系统的传递函数模型为阶次，为常数，物理可实现,2020/6/6,7,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,传递函数的引入,Pierre-SimonLaplace(1749-1827)，法国数学家Laplace变换Laplace变换的一条重要性质：若则,2020/6/6,8,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,传递函数表示,数学方式MATLAB输入语句,2020/6/6,9,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,传递函数输入举例,例3-1输入传递函数模型MATLAB输入语句在MATLAB环境中建立一个变量G,2020/6/6,10,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,另外一种传递函数输入方法,例3-2如何处理如下的传递函数？定义算子，再输入传递函数,2020/6/6,11,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,应该根据给出传递函数形式选择输入方法例3-3输入混合运算的传递函数模型显然用第一种方法麻烦，所以,2020/6/6,12,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB的传递函数对象,2020/6/6,13,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,传递函数属性修改,例3-4延迟传递函数，即若假设复域变量为，则,2020/6/6,14,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,传递函数参数提取,由于使用单元数组，直接用不行有两种方法可以提取参数这样定义的优点：可以直接描述多变量系统第i输入对第j输入的传递函数,2020/6/6,15,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.1.2线性系统的状态方程模型,状态方程模型状态变量，阶次n，输入和输出非线性函数：一般非线性系统的状态方程描述,2020/6/6,16,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,线性状态方程,时变模型线性时不变模型(lineartimeinvariant,LTI),2020/6/6,17,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,线性时不变模型的MATLAB描述,MATLAB输入方法矩阵是方阵，为矩阵为矩阵，为矩阵可以直接处理多变量模型给出矩阵即可注意维数的兼容性,2020/6/6,18,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-5,2020/6/6,19,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,带时间延迟的状态方程,数学模型MATLAB输入语句其他延迟属性：ioDelay,2020/6/6,20,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.1.3线性系统的零极点模型,零极点模型是因式型传递函数模型零点、极点和增益零极点模型的MATLAB表示,2020/6/6,21,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-5零极点模型MATLAB输入方法另一种输入方法,2020/6/6,22,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.1.4多变量系统传递函数矩阵模型,传递函数矩阵为第i输出对第j输入的传递函数可以先定义子传递函数，再由矩阵定义,2020/6/6,23,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-7多变量模型,2020/6/6,24,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.2线性离散时间系统的数学模型,单变量系统：差分方程取代微分方程主要内容离散传递函数离散状态方程,2020/6/6,25,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.2.1离散传递函数模型,数学表示（Z变换代替Laplace变换）MATLAB表示（采样周期）算子输入方法：,2020/6/6,26,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-8离散传递函数，采样周期MATLAB输入方法另一种输入方法,2020/6/6,27,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散延迟系统与输入,数学模型延迟为采样周期的整数倍MATLAB输入方法,2020/6/6,28,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,滤波器型描述方法,滤波器型离散模型分子、分母除以记，则,2020/6/6,29,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB表示方法例3-9,2020/6/6,30,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.2.2离散状态方程模型,数学形式注意兼容性MATLAB表示方法,2020/6/6,31,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散延迟系统的状态方程,数学模型MATLAB表示方法,2020/6/6,32,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.3方框图描述系统的化简,单环节模型前面已经介绍了实际系统为多个环节互连，如何解决互连问题，获得等效模型？主要内容控制系统的典型连接结构节点移动时的等效变换复杂系统模型的简化,2020/6/6,33,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.3.1控制系统的典型连接结构,系统串、并联串联传递函数并联传递函数,2020/6/6,34,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,串、并联状态方程模型,串联系统的状态方程并联系统的状态方程,2020/6/6,35,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,串、并联系统的MATLAB求解,若一个模型为传递函数、另一个为状态方程，如何处理？将二者变换成同样结构再计算基于MATLAB的计算方法串联注意次序：多变量系统并联优点，无需实现转换,2020/6/6,36,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统的反馈连接,反馈连接正反馈负反馈,2020/6/6,37,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,状态方程的反馈等效方法,其中若,2020/6/6,38,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,反馈连接的MATLAB求解,LTI模型符号运算(置于sym目录),2020/6/6,39,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-10,2020/6/6,40,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-11控制器为对角矩阵,2020/6/6,41,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,2020/6/6,42,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.3.2节点移动时的等效变换,考虑模型难点：A点在回路间，移至输出端,2020/6/6,43,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,节点移动,2020/6/6,44,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.3.3复杂系统模型的简化,例3-12原系统可以移动新支路模型,2020/6/6,45,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出,2020/6/6,46,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-13电机拖动模型,2020/6/6,47,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,信号单独输入得出另一个传递函数,2020/6/6,48,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,最终得出传递函数矩阵,2020/6/6,49,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4系统模型的相互转换,前面介绍的各种模型之间的相互等效变换主要内容连续模型和离散模型的相互转换系统传递函数的获取控制系统的状态方程实现状态方程的最小实现传递函数与符号表达式的相互转换,2020/6/6,50,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4.1连续模型和离散模型的相互转换,连续状态方程的解析阶采样周期选择,2020/6/6,51,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,这样可以得出离散模型记则可以得出离散状态方程模型MATLAB函数直接求解,2020/6/6,52,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,还可以采用Tustin变换（双线性变换）例3-14双输入模型，,2020/6/6,53,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,输入模型、变换,2020/6/6,54,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,模型,2020/6/6,55,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-15时间延迟系统的离散化MATLAB求解零阶保持器变换变换结果,2020/6/6,56,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Tustin变换数学表示其他转换方法FOH一阶保持器matched单变量系统零极点不变imp脉冲响应不变准则,2020/6/6,57,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散模型连续化,对前面的变换求逆Tustin反变换MATLAB求解(无需),2020/6/6,58,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-16对前面的连续状态方程模型离散化，对结果再连续化，则可以基本上还原连续模型,2020/6/6,59,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4.2系统传递函数的获取,已知状态方程两端Laplace变换则,2020/6/6,60,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,因此可以得出传递函数难点基于Fadeev-Fadeeva算法能得出更好结果由零极点模型，直接展开分子分母用MATLAB统一求解,2020/6/6,61,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-17多变量模型，求传递函数矩阵,2020/6/6,62,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4.3控制系统的状态方程实现,由传递函数到状态方程的转换不同状态变量选择，结果不唯一默认变换方式，采用MATLAB函数G可以是传递函数、状态方程和零极点模型适用于有延迟的、离散的或多变量模型,2020/6/6,63,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-18连续多变量模型状态方程获取,2020/6/6,64,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出的状态方程模型ioDelay矩阵,2020/6/6,65,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,该模型可以转换回传递函数矩阵得出的转换结果,2020/6/6,66,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,均衡实现（banlancedrealization）,由一般状态方程输入输出关系显著程度不明显，需要进一步变换均衡实现是一种很有用的方式用MATLAB直接求解得出均衡实现的模型得出排序的Gram矩阵,2020/6/6,67,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4.4状态方程的最小实现,例3-19观察传递函数模型未见有何特殊求取零极点模型,2020/6/6,68,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出结果相同位置的零极点，可以对消问题：状态方程如何处理？MATLAB解决方法,2020/6/6,69,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-20多变量模型不能直接看出是否最小实现,2020/6/6,70,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解,2020/6/6,71,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.4.5传递函数与符号表达式的相互转换,传递函数到符号表达式表达式到传递函数置于sym目录下,2020/6/6,72,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.5线性系统模型降阶,用低阶模型近似高阶模型和最小实现不同最早由EdwardJ.Davison提出(1966)主要内容与Routh算法时间延迟模型的近似带有延迟的最优降阶算法状态空间的降阶算法,2020/6/6,73,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.5.1降阶算法与Routh降阶算法,原始模型寻求降阶模型假设,2020/6/6,74,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,展开原模型其中时间矩量可以递推求出若已知状态方程模型,2020/6/6,75,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,时间矩量的MATLAB求解降阶思想：保留前时间矩量,2020/6/6,76,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,对比系数，则,2020/6/6,77,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,这样可以得出,2020/6/6,78,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,降阶求解函数,2020/6/6,79,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-21原始模型Pad近似结果,2020/6/6,80,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-22反例零极点模型求取稳定模型,2020/6/6,81,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Pad近似不稳定降阶模型Pad不能保证降阶模型的稳定性不稳定降阶模型可能得出稳定降阶模型,2020/6/6,82,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Routh降阶方法与实例,Routh算法（较烦琐，从略）,2020/6/6,83,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Routh算法的最大特色：稳定系统降阶后能保证降阶模型稳定性例3-23仍考虑稳定模型,2020/6/6,84,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.5.3时间延迟模型的Pad近似,纯延迟的Pad近似方法近似函数纯滞后逼近,2020/6/6,85,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,编写MATLAB函数其中r/m任意选择可以选择0/m，以避免非最小相位模型,2020/6/6,86,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-24纯延迟模型MATLAB求解拟合结果,2020/6/6,87,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-25已知带有延迟的线性模型可以得出近似模型,2020/6/6,88,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.5.4带有时间延迟系统的次最优降阶算法,降阶模型的降阶效果误差定义ISE准则,2020/6/6,89,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,原模型降阶模型降阶误差定义,2020/6/6,90,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,参数向量误差MATLAB实现（从略）调用格式,2020/6/6,91,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-26对给出的传递函数进行降阶研究可以给出下面的语句得出的降阶模型为,2020/6/6,92,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-27已知高阶模型可以给出如下命令得出的降阶模型,2020/6/6,93,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.5.4状态方程模型的降阶算法,均衡实现模型的降阶算法,2020/6/6,94,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解函数例3-28,2020/6/6,95,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,基于Schur均衡实现模型的降阶算法,MATLAB求解函数例3-29高阶传递函数思路：先转换成状态方程，再降阶,2020/6/6,96,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,模型输入与降阶Schur降阶模型,2020/6/6,97,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,最优Hankel范数的降阶模型近似,MATLAB求解函数例3-30仍采用前面模型,2020/6/6,98,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,降阶算法综述,状态方程方法不能任意选择分母分子阶次，而很多传递函数方法可以降阶效果比较，下章给出时域响应比较频域响应比较降阶模型的应用仿真应用（用途越来越小）控制器设计应用,2020/6/6,99,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.6线性系统的模型辨识,模型辨识由已知实测数据获得系统模型的方法实测数据时域响应数据、频率响应数据主要内容离散系统辨识方法辨识信号生成多变量系统辨识离散系统在线辨识,2020/6/6,100,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.6.1离散系统的模型辨识,离散传递函数模型对应的差分方程模型,2020/6/6,101,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,已知实测信号输入输出由数据可以得出,2020/6/6,102,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,矩阵形式定义残差最小指标最小二乘解,2020/6/6,103,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统辨识工具箱求解T为结构体变量，T.a,T.b,tf(T)当然由前面的公式也能直接求解,2020/6/6,104,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-31实测数据,2020/6/6,105,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,基于MATLAB的求解,2020/6/6,106,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,数学形式辨识模型的提取还可以写成,2020/6/6,107,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,还可以由下面语句求解辨识结果,2020/6/6,108,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,直接辨识方法辨识结果辨识界面：ident,2020/6/6,109,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.6.2离散系统辨识信号的生成,问题：什么样信号激励系统，辨识效果最好？有丰富频率信息的信号最好，如PRBS伪随机二进制序列pseudo-randombinarysequence频率丰富值为可重复构建MATLAB直接生成,2020/6/6,110,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例3-32生成63个点的PRBS信号辨识效果残差明显减小,2020/6/6,111,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,连续系统的辨识,可以考虑的方法连续频率拟合方法，不惟一离散方法，再转换成连续模型例3-33,2020/6/6,112,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,正弦信号激励辨识结果问题原因：输入频率单一,2020/6/6,113,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.6.3多变量离散系统的辨识,离散传递函数矩阵模型其中例3-34,2020/6/6,114,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解,2020/6/6,115,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出的高阶模型应该最小实现辨识结果,2020/6/6,116,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,3.6.4离散系统的递推最小二乘辨识,在系统运行过程中实时获取系统参数，而不是象前面介绍的方法那样一次性获得模型，适合于变参数模型的实时控制广泛应用于自适应控制渐近地逼近参数真值这里介绍算法，仿真研究将在后面介绍,2020/6/6,117,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,递推最小二乘辨识,传递函数模型差分方程模型待辨识参数,2020/6/6,118,控制系统