请问R中矩阵幂次怎么算？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问R中矩阵幂次怎么算？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：矩阵

A%*%A

自己写个函数用loop实现n次%*%。
或者安装个expm的package然后就可以调用%^%

ywh19860616 发表于 2012-5-19 08:43
A%*%A

那如果是-0.5次呢？？

zgu9 发表于 2012-5-19 08:52
自己写个函数用loop实现n次%*%。
或者安装个expm的package然后就可以调用%^%

感谢先> <..那如果是-0.5次呢？？

e <- eigen(A)
V <- e$vectors
g <- e$values
B <- V %*% diag(1/sqrt(g)) %*% t(V)
B %*% t(B) - solve(A)

每天都在放空 发表于 2012-5-19 09:01
感谢先>

X <- matrix(c(2,1,1,2),nrow=2)
  X
#      [,1] [,2]
# [1,]    2    1
# [2,]    1    2

  E <- eigen(X)
  V <- E$values
  Q <- E$vectors
  Y <- Q%*%diag(1/sqrt(V))%*%t(Q)
  Y
#            [,1]       [,2]
# [1,]  0.7886751 -0.2113249
# [2,] -0.2113249  0.7886751

  solve(Y%*%Y)    ## i.e. find its inverse
#      [,1] [,2]
# [1,]    2    1
# [2,]    1    2

# Hence (Y%*%Y)^(-1) = X, or Y = X^(-1/2)

zgu9 发表于 2012-5-19 10:40
X

你这个只能对对称矩阵做吧。。