12/01/17 22:19:59.38
>>307
R-Sources URLﾘﾝｸ(cse.naro.affrc.go.jp)
見たけど、これで行けるんじゃね？
> A <- matrix(c(1,2,0,3), 2, 2)
> B <- eigen(A)$vectors %*% diag(1/sqrt(eigen(A)$values)) %*% solve(eigen(A)$vectors)
> solve(A)
[,1] [,2]
[1,] 1.0000000 0.0000000
[2,] -0.6666667 0.3333333
> B %*% B
[,1] [,2]
[1,] 1.0000000 0.0000000
[2,] -0.6666667 0.3333333
> A %*% B %*% B
[,1] [,2]
[1,] 1.000000e+00 0
[2,] 2.217736e-16 1
最後のがキッチリ単位行列にならないのは、
丸め誤差の所為だから仕方ないと思われ。