纵横向拉开档次法的MATLAB实现(共6页).doc

《纵横向拉开档次法的MATLAB实现(共6页).doc》由会员分享，可在线阅读，更多相关《纵横向拉开档次法的MATLAB实现(共6页).doc（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上简介：本文档为纵横向拉开档次法的MATLAB实现，可适用于工程科技领域，主题内容包含globalxystdszxystdxy定义全局变量loadshuju原始数据xystd=zscore(shuju)数据无量纲处理xystdrow,符等。global xystdsz xystd x y %定义全局变量load shuju %原始数据xystd= zscore (shuju); %数据无量纲处理 xystdrow,xystdcol=size(xystd);%-区域知识创造能力评价-for tt=1:xystdcol xystdsztt(:,:)=xystdtt(:,1:

2、10); %提取区域知识创造能力指标无量纲值endxystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 优化初始值endAeq=; beq=;lb=zeros(1,xycol);ub=ones(1,xycol); %zeros生成零矩阵；ones生成全1阵。options =optimset(largescale,off); %优化函数，largescale大规模算法w,faval=fmincon(YHQU,w0,Aeq,

3、beq,lb,ub,fun,options ); %优化求权重；fmincon用来求解非线性多元函数最小值。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 ;clear global xystdsz;%-区域知识流动能力评价-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,11:16); %提取区域知识流动能力指标无量纲值e

4、ndglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 优化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 );wqz2=w./sum(w); %权重归一化for t

5、t=1:xystdszcol ztt(:,2)=xystdtt(:,11:16)*wqz2; % 求评价值 pxald(:,tt)=px(ztt(:,2) ; % 对评价值排序endclear w0 w lb ub faval ;clear global xystdsz;%-企业技术创新能力评价-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,17:28); %提取企业技术创新能力指标无量纲值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1)

6、;w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 优化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 );wqz3=w./sum(w); %权重归一化for tt=1:xystdszcol ztt(:,3)=xystdtt(:,17:28)*wqz3; % 求评价值 pxaqy(:,tt)=px(ztt(:,3) ; %

7、 对评价值排序endclear w0 w lb ub faval ;clear global xystdsz;%-创新环境能力评价-for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,29:35); %提取创新环境能力指标无量纲值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 优化w初始值endAeq=; beq=;lb=zeros(1,xy

8、col);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(:,tt)=px(ztt(:,4) ; % 对评价值排序endclear w0 w lb ub faval ;clear global xystdsz;%-创新经济绩效能力评价-for tt=1:xystdsz

9、col xystdsztt(:,:)=xystdtt(:,36:42); %提取创新绩效能力指标无量纲值endglobal xystdsz; xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz1);w0=zeros(1,xycol);for i=1:xycol w0(1,i)=1/xycol; % 优化w初始值endAeq=;beq=;lb=zeros(1,xycol);ub=ones(1,xycol);options =optimset(largescale,off); w,faval=fmincon(YHQU,w0,Ae

10、q,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 global xystdsz;xystdsz=z;global xystdsz;xystdszrow,xystdszcol=size(xystdsz);xyrow,xycol=size(xystdsz 1);w0=zero

11、s(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; % 求评价值 pxazz=px(zz) ; % 对评价值排序end%- zzfxcapcity表示各项能力得分-结果-for i=1:7 for j=1:5 zzfxcapcityj(:,i)=zi(:,j)

12、;endendzzfxcapcity6(:,:)=zz;%- pxafxcapcity表示各项能力排序结果-pxafxcapcity1=pxacz;pxafxcapcity2=pxald;pxafxcapcity3=pxaqy;pxafxcapcity4=pxahj;pxafxcapcity5=pxajx;pxafxcapcity6=pxazz; %-优化程序 YHQU.m function xysum=YHQU(w,xystdsz)global xystdsz xystdrow,xystdcol=size(xystdsz);xystdrow1,xystdcol1=size(xystdsz1)

13、;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; endendxysum=-sum(sum(xyvar);%-排序px.m function pxa=px(gyhjg) % px排序.pxa是排序结果,gyhjg是要进行排序的评价值列向量jga=gyhjg;m

14、,n=size(jga);for col=1:n for i=1:m-1 for j=i+1:m if jga(i,col)0 m=m1; %表示待评价对象数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=1:n if xyflag(i)=0 for j=1:m xyt(j,i)=xyt(j,i)-xymin(i); %对xy(j,i)0的指标进行平移处理,使所有指标非负 end endendxysumt=sum(xyt); %列求和xystdt=xyt./xysumt(ones(m,1),:); %进行列和等于11归一化处理xystdt=m*xystdt; %避免由于待评价对象m很大时,xystd值很小,所以乘以mend专心-专注-专业