我正在尝试寻找光线和3D曲面之间的交点,以球面坐标r=f(θ)给出。
由于不确定如何简单地使用Python,我首先将曲面方程转换为笛卡尔方程。我用渐近模,来求解关于t的二阶参数方程。
import numpy as np
from sympy.solvers import solve
from sympy import Symbol
def mpn(Pd,Bz,X,Y,Z,px,py,pz):
r0=(10.22+1.29*np.tanh(0.184*(Bz+8.14)))*Pd**(-1/6.6)
alpha=(0.58-0.007*Bz)*(1+0.024*np.log(Pd))
t=Symbol('t',positive="True")
x=X+px*t
y=Y+py*t
z=Z+pz*t
sol=solve(x**2+y**2+z**2-r0**2*(2/(1+x/np.sqrt(x**2+y**2+z**2)))**(2*alpha),t)[0]
return (X+px*sol,Y+py*sol,Z+pz*sol)
sol=mpn(3.69,-4.20,-6.71,-10.3,-5.05,0.0582,0.0012,-0.9983)
但这会产生一个奇怪的错误,'Add‘对象没有'sqrt’属性。我不知道为什么会这样,因为我用的是np.sqrt,而不是Add.sqrt或+.sqrt之类的东西。有人能帮我吗?
发布于 2018-07-10 08:53:20
您应该使用sympy.sqrt()而不是numpy.sqrt()。
https://stackoverflow.com/questions/51255742
复制相似问题