随机序列的Markov链预测

1978年至2004年城镇居民蔬菜人均消费数据，请给出城镇居民蔬菜消费量的长期趋势。

 

1.状态划分及构造

A=[130 137 143 152.34 159.08 165.48 149.04 144.36 148.32 142.58 147.02 144.56 138.70 132.18 124.91 ...
     120.64 120.74 116.47 118.51 113.34 113.76 114.94 114.74 115.86 116.52 118.34 122.30];
b=diff(A);
A1=A(1:end-1);
tt=b./A1*100;
z=[-3 0 3];
for i=1:length(tt)
    if tt(i)<z(1)
        aa(i)=1;
    elseif tt(i)<z(2)
        aa(i)=2;
    elseif tt(i)<z(3)
        aa(i)=3;
    else
        aa(i)=4;
    end
end

 

2.计算各状态转移次数及转移概率

for i=1:4
    for j=1:4
        k=[i,j];
        c(i,j)=length(strfind(aa,k));
    end
end
c
cs=sum(c,2);
p=zeros(4,4);
for i=1:4
    p(i,:)=c(i,:)./cs(i);
end
p

c =

     4     0     4     1
     1     0     1     0
     3     1     3     1
     1     1     0     4


p =

    0.4444         0    0.4444    0.1111
    0.5000         0    0.5000         0
    0.3750    0.1250    0.3750    0.1250
    0.1667    0.1667         0    0.6667

 

3.计算城镇居民年人均蔬菜消费量的长期趋势

p1=p'-eye(size(p'));
p2=[p1;ones(1,4)];
b11=[zeros(1,4),1]';
x=linsolve(p2,b11)

x =

    0.3600
    0.0800
    0.3200
    0.2400

环比增长率<-3%的概率:0.36

-3%<环比增长率<0%的概率:0.08

0%<环比增长率<3%的概率:0.32

环比增长率>3%的概率:0.24