나를 해결하기 위해 노력하고 다음과 같은 문제에 CVXPY.
문제를 혼합 정수 SDP 으로 인해 PSD 행렬 우리는 설명합니다. 그러나 그에 따라, 이 목록은 그것은 아무도 해법을 처리할 수 있 이러한 문제입니다.
우리가 사용할 수 있습니다 사실 A
은 2x2 매트릭스를 어떻게든 변환이 혼합된 정수 SOCP 문제입니까?
import cvxpy as cp
import matplotlib.pyplot as plt
import numpy as np
np.random.seed(271828)
m = 2; n = 50
x = np.random.randn(m,n)
off = cp.Variable(boolean=True)
A = cp.Variable((2,2), PSD=True)
b = cp.Variable(2)
obj = cp.Maximize(cp.log_det(A))
constraints = [ cp.norm(A@x[:,i] + b) <= 1 + 20*off for i in range(n) ]
constraints += [cp.sum(off) <= 20]
prob = cp.Problem(obj, constraints)
optval = prob.solve(solver='XPRESS', verbose=False) # seems to work, although it's not super accurate
print(f"Optimum value: {optval}")
# plot the ellipse and data
angles = np.linspace(0, 2*np.pi, 200)
rhs = np.row_stack((np.cos(angles) - b.value[0], np.sin(angles) - b.value[1]))
ellipse = np.linalg.solve(A.value, rhs)
plt.scatter(x[0,:], x[1,:])
plt.plot(ellipse[0,:].T, ellipse[1,:].T)
plt.xlabel('Dimension 1'); plt.ylabel('Dimension 2')
plt.title('Minimum Volume Ellipsoid')
plt.show()