2010-03-02 51 views
1

我试图解决在其结束是矩形结构域的泊松方程像 斧一个线性问题= B 但因为我知道的边界条件,有节点,其中我有解决方案的价值。我想我的问题是... 我怎样才能解决稀疏系统Ax = b的,如果我知道一些的x坐标是与未决定值取决于这些呢?除了我知道一些开始的解决方案之外,它与普通解决方案相同。求解稀疏线性问题的一些公知的边界值

谢谢!

回答

1

如果我理解正确,x的一些元素是已知的,有些则不是,你要解决Ax = b的为x的未知值,是否正确?

令组Ax = [A1 A2] [X1; x2] = b,其中向量x = [x1; x2],矢量x1具有x的未知值,而矢量x2具有x的已知值。然后,A1x1 = b-A2x2。因此,使用scipy.linalg.solve或任何其他所需的求解器求解x1。

+0

这可能是愚蠢的问,但我不太明白你的符号[A C] [X,Y] = B。你能解释一下吗?你是我的问题的重述听起来像它解决我的问题=] – 2010-03-02 06:16:58

+1

对不起,我用混乱的变量名。请参阅编辑。我将矢量x分成x1和x2,并将矩阵A划分为A1和A2。 – 2010-03-02 06:35:18

+0

如何将A分区为A1和A2? – 2010-03-02 17:14:15