[论文解读] 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

一句话总结

本文建立了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)的时间依赖性电网络，其中σ(t)通过dσij/dt + σij = ϕij进行更新。
• 利用基尔霍夫定律，从σ(t)和固定外部源向量b计算电流ϕ(t)和电势p(t)。
• 通过在子图Hr上基于归纳法的证明，展示电势在支持集逐步扩展过程中收敛至p*。
• 运用最大值原理和离散∞-调和函数理论，分析在扩展子图上的电势收敛性。
• 采用对数变换ξ(t) = log(σ(t)/σ(0))，分析不同弧集上电导率的渐近行为。
• 利用引理5比较线性拉普拉斯系统解，建立电势误差的指数衰减。

实验结果

研究问题

• RQ1 Physarum求解器是否在一般有向图上收敛至线性转送问题的最优解？
• RQ2 流量ϕ(t)和电导率σ(t)收敛至最优解的速率如何？
• RQ3 电势p(t)是否收敛至对偶问题的解？若是，其数学特征是什么？
• RQ4 该收敛性是否可证明适用于一般有向图，而不仅限于平面图或最短路径情形？
• RQ5 最优支撑集H*的结构如何影响电势和流量的收敛行为？

主要发现

• 电流ϕ(t)以指数速率收敛至最小化∑ℓijϕij的最优流量ϕ̂。
• 电导率σ(t)收敛至一个正极限，且其收敛也是指数的。
• 电势p(t)收敛至对偶问题的解p*，该解是在最优支撑H*上的离散∞-调和函数。
• 若最优支撑H*连通，则电势收敛至在图G的生成子图上定义的规范对偶解。
• 电势的收敛速率是指数的，其速率满足递推关系a_{k+1} = min{a_k, r - r'}，适用于连续的子图扩展。
• 分析表明，该算法能自然识别最优子图，并在扩展至新弧时仍保持指数收敛。