使用 numpy.piecewise 生成分段周期图时的条件检查

Question

我正在尝试使用 Numpy 和 matplotlib 生成分段周期性图，如下所示：

import numpy as np
import matplotlib.pyplot as plt

Q1 = lambda t, f, Q_max: Q_max * np.sin(2 * np.pi *f * t)
Q2 = lambda t, f, Q_max: 0

def Q_true(t, f, stat):
    while(t >= 1/f):
        t -= 1/f 
    while(t < 0): 
        t += 1/f 
    return ((t <= 1/(2*f)) == stat)

Q = lambda t, f, Q_max: np.piecewise(t, [Q_true(t, f, True) , Q_true(t,f, False)], [Q1, Q2], f, Q_max)

Q_max = 225 # mL/sec
f = 1.25 # Hz
t = np.linspace(0,4,101) # secs
plt.plot(t, Q(t, f, Q_max))

问题是 Q_true 正在接收整个 t 数组而不是单个点。如果我只是在 numpy.piecewise condlist 中使用小于/大于语句，这不是问题，但使用 Q_true 确定它是真还是假要容易得多。

情节应如下所示：

有什么想法吗？

谢谢！

Answer 1

以下版本的 Q_true 有效：

def Q_true(t, f, stat):
    period = 1/f
    return (t % period < period/2) == stat

请注意，您正在命名匿名函数 (Q1 = lambda ...)。在这种情况下，您应该只使用 def 定义函数。