Poisson过程模拟和检验.docVIP

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Poisson过程模拟和检验

Poisson过程的模拟和检验 实验目的:理解掌握Poisson过程的理论,了解随机过程的模拟实现技术,学习并掌握在实际中如何检验给定的随机过程是否为Poisson过程。 实验内容:利用C语言、MATLAB等工具,结合Poisson过程等相关结论,模拟Poisson过程(还可选:非齐次Poisson过程等);查找资料、学习关于Poisson过程假设检验的相关知识,检验上述模拟实现的到达过程是否满足Poisson过程的定义(编程或利用统计软件,如SPSS、SAS等作为辅助工具)。 作业要求:提交实验报告电子版,说明模拟实现的过程,检验原理、步骤等以及实现过程;提交程序源代码。 一、泊松过程的模拟 1.基本原理 根据服务系统接受服务顾客数服从泊松分布这一模型可知,{X(n),t QUOTE }是一个计数过程,{ QUOTE ,n QUOTE 是对应的时间间隔序列,若 QUOTE (n)(n=1,2,...)是独立同分布的均值为 QUOTE 的指数分布,则{X(n),t QUOTE }是具有参数为λ的泊松。 2.具休实现过程 思路:本实验从用MATLAB编程软件,从构造服从指数分布的时间间隔 QUOTE 入手,计算每个事件的发生时刻,最后得到X(t),也就模拟了泊松过程。 实现步骤如下: (1).由函数random(‘exponential’,lamda)构造服从指数分布的 QUOTE 序列。 (2).根据服务系统模型, QUOTE = QUOTE + QUOTE 。 (3).对任意t QUOTE ( QUOTE , QUOTE ),X(t)=n,由此得到泊松过程的模拟。 3.过程模拟验证 (1)设定t=0时刻,计数为0,满足X(0)=0这一条件。 (2) QUOTE 是由random(‘exponential’,lamda)生成, QUOTE 间相互独立。 (3)由实验结果图可以很清楚地看出,在充分小的时间间隔内,最多有一个事情发生,而不可能有两个或两个以上事件同时发生,同时可以看出X(t)是一个平稳增量过程,结合条件(2)可知,X(t)是独立平稳增量过程。 图1:模拟泊松过程图 由此可知,根据服务系统模型,由具有指数分布的时间间隔序列模拟泊松过程可行。 二、泊松过程的检验 1.检验方法 Kolmogorov-Smirnov检验(柯尔莫哥洛夫-斯摩洛夫),亦称拟合优度检验法,用来检用来检验模拟所得的数据的分布是不是符合一个理论的已知分布。 检验步骤及过程: (1)条件设定: H1:实验产生模拟泊松分布数据的总体分布服从泊松分布。 H0:实验产生模拟泊松分布数据的总体分布不服从泊松分布。 (2)检验准备: 对于H1,已经假定所产生模拟泊松过程数据服从泊松分布,而强度未知,利用函数poissfit(x,alpha)估算出模拟泊松过程的强度,再利用函数poisscdf(x,lamda)得到泊松分布的累积分布函数。 (3) Kolmogorov-Smirnov检验 直接调用Kolmogorov-Smirnov检验函数kstest(x,[x,p],alpha),其中,x为输入模拟泊松序列,P为累积分布函数,1- alpha为置信区间,当结果时,则输入数据位泊松分布,否则,不是泊松分布。 三、程序代码 clear lamda=2;Tmax=50; delta_t=0.1;%时间精度 i=1;a=random(exponential,lamda); T(1)=round(a*10)/10; w(1)=T(1);%初始化 %%%%%%%%%%%%%%泊松过程模拟%%%%%%%%%%%%%%% while(w(i)Tmax) T(i)=random(exponential,lamda);%构造服从指数分布的时间间隔序列Tn T(i)=round(T(i)*10)/10; w(i+1)=w(i)+T(i);%计算等待时间 i=i+1; end w=w; x=zeros(w(1)/delta_t,1); for k=1:size(w,1)-1 length=w(k+1)/delta_t-w(k)/delta_t; x=[x;ones(length,1)*k];%得到泊松分布X(t)序列 end %%%%%%%%%%%泊松过程检验%%%%%%%%%%%%%%%%% alpha=0.05; lamda1=poissfit(x,alpha);%用MLE算法计算出泊松分布的强度lamda,置信区间为1-lamda p=poisscdf(x,lamda1);%计算累计分布 [H,s]=kstest(x,[x,p],alpha)%利用Kolmogorov-Smirnov检验,置信区间为1-

文档评论(0)

jgx3536 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:6111134150000003

1亿VIP精品文档

相关文档