## LMI ConstraintsThere are a number of constraints arising from different requirements that needs to be computed as part of this algorithm: The input, , is inside the input polytope The derivative at the vertives , is pointing in/outwards wrt. a facet. All derivatives in the same direction.
Several of these LMIs require knowning the normal to a facet and this
can be computed as the null-space of a basis for the co-dimensional 1
space containing the facet.
It is also important to know if the computed normal points inwards or
outwards with respect to the polytope the facet is part of. This can
be computed by knowing a point inside the polytope eg. a vertex that
is not part of the facet in question.
where is negative if is pointing inwards. The outwards pointing normal can then be computed as . The outward pointing normal for facet will hereafter be denoted . ## Input inside the polytopeFor the input, , to be inside the polytope must be inside each of the polytopes facets. Computing a vector from a point on the facet to and then finding
that vectors inner product with the facets normal, yields a
scalar which sign determines if the point is on the inside or outside
of the facet.
## Derivatives pointing inwardsFor the derivative to point inwards the inner product between the derivative and the outwards normal of the facet must be negative.
## All derivatives in the same directionFor all the derivatives to be in the same direction there must exist
some vector , such that.
This is not an LMI as and are both variable, but for fixed , or the constraint can be tested. It is according to [Habets06a] possible to compute a set of LMIs testing this constraint, but it requires knowledge of the vertices of . It might seem straigth forward to obtain the vertices as we know the input polytope, but the input at each vertice is further constrained by the requirement to block certain facets. There exists algorithms which should be able to compute these vertices, but they are not immidiatly accesible from matlab. The toolbox uses both normals of all facets as testvectors, , along with any user defined testvectors specified in the configuration structure. |