矩阵逆的模(dropping strategy)

2011-10-17

下面的文本主要用于算法设计和算法实现的需要。 主要内容参考了Golub的书籍。

适合任何三角矩阵（无论对角元是否为1，属于通用算法）

1

估计U1.1 1.2 2

1 ................................................................................................................................... 1

方法1： ........................................................................................................................ 1 方法2： ........................................................................................................................ 2

1估计L .................................................................................................................................... 2 2.1 2.2

方法1：estfullnrm_forward(T,D) ................................................................................ 2 方法2： ........................................................................................................................ 3

1 估计U

1见Golub的书：

1.1 方法1：

1.2 方法2：

2 估计L

1参考Golub的书，在估计U1的模算法基础上得到的对于下三角矩阵的算法：

2.1 方法1：

p(1:n)0

for k1:n

Choose d(k)

y(k)(d(k)p(k))/T(k,k)

p(k1:n)p(k1:n)y(k)T(k1:n,k)

end

2.2 方法2：

p(1:n)0

for k1:n

y(k)(1p(k))/T(k,k) y(k)(1p(k))/T(k,k)

p(k)p(k1:n)y(k)T(k1:n,k) p(k)p(k1:n)y(k)T(k1:n,k)

if |y(k)|||p(k)||1|y(k)|||p(k)||1

else end end

y(k)y(k) p(k1:n)p(k)

y(k)y(k) p(k1:n)p(k)