%矩阵的满秩分解
clc;clear;
%设输入矩阵为A
A = [1 1 2 1 0 0;0 1 1 0 1 0;1 2 3 0 0 1];%这里输入待分解的矩阵
disp("original matrix(m×n) is:");
disp(A);
A1 = rref(A); %将矩阵A化成行最简形式保存在A1中
disp("The simplest shape is:")
disp(A1)
[m,n]=size(A); %获取矩阵A的大小:m行n列
B0= [];%生成一个空向量
C0= [];%生成一个空向量
for i=1:m %依次扫描矩阵m行
flag=1;
for j=1:n%依次扫描矩阵n列
if A1(i,j)==1 %若A1(i, j)等于1
for k=1:i-1%固定j列,扫描此列的第1行到i-1行元素
if A1(k,j)~=0 %判断是否全为0
flag=0; %若不全为0,则将flag置为0(说明此列不是单位矩阵的列)
break;
end
end
for k=i+1:m%固定j列,扫描此列的第i+1行到m行(即最后一行)元素
if A1(k,j)~=0 %判断是否全为0
flag=0; %若不全为0,则将flag置为0(说明此列不是单位矩阵的列)
break;
end
end
if flag==1%若flag为1(不为0),则说明此列是【矩阵的行最简形式矩阵】的单位矩阵的列
B0=[B0,A(:,j)];%将矩阵A的j列加到B0列向量之后
C0=[C0;A1(i,:)];%将矩阵A1的i行加到C0行向量之后,
end
end
end
end
[m1,n1]=size(B0);%获取矩阵B0的大小:m1行n1列
[m2,n2]=size(C0); %获取矩阵C0的大小:m2行n2列
disp("m×r矩阵:");
B=B0(:,1:n1);%将矩阵B0的第1列到最后一列赋值给矩阵B (m×r)
disp(B);
disp("r×n矩阵:");
C=C0(1:m2,:);%将矩阵C0的第1行到最后一行赋值给矩阵C (r×n)
disp(C);
%验证:BC=A
disp("验证B×C:");
A_1= B*C;
disp(A_1);
满秩分解例题:
满秩分解例题2: