马尔科夫链转移矩阵

clear all;
%  [1] 华林香等. 马尔可夫模型在一次能源消费预测中的应用——以福建省为例[J]. 福建师范大学学报:自然科学版 2013 29（5）:78-86.
%  [2] 王锋. 中国碳排放增长的驱动因素及减排政策评价[M]. 经济科学出版社 2011.
%  论文[1]中的数据
y=[63.8	25.1	0.2	10.9;
    62	22.3	0.1	15.6;
    61.8	21.6	0.1	16.5;
    65.2	21.4	0.1	13.3;
    64.8	18.9	0.3	16;
    67.6	18.3	1.3	12.8;
    57.8	23.6	4	14.6;
    64.4	22.5	4.3	8.8;];
%  论文[2]中的数据
% y= [67.75 23.21 2.35 6.69;
%     66.68 22.87 2.55 7.9;
%     66.32 23.41 2.56 7.71;
%     68.38 22.21 2.58 6.83;
%     67.99 22.33 2.60 7.08;
%     69.10 21.00 2.80 7.10;
%     69.40 20.40 3.03 7.20;
%     69.50 19.70 3.50 7.30;
%     68.67 18.68 3.77 8.89;
%     ] ;

for m = 1:（size（y1）-1）
    y1 = y（m :）;
    y2 = y（m+1 :）;
    y3 = y2 - y1;             %差值矩阵
    FF（::m）=zeros（size（y12）size（y12））; %每步概率转移矩阵
    P1=0;    
    for i = 1:size（y12）
        for j = 1:size（y12）
            if y3（i） < 0             % 减少
                FF（:im） = 0;
                FF（iim） = y2（i）/y1（i）;
                if y3（j） > 0
                    FF（ijm） = （（1- FF（iim））*y3（j））/（sum（y3） - y3（i））;
                    P1 = P1 + y3（j）;
                end
            else                     %增加
                FF（i:m） = 0;
                FF（iim） = 1;
            end
        end
        % 修正转移概率矩阵
        for n = 1:size（y12） 
            if y3（n） > 0 && y3（n） == P1 && sum（FF（i:m）） > 1
                FF（inm） = 1 - FF（iim）;
            end
            if y3（n） > 0 && y3（n） < P1
                FF（inm） = （1 - FF（iim））*（y3（n）/P1）;
            end
            
        end
        P1=0;
    end    
end
% 计算平均转移矩阵P1
P = FF（::1）;
for n = 2:（size（y1）-1） 
    P = P * FF（::n）;
end
P1 = P^（1/（size（y1）-1））;  

%  A = B .^（1/n）; % .^表示对矩阵B中每个元素进行操作，所以B.^（1/n）即表示对矩阵B中每个元素开n次方
%  a = b ^（1/n）; % 对变量b开n次方
   


        
        
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2043  2016-06-22 07:35  MarkovUP.m

     文件    1418392  2016-06-19 19:08  [1]中国碳排放增长的驱动因素及减排政策评价.pdf

     文件     345999  2016-06-20 03:29  [2]马尔可夫模型在一次能源消费预测中的应用_以福建省为例.pdf

----------- ---------  ---------- -----  ----

              1766434                    3