常微分方程欧拉折线法的Picard迭代改进.pdf
《常微分方程欧拉折线法的Picard迭代改进.pdf》由会员分享,可在线阅读,更多相关《常微分方程欧拉折线法的Picard迭代改进.pdf(7页珍藏版)》请在咨信网上搜索。
1、Advances in Applied Mathematics 应用数学进展应用数学进展,2023,12(8),3763-3769 Published Online August 2023 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2023.128370 文章引用文章引用:孙波,熊常鹏.常微分方程欧拉折线法的 Picard 迭代改进J.应用数学进展,2023,12(8):3763-3769.DOI:10.12677/aam.2023.128370 常微分方程欧拉折线法的常微分方程欧拉折线法的Pi
2、card迭代改进迭代改进 孙孙 波波,熊常鹏熊常鹏 湖南科技学院理学院,湖南 永州 收稿日期:2023年7月26日;录用日期:2023年8月18日;发布日期:2023年8月28日 摘摘 要要 考虑常微分方程数值解法,对欧拉折线法加以改进。主要方法是考虑常微分方程数值解法,对欧拉折线法加以改进。主要方法是,对每个节点处局部线性近似解进行一对每个节点处局部线性近似解进行一次次Picard迭代修正,用局部修正解估算下一节点处解值。证明了算法局部截断误差比欧拉折线法高一阶,迭代修正,用局部修正解估算下一节点处解值。证明了算法局部截断误差比欧拉折线法高一阶,分析了算法的收敛性与稳定性。通过方程实例数值计
3、算进一步验证了本文算法所得近似解比欧拉折线法分析了算法的收敛性与稳定性。通过方程实例数值计算进一步验证了本文算法所得近似解比欧拉折线法近似解更接近精确解近似解更接近精确解。关键词关键词 常微分方程,欧拉折线法,常微分方程,欧拉折线法,Picard迭代,修正迭代,修正 Improvement of Euler Method for Ordinary Differential Equations by Picard Iteration Bo Sun,Changpeng Xiong College of Science,Hunan University of Science and Engineer
4、ing,Yongzhou Hunan Received:Jul.26th,2023;accepted:Aug.18th,2023;published:Aug.28th,2023 Abstract In this paper we study numerical algorithms for ordinary differential equations,and improve the Euler method.We refine local linear approximate solution at each node by Picard iteration,and then estimat
5、e the numerical solution at next node by the refined local approximate solution.The local truncation error of our algorithm is smaller than that of Euler method by one order.The convergence and stability are proven.Finally,we show that our algorithm gives approximate solu-tion closer to the exact so
6、lution than Euler method by an example.Keywords Ordinary Differential Equation,Euler Method,Picard Iteration,Correct 孙波,熊常鹏 DOI:10.12677/aam.2023.128370 3764 应用数学进展 Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(C
7、C BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 一阶常微分方程初值问题形如()00d,d.x xyf x yxyy=(1.1)为保证解的存在唯一性,(),f x y在含点()00,xy的某区域连续,且关于 y 局部 Lipschitz 连续。当(),f x y表达式较复杂时,难以用分离变量积分等办法求出解析解,只能通过数值方法计算近似解了。数值方法主要思想是,用差商近似代替微商,将方程差分化,离散化。经典方法主要有 Euler 方法和Runge-Kutta 方法。Euler 方法的简单体现形式为显式 Euler 公式,即 E
8、uler 折线法:()1,nnnnyyhf xy+=+,0,1,n=.其局部截断误差为()2O h,其中 h 为离散化步长。若用向后差商()()1nny xy xh+近似代替导数()1nyx+,还可得隐式或后退 Euler 公式:()111,nnnnyyhf xy+=+.其局部截断误差也是()2O h。具体计算时要先用显式 Euler 公式给出初步近似解,再由隐式公式进行迭代修正。Euler 方法的高级体现形式主要为梯形公式和预估校正公式。梯形公式可表述为()()111,2nnnnnnhyyf xyf xy+=+,其局部截断误差为()3O h。由于是隐式,具体计算时也需先用显式 Euler 公
9、式给出初步近似解,再由隐式公式进行迭代修正。预估校正公式即对显式 Euler 公式所给初步近似解进行一次迭代修正,可表述为()1,nnnnyyhf xy+=+,()()111,2nnnnnnhyyf xyf xy+=+.其局部截断误差为()3O h。Runge-Kutta算法思想是,用几个点处f 预估值加权平均作为当前节点处y估计值,以提高计算精度。为了进一步提高计算精度,还可利用当前节点1nx+前面多个节点1,nnn rxxx处的数值结果1,nnn ryyy估算1ny+,即线性多步法。据我们所查阅数值分析教材和参考书所述(1-6),一阶常微分方程初值问题经典数值解法主要就是这些。期刊文献上也
10、发表了一些探索性工作。如 2011 年,钟巍提出的一种数值解法,先将方程变形为(),ddf x yyxyy=,Open AccessOpen Access孙波,熊常鹏 DOI:10.12677/aam.2023.128370 3765 应用数学进展 再根据右端在当前节点值的正负给出下一节点数值解。文章举例说明了算法的近似效果,但未给出收敛性证明和误差估计(7)。现代人工神经网络理论和应用的发展也丰富了微分方程数值解法,通过训练带权重参数的试验解(trial solution),使其逐渐适应方程,这便是微分方程的神经网络数值解法(8)。应该说,常微分方程人工神经网络数值解法具有一定开创性。202
11、2 年,孙波、陈慧雄对一阶常微分方程初值问题解的存在唯一性 Picard 迭代证明法加以改进,提高了迭代收敛速度9。微分方程初值问题(1.1)可转化为等价积分方程()00,dxxyyf x yx+=.(1.2)通过构造迭代函数序列()()()010,dnnxxxyf xxx+=(1.3)可推导解的存在性。一般教科书拿()00 xy起始迭代,只要(),f x y在()00,xy的某矩形邻域内连续,且关于 y 李普希兹连续,则()nx在0 x的某邻域内一致收敛,其极限函数()x便是积分方程(1.2)即初值问题(1.1)的解。Picard 迭代的作用是改进局部近似解,每迭代一步使近似解更接近精确解。
12、孙波、陈慧雄将迭代起始函数改为()()()00000,xyf xyxx+,并证明了所得迭代函数列也一致收敛,且比以()00 xy起始迭代生成的函数列更快收敛到精确解。考虑到 Picard 迭代具有修正局部近似解作用,而欧拉折线法实质上就是分段局部线性近似求解。若将每个节点处局部线性近似解用 Picard 迭代修正一次,应更接近精确解。基于这一考虑,我们用 Picard迭代对欧拉折线法每个节点处局部线性近似解修正一次。第二节陈述算法,证明其收敛性与稳定性,并给出误差估计。第三节通过计算实例说明算法效果,并和欧拉折线法加以比较。2.基于基于 Picard 迭代修正的数值解法迭代修正的数值解法 本文
13、算法可表述如下:第一步:将初始点0 x处一阶近似解()()0000,yyf xyxx=+代入(1.2)或(1.3)右边,得修正的局部解()()()000000,dxxyyf x yf xyxxx+=+.(2.1)然后用(2.1)估算节点10 xxh=+处的解值:()()()10100000,dxxyyf x yf xyxxx+=.第二步:对1x处线性近似:()()1111,yyf x yxx=+进行 Picard 修正:()()()111111,dxxyyf x yf x yxxx+=+,并用它估算2xx=处的解值。孙波,熊常鹏 DOI:10.12677/aam.2023.128370 376
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微分方程 折线 Picard 改进
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。