Line | |
---|
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.