用lingo求状态转移矩阵的稳态概率分布

用lingo解决的代码，如下：

model:
 !求转移矩阵的稳态概率分布;
 sets:
  prod /1..3/:P;
  link(prod,prod):T;
 endsets
 
data: !转移矩阵;
 T= 0.368 0     0.632
    0.368 0.368 0.264
    0.184 0.368 0.448;
 enddata
 @for(prod(i):
      !去掉一个冗余约束;    
      P(i)=@sum(link(j,i):P(j)*T(j,i)));
 @sum(prod: P)=1;
 @for(prod(i):
   @warn('输入矩阵的每行之和必须是1',
      @abs(1-@sum(link(i,j): T(i,j)))
       #gt# 0.000001););
 end

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

结果为：

解释如下：

    即矩阵P的稳态概率分布为：w=[0.285, 0.263, 0.452]（保留三位有效小数）。