Skip to content

线性常微分方程组

Linear system of ordinary differential equations

形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。

条目类型
模型

形式陈述

IR 为区间,A:IRn×nb:IRn 连续。一阶线性常微分方程组为

x(t)=A(t)x(t)+b(t),tI,

b0 时称为齐次。对任意 t0I 与任意初值 x0Rn,初值问题在整个区间 I 上存在唯一解;齐次方程的解集是 n 维向量空间。以 t0 为基点的基本矩阵 Φ 是矩阵初值问题

Φ(t)=A(t)Φ(t),Φ(t0)=I

的解,其各列构成齐次解空间的一组基,且 Φ(t) 对每个 tI 可逆。非齐次问题的解由常数变易(变参数)公式给出:

x(t)=Φ(t)(x0+t0tΦ(s)1b(s)ds).

A 为常矩阵时 Φ(t)=e(tt0)A,其中 etA=k0tkAk/k! 为矩阵指数。

直觉

线性是全部结构的来源:叠加原理使齐次解集成为向量空间,于是不必对每个初值单独求解,只需一次性弄清一组基如何演化——基本矩阵就是这本总账,第 j 列记录初始方向 ej 随时间的去向,任意初值的齐次解就是 Φ(t)x0。非齐次项按 Duhamel 原理理解:时刻 s 注入的微小冲击 b(s)ds 此后按齐次动力学演化,到时刻 t 变成 Φ(t)Φ(s)1b(s)ds;把所有时刻的贡献积分起来,再叠加初值的自由演化,正是常数变易公式的内容。常系数情形下"演化"就是矩阵指数:把 A 化成对角形或标准形后,耦合的方程组解耦为一组标量的指数增长、衰减与旋转,矩阵的谱数据直接翻译成动力学行为。

例子与边界

谐振子 u=u 引入 x=(u,u) 后化为 x=AxA=(0110)。由 A2=I 可直接求和得

etA=(cost)I+(sint)A=(costsintsintcost),

解沿圆周旋转,对应一对纯虚特征值 ±i。与之对照,Jordan 块 A=(λ10λ) 给出 etA=eλt(1t01):非平凡的 Jordan 标准形块使指数因子之外又长出多项式因子 t,这正是重根"共振"项的来源。

变系数情形有一个高频错误:把解写成 exp(t0tA(s)ds)。这仅当不同时刻的矩阵两两交换(例如 A(t) 恒为对角,或 A(t)=a(t)B)时才成立;一般非交换情形必须使用时间有序指数(Peano–Baker 级数),直接把积分塞进指数通常是错的。另一个普遍成立的事实是 Liouville 公式:基本矩阵的行列式(Wronskian)满足

detΦ(t)=exp(t0ttrA(s)ds),

因而永不为零——一组解要么处处线性无关,要么处处相关,不会中途退化。最后,"解在整个系数区间上存在"是线性特有的礼遇:非线性方程如 x=x2 的解可以在有限时间内爆破,线性系统的解则总能活到系数定义区间的尽头。

推论与应用

线性系统是矩阵谱与连续时间动力学之间的词典。连续系数使右端对状态局部 Lipschitz,Picard–Lindelöf 定理先给出唯一局部解,线性增长再让解延拓到整个系数区间;Grönwall 不等式则量化初值与扰动误差。

常系数情形,全部特征值实部为负当且仅当原点渐近稳定,这一谱判据是平衡点稳定性分析的基石;在非线性动力系统的平衡点处线性化,双曲情形的局部行为由线性部分主导——但线性化只是局部近似,并不意味着非线性系统可以全局线性化。工程上,控制理论用状态转移矩阵 Φ 描述系统对输入的响应。

矩阵指数仍是本页的精确演化结构,离散时间方法怎样用放大因子近似它由绝对稳定域分析。若快速衰减谱尺度迫使显式方法采用远小于精度需求的步长,则进入刚性 ODE的计算情形。高阶标量方程通过引入各阶导数作为分量化为一阶系统,因此本页理论同时覆盖振动方程等经典对象。

参考资料
  • Gerald Teschl, Ordinary Differential Equations and Dynamical Systems, American Mathematical Society, 2012,Ch. 3, linear systems and fundamental matrices。
  • Morris W. Hirsch, Stephen Smale, and Robert L. Devaney, Differential Equations, Dynamical Systems, and an Introduction to Chaos, 3rd ed., Academic Press, 2013,Chs. 4–6, linear systems, matrix exponentials and phase portraits。
关系图谱14 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组