问题描述
我有一个圆锥和一个环,正如你在这里看到的,我怎样才能找到这两个表面的交点?非常感谢,
这是绘制这些曲面的代码,
def Surfaces(r,T1,T2):
t=(32*np.pi)/180. #theta=constant
r1,f=np.mgrid[0:1.1:40j,0:2.*np.pi:40j]
x=r1*np.sin(t)*np.cos(f)
y=r1*np.sin(t)*np.sin(f)
z=r1*np.cos(t)
a = (20*np.pi)/180 # deg
t = np.transpose(np.array([x,y,z]),(1,2,0))
m = [[np.cos(a),np.sin(a)],[0,1,0],[-np.sin(a),np.cos(a)]]
x,z = np.transpose(np.dot(t,m),(2,1))
ax.plot_surface(x,z,alpha=0.3,color='b')
f,t= np.mgrid[0:2*np.pi:40j,(T1*np.pi)/180:(T2*np.pi)/180:40j]
x = r*np.cos(f)*np.sin(t)
y = r*np.sin(f)*np.sin(t)
z = r*np.cos(t)
ax.plot_surface(x,alpha=0.5,color='g')
for angle in range(0,90):
ax.view_init(0,angle)
plt.show()
Surfaces(1,35,37)
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)