已授權您依據自由軟體基金會發行的無固定段落、封面文字和封底文字GNU自由文件授權條款1.2版或任意後續版本，對本檔進行複製、傳播和／或修改。該協議的副本列在GNU自由文件授權條款中。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	此檔案採用創用CC 姓名標示-相同方式分享 3.0 未在地化版本授權條款。
	您可以自由： 分享 – 複製、發佈和傳播本作品 重新修改 – 創作演繹作品 惟需遵照下列條件： 姓名標示 – 您必須指名出正確的製作者，和提供授權條款的連結，以及表示是否有對內容上做出變更。您可以用任何合理的方式來行動，但不得以任何方式表明授權條款是對您許可或是由您所使用。 相同方式分享 – 若要根據本素材進行再混合、轉換或創作，則必須以與原作相同或相容的授權來發布您的作品。
	已新增授權條款標題至此檔案，作為GFDL授權更新的一部份。http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

R-Quelltext:

#Orenstein-Uhlenbeck Prozess mit normalverteiltem Anfangswert
ou<-function(theta=1, mu=1.2, sigma=0.3, init=c(mu, sigma^2/(2*theta)), from=0, to=2, steps=500*(to-from)) {
  t<-seq(from,to,length=steps)
  oup<-rnorm(n=1,mean=init[1],sd=sqrt(init[2]))
  for (i in (1:(steps-1))) {
    oup[i+1]=oup[i] + theta*(mu-oup[i])*(to-from)/steps+sigma*rnorm(n=1,mean=0,sd=sqrt(sigma^2*(to-from)/steps))
  }
  return(oup)
}

cols=c("navy", "darkolivegreen4", "red")
t<-seq(0,2,length=1000)
set.seed(311051)

png(filename="OrnsteinUhlenbeck3.png", width=1200, height=900, pointsize=12)
par(bg="lightgrey")
matplot(t,ou(init=c(0,0)),type="l", xlab="", ylab="",col=cols[1], lwd=2, ylim=c(0,2))
lines(t,ou(init=c(2,0)),col=cols[2], lwd=2)
lines(t,ou(),col=cols[3], lwd=2)
title(main="Ornstein - Uhlenbeck", cex.main=2.5)
text(x=1.85,y=0,expression(list(theta==1, mu==1.2, sigma == 0.3)) )
text(x=1.55,y=1.75,expression(paste("d",X[t]==theta(mu-X[t]),"dt+",sigma,"d", W[t])),cex=2.5)
text(x=0,y=0.07, expression(X[0]==0),col=cols[1])
text(x=0,y=1.92, expression(X[0]==2),col=cols[2])
text(x=0.03,y=1.29, expression(paste(X[0], "~N(", mu,",", frac(sigma^2,2*theta), ")") ),col=cols[3])
abline(h=1.2, lty=3, col="grey25")
dev.off()

 

本PNG 點陣圖使用R創作。