English: 2D Ornstein-Uhlenbeck process with time step of .0001, while theta = 1.0 and sigma = 300. The initial position is (10, 10). Finally the point wanders around the central point (0, 0).
# A simulation of 2D Ornstein-Uhlenbeck process with time step dt = .0001importmatplotlib.pyplotasplimportnumpyasnpt0=0.0dt=0.0001t_final=2T=np.arange(t0,t_final,dt)ax=pl.figure().add_subplot(111)ax.set_xlabel('X')ax.set_ylabel('Y')x,y=10.0,10.0ux,uy=0.0,0.0theta=1.0sigma=300.0np.random.seed(1)fortinT:new_x=x+theta*(ux-x)*dt+sigma*np.random.normal(0,dt)new_y=y+theta*(uy-y)*dt+sigma*np.random.normal(0,dt)ax.plot([x,new_x],[y,new_y],'b-',linewidth=0.5)x,y=new_x,new_ypl.show()
Licenciranje
Ja, nosilac autorskog prava nad ovim delom, objavljujem isto pod sledećom licencom:
da delite – da umnožavate, raspodeljujete i prenosite delo
da prerađujete – da preradite delo
Pod sledećim uslovima:
autorstvo – Morate da date odgovarajuće zasluge, obezbedite vezu ka licenci i naznačite da li su izmene napravljene. Možete to uraditi na bilo koji razuman manir, ali ne na način koji predlaže da licencator odobrava vas ili vaše korišćenje.
deliti pod istim uslovima – Ako izmenite, preobrazite ili dogradite ovaj materijal, morate podeliti svoje doprinose pod istom ili kompatibilnom licencom kao original.