在 r 中导出具有正对角线的下三角矩阵
Deriving lower triangular matrix with positive diagonal in r
我正在创建优化算法,我需要插入 A 的初始值,它是一个具有正对角线的下三角矩阵。我的第一个问题是如何在r中导出一个具有正对角线的随机下三角矩阵作为初始值矩阵?选择随机矩阵是个好主意吗?如果不是,有什么方法可以更好地初步猜测这种类型的矩阵?
我可以想象根据您的具体应用会有更好的解决方案,但我们可以将对角线设置为半法线(即 |e|
其中 e
~ N(0,1))并将下三角非对角线元素设置为标准法线值。 ...
n <- 10
M <- diag(abs(rnorm(n)))
M[lower.tri(M, diag = FALSE)] <- rnorm(n*(n-1)/2)
我正在创建优化算法,我需要插入 A 的初始值,它是一个具有正对角线的下三角矩阵。我的第一个问题是如何在r中导出一个具有正对角线的随机下三角矩阵作为初始值矩阵?选择随机矩阵是个好主意吗?如果不是,有什么方法可以更好地初步猜测这种类型的矩阵?
我可以想象根据您的具体应用会有更好的解决方案,但我们可以将对角线设置为半法线(即 |e|
其中 e
~ N(0,1))并将下三角非对角线元素设置为标准法线值。 ...
n <- 10
M <- diag(abs(rnorm(n)))
M[lower.tri(M, diag = FALSE)] <- rnorm(n*(n-1)/2)