matlab - 计算二项分布 rv 函数的期望值

Question

我正在尝试计算X的函数的期望值，其中X是二项式分布的。所以我想以sum(Pr(X=k)*f(k),k=0,..,n). 现在我想看看值是否随着二叉树的增加而收敛，即如果n增加。然而，对于增加n， ( nover k) 项在Pr(X=k)some 中趋于无穷大k，而函数f(k)=0, 这给出了在 matlabInf*0中产生的乘法。NaN因此结果是NaN。

但是，通过以矩阵形式计算类似问题，我知道期望值应该收敛。

所以我的问题是，我该如何处理？

这是我到目前为止的代码（对于一些参数，u~d~0.5. K=s=100）

sum=0;
for k=0:N;
    tmp=exp(-r*T)*nchoosek(N,k)*q_u^k*(1-q_u)^(N-k)*max(s*u^k*d^(N-k)-K,0);
    sum=sum+tmp;
end