首页 > 其他分享 >分数阶导数的基本概念与介绍

分数阶导数的基本概念与介绍

时间:2023-05-15 20:47:21浏览次数:43  
标签:分数 tau aD 导数 int 积分 alpha 基本概念

重要参考

以下是我的整理

一、学习的思路

我们熟知的是整数阶的微积分定义,分数阶微积分的定义一般就是由已知的定义推导而来的。所以学习之前需要看看整数阶微积分的定义再推分数阶即可。

在粘弹性研究中一般使用的是 Riemann-Liouville 分数阶导数定义,我就先主要学习这个。如果用上述定义,那么就得先学习分数阶积分定义,再用整数解导数。所以看似学习分数阶积分,实际上学的是分数阶积分的定义。

二、前置知识

1.柯西重复积分公式

重复积分的定义:

\[\begin{align} I^{(0)}(x)&=f(x)\\[1.2ex] I^{(1)}(x)&=\int^xf(t)\ dt\\[1.2ex] I^{(2)}(x)&=\int^x\int^{t_1}f(t_2)\ dt_2\ dt_1\\[1.2ex] \cdots\\[1.2ex] I^{(n)}(x)&=\int^x\int^{t_1}\cdots \int^{t_{n-1}}f(t_{n-1})\ dt_{n-1}\ dt_{n-2} \cdots\ dt_2 \ dt_1\\[1.2ex] I^{(n)}(x)&=\int^xI^{(n-1)}(t)\ dt\\[1.2ex] \end{align} \]

简单说就是 \(I^{(n)}(x)\) 就是对 \(f(x)\) 求 \(n\) 次积分。

推导柯西重复积分:

截屏2023-05-15 19.51.37

2.\(\Gamma\)函数

在推导分数阶积分的公式的时候用到。在这里只需要知道Gamma函数是阶乘概念的推广即可。

\[\Gamma(n+1)=n! \]

因此上面推导的柯西重复积分定义

\[I^{(n)}(x)=\dfrac{1}{\Gamma(n)}\int^x (x-t)^{n-1}f(t)dt \]

三、R-L 分数阶积分

如前所述,学习 Riemann-Liouville 分数阶导数定义之前,先要学习分数阶积分的概念。

记导数算子为 \(D\),即

\[Df(x)\triangleq\dfrac{df}{dx},D^2f(x)\triangleq\dfrac{d^2f}{dx^2},\cdots,D^nf(x)\triangleq\dfrac{d^nf}{dx^n} \]

记 \(_aD_x^{-1}\) 表示变上限积分,即

\[\begin{align} _aD_x^{-1}f(x)&\triangleq\int_a^xf(t)dt \\[2ex] _aD_x^{-2}f(x)&\triangleq_aD_x^{-1}\left(_aD_x^{-1}f(x)\right)\\[2ex] &\cdots\\[2ex] _aD_x^{-n}f(x)&\triangleq_aD_x^{-1}\left(_aD_x^{n-1}f(x)\right),\ x>a \end{align} \]

易知道,对任意正整数 \(n\),都有

\[D^n\left(_aD^{-n}_x f(x)\right)=f(x) \]

下面讨论 \(n\) 重变上限积分 \(_aD_x^{-n}f(x)\) 的表达形式:

\[_aD_x^{-n}f(x)=\int_a^xd\tau_1\int_a^{\tau_1}d\tau_2\cdots\int_a^{\tau_{n-1}}f(\tau_n)d{\tau_n} \]

由重积分的Cauchy公式,可得:

\[_aD_x^{-n}f(x) =\dfrac{1}{(n-1)!}\int_a^x(x-\tau)^{(n-1)}f(\tau)d\tau =\dfrac{1}{\Gamma(n)}\int_a^x(x-\tau)^{(n-1)}f(\tau)d\tau \]

由于Gamma函数对任意正实数都有定义,因此我们可以其推广到任意正实数情形。

3.1分数阶积分定义

设 \(a>0\),\(f(x)\in L_1[a,b]\),则 \(f(x)\) 的 \(\alpha\) 阶微积分为:

\[_aD_x^{-\alpha}f(x)=\dfrac{1}{\Gamma(\alpha)}\int_a^x(x-\tau)^{(\alpha-1)}f(\tau)d\tau \]

虽然名称中是"分数阶",但事实上是对任意正实数都有定义

显然,当 \(\alpha\) 是正整数时, \(_aD_x^{-\alpha}\)就是通常意义下的整数阶积分。

如果 \(f(x)\) 在 \([a,b]\) 上连续,则可以证明:

\[\lim\limits_{a\to 0}{_aD_x^{-\alpha}f(x)}=f(x) \]

通常记为

\[{_aD_x^{0}f(x)}=f(x) \]

上面定义的是左分数阶导数,这是因为 \(f(x)\) 在 \(x\) 点的分数阶导数是通过 \(f(x)\) 在区间 \([a, x]\) 中的值来表示的。实际上还可以定义右分数阶导数。一个代表 “过去”, 而另一个代表 “未来”

3.2分数阶积分的存在性

参看 华东师范大学 潘建瑜 教授 第一讲 分数阶微分方程

3.3分数阶积分的性质

参看 华东师范大学 潘建瑜 教授 第一讲 分数阶微分方程

四、Riemann-Liouville分数阶导数

Riemann-Liouville 分数阶导数是历史上最早的分数阶导数定义, 也是目前理论研究相对较完善的分数阶导数.

设 \(\alpha>0\) 是任意正实数, \(n\) 是大于 \(\alpha\) 的最小正整数,即 \(0\le n-1\le\alpha<n\) ,则 R-L 分数阶导数定义为:

\[_aD_x^{\alpha}f(x)=D^n(_aD_x^{\alpha-n}f(x))=\dfrac{1}{\Gamma(n-a)}\dfrac{d^n\left(\int_a^x\dfrac{f(\tau)}{(x-\tau)^{(\alpha-n+1)}}d\tau\right)}{dx^n} \]

即先做 \(n-\alpha\) 次分数阶积分,然后再求 \(n\) 次导数。我们注意到 \(0<n-\alpha\le1\)。

截屏2023-05-15 20.24.03

引理

引理1

截屏2023-05-15 20.25.31

引理2

截屏2023-05-15 20.26.46

引理3

截屏2023-05-15 20.29.11

标签:分数,tau,aD,导数,int,积分,alpha,基本概念
From: https://www.cnblogs.com/PeterRabbi/p/17403031.html

相关文章

  • 先利用PCA做主成分分析,通过累计贡献率确定最佳主成分数,然后再进行BP回归预测分析,两个
    先利用PCA做主成分分析,通过累计贡献率确定最佳主成分数,然后再进行BP回归预测分析,两个算法已经都写在一起了,可以直接运行,不用分段运行。ID:2635668336340841......
  • 导数最终章
    思考以前总是分不清,导数和导函数,首先明白什么是导数,导数描述某点的领域的变换\[\lim_{x->x_0}\frac{\Deltay}{\Deltax}=\lim\frac{f(x_0+\Deltax)-f(x_0)}{\Deltax}\]导函数,其实也相当于函数,只是由无数个导数组成的函数题型题型1对应求某点的值,一般就是......
  • IO多路复用基本概念:select,epoll(转)
    原文:https://blog.csdn.net/weixin_45743893/article/details/122970342IO多路复用概念IO多路复用是指内核一旦发现进程指定的一个或者多个IO条件准备读取,它就通知该进程。其实就是在单个线程中通过记录跟踪每一个I/O流的状态来管理多个I/O流。它们三的作用提供一种IO复用的方式......
  • 类和对象的基本概念
    封装:属性和行为作为一个整体,来表现各种事物。将属性和行为加以权限控制(private,public,protected)一些术语:属性(成员属性,成员变量),行为(成员函数,成员方法);统称为成员。实例化(通过一个类,创建一个对象的过程)访问权限:当一个形参是函数对象时,该函数体内部不可以访问该对象的私有属性......
  • Java | 一分钟掌握定时任务 | 1 - 基本概念
    作者:Mars酱声明:本文章由Mars酱原创,部分内容来源于网络,如有疑问请联系本人。转载:欢迎转载,转载前先请联系我!什么是定时任务定时任务就是在指定时间执行的一个或一串动作的行为。现实世界中,比如我们向自己心爱的女友送上祝福;给自己的好友送上祝福。那么,我们为了准时送祝福,会怎么做?当然......
  • 分数比较
    分数比较:比较两个分数的大小。解题思路:输入两个分数,构造可以求出最小公倍数的函数,利用函数求出分母相同时两分子的大小,比较大小,根据要求输出结果。代码:#include<iostream>usingnamespacestd;intbei(intx,inty){   intt,temp,k=x*y;   if(x<y)   {......
  • [NOIP2009 普及组] 分数线划定
    [NOIP2009普及组]分数线划定题目描述世博会志愿者的选拔工作正在A市如火如荼的进行。为了选拔最合适的人才,A市对所有报名的选手进行了笔试,笔试分数达到面试分数线的选手方可进入面试。面试分数线根据计划录取人数的\(150\%\)划定,即如果计划录取\(m\)名志愿者,则面试分数......
  • [SWPUCTF 2022 新生赛]base64 已解决 题目分数:213
    查壳:64位,操作系统是ubantu的,可能会有所不同稍加留意一下,进IDA:依旧是比较题,我们先看看s2里的内容:‘TlNTQ1RGe2Jhc2VfNjRfTlRXUTRaR0ROQzdOfQ==’目标是v3,看看v3调用的函数sub_124C:base64?看看是不是标准码:巨标准,那么直接base64解码就好了:得到NSSCTF{base_64_NTWQ4ZGDNC7N}......
  • Phred质量分数
    Phred质量分数是DNA测序数据中,用来评估碱基质量的一种标准化表示方法。它最初是由美国华盛顿大学的高通量测序专家PhredJ.在1997年发表的论文中提出的。Phred质量分数使用logarithmicscale(即对数刻度)来表示碱基质量值。一般取值范围为0-40之间,其中0表示最差的质量,40表示......
  • 列出真分数序列
    列出真分数序列:按递增顺序依次列出所有分母为40.分子小于40的最简分数。解题思路:利用for循环遍历所有分子的可能,在for循环中利用辗转相除法求出最大公约数,若最大公约数为1则输出,否则继续循环。代码:#include<iostream>#include<iomanip>usingnamespacestd;intmain(){......