一個(gè)時(shí)間序列如下：
0.19
0.19
0.13
0.01
-0.03
-0.03
0.03
0.09
-0.13
-0.13
0.05
-0.03
0.03
0.09
-0.07
0.11
0.05
0.01
-0.05
-0.01
0.07
0.01
-0.15
-0.01
······這個(gè)時(shí)間序列一共有4032個(gè)值，這些值的取值空間有115個(gè)取法a=[-0.49，-0.47,-0.45······1.75,1.77,1.79]，現(xiàn)在想求出來這個(gè)115*115的狀態(tài)轉(zhuǎn)移概率矩陣，寫了一個(gè)matlab程序但是計(jì)算出來有問題，想問問大神們?nèi)绻肞ython做應(yīng)該怎么做，作為剛剛?cè)腴T的小白，只知道用if語句來實(shí)現(xiàn)，但是這樣的話，就要寫115個(gè)if語句了，怎么做比較簡(jiǎn)潔呢，謝謝大家
1
0
0
不太清楚馬爾可夫鏈如何計(jì)算，猜測(cè)是狀態(tài)轉(zhuǎn)移，看看下面的有沒有用。
a = [1,2,3,0,0,0,3,3,2,3,3,2,3,2,1,2,3]
l = 4
N = [[0] * l for i in range(l)]
amount = 0
for i in data:
print i
amount += 1
N[i[0]][i[1]] += 1
print amount
for i in range(l):
for j in range(l):
N[i][j] /= float(amount)
print N
