初等变换的一个应用矩阵的满秩分解

第25卷第5期2009年10月

大 学 数 学

COLL EGE MA T H EMA TICS

Vol. 25, №. 5

Oct. 2009

初等变换的一个应用:矩阵的满秩分解

靳全勤

(同济大学数学系, 上海200092)

  [摘 要]给出利用初等变换求矩阵满秩分解的一个简洁方法.

[关键词]矩阵; 初等变换; 满秩分解

[中图分类号]O151. 2  [文献标识码]C   [文章编号]167221454(2009) 0520195203

矩阵是线性代数中最为重要的核心内容, 很多问题(如线性方程组求解、换的最简矩阵表示等等) . 的最有效工具之一. 熟知, 通过矩阵的初等变换, , 化二次型为标准型, 判别方阵是否可逆以及求解逆矩阵法. 内容虽然浅显简单, . , 敬请同仁批评指正.

设A 是一个, B 和r ×n 矩阵C 满足A =BC, 则称A =BC 为矩阵A 矩阵r m 和列数n. 如果n 阶方阵A 的秩r 远比其阶数n 小, 则满秩分. 设n 阶方阵A 有满秩分解A =BC, 根据文献[1], 我们有

n -r

|λE n -A |=|λE n -BC |==λ|λE r -CB |. 计算左边的行列式简单的多!

设A 为m ×n 矩阵, 通过行的初等变换, 可以将矩阵A 化为如下形式的简化行阶梯形矩阵

0000000

(1)

如果n 阶方阵A 的秩r 远远小于其阶数n , 无疑通过上式右边的行列式来求方阵A 的特征值要比直接

…0…0…0…0…0…0…0

1000000

c 1i 1+1…c 1i 2-1………………

000000

0100000

c 1i 2+1c 2i 2+1

…c 1i r -1…c 2i r -1……………

00000

0000100

c 1i r +1c 2i r +1c 3i r +1c r -1i r +1c r i r +1

………

c 1n c 2n c 3n

000000

00000

…c r -1n ………

c r n

00

00

其中1≤i 1≤i 2≤…≤i r ≤max (m , n ) .

根据列向量组的性质, 行初等变换不改变矩阵的列向量组的线性关系. 从矩阵A 的简化行阶梯形矩阵可以读出下面结论:

(i ) 矩阵A 的秩为r ;

 [收稿日期]2007204209

 [基金项目]国家自然科学基金(10671142) ; 教育部高校博士点基金([1**********])

196大 学 数 学              第25卷

(ii ) 矩阵A 的第i 1, i 2, …, i r 个列向量组成的向量组{A i 1, A i 2, …, A i r }为其列向量组的极大线性无

关组;

(iii ) 矩阵A 的每一个列向量均可由(ii ) 中的极大线性无关组线性表示, 并且第j 个列向量的线性表达式为

A j =c 1j A i 1+c 2j A i 2+…+c rj A i r .

从而, 矩阵A 可以表示为

c 11

A =(A 1 A 2 … A n ) =(A i 1 A i 2 … A i r )

c 21c r 1

c 12c 22c r 2

…c

1n …c 2n …c rn

.

若记B =(A i 1 A i 2 … A i r ) , 则矩阵B 为m ×r 列满秩矩阵; 若记

c 11c 12…c 1n

C =

c 21c r 1

c 22c r 2

…c 2n …c rn

,

则C 为r ×n 矩阵. 因为矩阵C 的第i 1列, 第i 2列, …, 第i r 列构成r 阶单位矩阵, C 为行满秩

矩阵. 这样, A =BC 就是矩阵A 的一个满秩分解.

注 上面讨论说明, 任意矩阵都存在满秩分解. 例1 经过行的初等变换, 我们有

231121310-11-1

=

31

3

2

3842131

732500

1-1110

01

1-11-11

-1-1111

-2-22-12

000

100

100

100

200

,

于是矩阵A 的秩为2. 令B =, C =, 则A =BC 为矩阵A 的一个满秩分解.

该例中矩阵A 的列向量极大线性无关组不唯一, 导致其满秩分解也不唯一. 取不同的列向量极大

线性无关组, 则得不同的满秩分解. 例如, 取第1列和第3列为极大列无关组, 则

[**************]

A =

131112253

12210

38411

3701

000

-1-111-1-11

-1-11

-2-22

000

100

100

100

20.

2

此时, 令B 1=

131

, C 1=

21

, 则A =B 1C 1为矩阵A 的又一个满秩分解. 再如, 取第1列

和第4列为极大列无关组, 则23151213A =

35281324

43712-1-11

1-1-11

3-1-11

3-2-22

1-1100

-2100

0100

-3200

.

000

000

第5期          靳全勤:初等变换的一个应用:矩阵的满秩分解

2

5381

111

, C 2=

197

此时, 令B 2=

131

10

-11

-21

0111

-32

, 则A =B 2C 2为矩阵A 的另一个满秩分解.

例2 求n 阶矩阵A =

…ω…

……11的特征值.

1

解 

容易看到, 矩阵A 的秩为1, 其满秩分解为A =

1

(1 1 … 1) , 由等式(1) 得|λE n -A |

n -1

(λ-n ) , 所以矩阵A 的特征值为λλ=λ1=λ2=…=λn -1=0, n =n.

本文只讨论了行初等变换的情形. 对于列初等变换, 有类似结果.

[参 考 文 献]

[1] 同济大学应用数学系. 高等代数与解析几何[M ].北京:of Two Full R ank Matrices :

of Elementary Matrix Operations

J I N Quan 2qi n

(Department of Mathematics , Tongji University , Shanghai 200092, China )

Abstract :We provide a simple method for decomposing a matrix into a product of two f ull rank matrices using elementary matrix operations.

K ey w ords :matrix ; elementary operation ; f ull rank decomposition

第25卷第5期2009年10月

大 学 数 学

COLL EGE MA T H EMA TICS

Vol. 25, №. 5

Oct. 2009

初等变换的一个应用:矩阵的满秩分解

靳全勤

(同济大学数学系, 上海200092)

  [摘 要]给出利用初等变换求矩阵满秩分解的一个简洁方法.

[关键词]矩阵; 初等变换; 满秩分解

[中图分类号]O151. 2  [文献标识码]C   [文章编号]167221454(2009) 0520195203

矩阵是线性代数中最为重要的核心内容, 很多问题(如线性方程组求解、换的最简矩阵表示等等) . 的最有效工具之一. 熟知, 通过矩阵的初等变换, , 化二次型为标准型, 判别方阵是否可逆以及求解逆矩阵法. 内容虽然浅显简单, . , 敬请同仁批评指正.

设A 是一个, B 和r ×n 矩阵C 满足A =BC, 则称A =BC 为矩阵A 矩阵r m 和列数n. 如果n 阶方阵A 的秩r 远比其阶数n 小, 则满秩分. 设n 阶方阵A 有满秩分解A =BC, 根据文献[1], 我们有

n -r

|λE n -A |=|λE n -BC |==λ|λE r -CB |. 计算左边的行列式简单的多!

设A 为m ×n 矩阵, 通过行的初等变换, 可以将矩阵A 化为如下形式的简化行阶梯形矩阵

0000000

(1)

如果n 阶方阵A 的秩r 远远小于其阶数n , 无疑通过上式右边的行列式来求方阵A 的特征值要比直接

…0…0…0…0…0…0…0

1000000

c 1i 1+1…c 1i 2-1………………

000000

0100000

c 1i 2+1c 2i 2+1

…c 1i r -1…c 2i r -1……………

00000

0000100

c 1i r +1c 2i r +1c 3i r +1c r -1i r +1c r i r +1

………

c 1n c 2n c 3n

000000

00000

…c r -1n ………

c r n

00

00

其中1≤i 1≤i 2≤…≤i r ≤max (m , n ) .

根据列向量组的性质, 行初等变换不改变矩阵的列向量组的线性关系. 从矩阵A 的简化行阶梯形矩阵可以读出下面结论:

(i ) 矩阵A 的秩为r ;

 [收稿日期]2007204209

 [基金项目]国家自然科学基金(10671142) ; 教育部高校博士点基金([1**********])

196大 学 数 学              第25卷

(ii ) 矩阵A 的第i 1, i 2, …, i r 个列向量组成的向量组{A i 1, A i 2, …, A i r }为其列向量组的极大线性无

关组;

(iii ) 矩阵A 的每一个列向量均可由(ii ) 中的极大线性无关组线性表示, 并且第j 个列向量的线性表达式为

A j =c 1j A i 1+c 2j A i 2+…+c rj A i r .

从而, 矩阵A 可以表示为

c 11

A =(A 1 A 2 … A n ) =(A i 1 A i 2 … A i r )

c 21c r 1

c 12c 22c r 2

…c

1n …c 2n …c rn

.

若记B =(A i 1 A i 2 … A i r ) , 则矩阵B 为m ×r 列满秩矩阵; 若记

c 11c 12…c 1n

C =

c 21c r 1

c 22c r 2

…c 2n …c rn

,

则C 为r ×n 矩阵. 因为矩阵C 的第i 1列, 第i 2列, …, 第i r 列构成r 阶单位矩阵, C 为行满秩

矩阵. 这样, A =BC 就是矩阵A 的一个满秩分解.

注 上面讨论说明, 任意矩阵都存在满秩分解. 例1 经过行的初等变换, 我们有

231121310-11-1

=

31

3

2

3842131

732500

1-1110

01

1-11-11

-1-1111

-2-22-12

000

100

100

100

200

,

于是矩阵A 的秩为2. 令B =, C =, 则A =BC 为矩阵A 的一个满秩分解.

该例中矩阵A 的列向量极大线性无关组不唯一, 导致其满秩分解也不唯一. 取不同的列向量极大

线性无关组, 则得不同的满秩分解. 例如, 取第1列和第3列为极大列无关组, 则

[**************]

A =

131112253

12210

38411

3701

000

-1-111-1-11

-1-11

-2-22

000

100

100

100

20.

2

此时, 令B 1=

131

, C 1=

21

, 则A =B 1C 1为矩阵A 的又一个满秩分解. 再如, 取第1列

和第4列为极大列无关组, 则23151213A =

35281324

43712-1-11

1-1-11

3-1-11

3-2-22

1-1100

-2100

0100

-3200

.

000

000

第5期          靳全勤:初等变换的一个应用:矩阵的满秩分解

2

5381

111

, C 2=

197

此时, 令B 2=

131

10

-11

-21

0111

-32

, 则A =B 2C 2为矩阵A 的另一个满秩分解.

例2 求n 阶矩阵A =

…ω…

……11的特征值.

1

解 

容易看到, 矩阵A 的秩为1, 其满秩分解为A =

1

(1 1 … 1) , 由等式(1) 得|λE n -A |

n -1

(λ-n ) , 所以矩阵A 的特征值为λλ=λ1=λ2=…=λn -1=0, n =n.

本文只讨论了行初等变换的情形. 对于列初等变换, 有类似结果.

[参 考 文 献]

[1] 同济大学应用数学系. 高等代数与解析几何[M ].北京:of Two Full R ank Matrices :

of Elementary Matrix Operations

J I N Quan 2qi n

(Department of Mathematics , Tongji University , Shanghai 200092, China )

Abstract :We provide a simple method for decomposing a matrix into a product of two f ull rank matrices using elementary matrix operations.

K ey w ords :matrix ; elementary operation ; f ull rank decomposition


相关文章

  • 北航的考博矩阵真题
  • 北航2009年矩阵考博试题 试题共分三部分,第一题是选择题,5-6个:第二题是填空题,5-6个:第三题是大题(求解和证明) ,8个. 题量比较大,今年的题比较难,但题型以及考察的内容,跟往年差不多. 选择和填空题因为比较凌乱,所以没有回顾, ...查看


  • 811高等代数
  • 考试科目:811高等代数 复习要求: 要求考生熟练掌握高等代数的基本理论以及常用的技巧和方法,能够熟练地综合运用高等代数的理论和方法去求解和证明有关问题 二.主要复习内容: 1. 行列式 行列式的定义.性质和常用计算方法(如:三角化法.加边 ...查看


  • 毕业论文 关于矩阵求逆的几种方法
  • 矩阵求逆的几种方法 曹远军 (陕西理工学院数学与计算机科学学院数学与应用数学专业10级2班,陕西 汉中 723100) 指导老师:陈露 [摘要] 矩阵求逆是高等代数的一个重要内容.本文对矩阵求逆的方法进行了系统的归纳和总结方法, 这些求逆方 ...查看


  • 617 数学分析
  • 617 数学分析 三.考试形式一)试卷满分及考试时间 本试卷满分为150分,考试时间为180分钟. (二)答题方式 答题方式为闭卷.笔试.试卷由试题和答题纸组成,所有题目的答案必须写在答题纸相应的位置上.考生不得携带具有存储功能的计算器. ...查看


  • 线性方程组的解法
  • 线性方程组的解法 1 引言 在科学研究和大型工程设计中出现了越来越多的数学问题,而这些问题往往需要求数值解.在进行数值求解时,经离散后,常常归结为求解形如Ax= b的大型线性方程组.而如插值公式,拟合公式等的建立,微分方程差分格式的构造等, ...查看


  • 矩阵论考试要求
  • 矩阵论考试要求 第一章线性空间与内积空间 1.2 线性空的定义 1.3 基与坐标:求线性空间在指定基下的坐标,判别线性的线性关系,求极大线性无关组,秩,不同基下的坐标 1.4 线性子空间(定义,基与维数,span,子空间的和与交,直和) 1 ...查看


  • 新东方考研线性代数(第二讲 矩阵的乘法及逆运算)
  • 第二讲 矩阵的乘法及逆运算 一.乘法 1. 定义与规律 ① 乘法三要素 ⅰ)条件, ⅱ)AB 的行列数, ⅲ)AB 的元素. c ij =a i 1b 1j +a i 2b 2j +⋅⋅⋅+a in b nj ② 规律 (A 1+A 2) ...查看


  • 可逆矩阵在通信中的应用 1
  • 可逆矩阵及其在保密通信中的应用 摘 要 本文在可逆矩阵的定义.性质及求法的基础上,讨论了判断可逆矩阵的方法.分块可逆矩阵的求法以及可逆矩阵的一类求法,并通过实例给出了具体应用.介绍了保密通信及可逆矩阵在其中的应用. 关键词 矩阵理论:可逆矩 ...查看


  • 矩阵对角化及应用论文
  • 矩阵对角化及应用 理学院 数学082 缪仁东 指导师:陈巧云 摘 要:本文是关于矩阵对角化问题的初步研究, 对矩阵对角化充要条件的归纳, 总结, 通过对实对称矩阵, 循环矩阵, 特殊矩阵对角化方法的计算和研究, 让读者对矩阵对角化问题中求特 ...查看


热门内容