Python-np.random.choice

Question

我正在使用 numpy.random.choice 模块根据函数数组生成选择“数组”：

def f(x):
    return np.sin(x)

def g(x):
    return np.cos(x)


base=[f, g]

funcs=np.random.choice(base,size=2)

此代码将生成一个包含 2 个项目的“数组”，引用基本数组中的一个函数。

这个帖子的原因是，我打印了funcs的结果并收到了：

[<function f at 0x00000225AC94F0D0> <function f at 0x00000225AC94F0D0>]

显然，这以某种形式返回对函数的引用，而不是我理解该形式是什么或如何操作它，这就是问题所在。我想更改函数的选择，以便不再是随机的，而是取决于某些条件，所以它可能是：

for i in range(2):
    if testvar=='true':
        choice[i] = 0 
    if testvar== 'false':
        choice[i] = 1

这将返回一个索引数组以放入以后的函数中

问题是，代码的进一步操作（我认为）需要以前形式的函数引用：[ ] 作为输入，而不是简单的 0,1 索引数组，我不知道该怎么做使用 if 语句获取表单 [ ] 的数组。

对于需要此输入的其余代码，我可能完全错了，但我不知道如何修改它，因此将其发布在这里。完整代码如下：（它是@Attack68 在Evolving functions in python 上提供的代码的轻微变体）它旨在存储一个函数，该函数在每次迭代时乘以随机函数并相应地积分。 （我已经对导致问题的函数上方的代码进行了注释）

import numpy as np
import scipy.integrate as int

def f(x):
    return np.sin(x)

def g(x):
    return np.cos(x)

base = [f, g]

funcs = np.random.choice(base, size=2)
print(funcs)
#The below function is where I believe the [<function...>] input to be required
def apply(x, funcs):
    y = 1
    for func in funcs:
        y *= func(x)
    return y

print('function value at 1.5 ', apply(1.5, funcs))

answer = int.quad(apply, 1, 2, args=(funcs,))
print('integration over [1,2]: ', answer)

这是我实现非随机事件的尝试：

import numpy as np
import scipy.integrate as int
import random
def f(x):
    return np.sin(x)

def g(x):
    return np.cos(x)

base = [f, g]

funcs = list()
for i in range(2):
    testvar=random.randint(0,100) #In my actual code, this would not be random but dependent on some other situation I have not accounted for here
    if testvar>50:
        func_idx = 0 # choose a np.random operation: 0=f, 1=g
    else:
        func_idx= 1
    funcs.append(func_idx)
#funcs = np.random.choice(base, size=10)
print(funcs)
def apply(x, funcs):
    y = 1
    for func in funcs:
        y *= func(x)
    return y

print('function value at 1.5 ', apply(1.5, funcs))

answer = int.quad(apply, 1, 2, args=(funcs,))
print('integration over [1,2]: ', answer)

这会返回以下错误：

TypeError: 'int' object is not callable

Answer 1

如果：您正在尝试将在随机选择的函数列表上运行的原始代码重构为使用与函数列表中的项目相对应的随机索引运行的版本。重构apply。

def apply(x,indices,base=base):
    y = 1
    for i in indices:
        f = base[i]
        y *= f(x)
    return y

---

...这以某种形式返回对函数的引用，而不是我了解该形式是什么或如何操作它...

函数是对象，列表包含对对象本身的引用。可以通过将它们分配给名称然后调用它们或索引列表并调用对象来使用它们：

>>> def f():
...     return 'f'

>>> def g():
...     return 'g'

>>> a = [f,g]
>>> q = a[0]
>>> q()
'f'
>>> a[1]()
'g'
>>> for thing in a:
    print(thing())

f
g

或者你可以传递它们：

>>> def h(thing):
...     return thing()

>>> h(a[1])
'g'
>>>

Answer 2

如果您仍想按原样使用函数apply，则需要将输入保留为函数列表。您可以使用这些索引来创建函数列表，而不是提供索引列表。

试试apply(1.5, funcs)，而不是：

apply(1.5, [base(n) for n in funcs])