控制系统仿真与CAD-chapter-4VIP

2020/6/7,1,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,第4章线性控制系统的计算机辅助分析,薛定宇著控制系统计算机辅助设计-MATLAB语言与应用第二版，清华大学出版社，2006CAI课件开发：张望舒哈尔滨工程大学薛定宇东北大学,2020/6/7,2,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统的分析方法,充分利用计算机对线性系统进行分析更新系统分析的观念求解传统方法难以求解的问题离散系统稳定性如何分析？Nyquist图、Nichols图没有频率信息，如何弥补？高阶系统的根轨迹如何绘制？多变量系统如何进行频域分析？,2020/6/7,3,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,本章主要内容,线性系统定性分析线性系统时域响应解析解法线性系统的数字仿真分析根轨迹分析线性系统频域分析,2020/6/7,4,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1线性系统性质分析,主要内容线性系统稳定性分析线性反馈系统内部稳定性分析线性系统的相似变换线性系统可控性分析线性系统可观测性分析Kalman分解系统状态方程的标准型系统的范数测度及求解,2020/6/7,5,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.1线性系统的稳定性分析,给定线性系统模型，如何分析稳定性？由控制理论可知，用Routh表格可以判定该系统稳定性。EdwardJohnRouth(1831-1907)历史局限性,2020/6/7,6,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,状态方程系统的稳定性,连续线性状态方程解析阶稳定性：矩阵的特征根均有负实部,2020/6/7,7,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散系统的稳定性,离散系统状态方程离散系统时域响应解析阶稳定性判定：所有特征根均在单位圆内,2020/6/7,8,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Routh判据的历史局限性,Routh判据提出时，没有求多项式根的方法现在求解矩阵特征根、求解多项式方程的根轻而易举，无需间接方法Routh判据只能得出是否稳定，进一步信息得不出来，如系统是否振荡离散系统无法由Routh方法直接判定，得借助于Jury判据，更复杂稳定性分析方法不统一,2020/6/7,9,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,基于MATLAB的稳定性判定方法,直接判定状态方程模型由可以求出所有特征根离散系统：传递函数模型：完全同样方法图解判定法连续系统：离散系统：，同时画出单位圆,2020/6/7,10,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-1高阶系统稳定性判定直接分析方法零极点模型,2020/6/7,11,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-2高阶离散单位负反馈系统模型MATLAB求解,2020/6/7,12,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.2线性反馈系统的内部稳定性,输入、输出稳定是不够的，因为若内部信号可能过大，对系统作硬件破坏应该引入内部稳定性概念，保证内部信号也是稳定的。,2020/6/7,13,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,由给定稳定输入到内部信号都稳定的系统称为内部稳定系统传递函数矩阵其中逐一判定每个子传递函数的稳定性很烦琐内部稳定性定理,2020/6/7,14,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,内部稳定性定理,闭环系统内部稳定的充要条件为没有不稳定零点没有不稳定零极点对消第一个条件等效于输入输出稳定性判定第2条件即可可以编写MATLAB函数判定内部稳定性,2020/6/7,15,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,判定的MATLAB函数,2020/6/7,16,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.3线性系统的线性相似变换,系统的状态方程表示称为系统实现不同状态选择下，状态方程不惟一相似变换非奇异矩阵状态变换新状态方程模型,2020/6/7,17,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,状态变换公式MATLAB求解方法,2020/6/7,18,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-3已知系统和转换矩阵MATLAB求解,2020/6/7,19,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,变换结果可见，相似变换能改变系统的结构引入相似变换矩阵，可以将已知系统转换成其他的形式,2020/6/7,20,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.4线性系统的可控性分析,可控性定义系统的可控性就是指系统内部的状态是不是可以由外部输出信号控制的性质，,2020/6/7,21,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,线性系统的可控性判定,可控性判定矩阵基于MATLAB的判定方法构造可控性判定矩阵,2020/6/7,22,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-4离散状态方程的可控性MATLAB求解,2020/6/7,23,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,判定矩阵判定矩阵构造方法这样的判定方法同样适合于连续系统和离散系统。也适用于多变量模型,2020/6/7,24,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,由Gram矩阵判定可控性,引入可控Gram矩阵该矩阵满足Lyapunov方程MATLAB求解矩阵构造,2020/6/7,25,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-5求Gram矩阵MATLAB命令Gram矩阵,2020/6/7,26,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,可控性阶梯分解,对于不完全可控的系统阶梯分解阶梯标准型MATLAB函数调用若原系统状态方程完全可控，则不必分解,2020/6/7,27,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-6不完全可控系统,2020/6/7,28,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.5线性系统的可观测性分析,可观测性定义系统的可观测性就是指系统内部的状态是不是可以由系统输出信号重建起来的性质,2020/6/7,29,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,可观测性判定,判定矩阵等同于系统可控性判定Gram矩阵MATLAB求解,2020/6/7,30,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Gram矩阵满足Lyapunov方程对偶问题,2020/6/7,31,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.6Kalman规范分解,Kalman规范分解,2020/6/7,32,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,子空间示意图,2020/6/7,33,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.6系统状态方程标准型的MATLAB求解,常用标准型单变量系统的标准型MATLAB默认的标准型可控标准型实现可观测标准型实现和Jordan标准型实现多变量系统Leunberge标准型侧重点：如何用MATLAB直接获取标准型,2020/6/7,34,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,单变量系统的标准型,可控标准型可观测标准型,2020/6/7,35,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,可控可观测标准型转换,2020/6/7,36,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,可控标准型和可观测标准型，对偶关系,2020/6/7,37,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Jordan标准型,MATLAB变换,2020/6/7,38,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,多变量系统的Leunberge标准型,由可控性判定矩阵构造矩阵,2020/6/7,39,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出Leunberge变换矩阵编写leunberge.m函数,2020/6/7,40,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB函数清单,2020/6/7,41,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,2020/6/7,42,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,标准型的变换方法总结,可控标准型可观测标准型Jordan标准型Leunberge标准型,2020/6/7,43,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-7求解可观测标准型标准型,2020/6/7,44,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-8已知模型,2020/6/7,45,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.1.7系统的范数测度及求解,系统也有范数,2020/6/7,46,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散系统的范数定义范数的MATLAB求解,2020/6/7,47,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-9已知离散系统模型,2020/6/7,48,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.2线性系统时域响应解析解法,给线性系统一个激励信号，输出是什么？有两大类方法解析解方法求解微分方程、差分方程解析解数值解方法主要内容基于状态方程的解析解方法基于传递函数部分方式展开的解析解方法二阶系统的解析解方法,2020/6/7,49,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.2.1基于状态方程的解析解方法,状态方程模型解析解求解难点,2020/6/7,50,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,状态增广方法,消除B矩阵，变成自治系统增广状态方程自治系统可以直接求解析解,2020/6/7,51,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,一般输入信号的系统增广,一般输入信号模型引入增广状态变量,2020/6/7,52,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,增广状态方程模型其中解析解,2020/6/7,53,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB实现函数,2020/6/7,54,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,调用格式信号描述,2020/6/7,55,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-10连续系统模型初值输入信号求解析解,2020/6/7,56,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统增广增广模型,2020/6/7,57,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,解析解求解解析解求解结果稳定性,2020/6/7,58,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.2.2基于部分分式展开方法求解,连续系统的解析解法无重根时部分方式展开,2020/6/7,59,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,由Laplace反变换求解析解有重根时相应项的解析解为,2020/6/7,60,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,部分分式的MATLAB求解例4-11输入信号为阶跃信号输出信号计算,2020/6/7,61,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解解析解解析解精确值,2020/6/7,62,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-12带有复数极点的系统阶跃响应解析解解析解,2020/6/7,63,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,解析解的进一步化简,基于Euler公式的化简其中新MATLAB函数,2020/6/7,64,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,新MATLAB函数清单,2020/6/7,65,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-13仍考虑MATLAB求解解析解,2020/6/7,66,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,基于Laplace变换的求解,参附录A步骤：定义符号变量描述原函数表达式调用laplace()函数或ilaplace()函数求解结果化简，如simple()函数求解举例,2020/6/7,67,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例1MATLAB求解解析解,2020/6/7,68,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例2MATLAB求解解析解,2020/6/7,69,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,离散系统的解析解法,Z变换无重根时部分分式展开解析解,2020/6/7,70,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,考虑采样周期例4-14,2020/6/7,71,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,输出信号解析解Z变换求解步骤定义符号变量调用iztrans()函数求解化简,2020/6/7,72,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,利用符号运算工具箱求解求解结果方法更规范，结果更简单,2020/6/7,73,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,有重根问题的解析解,部分分式表达式的Z反变换例4-15部分分式展开,2020/6/7,74,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,部分分式展开解析解,2020/6/7,75,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,符号运算求解解析解更直观，不建议用前者求解，而直接采用Z变换的符号运算方法求解,2020/6/7,76,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,时间延迟系统的解析解法,例4-16,2020/6/7,77,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,无延迟解析解有延迟解析解,2020/6/7,78,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.2.3二阶系统的阶跃响应及阶跃响应指标,二阶系统模型闭环模型记则,2020/6/7,79,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,阶跃响应的解析解,无阻尼振荡欠阻尼振荡临界阻尼振荡过阻尼振荡,2020/6/7,80,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,二阶系统阶跃响应曲线,2020/6/7,81,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,利用图形绘制功能，从新角度研究同样的问题,三维曲面绘制,2020/6/7,82,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,阶跃响应指标,超调量稳态值上升时间调节时间好的伺服控制系统，应该具有稳态误差小或没有稳态误差、超调量小或没有超调量、上升时间短、调节时间短等性能,2020/6/7,83,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.3线性系统的数字仿真分析,线性系统的解析解可以求解的条件4阶以上的系统需要求解4阶以上的多项式方程，根据Abel定理，无解析解。解析解和数值解结合实际应用需要数值解，需要阶跃响应曲线主要内容线性系统的阶跃响应与脉冲响应任意输入下系统的响应降阶模型的时域分析及比较,2020/6/7,84,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.3.1线性系统的阶跃响应与脉冲响应,阶跃响应曲线绘制函数多系统曲线绘制,2020/6/7,85,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-17延迟系统MATLAB语句利用MATLAB提供的功能，可以从曲线上得到更多的信息，如超调量等,2020/6/7,86,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解解析解解析解数值解精度比较,2020/6/7,87,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-18离散化采样周期求解得出的曲线可以比较,2020/6/7,88,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,ZOH变换Tustin变换，不同采样周期,2020/6/7,89,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-19多变量系统，阶跃响应MATLAB求解语句,2020/6/7,90,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统藕合的概念静态前置补偿矩阵不能直接乘法运算Pade近似,2020/6/7,91,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,补偿后系统的模型解藕效果还可以使得多变量系统能直接设计，在设计前必须解藕。后面将介绍解藕的频域方法,2020/6/7,92,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,系统的脉冲响应曲线,MATLAB下的impulse()函数与step()函数调用结构完全一致MATLAB求解可以容易地研究系统的脉冲响应曲线,2020/6/7,93,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.3.2任意输入下系统的响应,可以利用step()和impulse()函数求解输出信号计算如R(s)已知，则可以直接求解例4-20斜坡响应,2020/6/7,94,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解其他输入的响应可以由lsim()函数求取,2020/6/7,95,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-21多变量系统输入信号MATLAB求解,2020/6/7,96,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,多变量系统的时域响应可以这样求解比较容易理解曲线含义,2020/6/7,97,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.3.3降阶模型的时域分析及比较,前面介绍了降阶方法，但未比较效果例4-22MATLAB求解,2020/6/7,98,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-23非最小相位系统MATLAB求解降阶模型,2020/6/7,99,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.4根轨迹分析,单位负反馈闭环系统特征方程对K的不同取值，则可能绘制出每个特征根变化的曲线，这样的曲线称为系统的根轨迹。根轨迹用开环信息研究闭环特性,2020/6/7,100,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解该函数可以用于单变量不含有时间延迟的连续、离散系统的根轨迹绘制，也可以用于带有时间延迟的单变量离散系统的根轨迹绘制。,2020/6/7,101,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-24开环系统MATLAB求解如何求解临界增益？闭环系统稳定性如何变化,2020/6/7,102,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-25根轨迹求解求出阻尼在处的增益临界增益处阶跃响应,2020/6/7,103,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-26离散系统根轨迹根轨迹绘制,2020/6/7,104,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-27离散系统模型MATLAB求解临界增益求取,2020/6/7,105,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,带延迟的离散系统根轨迹,假设延迟为6步，则可以求临界增益延迟系统临界增益减小,2020/6/7,106,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-28延迟状态方程无延迟系统根轨迹,2020/6/7,107,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,延迟系统用Pade近似处理例4-29正反馈系统MATLAB求解,2020/6/7,108,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.5线性系统频域分析,频域分析Nyquist1932Bode，Nichols提出的新图形方法主要内容单变量系统的频域分析利用频率特性分析系统的稳定性系统的幅值裕度和相位裕度多变量系统的频域分析,2020/6/7,109,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.5.1单变量系统的频域分析,三种表示方法实部与虚部关系曲线即为Nyquist图Nyquist图的缺陷：无对应频率信息横轴对数坐标rad/s，纵轴分贝、度，Bode图幅值与相位关系，Nichols图，无频率信息,2020/6/7,110,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Nyquist曲线绘制grid命令绘制等M和等N圆,2020/6/7,111,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Bode图绘制Nichols图由nichols()函数绘制可以同样处理连续、离散、延迟、多变量系统，格式不变,2020/6/7,112,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-30开环传递函数Nyquist曲线绘制MATLAB曲线特色读取频率信息；频率范围,2020/6/7,113,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Bode图绘制快捷菜单读取特性Nichols图的绘制用鼠标读取频率信息弥补了传统Nichols图的不同,其他频域响应曲线,2020/6/7,114,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-31对下面模型离散化，MATLAB求解不同采样周期的离散模型Bode图,2020/6/7,115,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-32离散系统Nyquist图与Nichols图,2020/6/7,116,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-33延迟系统模型MATLAB求解,2020/6/7,117,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.5.2利用频率特性分析系统的稳定性,Nyquist定理可以进一步解释为,2020/6/7,118,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,可以用开环的系统模型，绘制Nyquist图并以此分析闭环系统的稳定性。,2020/6/7,119,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-34Nyquist图闭环阶跃响应,2020/6/7,120,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.5.3系统的幅值裕度和相位裕度,幅值裕度和相位裕度,相位裕度,幅值裕度,2020/6/7,121,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,稳定性裕度分析,如果系统的Nyquist图不与负实轴相交，则系统的幅值裕度为无穷大。如果系统的Nyquist图不与单位圆相交，则系统的相位裕度为无穷大。,2020/6/7,122,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,如果系统的Nyquist图在第三象限与单位圆有若干个交点，则系统的相位裕度以与离负实轴最近的为准。MATLAB求解方法如果某个裕度为无穷大，则返回Inf，相应的频率值为NaN。,2020/6/7,123,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-35MATLAB求解由于幅相裕度小，系统闭环响应有强振荡,2020/6/7,124,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,4.5.4多变量系统的频域分析,例4-36多变量系统的Nyquist图用nyquist函数直接求解,2020/6/7,125,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,多变量系统分析概述,前面的Nyquist图对多变量系统分析没有太大帮助，所以一般不采用这样的方法英国学派的频域方法SirHowardHRosenbrock教授提出的逆Nyquist阵列的方法(INA方法)剑桥大学SirMacFarlane教授特征轨迹方法帝国理工SirDQMayne教授序贯设计方法Sheffield大学的Owens教授的并矢算法,2020/6/7,126,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MFD工具箱,英国剑桥大学的Maciejowski教授开发基于MATLAB的工具箱多变量系统的描述还可以用传递函数描述，但需要已知公分母。所以过程烦琐。可以求出系统的传递函数矩阵模型,2020/6/7,127,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-37多变量模型传递函数矩阵变换,2020/6/7,128,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,得出公分母分子矩阵用这样的方法可以得出传递函数矩阵模型可以得出MFD能使用的模型,2020/6/7,129,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,对角优势分析,多变量频域分析的最重要内容是系统模型是不是解藕的模型，如果不是则需要变换如何判定是否解藕？闭环系统传递函数矩阵回差矩阵,2020/6/7,130,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,利用回差矩阵的逆矩阵性质，所以在频域分析中用逆的Nyquist矩阵分析更方便Rosenbrock教授采用逆Nyquist阵列方法单变量系统，Nyquist图是研究包围(-1,j0)点的周数来研究稳定性的多变量回差矩阵，研究包围(0,j0)点的情形Gershgorin定理可以分析对角占优性质，从而对系统的藕合进行分析，可以用于多变量系统的分析,2020/6/7,131,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,Gershgorin定理,复数矩阵对角占优矩阵,2020/6/7,132,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,进一步减小半径,2020/6/7,133,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,假设在w下，多变量系统前向回路INA为Gershgorin带，对不同的w值若对全部的w来说，各个对角元素的Gershgorin带均不包含圆心，则称原系统为对角占优系统。显而易见，对角优势矩阵的特征根不位于原点处，则单位反馈的闭环系统是稳定的。,2020/6/7,134,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB求解,MFD工具箱的频域响应数据INA绘制,2020/6/7,135,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,MATLAB函数编写,2020/6/7,136,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,由该函数可以直接绘制多变量系统的INA图，并分析其对角占优性质。采用更小的半径，非传统半径。,2020/6/7,137,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-38多变量模型MATLAB求解,2020/6/7,138,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,前置补偿矩阵MATLAB求解对角占优性质明显,2020/6/7,139,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,其他频域响应数据生成方法,其他,2020/6/7,140,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,例4-39多变量延迟系统MATLAB求解,2020/6/7,141,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,校正后特性绘制Gershgorin带的稳定性判定定理,2020/6/7,142,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,多变量系统的稳定性,2020/6/7,143,控制系统计算机辅助设计-MATLAB语言与应用东北大学信息学院,多变量系统的奇异值曲线绘制,单变量系统有Bode图，多变量系统能否采用这样的方法分析？传递函数矩阵的奇异值可以作为轨迹绘制出来，称为奇异值曲线奇异值曲线是多变量系统鲁棒控制