实验二矩阵的秩分解.docVIP

实验二 矩阵的秩分解 1、原理 设A是m×n复矩阵，rankA r≥1，则必存在秩为r的两个矩阵P1∈Cm×r，Q1∈Cr×n，使得A P1Q1 2、算法 在线性代数中，我们知道行最简形矩阵是一种特殊形式的矩阵，这种矩阵具有如下特征： （1）元素全为0的行在非零行的下方。 （2）非零行中从左到右的第一个非零元是1，称为该行的首1元素。首1元是所在列的惟一非零元。 （3）首1元左下方的元素全为0. 任一矩阵必可以通过有限次行的初等变换化成最简形。 设A是m×n矩阵，rankA r，如果A的最简行矩阵H的首1元分别在H的第j1，j2，…，jr列，取H的前r行构成的矩阵为C，A的第j1，j2，…，jr列构成的矩阵为B，则A BC就是矩阵A的一个秩分解。 3、程序 在matlab中并没有像LU分解（尽管和咱们学的不太一样）那样的现成的程序，因此在做实验室借助于网络，搜索到一个程序，在matlab中建立一个M文件，经试验，确实实现了矩阵的秩分解。程序如下： function zhifenjie A B rref A ; [m,n] size A ; P0 1:m,: 0; Q0 :,1:n 0; for i 1:m flag 1; for j 1:n if B i,j 1 for k 1:i-1 if B k,j ~ 0 flag 0; break; end end for k i+1:m if B k,j ~ 0 flag 0; break; end end if flag 1 P0 [P0,A :,j ]; Q0 [Q0;B i,: ]; end end end end [m1,n1] size P0 ; [m2,n2] size Q0 ; P P0 :,2:n1 Q Q0 2:m2,: 4、例子 以教材上的矩阵A [1 -1 2 3；-1 0 -1 0；0 -1 1 3]为例说明上述程序的正确性，在matlab中输入以下指令： 5、结果 输入上述指令以后，得到如下的结果： 此时可以说，上述程序是正确的。