[交流] 【求助】使用MATLAB拟合出按照公式的系数

function plottu
clc;
B=xlsread('E:\CT.xls','sheet1')
fun=inline('a(1)*exp(-b(1)*B(i,1))+a(2)*exp(-b(2)*B(i,1))+a(3)*exp(-b(3)*B(i,1))+y(0)','a','b','y(0)','x');
a(0)=[10,10,10];
b(0)=[30,30,30];
y(0)=36;
[x,resnorm]=lsqcurvefit(fun,a(1),b(1),y(1),B(i,1),B(i,2));
end

Originally posted by 信彼南山 at 2011-01-11 22:28:03:
你的y(0)为什么是36呢？

Originally posted by 信彼南山 at 2011-01-11 22:35:15:
你那个B(i,1)是啥东西？是不是想说B(:,1)啊？

你的函数是不是应该这样啊？
y(i)=a1*exp(-b1*x(i))+a2*exp(-b2*x(i))+a3*exp(-b3*x(i))+y(0)

Originally posted by 信彼南山 at 2011-01-11 22:35:15:
你那个B(i,1)是啥东西？是不是想说B(:,1)啊？

你的函数是不是应该这样啊？
y(i)=a1*exp(-b1*x(i))+a2*exp(-b2*x(i))+a3*exp(-b3*x(i))+y(0)

Originally posted by zyj8119 at 2011-01-11 22:46:51:

文件的第一列数据。

Originally posted by 信彼南山 at 2011-01-11 23:01:29:


大哥，你的第一列是编号，不是数据

Originally posted by 信彼南山 at 2011-01-11 23:00:46:
% zyj_main.m
clear;clc;
B=xlsread('CT.xls','sheet1');
xdata=B(:,2);
ydata=B(:,3);
a=[10,10,10];
b=[30,30,30];
y0=36;
x=lsqcurvefit(@zyj_fun1,[a,b,y0],xdata,ydata);


-------------------- ...

Originally posted by 信彼南山 at 2011-01-13 02:45:53:
我还真是这么写的程序

不过，按照你的这些数据，我觉得根本拟合不出来你那些参数，或者说可以有无穷多组结果能满足条件，建议换一个算法试试，或者在公示中多加一些约束条件来判别

Originally posted by zyj8119 at 2011-01-13 17:08:13:

但是别人拟合出来啊，是以前人。