[討論] 矩陣解聯立方程式機率算出來有負的 ><已回收

看板MATLAB作者vigirl (小潔寶)時間14年前 (2011/04/20 10:10)推噓2(2推 0噓 12→)

留言14則, 5人參與討論串1/2 (看更多)

最近在跑一個演算法跑出來可以得到一個N*N的矩陣裡面都是機率代表狀態之間的轉移機率接著我要利用解聯立的方法求出各個狀態的穩定機率 ex. 我的轉移機率矩陣 a=[0.1, 0.2, 0.3; 0.4, 0.5, 0.6; 0.7, 0.8, 0.9] 所以我有一些聯立方程式： 0.1*x+0.4*y+0.7*z=x 0.2*x+0.5*y+0.8*z=y 0.3*x+0.6*y+0,9*z=z 另外多加一個條件（機率總合為1）：x+y+z=1 我假設h=[x;y;z]是我欲求的穩定機率接著我把等號右邊x,y,z移項到左邊： e=eye(3) f=[1,1,1] d=[0,0,0,1] %最右邊的1是x+y+z="1"的1,前三個0是移項後等號右邊的0 於是c=[a'-e;f]; c=[0.1-1, 0.4, 0.7; 0.2, 0.5-1, 0.8; 0.3 0.6 0.9-1; 1 1 1 ] 然後解出 h=c\d; 但我發現我怎麼算因為我實際上矩陣a裡有幾千個機率值所以很難每一個每一個式子去檢查 >< 我算出來的h矩陣，也就是各個狀態的穩定機率竟然都會有負的?! 但機率根本不可能有負的呀... 請問有沒有人可以發現到我哪裡有不對的地方呢？先謝謝大家！ -- ╭。☆ http://www.wretch.cc/album/pinkyvigirl ☆。╮ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.117.164.26

→

CCWck

04/20 10:30, , 1^F

04/20 10:30, 1^F

推

Haderson

04/20 12:30, , 2^F

04/20 12:30, 2^F

→

verb

04/20 13:13, , 3^F

04/20 13:13, 3^F