列主消元法（Matlab实现）

目录

​​1、概述​​

​​2、代码​​

​​3、结果​​

1、概述

在前面一章知识点我们已经梳理清晰啦，这里直接上代码和运行结果。

2、代码

function [x]=abc(A,b,flag) 
if nargin<3,flag=0;end 
n=length(b); 
for k=1:(n-1) 
  [ap,p]=max(abs(A(k:n,k))); 
  p=p+k-1; 
  if p>k 
    A([k p],:)=A([p k],:); 
    b([k p],:)=b([p k],:); 
  end 
  m=A(k+1:n,k)/A(k,k); 
  A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-m*A(k,k+1:n); 
  b(k+1:n)=b(k+1:n)-m*b(k); 
  A(k+1:n,k)=zeros(n-k,1); 
  if flag~=0, Ab=[A,b],end 
end 
x=zeros(n,1); 
x(n)=b(n)/A(n,n); 
for k=n-1:-1:1 
  x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); 
end

3、结果