参考代码%生成测试数据a0=1;a1=2;a2=3;x1=rand(10,10);x2=rand(10,10);Y=a0+a1*x1+a2*x2;%加入随机噪声Y=Y+0.1*randn(size(y));%拟合函数形式f=@(k,x)k(1)+k(2)*x(:,1)+k(3)*x(:,2);x=[x1(:)x2(:)];k=lsqcurvefit(f,[1111],x,Y(:))n=size(x,1);plot(1:n,Y(:),'bo',1:n,f(k,x),'r:.')运行结果k=1.03091.97822.9800分别对应a0~a2，可见与原始系数比较吻合。
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024