Exemplo n.º 1
0
def run():
#--------------------------------------------------------------------------------------------------
    from animate import launch
    syst = Pendulum(1,9.81)
    animate=dict(
        repeat=False,
        interval=40,
        )
    launch(syst,animate,
        srate=25.,
        taild=1.,
        ini=syst.makestate(179.),
        )
Exemplo n.º 2
0
def run():
#--------------------------------------------------------------------------------------------------
    from animate import launch
    syst = DoublePendulum(1.,.5,4.,1.,9.81)
    animate=dict(
        repeat=False,
        interval=40,
        )
    launch(syst,animate,
        srate=25.,
        taild=1.,
        ini=syst.makestate(theta1=180.,theta2=1.),
        )
Exemplo n.º 3
0
def run():
#--------------------------------------------------------------------------------------------------
    from animate import launch
    syst = RandomWalk(
        stepg=bvn(vx=3.1,vy=2.8,cor=.5),
        N=400,
        )
    animate=dict(
        repeat=False,
        interval=50,
        )
    def axes(fig):
        fig.set_size_inches((8,12),forward=True)
        return fig.add_subplot(2,1,1),fig.add_subplot(2,1,2)
    launch(syst,animate,axes=axes)
Exemplo n.º 4
0
def run():
#--------------------------------------------------------------------------------------------------
    from animate import launch
    syst = Kmeans(
        data=mixture(
            N=1000,
            weights=array((.2,.4,.4)),
            comp=(
                bvn(10.,0.,5.,3.,.2),
                bvn(0.,0.,6.,7.,.5),
                bvn(0.,10.,3.,8.,.7),
                ),
            ),
        K=5,
        )
    animate=dict(
        repeat=False,
        interval=300,
        )
    launch(syst,animate)
Exemplo n.º 5
0
def run():
#--------------------------------------------------------------------------------------------------
    from animate import launch
    syst = GP(
        datagen(
            func=0,
            nsample=1000,
            ),
        regr='quadratic',
        corr='squared_exponential',
        theta0=1.,
        thetaL_r=.001,
        thetaU_r=100.,
        nugget=.01,
        random_start=3
        )
    launch(syst,animate={},
        logscale=False,
        nsample=100,
        clevel=.95,
        gpmean=False,
        )