如何将 R 中的矩阵化为阶梯形？

矩阵的阶梯形是指具有以下特征的矩阵

1. 每一行中的第一个非零元素，称为前导项，为 1。

2. 每个前导项所在列都在上一行的前导项之后。

3. 全部元素都为零的行（如果有）位于带有非零元素的行下面。

在 R 中，可以使用 matlib 包的梯度函数来查找矩阵的梯度形式。

示例

在线演示

> M<-matrix(rpois(25,10),ncol=5)
> M

输出

   [,1] [,2] [,3] [,4] [,5]
[1,] 8 11 3 10 13
[2,] 9 9 7 15 11
[3,] 10 13 10 14 13
[4,] 7 11 11 12 17
[5,] 13 10 9 20 13

示例

在线演示

> V<-rpois(5,2)
> V

输出

[1] 4 2 4 1 2

加载 matlib 包并找出矩阵 M 的梯度形式

示例

> library(matlib)
> echelon(M,V,verbose=TRUE,fractions=TRUE)

初始矩阵

输出

[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 8 11 3 10 13 4
[2,] 9 9 7 15 11 2
[3,] 10 13 10 14 13 4
[4,] 7 11 11 12 17 1
[5,] 13 10 9 20 13 2

行：1

输出

exchange rows 1 and 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 13 10 9 20 13 2
[2,] 9 9 7 15 11 2
[3,] 10 13 10 14 13 4
[4,] 7 11 11 12 17 1
[5,] 8 11 3 10 13 4

multiply row 1 by 1/13
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 9 9 7 15 11 2
[3,] 10 13 10 14 13 4
[4,] 7 11 11 12 17 1
[5,] 8 11 3 10 13 4

multiply row 1 by 9 and subtract from row 2
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 27/13 10/13 15/13 2 8/13
[3,] 10 13 10 14 13 4
[4,] 7 11 11 12 17 1
[5,] 8 11 3 10 13 4

multiply row 1 by 10 and subtract from row 3
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 27/13 10/13 15/13 2 8/13
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 7 11 11 12 17 1
[5,] 8 11 3 10 13 4

multiply row 1 by 7 and subtract from row 4
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 27/13 10/13 15/13 2 8/13
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 0 73/13 80/13 16/13 10 -1/13
[5,] 8 11 3 10 13 4

multiply row 1 by 8 and subtract from row 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 27/13 10/13 15/13 2 8/13
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 0 73/13 80/13 16/13 10 -1/13
[5,] 0 63/13 -33/13 -30/13 5 36/13

行：2

输出

exchange rows 2 and 4
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 73/13 80/13 16/13 10 -1/13
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 0 27/13 10/13 15/13 2 8/13
[5,] 0 63/13 -33/13 -30/13 5 36/13

multiply row 2 by 13/73
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 10/13 9/13 20/13 1 2/13
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 0 27/13 10/13 15/13 2 8/13
[5,] 0 63/13 -33/13 -30/13 5 36/13

multiply row 2 by 10/13 and subtract from row 1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 69/13 40/13 -18/13 3 32/13
[4,] 0 27/13 10/13 15/13 2 8/13
[5,] 0 63/13 -33/13 -30/13 5 36/13

multiply row 2 by 69/13 and subtract from row 3
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 -200/73 -186/73 -471/73 185/73
[4,] 0 27/13 10/13 15/13 2 8/13
[5,] 0 63/13 -33/13 -30/13 5 36/13

multiply row 2 by 27/13 and subtract from row 4
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 -200/73 -186/73 -471/73 185/73
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 63/13 -33/13 -30/13 5 36/13

multiply row 2 by 63/13 and subtract from row 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 -200/73 -186/73 -471/73 185/73
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 0 -573/73 -246/73 -265/73 207/73

行：3

输出

exchange rows 3 and 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 -573/73 -246/73 -265/73 207/73
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 0 -200/73 -186/73 -471/73 185/73

multiply row 3 by -73/573
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 -11/73 100/73 -27/73 12/73
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 0 -200/73 -186/73 -471/73 185/73

multiply row 3 by 11/73 and add to row 1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 80/73 16/73 130/73 -1/73
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 0 -200/73 -186/73 -471/73 185/73

multiply row 3 by 80/73 and subtract from row 2
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 -110/73 51/73 -124/73 47/73
[5,] 0 0 -200/73 -186/73 -471/73 185/73

multiply row 3 by 110/73 and add to row 4
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 257/191 -574/573 19/191
[5,] 0 0 -200/73 -186/73 -471/73 185/73

multiply row 3 by 200/73 and add to row 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 257/191 -574/573 19/191
[5,] 0 0 0 -262/191 -2971/573 295/191

行：4

输出

exchange rows 4 and 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 -262/191 -2971/573 295/191
[5,] 0 0 0 257/191 -574/573 19/191

multiply row 4 by -191/262
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 274/191 -172/573 21/191
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 257/191 -574/573 19/191

multiply row 4 by 274/191 and subtract from row 1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 -2249/393 226/131
[2,] 0 1 0 -48/191 730/573 73/191
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 257/191 -574/573 19/191

multiply row 4 by 48/191 and add to row 2
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 -2249/393 226/131
[2,] 0 1 0 0 874/393 13/131
[3,] 0 0 1 82/191 265/573 -69/191
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 257/191 -574/573 19/191

multiply row 4 by 82/191 and subtract from row 3
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 -2249/393 226/131
[2,] 0 1 0 0 874/393 13/131
[3,] 0 0 1 0 -152/131 16/131
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 257/191 -574/573 19/191

multiply row 4 by 257/191 and subtract from row 5
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 -2249/393 226/131
[2,] 0 1 0 0 874/393 13/131
[3,] 0 0 1 0 -152/131 16/131
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 0 -1595/262 423/262

行：5

输出

multiply row 5 by -262/1595
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 -2249/393 226/131
[2,] 0 1 0 0 874/393 13/131
[3,] 0 0 1 0 -152/131 16/131
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 0 1 -266/1003

multiply row 5 by 2249/393 and add to row 1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 0 331/1595
[2,] 0 1 0 0 874/393 13/131
[3,] 0 0 1 0 -152/131 16/131
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 0 1 -266/1003

multiply row 5 by 874/393 and subtract from row 2
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 0 331/1595
[2,] 0 1 0 0 0 1099/1595
[3,] 0 0 1 0 -152/131 16/131
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 0 1 -266/1003

multiply row 5 by 152/131 and add to row 3
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 0 331/1595
[2,] 0 1 0 0 0 1099/1595
[3,] 0 0 1 0 0 -296/1595
[4,] 0 0 0 1 2971/786 -295/262
[5,] 0 0 0 0 1 -266/1003

multiply row 5 by 2971/786 and subtract from row 4
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 0 331/1595
[2,] 0 1 0 0 0 1099/1595
[3,] 0 0 1 0 0 -296/1595
[4,] 0 0 0 1 0 -197/1595
[5,] 0 0 0 0 1 -266/1003

Nizamuddin Siddiqui

更新于：2020 年 11 月 23 日

1K+ 查看次数

启动你的事业

通过完成课程获得认证

开始