LeslieModel.zip

莱斯利矩阵是科学家LesliePH.于1945年引进一种数学方法，利用某一初始时刻种群的年龄结构现状，动态地预测种群年龄结构及数量随时间的演变过程。 LESLIE模型的matlab代码，可用于数学建模中的人口预测。

function [Y T Z] = LeslieModel（data b year number ratio migration）
%b是15-49岁妇女生育率
%data表示输入的4列分别为男性比例，男性死亡率，女性比例，女性死亡率
%number表示人口数，ratio表示男性的人口数/女性的人口数year表示要预测的年份
%输出T表示青壮年人口，Y表示未来年份所有的女性人口总数Z表示老龄人
L = Leslie（data b / （1 + ratio） / 1000）;
r = data（: 3） ./sum（data（: 3））; %r表示各个年龄段女性的比例
X0 = number * 1/（1 + ratio） .* r; %X0是初始值
X = X0;
Y = ones（1 year）;
for i = 1 : year
    X = L * X + migration * data（: 3） / 100;%如果不考虑迁移的情况，那么migration，那么就取migration为0向量
    Y（i） = （ratio + 1）* sum（X）;
    Z（i） = （ratio + 1） * sum（X（66:end））;
    T（i） = （ratio + 1） * sum（X（16:65））;    
end

%对城镇乡的leslie矩阵的稳定性分析程序
load bear_rate
load city
load town
load country
LC = Leslie（city_05 1/（1+1.14）*bear_rate（: 1））;
LT = Leslie（town_05 1/（1+1.17）*bear_rate（: 2））;
LR = Leslie（country_05 1/（1+1.21）*bear_rate（: 3））;
[V1 L1] = eigs（LC）;
[V2 L2] = eigs（LT）;
[V3 L3] = eigs（LR）;
L1L2L3

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1067  2015-05-15 09:49  LeslieModel.m