data1=xlsread('C:\Users\Li-PC\Desktop\dem4.xlsx',1)
data2=xlsread('C:\Users\Li-PC\Desktop\dem4.xlsx',2)
data1
data2
X=[data1(1:3,1:90)]
Y=[data1(4,1:90)]
Z=[data1(5,1:90)]
D=[data2(1:90,1:3)]
E=[data2(1:90,4)]
F=[data2(1:90,5)]
n=size(X',1);m=size(X,1);s=size(Y,1);q=size(Z,1);
A=zeros(1,n+m+s+q+1);b=0;
LB=zeros(n+m+s+q+1);UB=[];
LB(n+m+s+q+1)=-Inf;
for i=1:n;
f=[zeros(1,n)-D(i,:) zeros(1,s) zeros(1,q) 1];
Aeq=[Xeye(m) zeros(m,s) zeros(m,q) -X(:,i)
Y zeros(s,m) -eye(s) zeros(s,q) -Y(:,i)
Z zeros(q,m) zeros(q,s) eye(q) -Z(:,i)
zeros(1,n) zeros(1,m) E(i,:) F(i,:) 1];
beq=[zeros(m,1)
zeros(s,1)
zeros(q,1)
1];
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
end
w
lambda=w(1:n,:)
s_minus=w(n+1:n+m,:)
s_plusg=w(n+m+1:n+m+s,:)
s_plusb=w(n+m+s+1:n+m+s+q,:)
theta=w(n+m+s+q+1,:)
求大神解答,谢谢谢谢!


雷达卡


京公网安备 11010802022788号







