VIP20
本课程属于：006-微分方程与非线性动力系统 数学计算科学计算与计算物理工程 · 共 20 门相关课程，组合学习更完整。
查看所属VIP全部课程 →
← 返回所属目录 咨询 / 购买

微分方程求解排障实战45例

共计 45 章节
可可资料掘金营 · 点击在线试看
01
常微分方程初值问题求解失败：从报错定位到方法选择
针对一个给定的一阶常微分方程初值问题 · 在求解器报错或结果明显异常时 · 通过检查方程右端函数连续性
常微分方程初值问…从报错定位到方法…
02
刚性方程识别与求解器切换排障
对一个疑似刚性的常微分方程组 · 通过计算雅可比矩阵特征值量级比 · 观察显式方法步长被压缩的现象
刚性方程识别求解器切换排障
03
步长过大导致数值解振荡的排查与收敛性验证
在固定步长求解一个衰减型常微分方程时… · 通过稳定性区域分析定位步长越界问题 · 缩小步长并做步长减半收敛性验证
步长过大导致数值…收敛性验证
04
初值敏感与病态初值问题的求解诊断
对一个对初值高度敏感的常微分方程 · 观察不同初值下解轨线的分离现象 · 区分方程本身的不稳定性与求解器误差放…
初值敏感病态初值问题的求…
05
边值问题打靶法不收敛的调试
使用打靶法求解一个二阶常微分方程边值… · 当射击参数迭代不收敛时 · 通过分析残差随射击参数的变化曲线
边值问题打靶法不…使用打靶法求解一…
06
边值问题有限差分离散后线性系统求解异常排查
对二阶边值问题做有限差分离散 · 当所得线性系统求解出现矩阵奇异或解剧… · 检查网格步长与差分格式的一致性
边值问题有限差分…对二阶边值问题做…
07
热传导方程显式格式稳定性越界排障
对一维热传导方程采用显式差分格式 · 当时间步长超过稳定性限制导致解发散时 · 通过计算网格傅里叶数定位越界原因
热传导方程显式格…对一维热传导方程…
08
热传导方程隐式格式求解线性系统效率诊断
对一维热传导方程采用隐式差分格式 · 当每步线性系统求解过慢时 · 利用三对角结构改用追赶法
热传导方程隐式格…对一维热传导方程…
09
波动方程数值色散与耗散误差的识别与修正
对一维波动方程采用显式差分格式 · 观察到波包传播出现非物理色散或幅值衰… · 通过对比不同库朗数和空间差分阶数下的…
波动方程数值色散耗散误差的识别
10
对流方程迎风格式与中心格式的数值振荡排障
对一维对流方程分别采用中心差分和迎风… · 当中心格式产生非物理振荡 · 迎风格式过度耗散时
对流方程迎风格式心格式的数值振荡…
11
扩散-对流方程边界层分辨率不足的网格加密排障
对含小扩散系数的一维对流扩散方程 · 当数值解在边界层附近出现严重误差时 · 通过估计边界层厚度检查网格分辨率
扩散对流方程边界层分…
12
非线性方程牛顿迭代不收敛的求解排障
对非线性微分方程离散后得到的非线性代… · 当牛顿迭代发散或停滞时 · 通过检查雅可比矩阵条件数
非线性方程牛顿迭…对非线性微分方程…
13
常微分方程组长时间积分能量漂移的排查
对一个保守型常微分方程组进行长时间积… · 观察到能量随时间漂移 · 通过对比不同积分器的能量误差增长趋势
常微分方程组长时…对一个保守型常微…
14
刚性化学反应动力学方程求解中的负浓度排障
对化学反应动力学常微分方程组求解时出… · 通过检查求解器容差 · 方法阶数和反应速率常数尺度
刚性化学反应动力…负浓度排障
15
微分代数方程约束违约的检测与修正
对一个指标-1微分代数方程进行数值求… · 当约束方程随时间逐渐违约时 · 通过监测约束残差
微分代数方程约束…修正
16
延迟微分方程求解步长与历史插值排障
对一阶延迟微分方程进行数值求解 · 当延迟项历史数据插值不当导致解出现异… · 检查步长与延迟量的关系
延迟微分方程求解…历史插值排障
17
偏微分方程初边值条件不一致导致的求解异常排查
对抛物型方程给定初值和边界条件 · 当求解初期出现剧烈振荡或报错时 · 检查初值与边界条件在角点处的一致性
偏微分方程初边值…对抛物型方程给定…
18
二维泊松方程迭代求解收敛过慢的预条件排障
对二维泊松方程有限差分离散后使用迭代… · 当收敛过慢时 · 通过分析谱半径与网格规模的关系
二维泊松方程迭代…对二维泊松方程有…
19
稀疏线性系统求解器选择不当导致内存溢出的排障
对大型偏微分方程离散得到的稀疏线性系… · 当直接求解器因填充导致内存溢出时 · 通过分析矩阵稀疏结构
稀疏线性系统求解…对大型偏微分方程…
20
时间步长自适应控制失效的排查与参数调整
在使用自适应步长求解器时出现步长被反… · 计算停滞的现象 · 通过检查误差估计器
时间步长自适应控…参数调整
21
求解器容差设置与局部误差估计的验证排障
对同一常微分方程分别设置不同绝对/相… · 观察全局误差与容差的关系 · 当误差未随容差收紧而下降时
求解器容差设置局部误差估计的验…
22
事件检测漏检与过检的调试
在常微分方程积分中设置事件函数检测状… · 当事件被漏检或重复触发时 · 通过检查事件函数符号变化
事件检测漏检过检的调试
23
刚性常微分方程雅可比矩阵提供错误的诊断
对刚性常微分方程向隐式求解器提供解析… · 当求解结果异常或迭代不收敛时 · 通过有限差分近似对比解析雅可比
刚性常微分方程雅…对刚性常微分方程…
24
高阶常微分方程化为一阶方程组时的变量顺序排障
将一个高阶常微分方程化为一阶方程组后… · 当结果与预期不符时 · 通过检查变量排列顺序
高阶常微分方程化…将一个高阶常微分…
25
符号求解与数值求解结果不一致的交叉验证排障
对一个可解析求解的微分方程同时进行符… · 当两者结果不一致时 · 通过检查符号解的定义域
符号求解数值求解结果不一…
26
微分方程求解脚本的单元测试与回归验证设计
为一个包含多个求解器调用的微分方程求… · 覆盖典型方程 · 边界情况和已知解析解
微分方程求解脚本…回归验证设计
27
求解结果可视化异常：坐标轴、量纲与采样排障
对微分方程数值解进行可视化时出现曲线… · 量纲混乱或采样过疏 · 通过检查数据存储格式
求解结果可视化异…坐标轴
28
并行求解多个初值问题时结果串扰的排查
并行求解同一方程的不同初值问题 · 当部分结果相互污染或随机变化时 · 通过检查全局变量
并行求解多个初值…并行求解同一方程…
29
微分方程求解程序性能瓶颈定位与优化
对一个计算耗时较长的微分方程求解程序… · 通过计时右端函数调用 · 雅可比计算和线性求解各环节
微分方程求解程序…优化
30
右端函数不连续导致求解器反复失败的排障
当常微分方程右端函数含分段不连续项时 · 求解器频繁失败或步长骤减 · 通过定位不连续点
右端函数不连续导…当常微分方程右端…
31
刚性方程隐式求解中线性求解器失败的回退策略
在隐式求解刚性方程时 · 牛顿迭代中的线性系统求解失败导致积分… · 通过检查线性系统条件数
刚性方程隐式求解线性求解器失败的…
32
偏微分方程离散格式截断误差阶验证排障
对给定偏微分方程离散格式 · 通过网格加密计算误差收敛阶 · 当实测阶数低于理论阶数时
偏微分方程离散格…对给定偏微分方程…
33
守恒律方程数值通量选择不当导致解失真的排障
对一维守恒律方程采用不同数值通量 · 当解在激波处过度抹平或产生振荡时 · 通过对比通量函数与限制器组合
守恒律方程数值通…对一维守恒律方程…
34
多维偏微分方程边界条件实现错误的定位与修正
对二维椭圆方程实现狄利克雷和诺伊曼混… · 当解在边界附近异常时 · 通过检查边界节点索引
多维偏微分方程边…修正
35
时间推进格式与空间离散匹配不当的稳定性排障
对空间高阶离散配合显式时间推进 · 当出现不稳定增长时 · 通过分析空间算子特征值与时间推进稳定…
时间推进格式空间离散匹配不当…
36
微分方程反问题参数辨识不收敛的排障
对含未知参数的微分方程进行参数辨识 · 当优化迭代不收敛或参数漂移时 · 通过检查灵敏度矩阵
微分方程反问题参…对含未知参数的微…
37
伴随方程求解与梯度校验排障
为微分方程参数反演实现伴随方程求梯度 · 当梯度与有限差分结果不符时 · 通过检查伴随方程终值条件
伴随方程求解梯度校验排障
38
随机微分方程数值求解强收敛性验证排障
对随机微分方程采用欧拉-丸山格式 · 当强收敛阶实测低于理论值时 · 通过固定布朗路径
随机微分方程数值…对随机微分方程采…
39
随机微分方程弱收敛性统计误差排障
对随机微分方程进行弱收敛性测试 · 当蒙特卡洛样本量不足导致统计误差掩盖… · 通过增加样本量
随机微分方程弱收…对随机微分方程进…
40
刚性随机微分方程隐式格式实现与稳定性排障
对含刚性漂移项的随机微分方程采用隐式… · 当解出现非物理爆炸时 · 通过检查隐式迭代收敛性
刚性随机微分方程…稳定性排障
41
微分方程求解结果文件读写格式错误的排障
将微分方程求解结果写入文件后再读取用… · 当数据错位或精度丢失时 · 通过检查写入格式
微分方程求解结果…将微分方程求解结…
42
跨平台求解器结果差异的排查与统一
同一微分方程在不同平台或库版本上求解… · 通过对比求解器默认容差 · 步长控制和浮点运算差异
跨平台求解器结果…统一
43
微分方程求解工作流中单位与量纲错误的排查
在包含多物理场的微分方程求解工作流中 · 当结果量级明显异常时 · 通过逐项检查方程中各单位制
微分方程求解工作…单位
44
求解器版本升级后结果回归异常的定位
升级微分方程求解库后 · 原有算例结果出现偏差 · 通过对比新旧版本默认参数
求解器版本升级后…升级微分方程求解…
45
微分方程求解任务的容错与断点续算排障
对长时间积分的微分方程任务 · 当运行中断后需要续算时 · 通过保存求解器状态
微分方程求解任务…断点续算排障