大批量数据处理

%程序2

%2001—2005年城镇、乡村的男、女人数程序代码

aa=[127627 128453 129227 129988 130756]*10000;

%计算2005年城镇男，城镇女，乡村男，乡村女数目

czn5=(b5(1)+b5(3))*aa(5)/sum(b5);%czn5表示2005年城镇男数目czm5=(b5(2)+b5(4))*aa(5)/sum(b5);%czm5表示2005年城镇女数目xcn5=(b5(5))*aa(5)/sum(b5);%xcn5表示2005年乡村男数目

xcm5=(b5(6))*aa(5)/sum(b5);%xcm5表示2005年乡村女的数目

f5=[czn5 xcn5;czm5 xcm5]

czn4=(b4(1)+b4(3))*aa(4)/sum(b4);

czm4=(b4(2)+b4(4))*aa(4)/sum(b4);

xcn4=(b4(5))*aa(4)/sum(b4);

xcm4=(b4(6))*aa(4)/sum(b4);

f4=[czn4 xcn4;czm4 xcm4]

czn3=(b3(1)+b3(3))*aa(3)/sum(b3);

czm3=(b3(2)+b3(4))*aa(3)/sum(b3);

xcn3=(b3(5))*aa(3)/sum(b3);

xcm3=(b3(6))*aa(3)/sum(b3);

f3=[czn3 xcn3;czm3 xcm3]

czn2=(b2(1)+b2(2))*aa(2)/sum(b2);

czm2=(b2(3)+b2(4))*aa(2)/sum(b2);

xcn2=(b2(5))*aa(2)/sum(b2);

xcm2=(b2(6))*aa(2)/sum(b2);

f2=[czn2 xcn2;czm2 xcm2]

czn1=(b1(1)+b1(2))*aa(1)/sum(b1);

czm1=(b1(3)+b1(4))*aa(1)/sum(b1);

xcn1=(b1(5))*aa(1)/sum(b1);

xcm1=(b1(6))*aa(1)/sum(b1);

f1=[czn1 xcn1;czm1 xcm1]

%程序3（计算2001——2005年总的平均死亡比率）

%死亡

clc

%计算2005年总的平均死亡率

nx=a5(:,2).*a5(:,3);%男性比例*男性死亡率=男性死亡比例

mx=a5(:,4).*a5(:,5);

sw5=sum(nx)+sum(mx);

nx=a5(:,6).*a5(:,7);

mx=a5(:,8).*a5(:,9);

sw52=sum(nx)+sum(mx);

nx=a5(:,10).*a5(:,11);

mx=a5(:,12).*a5(:,13);

sw53=sum(nx)+sum(mx);

swl5=sum(sw5+sw52+sw53)/300; nx=a4(:,2).*a4(:,3);

mx=a4(:,4).*a4(:,5);

sw4=sum(nx)+sum(mx);

nx=a4(:,6).*a4(:,7);

mx=a4(:,8).*a4(:,9);

sw42=sum(nx)+sum(mx);

nx=a4(:,10).*a4(:,11);

mx=a4(:,12).*a4(:,13);

sw43=sum(nx)+sum(mx);

swl4=sum(sw4+sw42+sw43)/300; nx=a1(:,2).*a1(:,3);

mx=a1(:,4).*a1(:,5);

sw1=sum(nx)+sum(mx);

nx=a1(:,6).*a1(:,7);

mx=a1(:,8).*a1(:,9);

sw12=sum(nx)+sum(mx);

nx=a1(:,10).*a1(:,11);

mx=a1(:,12).*a1(:,13);

sw13=sum(nx)+sum(mx);

swl1=sum(sw1+sw12+sw13)/300; nx=a3(:,2).*a3(:,3);

mx=a3(:,4).*a3(:,5);

sw3=sum(nx)+sum(mx);

nx=a3(:,6).*a3(:,7);

mx=a3(:,8).*a3(:,9);

sw32=sum(nx)+sum(mx);

nx=a3(:,10).*a3(:,11);

mx=a3(:,12).*a3(:,13);

sw33=sum(nx)+sum(mx);

swl3=sum(sw3+sw32+sw33)/300; nx=a2(:,2).*a2(:,3);

mx=a2(:,4).*a2(:,5);

sw2=sum(nx)+sum(mx);

nx=a2(:,6).*a2(:,7);

mx=a2(:,8).*a2(:,9);

sw22=sum(nx)+sum(mx);

nx=a2(:,10).*a2(:,11);

mx=a2(:,12).*a2(:,12);

sw23=sum(nx)+sum(mx);

swl2=sum(sw2+sw22+sw23)/300;

swl=[swl1,swl2,swl3,swl4,swl5]

%回归分析

clc

y1=(b1(1)+b1(2)).*(a1(:,2)+a1(:,4))+(b1(3)+b1(4)).*(a1(:,6) +a1(:,8))+(b1(5)+b1(6)).*(a1(:,10)+a1(:,12));

y2=(b2(1)+b2(2)).*(a2(:,2)+a2(:,4))+(b2(3)+b2(4)).*(a2(:,6) +a2(:,8))+(b2(5)+b2(6)).*(a2(:,10)+a2(:,12));

y3=(b3(1)+b3(2)).*(a3(:,2)+a3(:,4))+(b3(3)+b3(4)).*(a3(:,6) +a3(:,8))+(b3(5)+b3(6)).*(a3(:,10)+a3(:,12));

y4=(b4(1)+b4(2)).*(a4(:,2)+a4(:,4))+(b4(3)+b4(4)).*(a4(:,6) +a4(:,8))+(b4(5)+b4(6)).*(a4(:,10)+a4(:,12))+4.3956e+005; y5=(b5(1)+b5(2)).*(a5(:,2)+a5(:,4))+(b5(3)+b5(4)).*(a5(:,6) +a5(:,8))+(b5(5)+b5(6)).*(a5(:,10)+a5(:,12));

y=[y1;y2;y3;y4;y5];

t=0:90;t=t';

tt=ones(91,1);

x=[tt,t;2.*tt,t;3.*tt,t;4.*tt,t;5.*tt,t];

rstool(x,y,'purequadratic')

%程序4

%出生

clc

%新生儿

cs5=sum((b5(1)+b5(2)).*a5(16:50,4).*c5(:,2))/1000;%城市出生的新生儿数

sz5=sum((b5(3)+b5(4)).*a5(16:50,4).*c5(:,3))/1000;%乡镇出生的新生儿数

nc5=sum((b5(5)+b5(6)).*a5(16:50,4).*c5(:,4))/1000;%乡村出生的新生儿数

xs5=sum(cs5+sz5+nc5)%总新生儿数

nh5=cs5*e(8,2)/(e(8,2)+100)+sz5*e(8,3)/(e(8,3)+100)+nc5*e(8 ,4)/(e(8,4)+100)%2001年新生儿中（城市，城镇，农村）男孩的数

mh5=xs5-nh5%2001年新生儿中（城市，城镇，农村）女孩的数

nyc5=nh5/sum(b5)*10;%将千分数转化为百分数

myc5=mh5/sum(b5)*10;

cs4=sum((b4(1)+b4(2)).*a4(16:50,4).*c4(:,2))/1000;

sz4=sum((b4(3)+b4(4)).*a4(16:50,4).*c4(:,3))/1000;

nc4=sum((b4(5)+b4(6)).*a4(16:50,4).*c4(:,4))/1000;