词条 | 分数阶微积分 |
释义 | 分数阶微积分是微积分的一个分支,它对函数进行分数阶微分积分,如对函数求1/2阶导数。 例如: 对x^n求1/2阶导数:1 简单介绍1简介现在关于分数阶导数研究论文每年约1000篇,且正在快速增长。分数阶微积分理论与应用的交流与学术会议日益频繁。每年都有比较大型的国际会议,小型会议越来越多(学术关注度详见图1)。 分数阶导数主要具有以下优势:1.分数阶导数具有全局相关能较好地体现系统函数发展的历史依赖过程;而整数阶导数具有局部性,不适合描述有历史依赖过程。 2.分数阶导数模型克服了经典整数阶微分模型理论与实验结果吻合不好的严重缺点,使用较少几个参数就可获得很好的效果。3.在描述复杂物理力学问题时,与非线性模型比较,分数阶模型的物理意义更清晰,表述更简洁。 图1 2.2数值方法的研究现状近年来分数阶微积分被广泛的应用于反常扩散、信号处理与控制、流体力学、图像处理、软物质研究、地震分析、粘弹性阻尼器、电力分形网络、分数阶正弦振荡器、分形理论、分数阶PID控制器设计。但是由于分数阶微积分具有历史依赖性与全域相关性,增加了分数阶导数方程的数值计算复杂性。 在数值算法方面主要存在的问题有:(1)长时间历程问题一直没有找到一个满意的解决途径,在数值模拟中,随着时间历程的增加计算量成指数增长。同时一些学者提出的短期记忆方法只对很少一些情况有效,并不具有普适性。因而长时间历程问题的解决任重道远。(2)在原有算法基础上开发出时间-空间混合的分数阶导数方程的算法和软件。一种数学工具要在工程中有广泛的应用,那么就必须有成熟的算法与软件,像有限元的计算模拟软件就有很多,所以有限元才能在工程界有如此广泛的应用。(3)分数阶导数的定义还不完善,现在分数阶导数的定义有多种,至今还没有一个完善到大多数学者能够接受的定义。 现阶段,分数阶导数方程的数值算法主要包括:(1)有限差分法:显示格式,隐式格式,Crank-Nicholson格式,预估校正算法,线性算法等;(2)级数逼近法:变分迭代法,Adomian分解法,同伦摄动法,通论分析法,微分转换法等;(3)有限元法;(4)无网格方法;(5)一些新的算法:矩阵转化法,外推法等。 以上这些数值算法各有优缺点,不同的条件与方程适用于不同的算法,这需要我们对各种方法比较熟悉,能够比较灵活的应用。否则很容易得到错误的计算结果。另外一个难点是在数值计算中使用哪一个分数阶导数的定义,这就涉及到定义的选择问题。根据大量文献参考后,一般在时间分数阶导数的计算中一般使用Caputo定义,在空间分数阶导数方程的数值计算中较多的使用Riemann-Liouville定义和级数定义。 3. 分数阶导数的定义关于分数阶导数的定义,许多数学家各自从不同角度入手, 给分数阶导数分别以不同的定义。其定义的合理性与科学性已在实践中得以检验。这个数学分支的发展已在实际问题中, 得到了广泛的应用。本文这部分重点将分析各种不同的定义,也说明各种定义之间的区别与联系。为了区分整数阶导数的表示形式,对于分数阶的导数,本文引入新的记号(下文Riemann-Liouville积分中表示阶积分,在此申明以防混淆)。 3.1Riemann-Liouville定义及性质 下面先阐述在研究中应用得比较广泛的一种分数阶导数的定义:Riemann-Liouville分数阶导数。在定义分数阶导数之前,先来阐述下Riemann-Liouville分数阶微积分。 3.1.1Riemann-Liouville定义 定义 设在上逐段连续,且在的任何有限子区间上可积,对,称 (1) 为函数的阶Riemann-Liouville积分(简称R-L积分),并且记为。其中为Gamma函数。 结合上面的阶Riemann-Liouville分数阶积分的定义以及经典微积分中的整数阶微积分可以给出如下的阶Riemann-Liouville分数阶微分的定义: 定义 设,,是大于或等于的最小正整数(),记。则称 (2) 为函数的阶Riemann-Liouville微分。 应用定义1可得阶Riemann-Liouville微分如下: 3.1.2Riemann-Liouville分数阶导数的性质 设,是满足定义1的函数,为任一常数,为分数,,则有文献【1】得到以下两条线性的性质: 性质1:。 性质2:。 3.2分数阶导数的级数定义 分数阶导数的级数定义又称为Grunwald-Letnikov定义。我们先来看整数阶导数的定义。 一阶导数的定义: (4) 二阶导数的定义: 通过选择相同变量h,令, 则(5)式等价于 (6) 那么对于阶导数来说就是以下的(7)式: 从(7)式整数阶的导数我们可以从形式上得到分数阶导数的级数定义。从形式上式(7)中的可以推广到非整数,组合数可以用Gamma函数来描述。而求和的上限(并非整数)也变成是(其中和分别是微分的上极限和下极限)。所以我们得到了用级数定义的分数阶导数(如下(8)式),又称the Grunwald-Letnikov fractional derivative。 定义3 在具有阶连续导数,并且至少取的条件下, Riemann-Liouville定义与Grunwald-Letnikov定义等价。所以由(7)到(8)的推广也是合理的。 3.3分数阶导数的Caputo定义 定义4 对于正的非整数(在阶数为负实数时,Caputo定义与Riemann-Liouville定义等价),,,是大于或等于的最小正整数()。则称 为函数的阶Caputo分数导数。 3.4三种分数阶导数定义的关系 Riemann-Liouville定义是Grunwald-Letnikov定义的扩充,其应用范围也就更广泛。与Grunwald-Letnikov定义扩展到Riemann-Liouville定义的思维方式相似, Caputo定义也是对Grunwald-Letnikov定义的另一种改进。对于函数的正的非整数阶导数,先进行阶导数, 再进行阶积分。 Riemann-Liouville定义与Caputo定义都是对Grunwald-Letnikov定义的改进。在阶数为负实数和正整数时, 它们是等价的。由文献【15】分析可知,在条件:(1)函数有阶连续导数, 至少取。(2)之下, 它们也是等价的。否则, 它们不等价。引入Riemann-Liouville定义, 可以简化分数阶导数的计算;引入Caputo定义, 让其拉普拉斯变换式更简洁, 有利于分数阶微分方程的讨论。 现实应用中,具体使用哪一个分数阶导数的定义,还是要看具体的情况而定。 4. 分数导数的应用举例分数阶导数在很多领域都有应用,下面拿与生活联系比较紧密的气候研究、医学图像处理、地震分析为例进行进一步地阐述与说明。 4.1天气和气候的研究 我们都知道没有一天天气是一样的,而气候的预测也不可能提到日程上来研究。这说明天气和气候的研究是比较困难的。天气和气候虽然遵从流体力学规律,但是却显示出随机性,研究天气和气候之间的关系必须引入分数阶的导数和积分,从物理上讲不外乎说明天气和气候的随机程度是不相同的。为此提出气候的q (0 ≤q ≤1) 阶微商是天气。此时引入天气和气候之间的桥梁——分数阶导数,这为天气与气候的研究带来很大的方便。 4.2医学图像处理 医学图像一般是指为了清楚地看到病人内部的局部器官病变情况而通过一定的设备仪器得到的图片,例如CT、B超等图片。由于设备,技术等方面的原因,得到的医学图像有可能模糊不清。图像的不清晰对临床诊断带来很大的麻烦。所以要考虑怎样处理,可以得到更清晰的医学图像。 现在从分数阶微分基本定义出发,可以作用于二维医学图像的分数阶微分掩模,掩模可以根据对图像的需求进行增强。通过实验证明,这个方法可以有效完成对医学图像的处理,并且弥补了传统方法不能连续改变处理效果的缺点,是一种简单可行并且效果较好的图像增强方法。 所以说分数阶导数对医学图像的处理,帮助是很大的。 4.3地震奇异性分析 由文献【4】,我们知道传统的地震解释主要是观测地震资料的振幅及相位的变化, 而振幅往往并不能反映真实的地质情况。地震界面可能是岩性分界面也可能是岩性过渡带, 岩性过渡带的地震反射波是入射波的分数阶导数。 因此我们将分数阶导数引入地震属性计算中, 构建一种对波形敏感而对振幅变化不敏感的新属性——奇异性, 用以刻画反射界面的横向变化。 方法的基本原理是首先计算地震子波的不同分数阶导数,然后利用匹配追踪算法将地震数据分解成地震子波的不同分数阶导数,进而获得反射波同相轴的分数阶。对胜利油田某区块实际二维地震资料进行了试处理,结果表明分数阶导数剖面能很好地描述不整合面,反映实际界面的横向变化。 参考文献[1] 王小东.Riemann-Liouville分数阶微积分及其性质证明[D].太原理工大学硕士学位论文,2008. [2]邓伟华.分数阶微分方程的理论分析与数值计算[D].上海大学博士学位论文,2007. [3]晏祥玉,周激流.分数阶微积分在医学图像处理中的应用[J].成都信息工程学院学报, 2008, 23 (01):38~41. [4]宋建国,刘垒,李辉等.分数阶导数在地震奇异性分析中的应用[J].石油物探, 2009,48(01):72~75. [5]刘朝辉,张为民.含分数阶导数的粘弹性固体模型及其应用[J].株洲工学院学报,2002,16(04):23~25. [6]刘林超,黄学玉,闫启方.基于Lagrange变分法的分数导数粘弹性梁振动分析[J].噪声与振动控制,2007,(02):43~46. [7]胡卫兵,何建.粘弹性阻尼材料的分数导数模型的对比分析[J].西安建筑科技大学学报(自然科学版),2002,34(03):222~227. [8]王少伟.分数阶微积分理论在粘弹性流体力学及量子力学中的某些应用[D].山东大学博士学位论文,2007. [9]姚奎.分形函数与分数阶微积分-构造性方法的应用[D].浙江大学博士学位论文,2003. [10] 王振滨,曹广益,曾庆山.分数阶PID控制器及其数字实现[J].上海交通大学学报,2004,38(04):517~520. [11] 袁晓,张红雨,虞厥邦.分数导数与数字微分器设计[J].电子学报,2004, 32(10):1658~1665. [12] 孙洪广.分数阶微积分的研究[EB/OL] .(2010-07-30)[2010-08-03]. [13] 陈纪修.数学分析(第二版上册)[M].北京:高等教育出版社,2008:119~129. [14] Adam Loverro. Fractional Calculus: History, Definitions and Applications for the Engineer.Department of Aerospace and Mechanical Engineering University of Notre Dame,May 8, 2004. [15] 林孔容.关于分数阶导数的几种不同定义的分析与比较[J].闽江学院学报,2003,24(05):3~6. [16] 胡亦郑,刘发旺.一类分数阶控制系统的数值解法[J].厦门大学学报(自然科学版), 2005,44(03):313~317. [17] 刘式达,时少英,刘式适等.天气和气候之间的桥梁-分数阶导数[J].气象科技,2007,35(01):15~19. 首先 对x^n求1阶导数后为nx^(n-1) 2阶导数后为n(n-1)x^(n-2) ... 那么m<n时 m阶导数后为n(n-1)(n-2)..(n-m+1)x^(n-m) 也就是n!/(n-m)! x^(n-m) 由于阶乘的概念可以扩展: n!=Γ(n+1) 故对x^n求m阶导数后为Γ(n+1)/Γ(n-m+1)x^(n-m) 令m=1/2 则得 d(1/2)x^n/dx^(1/2)=Γ(n+1)/Γ(n+1/2)x^(n-1/2)=n!(n-1)!2^(2n-1)x^(n-1/2)/(2n-1)!π^(1/2) 特别x的1/2导数为2(x/π)^(1/2) 对于其他函数……由于任何连续f(x)都可以展开Taylor级数……因此,任何连续函数都存在1/2阶导数 例如sin x的1/2阶导数为sin(π/4+x) 引言分数阶导数,简单来讲就是对整数阶导数理论的拓展。例如,我们一般对某个性质较好的可导函数,可以求出它的一阶导数、二阶导数、……、阶导数。那么我们是否可以对函数求分数阶导数呢?比如导数。再如某个函数不满足求导条件,我们是否可以使用微积分理论对这个函数进行分析性质的研究呢?根据多方文献的参考得知,答案是肯定的,这也是分数阶导数产生的源动力。 早在1695年,Leibniz给L’Hospital写了一封信,问:“整数阶导数的概念能否自然地推广到非整数阶导数。”L’Hospital对这个问题感到很新奇,作为回信他反问了一个简单的问题:“如果求导的次数为,那么将会是怎样的情况呢?”在这一年的9月30号,Leibniz给L’Hospital回信写到:“这会导致悖论,不过总有一天会得到有用的结果。”这个特殊的日子1695年9月30号被认为分数阶微积分的诞生日。 自1695年,分数阶导数的研究已经经历了三百多年。但是早期分数阶导数的研究主要存在于理论数学领域。在很长的一段时间内,分数阶微积分的研究没有得到自然科学与工程科学研究人员的关注,基本上没有相关的应用文章发表。分数阶微积分的研究热潮是在二十世纪七十年代,主要原因是因为研究人员发现分形几何、幂律现象与记忆过程等相关现象或过程可以与分数阶微积分建立起密切的联系。分数阶微积分可以作为一种很好的描述与刻画手段。 |
随便看 |
|
百科全书收录4421916条中文百科知识,基本涵盖了大多数领域的百科知识,是一部内容开放、自由的电子版百科全书。