列主元消去法例题详解_MATLAB列主元消去法

functionx=gauss(A,b)intAB=[A,b];n=length(b);RA=rank(A);RB=rank(B);d=RB-RA;ifd>0disp('此方程组无解')endifRA==RB&RA~=ndisp('此方程组有无穷解')endifRA==RB&RA==ndisp('此方程组...

function x = gauss(A,b)

int A

B=[A,b];

n=length(b);

RA=rank(A);

RB=rank(B);

d=RB-RA;

if d>0

disp('此方程组无解')

end

if RA==RB&RA~=n

disp('此方程组有无穷解')

end

if RA==RB&RA==n

disp('此方程组有唯一解')

[n,n]=size(A);

x=zeros(n,1);

for k = 1:n-1

[piv,r] = max(abs(B(k:n,k))); %找列主元所在子矩阵的行r

r = r + k - 1; % 列主元所在大矩阵的行 ,

if r>k

temp=B(k,:);

B(k,:)=B(r,:);

B(r,:)=temp;

end

end

if