clc
清晰的
清除所有
want_to_do =“展示”;
表= '是的';
save_results = 0;
% % % % % % % % % % % % % % % % % % %新添加% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %加载数据% % % % % % % % % % % % % % % % % % % %
data = xlsread(“文件C: \ \戴尔\用户桌面\ \ Workr \ DEA_prac \ DEA_prac \ K011.xls ');
id =数据(:1);年=数据(:,2);Y =数据(:,3);
X_H =数据(:,6);xk =数据(:7);X_L =数据(:,5);
% % % % % % % % % % % % % % % %平方和变量
Y_sq = (y ^ 2)之和;
X_Hsq =总和(X_H。^ 2);
X_Ksq =总和(xk。^ 2);
X_Lsq =总和(X_L。^ 2);
% % % % % % % % % % % % % % %变量变量正常化% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
Y_norm = y / Y_sq;
X_Hnorm = X_H. / X_Hsq;
X_Knorm = X_K. / X_Ksq;
X_Lnorm = X_L. / X_Lsq;
% % % % % % % %输入和输出向量
输入= [X_Hnorm、X_Knorm X_Lnorm];
输出= [Y_norm];
data1 =(输入、输出);
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
n =大小(数据,1);
t = (min(数据(:,2)):1:max(数据(:,2)))';
no_time =大小(t, 1);
N =大小(data1、1);
时间=数据(:,2);
%产出距离函数估计了DEA (O ' donnell, 2011)
函数(OTE、exitflag count_OTE] = outputTE_DEA_unban(输出、输入、N, t,时间,1)
exitflag = [];
选择= optimset(“显示”,“关闭”);
J =大小(输出,2);
K =大小(输入,2);
投票= 0 (N, 1);
n_t =大小(t, 1);
i = 1: n_t % % % % t
Q = [];
X = [];
j = 1: N
如果时间(j) = = t(我)
Q = [Q;输出(j)):;
X = [X;输入(j)):;
结束
结束
eval([‘问’,num2str(我),' =问;']);
eval ([' X ', num2str(我),' = X; ']);
结束
count_OTE = 0;
i = 1: N
j = 1: n_t % % % % % t
如果t (j) = =时间(我)
eval ([Qm =问,num2str (j), '; '])
eval ([Xm = X, num2str (j), '; '])
结束
结束
Qm = Qm ';
Xm = Xm ';
太=大小(Qm, 2);
% % % % % % % %
f = [1; 0 (mt, 1)];
A =[[输出(我,:)qm]; [0 (K, 1) Xm]];
b = [0 (J, 1);输入(我:)];
% % % % % % % % % % % % % % % % % % % % %
A_eq =[0的(mt)];
b_eq = 1;
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
磅= [0 0 (mt)] ';
乌兰巴托= [];
(para lambda_minus exitflag0] = linprog (f, A、b、A_eq b_eq,磅,乌兰巴托,[],选项);
% Para =[对位;帕拉');
如果exitflag0 = = 1
投票(i) = 1 / lambda_minus;
如果投票(我)> 1
count_OTE = count_OTE + 1;
投票(i) = 1;
结束
其他的
投票(i) = 0.5;
结束
exitflag = [exitflag; exitflag0];
结束
结束
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
我计算技术效率
谁能请告诉上面的代码有什么问题吗
