列主元消去法MATLAB程序

列主元消去法MATLAB程序：

function [RA,RB,n,X]=liezhu(A,b)

B=[A b]; n=length(b); RA=rank(A);

RB=rank(B);zhica=RB-RA;

if zhica>0,

disp('Çë×¢Ò⣺ÒòΪRA~=RB£¬ËùÒÔ´Ë·½³Ì×éÎÞ½â.') return

end

if RA==RB

if RA==n

disp('Çë×¢Ò⣺ÒòΪRA=RB=n£¬ËùÒÔ´Ë·½³Ì×éÓÐΨһ½â.') X=zeros(n,1); C=zeros(1,n+1);

for p= 1:n-1

[Y,j]=max(abs(B(p:n,p))); C=B(p,:);

B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;

for k=p+1:n

m= B(k,p)/ B(p,p);

B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);

end

end

b=B(1:n,n+1);A=B(1:n,1:n);

X(n)=b(n)/A(n,n);

for q=n-1:-1:1

X(q)=(b(q)-

sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);

end

else

disp('Çë×¢Ò⣺ÒòΪRA=RB

´Ë·½³Ì×éÓÐÎÞÇî¶à½â.')

end

end

列主元消除法，M文件

a=[2.51 1.48 4.53;1.48 0.93 -1.30 ;2.68 3.04 -1.48];

b=[0.05;1.03;-0.53];

[RA,RB,n,X]=liezhu(a,b)

运行结果；