G=[1/16, 1/6,1/8;1/32,0,1/4;1/32,1/4,1/12];%轉移概率a=20;%H=zeros(a);%[r,s]=size(H);%TotalNum=r*s;%IndexNum=1+floor(rand(1,floor(TotalNum*0.1))*TotalNum);%H(IndexNum)=1;h2=zeros(a*a); %轉移概率矩陣分佈x=1;for x=1:1:a*aif x==1h2((x-1)*a*a+2)=3/7;h2((x-1)*a*a+a+1)=3/7;h2((x-1)*a*a+a+2)=1/7;elseif x==a*ah2((x-1)*a*a+a*a-1)=16/25;h2((x-1)*a*a+a*(a-1)-1)=6/25;h2((x-1)*a*a+a*(a-1))=3/25;elseif x==ah2((x-1)*a*a+a-1)=4/13;h2((x-1)*a*a+2*a-1)=3/13;h2((x-1)*a*a+2*a)=6/13;elseif x==a*(a-1)+1h2((x-1)*a*a+a*(a-2)+1)=1/10;h2((x-1)*a*a+a*(a-2)+2)=1/10;h2((x-1)*a*a+a*(a-1)+2)=4/5;elseif (1<x)&&(x<a)h2((x-1)*a*a+x-1)=4/21;h2((x-1)*a*a+x+1)=6/21;h2((x-1)*a*a+x+a-1)=1/7;h2((x-1)*a*a+x+a)=2/7;h2((x-1)*a*a+x+a+1)=2/21;elseif mod((x-1),a)==0 && x~=a*(a-1)+1h2((x-1)*a*a+x-a)=3/62;h2((x-1)*a*a+x-a+1)=3/62;h2((x-1)*a*a+x+1)=12/31;h2((x-1)*a*a+x+a)=12/31;h2((x-1)*a*a+x+a+1)=4/31;elseif mod(x,a)==0 && x~=a*ah2((x-1)*a*a+x-a)=6/61;h2((x-1)*a*a+x-a-1)=6/61;h2((x-1)*a*a+x-1)=16/61;h2((x-1)*a*a+x+a)=24/61;h2((x-1)*a*a+x+a-1)=12/61;elseif (a*(a-1)+1<x)&&(x<a*a)h2((x-1)*a*a+x-1)=16/52;h2((x-1)*a*a+x+1)=24/52;h2((x-1)*a*a+x-a)=3/52;h2((x-1)*a*a+x-a-1)=6/52;h2((x-1)*a*a+x-a+1)=3/52;elseh2((x-1)*a*a+x-1)=1/6;h2((x-1)*a*a+x+1)=1/4;h2((x-1)*a*a+x-a-1)=1/16;h2((x-1)*a*a+x-a)=1/32;h2((x-1)*a*a+x-a+1)=1/32;h2((x-1)*a*a+x+a-1)=1/8;h2((x-1)*a*a+x+a)=1/4;h2((x-1)*a*a+x+a+1)=1/12;endend%h0=zeros(1,a*a);%h0=zeros(a);%[r,s]=size(h0);%TotalNum=r*s;%IndexNum=1+floor(rand(1,floor(TotalNum*0.2))*TotalNum);%H(IndexNum)=1;h0_1=rand(1,a*a);H0_1=sum(h0_1);h0=h0_1/H0_1;k=20;for i=1:k; S=h0*h2^k; S0 = reshape(S,a,a); T=S'; S0_max=max(S0(:));%求出矩陣中最大數[x,y]=find(S0==max(S0(:)));%求出矩陣中最大數所在的位置goalIndex=(y-1)*a+x;%矩陣中最大自所在的柵格序號 endm=[h2'-eye(a*a);ones(1,a*a)];n=[zeros(a*a,1);1];h2_limit=m\n;plot(T,'m') hold on plot(h2_limit,'g')
馬爾科夫鏈-轉移概率
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章
「 網絡安全術語解讀 」通用平臺枚舉CPE詳解
https://blog.csdn.net/u013129300/article/details/129329786
規格嚴格-功夫到家
2024-05-05 13:43:19
深入學習和理解Django視圖層:處理請求與響應
Mifen
2024-05-05 13:38:28
如何閱讀 Paper
mghio
2024-05-05 13:32:48
REACT: SYNERGIZING REASONING AND ACTING IN LANGUAGE MODELS
initial_h
2024-05-05 13:32:27
sysbench的部分基準性能測試學習
濟南小老虎
2024-05-05 13:29:27
[轉帖]sysbench基準測試
濟南小老虎
2024-05-05 13:29:27
C++ lambda的重載
apocelipes
2024-05-05 13:22:37
關於農業種植的若干思考
LittleHann
2024-05-05 13:16:36
爲什麼有的人說技術不重要,有的人說重要?
killianxu
2024-05-05 13:10:35
博客性能優化筆記 | 99分
粥裏有勺糖
2024-05-05 13:06:15
LLM2Vec介紹和將Llama 3轉換爲嵌入模型代碼示例
deephub
2024-05-05 13:04:45
開發Android應用程序,在Android10的系統上提示網絡出錯?
LoveCoder
2024-05-05 13:02:54