2022年纵横向拉开档次法的MATLAB实现

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:10); %提取区域

2、知识创造能力指标无量纲值end xystdszrow,xystdszcol=size(xystdsz); xyrow,xycol=size(xystdsz1); 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,

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) ; % 对评价值排序end clear w0 w lb ub faval ; clear global xystdsz; %- 区域知识流动能力评价- for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,11:16); %提取区域知识流动能力

4、指标无量纲值end global xystdsz; xystdszrow,xystdszcol=size(xystdsz); xyrow,xycol=size(xystdsz1); 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 ); 精品学习资料 可选

5、择p d f - - - - - - - - - - - - - - 第 1 页，共 6 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 6 页 - - - - - - - - -wqz2=w./sum(w); % 权重归一化for tt=1:xystdszcol ztt(:,2)=xystdtt(:,11:16)*wqz2; % 求评价值pxald(:,tt)=px(ztt(:,2) ; % 对评价值排序end clear w0 w lb ub faval ; clear global xystds

6、z; %- 企业技术创新能力评价- for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,17:28); %提取企业技术创新能力指标无量纲值end global xystdsz; xystdszrow,xystdszcol=size(xystdsz); xyrow,xycol=size(xystdsz1); 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 =optimse

7、t(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) ; % 对评价值排序end clear w0 w lb ub faval ; clear global xystdsz; %- 创新环境能力评价- for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,29:3

8、5); %提取创新环境能力指标无量纲值end global xystdsz; xystdszrow,xystdszcol=size(xystdsz); xyrow,xycol=size(xystdsz1); 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

9、 ); 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 6 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 6 页 - - - - - - - - -wqz4=w./sum(w); % 权重归一化for tt=1:xystdszcol ztt(:,4)=xystdtt(:,29:35)*wqz4; % 求评价值pxahj(:,tt)=px(ztt(:,4) ; % 对评价值排序end clear w0 w lb ub faval ; clear

10、global xystdsz; %- 创新经济绩效能力评价- for tt=1:xystdszcol xystdsztt(:,:)=xystdtt(:,36:42); %提取创新绩效能力指标无量纲值end global xystdsz; xystdszrow,xystdszcol=size(xystdsz); xyrow,xycol=size(xystdsz1); 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); optio

11、ns =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 global xystdsz; xystdsz=z; global xystdsz; xystdsz

12、row,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; % 求评价值pxazz=px(zz) ; % 对评价值排序end

13、 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 6 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 6 页 - - - - - - - - -%- zzfxcapcity表示各项能力得分-结果 - for i=1:7 for j=1:5 zzfxcapcityj(:,i)=zi(:,j); end end zzfxcapcity6(:,:)=zz; %- pxafxcapcity表示各项能力排序结果- pxafxcapcity1=pxacz;

14、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); for i=1:xystdcol for j=1:xystdrow1 xyvalue(j,i)=xystdszi(j,:)*w; e

15、nd end xymean=mean(mean(xyvalue); for i=1:xystdrow1 %xystdrow1=30 for j=1:xystdcol %xystdcol=7 xyvar(i,j)=(xyvalue(i,j)-xymean).2; end end xysum=-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 j

16、ga(i,col)0 m=m1; %表示待评价对象数else if n20 m=m2; end end for t=1:ycol if n20 xyt(:,1:n2)=yt; end if 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 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 6 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 6 页 - - - - - - - - -end end xyflag=all(xyt=zeros(m,n); %判断是否存在xy(j,i)0 xymin=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 end end xysumt