矩阵的满秩分解

§4.3矩阵的满秩分解

本节讨论一个m ⨯n 复矩阵A 可以分解为两个与A 的秩相同的矩阵之积的问题。

定义4.3.1设m ⨯n 复矩阵A 的秩为r ,如果存在两个与A 的秩相同的复矩阵F 与G ,使得A =FG ,则称此式为复矩阵A 的满秩分解。

当A 是满秩矩阵时(行满秩或列满秩)A 可以分解为单位矩阵与A 自身的乘积,这个满秩分解叫做平凡分解。

定理4.3.1设m ⨯n 复矩阵A 的秩为r >0,则A 有满秩分解。

证:因为rankA =r >0,对A 施行初等行变换,可得到阶梯形矩阵B = 0⎪⎪, ⎝⎭

其中G 为r ⨯n 矩阵,并且rankG =r >0;因此存在着有限个m 阶初等矩阵之积,

-1-1记作P ,有PA =B ,或者A =P B ,将矩阵P 分块为P -1=(F S ) ,其中F 为m ⨯r ⎛G ⎫

矩阵,S 为m ⨯(n -r ) 矩阵,并且rankF =r ,rankS =n -r 。

则有A =P -1B =(F S ) ⎛G ⎫⎪=FG ,其中F 是列满秩矩阵,S 是行满秩矩阵。 ▌ ⎪⎝0⎭

但是,矩阵A 的满秩分解不唯一。这是因为若取任意一个r 阶非奇异矩阵D ,则有

~~A =FG =(FD )(D -1G ) =F G 。

2⎫⎛-101 ⎪例1、 求矩阵A = 12-11⎪的满秩分解。

22-2-1⎪⎝⎭

解:对矩阵A 进行初等行变换

2100⎫⎛-1012100⎫⎛-101 ⎪ ⎪G ⎫ (A I )= 12-11010⎪→ 0203110⎪=B =⎛ 0⎪⎪ ⎝⎭ 22-2-1001⎪ 00001-11⎪⎝⎭⎝⎭

⎛-1012⎫⎛100⎫ ⎪ ⎪-1012⎛⎫其中G = 0203⎪⎪所以B = 0203⎪,P = 110⎪;而⎝⎭ 0000⎪ 1-11⎪⎝⎭⎝⎭

P -1

⎛100⎫⎛10⎫ ⎪ ⎪= -110⎪=(F S ),其中F = -11⎪ -211⎪ -21⎪⎝⎭⎝⎭

⎛10⎫ ⎪⎛-1012⎫G ⎛⎫-1⎪=FG =-11由此可见,所以有A =P B =(F S ) ⎪ 0⎪ 0203⎪⎪。 ⎝⎭⎭ -21⎪⎝⎝⎭

109

定义4.3.2设m ⨯n 复矩阵H 的秩为r (r >0),并且满足以下条件:

1)矩阵H 的前r 行中的每一行至少含有一个不为零的元素,并且第一个不为零的元素是1,而后m -r 行的元素均为零;

2)如果矩阵H 的第i 行的第一个不为零的元素1在第j i 列(i =1, 2, , r ),

则j 1

3)矩阵H 的j 1, j 2, , j r 列是单位矩阵I m 的前r 列;

则称矩阵H 为Hermite 标准形(最简型) 。

由此定义可见,对于任意一个秩为r 的m ⨯n 复矩阵A ,均可以经过初等行变换将其化为Hermite 标准形H ,而且矩阵H 的前r 列元素组成的列向量组线性无关。

定义4.3.3以n 阶单位矩阵I n 的n 个列向量e 1, e 2, , e n 为列构成的n 阶矩阵P =e j 1, e j 2, , e j n 叫做置换矩阵。其中j 1, j 2, , j n 是1, 2, , n 的一个全排列。

定理4.3.2设m ⨯n 复矩阵A 的秩为r (r >0),矩阵A 的Hermite 标准形为H ,则在矩阵A 的满秩分解A =FG 中,可以取矩阵F 为A 的j 1, j 2, , j r 列构成的m ⨯r 列矩阵,G 为H 的前r 行构成的r ⨯n 列矩阵。 ()

2⎫⎛-101 ⎪例2、求矩阵A = 12-11⎪的满秩分解。

22-2-1⎪⎝⎭

解:先求出矩阵A 的Hermite 标准形

2⎫⎛10-12⎫⎛-101 ⎪ ⎪A = 12-11⎪→ 0102⎪=H ,H 的第1列与第2列构成I 3的前两

22-2-1⎪ 0000⎪⎝⎭⎝⎭

列,所以矩阵F 为A 的第1列与第2列构成的3⨯2矩阵,G 为H 的前2行构成的2⨯4矩

⎛-10⎫ ⎪⎛10-12⎫阵,即F = 12⎪,G = 0102⎪⎪, ⎝⎭ 22⎪⎝⎭

⎛-10⎫ ⎪⎛10-12⎫所以A =FG = 12⎪ 0102⎪⎪。

⎭ 22⎪⎝⎝⎭

对比例1,可以看出矩阵A 的满秩分解不唯一。

110

§4.3矩阵的满秩分解

本节讨论一个m ⨯n 复矩阵A 可以分解为两个与A 的秩相同的矩阵之积的问题。

定义4.3.1设m ⨯n 复矩阵A 的秩为r ,如果存在两个与A 的秩相同的复矩阵F 与G ,使得A =FG ,则称此式为复矩阵A 的满秩分解。

当A 是满秩矩阵时(行满秩或列满秩)A 可以分解为单位矩阵与A 自身的乘积,这个满秩分解叫做平凡分解。

定理4.3.1设m ⨯n 复矩阵A 的秩为r >0,则A 有满秩分解。

证:因为rankA =r >0,对A 施行初等行变换,可得到阶梯形矩阵B = 0⎪⎪, ⎝⎭

其中G 为r ⨯n 矩阵,并且rankG =r >0;因此存在着有限个m 阶初等矩阵之积,

-1-1记作P ,有PA =B ,或者A =P B ,将矩阵P 分块为P -1=(F S ) ,其中F 为m ⨯r ⎛G ⎫

矩阵,S 为m ⨯(n -r ) 矩阵,并且rankF =r ,rankS =n -r 。

则有A =P -1B =(F S ) ⎛G ⎫⎪=FG ,其中F 是列满秩矩阵,S 是行满秩矩阵。 ▌ ⎪⎝0⎭

但是,矩阵A 的满秩分解不唯一。这是因为若取任意一个r 阶非奇异矩阵D ,则有

~~A =FG =(FD )(D -1G ) =F G 。

2⎫⎛-101 ⎪例1、 求矩阵A = 12-11⎪的满秩分解。

22-2-1⎪⎝⎭

解:对矩阵A 进行初等行变换

2100⎫⎛-1012100⎫⎛-101 ⎪ ⎪G ⎫ (A I )= 12-11010⎪→ 0203110⎪=B =⎛ 0⎪⎪ ⎝⎭ 22-2-1001⎪ 00001-11⎪⎝⎭⎝⎭

⎛-1012⎫⎛100⎫ ⎪ ⎪-1012⎛⎫其中G = 0203⎪⎪所以B = 0203⎪,P = 110⎪;而⎝⎭ 0000⎪ 1-11⎪⎝⎭⎝⎭

P -1

⎛100⎫⎛10⎫ ⎪ ⎪= -110⎪=(F S ),其中F = -11⎪ -211⎪ -21⎪⎝⎭⎝⎭

⎛10⎫ ⎪⎛-1012⎫G ⎛⎫-1⎪=FG =-11由此可见,所以有A =P B =(F S ) ⎪ 0⎪ 0203⎪⎪。 ⎝⎭⎭ -21⎪⎝⎝⎭

109

定义4.3.2设m ⨯n 复矩阵H 的秩为r (r >0),并且满足以下条件:

1)矩阵H 的前r 行中的每一行至少含有一个不为零的元素,并且第一个不为零的元素是1,而后m -r 行的元素均为零;

2)如果矩阵H 的第i 行的第一个不为零的元素1在第j i 列(i =1, 2, , r ),

则j 1

3)矩阵H 的j 1, j 2, , j r 列是单位矩阵I m 的前r 列;

则称矩阵H 为Hermite 标准形(最简型) 。

由此定义可见,对于任意一个秩为r 的m ⨯n 复矩阵A ,均可以经过初等行变换将其化为Hermite 标准形H ,而且矩阵H 的前r 列元素组成的列向量组线性无关。

定义4.3.3以n 阶单位矩阵I n 的n 个列向量e 1, e 2, , e n 为列构成的n 阶矩阵P =e j 1, e j 2, , e j n 叫做置换矩阵。其中j 1, j 2, , j n 是1, 2, , n 的一个全排列。

定理4.3.2设m ⨯n 复矩阵A 的秩为r (r >0),矩阵A 的Hermite 标准形为H ,则在矩阵A 的满秩分解A =FG 中,可以取矩阵F 为A 的j 1, j 2, , j r 列构成的m ⨯r 列矩阵,G 为H 的前r 行构成的r ⨯n 列矩阵。 ()

2⎫⎛-101 ⎪例2、求矩阵A = 12-11⎪的满秩分解。

22-2-1⎪⎝⎭

解:先求出矩阵A 的Hermite 标准形

2⎫⎛10-12⎫⎛-101 ⎪ ⎪A = 12-11⎪→ 0102⎪=H ,H 的第1列与第2列构成I 3的前两

22-2-1⎪ 0000⎪⎝⎭⎝⎭

列,所以矩阵F 为A 的第1列与第2列构成的3⨯2矩阵,G 为H 的前2行构成的2⨯4矩

⎛-10⎫ ⎪⎛10-12⎫阵,即F = 12⎪,G = 0102⎪⎪, ⎝⎭ 22⎪⎝⎭

⎛-10⎫ ⎪⎛10-12⎫所以A =FG = 12⎪ 0102⎪⎪。

⎭ 22⎪⎝⎝⎭

对比例1,可以看出矩阵A 的满秩分解不唯一。

110


相关文章

  • 第三节矩阵的基本运算
  • 第三节 矩阵的基本运算 §3.1 加和减 §3.2矩阵乘法 §3.2.1 矩阵的普通乘法 §3.2.2 矩阵的Kronecker 乘法 §3.3 矩阵除法 §3.4矩阵乘方 §3.5 矩阵的超越函数 §3.6数组运算 §3.6.1数组的加和 ...查看


  • 矩阵的分解
  • 1:三角(LU)分解. 如果方阵A可以分解成为一个下三角矩阵L和一个上三角矩阵U 的乘积,则称A可作三角(LU)分解.当L为单位下三角矩阵而U为上三角矩阵时,此三角分解称为杜利特(Doolittle)分解.当L为下三角矩阵而U为单位上三角矩 ...查看


  • 大规模稀疏矩阵PARD解方法介绍
  • 高性能计算 13 大规模稀疏矩阵PARDISO求解方法介绍 于 超 英特尔亚太研发中心 上海 200241 chao.yu@intel.com 摘要: 大规模稀疏矩阵的求解是高性能计算中的一个常见问题.本文介绍了用直接法(Direct Sp ...查看


  • MATLAB求矩阵分解,特征值与特征向量
  • (1) LU分解 A是非奇异的,LU分解总是可以进行的. [L,U]=lu(X):产生一个上三角阵U和一个变换形式的下三角阵L(行交换),矩阵X必须是方阵. [L,U,P]=lu(X):产生一个上三角阵U和一个下三角阵L以及一个置换矩阵P, ...查看


  • 计算机科学计算答案 第二章 矩阵变换和计算
  • 第二章 矩阵变换和计算 一.内容提要 本章以矩阵的各种分解变换为主要内容,介绍数值线性代数中的两个基本问题:线性方程组的求解和特征系统的计算,属于算法中的直接法.基本思想为将计算复杂的一般矩阵分解为较容易计算的三角形矩阵. 要求掌握Gaus ...查看


  • MATLAB小波图像分解
  • 上一篇文章中我们实现了小波的一维.二维信号分解与重构,其中的二维信号分解与重构,只要稍作修改,就可以实现图像的分解和重构了.修改的工作,主要是对图像信号进行规范化处理.数据格式转换和绘图细节处理等. 简单起见,我们从黑白(灰度)图像的分解. ...查看


  • 初等变换的一个应用矩阵的满秩分解
  • 第25卷第5期2009年10月 大 学 数 学 COLL EGE MA T H EMA TICS Vol. 25, №. 5 Oct. 2009 初等变换的一个应用:矩阵的满秩分解 靳全勤 (同济大学数学系, 上海200092) [摘 要] ...查看


  • 矩阵分解及无约束最优化方法
  • 矩阵分解及无约束最优化 方法的原理和应用简介 --最优化方法课程实验报告 学 院:数学与统计学院 班 级:硕2041班 姓 名:王彭 学 号:指导教师:阮小娥 同 组 人:陈莹 钱东东 矩阵分解及无约束最优化方法的原理和应用简介 矩阵分解及 ...查看


  • 矩阵的三角分解的编程实现
  • 矩阵的三角分解的编程实现 一. 课程设计目的: 1. 将高斯消去法改写为紧凑形式,可以直接从矩阵A 的元素得到计算L,U 元素的递推公式,而不需要任何中间的步骤. 2. 理解三角分解的具体方法. 3. 掌握对于矩阵的运算的c 语言代码.了解 ...查看


热门内容