In this thesis we study the utility maximization problem for assets whose prices are cointegrated, which arises from the investment practice of convergence trading and its special forms, pairs trading and spread trading. The major theme in the first two chapters of the thesis, is to investigate the assumption of market-neutrality of the optimal convergence trading strategies, which is a ubiquitous assumption taken by practitioners and academics alike. This assumption lacks a theoretical justification and, to the best of our knowledge, the only relevant study is Liu and Timmermann (2013) which implies that the optimal convergence strategies are, in general, not market-neutral. We start by considering a minimalistic pairs-trading scenario with two cointegrated stocks and solve the Merton investment problem with power and logarithmic utilities. We pay special attention to when/if the stochastic control problem is well-posed, which is overlooked in the study done by Liu and Timmermann (2013). In particular, we show that the problem is ill-posed if and only if the agent’s risk-aversion is less than a constant which is an explicit function of the market parameters. This condition, in turn, yields the necessary and sufficient condition for well-posedness of the Merton problem for all possible values of agent’s risk-aversion. The resulting well-posedness condition is surprisingly strict and, in particular, is equivalent to assuming the optimal investment strategy in the stocks to be market-neutral. Furthermore, it is shown that the well-posedness condition is equivalent to applying Novikov’s condition to the market-price of risk, which is a ubiquitous sufficient condition for imposing absence of arbitrage. To the best of our knowledge, these are the only theoretical results for supporting the assumption of market-neutrality of convergence trading strategies. We then generalise the results to the more realistic setting of multiple cointegrated assets, assuming risk factors that effects the asset returns, and general utility functions for investor’s preference. In the process of generalising the bivariate results, we also obtained some well-posedness conditions for matrix Riccati differential equations which are, to the best of our knowledge, new. In the last chapter, we set up and justify a Merton problem that is related to spread-trading with two futures assets and assuming proportional transaction costs. The model possesses three characteristics whose combination makes it different from the existing literature on proportional transaction costs: 1) finite time horizon, 2) Multiple risky assets 3) stochastic opportunity set. We introduce the HJB equation and provide rigorous arguments showing that the corresponding value function is the viscosity solution of the HJB equation. We end the chapter by devising a numerical scheme, based on the penalty method of Forsyth and Vetzal (2002), to approximate the viscosity solution of the HJB equation.

Dans cette thèse, nous avons développé deux approches pour la construction de contrôleurs approchés semi-décentralisés. La thèse est partagée en deux parties, chaque partie décrivant une approche précise. Première Partie: elle traite de l'approximation semi-décentralisée d'un contrôle optimal pour des équations aux dérivées partielles (EDPs) dans un domaine borné. Dans cette partie on présente une méthode de calcul de contrôle optimal pour des systèmes distribués linéaires avec un opérateur d'entré borné ou non borné. Sa construction repose sur le calcul fonctionnel des opérateurs auto-adjoints et sur la formule de Dunford- Schwartz. Elle est conçue pour des architectures de calcul à très fine granularité, avec coordination semi-décentralisée. Enfin, elle est illustrée par des exemples portant en particulier sur la stabilisation interne de la chaleur, la stabilisation des vibrations d'une poutre, la stabilisation des vibrations dans une matrice de micro-cantilevers... Deuxième Partie: elle est consacrée à l'obtention de réalisations d'état, d'opérateurs linéaires solutions de quelques équations opératorielles différentielles linéaires dans des domaines bornés mono-dimensionnels. Nous proposons deux approches dans le cadre de réalisations diffusives. La première utilise des symboles complexes et la seconde des symboles réels sur l'axe réel. Puis, on illustre la théorie et on développe des méthodes numériques pour le contexte d'une application à l'équation de Lyapunov issue de la théorie du contrôle optimal pour l'équation de la chaleur. Un intérêt pratique pour cette approche est le calcul en temps réel sur des processeurs organisés pour une architecture semi-décentralisée.

This work is devoted to investigating the algebraic and analytic integrability of first order polynomial partial differential equations via an understanding of the well-developed area of local and global integrability of polynomial vector fields. In the view of characteristics method, the search of first integrals of the first order partial differential equations P(x,y,z)∂z(x,y) ∂x +Q(x,y,z)∂z(x,y) ∂y = R(x,y,z), (1) is equivalent to the search of first integrals of the system of the ordinary differential equations dx/dt= P(x,y,z), dy/dt= Q(x,y,z), dz/dt= R(x,y,z). (2) The trajectories of (2) will be found by representing these trajectories as the intersection of level surfaces of first integrals of (1). We would like to investigate the integrability of the partial differential equation (1) around a singularity. This is a case where understanding of ordinary differential equations will help understanding of partial differential equations. Clearly, first integrals of the partial differential equation (1), are first integrals of the ordinary differential equations (2). So, if (2) has two first integrals φ1(x,y,z) =C1and φ2(x,y,z) =C2, where C1and C2 are constants, then the general solution of (1) is F(φ1,φ2) = 0, where F is an arbitrary function of φ1and φ2. We choose for our investigation a system with quadratic nonlinearities and such that the axes planes are invariant for the characteristics: this gives three dimensional Lotka– Volterra systems x' =dx/dt= P = x(λ +ax+by+cz), y' =dy/dt= Q = y(µ +dx+ey+ fz), z' =dz/dt= R = z(ν +gx+hy+kz), where λ,µ,ν 6= 0. v Several problems have been investigated in this work such as the study of local integrability and linearizability of three dimensional Lotka–Volterra equations with (λ:µ:ν)–resonance. More precisely, we give a complete set of necessary and sufficient conditions for both integrability and linearizability for three dimensional Lotka-Volterra systems for (1:−1:1), (2:−1:1) and (1:−2:1)–resonance. To prove their sufficiency, we mainly use the method of Darboux with the existence of inverse Jacobi multipliers, and the linearizability of a node in two variables with power-series arguments in the third variable. Also, more general three dimensional system have been investigated and necessary and sufficient conditions are obtained. In another approach, we also consider the applicability of an entirely different method which based on the monodromy method to prove the sufficiency of integrability of these systems. These investigations, in fact, mean that we generalized the classical centre-focus problem in two dimensional vector fields to three dimensional vector fields. In three dimensions, the possible mechanisms underling integrability are more difficult and computationally much harder. We also give a generalization of Singer’s theorem about the existence of Liouvillian first integrals in codimension 1 foliations in Cnas well as to three dimensional vector fields. Finally, we characterize the centres of the quasi-homogeneous planar polynomial differential systems of degree three. We show that at most one limit cycle can bifurcate from the periodic orbits of a centre of a cubic homogeneous polynomial system using the averaging theory of first order.

This thesis deals with development of guidance laws for advanced applications. Two class of guidance problems, namely, impact angle constrained guidance and pulsed guidance for exo-atmospheric engagements, are considered here. Three impact angle constrained guidance schemes are developed using (i) Proportional navigation guidance (PNG), (ii) State Dependent Riccati Equation (SDRE) technique and (iii) geometric concepts, respectively. A collision course based pulsed guidance law is presented for exo-atmospheric interceptors. Proportional Navigation Guidance (PNG) law is the most widely used guidance law because of its ease of implementation and efficiency. However, in its original form, it achieves only a limited set of impact angles. A two stage PNG law is presented for achieving all impact angles against a stationary target. In the first phase of guidance, an orientation PNG command is used. The orientation navigation constant (N ) is a function of the initial engagement geometry and has a lower value (N less than 2). It is proved that following the orientation trajectory, the interceptor can switch to N = 2 and achieve the desired impact angle. Simulations, with a constant speed and with a realistic interceptor model, show successful interception of the target with all desired impact angles. Feedback implementation of the guidance law results in negligible errors in impact angle with uncompensated autopilot delays. The idea of a two-stage PNG law with impact angle constraint is further used to develop a guidance law for intercepting moving targets. Following the orientation trajectory, the interceptor can switch to N = 3 and achieve the desired impact angle. It is proved that the guidance achieves all impact angles in a surface-to-surface engagement scenario with receding and approaching targets, respectively. In a air-to-surface engagement scenario, it is proved that the guidance law achieves all impact angles in a deterministic set. Constant speed and realistic interceptor models are used for simulations. Results show negligible error in impact angle and miss distance for moving targets. The guidance law, in its feedback implementation form, achieves the desired impact angle for interceptors with delay and with a maneuvering target. The impact angle errors are low with negligible errors in miss distance. Next, the impact angle constrained guidance problem against a stationary target is solved as a non-linear regulator problem using the SDRE technique. The interceptor guidance problems are of finite time nature. As the main contribution of this part of the work, we solve a finite time interceptor guidance problem with infinite horizon SDRE formulation by choosing the state weighting matrix as a function of time-to-go. Numerical simulations are carried out both for a constant speed interceptor model and a realistic interceptor model. Simulations for both the models are carried out for various impact angles and firing angles. Robustness of the proposed guidance law with respect to autopilot lag is also verified by simulations. Results obtained show the efficiency of the SDRE approach for impact angle constrained missile guidance. A geometric guidance scheme is proposed for lateral interception of targets in a planar engagement scenario in the absence of line-of-sight rate information. A kill-band is defined for target initial positions capturable by an arc maneuver, followed by a straight line path by the interceptor. Guidance law for capturing targets inside the kill-band is presented and is further modified for targets outside the kill-band. Based on analytical studies on the kill-band, a guidance law is proposed for lateral interception of maneuvering targets. Simulations are carried with for typical low speed engagements. The concept of kill-band provides an inherent robustness to the proposed guidance law with respect to uncompensated system delays and target maneuver. As the final part of the work, an interceptor endgame pulsed guidance law for exoatmospheric engagements is derived by using the notion of collision heading. The proposed guidance law is derived in steps by (i) Obtaining the collision heading based on the collision triangle engagement geometry and then (ii) Computing the width of the pulse fired by the divert thruster to attain the collision heading. It is shown that this strategy is more effective than the existing zero effort miss (ZEM) based guidance laws for intercepting targets with higher heading angles off the nominal head-on collision course. A result on pulse firing sequence is also presented showing that firing pulses in quick succession results in minimum pulse widths and hence minimum control effort for a desired miss distance. Simulations are carried out for various engagement scenarios. Results show better miss-distance and divert thrust performance as compared to the existing ZEM based law.

A promising line of research attempts to bridge the gap between detector and tracker by means of considering jointly optimal parameter settings for both of these subsystems. Along this fruitful path, this thesis study focuses on the problem of detection threshold optimization in a tracker-aware manner so that a feedback from the tracker to the detector is established to maximize the overall system performance. Special emphasis is given to the optimization schemes based on two non-simulation performance prediction (NSPP) methodologies for the probabilistic data association filter (PDAF), namely, the modified Riccati equation (MRE) and the hybrid conditional averaging (HYCA) algorithm. The possible improvements are presented in two domains: Non-maneuvering and maneuvering target tracking. In the first domain, a number of algorithmic and experimental evaluation gaps are identified and newly proposed methods are compared with the existing ones in a unified theoretical and experimental framework. Furthermore, for the MRE based dynamic threshold optimization problem, a closed-form solution is proposed. This solution brings a theoretical lower bound on the operating signal-to-noise ratio (SNR) concerning when the tracking system should be switched to the track before detect (TBD) mode. As the improvements of the second domain, some of the ideas used in the first domain are extended to the maneuvering target tracking case. The primary contribution is made by extending the dynamic optimization schemes applicable to the PDAF to the interacting multiple model probabilistic data association filter (IMM-PDAF). Resulting in an online feedback from the filter to the detector, this extension makes the tracking system robust against track losses under low SNR values.

We present a new method for the computation of low rank approximations to the solution of large, sparse, stable Lyapunov equations. It is based on a generalization of the classical Smith method and profits by the usual low rank property of the right hand side matrix. The requirements of the method are moderate with respect to both computational cost and memory. Hence, it provides a possibility to tackle large scale control problems. Besides the efficient solution of the matrix equation itself, a thorough integration of the method into several control algorithms can improve their performance to a high degree. This is demonstrated for algorithms for model reduction and optimal control. Furthermore, we propose a heuristic for determining a set of suboptimal ADI shift parameters. This heuristic, which is based on a pair of Arnoldi processes, does not require any a priori knowledge on the spectrum of the coefficient matrix of the Lyapunov equation. Numerical experiments show the efficiency of the iterative scheme combined with the heuristic for the ADI parameters.

Orientador: Prof. Dr. André Fenili / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia Mecânica, Santo André, 2018.

This thesis addresses several aerospace guidance and decision making problems using both no cooperative and cooperative game theoretical solution concepts in the differential games framework. In the first part of the thesis, state dependent Riccati equation (SDRE) method has been extended to a zero-sum nonlinear differential games setting. This framework is used to study problems of intercepting a manoeuvring target, with and without terminal impact angle constraints, in the zero-sum differential games theory perspective. The guidance laws derived according to the proposed method are in closed from and online implementable. In the second part of the thesis, cooperative game theoretic concepts are applied to make a group of unmanned aerial vehicles (UAV) achieve rendezvous, in a given finite time horizon. An algorithm has been proposed that enables the UAVs to realize Nash bargaining solution. In this context, criteria for time consistency of a cooperative solution of nonzero-sum linear quadratic differential games have been studied. The problems where the UAVs try to achieve rendezvous by implementing cooperative game theoretic strategies, based on local information structure only, is also addressed.

We present a new method for the computation of low rank approximations to the solution of large, sparse, stable Lyapunov equations. It is based on a generalization of the classical Smith method and profits by the usual low rank property of the right hand side matrix. The requirements of the method are moderate with respect to both computational cost and memory. Hence, it provides a possibility to tackle large scale control problems. Besides the efficient solution of the matrix equation itself, a thorough integration of the method into several control algorithms can improve their performance to a high degree. This is demonstrated for algorithms for model reduction and optimal control. Furthermore, we propose a heuristic for determining a set of suboptimal ADI shift parameters. This heuristic, which is based on a pair of Arnoldi processes, does not require any a priori knowledge on the spectrum of the coefficient matrix of the Lyapunov equation. Numerical experiments show the efficiency of the iterative scheme combined with the heuristic for the ADI parameters.

Dans cette thèse sont montrés des développements de l'approximation de la phase aléatoire (RPA) dans le contexte de théories à séparation de portée. On présente des travaux sur le formalisme de la RPA en général, et en particulier sur le formalisme "matrice diélectrique" qui est exploré de manière systématique. On montre un résumé d'un travail sur les équations RPA dans le contexte d'orbitales localisées, notamment des développements des orbitales virtuelles localisées que sont les "orbitales oscillantes projetées" (POO). Un programme a été écrit pour calculer des fonctions telles que le trou de d'échange, la fonction de réponse, etc... sur des grilles de l'espace réel (grilles parallélépipédiques ou de type "DFT"). On montre certaines de ces visualisations. Dans l'espace réel, on expose une adaptation de l'approximation du dénominateur effectif (EED), développée originellement dans l'espace réciproque en physique du solide. Également, les gradients analytiques des énergies de corrélation RPA dans le contexte de la séparation de portée sont dérivés. Le formalisme développé ici à l'aide d'un lagrangien permet une dérivation tout-en-un des termes courte- et longue-portée qui émergent dans les expressions du gradient, et qui montrent un parallèle intéressant. Des applications sont montrées, telles que des optimisations de géométries aux niveaux RSH-dRPA-I et RSH-SOSEX d'un ensemble de 16 petites molécules, ou encore le calcul et la visualisation des densités corrélées au niveau RSH-dRPA-I / In this thesis are shown developments in the random phase approximation (RPA) in the context of range-separated theories. We present advances in the formalism of the RPA in general, and particularly in the "dielectric matrix" formulation of RPA, which is explored in details. We show a summary of a work on the RPA equations with localized orbitals, especially developments of the virtual localized orbitals that are the "projected oscillatory orbitals" (POO). A program has been written to calculate functions such as the exchange hole, the response function, etc... on real space grid (parallelepipedic or of the "DFT" type) ; some of those visualizations are shown here. In the real space, we offer an adaptation of the effective energy denominator approximation (EED), originally developed in the reciprocal space in solid physics. The analytical gradients of the RPA correlation energies in the context of range separation has been derived. The formalism developed here with a Lagrangian allows an all-in-one derivation of the short- and long-range terms that emerge in the expressions of the gradient. These terms show interesting parallels. Geometry optimizations at the RSH-dRPA-I and RSH-SOSEX levels on a set of 16 molecules are shown, as well as calculations and visualizations of correlated densities at the RSH-dRPA-I level