de: Pfade von zwei kompensierten zusammengesetzten Poisson-Prozessen. Wie in de:Poissonprozess.png ist die Intensität (Sprunghäufigkeit) des blauen Prozesses mit 2,4 genau viermal so hoch wie die des roten Prozesses. Im gezeichneten Intervall [0,35] springt der blaue Prozess 66-mal, der rote 16-mal, also etwa viermal so oft. Bei beiden Prozessen sind die Sprünge normalverteilt mit Mittel 0,25. Diese Sprünge nach oben werden durch den negativen Drift genau so ausgeglichen (kompensiert), dass beide Prozesse Martingale sind. Da der blaue Prozess öfter nach oben springt, ist sein negativer "Drift" stärker.

已授權您依據自由軟體基金會發行的無固定段落、封面文字和封底文字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 source code:

set.seed(29011980)

##plots a compound (normal), compensated poisson process.
poisson.process<-function(lambda, T, mu, sigma, drift=-lambda*mu) {
  N=rpois(n=1, lambda=lambda*T)                                   #number of jumps
  t=sort(array(data=runif(n=N, min=0,max=T),dim=N))               #jump times
  y=array(data=rnorm(n=N,mean=mu,sd=sigma), dim=N)                #jump sizes
  t=t-t[1]
  s=array(NA, dim=2*N)
  s[2*(1:N)-1]=diff(c(t,T))*drift
  s[2*(1:N)]  =y
  s=c(0,cumsum(s))
  ti=seq(from=1, to=N)
  return(list(s=s,t=t,ti=ti))
}

pp1=poisson.process(lambda=2.4,T=35,mu=0.25,sigma=1)
pp2=poisson.process(lambda=0.6,T=35,mu=0.25,sigma=1)

png(filename="Comp-Comp-Poisson-Pro.png", width=1600, height=1200, pointsize=12)
par(bg="whitesmoke")
plot(c(0,pp1$t),c(0,pp1$t),type="n",xlab="t",ylab="",ylim=range(pp1$s,pp2$s))
segments(pp1$t[pp1$ti], pp1$s[2*pp1$ti-1], pp1$t[pp1$ti+1], pp1$s[2*pp1$ti], col="blue", lwd=2)
segments(pp2$t[pp2$ti], pp2$s[2*pp2$ti-1], pp2$t[pp2$ti+1], pp2$s[2*pp2$ti], col="red", lwd=2)

title(main="Pfade von kompensierten zusammengesetzten Poisson-Prozessen",cex.main=1.7)
legend(x="bottomleft", inset=0.01, cex=1.4, title="Intensität", legend=c(expression(lambda[1]==2.4,lambda[2]==0.6)), col=c("blue","red"), lwd=3) 
text(x=par("usr")[1]*0.9, y=par("usr")[4]*0.95, labels="Sprunghöhe ~ N(0.25,1)", adj=c(0,1))
dev.off()

 

本PNG 點陣圖使用R創作。