2

10,000個の変数があります。それらの100については、私は正確な値を知っています。

他のものは次のように与えられます：

a = 0.x_1 * b + 0.y_2 * c+ 0.z_1 * d + (1 - 0.x_1 - 0.y_1 - 0.z_1) * a
b = 0.x_2 * c + 0.y_2 * d+ 0.z_2 * e + (1 - 0.x_2 - 0.y_2 - 0.z_2) * b

...

q = 0.x_10000 * p + 0.y_10000 * r+ 0.z_10000 * s + (1 - 0.x_10000 - 0.y_10000 - 0.z_10000) * q

はい私は0.x_n、0.y_n、0.z_nの正確な値を知っています...（接頭辞として0を持っているという点は、それが1未満であり、0より大きいことを意味します）

Pythonでこれを解決するにはどうすればよいですか？このような簡単な方程式を使って、いくつかの例を教えていただければ幸いです。

x - y + 2z =  5
    y -  z = -1
         z =  3
4

1 に答える 1

10

( numpyを使用) 一次方程式系を書き直すと

x - y + 2z =  5
    y -  z = -1
         z =  3

行列方程式として

A x = b

と

A = np.array([[ 1, -1,  2],
              [ 0,  1, -1],
              [ 0,  0,  1]])

と

b = np.array([5, -1, 3])

次に、次xを使用して見つけることができますnp.linalg.solve。

import numpy as np

A = np.array([(1, -1, 2), (0, 1, -1), (0, 0, 1)])
b = np.array([5, -1, 3])
x = np.linalg.solve(A, b)

収量

print(x)
# [ 1.  2.  3.]

次のことを確認できA x = bます。

print(np.dot(A,x))
# [ 5. -1.  3.]

assert np.allclose(np.dot(A,x), b)
于 2012-11-23T04:52:15.627 に答える