国家级精品课程
控制系统仿真与CAD
第五章 线性控制系统的计算机辅助分析
东北大学信息学院 薛定宇
控制系统仿真与CAD 国家级精品课程
2014-12-31
1/139
本章主要内容
线性系统定性分析 线性系统时域响应解析解法 线性系统的数字仿真分析 根轨迹分析 线性系统频域分析 多变量系统的频域分析
如果系统中所有的状态都是可控的，则称该系 统为完全可控的系统。 系统的可控性就是指系统内部的状态是不是可 以由外部输出信号控制的性质
控制系统仿真与CAD 国家级精品课程
2014-12-31


20/139
线性系统的可控性判定

可控性判定矩阵


若矩阵 为满秩矩阵，则系统完全可控 基于 MATLAB 的判定方法
控制系统仿真与CAD 国家级精品课程
2014-12-31
43/139

离散系统的范数定义

范数的 MATLAB 求解
控制系统仿真与CAD 国家级精品课程
2014-12-31
44/139
例5-9 已知离散系统模型
控制系统仿真与CAD 国家级精品课程
2014-12-31
45/139
5-1 系统性质分析小结
控制系统仿真与CAD 国家级精品课程
2014-12-31
14/139

判定的 MATLAB 函数

内部稳定返回0，内部不稳定但输入输出稳定返 回1，否则返回2
控制系统仿真与CAD 国家级精品课程
2014-12-31
15/139
5.1.3 线性系统的线性相似变换



系统的状态方程表示称为系统实现 不同状态选择下，状态方程不唯一 相似变换

若原系统状态方程完全可控，则不必分解
控制系统仿真与CAD 国家级精品课程
2014-12-31
25/139
例5-6 不完全可控系统
控制系统仿真与CAD 国家级精品课程
2014-12-31
26/139
4.1.5 线性系统的可观测性分析

可观测性定义
假设系统由状态方程 给出，对 任意的初始时刻 ，如果状态空间中任一状 态 在任意有限时刻 的状态 可以由 输出信号在这一时间区间内 的值精 确地确定出来，则称此状态是可观测的。 如果系统中所有的状态都是可观测的，则称 该系统为完全可观测的系统。 系统的可观测性就是指系统内部的状态是不 是可以由系统输出信号重建起来的性质

控制系统仿真与CAD 国家级精品课程
2014-12-31
4/139
5.1.1 线性系统的稳定性分析

给定线性系统模型，如何分析稳定性？
单位负反馈闭环系统 由控制理论可知，用Routh 表格可以判定该系统稳定性。 Edward John Routh (1831-1907) 历史局限性

控制系统仿真与CAD 国家级精品课程
2014-12-31
49/139
5.2.2 基于状态方程的解析解方法

状态方程模型 解析解


求解难点 想法：如果能用某种方法消去积分项最好
控制系统仿真与CAD 国家级精品课程
2014-12-31

50/139
状态增广方法

消除B 矩阵，变成自治系统


增广状态方程

常用标准型

单变量系统的标准型
MATLAB 默认的标准型 可控标准型实现 可观测标准型实现 和 Jordan 标准型实现


多变量系统 Luenberger 标准型

侧重点：如何用 MATLAB 直接获取标准型
控制系统仿真与CAD 国家级精品课程
2014-12-31
31/139
单变量系统的标准型

传递函数的一般形式

处理后的传递函数标准型

其中 由前面给出的传递函数标准型可以直接写 出系统的可控和可观测标准型
控制系统仿真与CAD 国家级精品课程
2014-12-31
32/139

传递函数一般形式

可控标准型
控制系统仿真与CAD 国家级精品课程
2014-12-31
33/139

传递函数一般形式
离散系统：pzmap(G)，同时画出单位圆
控制系统仿真与CAD 国家级精品课程
2014-12-31
9/139
例5-1 高阶系统稳定性判定

直接分析方法

零极点模型
控制系统仿真与CAD 国家级精品课程
2014-12-31
10/139
例5-2 高阶离散单位负反馈系统模型

MATLAB 求解
控制系统仿真与CAD 国家级精品课程
控制系统仿真与CAD 国家级精品课程
2014-12-31
46/139
5.2 线性系统时域响应解析解法


给线性系统一个激励信号，输出是什么？ 有两大类方法

解析解方法

求解微分方程、差分方程解析解

数值解方法 基于状态方程的解析解方法 基于传递函数的解析解方法 二阶系统的解析解方法
控制系统仿真与CAD 国家级精品课程

标准4-12-31
41/139
例5-8 已知模型，求Luenberger标准型

MATLAB处理
控制系统仿真与CAD 国家级精品课程
2014-12-31
42/139
5.1.8 系统的范数测度及求解


系统也有范数：范数即测度 范数

范数

范数是系统频域响应的峰值
控制系统仿真与CAD 国家级精品课程
2014-12-31
18/139

变换结果

可见，相似变换能改变系统的结构 引入相似变换矩阵，可以将已知系统转换 成其他的形式
控制系统仿真与CAD 国家级精品课程
2014-12-31
19/139
4.1.4 线性系统的可控性分析

可控性定义

假设系统由状态方程 给出，对任 意的初始时刻 ，如果状态中任一状态 可 以从初始状态 处，由有界的输入信号 的驱动下，在有限时间 内能够到达任意预先 指定的状态 ，则称此状态是可控的。
控制系统仿真与CAD 国家级精品课程
2014-12-31



8/139
基于 MATLAB 的稳定性判定方法

直接判定
状态方程模型、传递函数模型等
eig(G) 可以求出所有特征根 离散系统：abs(eig(G))
由

isstable(G) 也可以用于系统的稳定性判定

图解判定法
连续系统：pzmap(G)
非奇异矩阵 状态变换 新状态方程模型

且
控制系统仿真与CAD 国家级精品课程
2014-12-31
16/139

状态变换公式

MATLAB 直接求解方法

也可以根据上面公式进行变换
控制系统仿真与CAD 国家级精品课程
2014-12-31
17/139
例5-3 已知系统和转换矩阵

MATLAB 求解
2014-12-31
11/139
5.1.2 线性反馈系统的内部稳定性

输入、输出稳定是不够的，因为若内部信 号可能过大，对系统作硬件破坏
应该引入内部稳定性概念，保证内部信号 也是稳定的。

控制系统仿真与CAD 国家级精品课程
2014-12-31
12/139


由给定稳定输入 到内部信号 都稳定的系统称为内部稳定系统 传递函数矩阵

自治系统

可观测标准型

可控标准型和可观测标准型对偶 可以用转置表示对偶
控制系统仿真与CAD 国家级精品课程
2014-12-31

34/139

可控可观测标准型转换
控制系统仿真与CAD 国家级精品课程
2014-12-31
35/139
Jordan 标准型



MATLAB 变换
控制系统仿真与CAD 国家级精品课程

介绍了系统定性、定量分析的方法
统一框架下分析系统的稳定性——直接方法，并 介绍了内部稳定性的判定，eig/pzmap/isstable 线性系统的相似变换方法 ss2ss 系统的可控性、可观测性分析及标准型 rank 系统的Kalman分解 系统的范数计算 norm


有些内容理论上较深，如果不能理解，知道 如何用MATLAB求解就足够了
2014-12-31

主要内容


47/139
5.2.1 直接积分方法

状态方程的解析解 直接积分语句



得到结果后有必要用simple()函数化简结果 若只需状态变量，则不用乘C
控制系统仿真与CAD 国家级精品课程
2014-12-31
48/139

例5-10 系统的状态方程为

状态变量初值 输入信号

其中 逐一判定每个子传递函数的稳定性很烦琐 内部稳定性定理
控制系统仿真与CAD 国家级精品课程
2014-12-31
13/139
内部稳定性定理

闭环系统内部稳定的充要条件为

没有不稳定零点 没有不稳定零极点对消


第一个条件等效于输入输出稳定性 判定第2条件即可 可以编写MATLAB函数判定内部稳定性