R语言笔记--从矩阵的操作说起

# ----------- 矩阵工具-----------------
A <- array(1:6,dim=c(2,3))
B <- array(1:6,dim=c(3,2))
t(A) # 转置

# A
nrow(A)
ncol(A)

# 矩阵相乘 %*%
A * A # 对应元素相称
A %*% B # 矩阵相称
x <- array(1:2)
x %*% x # 二次型
y <- array(1:2)

# 下面的两种运算等价
corssprod(x,y)
t(x) %*% y # 更高效



# 函数diag() 的含义依赖于它的参数。当v 是一个向量时,diag(v)返回以该向
# 量元素为对角元素的对角矩阵。当M 是一个矩阵时,diag(M) 返回M的对角元素。
# 这和Matlab 中diag() 的用法完全一致。不过有点混乱的是,如果k 是单个值 4 ,那
# 么diag(k) 的结果就是k ×k 的方阵!