常微分方程的差分方法-欧拉法.docxVIP

常微分方程的差分方法-欧拉法一、摘要：人类社会已迈进电子计算机时代。在今天，熟练地运用计算机进行科学计算，已成为广大科技工作者和学者的一项基本技能，数值分析的基本内容是数值算法的设计与分析，科学技术当中常常需要求解常微分方程的定解问题，本文中主要以解决此问题最简单形式（一阶方程的初值问题）来求解微分方程。虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，求解从实际问题中归结出来的微分方程主要主要靠数值解法，本文就数值解法中的差分方法进行求解微分方程。二、关键词：差分方法、初值问题、数值解法、MATLAB三、引言：科学计算不应当将计算方法片面的理解为各种算法的简单罗列和堆积，它也是一门内容丰富、思想方法深刻而有着自身理论体系的数学学科。微积分的发明是人类智慧的伟大发展。求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，求解从实际问题中归结出来的微分方程主要主要靠数值解法。怎样应用数值解法求解从实际问题中归结出来的微分方程呢？四、正文 =f, (1) y()= (2)方程（1）中含有导数项(x),这是微分方程的本质特征，也正是它难以求解的症结所在。数值解法的第一步就是设法消除其导数项，这项手续称离散化。由于差分是微分的近似运算，实现离散化的基本途径是用差商替代导数。譬如，若在点列出方程（1）：()=f(,y())并用差商替代其中的导数项(),结果有：y()y()+hf(,y())设用y（）的近似值代入上式的右端，记所得结果为,这样导出的计算公式：y()y()+hf(,y())，n=0,1,2,… (3)这就是众所周知的欧拉（Euler）格式。若初值是已知的，则据式（3）可以逐步算出数值解,,…。为简化分析，人们常在为准确，即在=y()的前提下估计误差y()-。这种误差称为局部截断误差。称一种数值方法的精度是p阶的，如果其局部截断误差为O（）。对于欧拉格式（3），假定=y()，则据方程（1）有：=y()+hf(,y())=y()+h()而按泰勒公式y()=y()+h()+(),因此有y()-=()这说明欧拉格式仅为一阶方法。例 求解初值问题解：用MATLAB写出如下程序：% By lyqmathfunction main()clc;clear all;close all;% 3??μt0 = 0;y0 = [0; 0];xmax = 1;h = 0.01;% ????n = (xmax - t0)/h;for i = 1 :n+1if i == 1t(i) = t0;y(:, i) = y0;elset(i) = t0 + (i - 1)*h;y(:, i) = y(:, i - 1) + h*test_fun(t(i-1), y(:, i-1));endend% ??í?figure;hold on;box on;plot(t, y(1, :), r-, t, y(2, :), g-);grid on;legend(x, y);xlabel(t, FontWeight, Bold, Color, r);ylabel(data, FontWeight, Bold, Color, r);title(Euler·¨?????￠·?·?3ì×é By lyqmath, FontWeight, Bold, Color, r);% ?￠·?·?3ìfunction fxy = test_fun(t,y)fxy = zeros(2, 1)fxy(1) = exp(t) - 5*y(1) - y(2);fxy(2) = y(1) + 3*y(2);当h=0.01时，得出以下图像：当h=0. 1时，得到以下图像：取其中一部分由上h=0.001和h=0.1所得的图像，可看出精度h越精，图像拟合的越好，拟合度越高，效果越佳。五、结束语虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程。求解从实际问题当中归结出来的微分方程主要靠数值解法。欧拉方法是一类重要的数值解法。这类方法回避解y(x)的函数表达式，而是寻求它在一系列离散节点上的近似值，相邻的两个节点的间距称作步长。假定步长为定数。欧拉方法是一类离散化方法，这类方法将寻求解y(x)的分析问题转化为计算离散值值的代数问题，从而使问题获得了实质性的简化。然而随之带来的困难是，由于数据量往往很大，差分方法所归结出的可能是个大规模的代数方程组。