此graph图片可使用矢量图形重新创建为SVG文件。这具有很多好处；更多信息请参见共享资源:待清理媒体。如果存在此图片的SVG格式，请将其上传，然后将此模板替换为{{vector version available|新图片名称}}。 建议您将SVG命名为“HittingTimes1.svg”，这样在使用Vector version available（或Vva）模板时就不需要填写新图片名称参数。

(An SVG version should ideally not include the title text, which should be a separate caption as text instead to allow this to be useful in any language).

已授权您依据自由软件基金会发行的无固定段落及封面封底文字（Invariant Sections, Front-Cover Texts, and Back-Cover Texts）的GNU自由文件许可协议1.2版或任意后续版本的条款，复制、传播和／或修改本文件。该协议的副本请见“GNU Free Documentation License”。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	本文件采用知识共享署名-相同方式共享 3.0 未本地化版本许可协议授权。
	您可以自由地： 共享 – 复制、发行并传播本作品 修改 – 改编作品 惟须遵守下列条件： 署名 – 您必须对作品进行署名，提供授权条款的链接，并说明是否对原始内容进行了更改。您可以用任何合理的方式来署名，但不得以任何方式表明许可人认可您或您的使用。 相同方式共享 – 如果您再混合、转换或者基于本作品进行创作，您必须以与原先许可协议相同或相兼容的许可协议分发您贡献的作品。
	本许可协议标签作为GFDL许可协议更新的组成部分被添加至本文件。http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

R-source code:

T=7
N=100000
set.seed(303898)
a=0.9
cs=c("orange", "navy", "seagreen4", "grey31")

bm1=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau1=which(bm1>a)[1]
bm1[tau1:(N+1)]=a

bm2=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau2=which(bm2>a)[1]
bm2[tau2:(N+1)]=a

bm3=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau3=which(bm3>a)[1]
bm3[tau3:(N+1)]=a

png(filename = "HittingTimes1.png", width=1200, height=900, pointsize = 12)
par(bg="gray88")
t=seq(0,T,length=N+1)
plot(t,bm1,type="l",ylim=range(bm1,bm2,bm3,a*1.2),col=cs[1],lwd=2,ylab="",xlab="Zeit",main="Hitting time von drei Sampels der Brownschen Bewegung als Beispiel für eine Stoppzeit")
lines(t,bm2,col=cs[2],lwd=2)
lines(t,bm3,col=cs[3],lwd=2)
abline(h=a,col=cs[4],lty=1,lwd=3)

text(x=tau1*T/N,y=a*1.14, labels = substitute(tau[2]==t1,list(t1=tau1*T/N)),pos=4,col=cs[1],cex=1.5)
text(x=tau2*T/N,y=a*1.14, labels = substitute(tau[1]==t2,list(t2=tau2*T/N)),pos=4,col=cs[2],cex=1.5)
text(x=T,y=a*1.14, labels = substitute(tau[3]>T,list(T=T)),pos=2,col=cs[3],cex=1.5)
dev.off()

 

本PNG 位图使用R创作.