Aritalab:Lecture/NetworkBiology/Percolation on Graph

Revision as of 21:39, 15 June 2011

グラフ上のパーコレーション

グラフをある次数分布に従うツリーとみなすことで、パーコレーションを扱ってみます。この内容を読む前に母関数の概念を復習しておいてください。

次数分布の母関数

グラフの次数分布 $p(k)$ を与えられたとき、その母関数は

$G_v(x) = \sum p(k) x^k$

になります。しかし、頂点 v の隣にある頂点 w の次数分布はもはや $p(k)$ ではありません。選ばれた辺の先には、他の点が等確率でくることはないからです。辺をランダムに1本選んだとき、もう片方にある端点がどんな次数分布をとるか考えましょう。次数の大きい点のほうが、次数に比例して存在する確率が高くなります。つまり、もう片側にある端点の次数分布は、最初の次数分布に辺の数で重みをつけてから正規化した $\textstyle kp(k)/\Sigma_j jp(j)$ となります（分母は正規化定数）。よってある頂点 v の隣にくる頂点 w の次数分布の母関数は

$\sum_{k=0}^{\infty}x^k\frac{kp(k)}{\Sigma_jjp(j)} = x\frac{\Sigma_k kx^{(k-1)}p(k)}{\Sigma_jjp(j)} = x\frac{G'_v(x)}{G'_v(1)}$

母関数においては、次数 $k$ に対して $x^k$ が対応します。頂点 v の先に w がついているとき、w から v 以外に伸びる辺数に注目しましょう。 v から w への辺1本ぶんを除いた w に関する母関数は、上の式を $x$ で割ったものにあたります。

$G_w(x) = \frac{G'_v(x)}{G'_v(1)}$

占有される頂点の母関数

パーコレーションで相転移が起きるのは、頂点から伸びていく辺の期待値が1をちょうど超えるときです。ただし、このときの母関数はグラフ全体の次数分布を示す $G_v$ ではなく、占有された頂点のみによる部分ネットワークの母関数です。したがって、母関数を区別して $F_v$ と記しましょう。関数 F と G の関係は、各頂点が確率 q で占有されるため G と F の次数平均をそれぞれ $\langle k \rangle$ , $\langle l \rangle$ と書くと $\begin{align} \langle l \rangle &= q \langle k \rangle \\ \end{align}$ です。

この関係から　(式の導出は確率母関数のページを参照）

$\begin{align} F'_v(1) &= q G'_v(1) = q \langle k \rangle \\ F''_v(1) &= q^2G''_v(1) = q^2 (\langle k^2 \rangle - \langle k \rangle) \end{align}$

相転移が起きるのは $F_w'(1) = 1$ のときだから

$F'_w(1) = \frac{F''_v(1)}{F'_v(1)} = 1$

代入すると

$q_c = \frac{1}{\frac{\langle k^2 \rangle}{\langle k \rangle} - 1}$

すなわち、 $\langle k^2 \rangle$ が大きく $\langle k \rangle$ を上回るときに　 $q_c$ は小さくなります。ハブの k は大きくなりますから、 $k^2$ は更に大きくなります。次数の大きなハブがあるほど、パーコレーションが起きやすいということです。

パーコレーション

いかなる次数分布でも、 $\langle k^2 \rangle$ が $\langle k \rangle$ に比して大きいと、低い浸透確率で相転移が起こることをみてきました。

ランダムグラフの場合

ランダムグラフの場合は次数分布がポアソン分布になります。このとき $\langle k^2 \rangle = \langle k \rangle^2 + \langle k \rangle$ になるので

$q_c = 1 / \langle k \rangle$

すなわち、平均次数が大きくなるほど浸透しやすくなります。

スケールフリーネットワークの場合

べき分布 $p(k)=ck^{-\gamma}$ のとき

$\langle k^2 \rangle = \Sigma_kk^2p(k) = c\Sigma k^{2-\gamma}$

この値はFailed to parse (lexing error): \gamma　\leq 3 のときに発散します。有限のネットワークではその値も有限ですが、大きな値になるために非常に低い値で浸透すると考えてよさそうです。

@@ Line 39: / Line 39: @@
 です。
-この関係から
+この関係から　(式の導出は[[Aritalab:Lecture/Basic/Probability_Generating_Function|確率母関数のページ]]を参照）
 <math>
@@ Line 60: / Line 60: @@
 </math>
-すなわち、<math>\langle k^2 \rangle</math>が大きく<math>\langle k \rangle</math>を上回るとき、<math>q_c</math>は小さくなります。
+すなわち、<math>\langle k^2 \rangle</math>が大きく <math>\langle k \rangle</math> を上回るときに　<math>q_c</math> は小さくなります。
+ハブの ''k'' は大きくなりますから、<math>k^2</math> は更に大きくなります。次数の大きなハブがあるほど、パーコレーションが起きやすいということです。
 ==パーコレーション==

Aritalab:Lecture/NetworkBiology/Percolation on Graph

Revision as of 21:39, 15 June 2011

Contents

グラフ上のパーコレーション

次数分布の母関数

占有される頂点の母関数

パーコレーション

ランダムグラフの場合

スケールフリーネットワークの場合

Personal tools

Namespaces

Variants

Views

Actions

Search

Navigation

metabolites

Toolbox