Question

我刚开始用python编程。我想为各种函数规范化numpy arange。到目前为止，我的代码是

from numpy import arange
from scipy.integrate import quad
import matplotlib.pyplot as plt

def function_A(x):
    return -(x-1)*(x-5)

def normalize(functionName, numpyRange):
    print "Inside normalize"

    min = numpyRange.min()
    max = numpyRange.max()

    print "functionName: ", functionName
    print "min: ", min
    print "max: ", max
    print "quad(functionName, min, max)[0]: ", quad(functionName, min, max)[0]
    print "functionName(numpyRange): ", functionName(numpyRange) 

    print "RIGHT BEFORE RETURN"

    return functionName(numpyRange) / quad(functionName, min, max)[0]

x = arange(1,5.01,0.01)
plt.plot(x, normalize(function_A, x))
plt.show()

print "This should equal unity: ", quad(normalize(function_A, x), x.min(), x.max(), args=(function_A, x))

运行代码时出现错误消息：

Traceback (most recent call last):
  File "forStackExchange3.py", line 33, in <module>
    print "This should equal unity: ", quad(normalize(function_A, x), x.min(), x.max(), args=(function_A, x))
  File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 247, in quad
    retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points)
  File "/usr/lib/python2.7/dist-packages/scipy/integrate/quadpack.py", line 312, in _quad
    return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit)
quadpack.error: First argument must be a callable function.

错误消息让我感到困惑，因为即使错误发生在行

上

print "This should equal unity: ", quad(normalize(function_A, x), x.min(), x.max(), args=(function_A, x))

我仍然能够“输入”我的规范化函数，正如我的多个打印语句所证明的那样。我试过尝试不同的方式来编写提供错误信息的行，但没有任何运气。任何帮助，将不胜感激。也许我需要重写我的normalize函数，以便我可以正确使用quad？另外，如果在我的特定情况下有一个更简单的方法来在python中规范化函数，请告诉我。

Answer 1

documentation for quad州：

func：function要集成的Python函数或方法。

但是，您传递的第一个参数是ndarray，它不可调用。换句话说，您没有传递要集成的函数，而是传递其他一些此类操作的结果。

那是技术上错误的原因。至于如何修复它，我必须说我不确定你在问什么。给定范围y，可以将其标准化。范围y 可能是将某些功能应用于x的结果，但除此之外还有其他功能。因此，您可以构建范围规范化器，并将其应用于函数的结果。

def function_A(x):
    return -(x-1)*(x-5)

x = arange(1,5.01,0.01)

def normalize(y):
    return y / sum(y)

>> sum(normalize(function_A(x)))
1.0

你也可以构建一个元函数，它接受一个函数并返回一个规范化结果的函数：

def make_normalized(fn):
    return lambda x: fn(x) / sum(fn(x))

因此，你可以定义

normalized_function_A = make_normalized(function_A)

并且只是为了验证：

>> sum(normalized_function_A(x))
1.0

规范化Numpy Arange

1 个答案: