基于STM的红外遥控和触屏控制的MP播放器剖析.ppt

一种改进的求解分数阶微分方程的预估-校正法及其应用 指导教师：虞继敏教授 项目成员：王臣 武卫红 匡春华 参赛时间：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时分数阶