第二章 数学模型微分方程
线性系统的数学模型

描述控制系统输入、输出变量及内部变量之间关 系的数学表达式称为系统的数学模型。
★ 描述控制系统的输入-输出变量数学模型:
微分方程、传递函数、方框图、频率特性
★ 描述控制系统的内部变量数学模型: 状态空间
说明 ◆ 要分析自动控制系统的性能,必须先建立该系统 的数学模型; ◆ 一个物理系统,要处理的问题或要达到的精度不 同,得到的数学模型也不同。
3.反馈
R(S) E(S) + B(S) H(S) C(S)
G(S)
负反馈 正反馈 单位反馈:H(S)=1
主 要 内 容
§2-1 微分方程 §2-2 传递函数
§2-3 典型环节的传递函数及动态响应
§2-4 电气网络的运算阻抗与传递函数 §2-5 方框图 §2-5 反馈控制系统的传递函数
§2-1
微分方程
对于线性定常系统, 可以用线性常系数微分方程 作为其数学模型,如 a 0dnc (t)/dtn +a1dn-1c (t) /dtn-1+…+anc (t) =b0dmr(t)/dtm +b1dm-1r(t)/dtm-1+…+bmr(t) c(t): 系统的输出; r(t): 系统的输入; a0……an ; b0……bm 均为实数,均由系统本身的结
对电气网络,可以不列微分方程,仅利用运算电 路,经过简单的代数运算,就可以求得传递函数!
§2-5 控制系统的方框图
方框图是以图形表示系统的数学模型;
通过方框图,能够非常清楚地表示出信号在系统各
环节之间的传递过程;
方框图可以方便地求出复杂系统的传递函数; 方框图是分析控制系统的一个简明而有效的工具。
八.二阶振荡环节 1、传递函数
自动控制原理(第三版)第2章控制系统的数学模型(2)

自动控制原理
第二章 控制系统的数学模型
求取该电路在单位阶跃输入时的响应。 U c ( s) 1 G( s ) T RC U r ( s ) Ts 1
ur 1( t )
方法1
U c ( s ) G( s )U r ( s )
1
U r (s)
1 s
方法2
1 (Ts 1) s
1 t 1 g (t ) 1[G ( s)] e T T t uc (t ) g (t )ur ( )d
0 1 1 ( t ) t t 1 T 1 T e d e e T d 0T 0 T t
1 uc (t ) L [ ] (Ts 1) s T 1 1 1 L ( )L ( ) s Ts 1 1 e
大连民族学院机电信息工程学院
自动控制原理
第二章 控制系统的数学模型
传递函数的求法
例2-1 方法一 R-L-C串联电路
d 2 uc ( t ) R duc ( t ) 1 1 uc ( t ) ur ( t ) 2 dt L dt LC LC传递Fra bibliotek数: G( s)
U c ( s) 1 U r ( s) LCs 2 RCs 1
大连民族学院机电信息工程学院
自动控制原理
第二章 控制系统的数学模型
零、极点分布图
传递函数的零、极点分 布图: 将传递函数的零、 极点表示在复平面上的 图形。
零点用“o”表示 极点用“×”表示
j
1 -3 -2
-1
s2 G( s) = ( s 3)( s 2 2s 2)
大连民族学院机电信息工程学院
控制工程基础_第二章(2017)

时,
R F (s) s
18
例 求单位斜坡函数f(t)=t的拉氏变换。 f (t )
单位斜坡函数如图(b) 所示,定义为
0 t 0 f (t ) t t 0
解:利用定义式,可得
O
t
(b)单位斜坡函数
F (s)
0
1 1 st 1 1 st 1 st t e dt t ( e ) e dt 0 e 2 0 0 s s s s 0 s
12
二.举例
1.机械系统的微分方程式
机械系统设备大致分两类:平移的和旋转的。它们之间的区 别在于前者施加的力而产生的是位移,而后者施加的是扭矩产生 的是转角。
牛顿定律和虎克定律等物理定律是建立机械系统数学模型的基础
c1 m c2 xo xi
例1(1)如图所示机械系统。求其微分方程,图中Xi 表示输入位移,Xo 表示输出位移,假设输出端无负 载效应。(c、c1、c2为阻尼系数,k1、k2为弹性系数) 由牛顿定律有: 化为标准式得:
st
例 求单位脉冲函数的拉氏变换。 单位脉冲函数如图(c)所示。定义为
0 t 0 且 (t ) t 0
0
f (t )
(t )
O
0
(t )dt 1
0
t
F ( s) (t )e st dt (t )e st dt (t )e st dt f (0) e st
图c
14
(4)机械旋转系统 图中所示转动惯量为J的转子与弹性系数为k的弹性轴和阻尼 系数为B的阻尼器连接。假设外部施加扭矩m(t),则系统产生一个 偏离平衡位置的角位移(t) 。研究外扭矩m(t)和角位移(t)的关系。
第二章控制系统数学模型

有极限存在。
终值定理在分析研究系统的稳态性能时(例如分析系统的稳态误差,求取系统
输出量的稳态值等)有着很多的应用。因此终值定理也是一个经常用到的运算
定理。
7.初值定理: lim f (t) lim sF (s)
18
2
例2-1:写出RLC串联电路的微分方程。
ui
L
R
i
C
uo
ui 输入
uo 输出
[解]:据基尔霍夫电路定理:
L di dt
Ri
1 C
idt
ui
①
uo
1 C
idt
②
由②: i C d,uo代入①得: dt
LC
d 2uo dt 2
RC
duo dt
uo
ui
这是一个线性定常二阶微分方程。
3
例2-2 设一弹簧、质量块、阻尼器组成的系统如图所示,当外力 F(t)作用于系统时,系统将产生运动。试写出外力F(t)与质量块的 位移y(t)之间的微分方程。
uR uc Us
把 uR i R
和
ic
C
duc dt
代入电路,可得到电路的
微分方程:
RC
duc dt
uc
Us
23
现在对于上面的微分方程,我们用Laplace变换求解。
首先,利用Laplace变换中的微分定理,将微分方程变换成如下形式:
RC
duc dt
uc
Us
RCsU c (s) Uc (s) Us R(s)
利用待定系数法可求得:
A 1 ARC B 0
F (s) L[ f (t)] f (t)e st dt 0
第二章 控制系统的数学模型

= Ur (s)
传递函数为: di + u ur= R · + L i c dt Uc (s) 1 = duc G (s) = i = C dt Ur (s) LCs2 + RCs + 1
电气系统三要素:电阻、电容、电感
+ ί(t) R –
u(t)= ί(t)· R
u (t )
ί(t) C
–
u(t) ί(t)= R
图2-9 速度控制系统
+
R1 R2 R2 R1 k2
ui
R1
k1 u 1
c
u2
功 ua 放
m
SM
ω
负 载
ut
TG
运算放大器
uu+ ii+
_ +
+
Add
uo
差模输入电压等于零
u+= u-
运放同相输入端与反向输入端两点的电压相等,如同该 两点短路一样,称为虚短。
i+=i-=0
运放同相输入端与反向输入端的电流都等于零,如同该 两点被断开一样,称为虚断。
Tm s m ( s ) m (t ) K1U a ( s )
Tm s 1 m ( s) K1U a ( s)
m ( s) K1 G ( s) U a ( s) Tm s 1
m ( s) K2 G ( s) M c ( s) Tm s 1
传递函数的性质(续)
(5)传递函数与微分方程有相通性;
b1s b2 C (s) G ( s) R( s ) a0 s 2 a1s a2
对角线相乘
a0 s 2 a1s a2 C ( s ) b1s b2 R ( s )
自动控制原理-第二章 控制系统的数学模型

t
f (t)dt 0
t
f ( )d
n
ki .L[ f (t )]
i 1
sF (s) f (0 )
s2F (s) sf (0 ) f (0 )
snF (s) sn1 f (0 ) sn2 f (0 ) f (n1) (0 )
电枢回路方程为
La
dia (t) dt
Raia (t)
Ea (t)
ua (t)
电磁转矩方程 M m Cmia (t)
电动机轴上转矩平衡方程
Jm
dm (t)
dt
fmm (t)
Mm
MC
(t)
若以角速度 m 为输出量、电枢电压 ua 为输入量,
消去中间变量,直流电动机的微分方程为
(s2+s+1)Uc(s)= Ur(s)+0.1(s+2)
即 U S 1 U S 0.1S 2
C
S2 S 1 r
S2 S 1
通电瞬间, ur(t)=1 或 Ur(s)=L[ur(t)]=1/S
故 U S 1 1 0.1S 2
C
S2 S 1 S S2 S 1
再对上式两边求反拉氏变换:
u c
t
L1 U C
S
L1
S
2
1 S
1
1 S
S
2
1 S
1
=1+1.15e-0.5tSin(0.866t-120°)+ 0.2e-0.5tSin(0.866t+30°)
第2章 2-1物理系统的微分方程

共计58页
20
此外,在工程实践中,控制系统都有一 个额定的工作状态和工作点,当变量在 工作点附近作小范围的变化,且变量在给 定的区域间有各阶导数时,便可在给定 工作点的邻域将非线性函数展开为泰勒 级数,忽略级数中高阶无穷小项后,就 可得到只包含偏差的一次项的线性方程。 这种线性化方法称为小偏差法。
共计58页 21
例如,设非线性函数y=f(x)如图所示, 其输入量为x,输出量为y,如果在给定工 作点y0=f(x0)处各阶导数均存在,则在 y0=f(x0)附近将y展开成泰勒级数:
共计58页
22
如果偏差Δx=x-x0很小,则可忽略级数中高阶无穷小项, 上式可写为
K表示y=f(x)曲线在(x0,y0)处切线的斜率。因此非线性函数在工作点处可以用该点 的切线方程线性化。
N1
2
(2 - 9)
[ f f ( N 1 ) 2 ] ( N 1 )T T [J1 J 2 ( ) ] 1 1 2 1 L M N2 N2 N2 N2 N 2 N 2 2 ) J 1 J 2 ]2 [ f 1 ( 2 ) f 2 ]2 T L ( 2 ) T M N1 N1 N1
共计58页
24
附录: 用拉普拉斯变换求解线性微分方程
建立了系统的微分方程以后,对微分方程求解 就可以得到表示系统动态性能的时间响应。微分方 程的求解可以用经典方法或借助于计算机进行,也 可以采用拉普拉斯变换法。 一、拉普拉斯变换定义:设有函数f (t), t为实变量, s=ζ+jω为复变量。 如果线性积分
共计58页 1
第二章 控制系统的数学模型
§2.1物理系统的微分方程 §2.2传递函数 §2.3系统框图及其简化 §2.4典型信号分析 §2.5信号流程图 §2.6线性系统的状态方程
02 数学模型 - 03非线性微分方程的线性化

第二章控制系统的数学模型第3讲非线性微分方程的线性化王燕舞为什么要进行线性化?严格的说,几乎所有元件或系统的运动方程都是非线性方程,即输入、输出和扰动等之间的关系都是非线性的。
非线性微分方程的求解和控制系统性能研究非常复杂,而线性化后的模型可借助叠加原理的性质,简化系统分析。
因此,研究非线性微分方程的线性化具有较强的工程实用价值。
什么是非线性数学模型的线性化?在一定的条件下或在一定范围内把非线性的数学模型化为线性模型的处理方法。
符合什么条件的系统可以进行线性化呢?▪条件1: 小偏差理论或小信号理论。
在工程实践中,控制系统都有一个额定的工作状态和工作点,当变量在工作点附近作小范围的变化时,就满足这个条件。
▪条件2: 在工作点附近存在各阶导数或偏导数。
如何进行线性化呢?假设微分方程模型中包含非线性函数f(x)如图所示。
设y=f(x),假设系统在工作点(x 0, y 0), y 0=f(x 0) 附近变化,且在该工作点处各阶导数均存在,在(x 0, y 0)附近将y 展开成泰勒级数:)()()()(000xx x x x f x f x f y -⎥⎦⎤⎢⎣⎡∂∂+==若偏差Δx=x-x 0很小,可忽略级数中高阶无穷小项,上式化为)()()()()()(00000x xK x f x x x x x f x f x f y -+=-⎥⎦⎤⎢⎣⎡∂∂+≈=K 表示y=f(x)曲线在(x 0,y 0)处切线的斜率。
因此非线性函数在工作点处可以用该点的切线方程线性化。
yy=f(x)y 0x 0x ⋯+-⎥⎦⎤⎢⎣⎡∂∂+20022)()(!21x x x x x f xK x f x f y y y ∆=-=-=∆)()(00如何进行线性化呢?小偏差法:在给定工作点的邻域将非线性函数展开为泰勒级数,忽略级数中的高阶项,得到只包含偏差的一次项的线性方程。
液位流体过程。
如图,Q1为流入量,也是输入量;Q 2为流出量;h 为液位高度,为系统输出;C 为液缸的截面积。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解 列方程
duc C i( t ) dt L di Ri u u( t ) c dt
2、S域模型(传递函数)
包括:系统传递函数G(s)(经典控制论的核心模型)、 系统结构图(方块图)、信号流图
3 .信号流图
图2-1 S域模型及信号流图
4 . 频率特性(幅频特性,相频特性)也是一种数学 模型。
i
dQ 1 , u idt dt c
6.电位器
7 发电机
按长度计算,或者按旋转角度计算
e k f (t )
8. 电动机
+
La
+
Ra
-
if
m
SM
ua
-
ia
- Ea +
负 J ,f m m 载
图2-3 电枢控制直流电动机原理图 电枢控制直流电动机将输入的电能转换为机械能,拖动负 载运动: ua (t ) ia (t ) M m (t ) 直流电动机运动方程满足三个方程式: 1.)电枢回路电压平衡方程: ua (t ) La
x
(思考题:上述方程中为何没有受重力mg的影响?)
控制系统的数学模型
例2-3:速度控制系统
R2
微分放 大器
功放
ug
R2
R1
K1 +
R1
u1
C
K2 +
u2
加法器
K3 +
u3 SM
m
Mc
负 载
设定部件
uf
TG
传感器
分析步骤:
⑴该系统的组成和原理; ⑵该系统的输出量是 ,输入量是 ug ,扰动量是 M c ⑶速度控制系统方块图:
四. 数学模型的简化性和精确性
一个物理系统的数学模型不是唯一的。数学模型是在一 定的条件下,根据具体问题的要求建立的,是有一定针对性和 局限性的,都是有一定近似性的。因此,使用数学模型要注意 其前提条件和局限性。要合理选择和应用不同的模型。
例 电阻模型
图2-2 电阻模型
2.2 系统的微分方程的建立方法
第二章 自动控制系统的数学模型
基本要求
1.了解建立系统动态微分方程的一般方法。 2.熟悉拉氏变换的基本法则及典型函数的拉氏变换形式。 3.掌握用拉氏变换求解微分方程的方法。 4.掌握传递函数的概念及性质。 5.掌握典型环节的传递函数形式。 6.掌握由系统微分方程组建立动态结构图的方法。 7.掌握用动态结构图等效变换求传递函数和用梅森公式求 传递函数的方法。 8.掌握系统的开环传递函数、闭环传递函数,对参考输入 和对干扰的系统闭环传递函数及误差传递函数的概念。
kf
i 测速发电机: uf K f
⑸消去中间变量:推出 ~ ug ( Mc )之间的关系:
Kg u M c g Kg ug Kc Tm
显然,转速 既与输入量
有关,也与干扰 M有关。 ug c
左 右
上式可以用与研究在给定电压或有 负载扰动转矩时,速度控制系统的 动态性能
三 .数学建模方法 1 分析法(解析法(机理分析法))——对系统各部分 的 运动机理进行分析,根据它们所依据的物理规律或化学规 律分别写相应的运动方程(由工作模型导出数学模型), 即变量间的数学表达式,并实验验证。 2 实验法(系统辨识方法)——对系统加入适当的测试 (激励)信号,记录输出响应,再用适当的数学关系模拟 其特性。实验方法适用于复杂、非常见的系统。
例如: 微分方程
L
R C
d 2uc(t) duc(t) LC RC uc(t) ui 2 dt dt
例如状态方程 例 RLC电路 输入: u(t) 输出: uc(t)
u(t )
i(t )
uc ( t )
状态变量 : uc ( t )、i ( t )
改写为
1 duc i( t ) dt C di 1 u R i 1 u c L L L dt
控制系统的微分方程
三 微分方程举例 [例2-2] 求弹簧-质量-阻尼器的机械位移系统的微分方程。设输 入量为外力F,输出量为位移x。
F m
k
F
kx
m
阻尼器是一种产生粘性摩擦的装置,由活塞和 充满油液的缸体组成。活塞和缸体之间的任何 相对运动都将受到油液的阻滞。阻尼器用来吸 收系统的能量并转变为热量而散失掉。
四 非线性微分方程的线性化
在实际工程中,构成系统的元件都具有不同程 度的非线性,如下图所示。
于是,建立的动态方程就是非线性微分方程,对其求解 有诸多困难,因此,对非线性问题做线性化处理确有必要。
对弱非线性的线性化
如上图(a),当输入信号很小时,忽略非线性影响, 近似为放大特性。对(b)和(c),当死区或间隙很小时 (相对于输入信号)同样忽略其影响,也近似为放大特性, 如图中虚线所示。
[解]:图1和图2分别为系统原理结构图 f 和质量块受力分析图。图中,m为质量, m x fx f为粘滞阻尼系数,k为弹性系数。 图1 图2 根据牛顿定理,可列出质量块的力平衡方程如下: d 2x dx d 2x dx m 2 F f kx 整理得 m 2 f kx F dt dt dt dt 这也是一个二阶定常微分方程。X为输出量,F为输入量。 在国际单位制中,m、f和k的单位分别为: kg、N s / m、N / m
输出量项=输入量项
相似系统和相似量
[需要讨论的几个问题]: 1、相似系统和相似量: 我们注意到RLC无源网络和弹簧质量阻尼机械系统的微分方程 形式是完全 一样的。 这是因为:若令 q idt (电荷),则例2-1①式的结果变为:
d 2q dq 1 L 2 R q ui dt dt C
T1T2
d 2uc duc T 2 uc ui
2.2 系统数学模型的建立方法
二 、单个部件的数学模型(输入变量和输出变量之 间的关系)
1.比较器 2.放大器 3.电阻 4.电感
5.电容
y=v1-v2 y=k *v v=R*i di
eL
dt
u
1 c C
idt
②
duc d 2uc duc LC RC uc ui 由②:i C , 代入①得: 2 dt dt dt
2 d uc ( t ) R duc ( t ) 1 1 这是一个线性定常二阶微分方程。 u ( t ) ur ( t ) c dt 2 L dt LC LC L 令: T1 ,T2 RC ,T1、T2的量纲为时间,称为电 路的时间常数。 R
一 、 微分方程的建立的一般方法步骤:
1.分析各元件的工作原理,确定系统输入变量、 输出变量; 2. 根据基本定律列原始方程; 3. 消去中间变量, 得输入、输出微分方程; 4. 微分方程标准化 ,均按微分级阶次从大到小 排列。 左边(输出变量)=右边(输入变量)
微分方程的一般形式
d n c( t ) d n 1c ( t ) dc( t ) nm an a ... a a c ( t ) n 1 1 0 dt n dt n1 dt d m r (t ) d m 1 r ( t ) dr ( t ) bm bm 1 ... b1 b0 r ( t ) m m 1 dt dt dt
式中 r(t):输入量 c(t):输出量 ai ,bj :常量 m、n :输入量、输出量导数的最高阶次
控制系统的微分方程
系统的数学模型的建立方法举例 [例1]: 写出RLC串联电路的微分方程。 L R [解]:据基尔霍夫电路定理: di 1 i ( t ) u(t ) L Ri idt ui ① uc ( t ) C dt C
K1 Cm ( Ra f m Cm Ce ) , K2 Ra ( Ra f m Cm Ce )
K1 和k2是电动机 传递系数
如果电枢电阻 Ra 和 电动机的转动惯量 J m都很小,可忽略时, 还可以简化为
Cem (t ) ua (t )
可见,数学模型是由系统结构,元件参数以及基本 运动规律所决定的。
ug ue
-
运放Ⅰ
u1
运放Ⅱ
u2
功放
ua
Mc
电动机
uf
测速
⑷各环节微分方程:
运放Ⅰ: u1 k1 (ug u f ) k1ue 功率放大: ua k3u2
1 u1 ) u2 k2 (u , 运放Ⅱ:
f
,反馈环节: u
电动机环节:
齿轮系环节:
m
kn ua kc M c Tm
例: 设线性微分方程式为
d 2c(t ) dc(t ) c(t ) r (t ) dt dt
若 r (t ) r1 (t ) 时,方程有解 c1 (t ) ,而 r (t ) r2 (t )时, c 方程有解 ,分别代入上式且将两式相加,则显 2 (t ) r2 (t ) r (t ) + r1 (t ) 时,必存在解 然有,当 为 c(t ) c1 (t ) c2 (t ) ,即为可叠加性。
微分方程求解方法
复习拉普拉斯变换有关内容(1)
2-2
2-3
由以上三式消去中间变量: ia (t ), Ea , M m (t )
工程应用中,电枢电路电感 La 较小,通常忽略不计 d m (t ) Tm m (t ) K1ua (t ) K 2 M c (t ) dt
输出变量 输入变量
电动机数学模型
R J 式中 是电动机机电时间常数 Tm a m ( Ra f m Cm Ce )
第二章 控制系统的数学模型
控制系统数学模型的建立、分类、应用等有关问题作一般性 的(概括性)讨论,为后续各章的系统分析奠定基础。 2.1 引言