matlabではexpm(A)で指数行列の計算ができる
Rではexp(A)では行列Aのそれぞれの要素のExponentialを取ってしまう
そこでパッケージMatrixをインストール
install.packages("Matrix")
expm(A)でmatlabと同じように計算できる
要素が0となる場合は「・」のような表示になる
そもそも自分で関数を定義すればよかった
matexp <- function(x) {
y <- eigen(x)
y$vectors %*% diag( exp(y$values) ) %*% t(y$vectors)
}
getwd() #カレントディレクトリ
setwd("~/") #ディレクトリ変更
factor()