Tools for Control System Design : Stratification of Matrix Pairs and Periodic Riccati Differential Equation Solvers

Modern control theory is today an interdisciplinary area of research. Just as much as this can be problematic, it also provides a rich research environment where practice and theory meet. This Thesis is conducted in the borderline between computing science (numerical analysis) and applied control theory. The design and analysis of a modern control system is a complex problem that requires high qualitative software to accomplish. Ideally, such software should be based on robust methods and numerical stable algorithms that provide quantitative as well as qualitative information. The introduction of the Thesis is dedicated to the underlying control theory and to introduce the reader to the main subjects. Throughout the Thesis, the theory is illustrated with several examples, and similarities and differences between the terminology from mathematics, systems and control theory, and numerical linear algebra are highlighted. The main contributions of the Thesis are structured in two parts, dealing with two mainly unrelated subjects. Part I is devoted to the qualitative information which is provided by the stratification of orbits and bundles of matrices, matrix pencils and system pencils. Before the theory of stratification is established the reader is introduced to different canonical forms which reveal the system characteristics of the model under investigation. A stratification reveals which canonical structures of matrix (system) pencils are near each other in the sense of small perturbations of the data. Fundamental concepts in systems and control, like controllability and observability of linear continuous-time systems, are considered and it is shown how these system characteristics can be investigated using the stratification theory. New results are presented in the form of the cover relations (nearest neighbours) for controllability and observability pairs. Moreover, the permutation matrices which take a matrix pencil in the Kronecker canonical form to the corresponding system pencil in (generalized) Brunovsky canonical form are derived. Two novel algorithms for determining the permutation matrices are provided. Part II deals with numerical methods for solving periodic Riccati differential equations (PRDE:s). The PRDE:s under investigation arise when solving the linear quadratic regulator (LQR) problem for periodic linear time-varying (LTV) systems. These types of (periodic) LQR problems turn up for example in motion planning of underactuated mechanical systems, like a humanoid robot, the Furuta pendulum, and pendulums on carts. The constructions of the nonlinear controllers are based on linear versions found by stabilizing transverse dynamics of the systems along cycles. Three different methods explicitly designed for solving the PRDE are evaluated on both artificial systems and stabilizing problems originating from experimental control systems. The methods are the one-shot generator method and two recently proposed methods: the multi-shot method (two variants) and the SDP method. As these methods use different approaches to solve the PRDE, their numerical behavior and performance are dependent on the nature of the underlying control problem. Such method characteristics are investigated and summarized with respect to different user requirements (the need for accuracy and possible restrictions on the solution time).Modern reglerteknik är idag i högsta grad ett interdisciplinärt forskningsområde. Lika mycket som detta kan vara problematiskt, resulterar det i en stimulerande forskningsmiljö där både praktik och teori knyts samman. Denna avhandling är utförd i gränsområdet mellan datavetenskap (numerisk analys) och tillämpad reglerteknik. Att designa och analysera ett modernt styrsystem är ett komplext problem som erfordrar högkvalitativ mjukvara. Det ideala är att mjukvaran består av robusta metoder och numeriskt stabila algoritmer som kan leverera både kvantitativ och kvalitativ information.Introduktionen till avhandlingen beskriver grundläggande styr- och reglerteori samt ger en introduktion till de huvudsakliga problemställningarna. Genom hela avhandlingen illustreras teori med exempel. Vidare belyses likheter och skillnader i terminologin som används inom matematik, styr- och reglerteori samt numerisk linjär algebra. Avhandlingen är uppdelade i två delar som behandlar två i huvudsak orelaterade problemklasser. Del I ägnas åt den kvalitativa informationen som ges av stratifiering av mångfalder (orbits och bundles) av matriser, matrisknippen och systemknippen. Innan teorin för stratifiering introduceras beskrivs olika kanoniska former, vilka var och en avslöjar olika systemegenskaper hos den undersökta modellen. En stratifiering ger information om bl.a. vilka kanoniska strukturer av matrisknippen (systemknippen) som är nära varandra med avseende på små störningar i datat. Fundamentala koncept i styr- och reglerteori behandlas, så som styrbarhet och observerbarhet av linjära tidskontinuerliga system, och hur dessa systemegenskaper kan undersökas med hjälp av stratifiering. Nya resultat presenteras i form av relationerna för täckande (närmsta grannar) styrbarhets- och observerbarhets-par. Dessutom härleds permutationsmatriserna som tar ett matrisknippe i Kroneckers kanoniska form till motsvarande systemknippe i (generaliserade) Brunovskys kanoniska form. Två algoritmer för att bestämma dessa permutationsmatriser presenteras. Del II avhandlar numeriska metoder för att lösa periodiska Riccati differentialekvationer (PRDE:er). De undersökta PRDE:erna uppkommer när ett linjärt kvadratiskt regulatorproblem för periodiska linjära tidsvariabla (LTV) system löses. Dessa typer av (periodiska) regulatorproblem dyker upp till exempel när man planerar rörelser för understyrda (underactuated) mekaniska system, så som en humanoid (mänsklig) robot, Furuta-pendeln och en vagn med en inverterad (stående) pendel. Konstruktionen av det icke-linjära styrsystemet är baserat på en linjär variant som bestäms via stabilisering av systemets transversella dynamik längs med cirkulära banor. Tre metoder explicit konstruerade för att lösa PRDE:er evalueras på både artificiella system och stabiliseringsproblem av experimentella styrsystem. Metoderna är sk. en- och flerskotts metoder (one-shot, multi-shot) och SDP-metoden. Då dessa metoder använder olika tillvägagångssätt för att lösa en PRDE, beror dess numeriska egenskaper och effektivitet på det underliggande styrproblemet. Sådana metodegenskaper undersöks och sammanfattas med avseende på olika användares behov, t.ex. önskad noggrannhet och tänkbar begränsning i hur lång tid det får ta att hitta en lösning