用Pyx繪製花括號

Question

如何用Pyx在兩個任意點之間繪製「支撐」線？

這將是這個樣子：

Brace example http://tof.canardpc.com/view/d16770a8-0fc6-4e9d-b43c-a11eaa09304d

Answer 1

您可以繪製使用sigmoidals相當大括號。我沒有安裝Pyx，所以我只是使用matplotlib（這裏是pylab）繪製這些圖。這裏beta控制大括號中曲線的清晰度。

import numpy as nx 
import pylab as px 


def half_brace(x, beta): 
    x0, x1 = x[0], x[-1] 
    y = 1/(1.+nx.exp(-1*beta*(x-x0))) + 1/(1.+nx.exp(-1*beta*(x-x1))) 
    return y 

xmax, xstep = 20, .01 
xaxis = nx.arange(0, xmax/2, xstep) 
y0 = half_brace(xaxis, 10.) 
y = nx.concatenate((y0, y0[::-1])) 

px.plot(nx.arange(0, xmax, xstep), y) 
px.show() 

alt text http://i26.tinypic.com/23iyp76.png

我策劃了這沿x軸，以節省屏幕空間，但相處y軸支架只是交換x和y。最後，Pyx還有很多內置的路徑繪製功能，它們也可以滿足您的需求。

Answer 2

tom10提供了一個很好的解決方案，但可以使用一些改進。
關鍵是在範圍[0,1]，[0,1]上創建一個大括號，然後對其進行縮放。
這個版本還可以讓你調整一下形狀。對於獎勵積分，它使用二階導數來計算多少密集間隔點。

mid設置下部和上部之間的平衡。
beta1和beta2控制曲線（下和上）有多銳。
您可以更改height（或者僅將y乘以標量）。
使它垂直而不是水平只需要交換x和y。
initial_divisions和resolution_factor指導如何選擇x值，但通常應該是可以忽略的。

import numpy as NP 

def range_brace(x_min, x_max, mid=0.75, 
       beta1=50.0, beta2=100.0, height=1, 
       initial_divisions=11, resolution_factor=1.5): 
    # determine x0 adaptively values using second derivitive 
    # could be replaced with less snazzy: 
    # x0 = NP.arange(0, 0.5, .001) 
    x0 = NP.array(()) 
    tmpx = NP.linspace(0, 0.5, initial_divisions) 
    tmp = beta1**2 * (NP.exp(beta1*tmpx)) * (1-NP.exp(beta1*tmpx))/NP.power((1+NP.exp(beta1*tmpx)),3) 
    tmp += beta2**2 * (NP.exp(beta2*(tmpx-0.5))) * (1-NP.exp(beta2*(tmpx-0.5)))/NP.power((1+NP.exp(beta2*(tmpx-0.5))),3) 
    for i in range(0, len(tmpx)-1): 
     t = int(NP.ceil(resolution_factor*max(NP.abs(tmp[i:i+2]))/float(initial_divisions))) 
     x0 = NP.append(x0, NP.linspace(tmpx[i],tmpx[i+1],t)) 
    x0 = NP.sort(NP.unique(x0)) # sort and remove dups 
    # half brace using sum of two logistic functions 
    y0 = mid*2*((1/(1.+NP.exp(-1*beta1*x0)))-0.5) 
    y0 += (1-mid)*2*(1/(1.+NP.exp(-1*beta2*(x0-0.5)))) 
    # concat and scale x 
    x = NP.concatenate((x0, 1-x0[::-1])) * float((x_max-x_min)) + x_min 
    y = NP.concatenate((y0, y0[::-1])) * float(height) 
    return (x,y) 

用法很簡單：

import pylab as plt 

fig = plt.figure() 
ax = fig.add_subplot(111) 

x,y = range_brace(0, 100) 
ax.plot(x, y,'-') 

plt.show() 

PS：不要忘了，你可以通過clip_on=False到plot，並把它的軸線之外。