我有列数等于行数。对角线等于零。我怎样才能建立这个矩阵?茱莉亚的下三角矩阵
#mat
# [,1] [,2] [,3] [,4]
#[1,] 0 NA NA NA
#[2,] 1 0 NA NA
#[3,] 2 4 0 NA
#[4,] 3 5 6 0
我想这
x=rand(4,4)
4x4 Array{Float64,2}:
0.60064 0.917443 0.561744 0.135717
0.106728 0.72391 0.0894174 0.0656103
0.410262 0.953857 0.844697 0.0375045
0.476771 0.778106 0.469514 0.398846
c=LowerTriangular(x)
4x4 LowerTriangular{Float64,Array{Float64,2}}:
0.60064 0.0 0.0 0.0
0.106728 0.72391 0.0 0.0
0.410262 0.953857 0.844697 0.0
0.476771 0.778106 0.469514 0.398846
但我米寻找这样的事情
c=LowerTriangular(x)
4x4 LowerTriangular{Float64,Array{Float64,2}}:
0.0 NA NA NA
0.106728 0.0 NA NA
0.410262 0.953857 0.0 NA
0.476771 0.778106 0.469514 0
对角线应该等于零。
只是FYI茱莉亚不使用NA值。我们有NaN –
x-ref:https://groups.google.com/forum/#!topic/julia-users/VGbdlfbRfbc – StefanKarpinski
此外,请考虑仍使用零而不是NA(或NaN)。为什么? 1.线性代数的意义大部分是使用NA丢失的,这将阻止任何现有的线性代数代码的使用。已知新生感染有传染性和产生病毒。如果你更多地描述用例,这个选择可能会更清楚。 –