I have got a bunch of linear constraints for a positive semidefinite matrix rho like Tr(rho*Gamma{i})=gamma(i), where Gamma is a set of Povm operators, and gamma is the corresponding array of expectation values. I want to get a rho which is close enough to that which satisfies these constraints. And here y is the benchmark
of closeness.
function [rho,y]=closest(dim,Gamma,gamma)
n=size(Gamma,1);
cvx_begin sdp quiet
cvx_precision best
variable rho(dim,dim) hermitian
variable y=(n,1)
minimize norm(y,1)
trace(rho)==1
for i=1:n
abs(trace(rho*Gamma{i})-gamma(i))<=y(i)
end
cvx_end
end
Error using variable (line 57)
Invalid structure specification: y=(n,1)
Error in closest (line 6)
variable rho(dim,dim) hermitian y=(n,1)
for i=1:n
abs(trace(rho*Gamma{n})-gamma(n))<=y(n)
end
Do you mean to have
abs(trace(rho*Gamma{i})-gamma(i))<=y(i)
++++++
Also, I recommend you not use quiet until you are confident your program is working correctly and you trust the results. If quiet is not used, all solver and CVX output will be displayed on the screen; whereas if quiet is used, none of it will be.
Thank you very much for your nice reply. Yeah, I made a mistake for the index when typing. But it still doesn’t work even though I revise the indexing and add up another variable claim. I don’t know why.
You have not provided a mathematical description of the problem you wish to solve. You have not provided the input data. You have not said what you mean by “it doesn’t work”. You have not shown us the results of running it, and why that constitutes not working. Until you provide some of these items, you are unlikely to get useful assistance from readers of the forum, other than my guess in the next paragraph.
The thread title refers to SDP, but your program has no SDP constraints. Do you require rho to be semidefinite? If so, you must declare it as such, or add an explicit semidefinite constraint. The easiest way to do so is to change variable rho(dim,dim) hermitian
to variable rho(dim,dim) hermitian semidefinite
Edit: I see that your first post mentions “positive semidefinite matrix rho”. So it appears my guess above is correct. Either declare rho hermitian semidefinite or hermitian semidefinite, depending on whether you want to allow complex elements in rho.
I recommend you carefully read or re-read the entirety of the CVX Users’ Guide http://cvxr.com/cvx/doc/