求马尔可夫模型的转移矩阵
要在Matlab中解决:
clc,clear,format rat
a=[ 4 3 2 1 4 3 1 1 2 3
2 1 2 3 4 4 3 3 1 1
1 3 3 2 1 2 2 2 4 4
222312431 );
a=a’; a=a(3360 ); 将%矩阵a逐行展开为一个行向量
for i=1:4
for j=1:4
f(I,j )=Length(Findstr([I,j],a ); %统计子字符串’ij’的个数
结束
结束
Ni=sum(f,2 ); %计算矩阵f的行和
phat=f./repmat(Ni,1,size(f ) f,2 ) ) %求出状态迁移的频率
格式%将恢复为较短数量的显示格式
(() ) ) ) (() ) () ) ) ) ) ) ) ) ) ) ) ) ) ) ) 652
结果如下。
phat=
2/5 2/51/101/10
3/11 2/11 4/11 2/11
4/11 4/11 2/11 1/11
0 1/74/7 2/7
(1)或份额)我喜欢0