大数据矩阵计算基础(二):R、Excel、Matlab中矩阵操作

一、R语言

• 构造矩阵

  a=matrix(c(1,2,3,4),ncol=2,byrow=T)
  b=matrix(c(5,6,7,8),ncol=2,byrow=T)
  c=matrix(c(1,2,3,4),ncol=2,byrow=F) 
  function (data = NA, nrow = 1, ncol = 1, byrow = FALSE, dimnames = NULL) 
  data项为必要的矩阵元素，nrow为行数，ncol为列数，注意nrow与ncol的乘积应为矩阵元素个数，byrow项控制排列元素时是否按行进行，dimnames给定行和列的名称。

 

 

一、R语言 构造矩阵 a=matrix(c(1,2,3,4),ncol=2,byrow=T) b=matrix(c(5,6,7,8),ncol=2,byrow=T) c=matrix(c(1,2,3,4),ncol=2,byrow=F) function (data = NA, nrow = 1, ncol = 1, byrow = FALSE, dimnames = NULL)  data项为必要的矩阵元素，nrow为行数，ncol为列数，注意nrow与ncol的乘积应为矩阵元素个数，byrow项控制排列元素时是否按行进行，dimnames给定行和列的名称。 > rowname [1] "r1" "r2" "r3" > colname=c("c1","c2","c3","c4") > colname [1] "c1" "c2" "c3" "c4" > matrix(1:12,nrow=3,ncol=4,dimnames=list(rowname,colname))   c1 c2 c3 c4 r1 1 4 7 10 r2 2 5 8 11 矩阵线性运算 a+b：对应元素相加 a-b：对应元素相减 2*a：数与每个元素相乘 a*b： 简单的对应元素相乘，注意不是矩阵相乘 矩阵乘法 a%*%b：矩阵乘法（不懂矩阵乘法的百度） 若A为n×m矩阵，要得到A'B，可用函数crossprod()，该函数计算结果与t(A)%*%B相同，但是效率更高。 > A=matrix(1:12,nrow=4,ncol=3) > B=matrix(1:12,nrow=4,ncol=3) > t(A)%*%B     [,1] [,2] [,3] [1,]  30   70 110 [2,]  70  174 278 [3,] 110  278 446 > crossprod(A,B)     [,1] [,2] [,3] [1,]  30  70 110 [2,]  70 174 278 [3,] 110 278 446 矩阵转置 t(a)、t(b) 取方阵的对角线元素 diag(a) 求方阵的迹 在线性代数中，一个n×n矩阵A的主对角线（从左上方至右下方的对角线）上各个元素的总和被称为矩阵A的迹（或迹数），一般记作tr(A)。 命令：sum(diag(a)) 构造对角矩阵 c=diag(c(1,2,3,4)) > diag(diag(A))     [,1] [,2] [,3] [,4] [1,]   1   0   0   0 [2,]   0   6   0   0 [3,]   0   0   11   0 [4,]   0   0   0   16 矩阵求逆 solve(a) 矩阵求逆可用函数solve()，应用solve(a, b)运算结果是解线性方程组ax = b，若b缺省，则系统默认为单位矩阵，因此可用其进行矩阵求逆， > a=matrix(rnorm(16),4,4) > a             [,1]     [,2]     [,3]     [,4] [1,] 1.6986019   0.5239738 0.2332094 0.3174184 [2,] -0.2010667 1.0913013 -1.2093734   0.8096514 [3,] -0.1797628 -0.7573283 0.2864535 1.3679963 [4,] -0.2217916 -0.3754700 0.1696771 -1.2424030 > solve(a)               [,1]     [,2]     [,3]     [,4] [1,] 0.9096360 0.54057479 0.7234861 1.3813059 [2,] -0.6464172 -0.91849017 -1.7546836 -2.6957775 [3,] -0.7841661 -1.78780083 -1.5795262 -3.1046207 [4,] -0.0741260 -0.06308603 0.1854137 -0.6607851 > solve (a) %*%a                 [,1]       [,2]           [,3]       [,4] [1,] 1.000000e+00 2.748453e-17 -2.787755e-17 -8.023096e-17 [2,] 1.626303e-19 1.000000e+00 -4.960225e-18 6.977925e-16 [3,] 2.135878e-17 -4.629543e-17 1.000000e+00 6.201636e-17 [4,] 1.866183e-17 1.563962e-17 1.183813e-17 1.000000e+00 矩阵的行列式 det(a) det(b) det(c) 二、Excel(复杂不常用) 定义矩阵名称 矩阵加法 矩阵减法 矩阵数乘 矩阵转置 transpose Shift+Ctrl+Enter 矩阵对应元素相乘 矩阵乘法 mmult 矩阵的逆 minverse 矩阵的行列式 mdeterm 计算矩阵对应元素乘积之和 sumproduct Shift+Ctrl+Enter 由于Excel矩阵操作比较复杂，描述起来比较麻烦，大家可以百度一下 三、Matlab MATLAB（矩阵实验室）是MATrix LABoratory的缩写，是一款由美国The MathWorks公司出品的