二阶线性方程的结构
到这里,我们已经会处理不少一阶方程:能分离的、线性的、精确的、自治的,也知道没有解析公式时可以用数值方法看解。接下来课程进入二阶方程。二阶方程第一次让“当前位置”和“当前速度”同时成为状态的一部分。
本章不急着求常系数方程的特征根。我们先把二阶线性方程的骨架看清楚:为什么两个齐次解可以叠加,为什么要找两个基本解,Wronskian 到底在检查什么,非齐次方程的通解为什么总是“齐次通解加一个特解”。
二阶方程的解不仅由初始位置决定,还需要初始速度共同确定。
从一阶到二阶:状态多了一维
一阶初值问题通常写成
y′=f(t,y),y(t0)=y0.
只要方向场足够规则,一个初始高度 y0 就能选出一条解曲线。二阶方程不同。它含有 y′′,也就是速度 y′ 的变化率;要预测下一刻,光知道位置还不够,还要知道此刻朝哪里走、走得多快。
本章讨论的二阶线性方程一般写成
a2(t)y′′+a1(t)
在 a2(t)=0 的区间上,把它除以 a2(t),得到标准形式
y′′+p(t)y′+q(t)y=g(t).
对应的初值问题是
y′′+p(t)y′+q(
这里 y0 是初始位置,v0 是初始速度。如果 p,q,g 在同一个区间 I 上连续,且 ,那么这个初值问题在整个区间 上有唯一解。这个结论比一阶存在唯一性定理更贴近本章:二阶线性方程的两个初始数据正好决定一条解。
二阶方程需要两个初始条件,不是因为公式里会出现两个积分常数这么表面,而是因为系统状态本来就是一对数 (y,y′)。给定位置和速度,才给定了运动的完整起点。
把二阶方程改写成一阶系统会更直观。令
x1=y,x2=y′.
则
x1′=x2,x
也就是说,二阶方程其实在描述平面中的状态点 (y,y′) 如何运动。后面学习线性系统时,我们会把这种想法发展成矩阵方法。
线性算子与叠加原理
先把标准形式左边看成一个作用在函数上的机器:
L[y]=y′′+p(t)y′+q(t)y.
这个 L 叫线性微分算子。说它“线性”,意思是对任意两个足够可导的函数 u,v 和常数 α,β,都有
L[αu+βv]=αL[u]+βL[v].
这条式子来自求导的线性性。导数会把和拆成和,把常数保留下来;再加上 p(t)、q(t) 乘函数,也不会破坏加法和数乘。
线性算子与叠加原理:若两个解都满足齐次方程,则它们的任意线性组合仍是解。
齐次方程是
L[y]=0.
如果 y1 和 y2 都是齐次方程的解,那么
L[c1y1+c2y
所以
c1y1+c2y2
仍是齐次方程的解。这就是叠加原理。
叠加原理直接适用于齐次方程。若 L[y1]=g 且 L[y2]=g,一般不能说 也是同一个非齐次方程的解,因为 。非齐次方程可以叠加的是“输入”和“响应”:若 、,则 。
这个区别很重要。自由振动、电路自然响应、没有外部驱动的系统,常由齐次方程描述;外力、输入电压、供给项则出现在非齐次右端。齐次解能自由叠加,非齐次解要看右端是否也一起叠加。
齐次与非齐次:先分清零输入和外部输入
标准二阶线性方程
y′′+p(t)y′+q(t)y=g(t)
按右端分成两种情况。
当 g(t)=0 时,方程
y′′+p(t)y′+q(t)y=0
叫齐次方程。它描述没有外部输入时,系统凭借自身状态继续运动。机械系统里,这常对应自由振动;电路里,这常对应自然响应。
当 g(t)=0 时,方程
y′′+p(t)y′+q(t)y=g(t)
叫非齐次方程。右端 g(t) 是外部输入、外力、源项或驱动项。它把系统从“自由运动”推向某种受迫响应。
这两个词容易让人误会。“齐次”不是说解长得整齐,也不是说方程里所有项次数相同。这里的齐次只看右端是否为零。
不要把“线性方程”和“齐次方程”混成一件事。y′′+ty′+y=sint 是线性的,但不是齐次的; 右端是零,却不是线性的。
判断线性时,只看未知函数 y 及其导数是不是以一次幂出现,并且彼此没有相乘。允许系数依赖 t,例如
(1+t2)y′′+ety
是线性的。下面这些则不是线性的:
y′′+y2=0,y
基本解组:用两个解张成所有齐次解
二阶齐次线性方程的解集合有一个像二维平面一样的结构。只要找到两个“方向不同”的非零解 y1,y2,就可以把所有齐次解写成
yh=c1y1+c
这样的两个解叫基本解组。这里“方向不同”说得更准确一些,就是它们线性无关:若
c1y1(t)+c2y
在整个区间上恒成立,那么只能有
c1=0,c2=0.
基本解组中的两个线性无关解张成齐次方程的二维解空间,任意齐次解都可写成 c_1y_1 + c_2y_2。
为什么两个解就够?一个干净的理由来自初值。
设 y1,y2 是同一个二阶齐次线性方程的解。如果在某个 t0 处,初值向量
(y1(t
能张成整个平面,那么任何初始状态 (y0,v0) 都能写成它们的线性组合:
(y
于是函数
y=c1y1+c2y
满足这个初始位置和初始速度。由于二阶线性初值问题唯一,它就是对应初始状态的唯一齐次解。每个齐次解都有某个初始状态,所以都被这两个基本解张成。
下面的交互把这个想法放在最熟悉的例子里。对方程 y′′+y=0,两个基础解是 cost 和 sint。滑动 c 时,注意它们怎样同时控制曲线和初始数据。
基本解组不是“随便找两个解”。它要能表达任意初始位置和初始速度。对二阶齐次线性方程来说,这正是“两个线性无关解”的含义。
Wronskian:初值平面里的面积
怎样判断两个解是否真的方向不同?二阶方程里最常用的工具是 Wronskian:
W[y1
这个行列式不是凭空出现的。它就是两个初值向量
(y1(t),y1′(t)),(y
在状态平面中张开的有向面积。
Wronskian 可看作两个解向量在 (y,y′) 初值平面中张开的面积因子;面积为零表示两个方向线性相关。
如果在某个 t0 处
W[y1,y2](t0)=
那么这两个初值向量不平行,可以张成平面,因此 y1,y2 构成基本解组。若 W(t0)=0,这两个解在该点提供的初始信息塌缩到同一条直线上,不能在这里作为基本解组。
对同一个二阶齐次线性方程的两个解,还有一个特别方便的结论:只要系数 p,q 在区间 I 上连续,Wronskian 在 I 上要么处处为零,要么处处不为零。它不会在区间中间偶尔变成零又恢复。
这个结论可以从 Abel 公式看出。对
y′′+p(t)y′+q(t)y=0
的任意两个解 y1,y2,Wronskian 满足
W′(t)=−p(t)W(t).
于是
W(t)=W(t0)exp(−∫t
指数因子永远不为零,所以 W(t) 是否为零完全由 W(t0) 决定。
下面的交互把 Wronskian 当作面积来看。拖动两个向量,观察行列式如何随着平行或张开而变化。
“Wronskian 为零”和“函数线性相关”的关系有使用边界。若 y1,y2 是同一个二阶线性齐次方程在同一区间上的解,那么 W 在一点非零就足以说明它们线性无关;若只是任意两个函数,不能随便把这个结论反过来用。
非齐次通解:一个特解加所有自由响应
现在回到非齐次方程
L[y]=g(t).
假设我们已经找到一个特解 yp,满足
L[yp]=g(t).
再设 yh 是齐次方程 L[y]=0 的任意解。那么
L[yp+yh]=L[y
所以 yp+yh 仍是非齐次方程的解。
反过来,如果 Y 也是非齐次方程的一个解,那么
L[Y−yp]=L[Y]−L[yp
也就是说,任意两个非齐次解的差都是齐次解。因此所有非齐次解恰好是
y=yp+yh.
若 y1,y2 是对应齐次方程的基本解组,那么非齐次通解写成
y=yp+c1y1+
非齐次方程的通解可以看作一个特解 y_p 加上齐次方程的所有解 c_1y_1+c_2y_2。
这个结构比具体求法更重要。下一章会讲常系数齐次方程怎样找 y1,y2,再下一章会讲待定系数法怎样找 yp。不管具体方法怎么变,解的结构始终是同一个:
总响应=一个受迫响应+所有自由响应.
下面的交互用
y′′+y=1
做例子。一个特解是 yp=1,齐次解是 c1cost+c2。改变 时,右端输入没有变,改变的只是叠加在特解上的自由响应。
例题:先看结构,再算常数
设方程
y′′−y=et
已知对应齐次方程 y′′−y=0 的两个解是
y1=et,y2=e
又知道非齐次方程的一个特解是
yp=21tet.
求满足
y(0)=1,y′(0)=0
的解。
先写出结构,而不是急着代初值。由于 y1=et 与 y 是齐次方程的基本解组,一个非齐次特解是 ,所以通解为
这个例子有一个小陷阱:右端 et 恰好也是齐次解之一,所以特解不能简单猜成 Aet。本章暂时不展开待定系数法的规则,只要先记住:特解 yp 负责右端输入,齐次部分负责初始条件留下的自由度。
代初值时必须对整个通解求导,包括特解部分。常见错误是只让 c1y1+c2y2 去满足 和 ,忘了 本身也贡献初始位置和初始速度。
结构来自模型
二阶线性方程不是只为解题而设计的格式。很多模型自然会长成
my′′+by′+ky=F(t).
在弹簧-质量-阻尼系统中,y(t) 表示位移,my′′ 是惯性项,by′ 是阻尼项,ky 是回复力项, 是外力。把力平衡整理后,就得到这样的二阶线性方程。
弹簧-质量-阻尼模型展示了二阶线性方程中惯性、阻尼、回复力与外力输入的结构来源。
如果 F(t)=0,系统做自由运动;若还有阻尼,运动会逐渐衰减。若 F(t)=0,外力持续输入能量,系统的响应由“外力造成的部分”和“初始状态留下的部分”共同组成。
电路里也会出现同样的结构。串联 RLC 电路中,电荷 q(t) 常满足
Lq′′+Rq′+C1q
其中 L 是电感,R 是电阻,C 是电容,E(t) 是外加电压。形式和弹簧系统很像:
惯性或储能项+耗散项+回复项=外部输入.
这些模型会在后续章节反复出现。本章要留下的不是具体物理细节,而是同一个结构判断:
齐次部分看系统自身,非齐次右端看外部输入.
常见误区与练习
本章的公式并不多,真正容易出错的是判断范围和结构。
误区一:两个解只要不同就能组成基本解组
不同不等于线性无关。例如 et 和 2et 是两个写法不同的解,但它们只差一个常数倍,张不开二维解空间。它们的 Wronskian 为零。
误区二:忘记所有结论都在同一个区间上说
二阶线性方程的标准理论要求系数在同一个区间上连续,并且最高阶系数不为零。若方程写成
t2y′′+ty′+y=0,
就不能在跨过 t=0 的区间上直接把它当作标准二阶线性方程使用,因为最高阶系数 t2 在 0 处为零。
误区三:把一个特解当成通解
非齐次方程的一个特解只回答“外部输入至少造成一种响应”。通解还要加上所有齐次解:
y=yp+c1y1+
如果题目给了两个初始条件,却你的答案里没有两个可调常数,通常说明你还没有写出通解。
练习
- 判断下面方程是否为二阶线性方程,并说明是否齐次:
y′′+t2y′+sinty=0,
y′′+yy′=t,
(1+t2)y′′+y=et.
第一个是二阶线性齐次方程。第二个不是线性方程,因为出现了 yy′。第三个是二阶线性非齐次方程,右端是 et。
- 已知 y1=1、y2=t 是 y 的两个解。求它们的 Wronskian,并说明它们是否构成基本解组。
有
W[1,t](t)=10
- 设 L[y]=y′′+y。若 yp=1 是 的一个特解,而 是齐次方程的基本解组,写出 的通解。
通解为
y=1+c1cost+c2sint.这里 是一个特解, 是对应齐次方程 的通解。
- 若 Y1 和 Y2 都满足 L[y]=g(t),证明 满足对应齐次方程。
利用线性算子:
L[Y1−Y2]=L[Y
本章的结论可以压缩成一句话:二阶线性方程的齐次解空间是二维的,两个线性无关解构成基本解组;非齐次方程的全部解,是这个二维齐次解空间加上任意一个特解后的平移。下一章要做的事,就是在常系数情形下具体把那两个基本解找出来。