高数 #高数

7.常微分方程 I

Shane Lorien

微分方程实际上对学过物理的人而言,是相当熟悉的事情,但如何严格化则还应归属于数学,尽管就应用而言,直觉可以解决大部分问题(

从场论到微分方程

在场论中,我们讨论了所谓散度无旋,旋度无源:

以及一个函数可以漂亮地分解成某个散度和某个旋度:

Maxwell 方程

如雷贯耳的 方程组的微分形式,实际上就是用 的语言写成的。这也在场论中提到过。

正如上面所说,我们可以试着把 拆解:

物理上,大概就是静电场和涡旋电场,分别由电势的梯度和变化的磁场产生。 能不能消掉 ,得到一个只含 的方程呢?

波动方程的推导 (Derivation of Wave Equation)

1. 基础方程

推导始于真空(或无源项 )中的麦克斯韦方程组:

  • 法拉第感应定律: - 安培-麦克斯韦定律:(注:此处简化了系数,并假设无源
2. 数学推导步骤

对法拉第定律两边同时取旋度(Curl)

左侧(利用矢量算子恒等式):

关键点: 在无源区域(真空),由于 (由高斯定律得出),左侧简化为

右侧(代入安培定律):

3. 最终结论:波动方程

将左右两项合并并消去负号,得到电磁波在空间的演化规律:

展开即为:

物理意义: 这正是一个典型的波动方程(Wave Equation)。它表明电场(和磁场)不需要介质,可以通过自身场的变化在空间中以波的形式传播。

可以看到,这样一个方程全是 的导数所满足的方程,也就是物理直觉上的微分方程。数学上如何描述呢,我们先从简单的情形开始——常微分方程,即不带偏导数的微分方程。

常微分方程 (ODE)

基本定义

设定 为一个自变量, 为因变量,以及 的各阶导数

定义( 阶 ODE):

满足该等式的方程称为一个 阶常微分方程

定义(解):

若函数 在区间 上满足方程 ,则称 为该方程在区间 上的一个


2. 示例

一阶模型:指数衰减

  • 解的形式: - 特性: 解存在且不唯一(取决于常数 )。

二阶模型:简谐振动

  • 也许可以注意到的解: - 注意到方程的线性,于是找到比较一般的解:

这个比较一般的解,我们称为通解。

通解 (General Solution) 的定义

对于含有 个相互独立的任意常数 的解函数组 ,若满足以下两个条件:

  1. 独立性: 常数 在区间 上是相互独立的。

  2. 有效性: 对于每一组固定的 ,函数 上均满足方程

则称 为方程 的一个 通解。 固定常数后的函数也就被称为特解,例如简谐振动示例中的 就是取常数为


常数独立性 (Wronskian 行列式)

为了判断一组解是否能构成通解,需要引入线性无关性的判定:

定义(常数独立性):

对于函数组 ,若 上线性独立,其充要条件通常涉及 Wronskian 行列式不为零:

逻辑闭环:

  • 一阶方程有一个待定常数

  • 二阶方程有两个相互独立的待定常数

  • n 阶方程则需要 个线性独立的解来构造通解。

我们可以验证一下简谐振动的解是不是符合常数独立性。

验证:通解中常数的独立性

1. 设定函数形式

已知二阶方程的解函数族为:

对其求一阶导数(关于时间 ):

2. 构造行列式 (Wronskian 思想的变形)

为了验证 的独立性,计算函数向量对常数向量的偏导数行列式:

3. 代入计算

将偏导数结果代入矩阵:

利用行列式计算规则


结论

由于该行列式的值为 1(恒不为零),根据判定准则,常数 在解空间中是相互独立的。

这证明了 确实是该二阶微分方程的通解


一般的方程的解,我们也可以看成函数,那么也就有所谓隐函数,相应的,这里有所谓通积分。

通积分 (General Integral)

1. 定义

对于 阶常微分方程 ,如果其通解是由一个隐函数方程给出的:

若该等式在满足常数独立性的前提下,能够通过对 求导等代数运算还原为原方程 ,则称该隐式方程为原方程的 通积分

核心区别:

  • 通解:通常指显式函数形式

  • 通积分:指隐式方程形式,不一定能轻易解出显式的

实际上,也就是消除了 的导数的影响,也就是解掉了微分。


2. 案例演示:圆族与微分方程

板书通过一个几何例子展示了通积分与微分方程的转换逻辑:

  • 微分方程: - 通积分猜测: (这是一个以原点为圆心的圆族)

验证过程:

对隐式方程 关于 求导:

  1. 2.

  2. 化简得:

结论:

隐式方程 是微分方程 的通积分。如果写成显式通解,则为:


3. 几何直观

在坐标系中,这个通积分代表了一系列同心圆。微分方程 实际上描述了圆上任意一点的切线斜率必须与该点到原点的连线(半径)垂直。


会不会有解落在通解之外呢,或者说遍历通解中的常数也取不到某个解?

定义:奇解 (Singular Solution)

若函数 满足原方程 ,但它不包含在通解中(即无法通过给通解中的任意常数 赋值来得到),则称其为方程的一个奇解

典型案例分析:

1. 通解及其验证

给定微分方程:

通解形式为:

常数独立性验证:

通过对常数 求偏导来判断其是否提供了一个有效的自由度:

核心逻辑: 虽然在某些特定点该偏导数为零,但在区间内它不恒等于零,因此常数 是独立的。


2. 奇解 (Singular Solutions)

观察方程结构,可以发现两个不包含在上述通解族中的特殊恒等解:

验证:,则 。代入原方程:,等式成立。

由于这两个解无法通过给通解中的 赋值得到( 函数无法恒等于 1),它们被定义为该方程的奇解


包络与奇解 (Envelope and Singular Solution)

核心定义

针对 一阶常微分方程 (1st Order ODE)

若某个函数 在其曲线上:

  • 每一点 都与该方程通解族中的 某一个解 相切;

  • 本身也满足该微分方程。

为该方程的一个 奇解


为什么“相切”意味着“唯一性失效”?

从几何角度来看,包络线具有一种奇特的“重合”属性:

  1. 坐标相同: 在切点处,包络线与通解曲线的 值一致。

  2. 斜率相同: 由于相切,它们的导数 也完全一致。

推论:

对于同一个初值点 ,现在出现了两条不同的积分曲线(包络线本身和穿过该点的通解曲线),它们在该点具有相同的斜率。这说明在该轨迹上,微分方程解的存在唯一性定理 遭到了破坏。


判别逻辑
  • 通解族: - 包络线方程: 通常通过联立以下方程组并消去参数 获得:
  • 验证: 求出的包络线方程如果代入原 ODE 成立,且无法通过给 赋值得到,则它就是奇解。

实际上奇解有无穷多个,就上面的例子而言,我们可以在 走一段,然后沿着三角函数走到 ,然后再随便走……也许正如,当你看到一只蟑螂,屋子里已经满是蟑螂了(?)从另一个视角来看,通解甚至只是解的一小部分(


初值问题 (Cauchy 问题)

对于一个 阶常微分方程,如何通过初始条件锁定唯一轨迹。

定义

对于 阶 ODE:,寻找一个特解 使其满足如下 个初始条件:

其中 给定的初值

研究核心

对于 Cauchy 问题,数学家和物理学家最关心的三个维度:

  1. 存在性 (Existence): 是否真的能找到这样一个解?

  2. 唯一性 (Uniqueness): 是否有且仅有一个解?(对应了刚才讨论的“包络线/奇解”是否会破坏唯一性)。

  3. 长期行为 (Long-term Behavior): 随着自变量 (通常是时间 )的增加,解是趋于稳定、发散还是进入混沌状态?

怎么看存在性呢?

通解常数的独立性与逆映射 (Independence & Inverse Mapping)

1. 初值条件的方程组形式

对于一个包含 个待定常数 的通解 ,给定初值 ,我们可以建立如下方程组:

2. 核心判定:Jacobi 行列式

要使上述初值能够唯一地锁定一组常数 ,必须保证该映射在局部是可逆的。其判定条件为 Jacobi 行列式不为零:

  • 数学含义: 这一条件等价于常数的 独立性 (Independence)

  • 几何直观: 它保证了从“常数空间”到“初值空间”的映射是一个局部微分同胚。即对于每一组合理的初值,我们都能反向找到唯一对应的常数解。

3. 逆映射定理的视角

更抽象的映射描述:

  • 设映射 ,其中 是常数向量, 是初值向量。

  • 若 Jacobi 行列式非零,则根据 逆映射定理 (Inverse Function Theorem),存在逆映射:

  • 这意味着:通解之所以能通过初值确定特解,本质上是因为常数的独立性保证了映射的可逆性。

通解的局限性

通解只是我们预期的一般的规律,除了上述的奇解之外,我们可以考虑:

我们把两个方程相乘就得到一个新的微分方程,这两个方程可以风马牛不相及,但是他们的通解都是新的方程通解,有两个通解?!

具体计算

变量分离法 (Separable Variables)

1. 基本形式

对于一阶常微分方程(First-order ODE),若其形式如下:

2. 解题步骤

  • 寻找常数解(平衡解):

    观察 是否有零点。即若存在常数 使得 ,则 是原方程的一个解

  • 变量分离:

    若在某个区间上 ,则可将方程改写为:

3. 原理解析

(假设 )。

对上述分离后的方程两边关于 求导/积分,可得隐式解:

存在逆函数,则显式解为:

验证:

利用复合函数求导法则:


例题演示

方程: 这里

  1. 检查零点:

    ,得到 是方程的一个特解

  2. 分离变量(假设 ):

  1. 两边积分:

为了解出 ,对等式两边取指数:

(由于 ,故 ),则有:

去掉左边的绝对值符号,右边引入正负号:

此时,系数 可以取遍所有非零实数。

注意到我们在步骤 ① 中发现的常数解

  • 如果我们允许上述式子中的系数取 0,那么这个式子就能涵盖

  • 重新定义一个新的任意常数 ,使其包含正数、负数以及零()。

最终,方程的**通解(General Solution)**写为:

变量替换(Variable Substitution)

基本模型:

对于形如 的微分方程:

  1. 令: 2. 求导: 3. 转化: 得到关于 的可分离变量方程:

典型例题

方程:

1. 变量替换

,则有

代入原方程得:

2. 寻找奇异解(Equilibrium Solutions)

考察 的情况:

() 时,

反代回原变量,得到一组直线解

3. 通解计算(分离变量法)

时,进行积分:

  • 左侧积分推导:

    利用三角恒等式变换,积分为:

(注:利用了半角公式进行简化)

  • 得出结果:

进一步整理得:

**$ y = 2 \arctan(x + C) + \frac{\pi}{2} - x - 1 $**

注意到, 可以加上 而不影响方程成立,所以实际上的通解为:


关于解

该方程解的族群分布:

  • 特征: 解曲线被一系列斜率为 的平行直线(即奇异解 )所分割。

  • 趋势: 在这些“隔离带”之间,通解曲线呈现出周期性的波动形态。当 时,解会无限趋近于这些特定的直线。


二,齐次

1. 识别模型

当微分方程可以表示为自变量与因变量比值的函数时:

2. 引入中间变量

令:

3. 求导变换

关于 求导(应用乘积法则):

4. 建立新方程

将上述结果代入原方程

整理得:


逻辑拆解

  • 本质: 这种变换的精髓在于通过“比例化”将一个二元函数 塌缩为一元函数 ,从而将方程转化为可分离变量的形式:
  • 适用场景: 如果方程中 的每一项次数相同(齐次),那么它一定能化为 的形式。

  • 潜在风险: 在后续积分过程中,需额外注意分母 的情况,这通常对应于解空间中的射线(即平衡解)。


三,分式线性组合

一、 核心模型

方程形式为:

这种方程的处理方式取决于分子、分母所代表的两条直线的几何关系,即行列式 是否为零。


二、 情况 (a):相交直线(

当两条直线相交时,可以通过平移坐标原点到交点来消除常数项

  1. 求交点: 解线性方程组 ,得到唯一解

  2. 坐标平移:

  3. 方程转化: 代入原方程后,常数项消失,方程变为:

此时,方程转化为了标准的齐次方程,可用 进一步求解。


三、 情况 (b):平行直线(

当两条直线平行时,意味着系数成比例:

  1. 比例简化: 此时分子和分母的线性部分是相关的。
  1. 变量替换: 此时回到最简单的模式,令

  2. 结果: 转化为可分离变量的方程进行计算。


总结与观察

这展示了数学处理中“化归”的思想:

  • 如果有交点,就通过“平移”消除常数项,化归为齐次方程

  • 如果无交点(平行),就通过“整体替换”化归为一元函数方程

至此,我们已经完整收集了一阶微分方程变量替换法的三大招式:

  1. (整体替换)

  2. (齐次化替换)

  3. (坐标平移/平行替换)

一阶线性 ODE(关于 线性)

1. 基本定义与分类

方程的标准形式为:

  • 非齐次方程: —— 记作 - 齐次方程:—— 记作

2. 解的结构理论

线性方程最迷人的地方在于其解的可叠加性:

  • 齐次解的线性性质: 的两个解,则其线性组合 亦为解。其通解具有 的形式,且没有奇异解。这不难验证,代入即可。

  • 通解结构公式:

    是非齐次方程 的一个特解 是对应齐次方程 的一个非零解,则 通解可表示为:

直观理解: 这与线性代数中“非齐次线性方程组的通解 = 特解 + 齐次方程组通解”的逻辑完全一致。


齐次方程 的求解推导

利用积分因子法求解齐次方程的过程:

  1. 构造原函数: 满足 ,即

  2. 引入积分因子: 在方程 两端同时乘以 (显然 ):

  1. 逆用乘积导数法则:
  1. 积分得出结论:

即齐次方程的通解为:

这实际上很符合直觉,正如 我们会猜解 。 我们利用了指数求导不变来构造微分。又由于实际上每一步都是等价变换,所以没有奇解。

非齐次方程的求解——常数变易法 (Variation of Parameters)

1. 核心思想

对于齐次方程 的通解 ,我们大胆假设非齐次方程 的通解具有相同的形式,但将常数 替换为待定函数

2. 求解过程

继续沿用积分因子

  • 原方程: ——- 乘积分因子: 左右同乘 ,其中
  • 合并左端: 根据导数乘积法则的逆过程:
  • 两端积分:

二、 一阶线性 ODE 的通用公式

通过上述推导,也就得到了最终的通解公式

结构拆解:

  1. 特解部分: 2. 齐次解部分:

这同样步步等价,所以这并没有奇解,这是线性方程的优良性质。


示例:RL电路

1. 电路构成

  • 电源电动势: - 电感:- 电阻:

2. 建立方程

根据基尔霍夫电压定律(KVL),电路的电压平衡方程为:

设定初始条件:(即电路在 时刻接通,初始电流为零)。


数学转化与求解

1. 标准化处理

为了套用一阶线性 ODE 的公式,将方程两端同除以

令常数 ,方程简化为标准型:

2. 积分因子法求解

  • 积分因子: - 变换:- 积分:- 通解公式:

带入物理量与初始条件

1. 代回物理参数

由于 ,通解写为:

2. 确定待定常数

代入初始条件

3. 最终特解


物理意义评估

  • 暂态响应: 指数项 描述了电流从零开始增长的过程。电感 的存在产生了“反电动势”,阻碍了电流的突变。

  • 稳态电流: 时,指数项趋于 ,电流趋于稳定值:

此时电感相当于短路,遵循基本的欧姆定律。


积分因子

一、 问题的引入:全微分方程

1. 基本形式

将微分方程写成对称形式:

2. 判定条件(恰当性检验)

如果存在一个函数 使得 ,那么该方程称为全微分方程(或恰当方程)。其通解直接为:

根据全微分的性质,这要求满足:

(即


二、 非全微分方程与积分因子

如果 ,方程就不是全微分的。此时,我们寻找一个非零函数 ,使得乘以它之后方程变得“恰当”:

这里的 被称为积分因子


三、 积分因子的求解推导(以一元因子为例)

寻找通用的 通常涉及偏微分方程,非常困难。寻找一元积分因子(仅与 有关)的捷径:

1. 展开判定式:

根据乘积求导法则, 展开为:

整理得:

2. 假设

此时 ,方程简化为:

变形得到关于 的可分离变量方程:

3. 求解结论:

如果右端项 仅是 的函数,记为 ,那么积分因子为:


补充情况: 的判定

如果积分因子仅是 的函数,推导过程与 对称:

  1. 判定条件: 考察表达式

  2. 结论: 若该式仅与 有关,记为 ,则有:

  1. 积分因子公式:

二、 综合例题演练

方程: 其中

1. 恰当性检查

  • - 显然 ,方程不是全微分方程。

2. 寻找积分因子(试错过程)

  • 尝试

    计算

    结果含有 ,说明不存在仅与 有关的积分因子。

  • 尝试

    计算

    结果仅与 有关!满足条件,令

3. 计算积分因子

4. 转化与求解

方程两端同乘 ,得到新方程:

简化为:

验证可知,此时 ,已变为全微分方程。

5. 最终通解

关于 积分,对 关于 积分并合并(注意去重):

另外,注意到 时, 同时 ,那么这时候方程也成立。

齐次方程的特殊积分因子

核心定理:

对于齐次微分方程 (即 为同次数的齐次函数),如果 ,则该方程的一个积分因子为:

证明:

1. 核心工具:欧拉齐次函数定理

根据齐次函数的性质,若 次齐次的,则满足:

具体到我们的函数:

  • (1) - (2)

2. 判定目标

我们需要证明,乘以 后的方程 满足全微分条件:


3. 分步求偏导

左侧(对 求偏导):

应用商法则 ,设分母

展开分母的偏导项

消去 ,得到:

右侧(对 求偏导):

同样应用商法则:

展开分母的偏导项

消去


4. 利用欧拉定理进行等价变换

现在我们要证明 。将两式写在一起观察:

消去两边的 ,整理得:

此时,代入欧拉齐次函数定理

  • 左边 - 右边

左右相等!


二、 典型例题演练

方程: 其中 ,均为一阶齐次函数。

1. 恰当性检查

。由于 ,方程不是全微分方程。

2. 解法视角一:齐次化替换(回顾)

可以将方程写为:

这是典型的 形式,可用 求解。

3. 解法视角二:构造特殊积分因子

利用板书右侧的公式:

  • 分母计算:

  • 锁定因子:

4. 转化与积分

乘回原方程:

拆分项:

5. 最终通解

  • 第一部分是 - 第二部分是

    合并得出: