1

FsolveinScipyがこれの正しい候補のようです。方程式を動的に渡すのに助けが必要です。事前にご意見いただければ幸いです。

動的とは、方程式の数が実行ごとに異なることを意味します。たとえば、ある状況では次のようになります。

alpha*x + (1-alpha)*x*y - y = 0
beta*x  + (1- beta)*x*z - z = 0
A*x + B*y + C*z = D

そして私が持っている別の状況:

alpha*x + (1-alpha)*x*y - y = 0
beta*x  + (1- beta)*x*z - z = 0
gama*x  + (1 -gama)*x*w - w =0
A*x + B*y + C*z + D*w = E

alpha、beta、A、B、C、DおよびEはすべて定数です。x、、、yは変数ですz。w

4

1 に答える 1

1

私は Fsolve を自分で使用したことはありませんが、そのドキュメントによると、呼び出し可能な関数が必要です。このようなものは、変数の数が不明な複数の関数を処理します。ここでは引数を正しく並べる必要がありますが、各関数は単純にリストを取ることに注意してください。

def f1(argList):
    x = argList[0]
    return x**2
def f2(argList):
    x = argList[0]
    y = argList[1]
    return (x+y)**2
def f3(argList):
    x = argList[0]
    return x/3

fs = [f1,f2,f3]
args = [3,5]
for f in fs:
    print f(args)

Fsolve の場合、次のようなことを試すことができます (未テスト):

def func1(argList, constList):
    x = argList[0]
    y = argList[1]
    alpha = constList[0]
    return alpha*x + (1-alpha)*x*y - y
def func2(argList, constList):
    x = argList[0]
    y = argList[1]
    z = argList[2]
    beta = constList[1]
    return beta*x  + (1- beta)*x*z - z
def func3(argList, constList):
    x = argList[0]
    w = argList[1] ## or, if you want to pass the exact same list to each function, make w argList[4]
    gamma = constList[2]
    return gama*x  + (1 -gama)*x*w - w
def func4(argList, constList):

    return A*x + B*y + C*z + D*w -E ## note that I moved E to the left hand side


functions = []
functions.append((func1, argList1, constList1, args01))
# args here can be tailored to fit your  function structure
# Just make sure to align it with the way you call your function:
# args = [argList, constLit]
# args0 can be null.
functions.append((func1, argList2, constList2, args02))
functions.append((func1, argList3, constList3, args03))
functions.append((func1, argList4, constList4, args04))

for func,argList, constList, args0 in functions: ## argList is the (Vector) variable you're solving for.
    Fsolve(func = func, x0 = ..., args = constList, ...)
于 2012-04-05T16:57:24.433 に答える