全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
11444 30
2018-09-20
用Matlab测算GML指数,总是出现 带有下标的赋值维度不匹配  ,具体错在哪不怎么了解,请高手指教 ,谢谢  以下是程序和excel .          X 投入 , Y 产出 ,B 非期望产出 ,P  时期 , F 单位数。


X=xlsread('E:\matlab\matlab2016\bin\L.xlsx','1','C2:D151');

Y=xlsread('E:\matlab\matlab2016\bin\L.xlsx','1','E2:E151');

B=xlsread('E:\matlab\matlab2016\bin\L.xlsx','1','F2:F151');

F=30;P=5;

i=size(X,2);

a=zeros(F,1);k=zeros(P*F,1);

f1=[a;-1];f2=[k;-1];

g=zeros(i,1);

lb1=zeros(F+1,1);lb1(F+1,1)=-Inf; lb2=zeros(P*F+1,1);lb2(P*F+1,1)=-Inf;ub=[];

for n=1:P;

for m=1+n*F-F:n*F;

d=Y(m,:);

c=Y(1+n*F-F:n*F,:);

z=X(1+n*F-F:n*F,:);

h=[z',g];

e=[-c', d'];

A=[e;h];

j=X(m,:);

b=[-d';j'];

l=B(1+n*F-F:n*F,:);

o=B(m,:);

Aeq=[l',o'];

beq=[o'];

q(:,m)=linprog(f1,A,b,Aeq,beq,lb1,ub);

end

end

for m=1:F*P;

d= Y(m,:);

c=Y(1:F*P,:);

z=X(1:F*P,:);

h=[z',g];

e=[-c',d'];

A=[e;h];

j=X(m,:);

b=[-d';j'];

l=B(1:F*P,:);

o=B(m,:);

Aeq=[l',o'];

beq=[o'];

p(:,m)=linprog(f2,A,b,Aeq,beq,lb2,ub);

end

for m=1:P*F-F;

ec(:,m)=(q(F+1,m)+1)/(q(F+1,m+F)+1);

bpc(:,m)= ((p(P*F+1,m)+1)* (q(F+1,m+F)+1))/((p(P*F+1,m+F)+1)* (q(F+1,m)+1));

gml(:,m)=(p(P*F+1,m)+1)/ (p(P*F+1,m+F)+1);

end

GML=reshape(gml,F,P-1)

EC=reshape(ec,F,P-1)

BPC=reshape(bpc,F,P-1)


1537441248(1).png




二维码

扫码加我 拉你入群

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

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

全部回复
2018-11-2 10:46:49
你好,我最近也打算用MATLAB去测算GML,请问代码问题您解决了没有?
二维码

扫码加我 拉你入群

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

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

2018-11-16 22:36:58
ZwRichard丶 发表于 2018-11-2 10:46
你好,我最近也打算用MATLAB去测算GML,请问代码问题您解决了没有?
没有啊 还不知代码是否正确
二维码

扫码加我 拉你入群

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

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

2018-12-22 21:55:32
您好,请问现在解决了吗?现在写论文也遇到相似的问题了,求助。。。
二维码

扫码加我 拉你入群

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

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

2019-1-13 09:59:06
请问楼主这个问题解决了吗?求助
二维码

扫码加我 拉你入群

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

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

2019-3-18 14:54:35
我用你的这个代码是可以做出来的,这个代码准确吗
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群