第1个回答 2019-08-16
你有60组数据，就是有
x
，y，z，k
四个数据，每个数据长度都是60
首先设置一个匿名函数fun，这里以k=a1*x*y+a2*y*z+a3*z*x为例子
fun=@(a,in)
a(1).*in(:,1).*in(:,2)+a(2).*in(:,2).*in(:,3)+a(3).*in(:,3).*in(:,1);
用a(n)去代替需要拟合的第n个参数，用in(,:1)代替x，in(,:2)代替y，in(,:3)代替z
根据你自己的函数形式，写好这个函数，记住乘除和乘方要用，.*，.\，.^
加个点
然后
in=[x(:)
y(:)
z(:)];
options
=
statset('TolFun',1e-4);
a=nlinfit(in,k,fun,[1
1
....1],options);
options那句控制精度
nlinfit的第一个参数是函数输入数据，这里是
in，存着x
y
z的数据
第二个参数是函数输出，这里是
k
第三个参数是你要拟合的方程，这里填fun，是我们之前辛苦写好的匿名函数
第四个参数是拟合参数的初值[1
1
....1]这里要填20个数
第五个参数是拟合设置，填入上一句设置好的options
输出a就是拟合好之后的参数，是一个长度为20的向量
有时候拟合结果会跟初值有关系，如果你的函数形式很复杂，
初值设不好可能会拟合不好，只能自己慢慢调整了