代码语言:english复制
#矩阵
A<-matrix(c(7,4,-3,4,6,93,0,93,23,5),byrow=TRUE,nrow=5,ncol=2)
dim(A)
B<-diag(c(1,2,3,4,5)) #对角矩阵
I<-diag(c(1),nrow=5,ncol=5) #单位矩阵
矩阵运算
代码语言:english复制C<-A B #相同位置元素相加或相减
C<-A%*%B#矩阵相乘
C<-A%x%B #直积
tA<-t(A) #转置
trA<-sum(diag(A)) #迹和
r<-qr(D)$rank #秩(rank)
detA<-det(A) # 行列式
Ainv<-solve(A)
#奇异阵的广义逆矩阵
library(MASS)
ginv(D)
#正定矩阵的cholesky分解
T<-chol(A)
#特征值与特征向量
eigen(A)