数值分析课件ch常微分方程数值解.ppt

第九章常微分方程数值解;常微分方程初始问题;举例;常微分方程解的存在唯一性;常微分方程的数值解;初始问题的欧拉方法;欧拉方法的推导;欧拉方法的推导;几个简单的数值方法;几个简单的数值方法;单步法的误差分析;算法的精度;例;高精度方法;Runge-Kutta方法;Runge-Kutta方法;Runge-Kutta方法;二级Runge-Kutta方法;二级Runge-Kutta方法;二级Runge-Kutta方法;三阶与四阶显式R-K方法 ;三阶与四阶显式R-K方法 ;三阶与四阶显式R-K方法 ;注;注; 例;;F=y-2*x/y; a=0; b=1; h=0.1; n=(b-a)/h; X=a:h:b; Y=zeros(1,n+1); Y(1)=1; for i=1:n x=X(i); y=Y(i); K1=h*eval(F); ;; 求微分方程（组）的解析解命令:;单步法的收敛性;定理:对初值问题的单步法 , 若局部截断误差为 ,且函数 对y 满足Lipschitz条件,即存在L0,使得 ; 对于欧拉方法，由于其增量函??? 就是 ，故当 关于 满足利普希茨条件时它是收敛的. ;单步法的稳定性;单步法的稳定性;单步法的稳定性; 隐式欧拉法的稳定性 ; 梯形公式的稳定性; 经典4阶龙格-库塔方法的稳定性; 研究一般方程的稳定性时， 相当于模型方程中的λ ;线性多步法;线性多步法;线性多步方法的构造 ;线性多步方法的构造 ;线性多步方法的构造 ;一些特殊的线性多步法;一些特殊的线性多步法; Simpson 2步4阶隐式公式 ;线性多步法的特点