如何编写n圆的过渡矩阵

Question

我已经编写了创建空白n×n矩阵的代码，但是我不确定如何使用for循环或lapply函数来填充矩阵。

在n圈中，第i个状态的转移概率应为状态i-1和i + 1的一半，而其他状态为0。

C <- matrix(rep(0,n^2),ncol = n, nrow = n)

Answer 1

Yo可以创建矩阵，然后填充对角线：

n<-10
m<-matrix(0,n,n)
diag(m[-1,]) <-0.5 #fill the second diagonal
diag(m[,-1]) <-0.5 #fill the other second diagonal
m[1,n]<-0.5
m[n,1]<-0.5