function?multidistribution=test（）
Qmatrix=input（‘請輸入速率轉移矩陣例如[-2?1?1;1?-4?3;2?2?-4]\n‘）;
Pmatrix=input（‘請輸入初始概率分布矩陣例如[0.5;0.2;0.3]\n‘）;
Tmatrix=[];
%求出Tmatrix，用來累計i分布
tmp=0;
for?i=1:1:length（Qmatrix）
????tmp=0;
????for?j=1:1:length（Qmatrix）
????????if（i==j）
????????????tmp=tmp+0;
????????elseif（i~=j）
????????????tmp=tmp+Qmatrix（ij）;
????????end
????????Tmatrix（ij）=tmp/（-Qmatrix（ii））;
????end?????
end
multidistribution=[];
%alltimedistribution=[];獲取所有時間列表
for?n=1:1:6
????%隨機產生一個變量，確定從哪一個狀態開始
????r=unifrnd（01）;
????tmp=0;
????for?i=1:1:length（Pmatrix）
????????tmp=tmp+Pmatrix（i）;
????????if?r>tmp
????????????continue;
????????elseif?r????????????break;
????????end
????end
????%此時的i是初始狀態
????%用一個矩陣來記錄下狀態的變化和在某一狀態停留的時間
????totaltime=0;
????timedistribution=linspace（00length（Pmatrix））;%用來記錄各個狀態時間分布情況
????distribution=linspace（00length（Pmatrix））;%用來記錄比例分布情況
????for?t=1:1:10^n
????????r1=unifrnd（01）;
????????r2=unifrnd（01）;
????????time=1/Qmatrix（ii）*log（r1）;%求出狀態停留時間下面求出下一次要跳到的狀態
????????for?k=1:1:length（Tmatrix）
????????????if?r2>Tmatrix（ik）
????????????????continue;
????????????elseif?r2<=Tmatrix（ik）
????????????????break;
????????????end
????????end
????????timedistribution（i）=timedistribution（i）+time;
????????totaltime=totaltime+time;
????????%記錄上一次執行的數據到矩陣state中去
????????i=k;
????end
????for?i=1:1:length（Pmatrix）
????????????distribution（i）=timedistribution（i）/totaltime;
????end
????format?long
????alltime=totaltime
????format?short
????%alltimedistribution=[alltimedistribution;totaltime]
????%已經求得一次distribution分布
????multidistribution=[multidistribution;distribution];
end
return;
%一組測試數據
%????0.2570????0.6436????0.0994
%????0.3589????0.3158????0.3253
%????0.4068????0.2576????0.3357
%????0.4326????0.2584????0.3090
%????0.4350????0.2618????0.3031
%????0.4357????0.2602????0.3041
%????0.4349????0.2607????0.3044

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2098??2014-12-01?12:10??隨機過程\test.m

?????文件???????1595??2014-12-01?12:10??隨機過程\test2.m

?????目錄??????????0??2015-01-19?15:43??隨機過程

-----------?---------??----------?-----??----

?????????????????3693????????????????????3