diag()提取或替换一个矩阵的对角线,或构造一个对角矩阵。 语法： diag(x = 1, nrow, ncol) diag(x) x x [,1] [,2] [,3] [,4] [1,] 1 5 9 13 [2,] 2 6 10 14 [3,] 3 7 11 15 [4,] 4 8 12 16 > diag(x) [1] 1 6 11 16 #展示x的对角线数据 > diag(diag(x)) [,1] [,2] [,3] [,4] [1,] 1 0 0 0 [2,] 0 6 0 0 [3,] 0 0 11 0 [4,] 0 0 0 16 #由x的对角线数据生成对角矩阵 > diag(4) [,1] [,2] [,3] [,4] [1,] 1 0 0 0 [2,] 0 1 0 0 [3,] 0 0 1 0 [4,] 0 0 0 1 #对角线为4个元素每个元素默认为1的对角矩阵 > diag(10,3,4) [,1] [,2] [,3] [,4] [1,] 10 0 0 0 [2,] 0 10 0 0 [3,] 0 0 10 0 #对角线数据位10,三行四列的矩阵 > a a [,1] [,2] [,3] [,4] [1,] 2 1 1 1 [2,] 1 2 1 1 [3,] 1 1 2 1 [4,] 1 1 1 2

看几个较复杂的例子 例1

3、判断等式两边的矩阵数值是否相等

4、所以这条语句的意思是取步骤3中矩阵的对角线

1、rep(“”,3) 把”“重复3次

2、colnames(M) 显示M的列名

3、把M的行命名