数学实验求微分方程的解.pptxVIP

数学实验求微分方程的解会计学问题背景和实验目的第1页/共25页 自牛顿发明微积分以来，微分方程在描述事物运动规律上已发挥了重要的作用。实际应用问题通过数学建模所得到的方程，绝大多数是微分方程。 由于实际应用的需要，人们必须求解微分方程。然而能够求得解析解的微分方程十分有限，绝大多数微分方程需要利用数值方法来近似求解。 本实验主要研究如何用 Matlab 来计算微分方程（组）的数值解，并重点介绍一个求解微分方程的基本数值解法－－Euler折线法。第2页/共25页Euler 折线法 考虑一维经典初值问题 基本思想：用差商代替微商根据 Talyor 公式，y(x) 在点 xk 处有 为分割点等距剖分：步长：得方程组：初值问题的Euler折线法第3页/共25页 具体步骤：分割求解区间，差商代替微商，解代数方程 分割求解区间 差商代替微商k = 0, 1, 2, ..., n-1yk 是 y (xk) 的近似第4页/共25页Euler 折线法举例例：用 Euler 法解初值问题解：取步长 h = (2 - 0)/n = 2/n，得差分方程当 h=0.4，即 n=5 时，Matlab 源程序见 fulu1.mEuler 折线法源程序第5页/共25页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;s