matsueushi | Juliaの単位行列

何かと必要になる単位行列だが、Julia v1.0以降だと LinearAlgebra をインポートして Matrix と I の組み合わせで作成するのが楽なようだ。

using LinearAlgebra
using Distributions
Matrix(I, 2, 2)

2×2 Array{Bool,2}:
  true  false
 false   true

Matrix{Int64}(I, 2, 2)

2×2 Array{Int64,2}:
 1  0
 0  1

Matrix{Float64}(I, 2, 2)

2×2 Array{Float64,2}:
 1.0  0.0
 0.0  1.0

MvNormal の分散を単位行列にしたいときは、 I だけで大丈夫。

MvNormal([3, 4], I)

IsoNormal(
dim: 2
μ: [3.0, 4.0]
Σ: [1.0 0.0; 0.0 1.0]
)

参考