行列の指数関数というのがある 簡単に言うと、固有値を指数の肩に乗せて、それを固有ベクトルの行列でサンドイッチする exp.m <- function(A,n=1,eigen=FALSE){ # 固有値分解 eigen.out<-eigen(A) # P=V,P^{-1}=U V<-eigen.out[[2]] U<-solve(V) B<-diag(exp…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。