[功課] 機率上機

看板NCCU09_Math作者chenghe (路人)時間13年前 (2011/05/04 18:47)推噓8(8推 0噓 0→)

留言8則, 7人參與討論串1/1

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

推

ahnjcfhjnkma

05/04 22:45, , 1^F

05/04 22:45, 1^F

推

yenli0501

05/04 23:13, , 2^F

05/04 23:13, 2^F

推

davidhuang28

05/04 23:29, , 3^F

05/04 23:29, 3^F