0

我有以下无限系列需要评估。目前，根据输入参数（kappa、sigmav、thetav），这个函数需要很长时间来计算。有没有可能加快速度？耗时的表达式是 symsum(..., 0, Inf)

 j = sym('j', 'positive');
 c =  1/((1/2*kappav)*(sigmav^2)*(1-exp(-kappav*deltat)));
 lambda = 2*c*Vt*exp(-kappav*deltat);
 v = 4*thetav/sigmav^2;
 ez = 2*c*exp(Vtt);
 % 
 p = ((exp(-1/2*(ez+lambda))*ez^(1/2*v)) / ( 2^(0.5*v)))* ... 
   vpa( symsum(   ((1/4*ez*lambda)^j)/ (gamma(1/2*v+j)*factorial(j)) ,j, 0, Inf  ));

PS。这是非中心卡方分布的 pdf ，有一些小的修改，因此我不能使用ncx2pdf. 我已经尝试过https://stackoverflow.com/a/15966126/321749这只会导致无穷大。

4

0 回答 0