Matlab形成方程

Question

我将此等式写入Matlab时遇到问题。不知怎的，只是不能正确。有人可以帮我吗？

Answer 1

你可以在积分之外取总和，并将丑陋的前因定义为一个新变量：

i=0;
sum=0;
while true
   i=i+1;
%{
   fac=(2^i*u/W-1)^(2/alpha)
   sum=sum+int(1/(1+fac)/int([1/(1+x^alpha/2), 1/fac,1.0/0])...
           /(lambda*A)^i*factorial(i)*e^(-lambda*A-K*rho*A)...
           *rho^K*K^K/gamma(K)*A^{K-1},0,1.0/0);
%}
end

然后检查sum>=1-eta。

请注意，为方便起见，我将总和移动了1，因为最好在matlab中对1中的所有内容进行索引。