首先,制作一个100 x100矩阵稀疏是愚蠢的。平的,愚蠢的。系统需要解决的时间很小,即使一个完整的矩阵。
不再是一个单一的你已经形成了吗?
气孔导度(全(A))
ans =
1043.6
甚至没有接近奇异。这意味着任何列排列没有影响的解决方案。
时间(@ ()\ B)
ans =
7.2113 e-05
几毫秒才解决。完整的形式,它是不同的吗?
房颤=全(A);
时间(@()房颤\ B)
ans =
7.734 e-05
事实上,没有。所以除非你真正的问题是非常大,那么你就是在浪费你的时间与整个努力使稀疏,或者就此而言,使用列排列减少填写。即使你做了这个1000 x1000问题,直接全部解决仍会快得多。
它将帮助如果你告诉我们你真正的问题是,多大或者为什么你迫切想要解决这个问题,你这样做。
即使没有,显然你仍然有一个问题在如何使用permutaiton列。
b = b(:)的意思(b (:));
陆[L U] =(一个);x1 = U \ (L \ b);
规范(x1 - A \ b)
ans =
0
x1计算正确。现在让我看看你是如何计算x2,使用列排列。
一个列排列是什么意思?
想象我们想解决这个问题一个* x = b。为此,我们将一个列的排列A。但这并不涉及交换b的行! ! ! ! ! ! !
的permutaiton相当于重新排序的列的元素x这个问题。我们预计,如果x是一个解决方案* x = b,和p是一个列排列的,然后我们会
* x = = (: p) * x (p)
不需要交换b的元素,也不是你想这样做。
x2 = 0(大小(x1));
p = colamd(一个);
陆[L U] = ((:, p));%这个构建排列成L和U。
x2 (p) = U \ (L \ b);
规范(x1, x2)
ans =
1.1414 e-13
工作的很好。的细微差别是气孔导度(A) *每股收益。
电导率(完整的(A)) *每股收益
ans =
2.3173 e-13
这就是为什么你有一个问题。由于某些奇怪的原因,你决定你需要交换的元素b。你还浪费时间与排列在第一时间在这个小问题,但也许你真正的问题是严重的,大规模的。谁知道呢?
