In mathematical optimization theory, the linear complementarity problem (LCP) arises frequently in computational mechanics and encompasses the well-known quadratic programming as a special case. It was proposed by Cottle and Dantzig in 1968.
Formulation Given a real matrix M and vector q, the linear complementarity problem LCP(q, M) seeks vectors z and w which satisfy the following constraints:
w , z ⩾ 0 , {\displaystyle w,z\geqslant 0,} (that is, each component of these two vectors is non-negative)
z T w = 0 {\displaystyle z^{T}w=0} or equivalently ∑ i w i z i = 0. {\displaystyle \sum \nolimits _{i}w_{i}z_{i}=0.} This is the complementarity condition, since it implies that, for all i {\displaystyle i} , at most one of w i {\displaystyle w_{i}} and z i {\displaystyle z_{i}} can be positive.
w = M z + q {\displaystyle w=Mz+q}
A sufficient condition for existence and uniqueness of a solution to this problem is that M be symmetric positive-definite. If M is such that LCP(q, M) has a solution for every q, then M is a Q-matrix. If M is such that LCP(q, M) have a unique solution for every q, then M is a P-matrix. Both of these characterizations are sufficient and necessary. The vector w is a slack variable, and so is generally discarded after z is found. As such, the problem can also be formulated as:
M z + q ⩾ 0 {\displaystyle Mz+q\geqslant 0}
z ⩾ 0 {\displaystyle z\geqslant 0}
z T ( M z + q ) = 0 {\displaystyle z^{\mathrm {T} }(Mz+q)=0} (the complementarity condition)
Convex quadratic-minimization: Minimum conditions Finding a solution to the linear complementarity problem is associated with minimizing the quadratic function
f ( z ) = z T ( M z + q ) {\displaystyle f(z)=z^{T}(Mz+q)}
subject to the constraints
M z + q ⩾ 0 {\displaystyle {Mz}+q\geqslant 0}
z ⩾ 0 {\displaystyle z\geqslant 0}
These constraints ensure that f is always non-negative. The minimum of f is 0 at z if and only if z solves the linear complementarity problem. If M is positive definite, any algorithm for solving (strictly) convex QPs can solve the LCP. Specially designed basis-exchange pivoting algorithms, such as Lemke's algorithm and a variant of the simplex algorithm of Dantzig have been used for decades. Besides having polynomial time complexity, interior-point methods are also effective in practice. Also, a quadratic-programming problem stated as minimize f ( x ) = c T x + 1 2 x T Q x {\displaystyle f(x)=c^{T}x+{\tfrac {1}{2}}x^{T}Qx} subject to A x ⩾ b {\displaystyle Ax\geqslant b} as well as x ⩾ 0 {\displaystyle x\geqslant 0} with Q symmetric is the same as solving the LCP with
q = [ c − b ] , M = [ Q − A T A 0 ] {\displaystyle q={\begin{bmatrix}c\\-b\end{bmatrix}},\qquad M={\begin{bmatrix}Q&-A^{T}\\A&0\end{bmatrix}}}
This is because the Karush–Kuhn–Tucker conditions of the QP problem can be written as:
… excerpt ends here. Continue reading the full article.
