I am trying to do partial derivatives using sympy and I want to convert it to a function so that I can substitute values and estimate the derivatives at some values of t_1, t_2. The code I am using is as follows:
我正在尝试使用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")
When I try to substitute a value for t_1 using,
当我尝试用t_1替换值时,
f(np.pi/2)
I get the following error:
我收到以下错误:
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'
I referred to the following links:
我提到了以下链接:
What causes this error (AttributeError: 'Mul' object has no attribute 'cos') in Python?
导致此错误的原因(AttributeError:'Mul'对象在Python中没有属性'cos')?
Python AttributeError:cos
but I think my imports of numpy and sympy are not *ing unlike the cases mentioned in those links. Any help is appreciated.
但我认为,与那些链接中提到的情况不同,我的numpy和sympy的输入并没有发生冲突。任何帮助表示赞赏。
1 个解决方案
#1
3
This type of error occurs when you call np.cos(a_symbol)
, which apparently translates under-the-hood in numpy to a_symbol.cos()
.
当您调用np.cos(a_symbol)时会发生这种类型的错误,这显然会将numpy中的引擎转换为a_symbol.cos()。
lambdify
is for numeric calculations - it replaces all sp
calls with np
calls. But what you're doing is symbolic. This is enough for your problem:
lambdify用于数值计算 - 它用np调用替换所有sp调用。但你所做的是象征性的。这足以解决您的问题:
f1 = lambda t: y_1.subs({t_1: t})
f2 = lambda t: y_2.subs({t_2: t})
#1
3
This type of error occurs when you call np.cos(a_symbol)
, which apparently translates under-the-hood in numpy to a_symbol.cos()
.
当您调用np.cos(a_symbol)时会发生这种类型的错误,这显然会将numpy中的引擎转换为a_symbol.cos()。
lambdify
is for numeric calculations - it replaces all sp
calls with np
calls. But what you're doing is symbolic. This is enough for your problem:
lambdify用于数值计算 - 它用np调用替换所有sp调用。但你所做的是象征性的。这足以解决您的问题:
f1 = lambda t: y_1.subs({t_1: t})
f2 = lambda t: y_2.subs({t_2: t})