lsqcurvefit最小二乘拟合

%拟合公式$y=a+b\ast x+c\ast x^2+d\ast x^3$
x2=x.*x;
x3=x.*x.*x;
%X为参数矩阵
f=@(X)(X(1)+X(2)*x+X(3)*x2+X(4)*x3-y);
x0=[1,1,1,1];%赋予参数初值
[X,resnorm]=lsqcurvefit(f,x0,x,y);