R 矩阵的特征值，假设对称

我试图在Julia中找到以下2 X 2矩阵（等于

a

）的特征值：

2×2 Array{Float64,2}:
0.120066  0.956959
0.408367  0.422321

我在R中有相同的数组，运行以下R命令可以得到以下输出：

eigen(a, symmetric=T, only.values=T)
$values
[1]  0.706626 -0.164245

eigvals(LowerTriangular(a))
2-element Array{Float64,1}:
0.120066
0.422321

然而，在Julia中，当我运行此命令时，我得到以下输出：

eigen(a, symmetric=T, only.values=T)
$values
[1]  0.706626 -0.164245

eigvals(LowerTriangular(a))
2-element Array{Float64,1}:
0.120066
0.422321

有没有办法复制Julia中对称矩阵的R

eigen（）

函数，因为我使用的

LowerTriangular

函数不起作用

> x
         [,1]     [,2]
[1,] 0.120066 0.956959
[2,] 0.408367 0.422321

在Julia中，

eigvals（LowerTriangular（a））

计算

x

下三角部分的特征值（即，严格上三角部分的条目设置为0）：

---

使用对称函数如下：

julia> eigvals(Symmetric(x, :L))
2-element Array{Float64,1}:
 -0.164241
  0.706628

---

由于Julia 0.7，您必须使用

使用LinearAlgebra

来导入函数。

能否请您澄清您执行

特征（…，symmetric=TRUE）的意图

在一个非对称矩阵上？@ArtemSokolov我在处理一个对称矩阵问题/评论：你在一个显然是非对称的矩阵上计算对称特征值，这似乎很奇怪<代码>？eigen阐明了当使用

symmetric=TRUE

时，“仅使用[矩阵的]下三角形（包括对角线）”。因此，

Symmetric

将如@BogumiłKamiński下面解释的那样，通过简单地将下三角复制到上三角来实现对称。我如何才能使R输出的值与Julia代码输出的值相等？@L.Klum这在我的答案中显示。无论如何，三角形矩阵的特征值是其对角线上的条目。你不需要

eigen

…对不起，我不太清楚（我的英语不是很好）。我需要在Julia中创建R中的0.7066279和-0.1642409值。