1

scipy で root-finder を使用するとエラーが発生するようです。誰かが私が間違っていることを指摘できるかどうか疑問に思っていました。
ルートを見つけている関数は簡単な例にすぎず、特に重要ではありません。

このコードを scipy 0.9.0 で実行すると:

import numpy as np
from scipy.optimize import fsolve

tmpFunc = lambda xIn: (xIn[0]-4)**2 + (xIn[1]-5)**2 + (xIn[2]-7)**3

x0 = [3,4,5]
xFinal = fsolve(tmpFunc, x0 )

print xFinal

次のエラー メッセージが表示されます。

Traceback (most recent call last):
  File "tmpStack.py", line 7, in <module>
    xFinal = fsolve(tmpFunc, x0 )
  File "/usr/lib/python2.7/dist-packages/scipy/optimize/minpack.py", line 115, in fsolve
    _check_func('fsolve', 'func', func, x0, args, n, (n,))
  File "/usr/lib/python2.7/dist-packages/scipy/optimize/minpack.py", line 26, in _check_func
    raise TypeError(msg)
TypeError: fsolve: there is a mismatch between the input and output shape of the 'func' argument '<lambda>'.
4

1 に答える 1

3

このルーチンを間違って使用しようとしていたようです。このルーチンには、私が与えた 3 つの変数を持つ 1 つの方程式と同じ数の方程式と変数が必要です。したがって、最小化する関数への入力が 3 次元配列である場合、出力は 3 次元配列でなければなりません。このコードは機能します：

import numpy as np
from scipy.optimize import fsolve

tmpFunc = lambda xIn: np.array( [(xIn[0]-4)**2 + xIn[1], (xIn[1]-5)**2 - xIn[2]) \
, (xIn[2]-7)**3 + xIn[0] ] )

x0 = [3,4,5]
xFinal = fsolve(tmpFunc, x0 )

print xFinal

これは、3 つの方程式を同時に解くことを表します。

于 2012-05-08T20:02:55.923 に答える