根据参数编写DEA-CCR代码 matlab

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

根据参数编写DEA-CCR代码 matlab

这是编写代码最简单的方式，用参数的数量编写，在这就免费分享给大家，希望能帮助到你们！

有数据，有结果，需要的拿去！数据在附件里面。

用心编写即可，后面我会把产出导向CCR也分享给大家，还有投入、产出导向BCC模型的代码！

数据(caj储军飞-交叉效率表2.1)

[nm]=size(X);

[ns]=size(Y);

A=[Y -X];

b=[zeros(n,1)];

for i=1:n;% i就是k

f=[-Y(i,:) zeros(1,m)];%返回矩阵Y的第i行行向量

Aeq=[ zeros(1,s) X(i,:)];

beq=1;

LB=zeros(1,s+m)';

w=linprog(f,A,b,Aeq,beq,LB);

opti=-f*w;

aa(i)=opti;

end

aa

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MATLAB matla atlab CCR Mat

代码第一行和第二行字母之间有空格[n m]=size(X);

[n s]=size(Y);

thanks for sharing

# 设置线性规划问题的约束
A = np.hstack((Y, -X))
b = np.zeros(n)

# 循环计算每个决策单元的效率
for i in range(n):
    f = np.hstack((-Y[i], np.zeros(m)))
    Aeq = np.hstack((np.zeros(s), X[i]))
    beq = 1
    LB = [(0, None) for _ in range(s + m)]
    result = linprog(f, A_ub=A, b_ub=b, A_eq=Aeq.reshape(1, -1), b_eq=[beq], bounds=LB, method='highs')
    aa[i] = -f @ result.x

感谢作者，给大家改写成一个python版本的