1

分数階微積分入門[2] ~分数階微分の定義~

41
0

この記事では
分数階微積分入門[1] ~ 分数階積分の定義 ~
の続きになってます. 前回は分数階積分を定義したので次に微分の方を定義していきます. 実は分数階積分を定義すれば微分の方は簡単に定義できます.

分数階微分のアイデア

次の定理が分数階微分を定義する上でアイデアである.

m,nN, m>n とし, fC([a,b]) とする. このとき,Dn=DmJamn が成立する.

微分積分学の基本定理 から, f=DmnJmn である. この両辺に微分作用素 Dn 作用させ, DnDmn=Dmに注意すれば,
Dnf=DnDmnJamnf=DmJamn を得る.

分数階積分は前回定義したので, 上の定理を使えば上手く微分が定義できることがわかる. 分数階積分は以下のように定義した.

Riemann-Liouvilleの積分

nR+ とする. 作用素JanL1[a,b]において次のように定義する.
Janf(x)=1Γ(n)ax(xt)n1f(t)dt
Riemann-Liouvillen階積分とよぶ. n=0に対しては, Ja0=I とする. ただし, Iは恒等写像とする.

分数階微分

分数階微分を以下のように定義する.

Riemann-Liouvilleの分数階微分

nR+とする. 作用素 DanDanf=DamJamnf と定める. つまり,
Danf(x)=1Γ(mn)dmdxmax(xt)mn1f(t)dt
これをRiemann-Liouvilleの分数階微分とよぶ. ここで m1<n<m, mN. n=0 に対しては Da0=I, として恒等作用とする.

注意
定理1ではn<mを満たす自然数mならどれでもいいが, 今回はn<mを満たす最小の自然数を選択した. 実はどのようにm選択しても結果は変わらないことがわかっている.

具体的な計算例

1<β,n>0 とし, 自然数mm1<n<m 満たすものとする. f(x)=(xa)βとしたとき, Danf(x)を求めよ.

前回の例1から,
Danf(x)=DamJmnf(x)=Γ(β+1)Γ(mn+β+1)Dam[(xa)mn+β].
特に, nβN のときは右辺の(xa)mn+β の次数はm 以下なので, Dm[(xa)mn+β]=0. 一方, nβN なら,
Danf(x)=Γ(β+1)Γ(β+1n)(xa)βn となる.

なのでこれは古典論と一致しているのがよくわかる.

最後に

今回は簡単な分数階微分の計算をしてみました. このように上手く古典論と一致しているのがよくわかるのですが, 全てが上手くいくわけではありません. 例えば定数関数を微分したとき, 古典論では0になるはずですが, 1/2階微分したものは0にはなりません. 余裕がある人は手を動かしてしてみましょう.
さて, 分数階微分まで定義しましたが, ここで記事が終わるのは面白くないのでもう少し続けようと思います. 次回はexの微分などを計算してみようと思います. その計算は定義通りにしようとすると上手くいかない(原始関数がが初等関数で表せない)ので少し工夫が必要です. それは次回に回しましょう.

ちょっと今回は短くなりましたが, 次回は少し重く?なります. 次々回は分数階微分方程式を予定しています. それでは.

投稿日:21
更新日:21
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 分数階微分のアイデア
  2. 分数階微分
  3. 具体的な計算例
  4. 最後に