2018/12/9
8
传递函数输入举例

例3-1 输入传递函数模型

MATLAB输入语句

在MATLAB环境中建立一个变量 G
9
2018/12/9
另外一种传递函数输入方法

例 如何处理如下的传递函数？

定义算子
，再输入传递函数
2018/12/9
10


应该根据给出传递函数形式选择输入方法 例 输入混合运算的传递函数模型
2018/12/9
59

得出结果

相同位置的零极点，可以对消 问题：状态方程如何处理？ MATLAB解决方法
2018/12/9
60
例3-24 多变量模型

不能直接看出是否最小实现
61
2018/12/9

MATLAB求解
2018/12/9
62
3.5 线性系统模型降阶


用低阶模型近似高阶模型 和最小实现不同 最早由Edward J. Davison提出(1966) 主要内容
4
2018/12/9
3.1 连续线性系统的数学 模型与MATLAB表示


3.1.1线性系统的状态方程模型
3.1.2 线性系统的传递函数模型 3.1.3 线性系统的零极点模型 3.1.4 多变量系统的传递函数矩阵模型

2018/12/9
5
3.1.1 线性连续系统数学模型及 MATLAB 表示
55
2018/12/9
例 连续多变量模型

状态方程获取
2018/12/9
56

得出的状态方程模型

ioDelay矩阵
2018/12/9
57

该模型可以转换回传递函数矩阵

得出的转换结果
2018/12/9
58
3.4.4 状态方程的最小实现
例3-23 观察传递函数模型


未见有何特殊 求取零极点模型
21
2018/12/9
例3-6 多变量模型
2018/12/9
22
3.2 线性离散时间系统的数学模型

单变量系统：差分方程取代微分方程

主要内容
离散传递函数 离散状态方程

2018/12/9
23
3.2.1 离散传递函数模型

数学表示 （Z变换代替Laplace变换）
MATLAB表示 （采样周期 ）
当然由前面的公式也能直接求解
2018/12/9
89
例3-32 实测数据
2018/12/9
90

基于MATLAB的求解
2018/12/9
91

数学形式

辨识模型的提取

还可以写成
2018/12/9
92

还可以由下面语句求解

辨识结果
2018/12/9
93

直接辨识方法

辨识结果

辨识界面：ident

2018/12/9
74
3.5.3 时间延迟模型的 Padé近似

纯延迟的Padé 近似方法 近似函数 纯滞后逼近


2018/12/9
75

编写 MATLAB 函数

其中 r/m 任意选择 可以选择 0/m ，以避免非最小相位模型
76
2018/12/9
例3-24 纯延迟模型

MATLAB求解


系统数学模型的获取
建模方法：从已知的物理规律出发，用数学推 导的方式建立起系统的数学模型 辨识方法：由实验数据拟合系统的数学模型

2018/12/9
2
系统数学模型的分类
非线性 系统 模型 线性 连续 单变量
定常
时变
离散 混合
多变量
2018/12/9
3
主要内容



线性连续系统的数学模型与MATLAB表示 线性离散时间系统的数学模型 方框图描述系统的化简 系统模型的相互转换 线性系统的模型降阶 线性系统的模型辨识 本章要点简介
例3-12 原系统可以移动

新支路模型
41
2018/12/9

得出
2018/12/9
42
例3-13 电机拖动模型

2018/12/9
43

信号单独输入

得出另一个传递函数
2018/12/9
44

最终得出传递函数矩阵
2018/12/9
45
3.4 系统模型的相互转换


前面介绍的各种模型之间的相互等效变换

2018/12/9
85
3.6.1 离散系统的模型辨识

离散传递函数模型

对应的差分方程模型
2018/12/9
86

已知实测信号

输入
输出

由数据可以得出
2018/12/9
87

矩阵形式

定义残差最小指标 最小二乘解
88
2018/12/9

系统辨识工具箱求解


T 为结构体变量，T.a, T.b, tf(T)


算子输入方法：
24
2018/12/9
例3-7 离散传递函数，采样周期

MATLAB输入方法

另一种输入方法
2018/12/9
25
离散延迟系统与输入

数学模型

延迟为采样周期的整数倍 MATLAB输入方法
2018/12/9
26

MATLAB表示方法
例3-8
2018/12/9
27
3.2.2 离散状态方程模型
与Routh算法 时间延迟模型的 近似 带有延迟的最优降阶算法 状态空间的降阶算法

2018/12/9
63
3.5.1 降阶算法 与 Routh 降阶算法

原始模型

寻求降阶模型

假设
64
2018/12/9

展开原模型

其中时间矩量 可以递推求出

若已知状态方程模型
2018/12/9
65
3.1.3 线性系统的零极点模型

零极点模型是因式型传递函数模型

零点 、极点 零极点模型的 MATLAB表示
和增益
2018/12/9
19
例3-5 零极点模型

MATLAB输入方法

另一种输入方法
2018/12/9
20
3.1.4 多变量系统传递函数矩阵模型

传递函数矩阵

为第 i 输出对第 j 输入的传递函数 可以先定义子传递函数，再由矩阵定义

时间矩量的MATLAB求解

降阶思想：保留前
时间矩量
2018/12/9
66

对比系数，则
2018/12/9
67

这样可以得出
2018/12/9
68

降阶求解函数
2018/12/9
69
例3-25 原始模型

Padé近似

结果
2018/12/9
70
例3-26 反例
零极点模型求取

稳定模型
2018/12/9

线性系统的传递函数模型

为阶次，
为常数，
物理可实现
6
2018/12/9
传递函数的引入
Pierre-Simon Laplace (1749--1827)，法国数学家 Laplace变换 Laplace变换的一条重要性质： 若 则
2018/12/9
7
传递函数表示

数学方式

MATLAB输入语句
时域响应比较 频域响应比较


降阶模型的应用
仿真应用（用途越来越小） 控制器设计应用

2018/12/9
84
3.6 线性系统的模型辨识


模型辨识

由已知实测数据获得系统模型的方法 时域响应数据、频率响应数据
实测数据

主要内容
离散系统辨识方法 辨识信号生成 多变量系统辨识 离散系统在线辨识

MATLAB实现（从略） 调用格式
2018/12/9
81
例3-30 对给出的传递函数进行降阶研究

可以给出下面的语句

得出的降阶模型为
2018/12/9
82
例3-31 已知高阶模型

可以给出如下命令

得出的降阶模型
2018/12/9
83
降阶算法综述


状态方程方法不能任意选择分母分子阶次， 而很多传递函数方法可以 降阶效果比较，下章给出
71

Padé近似

不稳定降阶模型


Padé不能保证降阶模型的稳定性 不稳定降阶模型可能得出稳定降阶模型
72
2018/12/9
Routh 降阶方法与实例

Routh算法（较烦琐，从略）
2018/12/9
73
Routh算法的最大特色：稳定系统降阶后能 保证降阶模型稳定性 例3-23 仍考虑稳定模型
3.1.2 线性系统的状态方程模型

状态方程模型

状态变量 ， 阶次 n ，输入和输出 非线性函数： 一般非线性系统的状态方程描述