[問題] 常態pdf積分已回收

看板MATLAB作者 (hitori)時間15年前 (2010/04/08 19:43), 編輯推噓0(004)
留言4則, 1人參與, 最新討論串1/3 (看更多)
不好意思 我有積分的問題想請教 首先I1的部分是對常態分配的PDF作積分 積分下界為負無窮 而積分上界為(a+sqrt(ro)*z)/sqrt(1-ro) 為一變數z的function 而val1為I1的k次方 接著I2的部分則是val1*normpdf(z,0,1) 變數仍為z 積分上下界分別正負無窮 照理說這東西積出來的值應為0.93656 但答案卻不對 能否請板友看看問題出在哪裡呢?? 感激~ k =3; ro = 0.25; a = 2; syms z; I1 = int(normpdf(z,0,1),z,-1e+10,(a+sqrt(ro)*z)/sqrt(1-ro)); val1 = I1^k; I2 = int(val1*normpdf(z,0,1),z,-1e+10,1e+10); vpa(I2) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.116.96.50

04/08 20:38, , 1F
既然是常態分布 上限>0的話 從0開始積就好了 積完+0.5
04/08 20:38, 1F

04/08 20:39, , 2F
上限小於0的話 從0積到|上限| 然後1-基完的值
04/08 20:39, 2F

04/08 20:44, , 3F
另外 erf是個好指令
04/08 20:44, 3F

04/08 20:47, , 4F
修正 上限小於0的話 從0積到|上限| 然後0.5-積完的值
04/08 20:47, 4F
文章代碼(AID): #1BlS5NkN (MATLAB)
文章代碼(AID): #1BlS5NkN (MATLAB)