R语言应用于分析矩阵

R语言应用于分析矩阵

学习《统计建模与R软件》（薛毅2006）这本书，使用R语言对书中练习题解答。矩阵分析为书中第二章内容，练习题如下：

2.1 建立一个R文件，在文件中输入变量 \(x = (1, 2, 3)^T\)，\(y = (4, 5, 6)^T\)，并作以下运算。
（1）计算\(z = 2x + y + e\)，其中\(e = (1, 1, 1)^T\)；
（2）计算\(x\)和\(y\)的内积；
（3）计算\(x\)和\(y\)的外积。
解：

x <- c(1, 2, 3)
y <- c(4, 5, 6)
e <- c(1, 1, 1)
#建立向量x，y和e

z = 2*x + y + e
print(z) #第一问矩阵加法计算 

## [1]  7 10 13

x %*% y #第二问矩阵内积计算，或者使用crossprod(x, y) 

##      [,1]
## [1,]   32

x %o% y #第三问矩阵外积计算，或者使用tcrossprod(x, y)

##      [,1] [,2] [,3]
## [1,]    4    5    6
## [2,]    8   10   12
## [3,]   12   15   18

---

2.2 将1，2，…，20构成两个\(4 × 5\)阶的矩阵，其中矩阵A是按列输入，矩阵B是按行输入，并作如下运算。
（1）\(C = A + B\)；
（2）\(D = AB\)