Matlab将数据拟合到函数，其中函数是求和和周期性语法

Question

我正在尝试使用lsqcurvefit拟合美国过去50年的温度数据，但是我担心我的拟合函数语法错误。

拟合函数本身具有以下形式：

∑(ai+bi*t+ci*t^2)*(cos(vt+p))

从i = 1到N的总和。我希望t对应于称为TIME的时间向量。 理想情况下，拟合应返回描述温度数据的a，b，c，v和p值

我正在尝试拟合称为TMAX的每日最高温度。

到目前为止，我有：

lsqcurvefit(fitFn(1,2,TIME),5,TIME,TMAX)

在脚本中将fitFn定义为

function f = fitFn(i,N,t)
f=0;
for i=1:N
    f =@(i,N,a,b,c,v,p,t) f + (a+b*t+c*t^2)*(cos(v*t+p));
end

但是，每当我运行lsqcurvefit时，我都会收到错误消息

Error using
fitFn>@(i,N,a,b,c,v,p,t)f+(a+b*t+c*t^2)*(cos(v*t+p)) (line
4)
Not enough input arguments.

Error in lsqcurvefit (line 199)
            initVals.F =
            feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});

Caused by:
    Failure in initial user-supplied objective function
    evaluation. LSQCURVEFIT cannot continue.

我确定我犯了一个简单的错误，但我有些茫然。 我是MatLab的新手，不太了解如何使用匿名函数。 任何帮助将不胜感激。

谢谢阅读

Answer 1

您的fitFn非常错误。 我认为这是您想要的：

f = @(x,t) (x(1)+x(2)*t+x(3)*t.^2)*cos(x(4)*t+x(5))
lsqcurvefit(f,zeros(1,5),TIME,TMAX)