Line | |
---|
1 | function fail = starzak |
---|
2 | |
---|
3 | N = 15; % Number of balls |
---|
4 | n = 6; % ...in R^n |
---|
5 | |
---|
6 | % Generate som balls |
---|
7 | randn('seed',123456); |
---|
8 | for i = 1:N |
---|
9 | P = randn(n);P = P*P'+eye(n)*0.1; |
---|
10 | Pi{i}=P; |
---|
11 | end |
---|
12 | |
---|
13 | % Define LMIs from S-procedure |
---|
14 | tau = sdpvar(N,1); |
---|
15 | P = sdpvar(n,n); |
---|
16 | G = lmi('P>0'); |
---|
17 | F = lmi(P > 1e3*eye(n)*sqrt(eps)); % Numericl reasons |
---|
18 | for i = 1:N |
---|
19 | F = F+lmi(tau(i)*Pi{i} > P); |
---|
20 | end |
---|
21 | F = F+lmi(tau<1); |
---|
22 | |
---|
23 | % MAXDET |
---|
24 | sol = solvesdp(F,-logdet(P),sdpsettings('verbose',0)); |
---|
25 | fail=getfail(sol.problem,-log(det(double(P))),10.8514,checkset(F)); |
---|
Note: See
TracBrowser
for help on using the repository browser.