带状对角矩阵的LU分解及回代求解算法实现

算法名称： 带状对角矩阵的LU 分解及回代求解

 

算法描述：

       分解主要是使用笔者前面几篇文章提到过的Crout 方法。因为不可能把一个带状对角矩阵 A 的LU 分解也像其压缩形式本是一样紧凑的存储起来，因为分解产生了附加的非零元素填入。一种直接的存储方案是，坝上三角因子（U ）返回到以前占有的相同的空间中，把下三角因子（L ）返回到单独的 N×m₁ 压缩矩阵中。U 的对角线元素被存放在 A 的存储空间的第一列。

       一旦矩阵A 被分解，任意数目的右端项都可以通过重复调用回代函数求解，方法参考笔者前面的几篇文章。

 

运行示例：

 

The origin matrix:

3.0 1.0 0.0 0.0 0.0 0.0 0.0

4.0 1.0 5.0 0.0 0.0 0.0 0.0

9.0 2.0 6.0 5.0 0.0 0.0 0.0

0.0 3.0 5.0 8.0 9.0 0.0 0.0

0.0 0.0 7.0 9.0 3.0 2.0 0.0

0.0 0.0 0.0 3.0 8.0 4.0 6.0

0.0 0.0 0.0 0.0 2.0 4.0 4.0

----------------------------------

The matrix (after compression):

0.0 0.0 3.0 1.0

0.0 4.0 1.0 5.0

9.0 2.0 6.0 5.0

3.0 5.0 8.0 9.0

7.0 9.0 3.0 2.0

3.0 8.0 4.0 6.0

2.0 4.0 4.0 0.0

----------------------------------

The U matrix compressed:

9.0 2.0 6.0 5.0

3.0 5.0 8.0 9.0

7.0 9.0 3.0 2.0

-5.280423280423281 -1.2539682539682542 -0.6137566137566138 0.0

7.287575150300601 3.6513026052104207 6.0 0.0

2.9979375773408496 2.353361748934415 0.0 0.0

-0.4729407448174646 0.0 0.0 0.0

----------------------------------

The L matrix compressed:

0.4444444444444444 0.3333333333333333

0.11111111111111112 0.037037037037037056

0.3068783068783069 -0.36507936507936506

-0.13827655310621242 -0.5681362725450901

-0.010724597827581485 0.27443970851093086

0.2283067327095946 0.0

0.0 0.0

----------------------------------

The right vector:

1.0

2.0

3.0

4.0

5.0

6.0

7.0

----------------------------------

The solution vector:

-0.12296353762606659

1.3688906128782

0.22459270752521324

0.00426687354538399

-0.14041892940263748

1.905352986811482

-0.08514352211016335

----------------------------------

示例程序： 

package  com.nc4nr.chapter02.bandec;

public   class  BanDec  ... {

    double[][] a = ...{
            ...{ 3.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0 },
            ...{ 4.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0 },
            ...{ 9.0, 2.0, 6.0, 5.0, 0.0, 0.0, 0.0 },
            ...{ 0.0, 3.0, 5.0, 8.0, 9.0, 0.0, 0.0 },
            ...