Rev | Line | |
---|
[37] | 1 | function test_complexsdp |
---|
| 2 | |
---|
| 3 | i=sqrt(-1); |
---|
| 4 | P = [4 1+2*i 3-i;1-2*i 3.5 0.8+2.3*i;3+i 0.8-2.3*i 4]; |
---|
| 5 | Zmanual = toeplitz([4.2827,0.8079+1.7342*sqrt(-1) 2.5574-0.7938*i]) |
---|
| 6 | |
---|
| 7 | Z = sdpvar(3,3,'toeplitz','complex'); |
---|
| 8 | Z = Z-sqrt(-1)*diag(imag(diag(Z))); |
---|
| 9 | |
---|
| 10 | t = sdpvar(1,1); |
---|
| 11 | e = Z(:)-P(:); |
---|
| 12 | F = set(Z > 0); |
---|
| 13 | F = F+set('||e||<t'); |
---|
| 14 | sol = solvesdp(F,t) |
---|
| 15 | mbg_asserttrue(sol.problem == 0); |
---|
| 16 | mbg_asserttrue(norm(Zmanual-double(Z)) < 1e-4) |
---|
| 17 | |
---|
| 18 | sol = solvesdp(set(Z > 0),norm(e)) |
---|
| 19 | mbg_asserttrue(sol.problem == 0); |
---|
| 20 | mbg_asserttrue(norm(Zmanual-double(Z)) < 1e-4) |
---|
Note: See
TracBrowser
for help on using the repository browser.