正切函数的反三角函数,python 反正切函数
我试图复制柯蒂斯关于轨道力学的那集,但我不能完全理解它。不过我从np.arctan切换到np.arctan2还是有进步的。
我可能没有正确实现arctan2吗?进口实验室
进口编号为NP
e=NP.arange (0.0,1.0,0.15)(。整形(-1,1))))))))))))))).
nu=NP.linspace (0.001,2 * NP.pi-0.001,50000)))))))nu=NP . Lin space(0.001,2 * NP.pi-0.001,50000))
m2 evals=(2* NP.arctan2(1,1/)(1-e))0.5 * NP . tan)nu/2)-
e *(1-e * * 2)* 0.5 * NP . sin(nu)1e * NP . cos)nu))
图2=pylab.figure(
ax2=fig2.add_subplot(111))
对于Me2,_einzip(m2evals,e.ravel()):
ax2.plot(nu.ravel)、Me2、label=str ) _e))
pylab.legend(
pylab.xlim(),7.75))
py lab . ylim((0,2 * NP.pi)))))))))))))))))))).
pylab.show()).
下图有一处中断。这个函数是平滑的,应该在y范围(0,2pi)内连接0和2pi,而不是使0和2pi接触。
教科书故事和方程式:
cqdyf的要求,我说:
问题可能出在以“主值”为输出的反正切函数上。
所以,如果x是第二象限或第三象限的角度,Arctan(tan(x))不会生成x,如果Arctan(tan) x)从x=0到x=Pi作图,可以看出x=Pi/2时有不连续的跳跃。
在你的情况下,我想写arctan (1,1/arg)而不是arctan(arg)。其中arg是反正切函数的参数。因此,如果arg为负,arctan2将角度转向第二象限,而不是第四象限。"
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。