Home > mpt > extras > analysis > mpt_getReachSubset.m

mpt_getReachSubset

PURPOSE ^

MPT_GETREACHSUBSET Computes a subset of P which enters Pfin

SYNOPSIS ^

function [Pret,keptrows,feasible] = mpt_getReachSubset(P,Pfin,A,B,F,G,horizon,Options)

DESCRIPTION ^

MPT_GETREACHSUBSET Computes a subset of P which enters Pfin

 [Pret,keptrows,feasible] = mpt_getReachSubset(P,Pfin,A,B,F,G,horizon,Options)

 ---------------------------------------------------------------------------
 DESCRIPTION
 ---------------------------------------------------------------------------
 this function returns a subset Pret={x | H x <= K} of a polyhedron defined by P
 from which all points enter a final polytope Pfin in N steps. This function is 
 almost identical to "domain". However, here we assume the feedback 
 to be time-varying while "domain" assumes the feedback to be time-invariant.

 ---------------------------------------------------------------------------
 INPUT
 --------------------------------------------------------------------------- 
 P                 - polyhedral partition given by the polytopic array
 Pfin              - target polytope which should be reached in "horizon" steps
 Fi,Gi             - optimal feedback law given by u=Fi{i}*x+Gi{i}
 A, B              - system dynamics are given by x(k+1)=Ax(k)+Bu(k)
 horizon           - number of steps in which the target poly should be reached
 Options.abs_tol   - if abs(x)<tolerance x is considered to be zero
 Options.lpsolver  - LP solver to be used (see help mpt_solveLP for details)

 Note: If Options is missing or some of the fields are not defined, the default
       values from mptOptions will be used

 ---------------------------------------------------------------------------
   OUTPUT:
 ---------------------------------------------------------------------------
 Pret       - subset Pret.H*x<=Pret.K of the input poly P.H*x<=P.K whose
              points enter (P) in "horizon" steps
 keptrows   - is a vector containing the origin of each facet of Pret.H*x<=Pret.K
               1 means that the facet originated from active constraints in P
               0 means that the facet originated from reachability restrictions Pfin
 feasible   - is 1 if the resulting region is not empty

 see also DOMAIN, RANGE

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:
Generated on Thu 30-Mar-2006 10:26:47 by m2html © 2003