S = matrix(c(0,1/3,1/3,1/3,0,0,1/2,1/2,0,0,0,1,0,1,0,0),4)
U = matrix(rep(1,times=16),4)
a= 0.85
G = a * S + ((1-a) * 1/4) * U
q0 = rep(1,4)
q = q0
repeat {
q = G%*% q
d = dist(matrix(c(q,q0),2,byrow=TRUE), method = "euclidean")
if (d <1e-6 ) break
q0 = q
}
print(q)
Normal 0 7.8 磅 0 2 false false false EN-US ZH-CN X-NONE Normal 0 7.8 磅 0 2 false false false EN-US ZH-CN X-NONE Normal 0 7.8 磅 0 2 false false false EN-US ZH-CN X-NONE
来自 “ ITPUB博客 ” ,链接:https://blog.itpub.net/26857237/viewspace-745536/ ,如需转载,请注明出处,否则将追究法律责任。