![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
柴原?
PSO（Particle Swarm Optimization）技術は、自然界の群れ行動に触発された最適化手法です。この手法では、複数の「粒子」（鳥）によって解の探索が行われます。
各粒子の位置は候補解を示しており、時間ごとに各粒子は自分の現在の位置でコスト関数を評価します。その評価値と、粒子がこれまでに経験した中で最も良かった位置（自己ベスト）を比較します。
各粒子は、自分の自己ベストを保持しており、他の粒子からの情報を交換することで、群れ全体の最良の位置（グローバルベスト）を特定します。粒子の移動の方向は、以下の2つの要素に基づいています：
1. 社会的要素：グローバルベストに向かって引っ張る成分 2. 認知的要素：自己ベストに向かって引っ張る成分
各粒子の速度は、以下の3つの成分から構成されています：
PSOの更新方程式は次のように表されます：
速度の更新： v_i(k+1) = α_d * v_i(k) + β_d * (x_lbi - x_i(k)) + γ_d * (x_gb - x_i(k))
位置の更新： x_i(k+1) = x_i(k) + v_i(k+1)
ここで、
基本的なPSOアルゴリズムは離散時間の性質を持ち、各粒子は単純な動きの方程式を持っています。これらの方程式は、群れ全体をコンパクトに表現することも可能です。
問題の次元を d、実行可能領域を Ω ⊆ ℝ^d、鳥の数を n、最小化する関数を f: Ω → ℝ と定義しています。
このようにして、PSOは群れの知識を活用し、最適解に向かって効率的に探索を行う手法です。
PSO（Particle Swarm Optimization）の新しい表現を理解するために、以下の要素を見ていきます。
位置行列 X: X = [ x_1 ⋮ x_n ] ∈ (Ω × ℝ^n) これは全ての粒子の位置を含む行列です。
速度行列 V: V = [ v_1 ⋮ v_n ] ∈ (ℝ^d × ℝ^n) これは全ての粒子の速度を含む行列です。
自己ベスト位置行列 X_lb: X_lb = [ x_lb1 ⋮ x_lbn ] ∈ (Ω × ℝ^n) 各粒子の自己ベスト位置を保持します。
グローバルベスト位置 X_gb: X_gb ∈ ℝ^d これは群れ全体の最良の位置を示します。
目的関数ベクトル F: F = [ f(x_1) ⋮ f(x_n) ] : (Ω × ℝ^n) → ℝ^n 各粒子の位置に対する目的関数の値を保持します。
ベクトル T: T ∈ ℝ^d 全ての要素が1で構成される行ベクトルです。
ベクトル Q_i: Q_i ∈ ℝ^n すべての要素が0で、i番目の要素のみが1の列ベクトルです。
単位行列 I_n: I_n ∈ (ℝ^n × ℝ^n) サイズnの単位行列です。
これらの定義を用いて、PSOアルゴリズムの更新式は次のように表されます：
速度の更新:
V(k+1) = α_d * V(k) + β_d * (X_lb(k) - 
k)) + γ_d * (X_gb(k) * T - k))
位置の更新:
X(k+1) = k) + V(k+1)
自己ベストの更新: X_lb(k+1) = 1/2 * (X(k+1) + X_lb(k))
自己ベストの改善: X_lb(k+1) = X_lb(k) + [X(k+1) - X_lb(k)] * ζ
ここで、ζは次のように定義されます： ζ = diag[sgn(F(X_lb(k)) - F(X(k+1)))]
このdiag[y]は、ベクトルyの要素を対角成分とする対角行列です。また、sgn(⋅)は符号関数で、次のように定義されます： sgn(y) = { 1, y ≥ 0; -1, y < 0 }
これにより、PSOアルゴリズムの全体的なフレームワークが整い、効率的に最適解を探索することができます。
この部分は、PSO（Particle Swarm Optimization）アルゴリズムにおける自己ベスト位置の更新方法を説明しています。具体的には、粒子の位置が改善されたかどうかに基づいて自己ベストを更新するための計算です。
自己ベスト位置の更新 自己ベスト位置 X_lb(k+1) は、次のように更新されます： X_lb(k+1) = X_lb(k) + [X(k+1) - X_lb(k)] * ζ
ここで、ζ は以下のように定義されます： ζ = diag[sgn(F(X_lb(k)) - F(X(k+1)))]
各要素の解説
自己ベストの更新式:
X_lb(k+1) は、前の自己ベスト X_lb(k) に新しい位置 k+1) から前の自己ベスト X_lb(k) を引いたベクトルに、ζ を掛けたものを足すことで計算されます。この計算によって、自己ベストが新しい位置に基づいて更新されますが、更新の影響は ζ によって制御されます。
ζ の定義: ζ は対角行列で、各対角要素は次のように定義される符号関数 sgn によって決まります。
符号関数 sgn(y) は、次のように定義されます： sgn(y) = { 1, if y ≥ 0; -1, if y < 0 }
つまり、F(X_lb(k))（自己ベスト位置での目的関数の値）と F(X(k+1))（新しい位置での目的関数の値）の差が0以上であれば、その成分は1になり、負であれば-1になります。
更新の意義: ζ が1の場合、その位置は改善されているため、自己ベストが更新されます。逆に、ζ が-1の場合、その位置は自己ベストより悪いため、自己ベストはそのまま保持されます。この仕組みによって、粒子は自己ベストを賢く更新し、最適解に向かって進むことができます。
まとめ: この更新メカニズムは、粒子が探索する際に過去の良い経験を活かしつつ、最新の情報を反映させるための重要なステップです。PSOアルゴリズム全体の効率を高める役割を果たしています。
