練習一

# 向量正規化計算，所有元素總和=1 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(.8, .2, .3, .7), nrow=2, ncol=2) A.Eigen <- Eigenize(A.Matrix) OutputEigen(A.Eigen)