可编程地求解方程组？

Question

可能重复:
C++中的线性方程组？

我有以下2个方程组:

对于a,b,c,d:

0 = a * r1_x + b * r1_x * r1_y + c * r1_y + d
1 = a * r2_x + b * r2_x * r2_y + c * r2_y + d
0 = a * r3_x + b * r3_x * r3_y + c * r3_y + d
1 = a * r4_x + b * r4_x * r4_y + c * r4_y + d

对于e,f,g,h:

0 = e * r1_x + f * r1_x * r1_y + g * r1_y + h
0 = e * r2_x + f * r2_x * r2_y + g * r2_y + h
1 = e * r3_x + f * r3_x * r3_y + g * r3_y + h
1 = e * r4_x + f * r4_x * r4_y + g * r4_y + h

我知道r1_x,r1_y,r2_x,r2_y,r3_x,r3_y,r4_x,r4_y的值,需要求解第一个中的a,b,c,d,以及e,f,g,h中的值.

我知道如何用铅笔和纸来解决这些问题,但我真的不确定如何编程.我怎样才能用C或C++(或伪代码)解决上述方程式.

谢谢

Answer 1

您可以将其映射到矩阵系统,A x = b其中A系数矩阵b是解决方案向量,并且x是未知数.您可以实现高斯消除,也可以使用众所周知的库.如果你使用LAPACK,你想要它的例程dgesv.