实验报告书8_常微分方程数值解.docVIP

PAGE  PAGE 7 东南大学《数学实验》报告 实验内容：常微分方程数值解 一 实验目的 自己编写常微分方程初值问题的常用算法，包括折线法、改进欧拉法、4阶龙格-库塔法（不允许直接使用ode45）,并用于对ODE模型的研究。 二 预备知识 (1)熟悉各种常用ODE数值算法原理 (2)了解各种算法的精度，熟悉ode45的用法 三 实验内容与要求 1.分别编写欧拉折线法、改进欧拉法和4阶龙格-库塔法通用算法 命令欧拉法： function [t,x1,x2]=ForwardEuler(a,b,c,d,n) h=(b-a)/n; t=a:h:b; x1=[c zeros(1,n)]; x2=[d zeros(1,n)]; y=zeros(2,1); for i=1:n y=ODE(x1(i),x2(i)); x1(i+1)=x1(i)+h*y(1); x2(i+1)=x2(i)+h*y(2); end end改进欧拉法： function [t,x1,x2]=ModifiedEuler(a,b,c,d,n) h=(b-a)/n; t=a:h:b; x1=[c zeros(1,n)]; x2=[d zeros(1,n)]; for i=1:n y=ODE(x1(i),x2(i)); yn1=x1(i)+h*y(1)