从 Python 中的文本文件导入方程答案

【问题标题】：Import equation from text file in Python从 Python 中的文本文件导入方程
【发布时间】：2018-02-24 02:08:38
【问题描述】：

我有一个复杂的方程式，它生成到一个.txt 文件中。我想导入这个方程（这是.txt 文件中的所有文本）并从中创建一个函数，然后可以拟合。

有人知道我会怎么做吗？要拟合的方程位于最底部。下面是我微弱的导入尝试...

myfile1= open("dummyfile.txt", 'r') 
def fcn(J1,J2,T,k,g):
    return myfile1.read()

“dummyfile.txt”内容：

B**2*N*(12.0*g**2*sp.exp(2.0*J2/(T*k)) + 60.0*g**2*sp.exp(6.0*J2/(T*k)) + 168.0*g**2*sp.exp(12.0*J2/(T*k)) + 360.0*g**2*sp.exp(20.0*J2/(T*k)) + 30.0*g**2*sp.exp((2.0*J1 + 4.0*J2)/(T*k)) + 168.0*g**2*sp.exp((4.0*J1 + 8.0*J2)/(T*k)) + 360.0*g**2*sp.exp((6.0*J1 + 14.0*J2)/(T*k)) + 180.0*g**2*sp.exp((8.0*J1 + 12.0*J2)/(T*k)) + 660.0*g**2*sp.exp((8.0*J1 + 22.0*J2)/(T*k)) + 660.0*g**2*sp.exp((12.0*J1 + 18.0*J2)/(T*k)) + 1092.0*g**2*sp.exp((16.0*J1 + 26.0*J2)/(T*k)) + 546.0*g**2*sp.exp((18.0*J1 + 24.0*J2)/(T*k)) + 1680.0*g**2*sp.exp((24.0*J1 + 32.0*J2)/(T*k)) + 1224.0*g**2*sp.exp((32.0*J1 + 40.0*J2)/(T*k)))/(3*T*k*(6*sp.exp(2.0*J2/(T*k)) + 10*sp.exp(6.0*J2/(T*k)) + 14*sp.exp(12.0*J2/(T*k)) + 18*sp.exp(20.0*J2/(T*k)) + 5*sp.exp((2.0*J1 + 4.0*J2)/(T*k)) + 14*sp.exp((4.0*J1 + 8.0*J2)/(T*k)) + 18*sp.exp((6.0*J1 + 14.0*J2)/(T*k)) + 9*sp.exp((8.0*J1 + 12.0*J2)/(T*k)) + 22*sp.exp((8.0*J1 + 22.0*J2)/(T*k)) + 22*sp.exp((12.0*J1 + 18.0*J2)/(T*k)) + 26*sp.exp((16.0*J1 + 26.0*J2)/(T*k)) + 13*sp.exp((18.0*J1 + 24.0*J2)/(T*k)) + 30*sp.exp((24.0*J1 + 32.0*J2)/(T*k)) + 17*sp.exp((32.0*J1 + 40.0*J2)/(T*k)) + 1))

【问题讨论】：

你应该看看eval，然后你应该重新考虑这个。
为了保持理智，请不要使用eval。但 SymPy 值得探索：docs.sympy.org/latest/modules/parsing.html
如果您已经可以将其作为纯文本读取，则无需在程序中“读取”它；只需将其复制到源代码中！提供语法正确，即。
你好。仅复制和粘贴文本当然可以，但是方程式是在上面生成的，并且会根据我告诉程序有关我的系统的内容而有所不同。因此，跳过这一步并让程序直接执行会更有吸引力

标签： python function import

【解决方案1】：

您可以使用exec() 做到这一点。

代码：

def build_function(filename):
    with open(filename, 'rU') as f:
        eqn = f.read().strip()
        exec("def fcn(J1, J2, T, k, g):\n return ({})".format(eqn))
        return locals()['fcn']

测试代码：

fcn = build_function('file1')
print(fcn(1, 2, 3, 4, 5))

文件1：

J2 + T*k

结果：

【讨论】：