双正交转移矩阵重整群於广义马可夫随机场之数值模拟Numerical.PDF

南榮學報第十 八期民國 104 年 7 月 Journal of Nan Jeon, Vol. 1 8, 20 15 雙正交轉移矩陣重整群 於廣義馬可夫隨機場之數值模擬 Numerical Simulations of General Markov Random Fields Using Biorthonormal Transfer-Matrix Renormalization-Group Method 黃玉坤 彭志豪 楊竣欽 Yu-Kun Huang Chih-Hao Peng Chun-Chen Yang 南榮科技大學工程科技研究所 摘 要 本研究使用雙正交轉移矩陣重整群 (Biorthonormal Transfer-Matrix Renormalization Group; BTMRG)方法探討廣義馬可夫隨機場 (Markov random field; MRF配分函數的數值計算。) BTMRG 乃 是改良自傳統轉移矩陣重整群 (TMRG) ，而TMRG的基本概念則是繼承自密度矩陣重整群 (Density-Matrix Renormalization Group; DMRG)並應用於non-Hermitian轉移矩陣問題 ，可在大型晶 格系統中精確計算系統的配分函數，不但有自動截取(decimate) Hilbert 空間的特性，同時控制重整化 的轉移矩陣維度在適度範圍 。本文中，我們將應用BTMRG 來對廣義MRF的基本轉移矩陣加以截取 。 我們選擇四個二元二階MRF 進行 BTMRG的數值模擬。其模擬結果皆呈現出相變 (phase transition) 的現象。本研究證實我們的 BTMRG 不論在精確度或運算速度上均優於傳統 TMRG ，並且較傳統 TMRG適用於大型晶格系統。 關鍵字：馬可夫隨機場 、配分函數 、密度矩陣重整群、轉移矩陣重整群 。 一、前言 馬可夫隨機場 (Markov Random Field; MRF)[1]是著名的統計 力學晶格模型。它在 各種領域中都有得到成功應用，比如說統計物理學，編碼理論，機器學習和影像建模。 例如在統計力學上，著名的 Ising模型和 Potts模型乃是馬可夫隨機場中的一個特例。 儘管已獲得如此廣泛而成功的應用，在 MRF研究中，配分函數 (partition function)的數 值計算仍然具相當高的挑戰性。在統計力學上， MRF配分函數的計算是極其基本而 重要的。例如， MRF的自由能 (free energy)跟 MRF配分函數的對數值有關，而一旦有 了自由能，幾乎所有熱力學屬性都可以從自由能當中經由標準化的推導過程獲得。因 為想要獲得廣義 MRF配分函數的解析解形式非常困難，所以實務上大多採取數值計 算來 求得配分函數的近似值，而在這些數值方法中，最重要的、也是最受歡迎的就是 B2-1 雙正交轉移矩陣重整群於廣義馬可夫隨機場之數值模擬 所謂的轉移矩陣法了。在這方法中， MRF的自由能是由一個非負轉移矩陣之 Perron root （即最大本徵值）的對數值所推導出來的。不幸的是，轉移矩陣的維度會隨著晶 格的尺寸大小呈現出幾何級數方式成長，因此這種方法會嚴重地侷限於小型系統。在 通常的情況下，我們可以利用轉移矩陣重整群 (Transfer-Matrix Renormalization Group; TMRG)[2,11,12,22]來避開這個轉移矩陣維度無限大的難題。TMRG 乃是繼承密度矩陣 重整群 (Density-Matrix Renormalization Group; DMRG)[4,14,17,20]的基本概念，在 Hilbert 空間中作系統性的截取(truncation)