楼主: hmhn421001
6495 21

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

  • 0关注
  • 2粉丝

无

已卖：1份资源

副教授

25%

还不是VIP/贵宾

-

0%
威望
0 级
论坛币
207 个
通用积分
52.6293
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
498 点
帖子
241
精华
0
在线时间
827 小时
注册时间
2010-8-6
最后登录
2024-11-23
楼主
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
关键词：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
数学好就是要天天学

返回列表
123下一页
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

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

GMT+8, 2025-12-30 22:00