→ celestialgod: 看起來就是不能直接算出來... 07/05 11:12
→ celestialgod: 你自己動手做，或是找文獻吧 07/05 11:12
→ systematic: 可以先用機率密度轉換，然後用蒙地卡羅算出期望值 07/07 02:59
→ systematic: 隨便找本統計計算裡面有教 07/07 03:00
我需要算一個很少見的分布的期望值
叫做non-central chi distribution
而且我要算的積分範圍不是全部
是一個常數c到inf
這邊我採用常見的int指令去算
syms x,i,a;
int((x-a)^2*pdf,x,a*c,inf)
pdf就是我這分布的pdf,a是另一個變數
到這裡都還好解決,可是後面問題來了
因為這分布的pdf中夾了一個
symsum(f(x,i),i,0,inf)
f是一個很醜的函數
而且這個summation雖然收斂,但是沒有公式解
所以這邊matlab總是整理不出來sum的式子
導致int那邊也做不出多項式,只能回傳int式子給我
因為我後面還要對a作微分求微分等於零的a'
所以他int之後不給我多項式的話diff跟solve就會error
想請教板上的各大神們
是否有別的辦法可以算這期望值??
因為我這原本的思路int->diff->solve看來是行不通了
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.23.7
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1435751187.A.4DB.html