【正文】 998199920002001200220032004200520062007200820092010表4:014歲人口數(shù)的殘差檢驗(yàn)?zāi)攴?14歲人口數(shù)實(shí)際值模型計(jì)算值殘差相對殘差200120022003200420052006200720082009表5：對1560歲的人口數(shù)的殘差計(jì)算年份1560歲的人口數(shù)的實(shí)際值模型計(jì)算值殘差相對殘差200120022003200420052006200720082009表6：對60歲以上的人口數(shù)的殘差計(jì)算年份60歲以上的人口數(shù)的實(shí)際值模型計(jì)算值殘差相對殘差200120022003200420052006200720082009表7：對男性人口比例的殘差計(jì)算年份男性人口比例的實(shí)際值模型計(jì)算值殘差相對殘差1980198119821983198419851986198719881989199019911992199319941995199619971998199920002001200220032004200520062007200820092010附錄1一、問題一程序：1. 層次分析法程序：clear。clc。a=[1 1/2 3/5 32 1 3/2 65/3 2/3 1 51/3 1/6 1/5 1]。[b,wi,lammax,CI,RI,CR]=AHP(a)%計(jì)算程序段輸出以下量% b:歸一化矩陣；% wi:權(quán)重向量；% lammax：最大特征值；% CI：隨機(jī)數(shù)一致性指標(biāo)；% RI：平均隨機(jī)數(shù)一致性指標(biāo)；% CR：隨機(jī)數(shù)一致性指標(biāo)比；function [b,wi,lammax,CI,RI,CR]=AHP(a)%和積法計(jì)算wi和規(guī)范化后的wi%輸出a的規(guī)范化矩陣b和wii,wi,lammax,CI；n=size(a,1)。b=ones(n,n)。for i=1:n。 for j=1:n b(i,j)=a(i,j)/sum(a(:,j))。 endendwii=sum(b,2)。%計(jì)算wiwi=wii/sum(wii)。%規(guī)范wiww=a*wi。for i=1:n ww(i)=ww(i)/wi(i)。endlammax=sum(ww)/n。 CI=(lammaxn)/(n1)。 RI=randRI(n)。CR=CI/RI。%%計(jì)算平均隨機(jī)數(shù)一致性指標(biāo)RI%CR（一致性指標(biāo)CI與同階平均隨機(jī)數(shù)指標(biāo)RI之比）=,則判斷矩陣具有滿意的一致性；%CI=（lamdamaxn）/(n1)。%注意：因T較大且程序多為for語句，使循環(huán)次數(shù)較多，程序運(yùn)行時間可能較長，如有需要可取較小的T值；%和積法計(jì)算RI值；%rand(n):生成n為隨機(jī)矩陣；%ceil(X):取大于X的最大整數(shù)（X為數(shù)值時）或取大于X的最大整數(shù)矩陣（X為矩陣)。function RI=randRI(N)for n=1:10T=10000。 %T值可取10000~30000；一般取10000即可，T越大，精度不一定越高；A=[1 2 3 4 5 6 7 8 9 1/2 1/3 1/4 1/5 1/6 1/7 1/8 1/9]。X=zeros(n)。if n==1 RI(n)=0。elsefor num=1:T for i=1:n for j=i+1:n B=ceil(17*rand(n))。 %生成1~17的隨機(jī)數(shù)矩陣； X(i,j)=A(B(i,j))。 X(j,i)=1/X(i,j)。 X(i,i)=1。 %生成判斷矩陣X； end end nn=size(X,1)。 b=ones(nn,nn)。 %將X歸一化處理成b矩陣 for i=1:n。 for j=1:n b(i,j)=X(i,j)/sum(X(:,j))。 end end wii=sum(b,2)。 %計(jì)算wi wi=wii/sum(wii)。 %規(guī)范化wi（即層因素的權(quán)重） ww=X*wi。 for i=1:n ww(i)=ww(i)/wi(i)。 end lammax=sum(ww)/n。 %求lamdamax(即lammax) CI(num)=(lammaxn)/(n1)。 %隨機(jī)數(shù)一致性指標(biāo)；end RII(n)=sum(CI)/T。 %平均隨機(jī)數(shù)一致性指標(biāo)；endend %%計(jì)算平均隨機(jī)數(shù)一致性指標(biāo)RI%CR（一致性指標(biāo)CI與同階平均隨機(jī)數(shù)指標(biāo)RI之比）=,則判斷矩陣具有滿意的一致性；%CI=（lamdamaxn）/(n1)。%注意：因T較大且程序多為for語句，使循環(huán)次數(shù)較多，程序運(yùn)行時間可能較長，如有需要可取較小的T值；%和積法計(jì)算RI值；%rand(n):生成n為隨機(jī)矩陣；%ceil(X):取大于X的最大整數(shù)（X為數(shù)值時）或取大于X的最大整數(shù)矩陣（X為矩陣)。function RI=randRI(N)for n=1:20T=10000。 %T值可取10000~30000；一般取10000即可，T越大，精度不一定越高；A=[1 2 3 4 5 6 7 8 9 1/2 1/3 1/4 1/5 1/6 1/7 1/8 1/9]。X=zeros(n)。if n==1 RI(n)=0。elsefor num=1:T for i=1:n for j=i+1:n B=ceil(17*rand(n))。 %生成1~17的隨機(jī)數(shù)矩陣； X(i,j)=A(B(i,j))。 X(j,i)=1/X(i,j)。 X(i,i)=1。 %生成判斷矩陣X； end end nn=size(X,1)。 b=ones(nn,nn)。 %將X歸一化處理成b矩陣 for i=1:n。 for j=1:n b(i,j)=X(i,j)/sum(X(:,j))。 end end wii=sum(b,2)。 %計(jì)算wi wi=wii/sum(wii)。 %規(guī)范化wi（即層因素的權(quán)重） ww=X*wi。 for i=1:n ww(i)=ww(i)/wi(i)。 end lammax=sum(ww)/n。 %求lamdamax(即lammax) CI(num)=(lammaxn)/(n1)。 %隨機(jī)數(shù)一致性指標(biāo)；end RII(n)=sum(CI)/T。 %平均隨機(jī)數(shù)一致性指標(biāo)；endendRI=RII(N)。 %輸出N階的平均隨機(jī)數(shù)一致性指標(biāo)值 clear。clcx=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\固定資產(chǎn)投資39。,39。sheet239。,39。c2:o1639。)。y=xlsread(39。D:\0000第一次模擬訓(xùn)練\數(shù)據(jù)處理\固定資產(chǎn)投資39。,39。sheet239。,39。b2:b1639。)。[m,n]=size(x)。%計(jì)算樣本的相關(guān)系數(shù)矩陣Rfor i=1:m for j=1:n xx(i,j)=(x(i,j)mean(x(:,j)))/std(x(:,j))。 % xx為對原數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的矩陣 endend[m,n]=size(y)。%計(jì)算樣本的相關(guān)系數(shù)矩陣Rfor i=1:m for j=1:n yy(i,j)=(y(i,j)mean(y(:,j)))/std(y(:,j))。 % xx為對原數(shù)據(jù)矩陣標(biāo)準(zhǔn)化后的矩陣 endendxx1=[ones(15,1