饿死抬头
function F = myfun(x,xdata)
Area = 481;
xlsFile = 'test.xls';
R = xlsread(xlsFile,1,'B2:B88');
for jj = 1:length(xdata);
A = 1;
for kk = jj:-1:1
qn(1,A) = (R(A,1)/x)*exp(-xdata(kk,1)/x)*Area*10^4/3600;
A = A+1;
end
end
F = sum(qn,2); %模拟值
end
以上是小弟的方程式

补充 xdata就是时间 [1:1:10］之类的就是说 第一个ydata对应xdata=1所算出来的F 第二个ydata对应xdata=2的F 然后xdata=1、2...所要算的F的方程式有一个循环的关系刚刚弄出来了 原来输出的F可以是一个矩阵