Rev | Line | |
---|
[37] | 1 | function socp1 |
---|
| 2 | |
---|
| 3 | X = sdpvar(3,3); |
---|
| 4 | x = sdpvar(3,1); |
---|
| 5 | obj = trace(X)+sum(x); |
---|
| 6 | F = set(X>0) + set(cone(x(2:end),1+x(1))) + set(trace(X)==x(1)+2*x(2)+3*x(3)+4)+set(X(1,3)==8); |
---|
| 7 | |
---|
| 8 | sol1 = solvesdp(F,obj); |
---|
| 9 | obj1 = double(obj); |
---|
| 10 | p1 = checkset(F); |
---|
| 11 | |
---|
| 12 | sol2 = solvesdp(F,obj,sdpsettings('dualize',1)); |
---|
| 13 | obj2 = double(obj); |
---|
| 14 | p2 = checkset(F); |
---|
| 15 | |
---|
| 16 | mbg_asserttolequal(sol1.problem,0); |
---|
| 17 | mbg_asserttolequal(sol2.problem,0); |
---|
| 18 | mbg_asserttolequal(obj1,obj2, 1e-5); |
---|
| 19 | mbg_asserttolequal(min(p1),0, 1e-5); |
---|
| 20 | mbg_asserttolequal(min(p2),0, 1e-5); |
---|
Note: See
TracBrowser
for help on using the repository browser.