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()
Лиценцирање
Ја, носилац ауторског права над овим делом, објављујем исто под следећом лиценцом:
да делите – да умножавате, расподељујете и преносите дело
да прерађујете – да прерадите дело
Под следећим условима:
ауторство – Морате да дате одговарајуће заслуге, обезбедите везу ка лиценци и назначите да ли су измене направљене. Можете то урадити на било који разуман манир, али не на начин који предлаже да лиценцатор одобрава вас или ваше коришћење.
делити под истим условима – Ако измените, преобразите или доградите овај материјал, морате поделити своје доприносе под истом или компатибилном лиценцом као оригинал.