matlab列主元高斯消去法程序,高斯消元法,列主元素消元法及LU分解法的matlab程序...

高斯消元法,列主元素消元法及LU分解法的matlab程序解析及举例

§2.2.1高斯消元法的MATLAB程序

function [RA,RB,n,X]=gaus(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

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

运行命令及结果

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] =gaus (A,b)

请注意：因为RA=RB=n，所以此方程组有唯一解.

RA =

3

RB =

3

n =

3

X =