列主元消去法例题详解_Gauss列主元消去法函数

%列主元消去法解方程组Ax=b，实现PA=LU

function [x,detA] =gauss(A,b)

n=length(b);[p,q]=size(A);

if p~=q||p~=n      fprintf('方阵的维数不同，请重新输！'); %检错

end

%为提高运行速度，给L,U,x,c,d1赋初值

L=zeros(n,n);

U=zeros(n,n);

x=zeros(n,1);

c=zeros(1,n);

d1=0;

%按列选主元，并进行行交换，记录行信息

for i=1:n-1

max=abs(A(i,i));

m=i;

for j=i+1:n

if max

max=abs(A(j,i));

m=j;

end

end

if (m~=i)

for k=i:n

c(k)=A(i,k);

A(i,k)=A(m,k);

A(m,k)=c(k);

end

d1=b(i);

b(i)=b(m);

b(m)=d1;

end

%进行消元计算

for k=i+1:m

for j=i+1:n

A(k,j)=A(i,j)-A(i,j)*A(i,i);