矩阵的满秩分解知识讲解
- 格式:doc
- 大小:223.00 KB
- 文档页数:4
§矩阵的满秩分解本节讨论一个n m ⨯复矩阵A 可以分解为两个与A 的秩相同的矩阵之积的问题。
定义设n m ⨯复矩阵A 的秩为r ,如果存在两个与A 的秩相同的复矩阵F 与G ,使得FG A =,则称此式为复矩阵A 的满秩分解。
当A 是满秩矩阵时(行满秩或列满秩)A 可以分解为单位矩阵与A 自身的乘积,这个满秩分解叫做平凡分解。
定理设n m ⨯复矩阵A 的秩为r 0>,则A 有满秩分解。
证:因为0>=r rankA ,对A 施行初等行变换,可得到阶梯形矩阵⎪⎪⎭⎫ ⎝⎛=0G B , 其中G 为n r ⨯矩阵,并且0>=r rankG ;因此存在着有限个m 阶初等矩阵之积,记作P ,有B PA =,或者B P A 1-=,将矩阵1-P 分块为()S F P =-1 ,其中F 为r m ⨯矩阵,S 为)(r n m -⨯矩阵,并且r rankF =,r n rankS -=。
则有()FG G S F B P A =⎪⎪⎭⎫ ⎝⎛==-01 ,其中F 是列满秩矩阵,S 是行满秩矩阵。
▌但是,矩阵A 的满秩分解不唯一。
这是因为若取任意一个r 阶非奇异矩阵D ,则有G F G D FD FG A ~~))((1===-。
例1、 求矩阵⎪⎪⎪⎭⎫ ⎝⎛----=122211212101A 的满秩分解。
解:对矩阵A 进行初等行变换()⎪⎪⎭⎫ ⎝⎛==⎪⎪⎪⎭⎫ ⎝⎛--→⎪⎪⎪⎭⎫ ⎝⎛----=0111000001130200012101100122201011210012101G B I A 其中⎪⎪⎭⎫ ⎝⎛-=30202101G 所以⎪⎪⎪⎭⎫ ⎝⎛-=000030202101B ,⎪⎪⎪⎭⎫ ⎝⎛-=111011001P ;而()S F P =⎪⎪⎪⎭⎫ ⎝⎛--=-1120110011,其中⎪⎪⎪⎭⎫ ⎝⎛--=121101F 由此可见,所以有()⎪⎪⎪⎭⎫ ⎝⎛--==⎪⎪⎭⎫ ⎝⎛==-12110101FG G S F B P A ⎪⎪⎭⎫ ⎝⎛-30202101。
满秩矩阵及满秩矩阵的应用专业:通信与信息系统姓名:李娜学号:6120140151目录一、满秩矩阵及满秩矩阵在矩阵分解方面的应用 (2)1.1矩阵的秩 (2)1.2满秩矩阵 (2)1.3满秩矩阵的性质 (3)1.3.1行(列)矩阵的一些性质 (4)1.4 行(列) 满秩矩阵在矩阵分解方面的应用 (6)二、满秩矩阵在保密通信中的应用 (8)2.1 基于满秩矩阵的保密通信模型 (8)2.1.1加密保密通信模型 (8)2.2.2满秩矩阵的应用 (8)2.2密钥的生成 (10)2.2.1加密密钥的生成 (10)2.2.2解密密钥的生成 (10)2.3其它问题 (10)2.3.1明文矩阵的选择 (10)2.3.2加密矩阵的选择 (11)2.3.3算法优化 (11)一、满秩矩阵及满秩矩阵在矩阵分解方面的应用引言矩阵是数学中的一个重要的基本概念,是现代数学的一个主要研究对象,也是数学研究和应用的一个重要工具。
“矩阵”这个词是由西尔维斯特首先使用的,他是为了将数学的矩形阵列区别于行列式而发明了这个述语,而实际上,矩阵这个课题在诞生之前就已经发展的很好了。
1.1矩阵的秩设A是一组向量,定义A的最大无关组中向量的个数为A的秩。
定义1 在m n矩阵A中,任意决定k行和k列交叉点上的元素构成A的一个k阶子矩阵,此子矩阵的行列式,称为A的一个k阶子式。
例如,在阶梯形矩阵中,选定1,3行和3,4列,它们交叉点上的元素所组成的2阶子矩阵的行列式就是矩阵A的一个2阶子式。
定义2 A=(a ij)m×n的不为零的子式的最大阶数称为矩阵A的秩,记作r(A),或rank(A)或R(A)。
特别规定零矩阵的秩为零。
显R(A)≤min(m,n)易得:若A中至少有一个r阶子式不等于零,且在R(A)<min(m,n)时,A中所有的r+1阶子式全为零,则A的秩为r。
由定义直接可得n阶可逆矩阵的秩为n,通常又将可逆矩阵称为满秩矩阵,不满秩矩阵就是奇异矩阵,det(A)=0。
第十一讲 满秩分解与奇异值分解一、矩阵的满秩分解1. 定义:设m n r A C (r 0)⨯∈>,若存在矩阵m r r F C ⨯∈及r nrG C ⨯∈,使得 A FG =,则称其为A 的一个满秩分解。
说明:(1)F 为列满秩矩阵,即列数等于秩;G 为行满秩矩阵,即行数等于秩。
(2)满秩分解不唯一。
r rrD C ⨯∀∈(r 阶可逆方阵),则 1111A FG F(DD )G (FD)(D G)F G --====,且m r r n1r 1rF C ,G C ⨯⨯∈∈ 2. 存在性定理:任何非零矩阵均存在满秩矩阵证明:采用构造性证明方法。
设m nr A C ⨯∈,则存在初等变换矩阵m mmE C ⨯∈, 使 G r EA B .......O (m r)⎡⎤⎢⎥==⎢⎥-⎢⎥⎣⎦行行, 其中r nr G C ⨯∈ 将A 写成1A E B -=,并把1E -分块成[]1r (m r)E F |S --=列列,其中m rrF C ⨯∈ .G A F .S ....FG .O ⎡⎤⎡⎤⎢⎥⎢⎥∴==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦E 是满秩分解。
3. Hermite 标准形(行阶梯标准形)设m nr B C (r 0)⨯∈>,且满足(1) B 的前r 行中每一行至少含一个非零元素(称为非零行),且第一个非零元素为1,而后(m r)-行的元素全为零(称为零行);(2) 若B 中第i 行的第一个非零元素(即1)在第i j 列(i 1,2,...,r)=,则 12r j j ...j <<<;(3) 矩阵B 的第1j 列,第2j 列,…,第r j 列合起来恰为m 阶单位方阵m I 的前r 列(即12r j ,j ,...,j 列上除了前述的1外全为0)则称B 为Hermite 标准形。
例1 561356120013001022B C 000111000000000000⨯⨯-⎡⎤⎢⎥⎢⎥=∈-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 为Hermite 标准形452245010200013B C 0000000000⨯⨯⎡⎤⎢⎥⎢⎥=∈⎢⎥⎢⎥⎣⎦ 也是Hermite 标准形4. 满秩分解的一种求法设m nr A C ⨯∈,(1) 采用行初等变换将A 化成Hermite 标准形,其矩阵形式为EA B =,其中B 为Hermite 标准形定义中给出的形状;(2) 选取置换矩阵1 P 的第i 列为i j e ,即该列向量除第i j 个元素为1外,其余元素全为零(i 1,2,...,r)=,其中i j 为Hermite 标准形中每行第一个非零元素(即1)所在的列数;2 其它(n r)-列只需确保P 为置换矩阵即可(P 的每一行,每一列均只有一个非零元素,且为1);3 用P 右乘任何矩阵(可乘性得到满足时),即可得该矩阵的第i j 列置换到新矩阵(即乘积矩阵)的第i 列 4 令[]1r (n r)P P |*-=列列,即12r n r1j j j rn rP e e ...e C ⨯⨯⎡⎤=∈⎣⎦(3)令G B =的前r 行r n n C ⨯∈,m r1rF AP C ⨯=∈则A FG = 证明:G EA B O ⎡⎤==⎢⎥⎣⎦,[]1G A E B F |S FG O -⎡⎤===⎢⎥⎣⎦则m r r F C ⨯∈,r nrG C ⨯∈,G 已知,但F ?=,当然可以通过求出1E,E -再将1E -分块得到,但这样G 就没必要采用Hermite 标准形形式,注意到r 1I BP O ⎡⎤=⎢⎥⎣⎦,则[]1r 11I AP E BP F |S F O -⎡⎤===⎢⎥⎣⎦证毕例1 1230A 02111021⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦求其满秩分解解:(1)首先求出A 的秩。
满秩矩阵及满秩矩阵的应用专业:通信与信息系统姓名:李娜学号:6120140151目录一、满秩矩阵及满秩矩阵在矩阵分解方面的应用 (2)1.1矩阵的秩 (2)1.2满秩矩阵 (2)1.3满秩矩阵的性质 (3)1.3.1行(列)矩阵的一些性质 (4)1.4 行(列) 满秩矩阵在矩阵分解方面的应用 (6)二、满秩矩阵在保密通信中的应用 (8)2.1 基于满秩矩阵的保密通信模型 (8)2.1.1加密保密通信模型 (8)2.2.2满秩矩阵的应用 (8)2.2密钥的生成 (10)2.2.1加密密钥的生成 (10)2.2.2解密密钥的生成 (10)2.3其它问题 (10)2.3.1明文矩阵的选择 (10)2.3.2加密矩阵的选择 (11)2.3.3算法优化 (11)一、满秩矩阵及满秩矩阵在矩阵分解方面的应用引言矩阵是数学中的一个重要的基本概念,是现代数学的一个主要研究对象,也是数学研究和应用的一个重要工具。
“矩阵”这个词是由西尔维斯特首先使用的,他是为了将数学的矩形阵列区别于行列式而发明了这个述语,而实际上,矩阵这个课题在诞生之前就已经发展的很好了。
1.1矩阵的秩设A是一组向量,定义A的最大无关组中向量的个数为A的秩。
定义1 在m n矩阵A中,任意决定k行和k列交叉点上的元素构成A的一个k阶子矩阵,此子矩阵的行列式,称为A的一个k阶子式。
例如,在阶梯形矩阵中,选定1,3行和3,4列,它们交叉点上的元素所组成的2阶子矩阵的行列式就是矩阵A的一个2阶子式。
定义2 A=(a ij)m×n的不为零的子式的最大阶数称为矩阵A的秩,记作r(A),或rank(A)或R(A)。
特别规定零矩阵的秩为零。
显R(A)≤min(m,n)易得:若A中至少有一个r阶子式不等于零,且在R(A)<min(m,n)时,A中所有的r+1阶子式全为零,则A的秩为r。
由定义直接可得n阶可逆矩阵的秩为n,通常又将可逆矩阵称为满秩矩阵,不满秩矩阵就是奇异矩阵,det(A)=0。
矩阵的满秩分解
§4.3矩阵的满秩分解
本节讨论一个n m ⨯复矩阵A 可以分解为两个与A 的秩相同的矩阵之积的问题。
定义4.3.1设n m ⨯复矩阵A 的秩为r ,如果存在两个与A 的秩相同的复矩阵F 与G ,使得FG A =,则称此式为复矩阵A 的满秩分解。
当A 是满秩矩阵时(行满秩或列满秩)A 可以分解为单位矩阵与A 自身的乘积,这个满秩分解叫做平凡分解。
定理4.3.1设n m ⨯复矩阵A 的秩为r 0>,则A 有满秩分解。
证:因为0>=r rankA ,对A 施行初等行变换,可得到阶梯形矩阵⎪⎪⎭
⎫ ⎝⎛=0G B , 其中G 为n r ⨯矩阵,并且0>=r rankG ;因此存在着有限个m 阶初等矩阵之积,
记作P ,有B PA =,或者B P A 1-=,将矩阵1-P 分块为()S F P =-1 ,其中F 为r m ⨯矩阵,S 为)(r n m -⨯矩阵,并且r rankF =,r n rankS -=。
则有()FG G S F B P A =⎪⎪⎭
⎫ ⎝⎛==-01 ,其中F 是列满秩矩阵,S 是行满秩矩阵。
▌
但是,矩阵A 的满秩分解不唯一。
这是因为若取任意一个r 阶非奇异矩阵D ,则有
G F G D FD FG A ~~))((1===-。
例1、 求矩阵⎪⎪⎪⎭
⎫ ⎝⎛----=122211212101A 的满秩分解。
解:对矩阵A 进行初等行变换
()⎪⎪⎭⎫ ⎝⎛==⎪⎪⎪⎭
⎫ ⎝⎛--→⎪⎪⎪⎭⎫ ⎝⎛----=0111000001130200012101100122201011210012101G B I A 其中⎪⎪⎭⎫ ⎝⎛-=30202101G 所以⎪⎪⎪⎭⎫ ⎝⎛-=000030202101B ,⎪⎪⎪⎭
⎫ ⎝⎛-=111011001P ;而()S F P =⎪⎪⎪⎭⎫ ⎝⎛--=-1120110011,其中⎪⎪⎪⎭
⎫ ⎝⎛--=121101F
由此可见,所以有()⎪⎪⎪⎭
⎫ ⎝⎛--==⎪⎪⎭⎫ ⎝⎛==-12110101FG G S F B P A ⎪⎪⎭⎫ ⎝⎛-30202101。
定义4.3.2设n m ⨯复矩阵H 的秩为r ()0>r ,并且满足以下条件:
1)矩阵H 的前r 行中的每一行至少含有一个不为零的元素,并且第一个不为零的元素是1,而后r m -行的元素均为零;
2)如果矩阵H 的第i 行的第一个不为零的元素1在第i j 列()r i ,,2,1 =, 则r j j j <<< 21;
3)矩阵H 的r j j j ,,,21 列是单位矩阵m I 的前r 列; 则称矩阵H 为Hermite 标准形(最简型)。
由此定义可见,对于任意一个秩为r 的n m ⨯复矩阵A ,均可以经过初等行变换将其化为Hermite 标准形H ,而且矩阵H 的前r 列元素组成的列向量组线性无关。
定义4.3.3以n 阶单位矩阵n I 的n 个列向量n e e e ,,,21 为列构成的n 阶矩阵()n j j j e e e P ,,,21 =叫做置换矩阵。
其中n j j j ,,,21 是n ,,2,1 的一个全排列。
定理4.3.2设n m ⨯复矩阵A 的秩为r ()0>r ,矩阵A 的Hermite 标准形为H ,则在矩阵A 的满秩分解FG A =中,可以取矩阵F 为A 的r j j j ,,,21 列构成的r m ⨯列矩阵,G 为H 的前r 行构成的n r ⨯列矩阵。
例2、求矩阵⎪⎪⎪⎭
⎫ ⎝⎛----=122211212101A 的满秩分解。
解:先求出矩阵A 的Hermite 标准形
H A =⎪⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎪⎭⎫ ⎝⎛----=0000230102101122211212101,H 的第1列与第2列构成3I 的
前两列,所以矩阵F 为A 的第1列与第2列构成的23⨯矩阵,G 为H 的前2行
构成的42⨯矩阵,即⎪⎪⎪⎭
⎫ ⎝⎛-=222101F ,⎪⎪⎭⎫ ⎝⎛-=230102101G , 所以⎪⎪⎪⎭
⎫ ⎝⎛-==222101FG A ⎪⎪⎭⎫ ⎝⎛-230102101。
对比例1,可以看出矩阵A 的满秩分解不唯一。