MATLAB马尔科夫链预测法
发布时间：-06-10 09:10,
浏览次数：1195
, 标签：
MATLAB
1，确定初始概率分布p0
输入初始状态矩阵：
>>A=[11 19 36 37 33];
p0=A./sum(A)
2,计算一步转移概率矩阵
输入频数矩阵：
>>n1=[48 3 0 0 0 ];
>> n2=[2 20 6 0 0 ];
>> n3=[0 4 68 4 0 ];
>> n4=[0 0 1 21 0 ];
>> n5=[0 0 0 1 10];
>> p=[n1./sum(n1);n2./sum(n2);n3./sum(n3);n4./sum(n4);n5./sum(n5)]
3，预测
确定预测的时间(如下一个时间点、两个时间间隔等等)
输入：
>>p1=p0*p%或者p2=p0*p^2等等
4，讨论遍历性，求平稳分布
>>syms x1 x2 x3 x4 x5%定义符号变量，以便后续计算
eq1=x1-[x1,x2,x3,x4,x5]*p(:,1);%等式都为0
eq2=x2-[x1,x2,x3,x4,x5]*p(:,2);
eq3=x3-[x1,x2,x3,x4,x5]*p(:,3);
eq4=x4-[x1,x2,x3,x4,x5]*p(:,4);
eq5=x5-[x1,x2,x3,x4,x5]*p(:,5);
eq6=x1+x2+x3+x4+x5-1;
[x1 x2 x3 x4 x5]=solve(eq1,eq2,eq3,eq4,eq5,eq6)%解线性方程组，非线性用fsolve函数
完成。
