常微分方程数值解学习课件.pptVIP

相轨线是封闭曲线，结构不稳定——一旦离开某一条闭轨线，就进入另一条闭轨线，不恢复原状。 自然界存在的周期性平衡生态系统是结构稳定的，即偏离周期轨道后，内部制约使系统恢复原状。 食饵-捕食者模型(Volterra)的缺点与改进 r1=1, N1=20, ?1=0.1, w=0.2, r2=0.5, ?2=0.18 相轨线趋向极限环 结构稳定 欧拉和龙格-库塔方法的共同点:只用yn计算yn+1 单步法 步长h?0时数值解yn无限接近解析解y(xn) 收敛性 单步法有p阶精度(局部截断误差O(hp+1)) 向前欧拉公式 改进欧拉公式 4阶龙格-库塔公式 数值算法的收敛性和稳定性 (显式)单步法: 整体误差 en=O(h) en=O(h2) en=O(h4) 单步法收敛 整体误差 en=y(xn)-yn=O(hp)?0 (h?0) 数值算法的收敛性和稳定性 稳定性 计算中舍入误差不会随步数的增加无限增大 yn的误差?n ?0?微分方程稳定 向前欧拉公式 向后欧拉公式 经典龙格-库塔公式 算法稳定 (特征根 -?) h任意 刚性现象与刚性方程 振动系统或电路系统的数学模型 现象 k=2000.5, r=1000, a=1, b=-1999.5, f(t)=1 瞬态解与稳态解 e-2000t ~快瞬态解 e-t/2 ~慢瞬态解 计算到t=0.005时