作者vacuityhu (真空管)
看板MATLAB
标题[讨论] matlab算期望值
时间Wed Jul 1 19:46:24 2015
我需要算一个很少见的分布的期望值
叫做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://webptt.com/cn.aspx?n=bbs/MATLAB/M.1435751187.A.4DB.html
2F:→ celestialgod: 看起来就是不能直接算出来... 07/05 11:12
3F:→ celestialgod: 你自己动手做,或是找文献吧 07/05 11:12
4F:→ systematic: 可以先用机率密度转换,然後用蒙地卡罗算出期望值 07/07 02:59
5F:→ systematic: 随便找本统计计算里面有教 07/07 03:00