← 回到目录
高 等 数 学 · 第 七 章 7.7 / 7.8

常系数线性方程:
车身压过减速带之后,
晃几下才停

上一页说线性方程有结构。这一页把结构用在最重要的那一类上——它是弹簧、电路、地震和汽车悬架共用的同一个方程。

第 1 幕

压过减速带之后,车身接下来怎么动

前轮碾过那道坎,车身被顶起来一下。然后呢?

你坐过的车不会都一样。有的忽悠忽悠晃三四下才安定(老车、减振器快废了);有的一下就回位,干脆利落;还有的像陷在糖浆里,慢慢往回蹭

同一辆车、同一道减速带,只换减振器 —— 三种完全不同的下场

三条曲线来自同一个方程,只差一个数。把车身当质量 m、弹簧劲度 k、减振器阻尼 c,牛顿第二定律写出来就是:

m y″ + c y′ + k y = 0

弹簧往回拉(−ky),减振器按速度反抗(−cy′),两个力加起来等于 my″移到一边就是上面这行。

课本管它叫「二阶常系数线性齐次微分方程」。
修车师傅管它叫「你这减振器该换了」

这一页要回答的就是一个问题:那个 c 取多少,车身才是上面第二条曲线?
第 2 幕

为什么解法是「猜一个 erx

课本一上来就说:设 y = erx凭什么这么猜?

看方程要求什么:y″y′y 三个东西按固定比例加起来要恰好抵成 0。能互相抵消,前提是它们得是同一类东西——

也就是说,我们要找一个「求导之后形状不变」的函数。

这个函数第十六页点过名。
那一页第 7 幕说过:sin 求导会左移,多项式求导会降次,只有 ex 是求导的不动点——求完还是它自己。
当时那句话像个冷知识。它其实是整个第七章的入场券。
这才是常系数方程的全部秘密:
它把「解微分方程」这件难事,换成了「解一个初中就会的二次方程」

上一页说过第七章的策略是化归——把没见过的变成见过的。
这一次化归的落点,一路退到了初中。
照这个办法做一道:y″ + y′ − 6y = 0 的通解是什么?先自己写,再看答案。
特征方程 r² + r − 6 = 0(r+3)(r−2) = 0r₁ = 2, r₂ = −3
通解 y = C₁e2x + C₂e−3x
注意这里悄悄用了上一页两条结论:两个解可以叠加(线性),而且两个就够了——二阶方程的解空间恰好二维。
写完两个常数就可以停手,不是「差不多了」,是「满了」。
第 3 幕

二次方程有三种下场,车身也有三种

一元二次方程 mr² + cr + k = 0 的判别式 Δ = c² − 4mk。中学里 Δ 只是「有没有实根」,在这儿它是「车好不好开」。

Δ特征根通解车身表现
Δ > 0两个不同实根 r₁, r₂C₁er₁x + C₂er₂x过阻尼
慢慢蹭回去
Δ = 0重根 r(C₁ + C₂x)erx临界阻尼
不过冲,回得最快
Δ < 0共轭复根 α ± βieαx(C₁cos βx + C₂sin βx)欠阻尼
来回晃几下
先记住这两个数怎么读,第 4 幕全靠它:
复根 α ± βi 里,
α(实部)管衰减——它有多负,晃动就消失得有多快;
β(虚部)管频率——它有多大,晃得就有多密。
一个根里塞了两件事:多快停下来,和停下来之前晃多密。

课本把这三行当三个公式让你背。下一幕把这三行连成一条连续的路:让阻尼从 0 一点点加大,看它自己走过这三种情形。

第 4 幕

⚠️ 支点:拖阻尼,看那道分界线

下面是同一辆车(取 m = 1, k = 1),只有减振器阻尼 c 在变。车身在 t = 0 被压下 1 个单位后松手。

上格是车身位移随时间的曲线,灰虚线是 ±2% 的「算安定了」判据带。
下格是那两个特征根在复平面上的位置——拖滑块,盯着它们走。

阻尼 c
两个根从虚轴上出发,沿着一段圆弧往左靠,在实轴上撞头,然后分开各走各的。
撞头的那一瞬间,就是 Δ = 0,就是临界阻尼。

它不是三个孤立的公式,是同一对根在复平面上的一次连续旅行
课本那三行,只是给这趟旅行的三段路各起了个名字。
c 拖到 0 看看:两个根停在虚轴上,实部 α = 0,车身永远晃下去
这就是没有减振器的弹簧,也就是第三页那个单摆——上一页刚把它线性化成 θ″ + (g/L)θ = 0
那笔账今天结清:它的解就是这一档,纯 cos 与 sin,永不停歇,周期 2π√(L/g)。
上一页说「线性化后模型说周期与振幅无关」——现在你能看见这句话是从哪儿来的:周期只由 β 决定,而 β 里根本没有初始条件。
拖到 c ≈ 1.6 你会发现:它的安定时间比临界阻尼(c = 2)还短
那「临界阻尼回得最快」这句话,是不是错的?
不是错,是那句话被省略了半截。完整说法是:在「不许过冲」的前提下,临界阻尼最快。
实测(本页数值):c = 2.05.83 个时间单位安定、零过冲c = 1.6 只用 3.75,但它冲过了原位 1.5%
只要允许一点点过冲,就能换来明显更快的安定——这是控制工程里天天在做的交易。
⚠️ 还有一件事值得注意:安定时间随 c 变化不是平滑下降的,会一跳一跳。因为「最后一次跨出 ±2%」这个判据本身是跳变的——某一次晃动刚好压进带子里,时间就整段掉下来。
指标怎么定义,结论就长什么样。这不是数学的毛病,是你得知道自己在量什么。
第 5 幕

两个技术疑点:那个 x 和那对 cos、sin

上面那张表里有两处,课本直接给结论,学生直接背,但它们都不是硬凑的。

疑点一:重根为什么要乘一个 x

Δ = 0 时只解出一个根 r,于是只得到一个erx。可上一页刚说过:

二阶线性齐次方程的解空间,是二维的。

一维的东西填不满二维的空间。所以不是「课本想乘个 x」,是「还差一个方向」这件事逼着你去找第二个解。课本用常数变易法找,找出来正好是 xerx

这是上一页第 4 幕那张线性代数对照表的第一次实战。
当时说:解空间的维数 = 方程的阶数,一个不多一个不少。
那还只是一句漂亮话;今天它变成了一条硬约束——它告诉你「还没找完」,并且告诉你「找到两个就该停」。
会用的定理和背过的定理,差别就在这儿。

疑点二:复根怎么变成了 cos 和 sin

Δ < 0 时根是 α ± βi,按套路通解该是 e(α+βi)x 那一套。可车身位移是实实在在的位移,不该冒出虚数。桥梁是欧拉公式:

e = cos θ + i sin θ
这条公式的出处,是第十页。
excos xsin x 三个泰勒展开写出来,
你会发现 e 的展开式拆开正好就是另外两个——欧拉公式不是定义出来的,是三个级数摆在一起自己撞出来的
第十页那台「用多项式算 sin」的机器,今天用来解释汽车为什么会晃。

代进去、把共轭的两项配对,虚部自动抵消,剩下的就是表里那行 eαx(C₁cos βx + C₂sin βx)

于是那两个字母的物理含义彻底落地:
eαx包络——把晃动一路压下去的那只手(第 4 幕图上的灰虚线);
cos βx、sin βx晃动本身
一个复根,一半管「多久停」,一半管「晃多密」。
重根时的第二个解,为什么偏偏是 xerx换成 x²erx 或者 2erx 行不行?
2erx 不行——它和 erx 线性相关,还是同一个方向,张不成二维(上一页那道 cos x 与 2cos x 的闸门,一模一样)。
x²erx 也不行——代回 y″ − 2ry′ + r²y = 0 算一遍,会剩下一个 2erx ≠ 0它压根不是解
xerx 代回去恰好抵成 0,且与 erx 线性无关。两个条件都满足的,就是它。
⚠️ 注意这里的判断方式:不是「记得课本写的是 x」,而是代回去验一遍——第十二页那条规矩(微分方程和积分一样,答案对不对自己就能判)在这一章一直有效。
第 6 幕

7.8:那条「撞上就乘 x」的规则,就是共振

路面不会只有一道减速带。发动机在抖、轮胎不平衡在转、搓板路一道接一道——右边不再是 0:

m y″ + c y′ + k y = F cos ωx

上一页的结构定理说:通解 = 齐次通解 + 一个特解。齐次那半截第 4 幕已经算完(它会衰减掉),所以真正决定长期表现的是特解。

课本求特解的办法叫待定系数法:右边长什么样,就猜一个同样形状的。规则里跟着一句话——

⚠️ 「若 λ 恰是特征根,特解要乘 x(k 重根就乘 xk)。」
这句话在课本里像一条纯技术性的补丁。它其实是这一章最有物理分量的一句。

把驱动频率 ω 一点点推向车身的固有频率(这里取 ω₀ = 1,并暂时去掉阻尼看得清楚),看会发生什么:

驱动频率 ω
ω 撞上 ω₀ 的那一刻,特解公式里的分母 ω₀² − ω² 变成 0——公式失效。
接手的正是那条「乘 x」的规则,而乘上去的那个 x,物理上就是「振幅随时间线性长大」。

这就是共振。不是玄学,是特解形式在分母归零那一刻换了个长相。
你早就见过它,只是不知道它写在哪一行:
洗衣机脱水到某个转速突然整台机器共振着走位;
发动机怠速调到某个转速,方向盘开始麻;
队伍过桥要便步走,不许齐步;
所有减振设计的第一条:把结构固有频率,调离激励频率。
第 5 幕说「重根要乘 x」,这一幕说「撞上特征根要乘 x」。这两个 x 是同一件事吗?
是同一件事。
两处的病根完全一样:你想放进解里的那个形状,已经被齐次解占掉了。
重根时,erx 这个方向被用过一次,第二个解得换个花样;
共振时,你猜的特解 cos ω₀x 本身就是齐次解——代进去左边只会得到 0,永远凑不出右边那个 F。
两次都是「这个方向已经满了,必须乘个 x 挤出一个新方向」。
⚠️ 记住这个判断,比记住课本那张分情况的表管用得多:先看你猜的东西是不是齐次解,是就乘 x。
第 7 幕

那么,悬架到底该取多少阻尼

第 4 幕给出的答案好像很明确:临界阻尼最好。可真车没有一辆是按临界阻尼调的。

工程上用阻尼比 ζ = c / (2√(mk)) 来说话——它就是第 4 幕滑块下面那个数。ζ = 1 是临界。

ζ什么车为什么
0.2 – 0.3家用轿车舒适优先。允许车身晃一下,换来对细碎路面不那么硬
0.3 – 0.4运动型 / 性能车车身姿态收得快,代价是路感生硬
≈ 1几乎没有量产车回位确实不过冲,但减振器太「硬」,高频路面输入直接顶进车厢
< 0.1减振器坏了车身发飘、连续起伏、刹车点头严重 —— 直接影响制动距离
为什么不取临界?因为第 4 幕只问了一个问题:压下去松手,多久回位。
真实路面还在持续输入。阻尼越大,路面的每一颗石子就越是被「硬连接」到车身上。

ζ 的选择是一场交易:车身姿态 vs 路面隔离。没有最优解,只有取舍——
而这个取舍,是在这一页的三种情形之间挑位置。
第二十三页那个词又回来了:jerk。
那一页讲弯道,说侧向加速度的变化率决定乘客舒不舒服,缓和曲线就是为了它。
今天这一页是纵向的同一件事:ζ 太小,车身反复起伏,加速度一次次换向;ζ 太大,每颗石子都变成一记冲击。
横向靠曲率连续,纵向靠阻尼比——两页讲的都是「别让加速度突变」。
路边检查减振器的土办法:用力按下车头再松手,看车身弹几下。
正常应该弹几下?用这一页的语言说清楚为什么。
正常是「弹起来一次、再轻轻一下,就停」。
把第 4 幕拖到 c = 0.5ζ ≈ 0.25,家用车那一档):读数写着最大过冲 44.4%,明显来回 4 次
44.4% 这个数就是答案——幅度按 44% → 20% → 8.8% → 3.9% 一路掉,后面两次肉眼基本看不见了,所以你只会数出一到两下。
要是明显晃了三四下,说明 ζ 掉到了 0.1 附近:拖到 c = 0.2 看看,衰减比变成 72.9%α = −c/2 太接近 0,包络压不下去——减振器基本失效(漏油或阀件磨损)。
⚠️ 注意这里的判断方式:不是数「几下」,是看每一下剩多少。那个百分比就是 e−ζπ/√(1−ζ²),一个只由 ζ 决定的数。
⚠️ 这不只是舒适问题:车身持续起伏时轮胎的垂直载荷在剧烈变化,附着力跟着忽大忽小,制动距离和转向响应都会变差
一个二阶方程的实部,能在事故报告里出现。
顺带交代第七章剩下的两节,省得你翻到那儿发慌:
7.9 欧拉方程 x²y″ + pxy′ + qy = 0:令 x = et 换个自变量,它就变成这一页的常系数方程又是化归,不必单独当新知识记。
7.10 常系数线性方程组:多个自由度耦合在一起(比如车身的上下 + 俯仰)。方法还是特征根,只不过二次方程换成了矩阵的特征值——上一页那张线性代数对照表,到这里就完全长成了线性代数本身。
第 8 幕

整节话,就这四行

常系数的全部技巧是「猜 erx,理由是它是求导的不动点;解方程就此塌成解二次方程
判别式的三种情形 = 车身的三种脾气:过阻尼慢、临界不过冲最快、欠阻尼来回晃
复根一半管衰减一半管频率:实部 α 是包络,虚部 β 是晃动的密度
「撞上特征根就乘 x」就是共振:方向被齐次解占满了,只能挤出一个新的
第七章到此走完。
回头看这一章只做了两次分岔:能不能分离(第二十五页)、是不是线性(第二十六页),
然后在线性这条路上一直走到常系数——也就是今天这一页,唯一一类能把通解完整写出来的方程。
再往后(变系数、非线性、偏微分),人类基本就只能数值解和定性分析了。你现在站的地方,是「能算出闭式答案」的边界上。
今 天 的 任 务

八道,第 5、6 题必须并排贴在一起

前四题练手,后四题是这一页的命门。每道解出来都代回去验一遍——这一章不必问别人对不对。

  1. y″ + 3y′ + 2y = 0 的通解特征根 −1、−2。答案 C₁e^(−x) + C₂e^(−2x)
  2. y″ + 2y′ + y = 0 的通解重根 −1。答案 (C₁ + C₂x)e^(−x)。写完顺手回答:那个 x 是从哪条约束逼出来的
  3. ⚠️ 求 y″ + 4y = 0 的通解,并说明它为什么永远不会停,周期是多少根 ±2i,α = 0 所以包络恒为 1;周期 π。这是第 4 幕滑块拖到 c = 0 的那一档
  4. ⚠️ 一辆车 m = 400 kg,k = 16000 N/m求临界阻尼 c;若实际减振器 c = 4000,它属于哪一种情形?c临界 = 2√(mk) = 5059.6 N·s/m;ζ = 0.79,欠阻尼。再算一下固有频率,看看是不是接近 1 Hz
  5. ⚠️ 求 y″ + y = cos 2x 的一个特解猜 A cos 2x,代回得 A = −1/3
  6. ⚠️ 求 y″ + y = cos x 的一个特解这次撞上了:猜 A cos x 必然失败,要乘 x。答案 (x/2)sin x
  7. 把第 5、6 两题的答案并排抄在一页上,写一句:两道题只差一个数字,为什么答案的形状完全不同差别只在驱动频率有没有撞上固有频率。这句话写出来,共振就不用背了
  8. ⚠️ 单摆线性化后是 θ″ + (g/L)θ = 0。用今天的方法解出来,求 L = 1 m 时的周期;并指出答案里为什么找不到振幅T = 2π√(L/g) ≈ 2.007 s。周期只由 β 决定,而 β 里没有初始条件——第三页那笔账到此彻底结清
第 5、6 题必须两道都做,而且贴在一起
如果你在第 6 题里先按第 5 题的样子猜了 A cos x请把那个失败的过程也留在纸上——
「代进去左边直接变成 0,怎么也配不出右边」这一步,就是共振在纸面上的样子。看懂它,比记住「乘 x」值钱得多。
对完答案,每道错题旁只写一句:
我是特征根解错了重根忘了乘 x没看出驱动频率撞上了特征根复根没转成 cos/sin,还是纯算错
做 完 之 后