MatrixCQY
首页
数学笔记
代码笔记
金融笔记
GitHub
首页
数学笔记
代码笔记
金融笔记
GitHub
  • 数学笔记

    • 数学笔记
    • 抽象代数

      • 群 (Group)
      • 群论进阶 (Advanced Group Theory)
      • 环与域 (Rings and Fields)
      • 模与代数 (Modules and Algebras)
      • 伽罗瓦理论 (Galois Theory)
      • 有限群表示论 (Representation Theory)
    • 李群与李代数

      • 矩阵李群 (Matrix Lie Groups)
      • 李代数初步 (Introduction to Lie Algebras)
    • 拓扑学

      • 拓扑空间与连续映射 (Topological Spaces)
      • 常见拓扑性质 (Topological Properties)
      • 滤子和网 (Filters and Nets)
      • 纽结理论 (Knot Theory)
    • 常微分方程

      • 常微分方程初等解法 (Elementary ODEs)
      • 线性微分方程 (Linear ODEs)
    • 偏微分方程

      • 偏微分方程的导出 (Derivation of PDEs)
    • 复变函数

      • 复数与复变函数
      • 全纯函数 (Holomorphic Functions)
      • 复变函数的展开
      • 积分变换
    • 实变函数

      • 预备知识 (Preliminaries)
      • Lebesgue 测度
      • 可测函数
      • Lebesgue 积分
    • 概率论

      • 事件与概率
      • 随机变量
      • 极限定理
      • 常用分布
    • 数理统计

      • 数理统计的基本概念
      • 参数估计与假设检验
    • 数值分析

      • 插值与样条
      • 矩阵与线性方程组
      • 非线性方程数值解法
      • 常微分方程数值解法
      • 偏微分方程的 Ritz-Galerkin 方法

线性微分方程 (Linear ODEs)

1. 线性微分方程组理论

一般形式: x′=A(t)x+f(t)\mathbf{x}' = A(t)\mathbf{x} + \mathbf{f}(t)x′=A(t)x+f(t),其中 x∈Rn\mathbf{x} \in \mathbb{R}^nx∈Rn.

  • 解的存在唯一性: 若 A(t),f(t)A(t), \mathbf{f}(t)A(t),f(t) 连续,则初值问题有唯一解。
  • 齐次方程: x′=A(t)x\mathbf{x}' = A(t)\mathbf{x}x′=A(t)x.
    • 解空间是 nnn 维线性空间。
    • 基解矩阵 (Fundamental Matrix) Φ(t)\Phi(t)Φ(t): 由 nnn 个线性无关解组成的矩阵。
    • 通解: x(t)=Φ(t)c\mathbf{x}(t) = \Phi(t)\mathbf{c}x(t)=Φ(t)c.
  • 非齐次方程: x(t)=Φ(t)c+Φ(t)∫Φ−1(s)f(s)ds\mathbf{x}(t) = \Phi(t)\mathbf{c} + \Phi(t) \int \Phi^{-1}(s)\mathbf{f}(s) dsx(t)=Φ(t)c+Φ(t)∫Φ−1(s)f(s)ds (常数变易公式).

2. 常系数线性方程组

x′=Ax\mathbf{x}' = A\mathbf{x}x′=Ax,其中 AAA 是常数矩阵。

  • 矩阵指数: 基解矩阵为 eAte^{At}eAt.
  • 解: x(t)=eAtx(0)\mathbf{x}(t) = e^{At} \mathbf{x}(0)x(t)=eAtx(0).
  • 计算 eAte^{At}eAt:
    • 若 AAA 可对角化 A=PΛP−1A = P \Lambda P^{-1}A=PΛP−1,则 eAt=PeΛtP−1e^{At} = P e^{\Lambda t} P^{-1}eAt=PeΛtP−1.
    • 若 AAA 有 Jordan 标准型 J=P−1APJ = P^{-1} A PJ=P−1AP,则 eAt=PeJtP−1e^{At} = P e^{Jt} P^{-1}eAt=PeJtP−1.

3. 高阶常系数线性方程

y(n)+a1y(n−1)+⋯+any=f(x)y^{(n)} + a_1 y^{(n-1)} + \dots + a_n y = f(x)y(n)+a1​y(n−1)+⋯+an​y=f(x).

  • 特征方程: λn+a1λn−1+⋯+an=0\lambda^n + a_1 \lambda^{n-1} + \dots + a_n = 0λn+a1​λn−1+⋯+an​=0.
  • 根据特征根 λ\lambdaλ 的情况(单实根、复根、重根)写出基解。
    • 单实根 λ\lambdaλ: eλxe^{\lambda x}eλx.
    • kkk 重根 λ\lambdaλ: eλx,xeλx,…,xk−1eλxe^{\lambda x}, x e^{\lambda x}, \dots, x^{k-1} e^{\lambda x}eλx,xeλx,…,xk−1eλx.
    • 复根 α±iβ\alpha \pm i\betaα±iβ: eαxcos⁡βx,eαxsin⁡βxe^{\alpha x} \cos \beta x, e^{\alpha x} \sin \beta xeαxcosβx,eαxsinβx.

4. 稳定性理论 (Stability)

讨论自治系统 x′=f(x)\mathbf{x}' = f(\mathbf{x})x′=f(x) 的平衡点 x0\mathbf{x}_0x0​ (f(x0)=0f(\mathbf{x}_0)=0f(x0​)=0) 的稳定性。

  • Lyapunov 稳定性: 初始状态靠近,解永远靠近。
  • 渐近稳定性: 初始状态靠近,解趋于平衡点。
  • 线性化判别法: 考察 Jacobian 矩阵 J=Df(x0)J = Df(\mathbf{x}_0)J=Df(x0​) 的特征值。
    • 若所有特征值实部均 <0<0<0,则渐近稳定。
    • 若有一个特征值实部 >0>0>0,则不稳定。
最近更新: 2026/2/15 04:15
Contributors: MatrixCQY
Prev
常微分方程初等解法 (Elementary ODEs)