通过代码求解联立方程

时间:2012-07-23 08:43:34

标签: algorithm math equation

这似乎是一个令人难以置信的简单而愚蠢的问题,但我发现的一切都太复杂了,我无法理解。

我有这两个非常基本的联立方程式:

X = 2x + 2z
Y = z - x

鉴于我知道X和Y,我将如何找到x和z?手动完成它非常容易,但我不知道如何在代码中完成。

5 个答案:

答案 0 :(得分:13)

  

这似乎是一个非常简单和愚蠢的问题

完全没有。这是一个非常好的问题,不幸的是它是一个复杂的答案。我们来解决

a * x + b * y = u
c * x + d * y = v

我坚持这里的2x2案例。更复杂的案例将要求您使用库。

首先要注意的是Cramer formulas不适合使用。计算行列式时

a * d - b * c

只要您a * d ~ b * c,就有catastrophic cancellation。这种情况很典型,你必须防范它。

简单/稳定性之间的最佳折衷是partial pivoting。假设|a| > |c|。那么系统相当于

a * c/a * x + bc/a * y = uc/a
      c * x +    d * y = v 

cx + bc/a * y = uc/a
cx +       dy = v  

现在,减去第一个到第二个产量

cx +       bc/a * y = uc/a
     (d - bc/a) * y = v - uc/a

现在很容易解决:y = (v - uc/a) / (d - bc/a)x = (uc/a - bc/a * y) / c。计算d - bc/aad - bc稳定,因为我们除以最大数字(它不是很明显,但它保持 - 以非常接近的系数进行计算,你会明白为什么它有效)

现在,如果|c| > |a|,您只需交换行并继续进行。

在代码中(请检查Python语法):

def solve(a, b, c, d, u, v):
    if abs(a) > abs(c):
         f = u * c / a
         g = b * c / a
         y = (v - f) / (d - g)
         return ((f - g * y) / c, y)
    else
         f = v * a / c
         g = d * a / c
         x = (u - f) / (b - g)
         return (x, (f - g * x) / a)

你可以使用完全旋转(要求你交换x和y,这样第一个除法总是最大的系数),但是编写起来比较麻烦,而且2x2情况几乎不需要。

对于n x n的情况,所有的旋转内容都封装在LU decomposition中,你应该使用一个库。

答案 1 :(得分:1)

@Alexandre,你错过了一个条件.. 这是最终的代码

void SolveLinearEquations (float a,float b,float c,float d,float u,float v, float &x, float &y)
{
float f;
float g;
if (abs(a) > abs(c))
{
     f = u * c / a;
     g = b * c / a;
     y = (v - f) / (d - g);
     if(c != 0)
        x = (f - g * y) / c;
     else
        x = (u - b * y)/a;
}
else
{
    f = v * a / c;
     g = d * a / c;
     x = (u - f) / (b - g);
     if (a != 0)
        y = (f - g * x) / a ;
     else
        y = (v - d * x)/c;
}
}

答案 2 :(得分:0)

(1)  ax + by = c    
(2)  dx + dy = f

(3)1*d     adx + bdy  = cd
(4)2*b     abx + bdy  = fb

(3)-(4)   adx - abx  = cd - fb

          x(ad-ab) = cd - fb


          x = (c*d - f*b)/(a*d-a*b) //use this equation for x



ax + by = c
     by = c - ax

     y = (c - a*x)/b    //use this equation for y

答案 3 :(得分:0)

以下功能可能对某些人有用:

function solve(s1,s2){    //only works if coefficients > 0

   str=s1 + " " +s2
   str=str.replace(/[^0123456789.-]/g, ' ') //eliminate letters
   str=str.replace( /\s\s+/g, ' ' )         //no double spaces


   var n=str.split(" ");            //put into an array
   var a=0,b=1,c=2,d=3,e=4,f=5      //see what were doing
   var x = ( n[c]*n[e] -n[b]*n[f])/(n[a]*n[e] - n[b]*n[d])
   var y= (n[c]-n[a]*x)/n[b]

   return({x:x, y:y})
}

使用:

result=solve("12x +  2y =32", "9x -5y=55")
alert (result.x+" ----- "+result.y)

答案 4 :(得分:0)

编织的算法:Ax + By = C,Dx + Ey = F

START

1)  PRINT "Enter A"
2)  INPUT A
3)  PRINT "Enter B"
4)  INPUT B
5)  PRINT "Enter C"
6)  INPUT C
7)  PRINT "Enter D"
8)  INPUT D
9)  PRINT "Enter E"
10) INPUT E
11) PRINT "Enter F"
21) INPUT F

22) detS <-- AE - BC
22) detX <-- CE - BF
24) detY <-- AF - CD

25) x <-- detX / detS
26) y <-- detY / detS

27) PRINT x, y

END