欧拉法的应用.doc

ID:59153298

大小:28.00 KB

页数:3页

时间:2020-09-11

欧拉法的应用.doc_第1页
欧拉法的应用.doc_第2页
欧拉法的应用.doc_第3页
资源描述:

《欧拉法的应用.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、欧拉法的应用1.问题背景在工程当中,我们所建立的方程一般都是微分方程,对于其求解非常重要,但是其求解又很难求得精确解,精确解也并不是我们所需要的。我们可以求得近似解。我们采用欧拉法来求解一阶常微分方程。2.数学模型一阶常微分方程如下式所示:3.算法及流程a.对于初值问题,先将其离散化,即把[a,b]区间n等分,得各离散节点,b.由上一点推出下一点的值,逐渐推导到n点。公式如下:c.最终可以计算出终点的近似值matlab程序为:clc;clear;symsxsymsyfun=input('请输入常微分函数右侧表达

2、式=');xinit=input('请输入初始点x0=');yint=input('请输入初始点y0=');xfinal=input('请输入求值点xfianl=');n=input('请输入划分区间个数n=');f=inline(fun);h=(xfinal-xinit)/n;x=zeros(1,n+1);y=zeros(1,n+1);x(1)=xinit;y(1)=yintfori=1:nx(i+1)=x(i)+h;y(i+1)=y(i)+h*f(x(i),y(i));endT=[x',y']y(n+1)p

3、lot(x,y)1.计算结果及分析请输入常微分函数右侧表达式=y-2*x/y请输入初始点x0=0请输入初始点y0=1请输入求值点xfianl=1请输入划分区间个数n=10T=01.00000.10001.10000.20001.19180.30001.27740.40001.35820.50001.43510.60001.50900.70001.58030.80001.64980.90001.71781.00001.7848ans=1.7848从图形结果可以看出,红色的点是每一步计算出来的点。求解精度并不高,并

4、不能得到很好的解。这种算法还需要改进。所以后面出现了改进的欧拉法。

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

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

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

《欧拉法的应用.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、欧拉法的应用1.问题背景在工程当中,我们所建立的方程一般都是微分方程,对于其求解非常重要,但是其求解又很难求得精确解,精确解也并不是我们所需要的。我们可以求得近似解。我们采用欧拉法来求解一阶常微分方程。2.数学模型一阶常微分方程如下式所示:3.算法及流程a.对于初值问题,先将其离散化,即把[a,b]区间n等分,得各离散节点,b.由上一点推出下一点的值,逐渐推导到n点。公式如下:c.最终可以计算出终点的近似值matlab程序为:clc;clear;symsxsymsyfun=input('请输入常微分函数右侧表达

2、式=');xinit=input('请输入初始点x0=');yint=input('请输入初始点y0=');xfinal=input('请输入求值点xfianl=');n=input('请输入划分区间个数n=');f=inline(fun);h=(xfinal-xinit)/n;x=zeros(1,n+1);y=zeros(1,n+1);x(1)=xinit;y(1)=yintfori=1:nx(i+1)=x(i)+h;y(i+1)=y(i)+h*f(x(i),y(i));endT=[x',y']y(n+1)p

3、lot(x,y)1.计算结果及分析请输入常微分函数右侧表达式=y-2*x/y请输入初始点x0=0请输入初始点y0=1请输入求值点xfianl=1请输入划分区间个数n=10T=01.00000.10001.10000.20001.19180.30001.27740.40001.35820.50001.43510.60001.50900.70001.58030.80001.64980.90001.71781.00001.7848ans=1.7848从图形结果可以看出,红色的点是每一步计算出来的点。求解精度并不高,并

4、不能得到很好的解。这种算法还需要改进。所以后面出现了改进的欧拉法。

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