第-4-章-线性控制系统的计算机辅助分析-《控制系统仿真与CAD》课课件

第4章

线性控制系统的计算机辅助分析薛定宇著《控制系统计算机辅助设计---MATLAB语言与应用》第二版，清华大学出版社，2006CAI课件开发：张望舒哈尔滨工程大学薛定宇东北大学12/11/20221控制系统计算机辅助设计---MATLAB语言与应用第4章

线性控制系统的计算机辅助分析薛定宇著《控制系统计系统的分析方法充分利用计算机对线性系统进行分析更新系统分析的观念求解传统方法难以求解的问题离散系统稳定性如何分析？Nyquist图、Nichols图没有频率信息，如何弥补？高阶系统的根轨迹如何绘制？多变量系统如何进行频域分析？12/11/20222控制系统计算机辅助设计---MATLAB语言与应用系统的分析方法充分利用计算机对线性系统进行分析12/11/2本章主要内容线性系统定性分析线性系统时域响应解析解法线性系统的数字仿真分析根轨迹分析线性系统频域分析12/11/20223控制系统计算机辅助设计---MATLAB语言与应用本章主要内容线性系统定性分析12/11/20223控制系统计4.1线性系统性质分析主要内容线性系统稳定性分析线性反馈系统内部稳定性分析线性系统的相似变换线性系统可控性分析线性系统可观测性分析Kalman分解系统状态方程的标准型系统的范数测度及求解12/11/20224控制系统计算机辅助设计---MATLAB语言与应用4.1线性系统性质分析主要内容12/11/20224控制4.1.1线性系统的稳定性分析给定线性系统模型，如何分析稳定性？由控制理论可知，用Routh表格可以判定该系统稳定性。EdwardJohnRouth(1831-1907)历史局限性12/11/20225控制系统计算机辅助设计---MATLAB语言与应用4.1.1线性系统的稳定性分析给定线性系统模型，如何分析状态方程系统的稳定性连续线性状态方程解析阶稳定性：矩阵的特征根均有负实部12/11/20226控制系统计算机辅助设计---MATLAB语言与应用状态方程系统的稳定性连续线性状态方程12/11/20226控离散系统的稳定性离散系统状态方程离散系统时域响应解析阶稳定性判定：所有特征根均在单位圆内12/11/20227控制系统计算机辅助设计---MATLAB语言与应用离散系统的稳定性离散系统状态方程12/11/20227控制系Routh判据的历史局限性Routh判据提出时，没有求多项式根的方法现在求解矩阵特征根、求解多项式方程的根轻而易举，无需间接方法Routh判据只能得出是否稳定，进一步信息得不出来，如系统是否振荡离散系统无法由Routh方法直接判定，得借助于Jury判据，更复杂稳定性分析方法不统一12/11/20228控制系统计算机辅助设计---MATLAB语言与应用Routh判据的历史局限性Routh判据提出时，没有求多项基于MATLAB的稳定性判定方法直接判定状态方程模型由

可以求出所有特征根离散系统：传递函数模型：完全同样方法图解判定法连续系统：离散系统：，同时画出单位圆12/11/20229控制系统计算机辅助设计---MATLAB语言与应用基于MATLAB的稳定性判定方法直接判定12/11/20例4-1高阶系统稳定性判定直接分析方法零极点模型12/11/202210控制系统计算机辅助设计---MATLAB语言与应用例4-1高阶系统稳定性判定12/11/202210控制系统例4-2高阶离散单位负反馈系统模型MATLAB求解12/11/202211控制系统计算机辅助设计---MATLAB语言与应用例4-2高阶离散单位负反馈系统模型12/11/202214.1.2线性反馈系统的内部稳定性输入、输出稳定是不够的，因为若内部信号可能过大，对系统作硬件破坏应该引入内部稳定性概念，保证内部信号也是稳定的。12/11/202212控制系统计算机辅助设计---MATLAB语言与应用4.1.2线性反馈系统的内部稳定性输入、输出稳定是不够的由给定稳定输入到内部信号都稳定的系统称为内部稳定系统传递函数矩阵其中逐一判定每个子传递函数的稳定性很烦琐内部稳定性定理12/11/202213控制系统计算机辅助设计---MATLAB语言与应用由给定稳定输入到内部信号12/11/2内部稳定性定理闭环系统内部稳定的充要条件为没有不稳定零点没有不稳定零极点对消第一个条件等效于输入输出稳定性判定第2条件即可可以编写MATLAB函数判定内部稳定性12/11/202214控制系统计算机辅助设计---MATLAB语言与应用内部稳定性定理闭环系统内部稳定的充要条件为12/11/202判定的MATLAB函数12/11/202215控制系统计算机辅助设计---MATLAB语言与应用判定的MATLAB函数12/11/202215控制系统计4.1.3线性系统的线性相似变换系统的状态方程表示称为系统实现不同状态选择下，状态方程不惟一相似变换非奇异矩阵状态变换新状态方程模型12/11/202216控制系统计算机辅助设计---MATLAB语言与应用4.1.3线性系统的线性相似变换系统的状态方程表示称为系状态变换公式MATLAB求解方法12/11/202217控制系统计算机辅助设计---MATLAB语言与应用状态变换公式12/11/202217控制系统计算机辅助设计-例4-3已知系统和转换矩阵MATLAB求解12/11/202218控制系统计算机辅助设计---MATLAB语言与应用例4-3已知系统和转换矩阵12/11/202218控制系统变换结果可见，相似变换能改变系统的结构引入相似变换矩阵，可以将已知系统转换成其他的形式12/11/202219控制系统计算机辅助设计---MATLAB语言与应用变换结果12/11/202219控制系统计算机辅助设计---4.1.4线性系统的可控性分析可控性定义

系统的可控性就是指系统内部的状态是不是可以由外部输出信号控制的性质，12/11/202220控制系统计算机辅助设计---MATLAB语言与应用4.1.4线性系统的可控性分析可控性定义12/11/20线性系统的可控性判定可控性判定矩阵

基于MATLAB的判定方法构造可控性判定矩阵12/11/202221控制系统计算机辅助设计---MATLAB语言与应用线性系统的可控性判定可控性判定矩阵12/11/202221控例4-4离散状态方程的可控性MATLAB求解12/11/202222控制系统计算机辅助设计---MATLAB语言与应用例4-4离散状态方程的可控性12/11/202222控制系判定矩阵判定矩阵构造方法这样的判定方法同样适合于连续系统和离散系统。也适用于多变量模型12/11/202223控制系统计算机辅助设计---MATLAB语言与应用判定矩阵12/11/202223控制系统计算机辅助设计---由Gram矩阵判定可控性引入可控Gram矩阵该矩阵满足Lyapunov方程MATLAB求解矩阵构造12/11/202224控制系统计算机辅助设计---MATLAB语言与应用由Gram矩阵判定可控性引入可控Gram矩阵12/1例4-5求Gram矩阵MATLAB命令Gram矩阵12/11/202225控制系统计算机辅助设计---MATLAB语言与应用例4-5求Gram矩阵12/11/202225控制系统可控性阶梯分解对于不完全可控的系统阶梯分解阶梯标准型MATLAB函数调用若原系统状态方程完全可控，则不必分解12/11/202226控制系统计算机辅助设计---MATLAB语言与应用可控性阶梯分解对于不完全可控的系统阶梯分解12/11/202例4-6不完全可控系统12/11/202227控制系统计算机辅助设计---MATLAB语言与应用例4-6不完全可控系统12/11/202227控制系统计算4.1.5线性系统的可观测性分析可观测性定义

系统的可观测性就是指系统内部的状态是不是可以由系统输出信号重建起来的性质12/11/202228控制系统计算机辅助设计---MATLAB语言与应用4.1.5线性系统的可观测性分析可观测性定义12/11/可观测性判定判定矩阵等同于系统可控性判定Gram矩阵MATLAB求解12/11/202229控制系统计算机辅助设计---MATLAB语言与应用可观测性判定判定矩阵12/11/202229控制系统计算机辅Gram矩阵满足Lyapunov方程对偶问题12/11/202230控制系统计算机辅助设计---MATLAB语言与应用Gram矩阵满足Lyapunov方程12/11/2024.1.6Kalman规范分解Kalman规范分解12/11/202231控制系统计算机辅助设计---MATLAB语言与应用4.1.6Kalman规范分解Kalman规范分解子空间示意图12/11/202232控制系统计算机辅助设计---MATLAB语言与应用子空间12/11/202232控制系统计算机辅助设计---M4.1.6系统状态方程标准型的MATLAB求解常用标准型单变量系统的标准型MATLAB默认的标准型可控标准型实现可观测标准型实现和Jordan标准型实现多变量系统Leunberge标准型侧重点：如何用MATLAB直接获取标准型12/11/202233控制系统计算机辅助设计---MATLAB语言与应用4.1.6系统状态方程标准型的MATLAB求解常用标准单变量系统的标准型可控标准型可观测标准型12/11/202234控制系统计算机辅助设计---MATLAB语言与应用单变量系统的标准型可控标准型12/11/202234控制系统可控可观测标准型转换12/11/202235控制系统计算机辅助设计---MATLAB语言与应用可控可观测标准型转换12/11/202235控制系统计算机辅可控标准型和可观测标准型，对偶关系12/11/202236控制系统计算机辅助设计---MATLAB语言与应用12/11/202236控制系统计算机辅助设计---MATLJordan标准型

MATLAB变换12/11/202237控制系统计算机辅助设计---MATLAB语言与应用Jordan标准型12/11/202237控制系统计算机多变量系统的Leunberge标准型由可控性判定矩阵构造矩阵12/11/202238控制系统计算机辅助设计---MATLAB语言与应用多变量系统的Leunberge标准型由可控性判定矩阵12得出Leunberge变换矩阵编写leunberge.m函数12/11/202239控制系统计算机辅助设计---MATLAB语言与应用得出Leunberge变换矩阵12/11/202239控MATLAB函数清单12/11/202240控制系统计算机辅助设计---MATLAB语言与应用MATLAB函数清单12/11/202240控制系统计算机12/11/202241控制系统计算机辅助设计---MATLAB语言与应用12/11/202241控制系统计算机辅助设计---MATL标准型的变换方法总结可控标准型可观测标准型Jordan标准型Leunberge标准型12/11/202242控制系统计算机辅助设计---MATLAB语言与应用标准型的变换方法总结可控标准型12/11/202242控制系例4-7求解可观测标准型标准型12/11/202243控制系统计算机辅助设计---MATLAB语言与应用例4-712/11/202243控制系统计算机辅助设计---例4-8已知模型12/11/202244控制系统计算机辅助设计---MATLAB语言与应用例4-8已知模型12/11/202244控制系统计算机辅助4.1.7系统的范数测度及求解系统也有范数

12/11/202245控制系统计算机辅助设计---MATLAB语言与应用4.1.7系统的范数测度及求解系统也有范数12/11/2离散系统的范数定义范数的MATLAB求解12/11/202246控制系统计算机辅助设计---MATLAB语言与应用离散系统的范数定义12/11/202246控制系统计算机辅助例4-9已知离散系统模型12/11/202247控制系统计算机辅助设计---MATLAB语言与应用例4-9已知离散系统模型12/11/202247控制系统计4.2线性系统时域响应解析解法给线性系统一个激励信号，输出是什么？有两大类方法解析解方法求解微分方程、差分方程解析解数值解方法主要内容基于状态方程的解析解方法基于传递函数部分方式展开的解析解方法二阶系统的解析解方法12/11/202248控制系统计算机辅助设计---MATLAB语言与应用4.2线性系统时域响应解析解法给线性系统一个激励信号，输4.2.1基于状态方程的解析解方法状态方程模型解析解求解难点12/11/202249控制系统计算机辅助设计---MATLAB语言与应用4.2.1基于状态方程的解析解方法状态方程模型12/11状态增广方法消除B矩阵，变成自治系统

增广状态方程自治系统可以直接求解析解12/11/202250控制系统计算机辅助设计---MATLAB语言与应用状态增广方法消除B矩阵，变成自治系统12/11/20225一般输入信号的系统增广一般输入信号模型引入增广状态变量12/11/202251控制系统计算机辅助设计---MATLAB语言与应用一般输入信号的系统增广一般输入信号模型12/11/20225增广状态方程模型其中解析解12/11/202252控制系统计算机辅助设计---MATLAB语言与应用增广状态方程模型12/11/202252控制系统计算机辅助设MATLAB实现函数12/11/202253控制系统计算机辅助设计---MATLAB语言与应用MATLAB实现函数12/11/202253控制系统计算机调用格式信号描述12/11/202254控制系统计算机辅助设计---MATLAB语言与应用12/11/202254控制系统计算机辅助设计---MATL例4-10连续系统模型初值输入信号求解析解12/11/202255控制系统计算机辅助设计---MATLAB语言与应用例4-10连续系统模型12/11/202255控制系统计算系统增广增广模型12/11/202256控制系统计算机辅助设计---MATLAB语言与应用系统增广12/11/202256控制系统计算机辅助设计---解析解求解解析解求解结果稳定性12/11/202257控制系统计算机辅助设计---MATLAB语言与应用解析解求解12/11/202257控制系统计算机辅助设计--4.2.2基于部分分式展开方法求解连续系统的解析解法

无重根时部分方式展开12/11/202258控制系统计算机辅助设计---MATLAB语言与应用4.2.2基于部分分式展开方法求解连续系统的解析解法12由Laplace反变换求解析解有重根时相应项的解析解为12/11/202259控制系统计算机辅助设计---MATLAB语言与应用由Laplace反变换求解析解12/11/202259控部分分式的MATLAB求解例4-11输入信号为阶跃信号输出信号计算12/11/202260控制系统计算机辅助设计---MATLAB语言与应用部分分式的MATLAB求解12/11/202260控制系MATLAB求解解析解解析解精确值12/11/202261控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解12/11/202261控制系统计算机辅助例4-12带有复数极点的系统阶跃响应解析解解析解12/11/202262控制系统计算机辅助设计---MATLAB语言与应用例4-12带有复数极点的系统12/11/202262控制解析解的进一步化简基于Euler公式的化简其中新MATLAB函数

12/11/202263控制系统计算机辅助设计---MATLAB语言与应用解析解的进一步化简基于Euler公式的化简12/11/2新MATLAB函数清单12/11/202264控制系统计算机辅助设计---MATLAB语言与应用新MATLAB函数清单12/11/202264控制系统计例4-13仍考虑MATLAB求解解析解12/11/202265控制系统计算机辅助设计---MATLAB语言与应用例4-13仍考虑12/11/202265控制系统计算机辅基于Laplace变换的求解参附录A步骤：定义符号变量描述原函数表达式调用laplace()函数或ilaplace()函数求解结果化简，如simple()函数求解举例12/11/202266控制系统计算机辅助设计---MATLAB语言与应用基于Laplace变换的求解参附录A12/11/202例1MATLAB求解解析解12/11/202267控制系统计算机辅助设计---MATLAB语言与应用例112/11/202267控制系统计算机辅助设计---MA例2MATLAB求解解析解12/11/202268控制系统计算机辅助设计---MATLAB语言与应用例212/11/202268控制系统计算机辅助设计---MA离散系统的解析解法Z变换无重根时部分分式展开解析解12/11/202269控制系统计算机辅助设计---MATLAB语言与应用离散系统的解析解法Z变换12/11/202269控制系统计考虑采样周期例4-1412/11/202270控制系统计算机辅助设计---MATLAB语言与应用考虑采样周期12/11/202270控制系统计算机辅助设计-输出信号解析解Z变换求解步骤定义符号变量调用iztrans()函数求解化简12/11/202271控制系统计算机辅助设计---MATLAB语言与应用输出信号12/11/202271控制系统计算机辅助设计---利用符号运算工具箱求解求解结果方法更规范，结果更简单12/11/202272控制系统计算机辅助设计---MATLAB语言与应用利用符号运算工具箱求解12/11/202272控制系统计算机有重根问题的解析解部分分式表达式的Z反变换例4-15部分分式展开12/11/202273控制系统计算机辅助设计---MATLAB语言与应用有重根问题的解析解部分分式表达式的Z反变换12/11/202部分分式展开解析解12/11/202274控制系统计算机辅助设计---MATLAB语言与应用部分分式展开12/11/202274控制系统计算机辅助设计-符号运算求解解析解更直观，不建议用前者求解，而直接采用Z变换的符号运算方法求解12/11/202275控制系统计算机辅助设计---MATLAB语言与应用符号运算求解12/11/202275控制系统计算机辅助设计-时间延迟系统的解析解法

例4-1612/11/202276控制系统计算机辅助设计---MATLAB语言与应用时间延迟系统的解析解法12/11/202276控制系统计算无延迟解析解有延迟解析解12/11/202277控制系统计算机辅助设计---MATLAB语言与应用无延迟解析解12/11/202277控制系统计算机辅助设计-4.2.3二阶系统的阶跃响应及

阶跃响应指标二阶系统模型闭环模型记则12/11/202278控制系统计算机辅助设计---MATLAB语言与应用4.2.3二阶系统的阶跃响应及

阶跃响应指标二阶跃响应的解析解无阻尼振荡欠阻尼振荡临界阻尼振荡过阻尼振荡12/11/202279控制系统计算机辅助设计---MATLAB语言与应用阶跃响应的解析解无阻尼振荡12/11/202279控制系统计二阶系统阶跃响应曲线12/11/202280控制系统计算机辅助设计---MATLAB语言与应用二阶系统阶跃响应曲线12/11/202280控制系统计算机辅利用图形绘制功能，从新角度研究同样的问题三维曲面绘制12/11/202281控制系统计算机辅助设计---MATLAB语言与应用利用图形绘制功能，从新角度研究同样的问题三维曲面绘制12/1阶跃响应指标超调量稳态值上升时间调节时间好的伺服控制系统，应该具有稳态误差小或没有稳态误差、超调量小或没有超调量、上升时间短、调节时间短等性能12/11/202282控制系统计算机辅助设计---MATLAB语言与应用阶跃响应指标超调量12/11/202282控制系统计算机辅助4.3线性系统的数字仿真分析线性系统的解析解可以求解的条件4阶以上的系统需要求解4阶以上的多项式方程，根据Abel定理，无解析解。解析解和数值解结合实际应用需要数值解，需要阶跃响应曲线主要内容线性系统的阶跃响应与脉冲响应任意输入下系统的响应降阶模型的时域分析及比较12/11/202283控制系统计算机辅助设计---MATLAB语言与应用4.3线性系统的数字仿真分析线性系统的解析解可以求解的条4.3.1线性系统的阶跃响应与脉冲响应阶跃响应曲线绘制函数多系统曲线绘制12/11/202284控制系统计算机辅助设计---MATLAB语言与应用4.3.1线性系统的阶跃响应与脉冲响应阶跃响应曲线绘制函例4-17延迟系统MATLAB语句利用MATLAB提供的功能，可以从曲线上得到更多的信息，如超调量等12/11/202285控制系统计算机辅助设计---MATLAB语言与应用例4-17延迟系统12/11/202285控制系统计算机辅MATLAB求解解析解解析解数值解精度比较12/11/202286控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解解析解12/11/202286控制系统计算例4-18离散化采样周期求解得出的曲线可以比较12/11/202287控制系统计算机辅助设计---MATLAB语言与应用例4-18离散化12/11/202287控制系统计算机辅助ZOH变换Tustin变换，不同采样周期12/11/202288控制系统计算机辅助设计---MATLAB语言与应用ZOH变换12/11/202288控制系统计算机辅助设计-例4-19多变量系统，阶跃响应

MATLAB求解语句12/11/202289控制系统计算机辅助设计---MATLAB语言与应用例4-19多变量系统，阶跃响应12/11/202289控制系统藕合的概念静态前置补偿矩阵不能直接乘法运算Pade近似12/11/202290控制系统计算机辅助设计---MATLAB语言与应用系统藕合的概念12/11/202290控制系统计算机辅助设计补偿后系统的模型解藕效果还可以使得多变量系统能直接设计，在设计前必须解藕。后面将介绍解藕的频域方法12/11/202291控制系统计算机辅助设计---MATLAB语言与应用补偿后系统的模型12/11/202291控制系统计算机辅助设系统的脉冲响应曲线MATLAB下的impulse()函数与step()函数调用结构完全一致MATLAB求解可以容易地研究系统的脉冲响应曲线12/11/202292控制系统计算机辅助设计---MATLAB语言与应用系统的脉冲响应曲线MATLAB下的impulse()4.3.2任意输入下系统的响应可以利用step()和impulse()函数求解输出信号计算如R(s)已知，则可以直接求解例4-20斜坡响应12/11/202293控制系统计算机辅助设计---MATLAB语言与应用4.3.2任意输入下系统的响应可以利用step()

MATLAB求解其他输入的响应可以由lsim()函数求取12/11/202294控制系统计算机辅助设计---MATLAB语言与应用12/11/202294控制系统计算机辅助设计---MAT例4-21多变量系统输入信号MATLAB求解12/11/202295控制系统计算机辅助设计---MATLAB语言与应用例4-21多变量系统12/11/202295控制系统计算机多变量系统的时域响应可以这样求解比较容易理解曲线含义12/11/202296控制系统计算机辅助设计---MATLAB语言与应用12/11/202296控制系统计算机辅助设计---MATL4.3.3降阶模型的时域分析及比较前面介绍了降阶方法，但未比较效果例4-22MATLAB求解12/11/202297控制系统计算机辅助设计---MATLAB语言与应用4.3.3降阶模型的时域分析及比较前面介绍了降阶方法，但例4-23非最小相位系统MATLAB求解降阶模型12/11/202298控制系统计算机辅助设计---MATLAB语言与应用例4-23非最小相位系统12/11/202298控制系统4.4根轨迹分析

单位负反馈闭环系统特征方程对K的不同取值，则可能绘制出每个特征根变化的曲线，这样的曲线称为系统的根轨迹。根轨迹用开环信息研究闭环特性12/11/202299控制系统计算机辅助设计---MATLAB语言与应用4.4根轨迹分析12/11/202299控制系统计算机MATLAB求解该函数可以用于单变量不含有时间延迟的连续、离散系统的根轨迹绘制，也可以用于带有时间延迟的单变量离散系统的根轨迹绘制。12/11/2022100控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解12/11/2022100控制系统计算机辅例4-24开环系统MATLAB求解如何求解临界增益？闭环系统稳定性如何变化12/11/2022101控制系统计算机辅助设计---MATLAB语言与应用例4-24开环系统12/11/2022101控制系统计算例4-25根轨迹求解求出阻尼在处的增益临界增益处阶跃响应12/11/2022102控制系统计算机辅助设计---MATLAB语言与应用例4-2512/11/2022102控制系统计算机辅助设计例4-26离散系统根轨迹根轨迹绘制12/11/2022103控制系统计算机辅助设计---MATLAB语言与应用例4-26离散系统根轨迹12/11/2022103控制系例4-27离散系统模型MATLAB求解临界增益求取12/11/2022104控制系统计算机辅助设计---MATLAB语言与应用例4-27离散系统模型12/11/2022104控制系统带延迟的离散系统根轨迹假设延迟为6步，则可以求临界增益延迟系统临界增益减小12/11/2022105控制系统计算机辅助设计---MATLAB语言与应用带延迟的离散系统根轨迹假设延迟为6步，则12/11/20例4-28延迟状态方程无延迟系统根轨迹12/11/2022106控制系统计算机辅助设计---MATLAB语言与应用例4-28延迟状态方程12/11/2022106控制系统延迟系统用Pade近似处理例4-29正反馈系统MATLAB求解12/11/2022107控制系统计算机辅助设计---MATLAB语言与应用延迟系统用Pade近似处理12/11/2022107控制4.5线性系统频域分析频域分析Nyquist1932Bode，Nichols提出的新图形方法主要内容单变量系统的频域分析利用频率特性分析系统的稳定性系统的幅值裕度和相位裕度多变量系统的频域分析12/11/2022108控制系统计算机辅助设计---MATLAB语言与应用4.5线性系统频域分析频域分析12/11/2022108控4.5.1单变量系统的频域分析

三种表示方法

实部与虚部关系曲线即为Nyquist图Nyquist图的缺陷：无对应频率信息

横轴对数坐标rad/s，纵轴分贝、度，Bode图幅值与相位关系，Nichols图，无频率信息12/11/2022109控制系统计算机辅助设计---MATLAB语言与应用4.5.1单变量系统的频域分析12/11/202210Nyquist曲线绘制grid命令绘制等M和等N圆12/11/2022110控制系统计算机辅助设计---MATLAB语言与应用Nyquist曲线绘制12/11/2022110控制系统Bode图绘制Nichols图由nichols()函数绘制可以同样处理连续、离散、延迟、多变量系统，格式不变12/11/2022111控制系统计算机辅助设计---MATLAB语言与应用Bode图绘制12/11/2022111控制系统计算机辅助例4-30开环传递函数Nyquist曲线绘制MATLAB曲线特色读取频率信息；频率范围12/11/2022112控制系统计算机辅助设计---MATLAB语言与应用例4-30开环传递函数12/11/2022112控制系统Bode图绘制快捷菜单读取特性Nichols图的绘制用鼠标读取频率信息弥补了传统Nichols图的不同其他频域响应曲线12/11/2022113控制系统计算机辅助设计---MATLAB语言与应用Bode图绘制其他频域响应曲线12/11/2022113控例4-31对下面模型离散化，MATLAB求解不同采样周期的离散模型Bode图12/11/2022114控制系统计算机辅助设计---MATLAB语言与应用例4-31对下面模型离散化，12/11/2022114控例4-32离散系统Nyquist图与Nichols图12/11/2022115控制系统计算机辅助设计---MATLAB语言与应用例4-32离散系统12/11/2022115控制系统计算例4-33延迟系统模型MATLAB求解12/11/2022116控制系统计算机辅助设计---MATLAB语言与应用例4-33延迟系统模型12/11/2022116控制系统4.5.2利用频率特性分析系统

的稳定性

Nyquist定理可以进一步解释为

12/11/2022117控制系统计算机辅助设计---MATLAB语言与应用4.5.2利用频率特性分析系统

的稳定性12/11/2

可以用开环的系统模型，绘制Nyquist图并以此分析闭环系统的稳定性。12/11/2022118控制系统计算机辅助设计---MATLAB语言与应用12/11/2022118控制系统计算机辅助设计---MA例4-34Nyquist图闭环阶跃响应12/11/2022119控制系统计算机辅助设计---MATLAB语言与应用例4-3412/11/2022119控制系统计算机辅助设计4.5.3系统的幅值裕度和相位裕度幅值裕度和相位裕度相位裕度幅值裕度12/11/2022120控制系统计算机辅助设计---MATLAB语言与应用4.5.3系统的幅值裕度和相位裕度幅值裕度和相位裕度相稳定性裕度分析如果系统的Nyquist图不与负实轴相交，则系统的幅值裕度为无穷大。

如果系统的Nyquist图不与单位圆相交，则系统的相位裕度为无穷大。12/11/2022121控制系统计算机辅助设计---MATLAB语言与应用稳定性裕度分析如果系统的Nyquist图不与负实轴相交，如果系统的Nyquist图在第三象限与单位圆有若干个交点，则系统的相位裕度以与离负实轴最近的为准。MATLAB求解方法如果某个裕度为无穷大，则返回Inf，相应的频率值为NaN。12/11/2022122控制系统计算机辅助设计---MATLAB语言与应用如果系统的Nyquist图在第三象限与单位圆有若干个交点例4-35MATLAB求解

由于幅相裕度小，系统闭环响应有强振荡12/11/2022123控制系统计算机辅助设计---MATLAB语言与应用例4-3512/11/2022123控制系统计算机辅助设计-4.5.4多变量系统的频域分析例4-36多变量系统的Nyquist图用nyquist函数直接求解12/11/2022124控制系统计算机辅助设计---MATLAB语言与应用4.5.4多变量系统的频域分析例4-36多变量系统的N多变量系统分析概述前面的Nyquist图对多变量系统分析没有太大帮助，所以一般不采用这样的方法英国学派的频域方法SirHowardHRosenbrock教授提出的逆Nyquist阵列的方法(INA方法)剑桥大学SirMacFarlane教授特征轨迹方法帝国理工SirDQMayne教授序贯设计方法Sheffield大学的Owens教授的并矢算法12/11/2022125控制系统计算机辅助设计---MATLAB语言与应用多变量系统分析概述前面的Nyquist图对多变量系统分析MFD工具箱英国剑桥大学的Maciejowski教授开发基于MATLAB的工具箱多变量系统的描述还可以用传递函数描述，但需要已知公分母。所以过程烦琐。可以求出系统的传递函数矩阵模型12/11/2022126控制系统计算机辅助设计---MATLAB语言与应用MFD工具箱英国剑桥大学的Maciejowski教授开例4-37多变量模型传递函数矩阵变换12/11/2022127控制系统计算机辅助设计---MATLAB语言与应用例4-37多变量模型12/11/2022127控制系统计算得出公分母分子矩阵用这样的方法可以得出传递函数矩阵模型可以得出MFD能使用的模型12/11/2022128控制系统计算机辅助设计---MATLAB语言与应用得出公分母12/11/2022128控制系统计算机辅助设计-对角优势分析多变量频域分析的最重要内容是系统模型是不是解藕的模型，如果不是则需要变换如何判定是否解藕？

闭环系统传递函数矩阵回差矩阵12/11/2022129控制系统计算机辅助设计---MATLAB语言与应用对角优势分析多变量频域分析的最重要内容是系统模型是不是解藕的利用回差矩阵的逆矩阵性质，所以在频域分析中用逆的Nyquist矩阵分析更方便Rosenbrock教授采用逆Nyquist阵列方法单变量系统，Nyquist图是研究包围(-1,j0)点的周数来研究稳定性的多变量回差矩阵，研究包围(0,j0)点的情形Gershgorin定理可以分析对角占优性质，从而对系统的藕合进行分析，可以用于多变量系统的分析12/11/2022130控制系统计算机辅助设计---MATLAB语言与应用利用回差矩阵的逆矩阵性质，所以在频域分析中用逆的NyquiGershgorin定理复数矩阵

对角占优矩阵12/11/2022131控制系统计算机辅助设计---MATLAB语言与应用Gershgorin定理复数矩阵12/11/2022131

进一步减小半径12/11/2022132控制系统计算机辅助设计---MATLAB语言与应用12/11/2022132控制系统计算机辅助设计---MA假设在w下，多变量系统前向回路INA为Gershgorin带，对不同的w值若对全部的w来说，各个对角元素的Gershgorin带均不包含圆心，则称原系统为对角占优系统。显而易见，对角优势矩阵的特征根不位于原点处，则单位反馈的闭环系统是稳定的。12/11/2022133控制系统计算机辅助设计---MATLAB语言与应用假设在w下，多变量系统前向回路INA为12/11/2MATLAB求解MFD工具箱的频域响应数据INA绘制12/11/2022134控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解MFD工具箱的频域响应数据12/11/MATLAB函数编写12/11/2022135控制系统计算机辅助设计---MATLAB语言与应用MATLAB函数编写12/11/2022135控制系统计算由该函数可以直接绘制多变量系统的INA图，并分析其对角占优性质。采用更小的半径，非传统半径。12/11/2022136控制系统计算机辅助设计---MATLAB语言与应用由该函数可以直接绘制多变量系统的INA图，并分析其对角占优例4-38多变量模型MATLAB求解12/11/2022137控制系统计算机辅助设计---MATLAB语言与应用例4-38多变量模型12/11/2022137控制系统计算前置补偿矩阵MATLAB求解对角占优性质明显12/11/2022138控制系统计算机辅助设计---MATLAB语言与应用前置补偿矩阵12/11/2022138控制系统计算机辅助设计其他频域响应数据生成方法

其他12/11/2022139控制系统计算机辅助设计---MATLAB语言与应用其他频域响应数据生成方法12/11/2022139控制系统例4-39多变量延迟系统MATLAB求解12/11/2022140控制系统计算机辅助设计---MATLAB语言与应用例4-39多变量延迟系统12/11/2022140控制系统

校正后特性绘制Gershgorin带的稳定性判定定理12/11/2022141控制系统计算机辅助设计---MATLAB语言与应用12/11/2022141控制系统计算机辅助设计---MA多变量系统的稳定性12/11/2022142控制系统计算机辅助设计---MATLAB语言与应用多变量系统的稳定性12/11/2022142控制系统计算机辅多变量系统的奇异值曲线绘制单变量系统有Bode图，多变量系统能否采用这样的方法分析？

传递函数矩阵的奇异值可以作为轨迹绘制出来，称为奇异值曲线奇异值曲线是多变量系统鲁棒控制中的重要指标，由sigma()函数绘制12/11/2022143控制系统计算机辅助设计---MATLAB语言与应用多变量系统的奇异值曲线绘制单变量系统有Bode图，多变量例4-40多变量延迟模型MATLAB绘制奇异值曲线12/11/2022144控制系统计算机辅助设计---MATLAB语言与应用例4-40多变量延迟模型12/11/2022144控制系统4.5.5降阶模型的频域分析及比较降阶模型对原模型的频域拟合如何？例4-41可以由前面的方法求出Pade、Routh模型、最优降阶模型、均衡实现降阶模型、最优Hankel范数模型、Schur模型12/11/2022145控制系统计算机辅助设计---MATLAB语言与应用4.5.5降阶模型的频域分析及比较降阶模型对原模型的频域降阶模型与比较最优降阶模型的效果最好12/11/2022146控制系统计算机辅助设计---MATLAB语言与应用降阶模型与比较12/11/2022146控制系统计算机辅助设本章主要内容MATLAB的使用为控制系统的分析提供了有力的工具，在控制系统发展初期，由于没有这样的强有力工具，出现了很多间接的方法，例如控制系统的稳定性分析以往的Routh判据可以完全由直接求根的方法取代，对控制系统来说，用eig()就可以直接求出系统的特征根，并给出了反馈控制系统的内部稳定性概念与判定方法。12/11/2022147控制系统计算机辅助设计---MATLAB语言与应用本章主要内容MATLAB的使用为控制系统的分析提供了有力利用MATLAB这样的工具还可以直接对控制系统的可控性、可观测性等进行直接判定，还介绍了系统的可控性、可观测性阶梯分解、Kalman分解、Leunberge标准型转换等，并介绍系统的范数测度及计算。本章介绍了线性系统的解析解算法，包括基于状态方程的解析解方法和基于部分分式展开技术的解析解方法，分别就连续系统和离散系统等问题进行了探讨，还介绍了改进的部分分式展开方法，从而可以得出更可读的解析解。由典型二阶系统的阶跃响应定义了系统的一些响应指标，如超调量、调节时间等，还介绍了这些指标的求解方法。12/11/2022148控制系统计算机辅助设计---MATLAB语言与应用利用MATLAB这样的工具还可以直接对控制系统的可控性连续和离散系统的阶跃响应曲线可以直接由MATLAB给出的step()函数直接绘制出来，还可以用函数就可以绘制出系统的脉冲响应曲线，还可以用函数都可以用于系统在任意输入下的时域响应，这些函数均可以用于所有能用MATLAB下线性时不变对象描述的线性系统时域分析。根轨迹分析是单变量系统稳定性分析与控制系统校正的一种有用方法，用rlocus()函数就可以直接绘制出单变量连续与离散系统的根轨迹曲线，并可以直接从根轨迹上读取临界稳定增益值。这样的方法还可以直接应用于绘制带有时间延迟的离散系统根轨迹绘制中。12/11/2022149控制系统计算机辅助设计---MATLAB语言与应用连续和离散系统的阶跃响应曲线可以直接由MATLAB给出的从频域响应中复数的几种表示方法引入了Nyquist图、Bode图和Nichols图，并介绍了在MATLAB控制系统工具箱中如何绘制这些图形的方法，介绍了应用频域响应进行闭环系统稳定性分析的方法，还介绍了幅值裕度和相位裕度的求取函数margin()，这些方法可以直接求解连续和离散单变量系统的频域响应分析。介绍了多变量系统的逆Nyquist阵列、回差矩阵与Gershgorin定理，基于这些概念介绍了多变量系统频域响应分析方法及相应的MATLAB实现。12/11/2022150控制系统计算机辅助设计---MATLAB语言与应用从频域响应中复数的几种表示方法引入了Nyquist图、Bo第4章

线性控制系统的计算机辅助分析薛定宇著《控制系统计算机辅助设计---MATLAB语言与应用》第二版，清华大学出版社，2006CAI课件开发：张望舒哈尔滨工程大学薛定宇东北大学12/11/2022151控制系统计算机辅助设计---MATLAB语言与应用第4章

线性控制系统的计算机辅助分析薛定宇著《控制系统计系统的分析方法充分利用计算机对线性系统进行分析更新系统分析的观念求解传统方法难以求解的问题离散系统稳定性如何分析？Nyquist图、Nichols图没有频率信息，如何弥补？高阶系统的根轨迹如何绘制？多变量系统如何进行频域分析？12/11/2022152控制系统计算机辅助设计---MATLAB语言与应用系统的分析方法充分利用计算机对线性系统进行分析12/11/2本章主要内容线性系统定性分析线性系统时域响应解析解法线性系统的数字仿真分析根轨迹分析线性系统频域分析12/11/2022153控制系统计算机辅助设计---MATLAB语言与应用本章主要内容线性系统定性分析12/11/20223控制系统计4.1线性系统性质分析主要内容线性系统稳定性分析线性反馈系统内部稳定性分析线性系统的相似变换线性系统可控性分析线性系统可观测性分析Kalman分解系统状态方程的标准型系统的范数测度及求解12/11/2022154控制系统计算机辅助设计---MATLAB语言与应用4.1线性系统性质分析主要内容12/11/20224控制4.1.1线性系统的稳定性分析给定线性系统模型，如何分析稳定性？由控制理论可知，用Routh表格可以判定该系统稳定性。EdwardJohnRouth(1831-1907)历史局限性12/11/2022155控制系统计算机辅助设计---MATLAB语言与应用4.1.1线性系统的稳定性分析给定线性系统模型，如何分析状态方程系统的稳定性连续线性状态方程解析阶稳定性：矩阵的特征根均有负实部12/11/2022156控制系统计算机辅助设计---MATLAB语言与应用状态方程系统的稳定性连续线性状态方程12/11/20226控离散系统的稳定性离散系统状态方程离散系统时域响应解析阶稳定性判定：所有特征根均在单位圆内12/11/2022157控制系统计算机辅助设计---MATLAB语言与应用离散系统的稳定性离散系统状态方程12/11/20227控制系Routh判据的历史局限性Routh判据提出时，没有求多项式根的方法现在求解矩阵特征根、求解多项式方程的根轻而易举，无需间接方法Routh判据只能得出是否稳定，进一步信息得不出来，如系统是否振荡离散系统无法由Routh方法直接判定，得借助于Jury判据，更复杂稳定性分析方法不统一12/11/2022158控制系统计算机辅助设计---MATLAB语言与应用Routh判据的历史局限性Routh判据提出时，没有求多项基于MATLAB的稳定性判定方法直接判定状态方程模型由

可以求出所有特征根离散系统：传递函数模型：完全同样方法图解判定法连续系统：离散系统：，同时画出单位圆12/11/2022159控制系统计算机辅助设计---MATLAB语言与应用基于MATLAB的稳定性判定方法直接判定12/11/20例4-1高阶系统稳定性判定直接分析方法零极点模型12/11/2022160控制系统计算机辅助设计---MATLAB语言与应用例4-1高阶系统稳定性判定12/11/202210控制系统例4-2高阶离散单位负反馈系统模型MATLAB求解12/11/2022161控制系统计算机辅助设计---MATLAB语言与应用例4-2高阶离散单位负反馈系统模型12/11/202214.1.2线性反馈系统的内部稳定性输入、输出稳定是不够的，因为若内部信号可能过大，对系统作硬件破坏应该引入内部稳定性概念，保证内部信号也是稳定的。12/11/2022162控制系统计算机辅助设计---MATLAB语言与应用4.1.2线性反馈系统的内部稳定性输入、输出稳定是不够的由给定稳定输入到内部信号都稳定的系统称为内部稳定系统传递函数矩阵其中逐一判定每个子传递函数的稳定性很烦琐内部稳定性定理12/11/2022163控制系统计算机辅助设计---MATLAB语言与应用由给定稳定输入到内部信号12/11/2内部稳定性定理闭环系统内部稳定的充要条件为没有不稳定零点没有不稳定零极点对消第一个条件等效于输入输出稳定性判定第2条件即可可以编写MATLAB函数判定内部稳定性12/11/2022164控制系统计算机辅助设计---MATLAB语言与应用内部稳定性定理闭环系统内部稳定的充要条件为12/11/202判定的MATLAB函数12/11/2022165控制系统计算机辅助设计---MATLAB语言与应用判定的MATLAB函数12/11/202215控制系统计4.1.3线性系统的线性相似变换系统的状态方程表示称为系统实现不同状态选择下，状态方程不惟一相似变换非奇异矩阵状态变换新状态方程模型12/11/2022166控制系统计算机辅助设计---MATLAB语言与应用4.1.3线性系统的线性相似变换系统的状态方程表示称为系状态变换公式MATLAB求解方法12/11/2022167控制系统计算机辅助设计---MATLAB语言与应用状态变换公式12/11/202217控制系统计算机辅助设计-例4-3已知系统和转换矩阵MATLAB求解12/11/2022168控制系统计算机辅助设计---MATLAB语言与应用例4-3已知系统和转换矩阵12/11/202218控制系统变换结果可见，相似变换能改变系统的结构引入相似变换矩阵，可以将已知系统转换成其他的形式12/11/2022169控制系统计算机辅助设计---MATLAB语言与应用变换结果12/11/202219控制系统计算机辅助设计---4.1.4线性系统的可控性分析可控性定义

系统的可控性就是指系统内部的状态是不是可以由外部输出信号控制的性质，12/11/2022170控制系统计算机辅助设计---MATLAB语言与应用4.1.4线性系统的可控性分析可控性定义12/11/20线性系统的可控性判定可控性判定矩阵

基于MATLAB的判定方法构造可控性判定矩阵12/11/2022171控制系统计算机辅助设计---MATLAB语言与应用线性系统的可控性判定可控性判定矩阵12/11/202221控例4-4离散状态方程的可控性MATLAB求解12/11/2022172控制系统计算机辅助设计---MATLAB语言与应用例4-4离散状态方程的可控性12/11/202222控制系判定矩阵判定矩阵构造方法这样的判定方法同样适合于连续系统和离散系统。也适用于多变量模型12/11/2022173控制系统计算机辅助设计---MATLAB语言与应用判定矩阵12/11/202223控制系统计算机辅助设计---由Gram矩阵判定可控性引入可控Gram矩阵该矩阵满足Lyapunov方程MATLAB求解矩阵构造12/11/2022174控制系统计算机辅助设计---MATLAB语言与应用由Gram矩阵判定可控性引入可控Gram矩阵12/1例4-5求Gram矩阵MATLAB命令Gram矩阵12/11/2022175控制系统计算机辅助设计---MATLAB语言与应用例4-5求Gram矩阵12/11/202225控制系统可控性阶梯分解对于不完全可控的系统阶梯分解阶梯标准型MATLAB函数调用若原系统状态方程完全可控，则不必分解12/11/2022176控制系统计算机辅助设计---MATLAB语言与应用可控性阶梯分解对于不完全可控的系统阶梯分解12/11/202例4-6不完全可控系统12/11/2022177控制系统计算机辅助设计---MATLAB语言与应用例4-6不完全可控系统12/11/202227控制系统计算4.1.5线性系统的可观测性分析可观测性定义

系统的可观测性就是指系统内部的状态是不是可以由系统输出信号重建起来的性质12/11/2022178控制系统计算机辅助设计---MATLAB语言与应用4.1.5线性系统的可观测性分析可观测性定义12/11/可观测性判定判定矩阵等同于系统可控性判定Gram矩阵MATLAB求解12/11/2022179控制系统计算机辅助设计---MATLAB语言与应用可观测性判定判定矩阵12/11/202229控制系统计算机辅Gram矩阵满足Lyapunov方程对偶问题12/11/2022180控制系统计算机辅助设计---MATLAB语言与应用Gram矩阵满足Lyapunov方程12/11/2024.1.6Kalman规范分解Kalman规范分解12/11/2022181控制系统计算机辅助设计---MATLAB语言与应用4.1.6Kalman规范分解Kalman规范分解子空间示意图12/11/2022182控制系统计算机辅助设计---MATLAB语言与应用子空间12/11/202232控制系统计算机辅助设计---M4.1.6系统状态方程标准型的MATLAB求解常用标准型单变量系统的标准型MATLAB默认的标准型可控标准型实现可观测标准型实现和Jordan标准型实现多变量系统Leunberge标准型侧重点：如何用MATLAB直接获取标准型12/11/2022183控制系统计算机辅助设计---MATLAB语言与应用4.1.6系统状态方程标准型的MATLAB求解常用标准单变量系统的标准型可控标准型可观测标准型12/11/2022184控制系统计算机辅助设计---MATLAB语言与应用单变量系统的标准型可控标准型12/11/202234控制系统可控可观测标准型转换12/11/2022185控制系统计算机辅助设计---MATLAB语言与应用可控可观测标准型转换12/11/202235控制系统计算机辅可控标准型和可观测标准型，对偶关系12/11/2022186控制系统计算机辅助设计---MATLAB语言与应用12/11/202236控制系统计算机辅助设计---MATLJordan标准型

MATLAB变换12/11/2022187控制系统计算机辅助设计---MATLAB语言与应用Jordan标准型12/11/202237控制系统计算机多变量系统的Leunberge标准型由可控性判定矩阵构造矩阵12/11/2022188控制系统计算机辅助设计---MATLAB语言与应用多变量系统的Leunberge标准型由可控性判定矩阵12得出Leunberge变换矩阵编写leunberge.m函数12/11/2022189控制系统计算机辅助设计---MATLAB语言与应用得出Leunberge变换矩阵12/11/202239控MATLAB函数清单12/11/2022190控制系统计算机辅助设计---MATLAB语言与应用MATLAB函数清单12/11/202240控制系统计算机12/11/2022191控制系统计算机辅助设计---MATLAB语言与应用12/11/202241控制系统计算机辅助设计---MATL标准型的变换方法总结可控标准型可观测标准型Jordan标准型Leunberge标准型12/11/2022192控制系统计算机辅助设计---MATLAB语言与应用标准型的变换方法总结可控标准型12/11/202242控制系例4-7求解可观测标准型标准型12/11/2022193控制系统计算机辅助设计---MATLAB语言与应用例4-712/11/202243控制系统计算机辅助设计---例4-8已知模型12/11/2022194控制系统计算机辅助设计---MATLAB语言与应用例4-8已知模型12/11/202244控制系统计算机辅助4.1.7系统的范数测度及求解系统也有范数

12/11/2022195控制系统计算机辅助设计---MATLAB语言与应用4.1.7系统的范数测度及求解系统也有范数12/11/2离散系统的范数定义范数的MATLAB求解12/11/2022196控制系统计算机辅助设计---MATLAB语言与应用离散系统的范数定义12/11/202246控制系统计算机辅助例4-9已知离散系统模型12/11/2022197控制系统计算机辅助设计---MATLAB语言与应用例4-9已知离散系统模型12/11/202247控制系统计4.2线性系统时域响应解析解法给线性系统一个激励信号，输出是什么？有两大类方法解析解方法求解微分方程、差分方程解析解数值解方法主要内容基于状态方程的解析解方法基于传递函数部分方式展开的解析解方法二阶系统的解析解方法12/11/2022198控制系统计算机辅助设计---MATLAB语言与应用4.2线性系统时域响应解析解法给线性系统一个激励信号，输4.2.1基于状态方程的解析解方法状态方程模型解析解求解难点12/11/2022199控制系统计算机辅助设计---MATLAB语言与应用4.2.1基于状态方程的解析解方法状态方程模型12/11状态增广方法消除B矩阵，变成自治系统

增广状态方程自治系统可以直接求解析解12/11/2022200控制系统计算机辅助设计---MATLAB语言与应用状态增广方法消除B矩阵，变成自治系统12/11/20225一般输入信号的系统增广一般输入信号模型引入增广状态变量12/11/2022201控制系统计算机辅助设计---MATLAB语言与应用一般输入信号的系统增广一般输入信号模型12/11/20225增广状态方程模型其中解析解12/11/2022202控制系统计算机辅助设计---MATLAB语言与应用增广状态方程模型12/11/202252控制系统计算机辅助设MATLAB实现函数12/11/2022203控制系统计算机辅助设计---MATLAB语言与应用MATLAB实现函数12/11/202253控制系统计算机调用格式信号描述12/11/2022204控制系统计算机辅助设计---MATLAB语言与应用12/11/202254控制系统计算机辅助设计---MATL例4-10连续系统模型初值输入信号求解析解12/11/2022205控制系统计算机辅助设计---MATLAB语言与应用例4-10连续系统模型12/11/202255控制系统计算系统增广增广模型12/11/2022206控制系统计算机辅助设计---MATLAB语言与应用系统增广12/11/202256控制系统计算机辅助设计---解析解求解解析解求解结果稳定性12/11/2022207控制系统计算机辅助设计---MATLAB语言与应用解析解求解12/11/202257控制系统计算机辅助设计--4.2.2基于部分分式展开方法求解连续系统的解析解法

无重根时部分方式展开12/11/2022208控制系统计算机辅助设计---MATLAB语言与应用4.2.2基于部分分式展开方法求解连续系统的解析解法12由Laplace反变换求解析解有重根时相应项的解析解为12/11/2022209控制系统计算机辅助设计---MATLAB语言与应用由Laplace反变换求解析解12/11/202259控部分分式的MATLAB求解例4-11输入信号为阶跃信号输出信号计算12/11/2022210控制系统计算机辅助设计---MATLAB语言与应用部分分式的MATLAB求解12/11/202260控制系MATLAB求解解析解解析解精确值12/11/2022211控制系统计算机辅助设计---MATLAB语言与应用MATLAB求解12/11/202261控制系统计算机辅助例4-12带有复数极点的系统阶跃响应解析解解析解12/11/2022212控制系统计算机辅助设计---MATLAB语言与应用例4-12带有复数极点的系统12/11/202262控制解析解的进一步化简基于Euler公式的化简其中新MATLAB函数

12/11/2022213控制系统计算机辅助设计---MATLAB语言与应用解析解的进一步化简基于Euler公式的化简12/11/2新MATLAB函数清单12/11/2022214控制系统计算机辅助设计---MATLAB语言与应用新MATLAB函数清单12/11/202264控制系统计例4-13仍考虑MATLAB求解解析解12/11/2022215控制系统计算机辅助设计---MATLAB语言与应用例4-13仍考虑12/11/202265控制系统计算机辅基于Laplace变换的求解参附录A步骤：定义符号变量描述原函数表达式调用laplace()函数或ilaplace()函数求解结果化简，如simple()函数求解举例12/11/2022216控制系统计算机辅助设计---MATLAB语言与应用基于Laplace变换的求解参附录A12/11/202例1MATLAB求解解析解12/11/2022217控制系统计算机辅助设计---MATLAB语言与应用例112/11/202267控制系统计算机辅助设计---MA例2MATLAB求解解析解12/11/2022218控制系统计算机辅助设计---MATLAB语言与应用例212/11/202268控制系统计算机辅助设计---MA离散系统的解析解法Z变换无重根时部分分式展开解析解12/11/2022219控制系统计算机辅助设计---MATLAB语言与应用离散系统的解析解法Z变换12/11/202269控制系统计考虑采样周期例4-1412/11/2022220控制系统计算机辅助设计---MATLAB语言与应用考虑采样周期12/11/202270控制系统计算机辅助设计-输出信号解析解Z变换求解步骤定义符号变量调用iztrans()函数求解化简12/11