为了计算衍生物和其它表达式我用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 )
。
感谢您的帮助
您的示例代码不是自包含的,例如,你引用诸如''g_T_np''之类的东西使得它很难提供帮助。但是从我收集的内容来看,您可能会对“lambdify”感兴趣,它会根据符号表达式为数值评估创建一个回调。 –
@BjoernDahlgren:我知道,但其余的代码对我的问题并不重要,可能会更混乱。我想要的只是改变符号T并将其替换为数字值。但我现在正在检查lamdify,它可能会非常有帮助!非常感谢你 – Shaun