some_tau = 1 / 0.5 ** 2
vel = pm.Normal('vel', mu=0, tau=some_tau)
states = [pm.Normal('s0', mu=true_positions[0], tau=some_tau)]
for i in range(1, len(true_positions)):
    pass
Hidden Markov in PyMC3