この質問は、もともとhttps://github.com/RobotLocomotion/drake/issues/12484で提示されています。問題は関数の使い方です。
関数の説明: https://drake.mit.edu/pydrake/pydrake.forwarddiff.html#pydrake.forwarddiff.jacobian
私のコード: '''
from pydrake.forwarddiff import jacobian
import math
import numpy as np
def f(t):
return np.array([1.5 - 0.5 * sigmoid(t - 2.5), 0.5, 0.0])
print jacobian(f,0.)
''' 返されたエラー情報:
AttributeError Traceback (most recent call last)
in ()
48
49
---> 50 print jacobian(f,0.)
51 #print jacobian(f,np.array[0.])
52
/drake/lib/python2.7/site-packages/pydrake/forwarddiff.py in jacobian(function, x)
38 The function should be vector-input and vector-output.
39 """
---> 40 x_ad = np.empty(x.shape, dtype=np.object)
41 for i in range(x.size):
42 der = np.zeros(x.size)
AttributeError: 'float' object has no attribute 'shape'
If I changed that line into
print jacobian(f,np.array[0.])
then the error becomes:
The autoreload extension is already loaded. To reload it, use:
%reload_ext autoreload
TypeError Traceback (most recent call last)
in ()
49
50 #print jacobian(f,0.)
---> 51 print jacobian(f,np.array[0.])
52
53
TypeError: 'builtin_function_or_method' object has no attribute 'getitem'
If I changed that line into:
print jacobian(f,[0.])
the error:
AttributeError Traceback (most recent call last)
in ()
50 #print jacobian(f,0.)
51 #print jacobian(f,np.array[0.])
---> 52 print jacobian(f,[0.])
53
54 # Run the simulation. Parameters not described above
/drake/lib/python2.7/site-packages/pydrake/forwarddiff.py in jacobian(function, x)
38 The function should be vector-input and vector-output.
39 """
---> 40 x_ad = np.empty(x.shape, dtype=np.object)
41 for i in range(x.size):
42 der = np.zeros(x.size)
AttributeError: 'list' オブジェクトに属性 'shape' がありません
問題は、それを正しく使用するにはどうすればよいかということです。ありがとう
****************アップデートライン******************************** **********
エリックの発見に感謝します。コードを次のように変更します: '''
from pydrake.forwarddiff import jacobian
import math
import numpy as np
def f(t):
return np.array([1.5 - 0.5 *(t - 2.5)])
def g(t):
return np.array([1.5 - 0.5 *(t - 2.5), 0.5, 0.0])
x=0.
x = np.asarray(x)
print jacobian(f,x) # accepted
print jacobian(f,0.) # not accepted
print jacobian(g,x) # not accepted
''' 結果の出力: '''
[[-0.5]]
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
<ipython-input-15-1b94cb728eed> in <module>()
9 x = np.asarray(x)
10 print jacobian(f,x)
---> 11 print jacobian(f,0.)
/drake/lib/python2.7/site-packages/pydrake/forwarddiff.py in jacobian(function, x)
38 The function should be vector-input and vector-output.
39 """
---> 40 x_ad = np.empty(x.shape, dtype=np.object)
41 for i in range(x.size):
42 der = np.zeros(x.size)
AttributeError: 'float' object has no attribute 'shape'
''' したがって、ここには少なくとも 2 つの問題があります。これは厄介です。2. g(t) のヤコビアンに注意してください。その行のエラーは次のとおりです: ''' ----------------------------- ---------------------------------------------- AttributeError トレースバック (ほとんどのrecent call last) in () 11 print jacobian(f,x) # 受け付けた 12 #print jacobian(f,0.) # 受け付けなかった ---> 13 print jacobian(g,x) # 受け付けなかった
/drake/lib/python2.7/site-packages/pydrake/forwarddiff.py in jacobian(function, x)
45 y_ad = function(x_ad)
46 return np.vstack(
---> 47 [y.derivatives() for y in y_ad.flat]).reshape(y_ad.shape + (-1,))
48
49
AttributeError: 'float' object has no attribute 'derivatives'
''' 私の推測が正しければ、g(x) 内のすべてを np.asarray に変更する必要があります。これは厄介です...