2016-11-15 38 views
0

为了计算衍生物和其它表达式我用sympy包,并表示T = sy.Symbol('T')现在我已经计算出了正确的表达:蟒:改变符号变量并分配数值

E= -T**2*F_deriv_T(T,rho) 

其中

def F_deriv_rho(T,rho): 
    ret = 0 
    for n in range(5): 
     for m in range(4): 
      inner= c[n,m]*g_rho_deriv_rho_np*g_T_np 
     ret += inner 
    return ret 

,看起来像这样:

F_deriv_rho: [0.0 7.76971e-5*T 0.0001553942*T**2*rho T*(-5.14488e-5*log(rho) - 5.14488e-5)*log(T) + T*(1.22574e-5*log(rho)+1.22574e-5)*log(T) + T*(1.89488e-5*log(rho) + 1.89488e-5)*log(T) + T(2.29441e-5*log(rho) + 2.29441e-5)*log(T) + T*(7.49956e-5*log(rho) + 7.49956e-5)*log(T) T**2*(-0.0001028976*rho*log(rho) - 5.14488e-5*rho)*log(T) + T**2*(2.45148e-5*rho*log(rho) + 1.22574e-5*rho)*log(T) + T**2*(3.78976e-5*rho*log(rho) + 1.89488e-5*rho)*log(T) + T**2*(4.58882e-5*rho*log(rho) + 2.29441e-5*rho)*log(T) + T**2*(0.0001499912*rho*log(rho) + 7.49956e 5*rho)*log(T)]

与python我想将T(和rho)作为符号更改为值。我怎么能这样做?

所以,我想创建10个数字,如T_def = np.arange(2000, 10000, 800),并通过迭代我在数组中创建的10个值来交换我所有的sy.symbol( T )

感谢您的帮助

+0

您的示例代码不是自包含的,例如,你引用诸如''g_T_np''之类的东西使得它很难提供帮助。但是从我收集的内容来看,您可能会对“lambdify”感兴趣,它会根据符号表达式为数值评估创建一个回调。 –

+0

@BjoernDahlgren:我知道,但其余的代码对我的问题并不重要,可能会更混乱。我想要的只是改变符号T并将其替换为数字值。但我现在正在检查lamdify,它可能会非常有帮助!非常感谢你 – Shaun

回答