2.1.1 三角形解法

如前所述,直接法都归结为解三角形方程组,所以我们在这里首先讨论一般上三角形和下三角形方程组的解法。下面我们重点讨论四阶方程组,设:

Lxb为下三角形方程组,

假定(),解得:

上三角形方程Uxb的方程为:

同样假定(),解得:

对于n阶上三角形方程组,我们有以下公式:

求解一个n阶上三角形方程组的回代过程,用n次除法,系数和未知量相乘的次数为:,故回代过程所用的乘除次数总共为:

MTALAB程序

function x=backsub(A,b)

% 上三角方程组的回代

n=length(b);

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

例 A=[3 -2 1 -1;0 4 -1 2;0 0 2 3;0 0 0 5];b=[8 -3 11 15]';

backsub(A,b)

ans = 2

    -2

     1

     3

对于n阶下三角形方程组,我们有以下公式: