常微分方程每章总结指南：从一阶方程到数值方法的知识主线梳理

1. 为什么“每章总结”比刷题更重要常微分方程这门课很多人学完的感觉是公式一大堆题型一大堆考试前翻书发现每一章都见过但合上书就串不起来。我当年第一次学的时候也是这样第一章会解可分离变量第二章会解一阶线性第三章会二阶常系数但问我“这些方法之间是什么关系”我答不上来。后来带过几轮答疑才慢慢意识到问题不在于题做得少而在于没有做每章总结。所谓每章总结不是把课本目录抄一遍也不是把公式表贴出来而是回答三个问题这一章引入了什么新对象、解决了什么旧方法解决不了的问题、它的方法和上一章是什么关系。常微分方程这门课有一条非常清晰的主线从最简单的方程类型出发每遇到一种解不了的形式就引入一种新工具。一阶方程搞定了发现二阶线性解不了于是引入特征方程常系数搞定了发现变系数解不了于是引入级数解和特殊函数再往后发现大多数方程根本没有初等解于是转向定性理论和数值方法。这条主线如果你不总结学到的就是一堆孤立的技巧总结之后你会发现整本书其实是在讲一个故事。这篇内容我就按章节把这条线捋一遍每一章讲清楚它的核心对象、核心方法、方法背后的动机以及和其他章节的衔接关系。适合正在学这门课、准备考试、或者想重新梳理一遍的人看。我会尽量用大白话解释但该严谨的地方不会含糊。2. 第一章基本概念与一阶方程2.1 这一章到底在定义什么第一章的核心任务其实只有一件事把“微分方程”这个东西说清楚。什么是阶、什么是线性、什么是解、什么是通解和特解、什么是初值问题。这些概念看起来简单但后面每一章都在反复用。比如“线性”这个概念第一章只是给个定义到第二章讲一阶线性方程时才真正体现出它的威力到第四章讲高阶线性方程时你才会明白为什么线性这么重要。我见过很多同学第一章学得很随意觉得不就是几个定义吗。结果到后面判断一个方程是不是线性、能不能用叠加原理全靠猜。这里我建议你把线性的定义抄下来贴在书桌前方程中未知函数及其各阶导数都是一次的系数只与自变量有关。这个定义在整本书里反复出现值得花时间吃透。阶的概念也值得多说一句。阶就是方程中出现的最高阶导数。一阶方程只含 y二阶方程含 y。为什么阶这么重要因为阶决定了通解中任意常数的个数。n 阶方程的通解含 n 个任意常数这是后面所有章节都在用的基本事实。初值条件给几个就确定几个常数这个逻辑贯穿始终。2.2 一阶方程的四种解法及其适用边界第一章真正要掌握的解法有四种可分离变量、齐次方程、一阶线性方程、伯努利方程。这四种不是随便列的它们之间有明确的层次关系。可分离变量是最基础的形式是 y f(x)g(y)把 y 和 x 分到两边分别积分就行。这个方法之所以基础是因为它不需要任何技巧只要你能分离就能解。但它的局限也很明显大多数方程分离不了。齐次方程是 y f(y/x) 的形式通过令 u y/x 把它化成可分离变量。这里的关键是识别什么时候一个方程是齐次的看它能不能写成 y/x 的函数。我当年总是判断错后来总结了一个笨办法把方程里的 x 和 y 同时换成 tx 和 ty如果方程形式不变那就是齐次的。这个方法虽然笨但从来不会错。一阶线性方程是 y p(x)y q(x) 的形式用积分因子法解。积分因子是 e 的 p(x) 积分次方这个公式必须背下来因为后面高阶线性方程、方程组都会用到类似的思想。这里我要强调一点积分因子法的推导过程比公式本身更重要。它的核心思想是“把左边凑成一个乘积的导数”这个思想在后面反复出现。伯努利方程是 y p(x)y q(x)y^n 的形式通过令 z y^(1-n) 化成一阶线性。这个方法的动机很直接非线性项 y^n 太讨厌了想办法把它消掉。换元之后确实变成线性的了这就是换元法的威力。方程类型标准形式核心方法关键换元可分离变量y f(x)g(y)分离后两边积分无齐次方程y f(y/x)化为可分离变量u y/x一阶线性y p(x)y q(x)积分因子法无伯努利方程y p(x)y q(x)y^n化为一阶线性z y^(1-n)2.3 一阶方程的几何解释与存在唯一性第一章还有一个容易被忽略但很重要的内容方向场和存在唯一性定理。方向场是把方程 y f(x,y) 在每一点画一个小线段线段斜率等于 f(x,y)。把所有小线段画出来解曲线就是沿着这些线段走的曲线。这个几何图像非常直观它告诉你微分方程的解在几何上是什么样子。存在唯一性定理说的是如果 f(x,y) 和它对 y 的偏导数在某区域连续那么过区域内每一点有且只有一条解曲线。这个定理的意义在于它告诉你什么时候方程一定有解、什么时候解是唯一的。很多同学觉得这个定理没用考试也不怎么考。但我要说这个定理是后面所有定性理论的基础。没有它你连“解是否存在”都不知道更别说研究解的性质了。我个人的经验是第一章的概念一定要慢不要急着刷题。概念清楚了后面每一章都会轻松很多。概念不清楚后面每一章都是在硬记公式。3. 第二章高阶线性方程的理论框架3.1 从二阶线性方程说起第二章通常从二阶线性方程开始形式是 y p(x)y q(x)y f(x)。这里要区分齐次和非齐次f(x) 0 是齐次f(x) ≠ 0 是非齐次。为什么先讲二阶因为二阶是最简单的能体现所有线性理论特征的阶数而且物理上二阶方程对应牛顿第二定律应用最广。这一章的核心是线性方程解的结构。齐次方程的通解是两个线性无关解的线性组合非齐次方程的通解是齐次通解加上一个特解。这个结构定理是整章的灵魂后面所有方法都是在找这两个东西齐次的两个线性无关解和非齐次的一个特解。我当年学这一章的时候最大的困惑是为什么齐次方程的通解一定是两个解的线性组合后来想明白了这其实是线性代数中“解空间是向量空间”这个事实在微分方程中的体现。齐次线性方程的所有解构成一个向量空间这个空间的维数等于方程的阶数。二阶方程的解空间是二维的所以需要两个基也就是两个线性无关的解。3.2 朗斯基行列式与线性无关的判定判断两个解是否线性无关用的是朗斯基行列式。对于两个函数 y1 和 y2朗斯基行列式 W(y1,y2) y1y2 - y1y2。如果 W 在某点不为零那么这两个函数线性无关。这个判据在第二章反复使用必须熟练掌握。这里有一个容易混淆的点朗斯基行列式不为零是线性无关的充分条件但不是必要条件。对于一般的函数W 恒为零也可能线性无关。但对于线性方程的解W 要么恒为零要么处处不为零不会出现时零时不零的情况。这个性质是阿贝尔定理给出的它把朗斯基行列式和方程本身联系起来了。阿贝尔定理说的是对于方程 y p(x)y q(x)y 0朗斯基行列式满足 W p(x)W 0所以 W C·exp(-∫p(x)dx)。这个公式非常有用它告诉你只要知道一个解就能通过朗斯基行列式求出另一个解。这就是降阶法的理论基础。3.3 降阶法与常数变易法降阶法是说如果已知齐次方程的一个非零解 y1那么可以令 y y1·u代入方程后得到一个关于 u 的一阶方程解出 u 再积分就得到 u从而得到第二个解 y2 y1·u。这个方法的动机很自然既然已经有一个解了就利用它把方程降一阶。常数变易法是求非齐次方程特解的方法。它的思想是把齐次通解中的常数换成函数然后代入非齐次方程确定这些函数。对于二阶方程设 y u1y1 u2y2其中 y1 和 y2 是齐次的两个线性无关解。代入非齐次方程后加上一个附加条件 u1y1 u2y2 0可以解出 u1 和 u2再积分得到 u1 和 u2。这两个方法的关系是降阶法解决齐次方程求第二个解的问题常数变易法解决非齐次方程求特解的问题。它们都依赖于已知齐次方程的解所以齐次方程的解是整章的基础。注意常数变易法中附加条件 u1y1 u2y2 0 不是随便加的它的目的是让代入后的方程尽可能简单。如果不加这个条件你会得到两个未知函数但只有一个方程无法求解。这个条件的选取是经过精心设计的值得仔细体会。3.4 这一章和第一章的关系第一章讲一阶线性方程用的是积分因子法。第二章讲二阶线性方程用的是常数变易法。表面上看方法不同但本质是一样的都是想办法把方程化成可以积分的形式。积分因子法是把左边凑成乘积的导数常数变易法是把非齐次项的影响分离出来。如果你理解了第一章积分因子法的推导第二章常数变易法的推导会容易很多。另外第一章的存在唯一性定理在第二章也有对应高阶线性方程也有存在唯一性定理条件类似结论类似。这个定理保证了后面所有讨论的前提解是存在的而且是唯一的。4. 第三章常系数线性方程的解法4.1 特征方程法的来龙去脉第三章的核心是常系数线性方程也就是 p(x) 和 q(x) 都是常数的情形。这一章的主角是特征方程法。对于方程 y py qy 0假设解的形式是 y e^(rx)代入后得到 r^2 pr q 0这就是特征方程。为什么假设解是指数形式因为指数函数的导数还是指数函数代入后可以把微分方程变成代数方程。这个想法非常巧妙用指数函数把微分运算变成乘法运算。特征方程的根有三种情况两个不同实根、两个相同实根、一对共轭复根。每种情况对应不同的通解形式。两个不同实根 r1 和 r2通解是 C1e^(r1x) C2e^(r2x)。两个相同实根 r通解是 (C1 C2x)e^(rx)。共轭复根 α ± iβ通解是 e^(αx)(C1cosβx C2sinβx)。这三种情况必须背熟因为后面所有常系数方程都归结为这三种情况。我当年总是记不住第二种情况为什么多一个 x。后来想明白了如果只有 e^(rx) 一个解那解空间就是一维的但二阶方程的解空间是二维的所以需要第二个解。第二个解是 xe^(rx)可以通过降阶法验证。这个推导过程值得自己动手做一遍做过之后就不会忘了。4.2 非齐次方程的特解求法非齐次常系数方程的特解通常用待定系数法。根据非齐次项 f(x) 的形式假设特解的形式代入方程确定系数。f(x) 常见的形式有多项式、指数函数、正弦余弦函数以及它们的乘积。待定系数法的关键是正确假设特解的形式。如果 f(x) 是多项式特解也假设为多项式如果 f(x) 是指数函数特解也假设为指数函数如果 f(x) 是正弦余弦特解也假设为正弦余弦的组合。但如果 f(x) 的形式和齐次解重复特解需要乘以 x 的若干次幂。这里有一个容易出错的地方如果 f(x) 是多种形式的和比如 f(x) x e^x特解要分别假设再相加。这个原理是叠加原理非齐次方程的解可以分解为多个非齐次方程的解之和。非齐次项形式特解假设形式注意事项Pn(x)Qn(x)若 0 是特征根乘以 xe^(ax)Ae^(ax)若 a 是特征根乘以 x 或 x^2sinβx 或 cosβxAcosβx Bsinβx若 ±iβ 是特征根乘以 xe^(ax)Pn(x)e^(ax)Qn(x)若 a 是特征根乘以 x 或 x^24.3 欧拉方程与变量替换欧拉方程是 x^2y pxy qy 0 的形式系数是 x 的幂次。这种方程可以通过令 x e^t 或 t lnx 化成常系数方程。这个换元的动机是欧拉方程的系数虽然不是常数但和 x 的幂次有关换元后可以消掉这个幂次。欧拉方程的特征方程和常系数方程类似只是把 r 换成 r(r-1) 等。具体来说x^2y 对应 r(r-1)xy 对应 ry 对应 1。这个对应关系可以通过换元推导出来推导过程值得自己动手做一遍。欧拉方程在物理和工程中很常见比如某些变截面杆的振动问题。它的解法本质上是常系数方程解法的推广理解了常系数方程欧拉方程就不难。4.4 这一章的核心逻辑第三章的核心逻辑是常系数方程可以通过特征方程法完全解决。齐次方程的通解由特征根决定非齐次方程的特解由待定系数法求出。这一章的方法非常程序化几乎不需要什么技巧只要按部就班就行。但这一章的局限也很明显它只适用于常系数方程。变系数方程怎么办这就是下一章要解决的问题。5. 第四章级数解与特殊函数5.1 为什么需要级数解变系数线性方程一般没有初等解比如贝塞尔方程、勒让德方程它们的解不是初等函数。这时候就需要用级数来表示解。级数解的思想是假设解可以展开成幂级数代入方程后确定系数。这个方法的动机是幂级数可以表示很多非初等函数而且幂级数的导数很容易计算。代入方程后通过比较同次幂的系数可以得到系数的递推关系从而求出级数解。级数解的关键是确定展开点。如果方程在 x0 处解析也就是系数在该点可以展开成幂级数那么可以在 x0 处展开。如果 x0 是奇点需要分正则奇点和非正则奇点。正则奇点可以用弗罗贝尼乌斯方法非正则奇点一般很难处理。5.2 弗罗贝尼乌斯方法弗罗贝尼乌斯方法是级数解的推广适用于正则奇点。它的思想是假设解的形式是 x^r 乘以一个幂级数代入方程后得到关于 r 的指标方程解出 r 后确定级数系数。指标方程是弗罗贝尼乌斯方法的核心。对于方程 x^2y xp(x)y q(x)y 0其中 p(x) 和 q(x) 在 x 0 处解析指标方程是 r(r-1) p(0)r q(0) 0。解出 r 后根据两个根的关系分情况讨论。这里的情况比较多两个根之差不是整数、是整数、是零每种情况对应的解的形式不同。我当年学这一章的时候觉得情况太多记不住。后来总结了一个原则先求指标方程再根据根的关系决定解的形式。根的关系只有三种每种对应一种处理方式记住这个框架就不会乱。5.3 贝塞尔方程与勒让德方程贝塞尔方程是 x^2y xy (x^2 - ν^2)y 0它的解是贝塞尔函数。贝塞尔函数在物理和工程中非常常见比如圆柱对称的振动问题、热传导问题。贝塞尔方程在 x 0 处是正则奇点可以用弗罗贝尼乌斯方法求解。勒让德方程是 (1-x^2)y - 2xy n(n1)y 0它的解是勒让德多项式。勒让德多项式在球坐标下的拉普拉斯方程中出现比如球对称的电场、引力场问题。勒让德方程在 x ±1 处是奇点但在 x 0 处是解析的可以用幂级数法求解。这两个方程是特殊函数的代表它们的解不是初等函数但有很多重要性质。这一章的目的不是让你记住所有特殊函数的性质而是让你掌握级数解的方法遇到新的变系数方程时知道怎么处理。5.4 这一章和前面章节的关系第四章是第三章的推广第三章解决常系数方程第四章解决变系数方程。方法从特征方程法变成级数法但核心思想是一样的假设解的形式代入方程确定参数。这一章也是整本书的转折点从求初等解转向求级数解。再往后你会发现大多数方程连级数解都求不出来只能研究解的定性性质。这就是下一章的内容。6. 第五章定性理论与稳定性6.1 从求解析解到研究解的性质第五章的核心转变是不再追求求出解的表达式而是研究解的性质比如解是否稳定、是否有周期解、长期行为如何。这个转变的动机是大多数非线性方程根本求不出解析解但实际问题中我们往往只需要知道解的大致行为。定性理论的核心工具是相平面和相图。对于二阶自治方程可以把方程写成两个一阶方程组成的方程组然后在相平面上画出轨线。相图的形状可以告诉你解的所有定性性质平衡点在哪里、是否稳定、是否有闭轨。平衡点是方程组右端为零的点也就是 y 0 和 x 0 同时成立的点。平衡点的稳定性可以通过线性化来判断在平衡点处求雅可比矩阵根据特征值的实部符号判断稳定性。实部都为负则稳定有一个实部为正则不稳定实部为零则需要进一步分析。6.2 李雅普诺夫方法与稳定性判据李雅普诺夫方法是判断稳定性的另一种方法不需要求解方程只需要构造一个李雅普诺夫函数。如果存在一个函数 V(x,y)在平衡点附近正定且沿轨线的导数负定那么平衡点是渐近稳定的。这个方法的动机是能量类比如果系统的能量函数随时间递减那么系统最终会停在能量最低点也就是平衡点。李雅普诺夫函数就是这种能量函数的推广。李雅普诺夫方法的难点在于构造 V 函数。没有通用的构造方法需要根据具体问题具体分析。常见的构造方法有二次型、能量函数、以及一些经验公式。这一章的重点不是构造技巧而是理解稳定性的概念和判据的逻辑。6.3 周期解与极限环周期解是相平面上的闭轨对应系统的周期运动。极限环是孤立的闭轨附近的轨线要么趋近它要么远离它。极限环在物理和生物中很常见比如心脏跳动、神经元放电、捕食者-被捕食者系统的周期振荡。判断极限环存在性的方法有庞加莱-本迪克松定理如果相平面上存在一个有界闭区域轨线进入后不再出去且区域内没有平衡点那么区域内存在闭轨。这个定理是定性理论中最重要的存在性定理之一。极限环的稳定性分为稳定、不稳定和半稳定。稳定极限环吸引附近的轨线不稳定极限环排斥附近的轨线。判断极限环稳定性的方法有庞加莱映射和弗洛凯理论但这些内容通常超出基础课的范围。6.4 这一章的意义和局限第五章的意义在于它告诉你即使求不出解也能知道解的行为。这个思想在现代数学和物理中非常重要因为大多数方程都求不出解析解。但这一章的局限也很明显定性理论主要适用于二阶自治方程高阶方程和非自治方程的处理要复杂得多。而且定性理论通常只能给出定性结论不能给出定量结果。如果需要定量结果还是要靠数值方法。7. 第六章数值方法与近似解7.1 为什么需要数值方法数值方法是最后的手段当解析解求不出、定性分析也不够用时用数值方法求近似解。数值方法的核心思想是把连续的微分方程离散化用差分方程近似微分方程然后逐步计算。数值方法的重要性怎么强调都不过分。实际工程中遇到的微分方程绝大多数都是靠数值方法求解的。解析解和定性分析更多是提供理论指导和验证真正算结果还是要靠数值方法。7.2 欧拉法与改进欧拉法欧拉法是最简单的数值方法y_{n1} y_n h·f(x_n, y_n)。它的思想是用切线近似曲线步长 h 越小近似越好。欧拉法的局部截断误差是 O(h^2)全局误差是 O(h)所以它是一阶方法。改进欧拉法也叫休恩法是欧拉法的改进先用欧拉法预测一个值再用梯形公式校正。改进欧拉法的局部截断误差是 O(h^3)全局误差是 O(h^2)所以它是二阶方法。欧拉法虽然简单但稳定性差步长稍大就会发散。改进欧拉法稳定性好一些但仍然是显式方法对刚性方程不适用。7.3 龙格-库塔方法龙格-库塔方法是一类高阶数值方法最常用的是四阶龙格-库塔方法。它的思想是在每一步中多次计算 f 的值然后加权平均得到更高阶的近似。四阶龙格-库塔方法的局部截断误差是 O(h^5)全局误差是 O(h^4)所以它是四阶方法。四阶龙格-库塔方法是工程中最常用的数值方法因为它精度高、稳定性好、实现简单。它的公式看起来复杂但结构很规整很容易编程实现。方法阶数局部误差全局误差稳定性欧拉法1O(h^2)O(h)差改进欧拉法2O(h^3)O(h^2)一般四阶龙格-库塔4O(h^5)O(h^4)好7.4 数值方法的稳定性和刚性方程数值方法的稳定性是指当步长 h 固定时数值解是否保持有界。对于测试方程 y λy欧拉法的稳定条件是 |1 hλ| 1也就是 h 2/|λ|。这个条件限制了步长的上限。刚性方程是指方程中包含多个时间尺度有些分量变化很快有些变化很慢。刚性方程的数值求解需要特别小心显式方法需要非常小的步长隐式方法稳定性好但每步计算量大。常用的刚性方程求解器有后向欧拉法、梯形法、以及各种隐式龙格-库塔方法。7.5 这一章和前面章节的关系第六章是前面所有章节的补充当解析方法失效时数值方法提供了一条出路。数值方法不是替代解析方法而是和解析方法互补。解析方法提供理论洞察数值方法提供具体结果。我个人的经验是学数值方法一定要动手编程。看公式觉得懂了一编程就会发现各种问题步长怎么选、误差怎么估计、稳定性怎么判断。这些问题只有动手做过才能理解。8. 每章总结的实操方法8.1 怎么做一份有效的章节总结说了这么多章节内容最后回到方法本身怎么做一份有效的章节总结。我的做法是每章总结包含四个部分核心对象、核心方法、方法动机、章节衔接。核心对象是这一章研究什么类型的方程。核心方法是这一章用什么方法解方程。方法动机是为什么要用这个方法它解决了什么前一章解决不了的问题。章节衔接是这一章和上一章、下一章是什么关系。这四个部分写下来一章的内容就串起来了。我建议用一张 A4 纸正面写核心对象和方法背面写动机和衔接。写完之后合上书试着把这张纸复述出来。如果能复述出来说明这一章真的掌握了。8.2 常见问题和避坑指南我见过很多同学的章节总结常见的问题有只抄公式不写动机、只写方法不写适用条件、只总结一章不联系其他章。这些问题导致总结做了等于没做因为知识还是孤立的。避坑的关键是每写一个公式问自己这个公式是怎么来的、什么时候能用、什么时候不能用。每写一个方法问自己这个方法解决了什么问题、和前面的方法是什么关系。每总结一章问自己这一章在整本书中处于什么位置。提示章节总结不是一次性的工作。学完第一章做一次总结学完第二章再回头补充第一章的总结。随着学习的深入你对前面章节的理解会不断加深总结也要不断更新。8.3 从总结到应用的转化章节总结的最终目的是应用考试时能快速定位方法做研究时能快速选择工具工作中能快速判断问题类型。从总结到应用关键是多练、多反思。我的做法是每做完一道题回头看看这道题对应哪一章、哪个方法、哪个公式。如果对应不上说明总结没做到位。如果对应上了说明总结有效。长期坚持你会发现自己看到题目就能判断出用什么方法这就是总结的威力。总结要素具体内容检查标准核心对象方程类型能说出方程的一般形式核心方法解法名称能写出方法的步骤方法动机解决什么问题能说出前一章为什么解决不了章节衔接和前后章的关系能画出章节关系图8.4 我个人的总结习惯我个人的习惯是每学完一章用一张 A4 纸做总结正面写核心内容背面写问题和衔接。总结完之后找一道综合题练手看看能不能把这一章的方法和前面章节的方法结合起来用。如果能说明总结到位了如果不能说明总结还有漏洞。另外我建议把总结电子化用 Markdown 或者 LaTeX 写方便修改和检索。纸质总结适合初期整理思路电子总结适合长期积累和复习。两者结合效果最好。最后说一点章节总结不是目的而是手段。目的是把知识串起来形成自己的理解框架。框架形成了具体公式忘了可以查但框架不会忘。这就是为什么每章总结比刷题更重要。