Sympy到numpy會導致AttributeError：'Symbol'對象沒有屬性'cos'

我想使用sympy做偏導數，我想將它轉換爲一個函數，以便我可以替換值並估計一些導數t_1，t_2的值。我使用的代碼如下：Sympy到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不會碰撞不像這些鏈接中提到的情況。任何幫助表示讚賞。

來源

2016-09-28 Vishnu

表達'y_1'是許多變量的函數，但是，你的'lambdify'定義意味着只有一個輸入參數（'t_1'），導致錯誤。要麼將'lambdify'函數定義爲具有多個輸入參數，要麼用數字替換其他符號變量。 – Stelios

嗨Stelios，謝謝你的回覆。我們不可能計算偏導數作爲符號處理其他變量，只能替代t_1？例如：如果t_1 = np.pi/2，則表達式'cos（t_1）* X_c1'相對於t_1的偏導數等於'-X_c1'。 – Vishnu

「Sympy」中使用的「符號」的概念在「Numpy」中無效。 'Numpy'只能理解帶有（多個）數值參數的函數。如果你只需要用數字替換't_1'，將其他變量作爲符號來處理，就應該使用'Sympy'（例如，使用'subs'方法） – Stelios

當您撥打np.cos(a_symbol)時，會發生此類型的錯誤，該錯誤明顯可將numpy中的引擎翻譯爲a_symbol.cos()。

lambdify用於數字計算 - 它將用np調用替換所有sp調用。但是你在做什麼是象徵性的。這是足以讓你的問題：

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

來源

2016-09-28 17:21:13 Eric

謝謝，這回答正確的問題。你知道爲什麼如果我使用f1（np.pi/2）它將cos（t_1）替換爲6.123 * 10 ** - 17。我如何讓它替代零？ – Vishnu

您可以使用'sp.pi'而不是 – Eric

完美地工作。再次感謝。 :) – Vishnu

Sympy到numpy會導致AttributeError：'Symbol'對象沒有屬性'cos'

回答

相關問題