假设我有以下 3 x 3 矩阵：

p<-3
X<-matrix(1/p,p,p)

--只是一个 x矩阵，其中每个条目都是。现在我想解决系统$\pmb X$$p$$p$$1/p$

其中是长度为的 1 向量。我知道是这样一个系统的解决方案，但是当我这样做时，$\pmb 1_p$$p$$\pmb 1_p$R

solve(X,rep(1,p))

我明白了

Error in solve.default(X, rep(1, p)) : 
  Lapack routine dgesv: system is exactly singular: U[2,2] = 0

而不是。我的问题是为什么$\pmb 1_p$