r - 导出r中具有正对角线的下三角矩阵

Question

我正在创建优化算法，我需要插入 A 的初始值，它是具有正对角线的下三角矩阵。我的第一个问题是如何导出一个随机下三角矩阵，正对角线作为 r 中的初始值矩阵？选择随机矩阵是个好主意吗？如果不是，有什么方法可以对这种类型的矩阵进行更好的初始猜测？

标签： roptimizationlinear-algebranumerical-analysis

Answer 1

我可以想象根据您的特定应用程序有更好的解决方案，但我们可以将对角线设置为半法线（即|e|~ eN(0,1)）并将下三角非对角线元素设置为标准法线值。...

n <- 10
M <- diag(abs(rnorm(n)))
M[lower.tri(M, diag = FALSE)] <- rnorm(n*(n-1)/2)