Skip to content

定义Definition

Riemann–Liouville分数阶积分

Riemann–Liouville fractional integral · Riemann-Liouville fractional integral · 左侧分数阶积分

在固定记忆起点上定义归一化幂核积分,证明L1与连续输入的接口、阶数半群、幂函数迁移及端点迹,区分阶数复合与时间重启。

积分两次可以合成一个权重为t−s的积分;积分三次,权重变成(t−s)²/2。既然整数次数已经藏在核的幂与阶乘中,就可以让次数取正实数。不过这样得到的算子首先是带记忆的积分,不是把普通导数符号的整数任意改成小数。

形式陈述 ​

固定从零开始记忆 ​

给定有限T>0和实数α>0。对 f∈L1(0,T),定义

(1)(Iαf)(t)=1Γ(α)∫0t(t−s)α−1f(s)ds对几乎处处的 t∈(0,T).

这里 Γ 是正实轴上的Gamma函数,积分采用Lebesgue意义,函数按几乎处处相等识别。它仍在L1中,并且

(2)‖Iαf‖1≤TαΓ(1+α)‖f‖1.

不能仅凭f可积就保证式(1)在每个点有限。

若f在[0,T]连续,则式(1)在每个t>0存在,并有连续延拓 Iαf(0)=0,满足

(3)|(Iαf)(t)|≤tαΓ(1+α)‖f‖∞.

对于0<α<1,还成立显式的Hölder估计

(4)|Iαf(t+h)−Iαf(t)|≤2‖f‖∞Γ(1+α)hα(0≤t<t+h≤T).

这没有宣称输出必定绝对连续或有连续普通导数。

阶数可以相加 ​

对α,β>0,

(5)IαIβf=Iα+βf=IβIαf.

对L1输入,等式按几乎处处理解;对连续输入,它处处成立。另行约定 I0 是恒等算子,不能把α=0直接代入含Γ(α)和不可积核的式(1)。

对p>−1和t>0,

(6)Iα(tp)=Γ(p+1)Γ(p+α+1)tp+α.

若记 ϕp(t)=tp/Γ(p+1),它简化成 Iαϕp=ϕp+α。允许p<0是一个L1输入的扩展,不表示它在零点连续。

直觉

非整数次积分仍是因果卷积 ​

令 gα(r)=rα−1/Γ(α)(r>0)。在有限窗中将它截在(0,T),并把f在窗外补零,式(1)正好是连续变量卷积在(0,T)上的限制。所以它是因果卷积的一个具体幂核实例。截断后的核有有限L1质量;不需要错误地声称未截断的幂核在整个正半轴上可积。

α=1时核为1,只累计总量。α=2时核为t−s,过去的输入已经被后续累计过更久。0<α<1时核在最近的过去变大,但局部面积有限;这就是弱奇异,既不允许直接在s=t取一个有限核值,也不会自动令整个积分发散。

两个端点条件怎样证明半群 ​

先对非负绝对值使用Tonelli定理,有

∫0T∫0t(t−s)α−1Γ(α)|f(s)|dsdt=∫0T(T−s)αΓ(1+α)|f(s)|ds.

这证明式(2),也说明定义几乎处处有限。复合时,在0<s<r<t上出现

1Γ(α)Γ(β)∫st(t−r)α−1(r−s)β−1dr.

令r=s+(t−s)v,内层变成

(t−s)α+β−1Γ(α)Γ(β)B(β,α)=(t−s)α+β−1Γ(α+β).

两个参数都为正,才满足Beta积分的双端收敛条件。对t再积分的总绝对质量,由式(2)的同一估计控制,所以Fubini定理允许给有符号f换序,得到式(5)。连续输入时总绝对积分在每个固定t也有限,故可逐点使用。

同样令s=tv,把幂函数代入式(1),得到 tp+αB(p+1,α)/Γ(α),即式(6)。其中p>−1检查s=0,α>0检查s=t,二者不能互相替代。

连续性与Hölder常数不靠微分奇异核 ​

对0<α<1,幂核随正距离递减。把两个积分分成共同的[0,t]和新增的[t,t+h],取绝对值后得到

‖f‖∞Γ(1+α){tα−(t+h)α+2hα}≤2‖f‖∞Γ(1+α)hα.

这一证明甚至只要求f有界可测。α=1时由积分直接给Lipschitz界。α>1时先用式(5)写成 I1Iα−1f;内层属于L1,外层连续且零点取零。于是全部正阶数都保留了式(3)的连续接口,没有对一个在积分端点无界的核强行逐项求导。

例子与边界

半积分两次是一次,但常数不能漏 ​

取f(t)=1+t。利用Γ(1/2)=√π,

I1/2f(t)=2πt1/2+43πt3/2,

再作半积分得到

I1/2I1/2f(t)=t+t2/2=I1f(t).

若把每次半积分的 1/π 漏掉,复合结果多一个π。这个检查同时检测阶数、幂和归一化,远比只在一个点比较小数更有辨识力。

已有Abel径向反演采用右端积分 AQ(t)=∫tLQ(s)/s−tds,且没有Gamma因子。令q(r)=Q(L−r),则

(7)(AQ)(L−r)=πI1/2q(r).

时间反向与√π都不可省;旧页的两次复合π倍普通积分正好与式(5)吻合。新页不改变旧径向坐标及其反演负号。

零点迹与可积输入是两种接口 ​

连续f总给 Iαf(0)=0,所以常数1不在 Iα:C[0,T]→C[0,T] 的像中。但当0<α<1时,可积的奇异输入

f(t)=t−αΓ(1−α)

却满足 Iαf(t)=1 对每个t>0。它在零点没有连续输入接口;在式(1)空区间处机械填0也不会消掉右极限1。讨论初始值前,必须先说使用哪一种函数类。

单点发散也可能发生。给 α=1/4、f(s)=(1−s)−1/2,0<s<1。f属于L1,但在t=1的核乘积为 (1−s)−5/4,积分无限。这不违反式(2),因为L1结论允许零测集上的例外。

半群的参数是阶数,不是经过的时间 ​

式(5)只说从同一个零点开始的两个积分可以复合。若t>a>0,

(8)Iαf(t)=1Γ(α)∫0a(t−s)α−1f(s)ds+1Γ(α)∫at(t−s)α−1f(s)ds.

第一项随未来t继续变化。即使知道 Iαf(a),通常也不能用这个单数替代整段旧输入。以α=1/2、f=1为例,两段各自从零累计再相加会给 2(a+t−a)/π,而原值为 2t/π;a=1、t=4时二者明显不同。

推论与应用

在适当空间上可以确实撤销积分 ​

当0<α<1且f∈L1,式(5)给 I1−αIαf=I1f。后者绝对连续,其几乎处处导数为f。因此

(9)ddtI1−αIαf=f几乎处处.

对于连续f,微积分基本定理使最后一步逐点成立。特别地Iα在L1上单射:若Iαf=0,先再积分到一阶,然后对原函数求导,得f=0。

反方向不自动是恒等。先取一个一般函数、作分数阶导数、再积分,会受到初始迹约束;Caputo初值接口明确保留了这份常数。积分的半群也不能直接转写成Caputo导数的半群。

有限幂和给出完全符号化的接口 ​

若输入为 ∑j=1mcjϕpj,所有p_j>−1,则分数阶积分只需把每个p_j加α,系数c_j不变。表达式生成需要O(m)次指数加法和记录操作,数值求值还须计Gamma与实幂的精度成本。两个指数相同的项可合并;输入不是有限幂和时,这不是一个通用的有限求积算法。

分数阶记忆终点用这套归一化坐标核半积分复合、非零初始常数和时间重启。它把“阶数相加成立”与“过去可丢弃”分成两个要独立回答的问题。

参考资料
关系图谱17 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
分类位置

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系