计算结果精度与预期相差很大

caojihai

结果与预期相差很大，尤其是Epsilon1大的时候，比如100以上。高人指点迷津。回帖或者发信到 caojihai@hotmail.com  十分感谢。
clc; clear; %Epsilon1 = 0:10:50;
Epsilon1 = 100;
%sv= 0:0.4:3;
Ns=length(Epsilon1); Nu = 1; Lambda=0.2;
mv=0:1:153; Nm=length(mv); sum1=zeros(1,Nm); sum3=zeros(1,Nm);
y=zeros(1,Ns); y1=zeros(1,Ns); y2=zeros(1,Ns); y3=zeros(1,Ns);
%lower bound
for indexs=1:1:Ns, Epsilon=Epsilon1(indexs);
    for indexm=1:1:Nm
        m=mv(indexm); j=0:1:m;
        sum2=1./gamma(m-j+1)./gamma(j+1).*gamma(m+Nu-j).*Lambda.^j.*(1+Nu/Epsilon).^j;
        sum1(indexm)=exp(-Lambda)*Nu^Nu*Epsilon^m/(Epsilon+Nu)^(m+Nu)/gamma(Nu)*sum(sum2)*m;
    end
    y(indexs)=-sum(sum1);
end

[ 本帖最后由 ChaChing 于 2009-9-11 11:53 编辑 ]

ChaChing

个人水平专业有限, 建议楼主说清楚些!
结果与预期相差很大! 结果? 预期?
有图吗?

caojihai

正确的结果 y=Epsilon1＋Lambda； 但是算出来的y却与此相差很大，尤其是Epsilon1越大，差的越多。

ChaChing

个人水平专业有限, 还是不清楚楼主的问题!
待高人路过

VibrationMaster

把两个求和的公式思路写出来