Question

我从数据库中获得了一些模型

f(t)=(2.128795454425367)+(208.54359721863273)*t+(26.098128487929266)*t^2+(3.34369909584111)*t^3+(-0.3450228278737971)*t^4+(-0.018630757967458885)*t^5+(0.0015029038553239819)*t^6;

其中有些是字符串。

我需要为t in range(1, 13)

评估此功能

现在我必须手动复制这些功能并运行它们

   print [1.2381648958643592 + \
          153.55656654019816 * t +\
          22.99318731025164 * (t**2) +\
          11.060577906796075 * (t**3) +\
          -1.3465054084767891 * (t**4) + \
          0.016926765998876842 * (t**5) +\
          0.001500086893490721 * (t**6) for t in range(1, 13)]

在python中有更好的方法吗？

Answer 1

如果表现不是主要问题 - 如果你只是在12分进行评估，我怀疑它不是 - 那么你可以利用方便的sympy库做很多工作为了你。例如：

>>> import sympy
>>> sympy.sympify("t**5 - t + 3")
t**5 - t + 3
>>> sympy.sympify("t**5 - t + 3").subs({"t": 10})
99993

我们可以将它包装在一个返回函数的函数中：

import sympy

def definition_to_function(s):
    lhs, rhs = s.split("=", 1)
    rhs = rhs.rstrip('; ')
    args = sympy.sympify(lhs).args
    f = sympy.sympify(rhs)
    def f_func(*passed_args):
        argdict = dict(zip(args, passed_args))
        result = f.subs(argdict)
        return float(result)
    return f_func

然后我们可以应用，甚至更容易达到正则表达式的更复杂的情况：

>>> s = "f(t)=(2.128795454425367)+(208.54359721863273)*t+(26.098128487929266)*t^2+(3.34369909584111)*t^3+(-0.3450228278737971)*t^4+(-0.018630757967458885)*t^5+(0.0015029038553239819)*t^6;"
>>> f = definition_to_function(s)
>>> f(0)
2.128795454425367
>>> f(10)
4230.6764921149115
>>> f = definition_to_function("f(a,b,c) = sin(a)+3*b-4*c")
>>> f(1,2,3)
-5.158529015192103
>>> import math
>>> math.sin(1)+3*2-4*3
-5.158529015192103

Answer 2

如果你想解析'function'字符串，你可以这样做：

import re

s = "f(t)=(2.128795454425367)+(208.54359721863273)*t+(26.098128487929266)*t^2\
    +(3.34369909584111)*t^3+(-0.3450228278737971)*t^4+(-0.018630757967458885)*t^5\
    +(0.0015029038553239819)*t^6;"

def f(t):
    l = map(float, re.findall("-?\\d+\\.\\d+", s))
    return sum(b * t**a for a,b in enumerate(l))

print map(f, xrange(1,13))

[239.75206957484252, 544.337732955938, 921.544112756058, 1366.6221363666925, 1864.8848673959649, 2393.2591324279497, 2922.9192385578326, 3423.0027817028927, 3865.4085456893295, 4230.676492114911, 4514.949840987468, 4738.019242139209]

此方法假定函数字符串的格式始终为

c0 + c1 t + c2 t^2 + c3 t^4 + ... cn t^(n+1)

并通过从字符串中提取浮点数并使用它们来生成实际的Python函数来工作。

Answer 3

您可以将函数存储为数据库中的python expersion，当您获取字符串时，只需执行eval（funcstr.replace（'x'，'yvalue'））。

为您举例：

funcstr = '2*x+5'
evalpoint = funcstr.replace('x', '5')
val = eval(funcstr)

此时应将val评估为15

Answer 4

正如NPE所说，这里正确的答案是为你的表达语言编写一个解析器（和简单的解释器）。

或者，更好的是，如果可能的话，首先在Python中生成表达式，而不是使用几乎与Python的子集完全兼容的语言。

或者，更好的是，如果语言只是表示多项式系数列表的一种方式，只需将其表示为系数列表，这将比任何实际的通用语言更容易解析。例如，假设数据库保留了这个：

2.128795454425367, 208.54359721863273, 26.098128487929266, 3.34369909584111, -0.3450228278737971, -0.018630757967458885, 0.0015029038553239819

然后，要在Python中执行它，你会这样做：

def eval_polynomial(polynomial, value):
    coefficients = [float(x.strip()) for x in polynomial.split(',')]
    return sum(coefficient * (value**exponent) 
               for exponent, coefficient in enumerate(coefficients))

然后：

>>> [eval_polynomial(expr, t) for t in range(1, 13)]

但是，如果你真的，真的想在不改变数据库中的内容的情况下这样做，你可以将其转换为Python表达式并对其进行评估：

>>> expr = 'f(t)=(2.128795454425367)+(208.54359721863273)*t+(26.098128487929266)*t^2+(3.34369909584111)*t^3+(-0.3450228278737971)*t^4+(-0.018630757967458885)*t^5+(0.0015029038553239819)*t^6;'
>>> removef = re.sub(r'f\((\w+)\)=', 'lambda \1: ', expr)
>>> fixpower = re.sub(r'(\w+)\^(\d+)', r'(\1**\2)', removef)
>>> nosemi = fixpower.replace(';', '')
>>> func = eval(nosemi)
>>> [func(t) for t in range(1, 13)]
[239.75206957484252, 544.337732955938, 921.544112756058, 1366.6221363666925, 1864.8848673959649, 2393.2591324279497, 2922.9192385578326, 3423.0027817028927, 3865.4085456893295, 4230.676492114911, 4514.949840987468, 4738.019242139209]

但是，你可能不想这样做。

而且，如果你这样做，你可能想要编写一个适用于你的实际语言的变换器，而不是根据一个例子对你的语言进行暗中猜测......

Answer 5

如果您信任您的来源，您可以使用正则表达式和eval执行此操作：

# deletes the simicolon and everything before the space
my_str = start_str.split('=')[1][:-1]
# change ^ to ** because that's the squared operator
my_str = re.sub('\^', '**', my_str)
# substitute the t for the numbers 1 to 13 and evaluate the string
results = [eval(re.sub('t', str(t), my_str)) for t in range(1,13)]

Python：如何评估一个字符串函数？

5 个答案: