Matlab
例子：
1.CCR 乘数形式
clear
load(‘biao.mat’);
X =[‘input1’; ‘input2’]; % 用户输入多指标输入矩阵X。行是输入变量，列是DMU
Y =[‘output1’; ‘output2’]; % 用户输入多指标输出矩阵Y。行是输出变量，列是DMU
n=size(X’,1); m=size(X,1); s=size(Y,1); %n 返回 DMU 的数目，m、s 分别为输入、输出指标数
A=[-X’ Y’]; %A 是不等式约束的系数矩阵，转置后行是 DMU，列是输入输出变量（符合模型
的式子）根据 MATLAB 解线性规划的标准形式，原式要变为：
u
′
y
―
w
′
x
≤
0
，因此 A 为[-X’
Y’]
b=zeros(n,1); %b 是不等约束式子的右边
LB=zeros(m+s,1);UB=[]; %%LB 和 UB 分别是变量的下界和上界，这里定义下界为 0，无上界
for i=1:n; %从第一个 DMU 开始到第 n 个 DMU，各 DMU 的目标函数系数和等式约束不一样，
因此可以用 for 循环结构。
f=[zeros(1,m) -Y(:,i)'];% f 是目标函数的系数向量，在这里目标函数ω的系数为 0，μ的系
数为 yi（即第 i 个被评价 DMU），因为要求最大值，而 matlab 的标准形式是最小值，因此要
在 Y(:,i)'前面加个负号“-”
Aeq=[X(:,i)' zeros(1,s)];beq=1;%Aeq 是等式约束的系数矩阵，beq 是等式约束的右边
w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB); %解线性规划，得 DMUi 的最佳权向量 wi，第 i 列为
DMUi 的权向量，其中 1：m 为 X 的权向量，m+1：m+s 为 Y 的权向量
E(i, i)=Y(:,i)'*w(m+1:m+s,i); %求出 DMUi 的相对效率值 Eii
end
w %输出最佳权向量
E %输出相对效率值 Eii
omega=w(1:m,:) %输出投入权向量 w
mu=w(m+1:m+s,:) %输出产出权向量 m
