大家请记得明天下午第七节 在系电考机率上机考 --以下附上考试有可能会用到的M-file-- 1. poissonarrivals.m 2. poissonprocess.m function s=poissonarrivals(lambda,T) function N=poissonprocess(lambda,t) n=ceil(1.1*lambda*T); s=poissonarrivals(lambda,max(t)); s=cumsum(exponentialrv(lambda,n)); N=count(s,t); while (s(length(s))<T), s_new=s(length(s))+... cumsum(exponentialrv(lambda,n)); s=[s;s_new]; end s=s(s<=T); 3. count.m 4. exponentialrv.m function n=count(x,y) function x=exponentialrv(lambda,m) [MX,MY]=ndgrid(x,y); x=-(1/lambda)*log(1-rand(m,1)); n=(sum((MX<=MY),1))'; 5. simswitch.m 6. brownian.m function M=simswitch(lambda,mu,t) function w=brownian(alpha,t) s=poissonarrivals(lambda,max(t)); t=t(:); y=s+exponentialrv(mu,length(s)); n=length(t); A=count(s,t); delta=t-[0;t(1:n-1)]; D=count(y,t); x=sqrt(alpha*delta).*gaussrv(0,1,n); M=A-D; w=cumsum(x); 7. gaussrv.m (brownian 会用到的M-file) function x=gaussrv(mu,sigma,m) x=mu +(sigma*randn(m,1)); --考试形式-- 会给定一些基本参数如 t、alpha等 然後要求你们画出相关的图 但因为是随机变数的关系，所以每个人的图都不一样 所以会以所打的程式码来当作评分重点 --重点提醒-- 1.brownian 2.poissonarrivals 输入： t：向量形式 输入： T：观察结束的时间 alpha：brownian 的 parameter lambda：arrival rate 3.poissonprocess 输入： t：向量形式 lambda：arrival rate 时间t 输入的形式为： t=0.01*(0:1000); 意思是从0到1000以每0.01的间格来取点 即t=0, 0.01, 0.02, 0.03, ......, 1000 也可写成 t=0:0.01:1000; -- 记得要画图形之前 要先输入变数值 然後给定函数为何 最後决定画什麽图 然後标上 x轴的参数，与 y轴的参数 这样大概就画得出图了 如果还不行，先检查看看自己程式码有没有打错 大概就这样了，学弟妹加油 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.119.232.211