[討論] matlab算期望值

作者: vacuityhu (真空管)   2015-07-01 19:46:24
我需要算一個很少見的分布的期望值
叫做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看來是行不通了
作者: celestialgod (天)   2015-07-05 11:08:00
http://goo.gl/1Jhuog看起來就是不能直接算出來...你自己動手做,或是找文獻吧
作者: systematic (卡卡獸)   2015-07-07 02:59:00
可以先用機率密度轉換,然後用蒙地卡羅算出期望值隨便找本統計計算裡面有教

Links booklink

Contact Us: admin [ a t ] ucptt.com