GL分数阶微积分

最新推荐文章于 2024-09-01 08:55:24 发布

小斜

最新推荐文章于 2024-09-01 08:55:24 发布

阅读量2.1k

点赞数 2

文章标签：线性代数概率论矩阵

本文链接：https://blog.csdn.net/xiexinning/article/details/121319442

版权

本文介绍了分数阶微积分的基本概念，包括预备公式如Γ函数、二项式系数等，并展示了如何将积分和导数通过变限积分求导统一起来。对于-p<0的情况表示积分，p>0表示导数，同时讨论了整数阶和分数阶的混合运算，给出了相关公式推导。这些理论和公式在分数阶微积分中具有重要意义。

摘要由CSDN通过智能技术生成

预备公式

$\Gamma(z+1)=z\Gamma(z)\tag{1}$

$\Gamma(z)=\int_{0}^{\infty} e^{-t}t^{z-1}\mathrm{d}t=\lim_{n\to\infty}\frac{n!n^{z}}{z(z+1)...(z+n)}=\lim_{n\to\infty}\frac{n^{z}}{\begin{bmatrix} p\\r \end{bmatrix}(p+r)}\tag{2}$

$B(z,w)=\int_0^1\tau^{z-1}(1-\tau)^{w-1}\mathrm{d}\tau=\frac{\Gamma(z)\Gamma(w)}{\Gamma(z+w)}\tag{3}$

$\begin{bmatrix} p\\r \end{bmatrix}=\frac{p(p+1)...(p+r-1)}{r!}=(-1)^r\begin{pmatrix}-p\\r \end{pmatrix}\tag{4}$

$\begin{pmatrix}p\\r \end{pmatrix}=\begin{pmatrix}p-1\\r \end{pmatrix}+\begin{pmatrix}p-1\\r-1 \end{pmatrix}\tag{5}$

$f^{n}(t)=\frac{d^{n}f}{dt^{n}}=\lim_{h\to 0} \frac{1}{h^n}\sum_{r=0}^{n}(-1)^r\begin{pmatrix}n\\r \end{pmatrix} f(t-rh)=\lim_{h\to 0} \frac{1}{h^n}\Delta^nf(t)\tag{6}$

变限积分求导
![avatar]（https://pic2.zhimg.com/aff0982565a44433f231be2a67ed263d_r.jpg?source=1940ef5c）

将积分和导数统一

我们把n阶导数的定义，即公式（6）极限号去掉，并且把p的范围从正整数改为实数，从而定义了下面这个函数
$f^{(p)}_h(t)=\frac{1}{h^p}\sum_{r=0}^{p}(-1)^r\begin{pmatrix}p\\r \end{pmatrix}f(t-rh) \tag{7}$
再定义极限
$\lim_{h\to 0 \atop nh=t-a }f^{(p)}_h(t)=_aD^{p}_tf(t)\tag{8}$

-p<0表示积分

此时是一个变上限积分，区间从a到t，证明利用定积分定义即可
当p是整数时，(8)的极限和p阶乘有关，即
$_aD^{-p}_tf(t)=\frac{1}{(p+1)!}\int_a^t(t-\tau)^{p-1}f(\tau)\mathrm{d}t\tag{9}$
p是任意阶时用gamma函数代替阶乘，证明过程在书2.2。

$_aD^{-p}_tf(t)=\frac{1}{\Gamma(p+1)}\int_a^t(t-\tau)^{p-1}f(\tau)\mathrm{d}t \tag{10}$
如果函数m+1阶连续，由分部积分可知
$_aD^{-p}_tf(t)=\sum_{k=0}^m\frac{f^{(k)}(a)(t-a)^{p+k}}{\Gamma(p+k+1)}+\frac{1}{\Gamma(p+k+1)}\int_a^t(t-\tau)^{p+m}f^{m+1}(\tau)\mathrm{d}\tau \tag{11}$

p>0表示导数

当p>0时上述定义函数具有性质（利用公式（5）证明）
$f^{(p)}_h(t)=\sum_{k=0}^{m}(-1)^{n-k}\begin{pmatrix}p-k-1\\rn-k\end{pmatrix}h^{-p}\Delta^kf(a+kh)+p^{-1}\sum_{r=0}^{n-m-1}(-1)^r\begin{pmatrix}p-m-1\\r \end{pmatrix}\Delta^{m+1}f(t-rh)$
对上式求积分估计，和积分形式是一样的。但是需要满足m<p<m+1
$_aD^{p}_tf(t)=\sum_{k=0}^m\frac{f^{(k)}(a)(t-a)^{-p+k}}{\Gamma(-p+k+1)}+\\\frac{1}{\Gamma(-p+m+1)}\int_a^t(t-\tau)^{-p+m}f^{m+1}(\tau)\mathrm{d}\tau \tag{12}$
特殊的函数 $(t-a)^{\nu}$ 有分数阶导数
$_aD^{p}_t(t-a)^{\nu}=\frac{\Gamma(\nu+1)}{\Gamma(-p+\nu+1)}(t-a)^{\nu-p}\tag{13}$

整数阶和分数阶混合运算

利用变上限积分求导
$\frac{d^n}{dt^n}(_aD^{p}_tf(t))=\sum_{k=0}^s\frac{f^{(k)}(a)(t-a)^{-p+n+k}}{\Gamma(-p+n+k+1)}+\frac{1}{\Gamma(-p+n+s+1)}\int_a^t(t-\tau)^{-p+s+n}f^{s+1}(\tau)\mathrm{d}\tau \\ =_aD^{p+n}_tf(t)\tag{14}$

$_aD^{p}_t(\frac{d^nf(t)}{dt^n})=\sum_{k=0}^s\frac{f^{(n+k)}(a)(t-a)^{-p+k}}{\Gamma(-p+k+1)}+\frac{1}{\Gamma(-p+s+1)}\int_a^t(t-\tau)^{-p+s}f^{n+s+1}(\tau)\mathrm{d}\tau\\ \tag{15}$
证明技巧：s取不同的值 P57 一式s=m+n-1，二式s=m-1，两式相减就有下式
$\frac{d^n}{dt^n}(_aD^{p}_tf(t))= _aD^{p}_t(\frac{d^nf(t)}{dt^n})+\sum_{k=0}^{n-1}\frac{f^{(k)}(a)(t-a)^{-p+n+k}}{\Gamma(-p+n+k+1)}\tag{16}$
显然当t=a，a导数为0有
$\frac{d^n}{dt^n}(_aD^{p}_tf(t))= _aD^{p}_t(\frac{d^nf(t)}{dt^n})=_aD^{p+n}_tf(t)$