caojihai 发表于 2009-9-9 09:07

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

结果与预期相差很大,尤其是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 发表于 2009-9-9 10:47

回复 楼主 caojihai 的帖子

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

caojihai 发表于 2009-9-10 23:22

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

ChaChing 发表于 2009-9-11 11:58

回复 板凳 caojihai 的帖子

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

VibrationMaster 发表于 2009-9-11 19:23

把两个求和的公式思路写出来
页: [1]
查看完整版本: 计算结果精度与预期相差很大