MATLAB最小二乘法拟合预测人口

已知进九年温州市人口，预测2024年温州市人口总量。

我们建立logistic人口模型，通过非线性最小二乘拟合，得到温州市人口的拟合曲线，进而求出下一年预测人口数量。

我们可以通过lsqcurvefit函数进行求解，由于此函数的工作原理是基于迭代法求最优解，因此对初始解十分敏感，为了求得较好的初始解，我们通过for循环不断改变初始参数值并记录残差平方和，最后寻找最小值的索引来求得较好的初始参数

（增长率函数中的r参数必须大于0，Xm必须大于往年的人口总量，并且根据常识一定不能超过现有人口的十倍，在此基础上拟合求参数）

求初始参数的程序如下：

clc,clear,close all
p = [911.7,917.5,921.5,925,930,957.29,964.5,967.9,976.1];
t = [2015:2023];
tt = [0:8];
fun =  @(x,tt)x(1)./(1+(x(1)/911.7-1).*exp(-x(2).*tt));
for i=1:200
    for j =1:300
options = optimoptions('lsqcurvefit','Algorithm','levenberg-marquardt');
[xishu,rnorm,r,exitflag] = lsqcurvefit(fun,[900+j,i],tt,p,[900,0],[9000,inf],[],[],[],[],[],options);
pr = fun(xishu,9);
rr(i,j)=rnorm;

end
end
minvalue = min(rr(:));
[row,col] = find(rr==mi