- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
一种改进的求解分数阶微分方程的预估-校正法及其应用 指导教师:虞继敏教授 项目成员:王臣 武卫红 匡春华 参赛时间:2014.11 1.引言 分数阶微积分理论的研究己有 300 多年, 但由于缺乏有效的应用背景,长期以来发展缓慢。在过去的近3个世纪里,分数阶微积分理论主要应用在纯数学领域,且其存在性也备受争议。1983 年 Mandelbort 首次报道了自然界存在大量分数维数的事实,此后分数阶微积分理论才得到真正的发展,并成为当前非线性科学的研究热点。因为分数阶系统可以更精准地描述出系统的特性,而且能真正地体现出系统的实际物理特性。在整数阶混沌系统的研究基础上, 众多学者将分数阶微分算子引入系统中, 发现当阶数为分数时, 系统仍表现出复杂的混沌行为由于分数阶混沌系统具有丰富的动力学特性和潜在的应用价值, 对其动力学特性及其应用研究引起了人们的广泛关注。 而Kai Diethelm等在2003年提出的预估-校正算法是目前比较成熟的求解分数阶微分方程数值求解方法。本文的工作主要就是对Kai Diethelm提出的预估-校正法做出改进,给出详细的证明过程,在第3节,给出一个数值例子来验证改进的预估-校正法在误差上的优势,在第4节,分别用改进的预估-校正法和Kai Diethelm的预估-校正法去寻找出现混沌现象的分数阶Lorenz系统的最小维数。 1.改进的预估-校正法 首先分数阶微积分的定义不尽相同,其中主流的定义有三种,分别Riemann- Liouville定义(RL),Grunwald-Letnikov定义(GL)和Caputo定义,本文采用的是Caputo定义,其表达形式如下 3.数值例子 在本节中,我们将给出一个具体的数值例子,计算并分析两种方法的误差,来证明改进的预估-校正法的精确性。数值例子如下: 其中有 且 。这个具有初始值问题的方程存在解析解,其解为 。我们选取一个正整数 ,令 。 我们有如下定义: ,为在 时的解析解,同理 为在 时用预估-校正法求得的数值解, 为 时用改进的预估-校正法求得的数值解。且两种方法求得的数值解的误差定义为: ,和 其中步长 的取值分别为1/100,1/200,1/400,1/800,1/1600和1/3200。下面表1是在不同的步长时两种方法的误差图。图1是 表1 步长不同时两种方法的误差绝对值 图1 当阶数为0.75,h=1/1000,N=900时两种方法误差的绝对值 在图1中,参数为 。可以看到,在初始值附近,两种方法的误差都是比较高的,之后随着步长增加都有下降。之后,预估-校正法的误差随步长增加增幅较大,而改进的预估-校正法的误差增幅较小,最后趋于稳定。 4.应用例子 在本节中,我们将采用 改进的预估-校正算法求解分数阶Lorenz混沌系统, 分析当阶数和系统参数变化时系统相图的复杂性变化,并与采用预估-校正法求解后得到的相图进行对比,以验证改进在的预估-校正法在求解上的优势。 广义的分数阶Lorenz系统如下 其中我们假设 参数设置为 系统的维数定义为 。2003年,Ilia Grigorenko在期刊《PHYSICAL REVIEW LETTERS》上发表文章《Chaotic Dynamics of the Fractional Lorenz System》,预测对于Lorenz系统,出现混沌现象的最小维数约为2.91 。如果假设 ,则有 。我们的工作就是,分别使用改进的预估-校正法和预估校正法对分数阶Lorenz系统进行仿真,并寻找出现混沌现象的最小维数。我们研究发现,使用预估-校正法找到的最小维数约为0.9935*3,而采用改进的预估-校正法,找到的最小维数约为0.9753*3。 当参数为 ,使用预估校正法求解分数阶Lorenz系统的数值解,其相图可表现出混沌现象;当 时,求解分数阶Lorenz系统后,得到的相图则表现不出混沌现象。图2和图3为阶数为0.9935和0.9934时分数阶Lorenz系统的相图。 图2 阶数为0.9935时分数阶
您可能关注的文档
- “新媒体”搭建青年成长新平台.pptx
- 基于net病历管理系统剖析.doc
- 集成信号发生器剖析.ppt
- 基于net的网上书店系统的设计与实现新剖析.doc
- 基于NGINX的代理服务器搭建毕业论文河源职业技术学院剖析.doc
- 集成运放的基础知识剖析.ppt
- 基于NiosII的数字钟设计剖析.doc
- 集成运放的设计与应用剖析.ppt
- “信息技术与学科教学有效整合的研究”课题工作汇报.ppt
- 集成运放电路及其应用剖析.ppt
- 运动风格趋势:轻量化户外 2025 S S.pdf
- 数字乡村建设村村通广播应用方案.pdf
- 高盛-日本市场策略 Japan Strategy Flash 1Q325 earnings summary Double-digit (%) profit growth as results largely surprised upwards 0820 2024.pdf
- CEEM《全球智库半月谈》(总第286期).pdf
- 智慧物流概论微课课件第4章智慧物流信息平台.pdf
- 数字乡村公共安全综合平台项目解决方案.pdf
- 智慧物流概论微课课件第7章智慧物流包装.pdf
- 【案例】”啤酒之王“的智慧物流探索之路.pdf
- 数字连江城市数字基座(一期)和数字乡村示范点建设项目.pdf
- 智慧物流概论微课课件第1章智慧物流概述.pdf
文档评论(0)