gaussl.m

来自「利用高斯消元法、雅克比消元法解决求解线性方程组的解的程序」· M 代码 · 共 19 行

function x=gaussl(A,b)
%本程序采用按列选主元方法消元。
n = length(b);
B = [A b];
for k = 1:n-1
   comp = abs(B(k:n,k));
   [M,I] = max(comp);   %%按列选主元，并记下行号（注意！是comp中的位置！！）。
   yhh = I+k-1;         %%将I换算为原行号yhh！！！ 
   h = B(k,:); B(k,:) = B(yhh,:); B(yhh,:) = h;   %%交换yhh行与k行。
   for i = (k+1):n
      for j = (k+1):(n+1)
         B(i,j) = B(i,j)-B(i,k)*B(k,j)/B(k,k);
      end
   end
end
%%此时B仅上三角元素改变，以下解方程也只用到了上三角部分。
%%以下解上三角方程。
x = uxb(B);