6533b85dfe1ef96bd12be91b
RESEARCH PRODUCT
Determinantal sets, singularities and application to optimal control in medical imagery
Bernard BonnardThibaut VerronMohab Safey El DinJean-charles FaugèreAlain Jacquemardsubject
FOS: Computer and information sciences[INFO.INFO-SC]Computer Science [cs]/Symbolic Computation [cs.SC]Computer Science - Symbolic Computation0209 industrial biotechnologyPolynomialRank (linear algebra)010102 general mathematicsBoundary (topology)Field (mathematics)02 engineering and technologySymbolic Computation (cs.SC)Optimal control01 natural sciencesPolynomial system solvingReal geometryPolynomial matrix[ INFO.INFO-SC ] Computer Science [cs]/Symbolic Computation [cs.SC]Set (abstract data type)Matrix (mathematics)020901 industrial engineering & automationApplications0101 mathematicsAlgorithmMathematicsdescription
International audience; Control theory has recently been involved in the field of nuclear magnetic resonance imagery. The goal is to control the magnetic field optimally in order to improve the contrast between two biological matters on the pictures. Geometric optimal control leads us here to analyze mero-morphic vector fields depending upon physical parameters , and having their singularities defined by a deter-minantal variety. The involved matrix has polynomial entries with respect to both the state variables and the parameters. Taking into account the physical constraints of the problem, one needs to classify, with respect to the parameters, the number of real singularities lying in some prescribed semi-algebraic set. We develop a dedicated algorithm for real root classification of the singularities of the rank defects of a polynomial matrix, cut with a given semi-algebraic set. The algorithm works under some genericity assumptions which are easy to check. These assumptions are not so restrictive and are satisfied in the aforementioned application. As more general strategies for real root classification do, our algorithm needs to compute the critical loci of some maps, intersections with the boundary of the semi-algebraic domain, etc. In order to compute these objects, the determinantal structure is exploited through a stratifi-cation by the rank of the polynomial matrix. This speeds up the computations by a factor 100. Furthermore, our implementation is able to solve the application in medical imagery, which was out of reach of more general algorithms for real root classification. For instance, computational results show that the contrast problem where one of the matters is water is partitioned into three distinct classes.
year | journal | country | edition | language |
---|---|---|---|---|
2016-07-20 |