分数阶导数

  • 格式:doc
  • 大小:384.50 KB
  • 文档页数:6

分数阶导数1引言我们都熟悉的导数的定义。

通常记作1()()df x D f x dx 或 222()()d f x D f x dx 或这些都是很容易理解的。

我们同样也熟悉一些有关导数的性质,例如[()()]()()D f x f y Df x Df y +=+但是像这样的记号1/21/21/2()D ()d f x f x dx 或者又代表什么意思呢大多数的读者之前肯定没有遇到过导数的阶数是1/2的。

因为几乎没有任何教科书会提到它。

然而,这个概念早在18世纪,Leibnitz 已经开始探讨。

在之后的岁月里,包括L’Hospital, Euler,Lagrange, Laplace, Riemann, Fourier, Liouville 等数学大家和其他一些数学家也出现过或者研究过的概念。

现在,关于“分数微积分”的文献已经大量存在。

近期关于“分数微积分”的两本研究生教材也出版了,就是参考文献[9]和[11]。

此外,两篇在会议上发表的论文[7]和[14]也被收录。

Wheeler 在文献[15]已编制了一些可读性较强,较易理解的资料,虽然这些都还没有正式出版。

本论文的目的是想用一种亲和的口吻去介绍分数阶微积分。

而不是像平常教科书里面的从定义-引理-定理的方法介绍它。

我们寻找了一个新的想法去介绍分数阶导数。

首先我们从熟悉的n 阶导数的例子开始,比如D n axn axe a e =。

然后用其他数字取代自然数字n 。

这种方式,感觉像是侦探一样,步步深入。

我们将寻求蕴含在这个构思里面的数学结构。

我们在探讨了各种思路,对分数阶导数的概念后,才对分数阶导数给出正式定义。

(如果想快速浏览它的正式定义,请参见米勒的优秀论文,参考文献[8]。

)随着探究的深入,我们会不时地让读者去思考一些问题。

对这些问题的答案将在本文的最后一节呈现。

那到底什么是一个分数阶导数呢让我们一起来看看吧……2指数函数的分数阶导数我们将首先研究指数函数ax e 的导数。

因为他们导数的形式,比较容易推广。

我们熟悉axe 的导数的表达式。

12233,,ax ax ax ax ax axD e ae D e a e D e a e ===,在一般情况下,当n 为整数时,n ax n ax D e a e =。

那么我们能不能用1/2取代n ,并记作1/21/2ax ax D e a e =呢我们何不尝试一下为什么不更进一步,让n 是一个无理数或者复数比如1+i 我们大胆地写作ax ax D e a e αα=, (1)对任意一个α,无论是整数,有理数,无理数,还是复数。

当α是负整数时,考虑(1)式的意义是很有趣的。

我们自然希望有1(())ax ax e D D e -=成立。

因为1(())axax e D e a=,所以我们有1()ax ax D e e dx -=⎰。

同理2()ax ax D e e dxdx -=⎰⎰。

当α是负整数时,我们将D α看作是n 次迭代的积分是合理。

当α是正实数,D α代表导数,当α是负实数,D α代表积分。

请注意,我们还没对一般函数给出分数阶导数的定义。

但是,如果这一定义被发现,我们期望指数函数的分数阶导数遵循关系式(1)。

我们注意到,刘维尔在他的论文[5]和[6]中就是采用这种方法去考虑微分的。

问题问题1:在上述情况下,12121212()a x a x a x a xD c e c e c De c De α+=+成立吗问题2:在上述情况下,ax ax D D e D e αβαβ+=成立吗问题3:上述1()ax ax D e e dx-=⎰和2()ax ax D e e dxdx-=⎰⎰,真的正确吗还是遗漏了一些东西问题4:用蕴含在(1)式的想法,怎样对一般性的函数求分数阶导数3三角函数:正弦函数和余弦函数我们对于正弦函数的导数很熟悉:012sin sin ,sin cos ,sin sin ,D x x D x x D x x ===-这些对于寻求1/2sin Dx ,并没有明显的规律。

但是,当我们画出这些函数的图形时,会挖掘出其中的规律。

即每当我们求一次微分,sin x 的图像向左平移/2π。

所以对sin x 求n 次微分,那么得到的图像就是sin x 向左平移/2n π,即得到sin sin()2n n D x x π=+。

如前,我们用任意数α替换正整数n 。

所以,我们得到正弦函数的任意α次导数的表达式,同理我们也得到余弦函数的:sin sin(),cos cos().22D x x D x x αααπαπ=+=+(2)在得到表达式(2)之后,我们自然想,这个猜测与指数函数的结果是否保持一致。

为了验证这个猜测,我们可以使用欧拉公式cos sin ixe x i x =+。

利用表达式(1),我们可以计算得到(/2)cos()sin()22ix ix i ix D e i e e e x i x ααπααπαπ===+++,这与(2)式是吻合的。

问题问题5:sin()D ax α是什么4px 的导数我们现在看看x 次方的导数。

我们以px 为例有:012,,(1),,(1)(2)(1).(3)p p p p p p n p p n D x x D x px D x p p x D x p p p p n x -===-=---+表达式(3)用连乘()!p n -的分子和分母去替换,则得到结果如下(1)(2)(1)()(1)1!(4)()(1)1()!p p n p n p p p p n p n p n p x x x p n p n p n -----+---==----上式就是n pD x 的一般表达式。

我们通过伽玛函数,用任意数α替换正整数n 。

当(4)式中的p 和n 是不是自然数时,伽玛函数使他们在替换后任然有意义。

伽马函数是欧拉在18世纪引进的概念。

当时是推广记号!z ,当z 不是整数时。

它的定义是1()d t z z e tt ∞--Γ=⎰,它具有这样的性质(+1)!z z Γ=。

那么我们可以将表达式(4)重新写作(1),(1)npp n p D x x p n -Γ+=Γ-+这使得当n 不是整数式,(4)式还是有意义的。

所以对于任意的α,我们写作(1)(5)(1)p p p D x x p ααα-Γ+=Γ-+利用(5)式,我们可以将分数阶导数延伸到很多的函数。

因为对于任意给定的函数,我们可以利用Taylor 级数展开成多项式的形式,0(),nn n f x a x∞==∑假设我们可以对()f x 进行任意次微分,那么我们得到(1)().(6)(1)nn n nn n n D f x a D x a x n αααα∞∞-==Γ+==Γ-+∑∑最终那个表达式(6)呈现出具有作为分数阶导数定义候选项的气质。

因为大量的函数都可以利用Taylor 公式展开成幂级数的形式。

然后,我们很快会发现它会导致矛盾的产生。

问题问题6:()D f x α是否有几何意义5一个神秘的矛盾我们将xe 的分数阶导数写为(7)x xD e e α=现在让我们拿它与(6)式进行对比,看看他们是否一致。

从Taylor 级数来看,01,!xnn e x n ∞==∑结合(6)式,我们得到如下表达式 0.(8)(1)nxn x D e n αα∞==Γ-+∑但是,(7)及(8)是不等价的,除非α是整数。

当α是整数时,(8)式的右侧是xe 的级数形式,只是用不同的表达方式。

但是当α不是整数时,我们得到两个完全不一样的函数。

我们发现了历史上引起大问题的矛盾。

这看起来好像我们,指数函数的分数阶导数的表达式(1)与次方函数的分数阶导数的公式(6)是相互矛盾。

正是因为有这样一个矛盾,所以分数阶微积分一般不会出现在初等阶段的教科书里面。

在传统的微积分中,导数的次数是整数次的,求导的函数是初等函数。

不幸的是,在分数阶微积分中,这是不正确的。

通常,一个初等函数的分数阶导数是较高级的超越函数。

关于分数阶导数的表格,请参阅文献[3]。

此时,您可能会问我们怎么继续探究呢这个谜团将在之后的部分中被解决。

敬请关注……6多重迭代积分我们一直在谈论导数。

积分也是反复被提及的。

我们可以写1()()D f x f x dx -=⎰,但是等式右边是不确定的。

我们可以写作1()()xD f x f t dt -=⎰。

第二次积分可以写成221120()()xt D f x f t dt dt -=⎰⎰。

积分区域是图1中的三角形。

如果我们交换积分的顺序,那么图1的右侧图可以表现出121210()()x xt D f x f t dt dt -=⎰⎰。

因为1()f t 不是一个关于2t 的函数,所以可以将里面的积分移到外面,即121211110()()()()x x xt D f x f t dt dt f t x t dt -==-⎰⎰⎰或者2()()()xDf x f t x t dt -=-⎰。

使用相同的过程步骤,我们可以写出32430011()()(),()()(),223x x D f x f t x t dt D f x f t x t dt --=-=-⋅⎰⎰ 在一般情况下,11()()().(1)!x n n D f x f t x t dt n --=--⎰现在,我们用先前做的方法,用任意数α替换n -,用伽玛函数替换阶乘,然后得到101()().(9)()()x f t dtD f x x t ααα+=Γ--⎰ 这个一般性的表达式(使用积分)的分数阶导数表达式,有成为定义的潜力。

但是存在一个问题。

如果1,α>-该积分是反常积分。

因为当,0.t x x t →-→对任意0α≥,积分是发散的。

当10,α-<<反常积分收敛。

所以当α是负数时,原表达式是正确的。

因此当α是负数时(9)式收敛,即它是一个分数阶次积分。

在我们结束这一部分之前,需要提下,趋于零的下极限是任意的。

可以简单的认为存在下极限b 。

但是会造成最后结果表达式的不同。

正因为如此,很多这个领域的研究人员使用符号()bx D f x α。

这个符号说明了极限过程是从b 到x 的。

这样我们从(9)式得到11()().(10)()()x b xb f t dtD f x x t ααα+=Γ--⎰问题问题7:如下分数阶微分b 的下极限是什么(1)()()(1)p p b x p D x c x c p ααα-Γ+-=-Γ-+7解秘现在,你可以开始去发现前面哪些地方出错了。

我们对于分数阶积分包含极限,并不感到惊讶。

因为积分是涉及到极限的。

然而普通的导数不涉及积分的极限,没有人希望分数阶导数包含这样的极限。

我们认为,导数是函数的局部性质。