常微分方程课程设计.pptxVIP

常微分方程课程设计 自牛顿发明微积分以来,微分方程在描述事物运动规律上已发挥了重要的作用。实际应用问题通过数学建模所得到的方程,绝大多数是微分方程。 由于实际应用的需要,人们必须求解微分方程。然而能够求得解析解的微分方程十分有限,绝大多数微分方程需要利用数值方法来近似求解。 本实验主要研究如何用 Matlab 来计算微分方程(组)的数值解,并重点介绍一个求解微分方程的基本数值解法－－Euler折线法。问题背景和实验目的 考虑一维经典初值问题 基本思想:用差商代替微商依照 Talyor 公式,y(x) 在点 xk 处有Euler 折线法 初值问题的Euler折线法 具体步骤:等距剖分：步长： 分割求解区间 差商代替微商得方程组：分割求解区间,差商代替微商,解代数方程 为分割点k = 0, 1, 2, 、、、, n-1yk 是 y (xk) 的近似 Euler 折线法举例例:用 Euler 法解初值问题取步长 h = (2 - 0)/n = 2/n,得差分方程当 h=0、4,即 n=5 时,Matlab 源程序见 fulu1、m解: Euler 折线法源程序clearf=sym(y+2*x/y^2);a=0; b=2;h=0、4;n=(b-a)/h+1; % n=(b-a)/h;x=0; y=1;szj=[x,y];for i=1:n-1 % i=