Hello, I am working on a small SDP program. From what I understand, this is what is called a matrix pencil because I only have one decision variable. Unfortunately, I am receiving the infeasible status.

I have a feeling that my problem may be impossible, but it is not very clear because of the simple formulation. The feasibility problem is

It is â€śstructurallyâ€ť infeasible with those input values.

First of all, the diagonal needs to be nonnegative (actually positive, if the SDP constraint is strict, as in the image). But b22 is negative, although it is of sufficiently small magnitude that you might get away with it. But just looking at the lower right 2 by 2 submatrix, which consists entirely of input data: in order for the problem to be feasible, require b22 >= 0, b33 >= 0, and b22*b33 >= b23^2. which for your sample inputs, means that not only must b22 be >= 0, but if b23 and b33are taken as â€śgivenâ€ť, would need b22 >= 3.1462 to get that block to be psd, which is necessary, but not sufficient for the whole matrix to be psd.

Actually, the following program shows that with b11, b12, b23, b33 given as in your example, the minimum value of b22 to make the problem feasible is the cvx_optval of 314.616.

Note that you should never use cvx_precision best Just stick with the default, no matter what solver you use. cvx_precision is a well-intentioned , but ill-advised option.

Another way of looking at it is determining how far R is from being psd, using the most â€śfavorableâ€ť value of x.

b11 = 26.2491;
b12 = -7.39239;
b22 = -3.12284e-12;
b23 = -7.30907;
b33 = 0.169802;
cvx_begin sdp
% define variables
variables x t
% define R
R = [b11 b12 0;b12 b22 b23;0 b23 b33] + x*[2 1 1;1 0 0;1 0 0];
% goal
minimize t
% constraints
subject to:
R + t*eye(3) >= 0
cvx_end

has optimal t = cvx_optval = 7.22466, which shows that the maximum achievable minimum eigenvalue of R is - 7.22466, so that is how far R is from being psd with the most favorable value of x. Or put another way if a scalar multiple of the identity matrix is added to R, that multiple would have to be at least 7.22466 in order to make the problem feasible.

On a different note, what @Mark_L_Stone just did here has intrigued me for my current problem. I noticed that Stephen Boydâ€™s book Linear Matrix Inequalities in System and Control Theory references EVPs. Do you recommend other resources to learn about EVPs and GEVPs?

I donâ€™t think it has any material directly on generalized eigenvalue problems, but if you can make it through Ben-Tal and Nemirovski â€śLECTURES ON MODERN CONVEX OPTIMIZATIONâ€ť, you should be mentally toughened up to make the literature on generalized eigenvalue problems seem not so difficult.