小编tco*_*sar的帖子

我知道这被问了很多次，但是，我仍然遇到以下问题。我为 mean 和 stdev 定义了自己的函数，但是 stdev 计算 std(Wapproxlist) 花费的时间太长。所以，我需要解决这个问题。

import numpy as np
def Taylor_Integration(a, b, mu):
    import sympy as sy
    A, B, rho = sy.symbols('A B rho', real=True)
    Wapp = (A + B*rho - rho/(2*mu*(1 - rho)))**2
    eq1 = sy.diff(sy.integrate(Wapp, (rho, a, b)),A)
    eq2 = sy.diff(sy.integrate(Wapp, (rho, a, b)),B)
    sol = sy.solve([eq1,eq2], [A,B])
    return sol[A], sol[B]

def Wapprox(rho, A, B):
    return A + B*rho

def W(mu, rho):
    return rho/(2*mu*(1-rho))

Wapproxlist = []
Wlist = []
alist = np.linspace(0, 0.98, …