[論文レビュー] On fractional Bessel equation and the description of corneal topography

W. Okrasiński, Łukasz Płociniczak|arXiv (Cornell University)|Jan 12, 2012

Fractional Differential Equations Solutions参考文献 2被引用数 3

ひとこと要約

本稿では、リーマン＝リウヴィルの分数階微分を用いて古典的な修正ベッセル関数を一般化することで、人間の角膜トポグラフィーを記述する分数階ベッセル微分方程式モデルを提案する。解は収束するべき級数として表現され、実際の角膜データに0.014–0.03 mmの誤差で適合し、分数階数α ∈ (0,1)に対して物理的整合性と漸近的挙動が向上していることが示された。

ABSTRACT

In this note we apply a modified fractional Bessel differential equation to the problem of describing corneal topography. We find the solution in terms of the power series. This solution has an interesting behavior at infinity which is a generalization of the classical results for modified Bessel function of order 0. Our model fits the real corneal geometry data with an error of order of a few per cent.

研究の動機と目的

• 古典的な円錐型またはゼルニケ基底モデルを一般化する、物理的根拠を持つ角膜トポグラフィーの数学的モデルの構築を目的とする。
• 分数階微分を角膜幾何学のモデリングに組み込み、物理的整合性と柔軟性を向上させることを目的とする。
• 0階数の修正分数階ベッセル微分方程式の解を導出し、x → ∞におけるその漸近的挙動を分析することを目的とする。
• 実際の角膜トポグラフィー測定データに対してモデルを検証し、外側および内側の角膜表面を最小限の誤差で適合させることを目的とする。

提案手法

• リーマン＝リウヴィルの分数階微分を用いた、修正分数階ベッセル微分方程式に基づく境界値問題を定式化する。
• 変数変換を用いて、軸対称な角膜表面方程式を非同次分数階ベッセル方程式に変換する。
• 同次部分をべき級数展開を用いて解き、古典的I₀がα = 1のときに還元される一般化関数I₀^α(x)を導出する。
• 特解（定数）と同次解の和として一般解を導出し、境界条件を適用する。
• ラプラス法と漸近的解析を用いて、I₀^α(x)のxが大きいときの挙動を特定し、指数関数的増加に累乗則係数が乗じられた形であることを示す。
• 最小二乗法を用いて123×123の実際の角膜測定データにモデルを適合させ、前後表面それぞれについてパラメータa, b, αを推定する。

実験結果

リサーチクエスチョン

• RQ1 分数階ベッセル微分方程式は、古典的モデルに比べて、角膜トポグラフィーの記述においてより物理的整合性があり、より高い正確性を提供できるか？
• RQ2 径方向変数x → ∞のとき、修正分数階ベッセル方程式の解はどのように漸近的挙動を示すか？
• RQ3 分数階数αを変化させた場合、角膜表面の形状および実データへのモデル適合度にどのような影響を与えるか？
• RQ4 モデルのパラメータa, b, αは角膜の物理的性質とどのように関係しており、期待されるオーダーの大きさであるか？
• RQ5 提案されたモデルは、実際の角膜トポグラフィー測定データに対してミリメートル未満の適合精度を達成できるか？

主な発見

• 分数階ベッセル方程式の解は絶対収束べき級数として表現され、古典的修正ベッセル関数I₀を一般化する新しい関数I₀^α(x)を定義する。
• x → ∞におけるI₀^α(x)の漸近的挙動は、exp(x^{2/(1+α)} / (2/(1+α)))に累乗則係数を乗じた形であることが示され、α = 0およびα = 1のときには既知の漸近的形に還元される。
• 外側角膜表面に対して、最適パラメータa = 0.580404, b = 1.19734, α = 0.421345を用いた場合、平均絶対誤差は0.014 mmに達する。
• 内側角膜表面に対しては、平均絶対誤差が0.03 mmであり、パラメータはa = 0.818763, b = 1.66664, α = 0.503431である。
• 分数階数αは0と1の間にあるため、非整数階微分の順序が、古典的モデルでは捉えきれない中間的な物理的挙動を記述できることを示している。
• モデルのパラメータaとbはオーダー1の大きさであり、事前に提唱された膜モデルからの理論的期待と整合しており、本手法の物理的妥当性を裏付けている。