%30个省份综合政策投入评价
pzc=［］；%手动在工作区内输入数据
［pzcb，MU，SIGMA］=zscore（pzc）；
a=zeros（size（pzcb（1∶30，：）'*pzcb（1∶30，：）））；
for i=1∶30∶481
t=pzcb（i∶i+29，：）；
a=a+
*t；
end%计算H，H=a
［x，y］=eig（a）；%特征值和特征向量
r=diag（y）；%或者写成r=abs（sum（y））
n=find（r==max（r））；
max_y=y（n，n）；%最大特征值
max_x=x（：，n）；%最大特征值对应的特征向量
zc_q=bsxfun（@rdivide，max_x，sum（max_x））；
pzcdf1=pzcb*zc_q；
%30个省份综合科技投入评价
pkj=［］；%输入数据
［pkjb，MU，SIGMA］=zscore（pkj）；
b=zeros（size（pkjb（1∶30，：）'*pkjb（1∶30，：）））；
for i=1：30∶481
t=pkjb（i∶i+29，：）；
b=b+
*t；
end%计算H，H=b
［x，y］=eig（b）；%特征值和特征向量
r=diag（y）；%或者写成r=abs（sum（y））
n=find（r==ma x（r））；
max_y=y（n，n）；%最大特征值
max_x=x（：，n）；%最大特征值对应的特征向量
kj_q=bsxfun（@rdivide，max_x，sum（max_x））；
pkjdfl=pkjb*kj_q；
%30个省份劳动者积极性投入评价(www.xing528.com)
pldz=［］；%输入数据
［pldzb，MU，SIGMA］=zscore（pldz）；
c=zeros（size（pldzb（1∶31，：）'*pldzb（1∶31，：）））；
for i=1∶30∶481
t=pldzb（i∶i+29，：）；
c=c+
*t；
end%计算H，H=c
［x，y］=eig（c）；%特征值和特征向量
r=diag（y）；%或者写成r=abs（sum（y））
n=find（r==ma x（r））；
max_y=y（n，n）；%最大特征值
max_x=x（：，n）；%最大特征值对应的特征向量
ldz_q=bsxfun（@rdivide，max_x，sum（max_x））；
pldzdfl=pldzb*ldz_q；
%30个省份软投入综合评价
%标准化处理
pp=［］；%输入数据
［Z，MU，SIGMA］=zscore（pp）；
d=zeros（size（Z（1∶30，：）'*Z（1∶30，：）））；
for i=1∶30∶481
t=Z（i∶i+29，：）；
d=d+
*t；
end%计算H，H=a
［x，y］=eig（d）；%特征值和特征向量
r=diag（y）；%或者写成r=abs（sum（y））
n=find（r==max（r））；
max_y=y（n，n）；%最大特征值
max_x=x（：，n）；%最大特征值对应的特征向量
pp_q=bsxfun（@rdivide，max_x，sum（max_x））；
pdf l=Z*pp_q；
免责声明：以上内容源自网络，版权归原作者所有，如有侵犯您的原创版权请告知，我们将尽快删除相关内容。
