写在前面
做了一个作业,觉得挺有意思,有空就发上来了。
阻滞增长模型–Logistic模型的部分推导、思路,并给出参考代码。
参考书籍:《数学模型》第四版,姜启源 谢金星 叶俊 编
MATLAB版本2021b
点我下载源代码及原始数据
提取码:p3uu
一、思路
阻滞作用体现在对人口增长率r的影响上,使得r随着人口数量x的增加而下降。于是有:
对r(x)做假设,设r(x)是x的线性函数,即,
这里r称固有增长率;引入人口容量xm,当x=xm时人口不在增长,于是有,
结合上述三式,有,
可解得,
式5即为Logistic阻滞增长模型。
取数据得第一个时间点得值为x0,即x0=3.9。
由式4可得,
等号左边项可由实际数据计算得出,左端项与右端x呈线性关系,可通过线性拟合得r,s,进而求得x_m。
至此,式5中固有人口增长率r、人口容量x_m、人口初始值x_0均已求得。
二、原始数据
点我下载mat格式原始数据
提取码:75d5
数据在参考书P164
三、源代码
// function_xt.m
function output = function_xt(xm, x0, r, t )
%UNTITLED2 此处提供此函数的摘要
% Logistic模型
% 输入:
% xm:人口容量
% r :<