第二章 线性系统的数学模型
2.1 列写系统微分方程
人们常将描述系统工作状态的各物理量随时间变化的规律用数学表达式或图形表示出来,这种描述系统各个物理量之间关系的数学表达式或图形称为系统的数学模型。
建立数学模型有两种方法:机理分析法和实验辨识法。机理分析法是通过理论推导得出,这种方法是根据各环节所遵循的物理规律来编写;实验辨识法是由实验求取,即根据实验数据通过整理编写出来。 本章着重讨论机理分析法。
一、列写部件微分方程的目的、方法与步骤
目的:通过该方程确定被控量与给定量及扰动量之间的函数关系。 (1)根据实际情况,确定系统的输入、输出变量。
(2)从输入端开始,按信号传递遵循的有关规律列出元件微分方程。 (3)消去中间变量,写出输入、输出变量的微分方程。 (4)整理,输入量项=输出量项。
举例1 编写RC 电路微分方程
(1)确定输入、输出量为ui 、u0 (2)根据电路原理列微分方程
0i u Ri u =+
0du i C dt
=
(3)消去中间变量,可得电路微分方程
00i du
RC u u dt
+=
举例2
编写电枢控制的他励直流电动机的微分方程
解:(1)确定输入、输出量为ud 、n (2)根据电路原理列微分方程
d d d d d di
e i R L u dt
++=
d e e c n =
根据电动机力矩平衡原理列微分方程
2375GD dn M dt
=
m d M c i =
(3)消去中间变量,可得电路微分方程
2222
375375d d d
d m
e e L GD R u GD d n dn n R dt c c dt c ++= 令:d
d d
L T R =
2375d
m m e
R GD T c c =
则得:22d m d m e
u d n dn
T T T n dt dt c ++=
举例3 具有质量弹簧阻尼器的机械位移系统
(1)确定输入、输出量为F 、y
(2)根据力学、运动学原理列微分方程
s f ma F F F =--
22d y a dt
=
f dy
F f dt
=
(3)消去中间变量,可得电路微分方程
22d y dy m f ky F dt dt
++= 以上两例中的物理系统不尽相同,但它们的数学模型却是相同的,我们把具有相同数学模型的不同物理系统称之为相似系统。在相似系统中,占据相应位置的物理量称为相似量。
对于同一个物理系统,当输入量、输出量改变时,所求出的数学模型却是不同的。利用相似系统的概念,我们可以用一个易于实现的系统来研究与其相似的复杂系统,并根据相似系统的理论出现了仿真研究法。 二、系统动态微分方程的编写 (1)确定系统输入量、输出量;
(2)从输入端开始将系统划分为若干个元部件,依有关定理列写各个部件的方程组; (3)消去中间变量; (4)整理。
举例4 列写直流调速系统的微分方程
(1)确定输入、输出量为Ug 、n
(2)根据力电路、电动机力矩平衡原理列微分方程
1()k g f u K u u =- d s k u K u =
22d d m m e
u d n dn T T T n dt dt c ++=
f f u K n =
(3)消去中间变量得直流调速系统的动态微分方程
211(1)d m m r
g k k k e T T T K d n dn n U K dt K dt K C ++=+++
其中1r s K K K =为正向通道电压放大系数
1s f
k e
K K K K C =
为系统开环放大系数
2.2 非线性数学模型的线性化
对于部分的非线性系统来说,是在一定的条件下可近似地视作线性系统,这种有条件地把非线性系统数学模型化为线性数学模型来处理的方法,称为非线性数学模型的线性化。
这样做会使问题简化,给控制系统的研究工作带来很大的方便,是工程中一种常见的、比较有效的方法。 线性化原理:
设非线性方程为 ()y f r =,工作点为00()y f r =,其各阶导数均存在,则可在工作点附近展开成泰勒级数
0020002
()1()
()()()()()2!r r df r d f r y f r r r r r dr dr
=+-+-
+
当 0()r r -很小时,忽略二次以上导数项0
00()
()()()r df r y f r r r dr =+- 或可表达为:00()y y k r r -=-
y k r ?=?
简写为 y kr = 式中0
()
(
)r df r k dr = 这就是非线性化方程,这种线性化方法叫做小偏差方法。 注意:
1.非线性方程必为连续。
原因:断续的方程不能用台劳级数展开,因此不能采用此方法。这类非线性称为本质非线性。
2.K 值与工作点的位置有关。
3.考虑增量ΔX 较小。 2.3 传递函数
求解微分方程可求出系统的输出响应,但如果方程阶次较高,则计算非常繁琐,因此对系统的设计分析不便,所以应用传递函数将实数中的微分运算变成复数中的代数运算,可使问题分析大大简化。 一、传递函数的概念及意义 (1)传递函数的定义:
线性系统在零初始条件下,输出信号的拉氏变换与输入信号的拉氏变换之比。 线性定常系统微分方程的一般表达式:
110110n n m c c c r
n n m
n n c m r d x d x dx d x dt dt dt dt a a a a x b b x ---++
++=+
+
1010()()()()()n n m c c n c r m r a s x s a s x s a x s b s x s b x s -+++=++
其中c x 为系统输出量,r x 为系统输入量
在初始情况为零时,两端取拉氏变换:
1010()()()()()n n m c c n c r m r a s x s a s x s a x s b s x s b x s -+++=++ 移项后得:
()()/()C r W s X s X s = 11011011/()m m n n m m n n b s b s b s b a s a s a s a ++--=++++++++
上式中Xc(s)输出量的拉氏变换;Xr(s)输入量的 拉氏变换; W(s) 为系统或环节的传递系数。
(2)传递函数的两种表达形式 a.传递函数的零极点表示形式
1110
11110
1
()
()()
m
m m i
m m i g
n
n n n n j
j s z b s d s d s d G s K a s c s c s c s p --=--=+++++=
=++
+++∏∏
b.传递函数的时间常数表示形式
1
011110111
(1)
1()1
(1)
m
m m i
m m i n
n n n n i
j s b f s f s f s G s K
a e s e s e s T s τ--=--=+++++==+++++∏∏
(3)关于传递函数的几点说明
a.传递函数的概念只适应于线性定常系统。
b.传递函数只与系统本身的特性参数有关,而与输入量变化无关。
c.传递函数不能反映非零初始条件下系统的运动规律。
d.传递函数分子多项式阶次低于或至多等于分母多项式的阶次。
二、典型环节的传递函数及其暂态特性
无论什么样的系统,它的传递函数都是一些基本因子相乘积而得到的。这些基本因子就是典型环节对应的传递函数。把复杂的物理系统划分为若干个典型环节,利用传递函数和框图来进行研究,这是研究系统的一种重要方法。
(1)比例环节(放大环节/无惯性环节)
特点:输入量与输出量的关系为一种固定的比例关系。
(2)惯性环节
特点:只包含一个储能元件,使其输出量不能立即跟随输入量的变化,存在时间上的延迟。
(3)积分环节
特点:输出量随时间成正比地无限增加。
(4)振荡环节
特点:振荡的程度与阻尼系数有关。
(5)微分环节
特点:是积分环节的逆运算,其输出量反映了输入信号的变化趁势。实践中,理想的微分环节难以实现。
(6)延迟环节(时滞环节、滞后环节)
特点:输出信号经过一段延迟时间τ后,可完全复现输入信号。
2.4 系统动态结构图
一、概念
1.动态结构图:是描述系统各组成元件之间信号传递关系的数学图形,它表示了系统的输入输出之间的关系。
2.结构图的组成:
(1)信号线:带箭头的直线,箭头表示信号传递方向。
(2)引出点(分离点):表示信号引出或测量的位置。
(3)比较点(相加点):对两个以上信号加减运算。
(4)方框:方框图内输入环节的传递函数。
3.动态结构图的绘制步骤:
(1)确定系统输入量与输出量。
(2)将复杂系统划分为若干个典型环节。
(3)求出各典型环节对应的传递函数。
(4)作出相应的结构图。
(5)按系统各变量的传递顺序,依次将各元件的结构图连接起来。
二、结构图的简化法则
常用的结构图变换方法可归纳为两类:一类是环节的合并,另一类是信号的分支点或相加点的移动。
结构图的变换必须遵循的原则是:变换前后的数学关系保持不变,因而也称为结构图的等效变换。
(一)环节的合并
法则一环节串联,传递函数相乘。
法则二环节并联,传递函数相加。
法则三反馈连接的等效传递函数。
法则四相加点从环节输入端移到输出端
法则五相加点从环节输出端移到输入端
法则六分支点从环节输入端移到输出端
法则七分支点从环节输出端移到输入端
法则八两个分支点、相加点间可以相互换位。相加点和分支点之间一般不能换位。
三、开环传递函数
1.开环传递函数:是闭环系统反馈信号的拉氏变换与偏差信号拉氏变换之比。
2.开环传递函数的求法
(1)单回路系统
(2)多回路系统
a.无交错局部反馈
结论: b.有交错局部反馈
1232312()1()
k
G G G H
G s G G H H =++
结论:
四、闭环传递函数 1.闭环传递函数:在初始条件为零的情况下,系统输出量的拉氏变换与输入量的拉氏变换之比。
2.闭环传递函数的求法
典型结构的闭环控制系统中有两个输入量——给定输入和扰动输入,它们同时作用于系统。对于线性系统而言,可通过分别讨论在不同输入信号作用下的系统传递函数和相应的输出,再通过叠加得到总的输出。 (1)给定输入单独作用下的闭环系统
(2)扰动输入单独作用下的闭环系统
3.误差传递函数:误差信号的拉氏变换与输入信号的拉氏变换之比。 (1)给定输入单独作用下的闭环系统
(2)扰动输入单独作用下的闭环系统
12343
232341
()1k G G G G H G s G G H G G H =
++1212
12()11k
G G G G
s G G H G Φ=
=++22
12()11k
G G s G G H G Φ=
=
++1211
()11r k
E s G G H G ==
++
4.给定输入和扰动输入作用下的闭环系统的总的输出量和偏差输出量 由线性系统叠加性可知:
2.5 信号流程图
一、基本概念及常用术语
信号流程图简称信号流图,是一种用图线表示线性方程组的方法。 设有线性方程组
11,2,n
i ij j
i x a x i n ===∑
经拉氏变换将微分方程转化成为线性代数方程
1
()()()
1,2,n
i ij j j X s G s X s i n ===∑
可以用信号流图来表示以上方程组各变量之间的关系。 信号流图的绘制方法:
设系统的描述方程为 12x ax =,其中1x 为输入变量,2x 为输出变量,a 为两个变量
间的传输。该系统可以用信号流图表示为:
常用术语:
1.节点 表示信号/变量。
(1)源节点 只有输出量的节点。 (2)汇节点 只有输入量的节点。
(3)混合节点 既有输入又有输出的节点。 2.支路 连接两节点间的有向线段。 3.支路增益 表示信号间的因果关系。
4.通路 从一个节点到另一节点的路径,其间每个节点只通过一次。 (1)前向通路 从源点到汇点 (2)闭通路 起点与终点为同一点
5.通路增益 通路中所有支路增益之积
6.不接触回路 回路之间没有公共节点 二、如何作信号流图
2212()11n
k G H G H E s G G H G ==++()()()()()()()()()()
n r n Y s s R s s N s E s E s R s E s N s =Φ+Φ=+
方法一:将系统微分方程作拉氏变换后,按所得代数方程作图。 例1 绘制二级RC 滤波电路的信号流图。
解(1)列写系统微分方程组 1311
u u i R -=
213i i i =-
331
1
u i dt c =?
3232u u i R -=
232
1
u i dt c =?
(2)对上述微分方程作拉氏变换
1131
1
()[()()]I s U s U s R =-
213()()()I s I s I s =-
3211
()()U s I s c s
=
3322
1
()[()()]I s U s U s R =-
2321
()()U s I s c s
=
(3)对上述代数方程作信号流图
(4)综合作出系统信号流图
方法二:由系统动态结构图变形得来。
变形原则:原信号线变为节点,传递函数变为支路增益。
三、由梅逊公式求传递函数
梅逊公式:1
1n
k k k T T ==??∑
T ——闭环传递函数 Δ——特征式
1231(1)m m L L L L ?=-+-++-∑∑∑∑
其中
1L ∑——为所有不同回路的增益之和
2
L ∑——每两个互不接触回路增益乘积之和 3
L ∑——每三个互不接触回路增益乘积之和 m
L ∑——每m 个互不接触回路增益乘积之和
n ——前向通路的条数
k T ——第K 条前向通路的增益
k ?——第K 条前向通路的余因子,在Δ中除去与第K 条前向通路相接触的回路增益后剩余的Δ值。
例2 用梅逊公式求下图中信号流图的传递函数。
解:(1)找出上图中所有的前向通路只有一条前向通路
11234T G G G G =
(2)找出系统中存在的所有的回路共有三个回路,三个回路的传输之和为
23634512347G G G G G G GG G G G ?=---
(3)这三个回路都存在公共节点,即不存在不接触回路。故系统的特征方程式为:
11L ?=-∑
236345123471G G G G G G GG G G G =+++
(4)由于这三个回路都与前向通路相接触,故其余因子Δ1=1。 (5)故该系统的传递函数为:
11()
()T Y s T R s ?==?
1234
236345123471G G G G G G G G G G G G G G G =
+++
2.6 脉冲响应函数
单位脉冲信号:
分析:()()r t t δ=其拉氏变换R(s)=1。
由 ()
()()
Y s s R s φ= 可知()()()()Y s s R s s φφ==
结论:系统或环节的单位脉冲响应函数的拉氏变换即为系统或环节的传递函数。
)Q c
r