7.常微分方程 I
微分方程实际上对学过物理的人而言,是相当熟悉的事情,但如何严格化则还应归属于数学,尽管就应用而言,直觉可以解决大部分问题(
从场论到微分方程
在场论中,我们讨论了所谓散度无旋,旋度无源:
以及一个函数可以漂亮地分解成某个散度和某个旋度:
Maxwell 方程
如雷贯耳的
正如上面所说,我们可以试着把
物理上,大概就是静电场和涡旋电场,分别由电势的梯度和变化的磁场产生。
能不能消掉
波动方程的推导 (Derivation of Wave Equation)
1. 基础方程
推导始于真空(或无源项
- 法拉第感应定律:
- 安培-麦克斯韦定律: (注:此处简化了系数,并假设无源 )
2. 数学推导步骤
对法拉第定律两边同时取旋度(Curl):
左侧(利用矢量算子恒等式):
关键点: 在无源区域(真空),由于
(由高斯定律得出),左侧简化为 。
右侧(代入安培定律):
3. 最终结论:波动方程
将左右两项合并并消去负号,得到电磁波在空间的演化规律:
展开即为:
物理意义: 这正是一个典型的波动方程(Wave Equation)。它表明电场(和磁场)不需要介质,可以通过自身场的变化在空间中以波的形式传播。
可以看到,这样一个方程全是
常微分方程 (ODE)
基本定义
设定
定义(
满足该等式的方程称为一个
定义(解):
若函数
2. 示例
一阶模型:指数衰减
- 解的形式:
- 特性: 解存在且不唯一(取决于常数 )。
二阶模型:简谐振动
- 也许可以注意到的解:
或 - 注意到方程的线性,于是找到比较一般的解:
这个比较一般的解,我们称为通解。
通解 (General Solution) 的定义
对于含有
-
独立性: 常数
在区间 上是相互独立的。 -
有效性: 对于每一组固定的
,函数 在 上均满足方程 。
则称
常数独立性 (Wronskian 行列式)
为了判断一组解是否能构成通解,需要引入线性无关性的判定:
定义(常数独立性):
对于函数组
逻辑闭环:
一阶方程有一个待定常数
。 二阶方程有两个相互独立的待定常数
。 n 阶方程则需要
个线性独立的解来构造通解。
我们可以验证一下简谐振动的解是不是符合常数独立性。
验证:通解中常数的独立性
1. 设定函数形式
已知二阶方程的解函数族为:
对其求一阶导数(关于时间
2. 构造行列式 (Wronskian 思想的变形)
为了验证
3. 代入计算
将偏导数结果代入矩阵:
利用行列式计算规则
结论
由于该行列式的值为 1(恒不为零),根据判定准则,常数
这证明了
一般的方程的解,我们也可以看成函数,那么也就有所谓隐函数,相应的,这里有所谓通积分。
通积分 (General Integral)
1. 定义
对于
若该等式在满足常数独立性的前提下,能够通过对
核心区别:
通解:通常指显式函数形式
。 通积分:指隐式方程形式,不一定能轻易解出显式的
。
实际上,也就是消除了
2. 案例演示:圆族与微分方程
板书通过一个几何例子展示了通积分与微分方程的转换逻辑:
- 微分方程:
- 通积分猜测: (这是一个以原点为圆心的圆族)
验证过程:
对隐式方程
-
2. -
化简得:
结论:
隐式方程
3. 几何直观
在坐标系中,这个通积分代表了一系列同心圆。微分方程
会不会有解落在通解之外呢,或者说遍历通解中的常数也取不到某个解?
定义:奇解 (Singular Solution)
若函数
典型案例分析:
1. 通解及其验证
给定微分方程:
其通解形式为:
常数独立性验证:
通过对常数
核心逻辑: 虽然在某些特定点该偏导数为零,但在区间内它不恒等于零,因此常数
是独立的。
2. 奇解 (Singular Solutions)
观察方程结构,可以发现两个不包含在上述通解族中的特殊恒等解:
验证: 若
由于这两个解无法通过给通解中的
包络与奇解 (Envelope and Singular Solution)
核心定义
针对 一阶常微分方程 (1st Order ODE):
若某个函数
-
每一点 都与该方程通解族中的 某一个解 相切;
-
且
本身也满足该微分方程。
则
为什么“相切”意味着“唯一性失效”?
从几何角度来看,包络线具有一种奇特的“重合”属性:
-
坐标相同: 在切点处,包络线与通解曲线的
值一致。 -
斜率相同: 由于相切,它们的导数
也完全一致。
推论:
对于同一个初值点
判别逻辑
- 通解族:
- 包络线方程: 通常通过联立以下方程组并消去参数 获得:
- 验证: 求出的包络线方程如果代入原 ODE 成立,且无法通过给
赋值得到,则它就是奇解。
实际上奇解有无穷多个,就上面的例子而言,我们可以在
初值问题 (Cauchy 问题)
对于一个
定义
对于
其中
研究核心
对于 Cauchy 问题,数学家和物理学家最关心的三个维度:
-
存在性 (Existence): 是否真的能找到这样一个解?
-
唯一性 (Uniqueness): 是否有且仅有一个解?(对应了刚才讨论的“包络线/奇解”是否会破坏唯一性)。
-
长期行为 (Long-term Behavior): 随着自变量
(通常是时间 )的增加,解是趋于稳定、发散还是进入混沌状态?
怎么看存在性呢?
通解常数的独立性与逆映射 (Independence & Inverse Mapping)
1. 初值条件的方程组形式
对于一个包含
2. 核心判定:Jacobi 行列式
要使上述初值能够唯一地锁定一组常数
-
数学含义: 这一条件等价于常数的 独立性 (Independence)。
-
几何直观: 它保证了从“常数空间”到“初值空间”的映射是一个局部微分同胚。即对于每一组合理的初值,我们都能反向找到唯一对应的常数解。
3. 逆映射定理的视角
更抽象的映射描述:
-
设映射
,其中 是常数向量, 是初值向量。 -
若 Jacobi 行列式非零,则根据 逆映射定理 (Inverse Function Theorem),存在逆映射:
- 这意味着:通解之所以能通过初值确定特解,本质上是因为常数的独立性保证了映射的可逆性。
通解的局限性
通解只是我们预期的一般的规律,除了上述的奇解之外,我们可以考虑:
我们把两个方程相乘就得到一个新的微分方程,这两个方程可以风马牛不相及,但是他们的通解都是新的方程通解,有两个通解?!
具体计算
变量分离法 (Separable Variables)
1. 基本形式
对于一阶常微分方程(First-order ODE),若其形式如下:
2. 解题步骤
-
寻找常数解(平衡解):
观察
是否有零点。即若存在常数 使得 ,则 是原方程的一个解。 -
变量分离:
若在某个区间上
,则可将方程改写为:
3. 原理解析
设
对上述分离后的方程两边关于
若
验证:
利用复合函数求导法则:
例题演示
方程:
-
检查零点:
令
,得到 是方程的一个特解。 -
分离变量(假设
):
- 两边积分:
为了解出
令
去掉左边的绝对值符号,右边引入正负号:
此时,系数
注意到我们在步骤 ① 中发现的常数解
-
如果我们允许上述式子中的系数取 0,那么这个式子就能涵盖
。 -
重新定义一个新的任意常数
,使其包含正数、负数以及零( )。
最终,方程的**通解(General Solution)**写为:
变量替换(Variable Substitution)
一
基本模型:
对于形如
- 令:
2. 求导: 3. 转化: 得到关于 的可分离变量方程:
典型例题
方程:
1. 变量替换
令
代入原方程得:
2. 寻找奇异解(Equilibrium Solutions)
考察
当
反代回原变量,得到一组直线解:
3. 通解计算(分离变量法)
当
-
左侧积分推导:
利用三角恒等式变换,积分为:
(注:利用了半角公式进行简化)
- 得出结果:
进一步整理得:
**$ y = 2 \arctan(x + C) + \frac{\pi}{2} - x - 1 $**
注意到,
关于解
该方程解的族群分布:
-
特征: 解曲线被一系列斜率为
的平行直线(即奇异解 )所分割。 -
趋势: 在这些“隔离带”之间,通解曲线呈现出周期性的波动形态。当
时,解会无限趋近于这些特定的直线。
二,齐次
1. 识别模型
当微分方程可以表示为自变量与因变量比值的函数时:
2. 引入中间变量
令:
3. 求导变换
对
4. 建立新方程
将上述结果代入原方程
整理得:
逻辑拆解
- 本质: 这种变换的精髓在于通过“比例化”将一个二元函数
塌缩为一元函数 ,从而将方程转化为可分离变量的形式:
-
适用场景: 如果方程中
和 的每一项次数相同(齐次),那么它一定能化为 的形式。 -
潜在风险: 在后续积分过程中,需额外注意分母
的情况,这通常对应于解空间中的射线(即平衡解)。
三,分式线性组合
一、 核心模型
方程形式为:
这种方程的处理方式取决于分子、分母所代表的两条直线的几何关系,即行列式
二、 情况 (a):相交直线( )
当两条直线相交时,可以通过平移坐标原点到交点来消除常数项
-
求交点: 解线性方程组
,得到唯一解 。 -
坐标平移: 令
。 -
方程转化: 代入原方程后,常数项消失,方程变为:
此时,方程转化为了标准的齐次方程,可用
三、 情况 (b):平行直线( )
当两条直线平行时,意味着系数成比例:
- 比例简化: 此时分子和分母的线性部分是相关的。
-
变量替换: 此时回到最简单的模式,令
。 -
结果: 转化为可分离变量的方程进行计算。
总结与观察
这展示了数学处理中“化归”的思想:
-
如果有交点,就通过“平移”消除常数项,化归为齐次方程。
-
如果无交点(平行),就通过“整体替换”化归为一元函数方程。
至此,我们已经完整收集了一阶微分方程变量替换法的三大招式:
-
(整体替换) -
(齐次化替换) -
(坐标平移/平行替换)
一阶线性 ODE(关于 线性)
1. 基本定义与分类
方程的标准形式为:
- 非齐次方程:
—— 记作 - 齐次方程: —— 记作
2. 解的结构理论
线性方程最迷人的地方在于其解的可叠加性:
-
齐次解的线性性质: 若
为 的两个解,则其线性组合 亦为解。其通解具有 的形式,且没有奇异解。这不难验证,代入即可。 -
通解结构公式:
若
是非齐次方程 的一个特解, 是对应齐次方程 的一个非零解,则 的通解可表示为:
直观理解: 这与线性代数中“非齐次线性方程组的通解 = 特解 + 齐次方程组通解”的逻辑完全一致。
齐次方程 的求解推导
利用积分因子法求解齐次方程的过程:
-
构造原函数: 令
满足 ,即 。 -
引入积分因子: 在方程
两端同时乘以 (显然 ):
- 逆用乘积导数法则:
- 积分得出结论:
即齐次方程的通解为:
这实际上很符合直觉,正如
非齐次方程的求解——常数变易法 (Variation of Parameters)
1. 核心思想
对于齐次方程
2. 求解过程
继续沿用积分因子
- 原方程:
—— - 乘积分因子: 左右同乘 ,其中 :
- 合并左端: 根据导数乘积法则的逆过程:
- 两端积分:
二、 一阶线性 ODE 的通用公式
通过上述推导,也就得到了最终的通解公式:
结构拆解:
- 特解部分:
2. 齐次解部分:
这同样步步等价,所以这并没有奇解,这是线性方程的优良性质。
示例:RL电路
1. 电路构成
- 电源电动势:
- 电感: - 电阻:
2. 建立方程
根据基尔霍夫电压定律(KVL),电路的电压平衡方程为:
设定初始条件:
数学转化与求解
1. 标准化处理
为了套用一阶线性 ODE 的公式,将方程两端同除以
令常数
2. 积分因子法求解
- 积分因子:
- 变换: - 积分: - 通解公式:
带入物理量与初始条件
1. 代回物理参数
由于
2. 确定待定常数
代入初始条件
3. 最终特解
物理意义评估
-
暂态响应: 指数项
描述了电流从零开始增长的过程。电感 的存在产生了“反电动势”,阻碍了电流的突变。 -
稳态电流: 当
时,指数项趋于 ,电流趋于稳定值:
此时电感相当于短路,遵循基本的欧姆定律。
积分因子
一、 问题的引入:全微分方程
1. 基本形式
将微分方程写成对称形式:
2. 判定条件(恰当性检验)
如果存在一个函数
根据全微分的性质,这要求满足:
(即
二、 非全微分方程与积分因子
如果
这里的
三、 积分因子的求解推导(以一元因子为例)
寻找通用的
1. 展开判定式:
根据乘积求导法则,
整理得:
2. 假设
此时
变形得到关于
3. 求解结论:
如果右端项
补充情况: 的判定
如果积分因子仅是
-
判定条件: 考察表达式
。 -
结论: 若该式仅与
有关,记为 ,则有:
- 积分因子公式:
二、 综合例题演练
方程:
1. 恰当性检查
- 显然 ,方程不是全微分方程。
2. 寻找积分因子(试错过程)
-
尝试
:计算
。结果含有
,说明不存在仅与 有关的积分因子。 -
尝试
:计算
。结果仅与
有关!满足条件,令 。
3. 计算积分因子
4. 转化与求解
方程两端同乘
简化为:
验证可知,此时
5. 最终通解
对
另外,注意到
齐次方程的特殊积分因子
核心定理:
对于齐次微分方程
证明:
1. 核心工具:欧拉齐次函数定理
根据齐次函数的性质,若
具体到我们的函数:
- (1)
- (2)
2. 判定目标
我们需要证明,乘以
3. 分步求偏导
左侧(对
应用商法则
展开分母的偏导项
消去
右侧(对
同样应用商法则:
展开分母的偏导项
消去
4. 利用欧拉定理进行等价变换
现在我们要证明
消去两边的
此时,代入欧拉齐次函数定理:
- 左边
- 右边
左右相等!
二、 典型例题演练
方程:
1. 恰当性检查
2. 解法视角一:齐次化替换(回顾)
可以将方程写为:
这是典型的
3. 解法视角二:构造特殊积分因子
利用板书右侧的公式:
-
分母计算:
-
锁定因子:
4. 转化与积分
将
拆分项:
5. 最终通解
-
第一部分是
- 第二部分是合并得出: