我目前正在尝试使用scipy.optimize.leastsq将一些数据与python进行拟合。我想要拟合的数据是这样的:Sigmap是一个数值积分,它是H的函数,需要相当长的时间来计算。我不希望将积分作为拟合例程的一部分,否则积分将重复执行,并显着增加拟合例程的时间。因此,我想从其他地方查找积分值。我用来实现这个的代码是:
integral = np.loadtxt(text file of form: H
我试图用Sympy的符号积分来找到一个定积分的封闭形式。symbols('x s H', real=True, nonnegative=True)不幸的是,使用Python2.7.11注:我还使用Mathematica的符号积分功能来实现它,并提出了一个Gauss超几何函数。不幸的是,计算这个函数会返回一个复数,这在计算一个实积分时
我想用Python语言实现这些积分,其中每个积分变量都是list enter image description here enter image description here 这是我在执行之前写每个积分的内容的尝试br) )]) 我记下H方程的误差 unsupported operand type(s) for +: 'generator' and 'generator' 我的第二个问题是如何对一个数组进行闭环积分