💥1 概述
所谓参数拟合，就是已知试验或者真实数据，然后寻找一个模型对其规律进行模拟，求取模型中未知参数的一个过程。
📚2 运行结果
🎉3 参考文献
[1]胡宗庆.基于布谷鸟搜索的地位指数曲线模型参数拟合[J].西南林业大学学报,2013,33(01):66-70.
👨💻4 Matlab代码
主函数部分代码：
%Liberia %Pop = 4.294e6; clear;clc; g1 = 0.1; psi = 0.486; b1 = 0.718; b2 = 1.34; b3 = 1.05e-5; g2 = 0.0524; r1 = 0.77; R0 = (1/(g1+psi))*(b1 + (b3*psi/g2) + (b2*r1*g1)) R0pconst = b1 + (b2*r1*g1) - (g1*b3/g2); R0primepsi = - R0pconst/(g1 + psi)^2 R0at0 = (b1/g1) + (b2*r1) R0atinf = b3/g2 x = linspace(0,2,100); for i = 1:length(x) yL(i) = (1/(g1+x(i)))*(b1 + (b3*x(i)/g2) + (b2*r1*g1)); yLp(i) = - R0pconst/(g1 + x(i))^2; end figure; plot(x,yL, x, yLp,'--'); legend('R_0(\psi)', 'R_0\prime (\psi)'); psiL = psi; R0L = R0; R0pL = R0primepsi; %% %Guinea %Pop = 11.76e6; %clear;clc; g1 = 0.1; psi = 1.11; b1 = 0.316; b2 = 0.446; b3 = 0.0325; g2 = 0.0269; r1 = 0.649; R0 = (1/(g1+psi))*(b1 + (b3*psi/g2) + (b2*r1*g1)) R0pconst = b1 + (b2*r1*g1) - (g1*b3/g2); R0primepsi = - R0pconst/(g1 + psi)^2 R0at0 = (b1/g1) + (b2*r1) R0atinf = b3/g2 x = linspace(0,2,100); for i = 1:length(x) yG(i) = (1/(g1+x(i)))*(b1 + (b3*x(i)/g2) + (b2*r1*g1)); yGp(i) = - R0pconst/(g1 + x(i))^2; end
