大家請記得明天下午第七節 在系電考機率上機考 --以下附上考試有可能會用到的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