source:
proiecte/pmake3d/make3d_original/Make3dSingleImageStanford_version0.1/third_party/opt/yalmip/dev/tests-mbg/sos/test_sos_matrix_1.m
@
37
Last change on this file since 37 was 37, checked in by , 14 years ago | |
---|---|
File size: 410 bytes |
Line | |
---|---|
1 | function massive |
2 | |
3 | sdpvar x y |
4 | P = [1+x^2 -x+y+x^2;-x+y+x^2 2*x^2-2*x*y+y^2]; |
5 | m = size(P,1); |
6 | v = monolist([x y],degree(P)/2); |
7 | Q = sdpvar(length(v)*m); |
8 | R = kron(eye(m),v)'*Q*kron(eye(m),v)-P; |
9 | s = coefficients(R(find(triu(R))),[x y]); |
10 | sol = solvesdp(set(Q > 0) + set(s==0)); |
11 | diff = (clean(P - kron(eye(m),v)'*double(Q)*kron(eye(m),v),1e-6)); |
12 | |
13 | mbg_asserttolequal(sol.problem,0); |
14 | mbg_asserttolequal(diff,[0 0;0 0]); |
Note: See TracBrowser
for help on using the repository browser.