常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc

常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc

ID:58450475

大小:111.00 KB

页数:3页

时间:2020-09-03

常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc_第1页
常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc_第2页
常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc_第3页
资源描述:

《常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、常微分方程数值解的Euler方法和改进的Euler法机械工程学院孙志利一、简介在数学和计算机科学中,欧拉方法(Eulermethod)命名自它的发明者莱昂哈德·欧拉,是一种一阶数值方法,用以对给定初值的常微分方程(即初值问题)求解。它是一种解决常微分方程数值积分的最基本的一类显型方法(Explicitmethod)。欧拉法(Eulermethod)是以流体质点流经流场中各空间点的运动即以流场作为描述对象研究流动的方法。它不直接追究质点的运动过程,而是以充满运动液体质点的空间——流场为对象。研究各时刻质点在流场中的变化规

2、律。将个别流体质点运动过程置之不理,而固守于流场各空间点。通过观察在流动空间中的每一个空间点上运动要素随时间的变化,把足够多的空间点综合起来而得出的整个流体的运动情况。其基本思想是迭代,其中分为前进的EULER法、后退的EULER法、改进的EULER法。所谓迭代,就是逐次替代,最后求出所要求的解,并达到一定的精度。误差可以很容易的计算出来。欧拉法简单地取切线的端点作为下一步的起点进行计算,当步数增多时,误差会因积累而越来越大。因此欧拉格式一般不用于实际计算。为提高精度,需要在欧拉格式的基础上进行改进。采用区间两端的函数

3、值的平均值作为直线方程的斜率,改进欧拉法的精度为二阶。二、基本知识与内容(一)Euler格式:对于方程组容易想到,可从y’着手解决它的数值计算问题。由数值微分向前差商公式得方程实际上给出故得。再由导出差分格式:称为Euler格式(二)改进的Euler格式:利用数值积分法导出新的差分格式,在此基础上提高Euler格式的阶数。对Euler公式的积分运算应用梯形积分公式,可得到如下求解公式:该公式称为梯形公式。Euler法是显示算法,计算量小,但精度低;梯形法是隐式算法,须借助迭代求解,计算量大,但提高了精度。综合使用上述两

4、种算法,可得到如下改进的Euler格式,又称Euler预测——校正公式:三、问题求解方法及步骤取h=0.1,分别用Euler格式、改进的Euler格式和Runge-Kutta格式解微分方程解:根据Euler公式根据改进Euler格式可得到用Euler法的程序如下:%欧拉方法求节点处的近似解x(1)=0;y(1)=1;h=0.1;forn=1:10y(n+1)=y(n)+h*(y(n)-2*x(n)/y(n));x(n+1)=x(n)+h;endplot(x,y,'r')%用红色线作近似解的图像holdonu=0:0.1

5、:1;v=sqrt(1+2*u);plot(u,v,'b')%用蓝色线作精确解的图像%作近似解得散点图,并标出点的坐标forn=1:11plot(x(n),y(n),'ro')%数值解的坐标点str=['(',num2str(x(n)),',',num2str(y(n)),')'];text(x(n),y(n)+0.2,str)%标出数值解坐标点endholdoff%改进的欧拉方法求节点处的近似解x(1)=0;y(1)=1;h=0.1;forn=1:10yp=y(n)+h*(y(n)-2*x(n)/y(n));%显式欧

6、拉法求预测值x(n+1)=x(n)+h;yc=y(n)+h*(yp-2*x(n+1)/yp);%隐式欧拉法求校正值y(n+1)=(yp+yc)/2;%求平均值endplot(x,y,'r')%用红色线作近似解的图像holdonu=0:0.1:1;v=sqrt(1+2*u);plot(u,v,'b')%用蓝色线作精确解的图像%作近似解得散点图,并标出点的坐标forn=1:11plot(x(n),y(n),'ro')%数值解的坐标点str=['(',num2str(x(n)),',',num2str(y(n)),')'];

7、text(x(n)-0.1,y(n)+0.1,str)%标出数值解坐标点endholdoff欧拉方法改进后的欧拉法四、计算结果分析Euler法改进后的Euler法显然改进型Euler方法的精度更高。

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。