流れる空の中で数学を。

とある数学好きの「手作りすうがく」と「気ままな雑記」。

線形代数

【Python】5000×5000直交行列の対角化【線形代数】

巨大なランダム直交行列 各成分がの一様乱数で決定される5000×5000行列を近似的に対角化してみた。 手法 規格化したベクトルを用意し、各成分がという位相だけを持つとして、をモンテカルロ法により最適化し、対角化を行った。 プログラム github.com gistc1…