Sympy to numpy вызывает AttributeError: объект «Символ» не имеет атрибута «cos»

Question

Я пытаюсь сделать частные производные с помощью sympy, и я хочу преобразовать его в функцию, чтобы я мог подставлять значения и оценивать производные на некотором значения t_1, t_2. Код, я использую выглядит следующим образом:

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") 

Когда я пытаюсь заменить значение Т_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 не сталкивается в отличие от случаев, упомянутых в этих ссылках. Любая помощь приветствуется.

Answer 1

Этот тип ошибки возникает, когда вы вызываете np.cos(a_symbol), что, очевидно, переводит под капот в numpy на a_symbol.cos().

lambdify предназначен для цифровых вычислений - он заменяет все звонки sp с звонками np. Но то, что вы делаете, символично. Этого достаточно для вашей проблемы:

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