线性系统中矩阵的部分消除

时间:2017-06-06 21:05:39

标签: matlab

我想知道在matlab中消除(在线性系统求解意义上)矩阵的某些行的最有效方法是什么。基本上,给定行“s”和“l”(这样s和l是不同的,它们的联合是1:n),如何有效地计算

A是一个很大的n x n稀疏SPD矩阵。

R = chol(A(s,s))
Atemp = R' \ A(s,l)
All = A(l,l) - Atemp' * Atemp

基本上我正在计算部分高斯消元的尾随矩阵。 上面的代码,对于大而稀疏的A,几乎所有的时间花在第二行,特别是第三行。有关如何改善这一点的任何提示? s的大小通常是l的大小的10%。

非常感谢

0 个答案:

没有答案