用于参数拟合、灵敏度分析（Matlab代码实现）

👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

所谓参数拟合，就是已知试验或者真实数据，然后寻找一个模型对其规律进行模拟，求取模型中未知参数的一个过程。

📚2 运行结果

 

  

 

 部分代码：

%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

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]胡宗庆.基于布谷鸟搜索的地位指数曲线模型参数拟合[J].西南林业大学学报,2013,33(01):66-70.

🌈4 Matlab代码实现