分数阶微积分的应用研究

12

无线互联科技

·通信观察

分数阶微积分的应用研究

王福兴，蒲亦非，周激流

摘 要：随着信息技术的飞速发展，分数阶微积分作为数学的重要分支在信号分析与处理等领域得到了广泛的研究和运用。在很多问题的处理过程中，分数阶微积分所拥有的优点正逐渐显露出来。本文阐述了分数阶微积分的时域和频域定义，介绍了其在工程技术中的应用。关键词：分数阶微积分

1 引言

分数阶微积分与菲涅耳衍射、Wigner 分布、小波变换等有着密切的联系，已经被广泛应用于科学研究和工程技术的很多方面，如神经网络、解微分方程、量子力学、衍射理论、光学系统、光图像处理以及雷达、通信、声纳等领域，成为跨学科的研究方向。继1980年被提出之后很快在光学领域得到广泛应用，其在信号处理领域的潜力在20世纪90年代中期得到发掘。近十年来，随着人们对分数阶微积分的研究不断深入，分数阶微积分的理论体系被逐步建立起来，并且被信息领域的研究人员所广泛接受。

随着信息技术的飞速发展，分数阶微积分的应用越来越多地得到人们的重视，在很多问题的处理过程中所拥有的优点正逐渐显露出来。许多科学家从不同的角度进行了不同的尝试，得到了不同的定义，其中比较有名的有G-L定义、R-L定义和Caputo定义[1,2]。

分数阶微积分是整数阶微积分运算的推广，在信号分析与处理等领域得到了广泛应用，特别是在信号的奇异性检测和提取方面具有特殊的作用。分数阶微积分的实现算法分为解析算法、数值算法两种。解析算法如拉普拉斯变换法、傅里叶变换法以及FFT法；数值算法如Zhang and Shimizu法，L-1法和池田法等。另外，在信号处理方面，有学者尝试用数字滤波器和模拟分抗电路来实现分数阶微积分并取得了一些进展。本文阐述了分数阶微积分的时域和频域两种定义方法，介绍了其在工程技术领域的应用。

2 概念

2.1 时域定义

文献[1，2]对分数阶微积分进行了3种时域的定义。

Grumwald-Letnikov定义是从研究连续函数整数阶导数的经典定义出发，将微积分的阶数与因次由整数扩大到分数推衍而来的。Grumwald-Letnikov的υ阶导数定义为：

式中，υ∈R，用 表示υ的整数部分；信号s(t)∈[a,t](α

P∈Z，Z为整数域；当υ>0时，P至少取到[υ]。在定义中，上标G 表示Grumwald-Letnikov定义，上标υ表示求υ阶微分，下标α和t

表示积分式的下界和上界。α为时间t 的初值。当υ为负实数时，公式转化为积分公式。

Riemann-Liouville定义从分数阶微积分应满足的性质入手，对Grumwald-Letnikov定义进行了改进，使之计算简单化。Riemann-Liouville定义的分数阶微分公式如下：

Riemann-Liouville定义的分数阶微分是先进行(n-υ)阶积分，然后进行υ阶微分。Riemann-Liouville定义的分数阶积分为：

(四川大学电子信息学院，四川 成都 610065)

Research on Application of Fractional Calculus

WANG Fu-xing,PU Yi-fei,ZHOU Ji-liu

(School of Electronics and Info.,Sichuan Univ.,Chengdu 610065,China)

Abstract: With the rapid development of information technology, fractional calculus as an important branch

of mathematics in signal analysis and processing and other fields has been widely studied and applied. In the

course of handling the many problems, fractional calculus have advantages are gradually revealed. This paper expounds the fractional calculus in time domain and frequency domain definition, introduces its application in engineering.

Keywords: fractional calculus

(1)

(2)

(3)

13

通信观察·

Caputo定义是对Grumwald-Letnikov定义的另一种改进，其目的是为了让拉普拉斯变换更加简洁，从而便于分数阶微分方程，易于工程应用。Caputo定义是对GL定义的另一种改进，其分数阶微分定义为：

类似地，其分数阶积分定义为：

该定义先进行n阶微分，再进行n-υ阶积分。

在上述众多的定义之间存在着一定的区别和联系，或者在某些条件下它们之间可以

相互转换。可以证明，对于很广一类实际函数来说，GL定义和RL定义是完全等效的。RL定义是GL定义的推广，其应用范围更为广泛。Caputo定义和RL定义都是在Grünwald-Letniko定义的基础上进行改进的，二者主要区别表现在对常数求导的定义上，前者对常数的求导是有界的(为0)，而后者求导是无界的，Caputo定义让分数阶微积分的拉普拉氏变换式更为简洁，它更适合用于分数阶微分方程初值问题的描述[3]。

2.2 傅里叶变换域定义

根据傅里叶变换的性质，υ阶微分算子D υ=D υ

是υ阶微分

乘子函数 的乘性算子，在复数域中分数阶微积分指数形式为：

其中，sgn(*)为符号函数。从通信调制角度看，信号分数阶微分的物理意义可以理解为广义的调幅调相，其振幅随频率成分数阶幂指数变化，其相位是频率的广义Hilbert变换。从信号处理角度看，υ阶微积分运算其实是对信号的一个线性时不变滤波系统，其滤波函数为

3 应用

分数阶微积分运算应用在工程实践中，是近几十年来新兴的研究方向。随着计算技术的发展和对分数阶微积分运算应用研究探索的深入，分数阶微积分运算在多个领域中起到了越来越重要的作用，其中包括机械力学、电子学、化学、生物学、经济学、控制理论、机器人、材料科学、岩石力学、分形理论、电磁场理论、图像与信息处理等方面。特别是在信息科学领域中，一些新颖的应用被相继地实现，如系统建模、曲线拟合、信号滤波、模式识别、图像边界提取、系统辨识、系统稳定性分析等等。

3.1 指纹识别系统[4]

指纹识别的方法随着计算机技术、图像处理与模式识别

方法的不断发展正日臻完善，目前使用的指纹识别系统大多是数字处理指纹识别系统，具有可自动化、可编程等优点。分数阶微积分具有时频旋转特性[5]，可以展示出信号从时域逐步变化到频域的更多特征，这一优势有助于将分数阶微积分的应用研究扩展到了生物医学领域，例如可用于有效解决指纹识别领域中存在的问题和改进传统指纹识别算法的不足，特别是分数阶傅立叶变换算法对发生形状畸变的图像有很好的识别能力，可获得比基于经典傅立叶变换的算法性能更好的识别准确性和灵活性。

分数阶傅里叶变换具有典型的光学含义，它有两种基本的方式可以在光学上实现，一种是利用透镜等分立光学元件，另一种是利用渐变折射率介质(GRIN)。基于分数阶傅里叶变换这一新兴信号处理手段进行指纹等生物特征认证和识别具有研究价值和应用前景。在文献[4]中，作者给出了基于分数阶傅里叶联合变换相关器的指纹识别光电混合系统，如图1所示。CCD1 将采集的待检测的目标指纹图像传输到计算机，同时计算机从指纹数据库中调用作为参考的指纹图像，目标指纹图像和参考指纹图像在计算机中分别进行相位掩模调制，然后计算机将调制后的目标指纹图像输出到空间光调制器SLM1上，将调制后的参考指纹图像输出到空间光调制器SLM2上；再分别经过透镜L1

和L2进行分数阶傅里叶变换，它们的分数阶傅里叶变换谱在分光镜处叠加，由CCD2得到分数傅里叶联合变换功率谱，并经计算机输出到SLM2，经过透镜L2完成分数阶傅里叶变换，最后CCD2 接收到的光强分布即是目标指纹图像和参考指纹图像的分数傅里叶联合变换相关结果，输出到计算机可进行进一步的数据处理和检测统计。

3.2 数字水印系统[6]

随着多媒体技术和网络技术的迅速发展与广泛应用，对多

媒体数字产品的版权保护已成为迫切需要解决的问题。传统的加密技术已经不足以解决问题，而数字水印技术在这方面显示出了巨大的潜力。数字水印算法将一个版权识别代码序列(水印信号)嵌入到图像(空域或变换域)中，利用它可以跟踪数字产品拷贝的非法销售和使用。数字水印系统的实现分为两步：水印

()()()()()01

01

n t v t

v n f D f t d n v t τττ-+=

Γ--⎰ 1,n v n n N -<<∈ (4)

()()()()0

101t v t v

f D f t d v t τττ-=

Γ-⎰ 10,0v t -<<>

(5)

()()ˆv

d i ωω= ()()()()()()()()()()ˆˆˆˆexp ˆˆ,sgn 2v v v v v

v v d iw a i a p v a ωωθωωωπωωθωω⎧==∙=∙⎪

⎨==⎪

⎩ (6)()()()()ˆexp v

v v v d iw i ωωθω==∙ (7)

图1 基于分数阶傅里叶联合变换相关器的指纹识别光电混合系统示意图