数值分析第八章.pptVIP

例1.一阶常微分方程;例2. Logistic模型 ;一阶常微分方程初值问题:;取定步长: h,记 xn = x0 + nh, ( n = 1,2, ···, N ) Euler公式: yn+1 = yn + h f( xn, yn ) ;y’ = f (x, y);预-校方法又称为修正的Euler法,算法如下 k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) ,;设 yn= y(xn), 称 Rn+1=y(xn+1) - yn+1为局部截断误差.;若局部截断误差为: O(h p +1) ;利用台劳展开式 y(xn+1)=y(xn)+hy’(xn)+0.5h2y”(xn)+O(h3);三阶Range-Kutta公式一般形式;例4;MATLAB求解常微分方程初值问题命令： ;实验例题1 蛇形曲线的常微分方程初值问题 ;MATLAB解常微分方程初值问题命令;1. 创建罗伦茨模型右端函数的M文件 function z=flo(t,P) A=[-8./3 0 P(2);0 -10 10;-P(2) 28 -1]; z=A*P;;一阶常微分方程组的向量表示;;捕食者与被捕食者问题 ;平面向量场:;MATLAB命令求解: Y0=[100,20]; [t,Y]=ode2