三角指数拟合runge-kutta方法-runge - kutta method for triangular exponential fitting.docxVIP

三角指数拟合runge-kutta方法-runge - kutta method for triangular exponential fitting.docx

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
三角指数拟合runge-kutta方法-runge - kutta method for triangular exponential fitting

1绪论在微生物、电子网络、热核反应、航空航天以及化学动力学等许多重要科学技术领域和实际问题中,经常出现一类可用间歇性震荡或者刚性的常微分系统描述的物理化学过程。在这些问题中,往往会有一些可预先估计的特殊条件,比如震荡的频率、震荡的幅度及震荡的大体走向等。针对这些特殊的条件,在数值求解领域中出现了一类称为函数拟合的数值方法,这类方法首先应用于线性多步法CBettis[1]、Gautschi[2]),之后又在Runge-Kutta方法中得到了很大的发展。关于应用于常微分方程初值问题??y??f(x,y),C1.1)??y(x0)??y0的函数拟合Runge-Kutta方法C记为FRK方法),Ozawa在[3]中对包括FRK方法的构造、FRK方法系数的确定性、FRK方法的阶、提高阶的办法等一般形式的FRK方法可行性理论进行了详细的分析。ii?1在构造一个FRK方法之前,首先要确定一组线性无关的函数集合{u}k,称之为基础函数。而如果一个常微分方程初值问题的解可由这一组基础函数线性表出,那么应用此FRK方法所获得的解是精确的。另外,Hoang、Sidje和Cong[4]又在Ozawa的基础上,将FRK方法与配置条件、配置解及配置方法的理论结合起来讨论,指出了FRK方法其实就是配置Runge-Kutta方法在基础函数上的扩展和延伸。k然而在分析构造方法的过程中,Ozawa直接将方法的系数定位于与x无关,H没有指出具体的原因,也没有给出达到这样效果的方法。Sidje、Hoang和Cong在[5]中对这一问题做了补充说明,他们给出了可分的FRK方法的概念,即要求构造出来的FRK方法的系数与x是无关的,并且指出通过对函数集合的选择可使方法具有可分性。Sidje和Hoang又在[6]中对基础函数的选取做了进一步的分析,指出若将函数集合扩充为一个函数空间H,即H??span{{ui}i?1?{1}},那么FRK方法为独立的充ii?1分必要条件是H??H,其中H??span{{u}k?{1}}。ii?1基于上式,再区别于基础函数{u}k的不同,我们可以将FRK方法C此后的FRK方法若不作特殊说明均指可分的FRK方法)分为三大类:ii?1第一类是取{u(x)}k??{x,?,xk},称为多项式拟合,也就是经典的配置Runge-Kutta方法。第二类是取{u(x)}2r??{sin(??x),cos(??x),?,sin(r??x),cos(r??x)},或者取ii?1ii?1{u(x)}2r?1??{sin(??x),cos(??x),?,sin(r??x),cos(r??x)}?{x}等,其中????R,称为sine-cosine拟合或者三角函数拟合。Sidje、Hoang和Cong[5]对一般化的三角函数拟合方法进行了详细的研究,Calvoetal[7]在基于对方法的阶的提高的目的上,对偶数三角函数拟合方法进行了进一步的研究。ii?1第三类是取{u(x)}2rjjj?1??{xjexp(??x),xjexp(???x)}r,其中????R,称为exponetial拟合或者指数拟合。指数拟合方法在函数拟合方法的应用中最为广泛,分支也最为细致,从对一般化的指数拟合方法开始,如Simos[8],VandenBergheetal[9-10],VandenBerghe、DeMeyer和VanDaele[11]等,到对一些特殊情况或者特殊性能的指数拟合方法,如Vigo-Aguiar、Ramos[12],Simos、Vigo-Aguiar[13],VandeVyver[14],Calvoetal[15]等,都有很多研究成果。另外也有一些学者对上述的各种函数集合加以改动或者组合,构造出了一些新的方法,例如VandenBerghe、Ixaru和VanDaele[11]将多项式拟合与指数拟合相结合,ii?1即取函数集合{u(x)}k??{exp(??x),exp(???x),x,?,xk?2},????RColeman、Duxbury[16]k则是将多项式拟合与三角函数拟合相结合,即取函数集合{ui(x)}i?1??{sin(??x),cos(??x),x,?,xk?2},????R等。关于FRK方法CA,b,c)的系数,Ozawa在[3]中提出了向量c的各个分量必须是两两互异的要求,以此来保证方法系数的确定性。在更细致的研究中,越来越多的学者都不约而同的选择了对称的配置点,例如Gautschi[2],Calvoetal[7],VandenBergheetal[9]等,而Sidje、Hoang和Cong又在[5]指出了任意的二阶FRK方法如果使用了对称的配置点都将是A-定的C0v=pi,v???h,????R),这在一定程度上说明了对称配置点对定性的作用。同时一些学

您可能关注的文档

文档评论(0)

peili2018 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档