[論文レビュー] Analysis of a time-dependent problem of mixed migration and population dynamics

François Castella, Philippe Chartier+1 more|arXiv (Cornell University)|Dec 7, 2015

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 5被引用数 3

ひとこと要約

本稿では、速い空間的移動と遅い人口動態的相互作用を有する時間依存型の被捕食者-被捕食者モデルを分析するための数学的枠組みを構築する。一般化された中心多様体定理と平均化法を組み合わせた手法を用い、主な結果として、システムが平均化係数を有する非周期的ロトカ＝ヴォリエーラ系に還元されることを示す。さらに高次の解析により、速い振動が極小周期軌道の不安定化を引き起こすことが明確に特定される。

ABSTRACT

In this work, we consider a system of differential equations modeling the dynamics of some populations of preys and predators, moving in space according to rapidly oscillating time-dependent transport terms, and interacting with each other through a Lotka-Volterra term. These two contributions naturally induce two separated time-scales in the problem. A generalized center manifold theorem is derived to handle the situation where the linear terms are depending on the fast time in a periodic way. The resulting equations are then amenable to averaging methods. As a product of these combined techniques, one obtains an autonomous differential system in reduced dimension whose dynamics can be analyzed in a much simpler way as compared to original equations. Strikingly enough, this system is of Lotka-Volterra form with modified coefficients. Besides, a higher order perturbation analysis allows to show that the oscillations on the original model destabilize the cycles of the averaged Volterra system in a way that can be explicitely computed.

研究の動機と目的

• 速い移動と遅い人口動態的相互作用という2つの時間スケールを持つ、捕食者-被捕食者ダイナミクスをモデル化する時間依存型微分方程式系の解析。
• 速いつきの時間に依存する周期的で急激に変動する移動項が、標準的な平均化法の仮定を破るという課題に対処すること。
• 元の高周波数のモデルの本質的な長期的挙動を捉える、簡略化された非周期的力学系への還元。
• 高次の摂動解析を用いて、平均化系における極小周期軌道の安定性に及ぼす速いつきの振動の影響を定量化すること。
• 2つの空間的サイトを含む数値例を用いて理論的枠組みの妥当性を検証し、収束性と誤差スケーリングを示すこと。

提案手法

• 移動と人口動態の時間スケール比を表す小パラメータ ε を用いた常微分方程式系を定式化する。
• 一般化された中心多様体定理を適用し、速いつき時間変数 θ = t/ε に周期的に依存する遅い多様体を構築する。
• 偏微分方程式の解を用いて、ε の任意の次数まで正確に近似可能な中心多様体の再帰的公式を導出する。
• 速いつきの振動を除去するための周期的変数変換を実行し、平均化法を適用して簡略化された非周期的力学系を導出する。
• 数値積分（RK4）と陰的エイラー法を用いて近似の妥当性を検証し、正確な解と近似解を比較する。
• 対数-対数スケールでの誤差解析により、1次中心多様体近似の収束率が O(ε) であることを確認する。

実験結果

リサーチクエスチョン

• RQ1 速いつきの周期的移動と遅い人口動態を有する系は、どのように簡略化された非周期的力学系に還元できるか？
• RQ2 線形項が速いつき時間変数に周期的に依存する場合、中心多様体の構造はどのようなものか？
• RQ3 移動率の速いつきの振動は、平均化されたロトカ＝ヴォリエーラ系における極小周期軌道の安定性にどのように影響するか？
• RQ4 再帰的関係を用いて、ε の任意の次数まで中心多様体を明示的に計算できるか？
• RQ5 小パラメータ ε の観点から、1次中心多様体近似の収束率はどのようになるか？

主な発見

• 還元された系は、元のシステムのパラメータの空間的・時間的平均を係数に持つロトカ＝ヴォリエーラ形式をとる。
• 中心多様体は速いつき時間変数に周期的に依存し、偏微分方程式の再帰的解を用いて任意の次数まで近似可能である。
• 高次の摂動解析により、速いつきの振動が平均化系の極小周期軌道を不安定化することが判明し、その不安定化効果は明示的に計算可能である。
• 数値的検証により、1次中心多様体近似の誤差が O(ε) のオーダーであることが確認され、対数-対数誤差プロットにおける捕食者と被捕食者の勾配はそれぞれ 1.0386 と 1.0080 であった。
• 初期条件を中心多様体上に選んだ場合、誤差は t/ε に対して指数関数的に減少し、O(ε) に安定化するため、近似の漸近的妥当性が裏付けられる。
• 本手法は2サイトモデルのダイナミクスを効果的に捉えており、数値結果は解軌道および誤差挙動の両面で理論的予測と一致している。