1

さまざまな動作を確認するために変更できる2つのパラメーターを持つ関数をグラフ化しようとしています。スライダーを使用してパラメーターを変更したいと思います。

私の検索では、軸を変更するが数学関数の一部を変更しないスライダーに出くわしました。

したがって、2 つのパラメーター Gmax と Km が軸である場合に機能する次のコードがあります。

    from matplotlib.widgets import Slider
    import numpy as np

    Gmax=1
    Km= 1

    def f(S):
        s1 = Gmax*S  #G_max
        e1 = S + Km #K_m
        return divide(s1,e1)

    S=arange(0,100,0.1)

    ax = subplot(111)
    subplots_adjust(left=0.15, bottom=0.25)
    l = plot(f(S))
    grid(False)
    title('Playing with sliders')
    xlabel('time')
    ylabel('concentration')


    axcolor = 'lightgoldenrodyellow'
    axGmax = axes([0.15, 0.1, 0.65, 0.03], axisbg=axcolor)
    axKm = axes([0.15, 0.15, 0.65, 0.03], axisbg=axcolor)

    sGmax = Slider(axGmax, 'Gmax', 0.1, 3.0, valinit=1)
    sKm = Slider(axKm, 'Km', 0.01, 1.0, valinit=1)

    def update(val):
        s1 = Gmax*S * sGmax.val 
        e1 = S + Km * sKm.val
        l.set_ydata(y)    
        ax.set_ylim(y.min(), y.max())  
        draw()

    sGmax.on_changed(update)
    sKm.on_changed(update)

    show()

だから私の質問は、軸スライダーの ax コマンドの代わりにパラメーターのコマンドがあるかどうかだと思いますか? または、それを行う別の方法がある場合は？

4

1 に答える 1

3

あなたのコードはほぼ正しいですが、 plot() がリストを返すため、に変更l = plot(f(S))する必要があります。l, = plot(f(S))次に、呼び出しl.set_ydata(...)て新しい値を設定できます。

コードは次のとおりです。

import matplotlib.pyplot as plt
import pylab
import numpy as np


def f(S, Gmax, Km):
    s1 = Gmax*S   # G_max
    e1 = S + Km  # K_m
    return np.divide(s1, e1)


def update(val):
    l.set_ydata(f(S, sGmax.val, sKm.val))


S = np.arange(0, 100, 0.1)

ax = plt.subplot(111)
plt.subplots_adjust(left=0.15, bottom=0.25)
l, = plt.plot(f(S, 1.0, 1.0))
plt.grid(False)
plt.title('Playing with sliders')
plt.xlabel('time')
plt.ylabel('concentration')

axcolor = 'lightgoldenrodyellow'
axGmax = plt.axes([0.15, 0.1, 0.65, 0.03], facecolor=axcolor)
axKm = plt.axes([0.15, 0.15, 0.65, 0.03], facecolor=axcolor)

sGmax = pylab.Slider(axGmax, 'Gmax', 0.1, 3.0, valinit=1)
sKm = pylab.Slider(axKm, 'Km', 0.01, 1.0, valinit=1)

sGmax.on_changed(update)
sKm.on_changed(update)

plt.show()
于 2012-08-19T05:24:01.033 に答える