x=[6行2列的矩阵];
y=[6行25列矩阵];
n=size(x',1);
m=size(x,1);
s=size(y,1);
a=[-x' y'];b=zeros(n,1);
lb=zeros(m+s,1);ub=[];
options=optimset('maxiter',1e4,'tolfun',1e-6);
for i=1:n
aeq=[x(:,i)' zeros(1,s)];beq=1;
f=[zeros(1,m) -y(:,i)'];
w(:,i)=linprog(f,a,b,aeq,beq,lb,ub,[],options);
E(i,i)=y(:,i)'*w(m+1:m+s,i);
for k=1:n
f=[zeros(1,m) y(:,k)'];
aeq=[x(:,k)' zeros(1,s);E(i,i)*x(:,i)' -y(:,i)'];
beq=[1;0];
v=linprog(f,a,b,aeq,beq,lb,ub,[],options);
E(i,k)=(y(:,k)'*v(m+1:m+s))/(x(:,k)'*v(1:m));
end
end
该程序出现错误,请给予改正,您可设置出售帖,1-500币悬赏,请讲诚信
|