縱橫向拉開檔次法的MATLAB實現(xiàn)

1、clear allglobal xystdsz xystd x y %定義全局變量load data.mat %原始數(shù)據(jù)xystd=bzycl(x,y); %數(shù)據(jù)無量綱處理xystdrow,xystdcol=size(xystd);%-區(qū)域知識創(chuàng)造能力評價-for tt=1:xystdcol xystdsztt(:,:)=xystdtt(:,1:10); %提取區(qū)域知識創(chuàng)造能力指標無量綱值endxystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i

2、)=1/xycol; % 優(yōu)化初始值endAeq=; beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off);w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options ); %優(yōu)化求權重wqz1=w./sum(w); %權重歸一化for tt=1:xystdszcol ztt(:,1)=xystdtt(:,1:10)*wqz1; % 求評價值 pxacz(:,tt)=px(ztt(:,1) ; % 對評價值排序endclear w0 w lb ub faval

3、;clear global xystdsz;%-區(qū)域知識流動能力評價-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,11:16); %提取區(qū)域知識流動能力指標無量綱值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 優(yōu)化w初始值endAeq=; beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options

4、 =optimset(largescale,off);w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options );wqz2=w./sum(w); %權重歸一化for tt=1:xystdszcol ztt(:,2)=xystdtt(:,11:16)*wqz2; % 求評價值 pxald(:,tt)=px(ztt(:,2) ; % 對評價值排序endclear w0 w lb ub faval ;clear global xystdsz;%-企業(yè)技術創(chuàng)新能力評價-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,

5、17:28); %提取企業(yè)技術創(chuàng)新能力指標無量綱值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 優(yōu)化w初始值endAeq=;beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off);w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options );wq

6、z3=w./sum(w); %權重歸一化for tt=1:xystdszcol ztt(:,3)=xystdtt(:,17:28)*wqz3; % 求評價值 pxaqy(:,tt)=px(ztt(:,3) ; % 對評價值排序endclear w0 w lb ub faval ;clear global xystdsz;%-創(chuàng)新環(huán)境能力評價-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,29:35); %提取創(chuàng)新環(huán)境能力指標無量綱值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow

7、,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 優(yōu)化w初始值endAeq=; beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off);w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options );wqz4=w./sum(w); %權重歸一化for tt=1:xystdszcol ztt(:,4)=xystdtt(:,29:35)*wqz4; % 求評價值 pxahj

8、(:,tt)=px(ztt(:,4) ; % 對評價值排序endclear w0 w lb ub faval ;clear global xystdsz;%-創(chuàng)新經(jīng)濟績效能力評價-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,36:42); %提取創(chuàng)新績效能力指標無量綱值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 優(yōu)化w初始值end

9、Aeq=;beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off); w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options );wqz5=w./sum(w); %權重歸一化for tt=1:xystdszcol ztt(:,5)=xystdtt(:,36:42)*wqz5; % 求評價值 pxajx(:,tt)=px(ztt(:,5) ; % 對評價值排序end% -求整體綜合能力評價排序-clear w0 w lb ub faval ;clear globa

10、l xystdsz;xystdsz=z;global xystdsz;xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz 1);w0=zeros(1,xycol);Aeq=;beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off);w,faval=fmincon(YHQU,w0,Aeq,beq,lb,ub,fun,options );wqz6=w./sum(w); %權重歸一化for kk=1:zcol z(:,kk)=zkk*wqz6

11、; % 求評價值 pxazz=px(zz) ; % 對評價值排序end%- zzfxcapcity表示各項能力得分-結果-for i=1:7 for j=1:5 zzfxcapcityj(:,i)=zi(:,j);endendzzfxcapcity6(:,:)=zz;%- pxafxcapcity表示各項能力排序結果-pxafxcapcity1=pxacz;pxafxcapcity2=pxald;pxafxcapcity3=pxaqy;pxafxcapcity4=pxahj;pxafxcapcity5=pxajx;pxafxcapcity6=pxazz; %-優(yōu)化程序 YHQU.m funct

12、ion xysum=YHQU(w,xystdsz)global xystdsz xystdrow,xystdcol=size(xystdsz);xystdrow1,xystdcol1=size(xystdsz1);for i=1:xystdcol for j=1:xystdrow1 xyvalue(j,i)=xystdszi(j,:)*w; endendxymean=mean(mean(xyvalue);for i=1:xystdrow1 %xystdrow1=30 for j=1:xystdcol %xystdcol=7 xyvar(i,j)=(xyvalue(i,j)-xymean).2;

13、 endendxysum=-sum(sum(xyvar);%-排序px.m function pxa=px(gyhjg) % px排序.pxa是排序結果,gyhjg是要進行排序的評價值列向量jga=gyhjg;m,n=size(jga);for col=1:n for i=1:m-1 for j=i+1:m if jga(i,col)0 m=m1; %表示待評價對象數(shù)else if n20 m=m2; endendfor t=1:ycolif n20 xyt(:,1:n2)=yt;endif n10 xmin=min(xt); xmax=max(xt); for i=1:n1 for j=1:m xyt(j,n2+1:n1+n2)=xmax(i)+xmin(i)-xt(j,i); %對逆向指標進行處理； end endendxyflag=all(xyt=zeros(m,n); %判斷是否存在xy(j,i)0xymin=min(xyt);for i=