b_measured = unifrnd(1,5,200,1);
b_targetl = 5-5e-3;
b_targetu = 5+5e-3;
c = ones(100,1);
dmu = 0.025;
a = zeros(200,100);
for m= 1:200
a(m, = exprnd((m-1)*0.025, 1, 100);
end
b = zeros(200, 1);
bl = b+b_targetl;
bu = b+b_targetu;
aa = [ -a; a];
bb = [ -bl; bu ];
[x,y] = linprog(c,aa,bb,[],[],zeros(100,1));
Exiting: One or more of the residuals, duality gap, or total relative error
has stalled:
the primal appears to be infeasible (and the dual unbounded).
(The dual residual < TolFun=1.00e-008.) |