2016-09-28 42 views
0

sympyを使用して部分偏導関数を実行しようとしていますが、これを関数に変換して値を代入して、 t_1、t_2の値。次のように私が使用していたコードは次のとおりです。Sympy to numpy AttributeError: 'Symbol'オブジェクトに 'cos'という属性がありません

import sympy as sp 
import numpy as np 
from sympy import init_printing 
init_printing() 
t_1,t_2,X_1,X_2,Y_1,Y_2,X_c1,X_c2,Y_c1,Y_c2,a_1,a_2,psi_1,psi_2,b_1,b_2= sp.symbols('t_1 t_2 X_1 X_2 Y_1 Y_2 X_c1 X_c2 Y_c1 Y_c2 a_1 a_2 psi_1 psi_2 b_1 b_2') 

X_1=X_c1 + (a_1 * sp.cos(t_1) * sp.cos(psi_1)) - ((b_1) * sp.sin(t_1)* sp.sin(psi_1)) 

X_2=X_c2 + (a_2 * sp.cos(t_2) * sp.cos(psi_2)) - ((b_2) * sp.sin(t_2)* sp.sin(psi_2)) 

Y_1=Y_c1 + (a_1 * sp.cos(t_1) * sp.sin(psi_1)) + ((b_1) * sp.sin(t_1)* sp.cos(psi_1)) 

Y_2=Y_c2 + (a_2 * sp.cos(t_2) * sp.sin(psi_2)) + ((b_2) * sp.sin(t_2)* sp.sin(psi_2)) 

D=(((X_2-X_1)**2) + ((Y_2-Y_1)**2))**0.5 

y_1=sp.diff(D,t_1) 

y_2=sp.diff(D,t_2) 

f=sp.lambdify(t_1, y_1, "numpy") 

g=sp.lambdify(t_2, y_2, "numpy") 

私が使用してT_1のために値を代入しようとすると、

f(np.pi/2) 

私は次のエラーを取得:

AttributeError Traceback (most recent call last) 
<ipython-input-26-f37892b21c8b> in <module>() 
----> 1 f(np.pi/2) 

/users/vishnu/anaconda3/lib/python3.5/site-packages/numpy /__init__.py in <lambda>(_Dummy_23) 

AttributeError: 'Symbol' object has no attribute 'cos' 

を私が言及しました次のリンク:

What causes this error (AttributeError: 'Mul' object has no attribute 'cos') in Python?

Python AttributeError:cos

が、私はnumpyのとsympyの私の輸入は、それらのリンクに説明した場合とは異なり、衝突されていないと思います。どんな助けもありがとうございます。

+2

式 'y_1'は、多くの変数の関数である:これはあなたの問題のために十分です。 'lambdify'関数を複数の入力引数を持つものとして定義するか、他のシンボリック変数を数値で置き換えてください。 – Stelios

+0

こんにちはStelios、返信ありがとうございます。他の変数を記号として取り扱う偏微分を計算し、t_1を代用することはできませんか?たとえば、t_1 = np.pi/2の場合、t_1に対する式 'cos(t_1)* X_c1'の偏微分は '-X_c1'に等しくなります。 – Vishnu

+0

「シンピー」で使用されている「シンボル」という概念は、「ナンシー」では無効です。 'Numpy'は(複数の)数値引数を持つ関数しか理解しません。 't_1'を数字だけに置き換えて、他の変数を記号として扱う必要があるなら、' sympy'を使うべきです(例えば 'subs'メソッドを使って) – Stelios

答えて

3

np.cos(a_symbol)を呼び出すと、このタイプのエラーが発生します。これは明らかにnumpyのフードをa_symbol.cos()に変換します。

lambdifyは、数値計算のためのもので、npコールですべてのspコールを置き換えます。しかし、あなたがしていることは象徴的です。しかし、あなたの `lambdify`定義がエ​​ラーになり、単一の入力引数(` t_1`)を意味

f1 = lambda t: y_1.subs({t_1: t}) 
f2 = lambda t: y_2.subs({t_2: t}) 
+0

ありがとうございます。これは質問に正しく答えます。あなたがf1(np.pi/2)を使うと、cos(t_1)を6.123 * 10 ** - 17に置き換えた理由が分かりますか?どのように代用ゼロにするのですか? – Vishnu

+0

代わりに 'sp.pi'を使用します – Eric

+0

完全に動作します。再度、感謝します。 :) – Vishnu

関連する問題