|
提供一个思路,希望对你有用:求出原矩阵的行向量组的极大无关组,然后用极大无关组组成的矩阵代替原行不满秩的矩阵即得到您想要的结果。
+ A" S }5 d6 A" j" w' o1 FMATLAB代码如下: Q' T5 M, \2 N9 i
M=size(H,1);
! M2 Z7 B: F. O/ Dwhile rank(H)~=M %当构造出的H不是行满秩时,找出它行向量组的极大无关组组成newH代替H/ ^0 J7 S6 W7 ]8 j' J/ V. x/ Z% L4 N
[R,jb]=rref(H'); %结果中R为rref函数中目标矩阵的行阶梯最简矩阵,jb是一个向量,为目标矩阵的列极大无关组所在的列数。注意这里需要先将H转置,因为我们的目标是要求H的行向量组的极大无关组
& y K' d1 Y2 knewH=H(jb,: );
. j, Q) v. ]3 g) r1 E3 CH=newH;
/ R8 _$ k+ R2 GM=size(H,1);# W9 I" A' }3 T; b
end" ~" v. @# j( ~; l/ j& L, ~
亲测有效,望有用。
: O* |; t. g2 u& ^& z
; I6 R" ^$ V+ o4 G T# h4 ^ F |
|