MATLAB中有一个多元非线性拟合的功能是nlinfit 基本语法是：

beta = nlinfit(X,Y,modelfun,beta0)

式子左边的beta可以是一个向量，向量的元素就是要回归的模型中的参数。 式子右边，modelfun是要回归的函数形式。X是函数的自变量数据；Y是函数的因变量数据；beta0是待回归参数的初始值。

modelfun如果是MATLAB内置的函数形式，则直接使用，例如，内置的Hougen-Watson model，形式如下：

使用时，这么用

beta = nlinfit(X,Y,@hougen,beta0)

这里@是函数句柄的操作符。 因为Hougen-Watson model中有5个待回归的参数（ β 1 \beta_{1} β1​， β 2 \beta_{2} β2​， β 3 \beta_{3} β3​， β 4 \beta_{4} β4​， β 5 \beta_{5} β5​），因此beta0是一个1×5的向量。模型中有3个未知数，因此X是个n×3的矩阵，Y是n×1的矩阵。n是数据的条数。

如果modelfun不是MATLAB内置的函数，就需要自己建立一个函数的形式。例如，下面的语句：

modelfun = @(beta,x)((beta(1)*x1+(1-beta(1))*x2)*beta(2));