河海大学《数值分析》课件-第9章常微分方程数值解法.pptVIP

数值分析ComputationalMethod

Chapter9常微分方程数值解法

9.常微分方程数值解法

9．1引言一阶常微分方程初值问题可表为:若函数适当光滑，且关于y满足李普希兹(Lipschitz)条件则上述问题存在唯一连续可微的解y(x)。

数值方法基本思想：在初值问题的求解区间[a,b]上插入若干个节点：称为步长数值解法的特点:采用步进式或递推式算法.由初始值y(x0)推导出x1的近似解y1，再由y1计算出y2，y3，?，直到求出yn为止。沿节点次序一步一步地向前推进.总假定等步长数值解法:将常微分方程这一连续模型离散化,求处的近似值得初值问题的解函数y(x)在离散节点

由yi-k,yi-k+1,?,yi-2,yi-1得到yi称为多步法。或k步法仅由yi-1得到yi称为单步法。所求出的yi(i=1,2,?,n)称为初值问题的数值解。

9.2.简单的数值方法与基本概念1.微分法:欧拉公式(显式)向前差商向后差商欧拉公式(隐式)

2.积分法:欧拉公式(显式)欧拉公式(隐式)

梯形公式(隐式)取初值:(预报值)(校正值)改进欧拉公式

例设(1)分别写出求解此微分方程的欧拉格式和改进欧拉格式;(2)取步长h=0.1,求y(0.2)的近似值,并与精确解y(x)=e-x2作比较.解:(1)欧拉格式: