为什么这个公式总是计算不对？

furoo

计算一个公式（见附件），m取值是2<=m<=100，步长是1.程序如下：
for m=2:100;
   for i=1:m, s1(i)=log(r(i))-log(r(m+1)); end;
   for i=1:m, s2(m)=sum(s1(i)); end;
   alfa(i)=((1/m)*s2(m))^(-1)
end
为什么计算结果总是不对呢？我初学matlab,还望多多指教，非常感谢

运行结果：m=100,然后alfa是1*100个值，正常应该是m是2到100，alfa是99个值啊

[ 本帖最后由 ChaChing 于 2009-3-11 14:08 编辑 ]

friendchj

try：

clc; clear; r=rand(1,201);
alfa=zeros(1,99);
for m=2:100;
   for i=1:m;
   s1(i)=log(r(i))-log(r(m+1));
   end;
   alfa(m-1)=(sum(s1)/m)^(-1);
end

[ 本帖最后由 ChaChing 于 2009-3-11 14:23 编辑 ]

ChaChing

个人习惯, 参考下
clc; clear; r=rand(1,201); alfa=zeros(1,99);
for m=2:100, ttt=log(r(m+1));
   for i=1:m,  s1(i)=log(r(i))-ttt; end
   alfa(m-1)=1/(sum(s1)/m);
end