学习目标
学完本章后,你应该能够:
- 用输入-系统-输出的框架描述任何动态系统,区分动态系统和静态系统。
- 对简单机械系统(阻尼器并联/串联)和电气系统(LRC 电路)列写微分方程。
- 对一阶系统求自然响应,理解时间常数 τ 和直流增益 K 的物理含义。
- 对二阶系统区分欠阻尼、临界阻尼、过阻尼三种响应形态,理解阻尼比 ζ 和自然频率 ωn 的含义。
- 用拉普拉斯变换将常系数微分方程转化为代数方程。
- 定义传递函数,理解它的前提条件(零初始状态)和组合规则。
- 对典型输入(阶跃、斜坡、正弦)写出对应的拉普拉斯变换。
先用人话理解本章在讲什么
控制理论要回答的问题很简单:给一个系统加上输入,输出会怎么变?但”会怎么变”不是一眼能看出来的——温度不会瞬间到设定值,机械臂不会瞬间停住,电路里电流也不是瞬间稳定的。原因在于系统有”记忆”:过去的输入还残留在系统里,影响着当前的输出。
这一章的任务就是建模:把一个物理系统变成数学方程,然后找到一个叫”传递函数”的工具,把微分方程变成代数运算,让后面的设计和分析变得可算。
初学者最容易卡住的地方有三个:
- 分不清并联阻尼和串联阻尼的等效公式(加法 vs 倒数加法),后面列方程容易写错。
- 不理解传递函数为什么要求”零初始状态”——这是因为传递函数描述的是系统本身的性质,和初始条件无关。
- 记住了公式但不知道 τ 和 ζ 具体意味着什么,导致分析响应曲线时没有直觉。
本章是整门课的地基。后面学稳定性、根轨迹、频率响应、PID 控制,全都要从这里出发。
核心概念
1. 系统的基本结构:输入—系统—输出
任何一个系统都可以用三个部分描述:
- 输入(激励 / excitation):你施加给系统的东西,比如力、电压、温度指令。
- 系统本身:内部的物理结构和参数,比如质量、阻尼系数、电阻、电容。
- 输出(响应 / response):你关心的结果,比如位移、电流、温度。
输入系统输出
动态系统 vs 静态系统:
- 动态系统:输出同时依赖于过去和当前的输入。关键词是”记忆”——系统内部有储能元件(弹簧、电容、电感等),使得状态不能瞬间改变。
- 静态系统:输出只取决于当前输入,和过去无关。比如一个纯电阻电路,给电压就立刻有电流,没有”记忆”。
判断标准:看系统里有没有储能元件。弹簧存弹性势能、电容存电场能、电感存磁场能,这些都意味着系统是动态的。
2. 机械系统的建模:阻尼器
阻尼器(damper)的本构关系是:力和速度成正比。
F=Bx˙
其中 B 是阻尼系数(和摩擦有关),x˙ 是速度。阻尼器本身不存能量,但它是机械网络中的基本元件。
并联阻尼器:两个阻尼器承受相同的位移和速度,总力是各阻尼器的力之和。
F=B1x˙+B2x˙=(B1+B2)x˙
等效阻尼系数:Beq=B1+B2
直觉:两个并联的阻尼器”同时”作用于同一个物体,力叠加,所以等效阻尼是直接相加。
串联阻尼器:两个阻尼器承受相同的力,但各自的位移之和等于总位移。
x˙=x˙1+x˙2,F=B1x˙1=B2x˙2
等效阻尼系数:Beq1=B11+B21
直觉:串联时力一样但速度分开,类似电路中电阻并联的倒数关系。注意不要搞反了:阻尼器串联是倒数加法,不是直接加。
3. 电气系统的建模:LRC 电路
一个由电感 L、电阻 R、电容 C 串联组成的电路,输入是电压 v(t),回路电流为 I,电荷 q=∫Idt。根据基尔霍夫电压定律(KVL):
LdtdI+RI+C1∫Idt=v(t)
如果用电荷 q 表示(I=q˙),方程变成:
Lq¨+Rq˙+C1q=v(t)
这是一个二阶常系数线性微分方程。其中:
- Lq¨:电感产生的电压降(和加速度类比)。
- Rq˙:电阻产生的电压降(和速度成正比,类似阻尼)。
- C1q:电容产生的电压降(和位移成正比,类似弹簧)。
机械和电气之间的类比关系是控制理论的基本功:力 ↔ 电压,速度 ↔ 电流,质量 ↔ 电感,阻尼 ↔ 电阻,弹簧 ↔ 电容。
4. 运动系统的建模:牛顿第二定律
物体的运动遵循:
F=mx¨
其中 m 是质量,x¨ 是加速度。定义两个状态变量:
- 位置:x
- 速度:v=x˙
状态方程(标准写法):
dtdx=v,dtdv=mF
这种写法把一个二阶方程拆成两个一阶方程,是后续状态空间分析的基础。考试中经常要求列状态方程,格式要记住。
5. 一阶系统
标准形式:
x˙+τx=Ku
其中 τ 是时间常数(单位:秒),K 是直流增益,u 是输入。
传递函数: 对上式做拉普拉斯变换(零初始条件),得到:
G(s)=U(s)X(s)=τs+1K
自然响应(令 u=0):
x(t)=Ae−t/τ
A 由初始条件决定。
时间常数 τ 的物理含义:
- τ 越小,响应越快。经过 1τ 的时间,响应衰减到初始值的 36.8%;经过 3τ 衰减到约 5%;经过 5τ 基本到达稳态。
- τ 是一阶系统的”速度标尺”。
直流增益 K 的物理含义:
- 输入为常数 u0 时,输出最终稳定在 Ku0。
- K 决定了”输入变化一点,输出最终会变化多少”。
常见错误: 混淆 τ 和 1/τ。传递函数分母写 τs+1,不是 s+τ。
6. 二阶系统
标准形式:
x¨+2ζωnx˙+ωn2x=Ku
其中 ζ 是阻尼比(无量纲),ωn 是自然频率(rad/s)。
传递函数:
G(s)=U(s)X(s)=s2+2ζωns+ωn2Kωn2
三种响应情况:
| 情况 | 条件 | 响应特征 |
|---|
| 欠阻尼 | ζ<1 | 输出振荡后收敛,有超调 |
| 临界阻尼 | ζ=1 | 最快无振荡收敛 |
| 过阻尼 | ζ>1 | 缓慢收敛,无振荡 |
ζ 的物理含义: 衡量系统的”刹车能力”。ζ 越小,振荡越剧烈,超调越大。ζ=0.7 附近是工程中常用的折中值。
ωn 的物理含义: 系统无阻尼时的自由振荡频率。ωn 越大,系统响应越快。
常见错误: 分母特征方程 s2+2ζωns+ωn2=0 的判别式是 (2ζωn)2−4ωn2,不是 4ζ2−4ωn2。推导时要小心系数。
7. 高阶系统
任何线性时不变(LTI)系统都可以分解为若干一阶和二阶子系统的组合。因此,只要搞懂了一阶和二阶系统的行为,就能分析任意高阶系统。
8. 拉普拉斯变换基础
拉普拉斯变换把时间域的微分方程变成复频域的代数方程,是控制理论的核心数学工具。
定义:
F(s)=L{f(t)}=∫0∞e−stf(t)dt
关键性质:
| 性质 | 公式 | 用途 |
|---|
| 线性 | L{af(t)+bg(t)}=aF(s)+bG(s) | 拆开求和 |
| 微分 | L{f′(t)}=sF(s)−f(0) | 消去微分符号 |
| 积分 | L{∫fdt}=sF(s) | 消去积分符号 |
| 零初值微分 | L{f′(t)}=sF(s)(当 f(0)=0) | 写传递函数时用 |
常用变换对照表:
| 时域 f(t) | 复频域 F(s) | 常见错误提醒 |
|---|
| e−at | s+a1 | a 是正数时,极点在左半平面 |
| sin(ωt) | s2+ω2ω | 分子是 ω,不是 1 |
| cos(ωt) | s2+ω2s | 分子是 s,不是 1 |
| tn | sn+1n! | 别漏掉 n! |
逆变换:
f(t)=L−1{F(s)}
从 F(s) 回到时域,通常用部分分式展开法。
9. 传递函数
定义: 在零初始条件下,输出的拉普拉斯变换与输入的拉普拉斯变换之比:
G(s)=R(s)C(s)
其中 C(s) 是输出,R(s) 是输入。
关键前提:“零初始状态”。 传递函数描述的是系统本身的输入输出关系,和初始条件无关。这不是一个可有可无的说明——如果不满足零初始条件,G(s) 就不是简单的输出比输入。
组合规则:
- 串联:G(s)=G1(s)⋅G2(s)
- 并联:G(s)=G1(s)+G2(s)
- 反馈:G(s)=1+G1(s)H(s)G1(s)(负反馈)
传递函数是 s 的有理多项式之比,分子分母的阶次和系统阶次直接相关。
10. 典型输入信号
| 输入类型 | 时域表达式 r(t) | 拉普拉斯变换 R(s) | 用途 |
|---|
| 单位阶跃 | u(t)=1,t≥0 | s1 | 测试系统的瞬态响应 |
| 单位斜坡 | r(t)=t,t≥0 | s21 | 测试系统的跟踪能力 |
| 单位正弦 | sin(ωt) | s2+ω2ω | 测试系统的频率响应 |
用传递函数分析响应时,步骤是:C(s)=G(s)⋅R(s),然后对 C(s) 做逆变换得到 c(t)。
核心公式与推导
推导 1:从微分方程到传递函数(一阶系统)
给定一阶系统:
x˙+τx=Ku
两边做拉普拉斯变换,令 x(0)=0:
sX(s)+τ1X(s)=KU(s)
(s+τ1)X(s)=KU(s)
G(s)=U(s)X(s)=s+τ1K=τs+1Kτ
这一步在干什么: 把一阶微分方程转化成 s 域的代数关系。分母的根 s=−1/τ 决定了系统的衰减速度。
推导 2:从微分方程到传递函数(二阶系统)
给定二阶系统:
x¨+2ζωnx˙+ωn2x=Kωn2u
两边拉普拉斯变换,令初始条件为零:
s2X(s)+2ζωnsX(s)+ωn2X(s)=Kωn2U(s)
G(s)=U(s)X(s)=s2+2ζωns+ωn2Kωn2
这一步在干什么: 二阶系统的传递函数分母是一个二次多项式。通过分析这个多项式的根(特征方程的根),可以判断系统是欠阻尼、临界阻尼还是过阻尼。
特征方程的根:
s=−ζωn±ωnζ2−1
- ζ<1:根是复数,实部 −ζωn 决定衰减,虚部 ωn1−ζ2 决定振荡频率。
- ζ=1:根是实重根 −ωn。
- ζ>1:根是两个实数,较靠近虚轴的那个根主导响应速度。
推导 3:阻尼器等效公式
并联:两个阻尼器速度相同 x˙1=x˙2=x˙,力相加。
Ftotal=B1x˙+B2x˙=(B1+B2)x˙⟹Beq=B1+B2
串联:力相同 F=B1x˙1=B2x˙2,速度相加。
x˙=x˙1+x˙2=B1F+B2F=F(B11+B21)⟹Beq1=B11+B21
配套例题
例题 1:弹簧-质量-阻尼系统的传递函数
题目: 一个质量为 m 的物体受外力 F(t) 驱动,通过阻尼系数 B 的阻尼器和刚度为 Ks 的弹簧连接到固定墙面。求位移 x(t) 对力 F(t) 的传递函数。
解题步骤:
第 1 步:列牛顿第二定律方程。
mx¨+Bx˙+Ksx=F(t)
这里 Bx˙ 是阻尼力,Ksx 是弹簧力,方向都与位移相反。
第 2 步:零初始条件下做拉普拉斯变换。
(ms2+Bs+Ks)X(s)=F(s)
第 3 步:写出传递函数。
G(s)=F(s)X(s)=ms2+Bs+Ks1
回答什么问题: 这个传递函数描述了力到位置的动态关系。分母是系统的特征多项式,根的位置决定系统是否振荡以及振荡多久衰减。
易错点: 忘记阻尼力前面的符号。阻尼力总是阻碍运动,所以方程中是 +Bx˙(移项后),不是 −Bx˙。
例题 2:一阶系统的自然响应
题目: 一阶系统 x˙+2x=0,初始条件 x(0)=5。求 x(t)。
解题步骤:
第 1 步:识别时间常数。方程改写为 x˙+x/τ=0,对比得 1/τ=2,所以 τ=0.5 秒。
第 2 步:写出自然响应公式。x˙+x/τ=0 的解为:
x(t)=x(0)e−t/τ
第 3 步:代入初始条件。
x(t)=5e−t/0.5=5e−2t
验证: t=0 时 x=5(满足初始条件),t→∞ 时 x→0(系统衰减到平衡态)。τ=0.5 秒意味着约 2.5 秒后响应基本衰减完。
例题 3:LRC 电路的传递函数
题目: 串联 LRC 电路中,输入为电压源 V(s),输出为电容两端电压 VC(s)。已知电感 L=1H,电阻 R=3Ω,电容 C=0.5F。求传递函数 G(s)=VC(s)/V(s),并判断阻尼比。
解题步骤:
第 1 步:用电流 I 作为中间变量列 KVL 方程。
LsI(s)+RI(s)+Cs1I(s)=V(s)
第 2 步:电容电压和电流的关系为 VC(s)=Cs1I(s),所以 I(s)=Cs⋅VC(s)。
第 3 步:代入第 1 步的方程。
(Ls+R+Cs1)CsVC(s)=V(s)
(LCs2+RCs+1)VC(s)=V(s)
G(s)=V(s)VC(s)=LCs2+RCs+11
第 4 步:代入数值 L=1,R=3,C=0.5。
G(s)=0.5s2+1.5s+11=s2+3s+22
第 5 步:对比标准二阶形式 s2+2ζωns+ωn2=0。
- ωn2=2⟹ωn=2≈1.414rad/s
- 2ζωn=3⟹ζ=223≈1.06
ζ>1,系统是过阻尼的,响应不会振荡。
易错点: 在求 ωn 和 ζ 时,一定要先把分母写成 s2+2ζωns+ωn2 的形式(s2 系数为 1),再比较系数。很多人忘记归一化就直接比较,结果算错。
例题 4:二阶系统的阶跃响应
题目: 二阶系统传递函数 G(s)=s2+2s+44。输入为单位阶跃 R(s)=1/s。求输出 C(s) 并判断响应类型。
解题步骤:
第 1 步:写出 C(s)=G(s)⋅R(s)。
C(s)=s2+2s+44⋅s1=s(s2+2s+4)4
第 2 步:分析分母特征方程。
s2+2s+4=0
ωn2=4⟹ωn=2,2ζωn=2⟹ζ=0.5。
ζ<1,系统是欠阻尼的,阶跃响应会振荡。
第 3 步:写出超调量。对于欠阻尼系统,阶跃响应的超调量为:
Mp=e−ζπ/1−ζ2×100%
代入 ζ=0.5:
Mp=e−0.5π/1−0.25=e−π/3≈e−1.814≈16.3%
需要记住的结论: ζ=0.5 时超调约 16%,这在考试中经常直接用,记住可以节省推导时间。
重点难点总结
需要区分的概念
| 概念 | 含义 | 一句话区分 |
|---|
| 动态系统 | 输出依赖于过去和当前输入 | 有”记忆”(储能元件) |
| 静态系统 | 输出只依赖当前输入 | 没有”记忆” |
| 一阶系统 | 微分方程最高一阶 | 时间常数 τ 决定响应快慢 |
| 二阶系统 | 微分方程最高二阶 | ζ 和 ωn 共同决定响应形态 |
| 传递函数 | 零初始条件下输出/输入的拉普拉斯变换比 | 描述系统本身,和初始条件无关 |
高频易错点
-
阻尼器并联 vs 串联:并联直接加(B1+B2),串联倒数加(1/B1+1/B2)。和电阻的串并联规律做对比:电阻串联是直接加,并联是倒数加。阻尼器恰好反过来。
-
传递函数的零初始条件:传递函数只能在零初始条件下定义。如果题目给定了非零初始条件,不能直接用 G(s) 乘 R(s),要先用拉普拉斯变换处理初值项。
-
二阶系统特征方程归一化:求 ζ 和 ωn 之前,务必先把分母 s2 系数除到 1。
-
自然频率的单位:ωn 的单位是 rad/s,不是 Hz。如果需要频率 f,用 f=ωn/(2π)。
-
拉普拉斯变换对的分子:sin(ωt) 的变换分子是 ω,cos(ωt) 的变换分子是 s。很多人记混。
自测题与答案
自测题
- 一个系统里有弹簧和阻尼器,没有电容也没有电感。这个系统是动态系统还是静态系统?为什么?
- 两个阻尼器串联,B1=2N⋅s/m,B2=3N⋅s/m。等效阻尼系数是多少?
- 一阶系统 x˙+5x=10u 的时间常数和直流增益分别是多少?
- 二阶系统的传递函数分母为 s2+6s+9。求 ωn 和 ζ,并判断阻尼类型。
- 求 L{3e−2t}。
- 已知 G(s)=s+32,输入为单位阶跃。求输出 C(s)。
- 传递函数 G(s)=s2+2ζωns+ωn25ωn2 中,当 ζ=0 时,系统的阶跃响应会怎样?
- LRC 串联电路中,输入电压 V(s),输出电容电压 VC(s),写出传递函数的一般形式。
自测题答案
-
动态系统。弹簧是储能元件(存储弹性势能),所以输出依赖于系统的状态历史,系统有”记忆”。
-
串联等效:Beq1=21+31=65,所以 Beq=56=1.2N⋅s/m。
-
标准形式 x˙+τx=Ku。对比 x˙+5x=10u,得 1/τ=5,所以 τ=0.2 秒;K=10/5=2(因为标准形式应写成 x˙+x/τ=Ku,即 x˙+5x=5Ku,所以 5K=10,K=2)。
-
已是归一化形式。ωn2=9⟹ωn=3,2ζωn=6⟹ζ=1。系统是临界阻尼。
-
L{3e−2t}=3⋅s+21=s+23。
-
C(s)=G(s)⋅R(s)=s+32⋅s1=s(s+3)2。
-
ζ=0 时系统无阻尼,阶跃响应会持续等幅振荡,永不衰减。物理上对应一个无摩擦的弹簧-质量系统。
-
一般形式:G(s)=V(s)VC(s)=s2+(R/L)s+1/(LC)1/(LC),或写成 G(s)=LCs2+RCs+11。
学习路线
建议按下面顺序学习:
-
先理解第 1 节”输入—系统—输出”的框架。这是分析任何系统的起点,后面每一章都在这个框架里操作。
-
再学第 2-4 节的物理建模。阻尼器等效公式要亲手推一遍,不要只看结论。LRC 电路的推导和机械系统对照看,建立起跨领域的类比直觉。
-
第 5-6 节是一阶和二阶系统的标准形式。这两个是核心中的核心。先把传递函数写对,再理解 τ、ζ、ωn 的含义。最好能手画不同 ζ 值下的响应曲线草图。
-
第 8 节拉普拉斯变换。不需要死记所有变换对,重点记住微分性质和 τ 与 s 域极点的关系。后续章节会反复用到。
-
第 9 节传递函数。理解”零初始条件”的含义,然后练习从微分方程推传递函数。组合规则(串联、并联、反馈)先记住,后面第 5-6 章会大量使用。
-
最后做例题和自测题。如果第 4 题和第 8 题卡住,说明归一化和电路建模需要重新复习。
和后续章节的关系
第 1 章解决的是”怎么描述系统”。后续章节的核心推进方向:
- 第 2-3 章(特征方程与稳定性):用第 1 章的传递函数分母(特征多项式)分析系统是否稳定。第 1 章中 ζ 和 τ 的讨论是稳定性的前奏。
- 第 4-5 章(根轨迹与频域分析):传递函数 G(s) 是根轨迹和 Bode 图的分析对象。第 1 章的拉普拉斯变换基础是频域分析的数学工具。
- 第 6-7 章(控制器设计):PID 控制器设计的目标就是把第 1 章中描述的开环系统改造成期望的闭环响应。
- 后续实验/仿真:MATLAB/Simulink 中建立系统模型时,输入的就是传递函数的分子分母系数。
第 1 章的传递函数写错了,后面全部跟着错。建议在学完本章后,找三到五个不同的物理系统,练习”从物理系统列微分方程 → 做拉普拉斯变换 → 写传递函数”这个完整流程。这个动作做到熟练,后面的章节会轻松很多。