練習四

使用 GNU-R 求解：

Normalize <- function(V) { Base.Value <- 0 nrow <- length(V) for (i in seq(1:nrow)) { Base.Value <- Base.Value + abs(V[i]) } R <- V / Base.Value R } Eigenize <- function(A) { myEigen <- eigen(A) nrow <- length(myEigen$values) for (i in seq(1:nrow)) { myEigen$vector[, i] <- Normalize(myEigen$vector[, i]) } myEigen } OutputEigen <- function(myEigen) { print(myEigen$values) nrow <- length(myEigen$values) for (i in seq(1:nrow)) { print(myEigen$vector[, i]) } } A.Matrix <- matrix(c(5, -1, 3, -1, 2, -2, 3, -2, 3), nrow=3, ncol=3) A.Eigen <- Eigenize(A.Matrix) OutputEigen(A.Eigen)