[論文レビュー] Stability Analysis of Fractional Differential Equations with Unknown Parameters

Mehmet Emir Köksal|arXiv (Cornell University)|Sep 15, 2017

Fractional Differential Equations Solutions参考文献 43被引用数 18

ひとこと要約

本稿では、未知のパラメータを有する分数階微分方程式（FDEs）のパrametric安定性解析のための図式的D-分解法を提案する。複雑な解析的計算を回避し、パrametric空間における実根境界、複素根境界、無限大根境界を導出することで、視覚的に安定性領域を特定し、広範なパラメータ区間におけるロバストな安定性評価を可能にする。本手法はベンチマークFDEsを用いて検証され、分數階動的を有する工学的システムにおける信頼性の高い実用的安定性解析を示している。

ABSTRACT

In this paper, the stability of fractional differential equations (FDEs) with unknown parameters is studied. FDEs bring many advantages to model the physical systems in the nature or man-made systems in the industry. Because this representation has a property between linear differential equations and nonlinear differential equations. Therefore, the designer may use the FDEs to model complex systems instead of nonlinear differential equations which have hard mathematical background. Using the graphical based D-decomposition method, we investigate the parametric stability analysis of FDEs without complicated mathematical analysis. To achieve this, stability boundaries are obtained firstly, and then the stability region set depending on the unknown parameters is found. The applicability of the presented method is shown considering some benchmark equations which are often used to verify the results of a new method. Simulation examples shown that the method is simple and give reliable stability results.

研究の動機と目的

• 未知または変動するパラメータを有する分数階微分方程式（FDEs）における安定性解析の課題に取り組む。これは古典的安定性解析よりも複雑である。
• 既存手法がパラメータの不確実性を小さな区間に制限しているという限界を克服し、ゼロから無限大までの全範囲で解析を可能にする。
• エンジニアや研究者が複雑な解析的解に依存せずに、安定性を評価できる実用的で視覚的かつ計算効率の良い手法を開発する。
• 共通順序および非共通順序FDEsの両方について、パラメトリック空間における安定性領域を特定する体系的なフレームワークを提供する。
• パラメータの変動に伴う安定性の変化を視覚化することで、分數階システムの設計および制御を支援するロバストネス解析を可能にする。

提案手法

• D-分解法を適用し、s平面における安定性境界（虚軸）をパラメトリックドメインにおける3種類の境界へ写像する：実根境界、複素根境界、無限大根境界。
• FDE伝達関数の特性準多項式に基づき、CaputoまたはGrünwald-Letnikovの分數階微分定義を用いて、安定性境界の解析的表現を導出する。
• Matignonの安定性基準を基盤とし、すべての固有値λiについて|arg(−λi)| > απ/2を満たす条件を、パラメトリック空間における安定性条件として定義する。
• 可能な限りFDE伝達関数を共通順序形式に変換し、三角恒等式を用いて境界曲線の導出を簡略化する。
• パrameterを変化させながら導出した境界方程式（例：a = −bω−α sin(0.5πα)/sin(πα), c = −aω²α cos(πα) − bωα cos(0.5πα)）をプロットすることで、安定性領域を生成する。
• 2次元または3次元プロット（必要に応じて色を4番目の次元として使用）を用いて、パrameterの変化に伴う安定性領域の変化を可視化する。

実験結果

リサーチクエスチョン

• RQ1 複雑な解析的手法に依存せずに、未知のパラメータを有する分數階微分方程式の安定性をどのように解析できるか？
• RQ2 変動する係数を有するFDEsにおいて、安定領域と不安定領域を分かつパラメトリック境界は何か？
• RQ3 減衰係数や剛性係数（a, b, c）が広範囲（正負を含む）に変動する場合、安定性領域はどのように変化するか？
• RQ4 提案手法が、Basset方程式や産業用加熱炉モデルなどのベンチマークFDEsについて、既知の安定性結果を的確に検証および拡張できるか？
• RQ5 共通順序と非共通順序の分數階微分の違いが、安定性領域の形状およびロバストネスに与える影響は何か？

主な発見

• 提案されたD-分解法は、複雑な解析的解に依存せずに、図式的な可視化によって未知パラメータを有するFDEsの安定性領域を的確に特定できる。
• α = 0.5の共通順序FDEでは、複素根境界が対称性を示す。すなわち、(a, c)が曲線上にあるならば、(c, a)も同様に曲線上にある。これは一般の非共通順序では成り立たない性質である。
• 例2では、α = 0.2およびb ∈ [−5, −1]の安定性領域が図8に可視化されており、αが増加するにつれて領域が収縮し、α → 0のときには右上象限が埋まることが示されている。
• α = 0.8の場合、bを−bに置き換えると、安定性領域はα = 0.2の場合と原点に関して対称となる。これは本手法の整合性を裏付けている。
• 例3では、産業用加熱炉の名目パラメータ（a = 14994, b = 6009.5, c = 1.69）が計算された安定性領域内に位置しており、追加計算なしにSondhiとHoteの安定性主張を検証している。
• 図12の3次元安定性領域は、すべての正のaとc、およびすべての負のaとcが、bの符号に応じて楕円形領域内で安定性を保つことを示しており、広範なパラメータ範囲にわたるロバストネスを示している。