您所在位置：主题：请教！挑战matlab!!
dongsheng398
[专家分：0] 发布于 2008-07-22 16:17:00
function F(to)=Y1+Y2
t=[3514 3852 4004 4334 4515 4727 4805 5000 5206 5411 5520 5722 6027 6109 6327 6523 6717 7326 7743 8489];
Y1=1;m1=1;m2=1;Y2=1
for i=1:20
Y1=Y1*(log(t(i)-t0)/(t(i)-t0));
m1=m1*log(t(i)-t0);
m2=m2*(1/(t(i)-t0));
end;
Y2=Y2*m1*m2/20;
F(t0)=Y1+Y2;
t=[3514 3852 4004 4334 4515 4727 4805 5000 5206 5411 5520 5722 6027 6109 6327 6523 6717 7326 7743 8489];
Y1=1;m1=1;m2=1;Y2=1
for i=1:20
Y1=Y1*(log(t(i)-t0)/(t(i)-t0));
m1=m1*log(t(i)-t0);
m2=m2*(1/(t(i)-t0));
end;
Y2=Y2*m1*m2/20;
F(t0)=Y1+Y2;