人大经济论坛 › 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › MATLAB等数学软件专版 › matlab做DEA,程序出错了,急救啊,谁能帮我解答下

CDA数据分析研究院

商业数据分析与大数据领航教育品牌



经管云课堂

经管/金融/财会/社科/名师公开课



学术培训

Stata 空间计量 SSCI Python

贵宾：通行论坛特权+数据库权限
+案例库+下载特权 VIP：论坛特权+更多下载次数
+ccerdata数据库+更高阅读权限+……

12 3 下一页

发帖

楼主: hmhn421001

5535 21

[问答] matlab做DEA,程序出错了,急救啊,谁能帮我解答下 [推广有奖]

0关注
2粉丝

无

副教授

25%

还不是VIP/贵宾

威望: 0 级
论坛币: 207 个
通用积分: 52.6293
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 498 点
帖子: 241
精华: 0
在线时间: 819 小时
注册时间: 2010-8-6
最后登录: 2024-4-29

楼主

hmhn421001 发表于 2011-12-8 01:43:31 |只看作者 |坛友微信交流群|倒序 |AI写论文

相似文件

换一批

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

程序如下:
clear
X=[264806.95,3998,340244.15;
146823.04,4179,131544.52;
25401.48,929,7862.49;
215305.84,1674,279184.52;
417338.19,6656,240468.29;
475780.83,5210,164915.04;
163304.52,540,45263.70;
255716.76,1876,220412.48;
115885.38,869,108008.70;
232136.41,3325,158022.80;
161578.74,1392,127160.69;
];
Y=[20818.89,406987.07;
14479.76,174703.22;
904.18,12587.72;
3328.15,313829.38;
31047.86,327888.54;
21627.37,253739.73;
2956.68,62213.18;
1964.39,249315.82;
2165.19,117318.03;
9586.89,191067.60;
584.47,140788.65;
];
n=size(X', 1);
m=size(X,1);
s=size(Y,1);
epsilon=10^-10;
f=[zeros(1,n) -epsilon*ones(1,m+s) 1];
A=zeros(1,n+m+s+1);b=0;
LB=zeros(n+m+s+1,1);
UB=[];
LB(n+m+s+1)=-inf;
for i=1:n;
Aeq=[X eye(m) zeros(m,s) -X(:,i)
   Y zeros(s,m) -eye(s) zeros(s,1)];
beq=[zeroS(m,1)
   Y(:,i)];
  w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
w
Lambda=w(1:n,:)
s_minus=w(n+1:n+m,:)
s_plus=w(n+m+1:n+m+s,:)
theta=w(n+m+s+1,:)
运行出现:
Error in ==> CCR_DEA at 36
Aeq=[X eye(m) zeros(m,s) -X(:,i)

扫码加我拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB atlab matla Lab DEA matlab Error 程序

相关帖子

无

使用道具举报

沙发

liuxin9023 发表于 2011-12-8 10:10:09 |只看作者 |坛友微信交流群

把XY转置下

使用道具举报

藤椅

hmhn421001 发表于 2011-12-8 11:07:57 |只看作者 |坛友微信交流群

liuxin9023 发表于 2011-12-8 10:10
把XY转置下

谢谢了,还请问下哪个地方需把XY转置下?由于不是很懂,还麻烦详细指点一下

无

使用道具举报

板凳

hmhn421001 发表于 2011-12-8 11:25:48 |只看作者 |坛友微信交流群

liuxin9023 发表于 2011-12-8 10:10
把XY转置下

clear
X=[264806.95 146823.04 25401.48 215305.84 417338.19 475780.83 163304.52 255716.76  115885.38 232136.41 161578.74;
3998 4179 929    1674 6656 5210 540 1876 869    3325 1392;
340244.15 131544.52 7862.49 279184.52 240468.29 164915.04 45263.7 220412.48 108008.7 158022.8 127160.69;
];
Y=[20818.89  14479.76 904.18 3328.15 31047.86 21627.37 2956.68 1964.39 2165.19 9586.89 584.47;
406987.07 174703.22 12587.72 313829.38 327888.54 253739.73 62213.18 249315.82 117318.03 191067.6 140788.65;];

n=size(X',1);
m=size(X,1);
s=size(Y,1);
epsilon=10^-10;
f=[zeros(1,n) -epsilon*ones(1,m+s) 1];
A=zeros(1,n+m+s+1);b=0;
LB=zeros(n+m+s+1,1);
UB=[];
LB(n+m+s+1)=-inf;
for i=1:n;
Aeq=[X eye(m) zeros(m,s) -X(:,i)
   Y zeros(s,m) -eye(s) zeros(s,1)];
beq=[zeros(m,1) Y(:,i)];
  w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
w
Lambda=w(1:n,:)
s_minus=w(n+1:n+m,:)
s_plus=w(n+m+1:n+m+s,:)
theta=w(n+m+s+1,:)
是需要转置一下。不过还是出现了同样的错误？

无

使用道具举报

报纸

epoh 发表于 2011-12-8 19:04:44 |只看作者 |坛友微信交流群

hmhn421001 发表于 2011-12-8 11:25
clear
X=[264806.95 146823.04 25401.48 215305.84 417338.19 475780.83 163304.52 255716.76 115885.3 ...

% Input matrix
X=[264806.95,3998,340244.15;
146823.04,4179,131544.52;
25401.48,929,7862.49;
215305.84,1674,279184.52;
417338.19,6656,240468.29;
475780.83,5210,164915.04;
163304.52,540,45263.70;
255716.76,1876,220412.48;
115885.38,869,108008.70;
232136.41,3325,158022.80;
161578.74,1392,127160.69];
% Output matrix
Y=[20818.89,406987.07;
14479.76,174703.22;
904.18,12587.72;
3328.15,313829.38;
31047.86,327888.54;
21627.37,253739.73;
2956.68,62213.18;
1964.39,249315.82;
2165.19,117318.03;
9586.89,191067.60;
584.47,140788.65];
%extracts the number of inputs and outputs;
[n,m] = size(X);
[n,s] = size(Y);
epsilon=10^-3;
% Output oriented CCR model;
w = zeros(n,n+m+s+1);
f=-[zeros(1,n) epsilon*ones(1,m+s) 1];
lb=ones(n+s+m+1,1)*(-Inf);
for j=1:n
   Aeq = [-Y', eye(s,s), zeros(s,m), Y(j,:)'; ...
            X', zeros(m,s), eye(m,m), zeros(m,1)];
   beq = [zeros(s,1);X(j,:)'];
   w(j,:) = linprog(f,[],[],Aeq,beq,lb);
end
w

已有 1 人评分	学术水平	热心指数	信用等级	收起理由
zhangtao	+ 5	+ 5	+ 5	热心帮助其他会员

总评分: 学术水平 + 5 热心指数 + 5 信用等级 + 5 查看全部评分

使用道具举报

地板

liuxin9023 发表于 2011-12-9 15:07:11 |只看作者 |坛友微信交流群

clear
X=[264806.95,3998,340244.15;
146823.04,4179,131544.52;
25401.48,929,7862.49;
215305.84,1674,279184.52;
417338.19,6656,240468.29;
475780.83,5210,164915.04;
163304.52,540,45263.70;
255716.76,1876,220412.48;
115885.38,869,108008.70;
232136.41,3325,158022.80;
161578.74,1392,127160.69;
];
Y=[20818.89,406987.07;
14479.76,174703.22;
904.18,12587.72;
3328.15,313829.38;
31047.86,327888.54;
21627.37,253739.73;
2956.68,62213.18;
1964.39,249315.82;
2165.19,117318.03;
9586.89,191067.60;
584.47,140788.65;
];
X=X';
Y=Y';
n=size(X', 1);
m=size(X,1);
s=size(Y,1);
epsilon=10^-10;
f=[zeros(1,n) -epsilon*ones(1,m+s) 1];
A=zeros(1,n+m+s+1);b=0;
LB=zeros(n+m+s+1,1);
UB=[];
LB(n+m+s+1)=-inf;
for i=1:n;
Aeq=[X eye(m) zeros(m,s) -X(:,i)
   Y zeros(s,m) -eye(s) zeros(s,1)];
beq=[zeroS(m,1)
   Y(:,i)];
  w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
w
Lambda=w(1:n,:)
s_minus=w(n+1:n+m,:)
s_plus=w(n+m+1:n+m+s,:)
theta=w(n+m+s+1,:)

已有 2 人评分	学术水平	热心指数	信用等级	收起理由
zhangtao	+ 5	+ 5	+ 5	精彩帖子
hmhn421001	+ 1	+ 1	+ 1	热心帮助其他会员

总评分: 学术水平 + 6 热心指数 + 6 信用等级 + 6 查看全部评分

使用道具举报

7楼

hmhn421001 发表于 2011-12-9 17:00:53 |只看作者 |坛友微信交流群

liuxin9023 发表于 2011-12-9 15:07
clear
X=[264806.95,3998,340244.15;
146823.04,4179,131544.52;

谢谢你了,前面已经把数据转置过来,做好了

无

使用道具举报

8楼

zhangtao 发表于 2011-12-9 17:25:54 |只看作者 |坛友微信交流群

s_minus =

  1.0e+004 *

  Columns 1 through 6

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.3510 0.0000 0.0000 0.0000 0.0000 0.0000

  Columns 7 through 11

0.0000 0.1514 0.0000 0.1790 0.0130
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 2.2154 0.0000 0.0000 0.0000

s_plus =

  1.0e+003 *

  Columns 1 through 6

0.0991 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

  Columns 7 through 11

0.0000 6.7075 0.8214 3.4732 5.7518
0.0000 0.0000 0.0000 0.0000 0.0000

theta =

  Columns 1 through 6

1.0091 1.0000 1.0000 1.0000 1.0000 1.0000

  Columns 7 through 11

1.0000 1.0446 0.9297 0.9162 0.9100

数学好就是要天天学

使用道具举报

9楼

zhangtao 发表于 2011-12-9 17:29:24 |只看作者 |坛友微信交流群

w =

  1.0e+004 *

  Columns 1 through 6

0.0001 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0001 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0001 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0001 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0001 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0001
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.3510 0.0000 0.0000 0.0000 0.0000 0.0000
0.0099 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0001 0.0001 0.0001 0.0001 0.0001 0.0001

  Columns 7 through 11

0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0001 0.0001 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.1514 0.0000 0.1790 0.0130
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 2.2154 0.0000 0.0000 0.0000
0.0000 0.6707 0.0821 0.3473 0.5752
0.0000 0.0000 0.0000 0.0000 0.0000
0.0001 0.0001 0.0001 0.0001 0.0001

数学好就是要天天学

使用道具举报

10楼

zhangtao 发表于 2011-12-9 17:29:45 |只看作者 |坛友微信交流群

w =

  1.0e+006 *

  Columns 1 through 6

0.0001 -0.0002 -0.0009 0.0001 0.0000 0.0001
0.0002 -0.0001 -0.0009 0.0002 -0.0000 0.0002
0.0007 0.0006 -0.0000 0.0020 -0.0009 0.0021
0.0001 -0.0001 -0.0011 0.0001 0.0000 0.0001
0.0001 -0.0002 -0.0013 0.0001 0.0000 0.0001
0.0001 -0.0001 -0.0010 0.0001 0.0000 0.0001
0.0002 -0.0000 -0.0009 0.0005 -0.0002 0.0005
0.0001 -0.0001 -0.0009 0.0001 -0.0000 0.0001
0.0001 -0.0001 -0.0010 0.0002 -0.0001 0.0002
0.0002 -0.0001 -0.0011 0.0002 -0.0000 0.0002
0.0001 -0.0001 -0.0014 0.0002 -0.0001 0.0003

  Columns 7 through 12

0.0001 -0.0000 -0.0001 -0.0001 -0.0001 1.0824
0.0001 -0.0000 -0.0003 -0.0003 -0.0002 1.1131
-0.0003 0.0004 -0.0042 -0.0028 -0.0014 1.4276
0.0001 -0.0000 -0.0002 -0.0002 -0.0002 1.2666
0.0001 -0.0000 -0.0001 -0.0002 -0.0002 1.2082
0.0001 -0.0000 -0.0002 -0.0002 -0.0002 0.9979
0.0001 0.0001 -0.0009 -0.0006 -0.0004 1.3421
0.0001 -0.0000 -0.0002 -0.0002 -0.0002 1.0418
0.0001 0.0000 -0.0003 -0.0003 -0.0002 1.0934
0.0001 -0.0000 -0.0003 -0.0003 -0.0003 1.3107
0.0001 0.0000 -0.0004 -0.0003 -0.0002 1.2825

  Columns 13 through 17

1.1043 0.9693 1.0070 1.1040 0.0000
1.2941 1.7714 1.1146 1.2917 0.0000
1.0401 1.6688 1.0135 1.0381 0.0078
1.4352 1.1737 1.1392 1.4368 0.0000
1.3535 1.0801 1.2664 1.3538 0.0000
1.2725 1.0665 1.0202 1.2741 0.0000
1.3025 1.2773 1.1998 1.3026 0.0003
1.7167 1.1435 1.0171 1.7222 0.0000
1.1504 1.5029 1.1154 1.1492 0.0001
1.8822 1.9726 1.2494 1.8766 0.0000
1.3567 1.8439 1.3853 1.3547 0.0000

数学好就是要天天学

使用道具举报

返回列表

12 3 下一页

发帖

本版微信群

加好友,备注cda
拉您进交流群

手机版 |

意见反馈 |

帮助 |

新手入门 |

用户手册 |

友情链接 |

如有投资本站、合作意向或投放广告，请联系：13661292478（刘老师）

联系客服

邮箱：service@pinggu.org 投诉或不良信息处理：（010-68466864）

京ICP备16021002-2号京B2-20170662号京公网安备 11010802022788号论坛法律顾问：王进律师知识产权保护声明免责及隐私声明

[问答] matlab做DEA,程序出错了,急救啊,谁能帮我解答下 [推广有奖]

经管之家送您一份

经管之家联合CDA

感谢您参与论坛问题回答

扫码加我 拉你入群

相关帖子

本版微信群

扫码加我拉你入群