(完整版)纵横向拉开档次法的MATLAB实现.docxVIP

PAGE1 / NUMPAGES1 简介：本文档为《纵横向拉开档次法的MATLAB实现》，可适用于工程科技领域，主题内容包含globalxystdszxystdxy定义全局变量loadshuju原始数据xystd=zscore(shuju)数据无量纲处理xystdrow,符等。 globalxystdsz xystd x y%定义全局变量 load shuju%原始数据 xystd=zscore(shuju);%数据无量纲处理 [xystdrow,xystdcol]=size(xystd); %区域知识创造能力评价 for tt=1:xystdcol xystdsz{tt}(:,:)=xystd{tt}(:,1:10); %提取区域知识创造能力指标无量纲值 end [xystdszrow,xystdszcol]=size(xystdsz); [xyrow,xycol]=size(xystdsz{1}); w0=zeros(1,xycol); for i=1:xycol w0(1,i)=1/xycol; %优化初始值 end Aeq=[]; beq=[]; lb=zeros(1,xycol);ub=ones(1,xycol);%zeros生成零矩阵；ones生成全1阵。 options =optimset(largescale,off);%优化函数，largescale大规模算法 [w,faval]=fmincon(@YHQU,w0,[],[],Aeq,beq,lb,ub,@fun,options ); %优化求权重；fmincon用来求解非线性多元函数最小值。 wqz1=w./sum(w); %权重归一化 for tt=1:xystdszcol z{tt}(:,1)=xystd{tt}(:,1:10)*wqz1;%求评价值 pxacz(:,tt)=px(z{tt}(:,1)) ;%对评价值排序 end clear w0 w lb ub faval ; clear global xystdsz; %区域知识流动能力评价 for tt=1:xystdszcol xystdsz{tt}(:,:)=xystd{tt}(:,11:16); %提取区域知识流动能力指标无量纲值 end global xystdsz; [xystdszrow,xystdszcol]=size(xystdsz); [xyrow,xycol]=size(xystdsz{1}); w0=zeros(1,xycol); for i=1:xycol w0(1,i)=1/xycol; %优化w初始值 end 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 ); wqz2=w./sum(w); %权重归一化 for tt=1:xystdszcol z{tt}(:,2)=xystd{tt}(:,11:16)*wqz2;%求评价值 pxald(:,tt)=px(z{tt}(:,2)) ;%对评价值排序 end clearw0 w lb ub faval ; clear global xystdsz; %企业技术创新能力评价 for tt=1:xystdszcol xystdsz{tt}(:,:)=xystd{tt}(:,17:28); %提取企业技术创新能力指标无量纲值 end global xystdsz; [xystdszrow,xystdszcol]=size(xystdsz); [xyrow,xycol]=size(xystdsz{1}); w0=zeros(1,xycol); for i=1:xycol w0(1,i)=1/xycol; %优化w初始值 end 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 ); wqz3=w./sum(w); %权重归一化 for tt=1:xystdszcol z{tt}(:,3)=xystd{tt}(:,17:28)*wqz3;%求评价值 pxaqy(:,tt)=px(z{tt}(:,3)) ;%对评价值排序 end clearw0 w lb ub faval ; clear global xystdsz; %创新环境能力评价 for tt=