迭代计算
xdata=1:10;
ydata=xdata+(1+(rand()-0.5)); %输入x,y值
c0=[1 1]; %设置初值;
for i=1:100
c = lsqcurvefit ('fun', c0, xdata, ydata) ;
c0 = c; %以计算出的c为初值进行迭代;
end
c
y1=polyval(c,xdata);
plot(xdata, ydata,'o',xdata,y1,'r');
legend ('实验数据 ','拟合曲线');
设置过定点的函数
function f = fun( c,x )
%在这里设置定点
f=c(1).*(x-1)+2;
end
这个函数过定点(1,2)