[論文レビュー] Convergence Properties for the Physarum Solver

Kentaro Ito, Anders Johansson+2 more|arXiv (Cornell University)|Jan 27, 2011

Slime Mold and Myxomycetes Research参考文献 30被引用数 27

ひとこと要約

この論文は、線形輸送問題を解く生物学的にインspiredなアルゴリズムであるPhysarumソルバの収束を確立し、流量および伝導度が最適解へ指数関数的に収束することを示している。また、電気的ポテンシャルが∞-調和関数に収束することを証明し、このアルゴリズムのネットワーク最適化における性能に対する厳密な数学的基盤を提供している。

ABSTRACT

The Physarum solver is an intuitive mechanism for solving optimisation problems based on the idea of an electrical network, whereby the conductivity is reinforced by the current. We show that the Physarum solver obtains the solution to the linear transshipment problem on a digraph and that the electrical potential converges to an infinity-harmonic solution of the dual problem.

研究の動機と目的

• 有向グラフ上での線形輸送問題を解くためのPhysarumソルバの収束特性を厳密に分析すること。
• コスト関数∑ℓijϕijを最小化する最適解への流量および伝導度ベクトルの指数的収束を確立すること。
• 電気的ポテンシャルが最適化問題の双対解に対応する離散∞-調和関数に収束することを証明すること。
• 平面グラフにおける最短経路収束に関する先行結果を、任意の重みを持つ一般の有向グラフへ一般化すること。
• アルゴリズムの生物学的妥当性および分散型計算への応用可能性のための数学的基盤を提供すること。

提案手法

• Physarumソルバを時間依存の電気回路としてモデル化し、時間変動する伝導度σ(t)をdσij/dt + σij = ϕijに従って更新する。
• キルヒホッフの法則を用いて、固定された外部源ベクトルbとσ(t)から電流ϕ(t)およびポテンシャルp(t)を計算する。
• 部分グラフHr上で帰納的証明を適用し、増大するサポート上でのポテンシャルの指数的収束を示す。
• 最大原理および離散∞-調和関数論を用いて、拡張された部分グラフ上でのポテンシャル収束を分析する。
• 伝導度の漸近的挙動を分析するため、対数変換ξ(t) = log(σ(t)/σ(0))を用いる。
• 補題5を活用して線形ラプラシアン方程式系の解を比較し、ポテンシャル誤差の指数的減衰を確立する。

実験結果

リサーチクエスチョン

• RQ1 Physarumソルバは一般の有向グラフ上での線形輸送問題の最適解へ収束するか？
• RQ2 流量ϕ(t)および伝導度σ(t)の最適解への収束速度は何か？
• RQ3 電気的ポテンシャルp(t)は双対問題の解に収束するか？ もしそうならば、その数学的特徴は何か？
• RQ4 収束は平面グラフや最短経路のケースに限らず、一般の有向グラフに対しても証明可能か？
• RQ5 最適サポート集合H*の構造は、ポテンシャルおよび流量の収束挙動にどのように影響するか？

主な発見

• 電流ϕ(t)は、∑ℓijϕijを最小化する最適流量ϕ̂へ指数的レートで収束する。
• 伝導度σ(t)は正の極限値へ収束し、その収束も指数的である。
• 電気的ポテンシャルp(t)は、最適サポートH*上で定義された離散∞-調和関数である双対問題の解p*に収束する。
• 最適サポートH*が連結であれば、ポテンシャルはGのスパニング部分グラフ上で定義された標準的な双対解に収束する。
• ポテンシャルの収束レートは指数的であり、連続する部分グラフ拡張においてa_{k+1} = min{a_k, r - r'}と表される。
• 解析により、アルゴリズムが自然に最適部分グラフを特定し、新しいアークへ拡張しても指数的収束を維持することが示された。