# A simulation of 1D Ito process with time step dt = .0001# mu is zero. sigma is a Ricker wavelet.importmatplotlib.pyplotasplimportnumpyasnpt0=0.0dt=0.0001t_final=10.0T=np.arange(t0,t_final,dt)ax=pl.figure().add_subplot(111)ax.set_xlabel('t')ax.set_ylabel('X')x=0.0mu=0.0# Ito Integralnp.random.seed(1)defRicker(t):return2.0/(np.sqrt(3)*(np.pi**.25))*(1-t*t)*np.exp(-0.5*t*t)fortinT:sigma=Ricker(t-t_final/2.0)new_x=x+mu*dt+sigma*np.random.normal(0,dt)ax.plot([t,t+dt],[x,new_x],'b-',linewidth=0.5)x=new_xpl.show()
Лицензирование
Я, владелец авторских прав на это произведение, добровольно публикую его на условиях следующей лицензии:
делиться произведением – копировать, распространять и передавать данное произведение
создавать производные – переделывать данное произведение
При соблюдении следующих условий:
атрибуция – Вы должны указать авторство, предоставить ссылку на лицензию и указать, внёс ли автор какие-либо изменения. Это можно сделать любым разумным способом, но не создавая впечатление, что лицензиат поддерживает вас или использование вами данного произведения.
распространение на тех же условиях – Если вы изменяете, преобразуете или создаёте иное произведение на основе данного, то обязаны использовать лицензию исходного произведения или лицензию, совместимую с исходной.
Файл содержит дополнительные данные, обычно добавляемые цифровыми камерами или сканерами. Если файл после создания редактировался, то некоторые параметры могут не соответствовать текущему изображению.