你好,我也学习了你编的sbm程序,但是感觉有疑惑的地方。我的理解是只有Aeq,beq没有A和b;第一段为什么要分开求linprog?对于第一部分我的理解是:
clc,clear all;
output=[1;2;3;4;5;6;7;8;9];
input=[7 3;8 1;4 3;2 4;4 2;11 1;10 1;10 3];
[n,m]=size(output)
[j,k]=size(input)
for i=1:j
f1=[zeros(j,1);-1/(2*k);-1/(2*k)];
Aeq1=[input' [1 0;0 1];ones(j,1)' 0 0];
beq1=[input(i,:)';1];
lb1=[zeros(j+k,1)]
[x,fival]=linprog(f1,[],[],Aeq1,beq1,lb1)
end;
|