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