分数阶微积分的探讨 (2)

绪论...................................................1 1 分数阶微分的基本理论..............................................1 1.1 分数阶微积分..................................................2 1.2 分数阶微积分的定义............................................3 1.3 分数阶微积分的性质............................................5 1.4 各种定义之间的联系与区别......................................6 1.5 一些初等函数的分数阶微积分....................................8 1.6 分数阶微积分的物理意义.......................................10 1.7 分数阶微积分在自然中的存在...................................11 2 分数阶微积分的应用...............................................12 2.1 医学图像处理.................................................12 2.2 天气和气候的研究.............................................13 2.3 地震奇异性分析...............................................14 参考文献...........................................................15 致谢...............................................................16

I

湖北大学本科毕业论文（设计） 分数阶微积分及其应用

摘 要

分数阶微积分作为整数阶微积分的推广，其概念早已提出，近300年来，分数阶微积分这一重要数学分支渐成体系，它是研究分形分析的重要工具被应用于许多工程计算中。本文给出了分数阶微积分的一些性质及其推导过程，并给出一些初等函数的分数阶微积分,及其应用。

【关键词】分数阶微积分 分数阶微分 分数阶积分 图像增强 模板 应用

II

湖北大学本科毕业论文（设计） Fractional calculus and its applications

Abstract

Fractional Calculus as extention of integral calculus, its concept has long been proposed, for nearly 300 years, fractional calculus of this important branch of mathematics that had gradually become the system , it is the study of fractal analysis tools are used in many engineering

calculations .in his paper, some properties of the fractional calculus and the derivation process of the fractional calculus are given, Besides some elementary functions of fractional calculus and its applications.

【Key Words】Fractional Calculus Fractional derivatives Fractional integrals image

enhancement applications

III

湖北大学本科毕业论文（设计） 绪论

分数阶微积分是微积分的一个分支，它对函数进行分数阶微分积分，如对函数求1/2阶导数。分数阶微积分(Fractional Calculus)是相对传统意义上的整数阶微积分提出来的。普通的微积分运算，如一阶微分、二阶微分、一阶积分、高阶积分等都是在积分运算次数为整数情况下的微积分运算，而分数阶微分，顾名思义，就是将通常意义下的整数阶微积分运算推广到运算阶次为分数的情况。例如：对x^n求1/2阶导数; 分数阶微积分就可以看作是整数阶微积分运算的推广。值得注意的是，其实将上述这种非整数阶的微积分称之为“分数阶微积分”或“分数阶演算”是一种不严格的命名，其中的“分数”是一个不准确的归类。因为指数v可以被推广到有理分数、无理数甚至复数，所以从严格数学意义上讲，它应该被称为“非整数阶微积分”或“非整数阶演算”。但是由于历史的原因，“分数阶微积分”或“分数阶演算”的命名己经成为习惯用法。

分数阶微积分的概念虽然早就提出来了，但是由于对分数阶微分方程的求解缺乏相应的数学工具，所以它在工程中的应用一直到上世纪后期才有研究。实际系统中有许多是分数阶的而不是整数阶的，之所以将它们作为整数阶系统来考虑是由于其复杂性。然而随着分数阶微积分理论的发展，对于分数阶微积分在实际系统中的应用开始了研究。

分数阶微积分是一个古老而新鲜的概念。早在整数阶微积分创立的初期，就有一些数学家，如L‘hospital、Leibniz等开始考虑它的含义。然而，由于缺乏应用背景支撑等多方面的原因，它长期以来并没有得到较多的关注和研究。随着自然科学和社会科学的发展、复杂工程应用需求的增加，尤其是20世纪七八十年代以来对分形和各种复杂系统的深入研究，分数阶微积分理论及其应用开始受到广泛关注。

进入21世纪以来，分数阶微积分建模方法和理论在高能物理、反常扩散、复杂粘弹性材料力学本构关系、系统控制、流变性、地球物理、生物医学工程、经济学等诸多领域有了若干非常成功的应用，凸显了其独特优势和不可代替性，其理论和应用研究在国际上已成为一个热点。近年来分数阶微积分被广泛的应用于反常扩散、信号处理与控制、流体力学、图像处理、软物质研究、地震分析、粘弹性阻尼器、电力分形网络、分数阶正弦振荡器、分形理论、分数阶PID控制器设计。但是由于分数阶微积分具有历史依赖性与全域相关性，增加了分数阶导数方程的数值计算复杂性。特别是在信息科学领域中，一些新颖的应用被相继地实现，如系统建模、曲线拟合、信号滤波、模式识别、图像边界提取、系统辨识、系统稳定性分析等等。

分数阶微积分的非局域性质，导致分数阶导数控制方程数值模拟的计算量和存储量随问题规模的增大而增加得比相应整数阶方程快得多，一些计算整数阶方程十分有效的数值方法对分数阶方程也完全失效。而且，目前大多数的分数阶微积分方程模型还是唯象模型，其内在的物理和力学机理还不是很清楚，有待进一步的深入研究。同时一些学者提出的短期记忆方法只对很少一些情况有效，并不具有普适性。因而长时间历程问题的解决任重道远。⑵在原有算法基础上开发出时间－空间混合的分数阶导数方程的算法和软件。一种数学工具要在工程中有广泛的应用，那么就必须有成熟的算法与软件，像有限元的计算模拟软件就有很多，所以有限元才能在工程界有如此广泛的应用。⑶分数阶导数的定义还不完善，现在分数阶导数的定义有多种，至今还没有一个完善到大多数学者能够接受的定义。

1 分数阶微分的基本理论

1.1 分数阶微分

目前通常使用的分数阶微积分记法是：

1

湖北大学本科毕业论文（设计） 如果对于自变量为 x的函数 y = f ( x)，自变量取值范围为( a , b )区间，在该区间上，对函数 f ( x )的v阶分数阶微积分记为：

vDfxfabvx

其中微积分阶次v是分数，若 v > 0,v ∈ R则称为分数阶微分。如果v< 0, 则称为分数

阶积分。

对于连续函数yft,依据整数阶导数的定义，它的一阶导数定义为： fftfthdflim （1.1） dthoh依据相同的定义，可以推出二阶导数的定义式：

ftfth2fft2lim

h0th lim1ftfthfthft2h h0hhhft2fthft2h （1.2） 2h limh0根据（1.1）和（1.2）我们得知：

ft3fth3ft2hft3hd3fft3lim

h0dth3更一般的，n阶导数的一般定义可以记为：

f其中

ndnf1tlimdtnh0hn1ftrh，

nrr0nr为二项式系数，

nrpnrnn1n2...nr1， !r现在，让我们考虑一般化分数阶在表达式中的情况：

fh1tph1ronrftrh

pr公式中：p,n是两个任意的整数。 显然对于pn,我们有 limfhh0ptfpdpftp，因为在这种情况下，开始的分

dt子中所有的系数在p之后都等于零。

让我们考虑p为负值的情况。为方便起见，我们定义为：

2

湖北大学本科毕业论文（设计） N!NttlimNtt1t2...tN

因此我们有

p1rrpp1...pr1pr !r如果用-p代替p可以写作：

p1t1tfd， ap1!其中p是一个正整数。

如果n是一个固定的数，那么当h0时，fh数，并且考虑到fhppt的极限趋近于零。当时，为了渠

到一个非零的极限，我们必须假设n，我们可以令h=(t-a)/n，其中a是一个实常

t的极限可能为无限或有限值，我们可以表示为：

limfhptaDtpft

h0我们可以推导出一般的表达式为：

aDft=limhpth0pronprftrh

p1t1=tfd ap1!1.2分数阶微积分的定义

当分数阶微积分阶次v取自然数n时(n ∈ Z)，即v = n时，表示为通常意义下的整数阶微分；v = − n时表示整数阶积分。由此可见，分数阶微积分是整数阶微积分的推广，整数阶微积分可以看成是分数阶微积分的特例。

在分数阶微积分理论发展过程中，出现了很多种函数的分数阶微积分的定义，如由整数阶微积分直接扩展而来的 Cauchy 积分公式，Grünwald-Letnikov分数阶微积分定义、Riemann-Liouville 分数阶微积分定义以及 Caputo 定义等。

关于分数阶导数的定义，许多数学家各自从不同角度入手，给分数阶导数分别以不同的定义。其定义的合理性与科学性已在实践中得以检验。这个数学分支的发展已在实际问题中，得到了广泛的应用。本文这部分重点将分析各种不同的定义，也说明各种定义之间的区别与联系。

1.2.1 Cauchy 积分公式

该公式由整数积分直接扩展而来

3

湖北大学本科毕业论文（设计） Dvftv1fd v12jct其中c为包含ft 单值和解析开区域的光滑曲线， Γ (v+ 1)为 Gamma 函数，它也是第二类欧拉积分，在分数阶微积分的运算中有着非常重要的作用。完全的函数是以

Euler极限的形式给出，如下面所示：

N!Nttlim Ntt1t2...tN常用的函数的形式是积分变换形式： tt1yyedy 01.2.2 Grunwald-Letnikov 分数阶微积分定义

aDptft=limhh0pronprftrh

1.2.3 Riemann-Liouville 分数阶微积分定义

11=tfd Dftapaptt0 << 1，且为初值，特别地，D左右侧的下标分别表示积分的下、上界限。 还可以定义出分数阶微分，假设分数阶 n − 1< ≤ n，则定义其分数阶微分为

ndnn1dDftaDtftnatdtndtntffn1

at1.2.4 Caputo 分数阶微分定义

tfm11d 0Dtfv1v0t其中mv，m 为整数，0 < v< 1。类似地，Caputo 分数阶积分定义为

ty1rd，v<0 0Dtyt1vv0t对很广一类实际函数来说，Grunwald-Letnikov 分数阶微积分定义及 Riemann-Liouville 分数阶微积分定义是完全等效的，

1.3 分数阶微积分的性质

4

湖北大学本科毕业论文（设计） 对上述所有分数阶微积分的定义形式，都有如下几条共同性质: (1) 解析函数ft的分数阶0Dtft导数对t和都是解析的。

(2) =n为整数时，分数阶微分与整数阶微分与整数阶微分的值完全一样，且

0D0tft=ft。

x1ux(3) aDtt，1xum为整数时，Dxtu=0 .

xu1xut(4) 分数阶微积分算子是线性的，亦即对任意常数 x,y，有

Dtxftygtx0Dtfty0Dtgt .

0 （5） 分数阶微积分算子满足交换律

 0DtDftDDftDft . 0t0t0t0t（6） 当N>0且N为整数，u为分数时

0DtNuft0DtN0Dtuft .

（7）u0,1，且uN，NN*，则

0Dtu0Dtft0DtN0ftQNt,N .

证明：

（2）当=n为整数时，明显可得分数阶微分与整数阶微分与整数阶微分的值完全一样；0Dtft相当于整数阶微积分，也就相当于没有做运算，所以依然等于ft，（也

0可根据定义式带入求解）.

（3） Dxtu1dmxdt1dmxdtmtmx1tttz0mx1tzudz ，

mtty0yudy，

5

湖北大学本科毕业论文（设计） m1dmxuBu1,mx tmxdt y(u1)uxt m1xm;z .

t(ux1) （4）由0Dtft的微分算子是线性的，即有

0Dtxftygt=0Dtxft0Dtygt

=x0Dtfty0Dtgt . （5）由定义易得知.

（6）令vmu且m为大于v的最小整数则

v 0DtNuft0DtNmD0tft

0vDtN0Dtuft0DtN0Dtm0Dtvft0DtNm0Dtft ，

即证.

（7） 0Dtu0Dtft0Dtu0DtN0DtNft

0Dtu0DtN0DtNft0DtuQNt,NN

0Dtu0Dt 0DtNNNuDftD0t0tQNt,

0DtN0DtNft0DtuQNt,

0DtNft0Dtu0DtQNt,00DtNftQNt,N.

1.4 各种定义之间的联系与区别

从最初提出分数阶微积分的概念开始，许多数学家从不同的分析角度入手，结合实际问题的应用，分别给出了分数阶微积分的不同定义，这就是目前对于分数阶微积分存在多个运算定义的原因。在这众多的定义之间，存在着一定的区别和联系，或者在某些条件下它们之间可以进行相互转换。目前最为常用的分数阶微积分定义应当是 Riemann-Liouville 分数阶微积分定义。

用Grunwald-Letnikov分数阶导数定义分数阶向后分差2是不方便的。计算得到的表达式中因为有整数的原因，所以表面上看是比较好一些；单对于非整数的情况显然不适合计

6

湖北大学本科毕业论文（设计） 算。因为其表达式中所要求的函数必须是m次以上连续可微，其所能表达的分数阶导数的定义也包含在Riemann-Liouville 分数阶微积分定时的使用情况之中。因此，

dpDftatdt m1ttampfd

k0mtfkata1mpm1tf() pk1pk1apk aDtpft,mpm1

因此，如果我们考虑函数ft，当t>0时，才在m+1阶可微且连续，那么Grunwald-Letnikov分数阶导数定义：

fptndnf1tnlimh0hndtpnpr1ftrh

nrr0nrp1t1又因为 aDft=limhfd ftrh=p1!ath0ro由 n1n(n1)(n2)...21n!，可得

1 aDtft=limnh0hp1r0nrnr1=ftrhpt1fd

at由上式可知，对于整数阶微积分运算，即微积分阶次v = n ∈ Z时，Riemann-Liouville(R-L) 定义与 Grünwald-Letnikov(G-L) 定义计算所得到的结果一致。对于运算阶次为分数阶 v ( n − 1 < v < n)时，在假设函数 f ( x )满足具有 m+ 1阶连续导数，并且m 至少取 m = n− 1的条件下，Grünwald-Letnikov 定义与 Riemann-Liouville 定义等价。但他们在相当广的一类函数中还是相等的，Riemann-Liouville 分数阶微积分定义是Grunwald-Letnikov分数阶微积分定义的推广

而Caputo 分数阶微分定义：

tfn1d aDtft=n1nat（n1,n1n,t），式子中n为大于a的最小正整数：fn为函数

f的n阶导数。将其与Riemann-Liouville分数阶定义相比，Caputo定义将函数的整数

阶导数的运算转移到了积分内部，改为对的求导，运用分部积分公式，可以转变为如下形式：

tfd1aDtftn1

natn7

湖北大学本科毕业论文（设计） fnatan1na1nn1tfd n1at其式子的右端第一项为整数阶导数的初始条件，这在某一方面就增加了定义式的适用性。因此它在物理学的流体学和力学中对于问题的求解相比其他定义而言比较方便。当

n时，及分数阶导数靠近整数阶导数：

limaDtftlimnfnnatan1tna1nn1tfd n1atfnafn1d

a fnt （n为自然数）

比较Caputo定义式子与Riemann-Liouville定义式，两式的区别在于对微分与积分顺序

的不同，前面是积分在内而后面则是积分在外。从对函数的要求来看，Caputo定义式的更严格，它需要函数n阶可微。与Grunwald-Letnikov定义扩展到Riemann-Liouville定义的思维方式相似，Caputo定义也是对Grunwald-Letnikov定义的另一种改进。对于函数的正的非整数阶导数，先进行阶导数，再进行阶积分。前者对常数的求导是有界的(为 0)，而后者求导是无界的，Caputo 定义更适用于分数阶微分方程初值问题的描述。

Caputo定义式最主要的优点是分数阶微分方程的初始条件，对Caputo定义采取了与整数阶微分方程一样的形式，加强了在端点值时的。

1.5 一些初等函数的分数阶微积分

（i）常数C 的分数阶微积分：

11mp1dmdmCDCxtCdtxmp 0mmmpdx0mpdxmppxxCxp ,(m1pm)

1p(ii) 指数函数eλx的分数阶微积分:

1p1Dextetdtpex,(0) ppxxx(iii) (A)正弦函数sinx的分数阶微积分:

8

湖北大学本科毕业论文（设计） 1p1pDsinxxtsintdtsinxp,(0,p1) p2pxx (B)余弦函数cosx的分数阶微积分:

1p1pDcosxxtcostdtcosxp,(0,p1) p2pxx明显，当上面几个式子中的p为整数时，其积分结果可验证与整数阶微积分结果相同，只是积分的下界不一致。表1-5-1是一些函数半阶导数和半阶积分

半阶积分 函数 半阶微分 2Cx C 1 xC1 x x 24x 3320 x  2xx 2 xx(1x)arctan(1x)arctanx x 1x1 x1arctanarctanx 1x xx 322arcsinh(x) x(1x)2arcsin(x) x(1x)1 1x1xxarcsinh(x) x(1x)1 1x1xxarcsin(x)x(1x)32 9

湖北大学本科毕业论文（设计） 2x(ln(4x)2) lnx 表 1-5-1 ln(4x) x1.6 分数阶微积分的物理意义

对于整数阶微积分理论而言，有着明确的物理意义与几何意释，然而对于分数阶微积分，目前仍然没能找到一种明确的物理与几何解释。

分数阶微积分缺少物理与几何意义的解释这个问题从一开始就被人们所注意，直到 1974 年在美国 New Haven 召开的第一届分数阶微积分国际会议上再次被人们提出，随后的多次国际会议也被提出进行讨论，仍然没能得到满意的结果。许多科学家分别从不同的角度，对该问题进行了深入探讨研究，取得了众多的研究结果，对于分数阶微积分的物理与几何意义做出了部分解释。但是多数解释是从一些分数阶微积分特例中得到，普适性不强，并不能很好地完全解释分数阶微积分的物理意义与几何意义。该问题仍然是一个亟待解决的富有挑战性的问题。

以下将介绍一种关于分数阶微积分物理意义与几何意释，该观点是Igor Podlubny 教授在 2001 年提出的。其解释是基于多种分数阶微积分定义所做出，其中也包括了 Riemann-Liouville 分数阶微积分的定义，这里对该观点进行简单阐述。

首先介绍两个不同的时间的概念：个人时间τ 和宇宙时间T 。宇宙时间是均匀的、等间隔的流逝的时间，也叫做绝对时间，记为T ，个人时间τ 是非均匀的时间。两种时间表示可以如图 所示：

（i） 绝对时间T

（ii）个人时间

上图为两种时间示意图，图（i）为绝对时间，其每个时间间隔都是均匀的，等间隔流失的。图（ii）为个人时间，明显可看出每个时间间隔是非均匀的。

可以这样理解，如果有某个观测者甲对一运动物体进行观测，他使用时钟进行记时，但是这个时钟由于某种故障而越走越慢，每个时间间隔都比上一个时间间隔要长，那么观测者 甲拥有一个自己并不知道的错误的计时时间 。另一个观测者乙同样对该运动的物体

10

湖北大学本科毕业论文（设计） 进行观测，但是他的时间是宇宙时间T ，每个时间间隔都是均匀的。两个时间的转换关系为Tg。

设观测者 甲(拥有个人时间)对运动物体测出的速度为u ，那么他将测出该运动物体运行的距离为：

SAtud

0t那么，对于同一个运动物体，观测者乙(拥有宇宙时间T )测出的该物体运动的距离为：

SBtudg

0t实际上，该运动物体运动的距离应该是观测者B 测量出来的距离，即SBt

另外，考虑分数阶微积分的 Riemann-Liouville 分数阶积分定义：

tt1p1Dftftdfdg （1.6.1） 0tp00p如果将式中的函数 看作物体运动的速度，绝对时间T 与个人时间满足关系 T = g如式（1.6.1）所示中的函数f看作物体运动的速度u，绝对时间T 与个人时间 满足

关系Tg如式(1.6.1)所示。则 Riemann-Liouville 定义下的分数阶积分运算其物理意义就是：在个人时间测度上，所观测的运动物体以速度f 实际走过的距离 ，记为：

SBt0Dtpft

类似地对分数阶微分运算分析其物理意义，同样使用 Riemann-Liouville 分数阶微分定义。根据物体的运动速度是该物体移动距离的一阶导数关系，可得：

dduBtSBtSBt0Dtpft0Dt1pft

dtdt也就是说，观测者乙测量到物体的速度uBt就是函数ft的1p阶Riemann-Liouville 微分。这就是 Riemann-Liouville 定义下的分数阶微分运算的物理解释。

1.7 分数阶微积分在自然中的存在

分数阶微积分现象在自然界中是存在的，比如生物的神经电脉冲信号等。能否找到自

然界中的某种物质或材料，或者通过实际的电子电路实现分数阶微积分，一直是众多科研工作者研究的问题。能否在自然界中找到一种材料或器件，它的某种电气或物理特性与分数阶微积分特性相同，这样就可以利用该物质的这一特性进行分数阶微积分的有效运算，从而避免了较为繁琐的模拟或数字电路实现。另一方面，如果该材料存在，那么也能证明分数阶微积分确实蕴含在自然界中，人们有必要对分数阶微积分的独特特性进行深入研究和认识。如果自然界中确实存在着具有分数阶微积分特性的材料，那么同时也肯定了分数

11

湖北大学本科毕业论文（设计） 阶微积分的存在，其性质不同于通常的整数阶微积分特性，如果此时仍使用整数阶微积分特性进行分析，是达不到预期效果的。

2 分数阶微分的应用

分数阶导数在很多领域都有应用，下面拿与生活联系比较紧密的气候研究、医学图像处理、地震分析为例进行进一步地阐述与说明。

2.1 医学图像处理

医学图像一般是指为了清楚地看到病人内部的局部器官病变情况而通过一定的设备仪器得到的图片，例如CT、B超等图片。由于设备，技术等方面的原因，得到的医学图像有可能模糊不清。图像的不清晰对临床诊断带来很大的麻烦。所以要考虑怎样处理，可以得到更清晰的医学图像。由于分数阶微分的阶数是可以连续变化的，因此在图像处理的过程中可以通过调节微分算子的阶数，在适当增加高频信息的同时又保留一定的低频信息以达到图像增强的目的，这是整数阶微分算子不能实现的。

图像增强的方法分为两大类：空域增强法和频域增强法。空域增强法是对图像的像素直接进行处理。频域增强法主要是将图像在频域里进行图像平滑和锐化滤波等，然后进行反变换以达到图像增强的方法。目前图像增强的算法主要有：直方图均衡化，直方图规定化，直方图映射变换灰度变换，邻域平均法，低通滤波法，多图像平均法，高通滤波法，微分锐化法和同态图像增强法。文中用的图像增强的方法属于空域中的微分增强法。

从前面知道分数阶导数的G-L定义如式(1)。将一元函数的持续期按等间隔等分，由此可以推导出一元函数分数阶微分的差分表达式：

dafta1a1aft2ftft1a+…ftn adt2!n!an1将一元的分数阶微积分的表达式推到二元函数，令二元函数的分数阶微积分的表达式在x方向和y方向的形式如下：

a1afx2,yaffx,yafx1,y ax2a1afx,y1affx,yafx,y1 xa2如此，由此式，可以得到作用到二维图像上的分数阶微积分在x,y方向的模板。为了可以简化运算，使得掩模操作更加方便，可以将x,y方向的模板叠加在一起，并且从模板中可看出，x,y的作用中心点并不一致，可以将作用点选作为模板中心，然后又考虑到倾斜方向上其他因素的影响，可以加上适当的系数，得到了新的微分模板。同此也可以得到分数阶积分的模板，显然a的值便变为负数。

为了使处理之后的每个像素的灰度值不变化，一边保持原本图像的大概包络，将模板中的每个系数都除以4+2a-6a。

用模板对已有的图像进行滤波处理，令m(p,q)为模板中方位为（p，q）的元素，f

212

湖北大学本科毕业论文（设计） （x,y）是模板作用于图像的中心像素，f（x+p,y+q）是离中心f（x,y）距离（p，q）的点，将得到模板的映射规则：Dax,ymp,qfxp,yq

p1q111移动模板，使其可以作用到图片中的每一个元素，到达曾强图像的目的。当元素距离图片边缘小于两个点位的距离时，模板不能完全的覆盖到图片上，此时将其忽略，可以不进行运算。

采用的微分阶数都是正数，，图像的细节和边缘就变得更明显，原图比较模糊的地方得到了明显的改善，但是在当阶数大到一定程度时就会出现明显的锐化现象，图像质反而降低。由此可以看出，文中设计的微分增强掩模在阶数0.2—0.6时，增强效果比较明显。当微分阶数为负数的时候，设计的掩模就可以看成是积分掩模，其对图像具有平滑去噪的效果，当获得的图像具有一定噪声的时候，就可以通过对负阶数的选取，从而得到去噪后的图像。

下面分别运用空间滤波的分数阶微分算子(阶数=0．5)、巴特沃斯算子、拉普拉斯算子对人脑的CT图像进行增强，结果如图3所示。从图中的比较可以看出，经过拉普拉斯算子处理过的图像出现了明显的锐化现象，经过巴特沃斯滤波后图像增强效果并不理想，而通过设计的掩模对图像处理后明显对细节进行了加强，同时，可以根据不同需求，调节分数阶微分阶数以得到符合医学上面需要的图片。分数阶微分掩模有很好的图像增强功能，同时可以通过改变微分的阶数，来适当增强高频信息的同时保留一定的低频信息以满足不同的要求，需要注意的是，当微分阶数过大时将会出现明显的锐化现象。与其它空间滤波算子的效果比较，可以看出，分数阶微分算子弥补了传统图像增强算子不能通过改变参数来得到连续变化的增强效果。当微分阶数为负数时，分数阶微分算子就成了分数阶积分算子，也可以通过调节积分阶数，以得到平滑效果不同的平滑图像，因此，该算子可以将分数阶微分和积分进行统一，在结构简单的同时也具有很强的灵活性。同时作为一个新兴的图像处理方法，分数阶微积分还具有更大的发展空间。

2.2 天气和气候的研究

我们都知道没有一天天气是一样的，而气候的预测也不可能提到日程上来研究。这说明天气和气候的研究是比较困难的。天气和气候虽然遵从流体力学规律，但是却显示出随机性，研究天气和气候之间的关系必须引入分数阶的导数和积分，从物理上讲不外乎说明天气和气候的随机程度是不相同的。为此提出气候的q(0 ≤q≤1） 阶微商是天气。此时引入天气和气候之间的桥梁——分数阶导数，这为天气与气候的研究带来很大的方便。

现在从分数阶微分基本定义出发，可以作用于二维医学图像的分数阶微分掩模，掩模可以根据对图像的需求进行增强。通过实验证明，这个方法可以有效完成对医学图像的处理,并且弥补了传统方法不能连续改变处理效果的缺点，是一种简单可行并且效果较好的图像增强方法。

2.3地震奇异性分析

我们知道传统的地震解释主要是观测地震资料的振幅及相位的变化，而振幅往往并不能反映真实的地质情况。地震界面可能是岩性分界面也可能是岩性过渡带，岩性过渡带的地震反射波是入射波的分数阶导数。

13

湖北大学本科毕业论文（设计） 因此我们将分数阶导数引入地震属性计算中，构建一种对波形敏感而对振幅变化不敏感的新属性——奇异性，用以刻画反射界面的横向变化。

方法的基本原理是首先计算地震子波的不同分数阶导数，然后利用匹配追踪算法将地震数据分解成地震子波的不同分数阶导数，进而获得反射波同相轴的分数阶。对胜利油田某区块实际二维地震资料进行了试处理，结果表明分数阶导数剖面能很好地描述不整合面，反映实际界面的横向变化。

14

湖北大学本科毕业论文（设计） 参考文献

[1] 陈遵德, 陈富贵，《非整数阶微积分的滤波特性及数值算法.数值计算与计算机应用》

[J]，1999 [2] 王竹溪，郭敦仁, 《特殊函数概论》[M]，北京，北京大学出版社, 2000

[3] 林孔容,《关于分数阶导数的几种不同的定义的分析与比较》[J],闽江学院学报，2003 [4] 刘式达，刘式适，《特殊函数(第二版)》[M]，气象出版社, 2002

[5] 段俊生，天津轻工业学院学报，《含 Caputo 分数阶导数的分数阶微分方程》[J]，

2003 [6] 赵元英，袁晓，滕旭东等，《常用周期信号的分数微分运算》[M]，四川大学学报(工

程科学版)，2004 [7] 袁晓,张红雨,虞厥邦,《分数导数与数字微分器设计》[J] 电子学报 , 2004 [8] 陶然,邓兵,王越, 《分数阶傅里叶变换及其应用》[M], 北清华大学出版社, 2009 [9] 王在华，《分数阶微积分》[J], 科学中国人,2011

[10] Vinagre B M and YangQuan Chen, Fractional calculus applications in automatic control and robotics, Las Vegas, 41st IEEE CDC, Tutorial workshop 2, 2002

15

湖北大学本科毕业论文（设计） 致谢

本文是在陈老师精心指导和大力支持下完成的。陈娜老师以其严谨求实的治学态度、高度的敬业精神、兢兢业业、孜孜以求的工作作风和大胆创新的进取精神对我产生重要影响。她渊博的知识、开阔的视野和敏锐的思维给了我深深的启迪。同时，在此次毕业论文设计过程中我也学到了许多的关于分数阶微积分方面的知识，视野得到了极大的开阔。同时我还要感谢我们班的同学，是他们的督促与指导给了我好大的动力。

另外，我还要特别感谢学校为我完成这篇论文提供了巨大的帮助，使我得以顺利完成论文。最后，再次对关心、帮助我的老师和同学表示衷心地感谢

16

湖北大学本科毕业论文（设计）外文翻译

1