%%
A = [1,2,3;1,0,1;2,0,2;2,4,6];
b = [1;0;1;3];
Ab = [A,b]; % 构造增广矩阵
r1 = rank(A); % 求系数矩阵的秩
r2 = rank(Ab); % 求增广矩阵的秩
x = pinv(A)*b; % 求线性方程组的最小二乘解
%%
A = [1,-5,2,-3;5,3,6,-1;2,4,2,1];
b = [11;-1;-6];
c = [A,b];
% 化成行最简行
d = rref(c);
% 把数值矩阵化成符号矩阵,Matlab转成符号矩阵才能计算
e = sym(c);
f = rref(e);
%% 循环
clear,clc
N = 1e4;
result = [];
for n