-1

point から始まる一連のNランダム パスがありますX。初期ポイントを変更して新しいパスXを生成する Slider があります。Nパスを更新できるようにしたいのですが、set_ydata1D 配列しか受け入れられません。現時点では、更新ごとに軸をクリアしてプロットしていますが、これはあまり効率的ではありません。でそれを行う組み込みの方法はありmatplotlibますか？

xJ = arange(-10,10,0.1)
psinaive = zeros((xJ.shape[0]))
uapprox = zeros((xJ.shape[0],Nt))

wplot = []
wcondplot = []

for i,x in enumerate(xJ):
    WJ = sqrt(dt)*np.random.randn(Ntraj,Nt)
    WJ[:,0] = x
    WJ = np.cumsum(WJ,1)
    wplot.append(WJ)
    cond = V(WJ,limits)[0]
    wcondplot.append(WJ[cond,:])
    wa = 1.0/WJ.shape[0]*exp(-phi(WJ[:,-1],alpha)/lmbda)
    psinaive[i] = sum(wa[cond])
    uapprox[i,:] = 1.0/psinaive[i]*np.dot(wa[cond],WJ[cond,:]).flatten()
    if i % 10 == 0:
        print '..%.1f'%x,



J = -lmbda*log(psinaive)

#plot results

ax1=subplot(221)
plot(xJ,J)
plot(xJ,Jl(xJ,ti,tf,alpha,R,v,t1,limits))
title('$J(x,t)$')
plt.axvline(x=-10)

subplot(222)
plot(xJ,uapprox[:,0])

ax3 = subplot(223)
plot(t,wplot[0].T,alpha=0.2)
title('%d paths'%Ntraj)
ylim((-15,15))
ymin,ymax = ylim()
plt.axvline(x=t1, ymin=(d-ymin)/(ymax-ymin), linewidth=2, color='k')
plt.axvline(x=t1, ymin=(b-ymin)/(ymax-ymin), ymax = (c-ymin)/(ymax-ymin), linewidth=2, color='k')
plt.axvline(x=t1, ymax=(a-ymin)/(ymax-ymin), linewidth=2, color='k')

ax4 = subplot(224)
ax4.set_title('No alive paths')
if len(wcondplot[0])>0:
    ax4.plot(t,wcondplot[0].T,alpha=0.2)
    ax4.set_title('%d alive from %d paths'%(len(wcondplot[0]),Ntraj))
ylim((-15,15))
ymin,ymax = ylim()
plt.axvline(x=t1, ymin=(d-ymin)/(ymax-ymin), linewidth=2, color='k')
plt.axvline(x=t1, ymin=(b-ymin)/(ymax-ymin), ymax = (c-ymin)/(ymax-ymin), linewidth=2, color='k')
plt.axvline(x=t1, ymax=(a-ymin)/(ymax-ymin), linewidth=2, color='k')
subplots_adjust(0.15,0.25)


axsx = axes([0.15,0.1,0.75,0.1])
slx = Slider(axsx,'x',0,len(xJ),0,valfmt='%.0f')



def updatex(val):
    x = int(val)

    ax1.cla()
    ax1.plot(xJ,J)
    ax1.plot(xJ,Jl(xJ,ti,tf,alpha,R,v,t1,limits))
    ax1.set_title('$J(x,t)$')
    ax1.axvline(x=xJ[x])

    ax3.cla()
    ax3.plot(t,wplot[x].T,alpha=0.2)
    ax3.set_title('%d paths'%Ntraj)
    ax3.set_ylim((-15,15))
    ymin,ymax = ax3.get_ylim()
    ax3.axvline(x=t1, ymin=(d-ymin)/(ymax-ymin), linewidth=2, color='k')
    ax3.axvline(x=t1, ymin=(b-ymin)/(ymax-ymin), ymax = (c-ymin)/(ymax-ymin), linewidth=2, color='k')
    ax3.axvline(x=t1, ymax=(a-ymin)/(ymax-ymin), linewidth=2, color='k')

    ax4.cla()
    ax4.set_title('No alive paths')
    if len(wcondplot[x])>0:
        ax4.plot(t,wcondplot[x].T,alpha=0.4)
        ax4.set_title('%d alive from %d paths'%(len(wcondplot[x]),Ntraj))
    ax4.set_ylim((-15,15))
    ymin,ymax = ax4.get_ylim()
    ax4.axvline(x=t1, ymin=(d-ymin)/(ymax-ymin), linewidth=2, color='k')
    ax4.axvline(x=t1, ymin=(b-ymin)/(ymax-ymin), ymax = (c-ymin)/(ymax-ymin), linewidth=2, color='k')
    ax4.axvline(x=t1, ymax=(a-ymin)/(ymax-ymin), linewidth=2, color='k')
    draw()

slx.on_changed(updatex)

結果は次のとおりです。

結果

4

2 に答える 2

1

配列に nans が含まれている場合、配列の 1d-mensionality にもかかわらず、matplotlib は個別のセグメントを作成するという事実を使用できます。

x = (np.arange(5 * 4) % 4).reshape(5, 4) * 1.
x[x==3] = np.nan
y = x + (np.arange(5 * 4)/4).reshape(5,4)*1.
line2d = plt.plot(x.flatten(),y.flatten()) [0]
print x,y
> [[  0.,   1.,   2.,  nan],
> [  0.,   1.,   2.,  nan],
> [  0.,   1.,   2.,  nan],
> [  0.,   1.,   2.,  nan],
> [  0.,   1.,   2.,  nan]]
> [[  0.   1.   2.  nan]
> [  1.   2.   3.  nan]
> [  2.   3.   4.  nan]
> [  3.   4.   5.  nan]
> [  4.   5.   6.  nan]]

次に、問題なく set_data() メソッドを使用できます。次に例を示します。

line2d.set_data(y.flatten()+3)
于 2013-03-05T12:29:02.093 に答える
0

クリアして再プロットする必要がなくなるので、使用set_dataしてみてください（これは非常に遅いことに気づきました）。使用例:

fig = figure()
ax = fig.add_subplot(111)
p = ax.plot(x, y)

# you do something to your data and want to replot

p.set_data(x, y)

コードを投稿していただければ、正確な解決策を見つけるお手伝いをいたします。

于 2013-03-05T11:28:52.037 に答える