Python 为什么scipy中的mathieu函数不是周期性的?
我在scipy.special中发现Mathieu函数不是周期性的。以下是代码和绘图:Python 为什么scipy中的mathieu函数不是周期性的?,python,numpy,scipy,wolfram-mathematica,gsl,Python,Numpy,Scipy,Wolfram Mathematica,Gsl,我在scipy.special中发现Mathieu函数不是周期性的。以下是代码和绘图: import scipy.special as sp import matplotlib.pyplot as plt import numpy as np x = np.linspace(0,10*np.pi,100) plt.plot(x,sp.mathieu_cem(3,3,x)[0]) plt.show() 然而,在Mathematica中,它显然是周期性的: Plot[MathieuC[Mathi
import scipy.special as sp
import matplotlib.pyplot as plt
import numpy as np
x = np.linspace(0,10*np.pi,100)
plt.plot(x,sp.mathieu_cem(3,3,x)[0])
plt.show()
然而,在Mathematica中,它显然是周期性的:
Plot[MathieuC[MathieuCharacteristicA[3, 3], 3, x], {x, 0, 10 Pi}]
这两个函数不一样吗?由于scipy的文档没有太多细节,我不知道如何获得ce
的定期版本
顺便问一下,在gsl中有,与scipy有什么关系?他们是同一个吗
编辑
转换成学位后,我发现scipy中的循环比mathematica中的循环多,有什么解释吗
scipy.special.mathieu\u cem
以度为单位获取x
。Python代码上升到31.4度,Mathematica代码上升到31.4弧度时,此输出看起来差不多正确。scipy.special.mathieu_cem
以度为单位取x。此输出看起来与Python代码的角度上升到31.4度和Mathematica代码的角度上升到31.4弧度差不多。SciPy的SciPy.special.mathieu\u cem
采用x
度。re:edit。您没有正确地转换为度Cipy的scipy.special.mathieu_cem
以度为单位获取x
。re:edit。未正确转换为度要将弧度转换为度,请乘以(180/π)。在问题中更新的代码中,您没有除以π。要将弧度转换为度,请乘以(180/π)。在问题中更新的代码中,没有除以π。