clear;
clc;
warning?off;
close?all;
%%?從Excel表格中讀取數據
female?=?xlsread（‘A0301a.xls‘1‘D7:D126‘）;?%2010年各年齡段女性人口數量
rsex?=?xlsread（‘A0301a.xls‘1‘H7:H126‘）;?%2010年各年齡段性別比，女=100的男性數量
death?=?xlsread（‘A0604a.xls‘1‘J8:J127‘）;?%2010年各年齡段死亡率，千分之單位
%%?從數據中讀取各年齡女性人口數、女性人口比例和存活率
j?=?1;
for?i=1:1:length（female）
????if?mod（i-16）?~=?0
????????x（j1）?=?female（i）/1e8;
????????w（j1）?=?100/（rsex（i）+100）;
????????s（j1）?=?1-death（i）/1000;
????????j?=?j+1;
????end
end
%%?總和生育率
b1（1:151）?=?0;?%0-14歲
b1（16:501）?=?1.4/35;?%15-49歲，計劃生育政策
b1（51:1001）?=?0;?%50-100歲

b2（1:151）?=?0;?%0-14歲
b2（16:501）?=?1.8/35;?%15-49歲二胎政策
b2（51:1001）?=?0;?%50-100歲
%%?計算Leslie人口預測模型中Leslie矩陣
L1?=?zeros（100）;?%計劃生育政策
for?i=1:1:length（b1）
????L1（1i）?=?w（i）*b1（i）;
????if?i?~=?length（b1）
????????L1（i+1i）?=?s（i）;
????end
end

L2?=?zeros（100）;?%二胎政策
for?i=1:1:length（b2）
????L2（1i）?=?w（i）*b2（i）;
????if?i?~=?length（b2）
????????L2（i+1i）?=?s（i）;
????end
end
%%?預測每年的女性人口數量
t（11）?=?