zoukankan      html  css  js  c++  java
  • (四)分数阶微积分


    我们重点考察$R-L$型分数阶微积分的性质,简记${}_{0}^{RL}D_{t}^{eta}=D_{t}^{eta}$,若无特殊说明。
    a). 线性性
    $$D_{t}^{eta}[f(t)+g(t)]=D_{t}^{eta}f(t)+D_{t}^{eta}g(t)$$
    $$D_{t}^{eta}lambda f(t)=lambda D_{t}^{eta}f(t) $$
    证明:直接带入定义验算即可.设$m=[eta]+1$
    $${}_{0}^{RL}D_{t}^{eta}f(t)=frac{1}{Gamma(m-eta)}frac{d^{m}}{dt^{m}}int_{0}^{t}(t- au)^{m-eta-1}f( au)d au$$
    b). 积分的叠加性
    $$D_{t}^{-alpha}D_{t}^{-eta}f(t)=D_{t}^{-alpha-eta}f(t) (alpha,eta>0)$$
    证明:对整数阶积分结论是显然的,对于分数阶R-L积分仍然具有叠加性。
    由定义知
    $${}_{0}D_{t}^{-eta}f(t)=frac{1}{Gamma(eta)}int_{0}^{t}(t-x)^{eta-1}f(x):=g(t)$$
    那么
    egin{eqnarray*}
    {}_{0}^{}D_{t}^{-alpha}g(t)&=&frac{1}{Gamma(alpha)}int_{0}^{t}(t- au)^{alpha-1}g( au) d au\
    &=&frac{1}{Gamma(alpha)Gamma(eta)}int_{0}^{t}(t- au)^{alpha-1}d au int_{0}^{ au}( au-x)^{eta-1}f(x)dx\
    &=&frac{1}{Gamma(alpha)Gamma(eta)}int_{0}^{t}f(x)dx int_{x}^{t}(t- au)^{alpha-1}( au-x)^{eta-1}d au(交换积分次序)\
    &=&frac{1}{Gamma(alpha)Gamma(eta)}int_{0}^{t}f(x)dxint_{0}^{1}(t-x)^{alpha+eta-1}(1-xi)^{alpha-1}xi ^{eta-1}dxi (Let xi=frac{ au-x}{t-x})\
    &=&frac{B(alpha,eta)}{Gamma(alpha)Gamma(eta)}int_{0}^{t}(t-x)^{alpha+eta-1}f(x)dx\
    &=&frac{1}{Gamma(alpha+eta)}int_{0}^{t}(t-x)^{alpha+eta-1}f(x)dx\
    &=&{}_{0}D_{t}^{-alpha-eta}f(t)
    end{eqnarray*}
    由此我们也得到了积分满足交换性,即
    $$D_{t}^{-alpha}D_{t}^{-eta}f(t)=D_{t}^{-alpha-eta}f(t)=D_{t}^{-eta}D_{t}^{-alpha}f(t) (alpha,eta>0)$$
    c). 上式考虑了积分叠加的情形,对于连续函数$f(t)$考虑混合运算“先积分再微分”.(还记得R-L定义思路$D^{eta}=D^{m}D^{-(m-eta)}$)
    $${}_{0}^{}D_{t}^{alpha}{}_{0}D_{t}^{-eta}f(t)={}_{0}^{}D_{t}^{alpha-eta}f(t) (alpha>0,eta>0)$$
    证明:先探讨一种特殊的情形
    $$D^{lambda}D^{-lambda}f(t)=f(t) (lambda>0)$$
    当$lambda$为整数时结论显然成立。不妨设$k-1 leq lambda<k$,即$k=[lambda]+1$.故由定义有
    $$D^{lambda}=D^{k}D^{-(k-lambda)}$$
    代入下式
    $$D^{lambda}D^{-lambda}=D^{k}D^{-(k-lambda)}D^{-lambda}=D^{k}D^{-k}=I (use b).)$$
    $ullet$若$alpha<eta$,记$m=[alpha]+1,n=[alpha-eta]+1$,则
    egin{eqnarray*}
    D^{alpha}D^{-eta}f(t)&=&D^{m}D^{-m-alpha}D^{-eta}f(t)
    &=&D^{m}D^{-(m-alpha+eta)}f(t)\
    &=&D^{m}D^{-m}D^{-(eta-alpha)}f(t)\
    &=&D^{alpha-eta}f(t)\
    end{eqnarray*}
    $ullet$若$alpha>eta$,则
    egin{eqnarray*}
    D^{alpha}D^{-eta}f(t)&=&D^{m}D^{-m-alpha}D^{-eta}f(t)\
    &=&D^{m}D^{-(m-alpha+eta)}f(t)\
    &=&D^{n}D^{m-n}D^{-(eta-alpha)}f(t)\
    &=&D^{n}D^{-[n-(alpha-eta)]}\
    &=&D^{alpha-eta}f(t)\
    end{eqnarray*}
    综上有
    $${}_{0}^{}D_{t}^{alpha}{}_{0}D_{t}^{-eta}f(t)={}_{0}^{}D_{t}^{alpha-eta}f(t) (alpha>0,eta>0)$$
    d). 我们前面讲了一般$D^{alpha}D^{-eta}f(t) eq D^{-eta}D^{-alpha}f(t)$,也就是说积分算子与微分算子的交换性并不总是成立的,在经典的微积分里同样要满足一些特定的条件。
    例如$f(x)equiv C$
    先微分再积分:
    $$D^{-1}Df(x)=0$$
    先积分再微分:
    $$DD^{-1}f(x)=C$$
    下面我们研究混合运算中的“先微分再积分”情形,即$D^{-eta}D^{alpha},(eta>0,alpha>0)$形式.
    首先,我们有结论:
    $$D^{-lambda}D^{lambda}f(t)=f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+1)}t^{lambda-j}$$
    证明:回忆当$f(x,y)$满足连续性条件且可微时
    $$frac{d}{dx}int_{a}^{b}f(x,y)dy=int_{a}^{b}frac{partial f(x,y)}{partial x}dy$$
    $$frac{d}{dx}int_{alpha (x)}^{eta(x)}f(x,y)dy=f(x,eta(x))eta'(x)-f(x,alpha(x))alpha'(x)+int_{alpha(x)}^{eta(x)}frac{partial f(x,y)}{partial x}dy$$
    特殊地,
    $$frac{d}{dx}int_{0}^{x}f(x,y)dy=f(x,x)+int_{0}^{x}frac{partial f(x,y)}{partial x}dy$$
    设$g(t, au)=(t- au)^{lambda}D^{lambda}f( au)$, 则$g( au, au)=0$
    $$frac{partial g}{partial t}=lambda (t- au)^{lambda -1}D^{lambda}f( au)$$
    所以
    egin{eqnarray*}
    D^{-lambda}D^{lambda}f(t)&=&frac{1}{Gamma(lambda)}int_{0}^{t}(t- au)^{lambda-1}D^{lambda}f( au)d au\
    &=&frac{1}{Gamma(lambda+1)}frac{d}{dt}int_{0}^{t}(t- au)^{lambda}D^{lambda}f( au)d au\
    end{eqnarray*}
    注意上式中不能轻易使用分部积分法,因为尚未验证分数阶微分是否具有叠加性。
    不妨设$m-1 leq lambda <m$,即$m=[lambda]+1$,由定义和$ c).$知
    egin{eqnarray*}
    frac{1}{Gamma(lambda+1)}int_{0}^{t}(t- au)^{lambda}D^{lambda}f( au)d au&=&
    frac{1}{Gamma(lambda+1)}int_{0}^{t}(t- au)^{lambda}D^{m}D^{-(m-lambda)}f( au)d au\
    &=&frac{1}{Gamma(lambda+1)}int_{0}^{t}(t- au)^{lambda}dD^{m-1}D^{-(m-lambda)}f( au)\
    &=&frac{D^{lambda-1}f( au)}{Gamma(lambda+1)}(t- au)^{lambda}|_{0}^{t}
    +frac{1}{Gamma(lambda)}int_{0}^{t}(t- au)^{lambda-1}D^{lambda-1}f( au)d au\
    &=&frac{1}{Gamma(lambda)}int_{0}^{t}(t- au)^{lambda-1}D^{lambda-1}f( au)d au-frac{D^{lambda-1 }f(0)}{Gamma(lambda+1)}t^{lambda}\
    &=&frac{1}{Gamma(lambda-1)}int_{0}^{t}(t- au)^{lambda-2}D^{lambda-2}f( au)d au-frac{D^{lambda -2}f(0)}{Gamma(lambda)}t^{lambda-1}-frac{D^{lambda-1}f(0)}{Gamma(lambda+1)}t^{lambda}\
    &=&cdots cdots cdots \
    &=&frac{1}{Gamma(lambda-m+1)}int_{0}^{t}(t- au)^{lambda-m}D^{-m-lambda}f( au)d au-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+2)}t^{lambda-j+1}\
    &=&D^{-(lambda-m+1)}D^{-(m-lambda)}f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+2)}t^{lambda-j+1}\
    &=&D_{t}^{-1}f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+2)}t^{lambda-j+1}\
    end{eqnarray*}

    egin{eqnarray*}
    D^{-lambda}D^{lambda}f(t)&=&frac{d}{dt}[D_{t}^{-1}f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+2)}t^{lambda-j+1}]\
    &=&f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(lambda-j+1)}t^{lambda-j}\
    end{eqnarray*}
    下面我们考虑$D^{-eta}D^{alpha},(eta>0,alpha>0)$形式.
    $ullet$若$eta leq alpha ,$则
    egin{eqnarray*}
    D^{-eta}D^{alpha}f(t)&=&D^{alpha-eta}D^{-(alpha-eta)}D^{-eta}D^{alpha}f(t)\
    &=&D^{alpha-eta}D^{-alpha}D^{alpha}f(t) (use b). )\
    &=&D^{alpha-eta}[f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(alpha-j+1)}t^{alpha-j}]\
    &=&D^{alpha-eta}f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(alpha-j+1)}D^{alpha-eta}t^{alpha-j}\
    &=&D^{alpha-eta}f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(eta-j+1)}t^{eta-j}\
    end{eqnarray*}
    $ullet$若$eta>alpha$,则
    egin{eqnarray*}
    D^{-eta}D^{alpha}f(t)&=&D^{-(eta-alpha)}D^{-alpha}D^{alpha}f(t)\
    &=&D^{-(eta-alpha)}[f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(alpha-j+1)}t^{alpha-j}]\
    &=&D^{alpha-eta}f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(alpha-j+1)}D^{-(eta-alpha)}t^{alpha-j}\
    &=&D^{alpha-eta}f(t)-sum_{j=1}^{m}frac{D^{lambda-j}f(0)}{Gamma(eta-j+1)}t^{eta-j}\
    end{eqnarray*}
    总而言之
    $$D^{-eta}D^{alpha}f(t)=D^{alpha-eta}f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(eta-j+1)}t^{eta-j}$$
    e). 前面我们讨论了积分叠加的情形和“先积分再求导”和“先求导再积分”的情形,现在我们考虑微分是否具有叠加性,即
    $$D^{eta}D^{alpha}f(t)=D^{eta+alpha}f(t) (alpha>0,eta>0) ( extbf{?})$$
    一般结论并不是总成立的,例如$f(t)=1$,则有
    $$D^{eta}Df(t)=0$$
    $$DD^{eta}f(t)=Dfrac{t^{-eta}}{Gamma(1-eta)}=-eta frac{t^{-eta-1}}{Gamma(1-eta)} eq 0$$
    设$n-1leq eta <n$即$n=[eta]+1$利用$ d). $ 可得
    $$D^{eta}D^{alpha}=D^{n}D^{-(n-eta)}D^{a}=D^{n}[D^{alpha+eta-n}f(t)-sum_{j=1}^{m}frac{D^{alpha-j}f(0)}{Gamma(n-eta-j+1)}t^{n-eta-j}]$$
    f ). 回忆
    二项式定理
    $$(f(x)+g(x))^{n}=C_{n}^{k}f^{n-k}(x)g^{k}(x)$$
    Leibnitz公式
    $$(f(x)g(x))^{(n)}=C_{n}^{k}f^{(n-k)}(x)g^{(k)}(x)$$
    分数阶R-L积分的Leibniz公式为
    $$D^{- u}[f(x)+g(x)]=sum_{k=0}^{infty}C_{- u}^{k}D^{k}g(x)D^{- u-k}f(x)$$
    证明:根据定义有
    $$D^{- u}f(x)g(x)=frac{1}{Gamma( u)}int_{0}^{x}(x- au)^{ u-1}f( au)g( au)d au$$
    将$g( au)$在$ au=x$处展成$Taloy$公式
    $$g( au)=g(x)+sum_{k=1}^{infty}frac{D^{k}g(x)}{k!}( au-x)^k$$
    代入定义式得
    egin{eqnarray*}
    D^{- u}f(x)g(x)&=&frac{1}{Gamma( u)}g(x)int_{0}^{x}(x- au)^{ u-1}f( au)d au+
    sum_{k=1}^{infty}frac{D^{k}g(x)}{k!}frac{(-1)^{k}}{Gamma( u)}int_{0}^{x}(x- au)^{k+ u-1}f( au)d au \
    &=&g(x)D^{- u}f(x)+sum_{k=1}^{infty}C_{- u}^{k}D^{k}g(x)D^{- u-k}f(x)\
    &=&sum_{k=0}^{infty}C_{- u}^{k}D^{k}g(x)D^{- u-k}f(x)
    end{eqnarray*}
    其中$$C_{- u}^{k}=frac{(-1)^kGamma(k+ u)}{k!Gamma( u)}$$.
    再根据分数阶导数定义可得到分数阶导数的Leibniz公式
    egin{eqnarray*}
    D^{ u}f(x)g(x)&=&D^{m}D^{-(m- u)}f(x)g(x)\
    &=&D^{m}[sum_{k=0}^{infty}C_{ u-m}^{k}D^{k}g(x)D^{ u-m-k}f(x)]
    end{eqnarray*}
    g ). R-L型分数阶导数与Caputo型导数的关系
    一般两者并不相同,但在下面定理的条件下两者等价。
    定理:若$f(t)$具有$m$阶连续导($m=[eta]+1$),且$f^{k}(a)=0 (k=0,1,2,cdots,m-1)$,则有
    $$_{a}^{RL}D_{t}^{eta}f(t)=_{a}^{C}D_{t}^{eta}f(t)$$
    证明:反复利用分部积分可得
    egin{eqnarray*}
    {}_{a}^{C}D_{t}^{eta}f(t)&=&frac{1}{Gamma(m-eta)}int_{a}^{t}(t- au)^{m-eta-1}f^{(m)}( au)d au\
    &=&frac{(m-eta-1)(m-eta-2)cdots (-eta)}{Gamma(m-eta)}int_{a}^{t}(t- au)^{-eta-1}f( au)d au\
    end{eqnarray*}
    另一方面反复利用
    $$frac{d}{dt}int_{a}^{t}f(t, au)d au=f(t,t)+int_{a}^{t}frac{partial f}{partial x}$$
    egin{eqnarray*}
    {}_{a}^{RL}D_{t}^{eta}f(t)&=&frac{1}{Gamma(m-eta)}frac{d^{m}}{dt^{m}}int_{a}^{t}(t- au)^{m-eta-1}f( au)d au\
    &=&frac{1}{Gamma(m-eta)}frac{d^{m-1}}{dt^{m-1}}frac{d}{dt}int_{a}^{t}(t- au)^{m-eta-1}f( au)d au\
    &=&frac{m-eta-1}{Gamma(m-eta)}frac{d^{m-1}}{dt^{m-1}}int_{a}^{t}(t- au)^{m-eta-2}f( au)d au\
    &=&cdots cdots cdots \
    &=&frac{(m-eta-1)(m-eta-2)cdots (-eta)}{Gamma(m-eta)}int_{a}^{t}(t- au)^{-eta-1}f( au)d au\
    end{eqnarray*}
    所以此时
    $${}_{a}^{RL}D_{t}^{eta}f(t)={}_{a}^{C}D_{t}^{eta}f(t)$$

    h ).分数阶导数与整数阶导数的关系
    评注:分数阶导数是整数阶导数的推广,然而这种推广并不是唯一的,历史上出现了多种版本的分数阶导数,例如$R_L$ 型、$Caputo$型、$Resiz$ 型等各个版本计算结果都不尽相同,都有独立的应用,$Mill$ 和$Ross$ 在此基础上提出了序列微积分的概念在形式上统一了分数阶微积分。
    i). 序列分数阶导数实际上是基于观察所得出的定义,对整数阶的导数有
    $$D^{n}f(t)=DDD cdots Df(t)$$
    推广到更一般的形式
    $$D^{a}f(t)=D^{alpha_{1}}D^{alpha_{2}}cdots D^{alpha_{n}}f(t)$$
    其中$alpha=alpha_{1}+alpha_{2}+ cdots + alpha_{n}$.
    从上面的形式可以看到分解并不是唯一的
    当分解形式是
    $$ D^{alpha}=D^{n}D^{-(n-alpha)}$$
    为$Riemann-Liouville$型分数阶导数.
    当分解形式是
    $$ D^{alpha}=D^{-(n-p)}D^{n}$$
    为$Caputo$型分数阶导数.


    j ). Remark: 利用卷积工具可推广到广义函数的分数阶导数,经典分析理论到近代分析跨越.

  • 相关阅读:
    Win8系统 Python安装
    一些安卓开源框架整理
    Android 媒体键监听以及模拟媒体键盘的实现 demo
    android View 自动 GONE 问题
    Android 定时器TimerTask 简单使用
    关于Android studio 相对 eclipse 优点
    Java序列化与反序列化
    android shape的使用 边框
    Android Studio 修改 包名 package name
    Android WebView Long Press长按保存图片到手机
  • 原文地址:https://www.cnblogs.com/zhangwenbiao/p/3705488.html
Copyright © 2011-2022 走看看