markov_ergodic_test.m

来自「Markov分析的matlab工具包」· M 代码 · 共 45 行

%
%   随机数学(30250143, 07Fall) 大作业程序
%
%   作者： 王  音 2004011448 
%         葛元恺 2004011425 <[e-mail]>
%
%   All rights reserved. 2008.1

%一个测试遍历性的脚本

P = [   0   0   1   0; 
        0   0   1   0; 
        1/3 0   2/3 0; 
        1/2 1/2 0   0]   % P147的例子 
    
if Markov_Ergodic(P), % should be 1 (Yes)
    disp('该链是遍历链');
else
    disp('该链是非遍历链');
end


P2 = [  0   0   1   0; 
        0   0   1   0; 
        1/3 2/3 0   0; 
        1/2 1/2 0   0]  % P148的例子
    
if Markov_Ergodic(P2) % should be 0 (No)
    disp('该链是遍历链');
else
    disp('该链是非遍历链');
end


P3 = [  1/2 1/2 0   0;
        0   1/2 1/2 0;
        1/2 0   1/2 0;
        0   0   0   1]  %一个非遍历链的例子

if Markov_Ergodic(P3)  %should be 0 (No)
    disp('该链是遍历链');
else
    disp('该链是非遍历链');
end