0

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全体を受信して​​いることです。tcondlistでより小さい/より大きいステートメントを使用する場合、これは問題ではありませんが、を使用numpy.piecewiseして真か偽かを判断する方がはるかに簡単Q_trueです。

プロットは次のようになります。

プロット例

何か案は？

ありがとう！

4

1 に答える 1

2

次のバージョンのQ_true作品：

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

匿名関数に名前を付けていることに注意してください（Q1 = lambda ...）。この場合、を使用して関数を定義する必要がありますdef。

于 2011-03-22T16:30:37.747 に答える