R Data Il 矩阵运算例子 A <-array(1:9,dim=c(3,3)) >A<-t(array(c(1:8,10),dim=c(3,3))) B <-array(1:9,dim=c(3,3)) >b<-c(1,1,1) >C<-A B x <solve(A,b) >C >X [,1][,2][,3] [1]-1.000000e+00 1.000000e+003.330669e-16 [1,] 1 16 49 [2,] 4 25 64 B <solve(A) [3,] 9 36 81 >B [,1] [,2][,3] D <-A%*%B [1,]-0.6666667-1.333333 1 D [2,]-0.66666673.666667 -2 [3,] 1.0000000-2.000000 1 [,1] [,2] [,3] [1,] 30 66102 >A [2,] 36 81126 > Sm <-crossprod(A,A) [,1][,2][,3] [3,] 42 96150 ev eigen(Sm) [1,] 1 <- 2,] 4 6 ev [3,] 810 M <-array(1:9,dim=c(3,3)) Svalues >M [1]303.19533618 0.76590739 0.03875643 [,1][,2][,3] [1,] 1 4 7 Svectors [2,] 2 5 8 [,1 [,2] [,3] [3,] 3 6 9 [1,]-0.4646675 0.833286355 0.2995295 diag(M) [2,]-0.5537546-0.009499485-0.8326258 [1]159 [3,]-0.6909703-0.5527599940.4658502R Data II 矩阵运算例子