目前CSDN里面几乎没有一个正确的该模型代码,本人分享代码的测算结果已发表在C刊。
clc
clear
X =[];%投入指标数据
Y = [];%期望产出指标数据
Z = [];%非期望产出指标数据
[m,n]=size(X);%m为投入指标数,n是样本数
s=size(Y,1);%s为期望产出数
q=size(Z,1);%q为非期望产出数
D=1./(m*X');%目标函数的分母
E=1./((s+q)*Y');%等式约束期望产出的分母
F=1./((s+q)*Z');%等式约束非期望产出的分母
LB=zeros(n+m+s+q+1,1);UB=[];
theta=zeros(n,1);
for i=1:n
……
后续代码感兴趣的,可以联系留言或者私聊