行列式求解，高斯消元法

时间 2019-11-12

标签行列式求解高斯消元繁體版

原文原文链接

程序部分

gauss.m shell

function [x]=gauss(A,b)
  n=size(A,1);
  for k=1:n-1
    m=A(k+1:n,k)/A(k,k);
    A(k+1:n,k+1:n)-=m*A(k,k+1:n);
    A(k+1:n,k)=zeros(n-k,1);
    b(k+1:n)-=m*b(k);
  endfor
  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);
  endfor
endfunction

参数： A 输入系数矩阵 spa

b 右端项 .net

程序原理就是利用基础行列式变换构造一个三角形矩阵，而后三角形矩阵按照逆序分别求出x4 x3 x2 x1获得原始解。 code

示例

输入
A =

   1   1   1   1
  -1   2  -3   1
   3  -3   6  -2
  -4   5   2  -3

b =

   10
   -2
    7
    0

输出
ans =

   1
   2
   3
   4

局限性是 :高斯消去法没法处理输入矩阵A(i,i)=0的状况，并且使用过程当中不稳定因素也比较大，属于基本的求值法。另一种改进版的高斯消元方法是将主元进行排序而后再消元，能够解决传统方法的弊端。列主元高斯消元法 blog