Cette thèse porte sur le calcul d'une sur-approximation conservative pour les solutions d'équations différentielles ordinaires en présence d'incertitudes et sur son application à l'estimation et l'analyse de systèmes dynamiques à temps continu. L'avantage principal des méthodes et des algorithmes de calculs présentés dans cette thèse est qu'ils apportent une preuve numérique de résultats. Cette thèse est organisée en deux parties. La première partie est consacrée aux outils mathématiques et aux méthodes d'intégration numérique garantie des équations diff érentielles incertaines. Ces méthodes permettent de caractériser de manière garantie l'ensemble des trajectoires d'état engendrées par un système dynamique incertain dont les incertitudes sont naturellement représentées par des intervalles bornés. Dans cette optique, nous avons développé une méthode d'intégration hybride qui donne de meilleurs résultats que les méthodes d'intégration basées sur les modèles de Taylor intervalles. La seconde partie aborde les problèmes de l'identification et de l'observation dans un contexte à erreurs bornées ainsi que le problème d'atteignabilité continue pour la véri cation de propriétés des systèmes dynamiques hybrides. / This thesis addresses the computation of conservative over-approximation of the solutions of uncertain ordinary di erential equations and its application to the estimation and the analysis of uncertain continuous-time dynamical systems. The main feature of the methods and algorithms presented in this thesis is the fact that they are numerically veri ed and hence can be used to obtain numerical proof of properties. This thesis is organized in two parts. The first part is devoted to the mathematical tools and the guaranteed numerical integration methods for uncertain ordinary di erential equations. These methods make it possible to characterize in a guaranteed way all the state trajectories generated by an uncertain dynamical system whose uncertainties are in a natural way described by bounded boxes. Accordingly, we have developed a hybrid integration method which gives better results than the integration methods based on interval Taylor models. The second part is dedicated to the resolution of identi cation and observation issues in a bounded error context. It also deals with continuous reachability computation for the veri cation of the properties of hybrid dynamical systems.

Cette thèse traite de quatre sujets en rapport avec les dynamiques non-réversibles. Chacun fait l'objet d'un chapitre qui peut être lu indépendamment.Le premier chapitre est une introduction générale présentant les problématiques et quelques résultats majeurs de physique statistique computationnelle.Le second chapitre concerne la résolution numérique d'équations aux dérivées partielles hypoelliptiques, c'est-à-dire faisant intervenir un opérateur différentiel inversible mais non coercif. Nous prouvons la consistance de la méthode de Galerkin ainsi que des taux de convergence pour l'erreur. L'analyse est également conduite dans le cas d'une formulation point-selle, qui s'avère être la plus adaptée dans les cas qui nous intéressent. Nous démontrons que nos hypothèses sont satisfaites dans un cas simple et vérifions numériquement nos prédictions théoriques sur cet exemple.Dans le troisième chapitre nous proposons une stratégie générale permettant de construire des variables de contrôle pour des dynamiques hors-équilibre. Cette méthode permet en particulier de réduire la variance des estimateurs de coefficient de transport par moyenne ergodique. Cette réduction de variance est quantifiée dans un régime perturbatif. La variable de contrôle repose sur la solution d'une équation aux dérivées partielles. Dans le cas de l'équation de Langevin cette équation est hypoelliptique, ce qui motive le chapitre précédent. La méthode proposée est testée numériquement sur trois exemples.Le quatrième chapitre est connecté au troisième puisqu'il utilise la même idée de variable de contrôle. Il s'agit d'estimer la mobilité d'une particule dans le régime sous-amorti, où la dynamique est proche d'être Hamiltonienne. Ce travail a été effectué en collaboration avec G. Pavliotis durant un séjour à l'Imperial College London.Le dernier chapitre traite des processus de Markov déterministes par morceaux, qui permettent l'échantillonnage de mesure en grande dimension. Nous prouvons la convergence exponentielle vers l'équilibre de plusieurs dynamiques de ce type sous un formalisme général incluant le processus de Zig-Zag (ZZP), l'échantillonneur à particule rebondissante (BPS) et la dynamique de Monte Carlo hybride randomisée (RHMC). La dépendances des bornes sur le taux de convergence que nous démontrons sont explicites par rapport aux paramètres du problème. Cela permet en particulier de contrôler la taille des intervalles de confiance pour des moyennes empiriques lorsque la dimension de l'espace des phases sous-jacent est grande. Ce travail a été fait en collaboration avec C. Andrieu, A. Durmus et N. Nüsken. / This thesis deals with four topics related to non-reversible dynamics. Each is the subject of a chapter which can be read independently. The first chapter is a general introduction presenting the problematics and some major results of computational statistical physics. The second chapter concerns the numerical resolution of hypoelliptic partial differential equations, i.e. involving an invertible but non-coercive differential operator. We prove the consistency of the Galerkin method as well as convergence rates for the error. The analysis is also carried out in the case of a saddle-point formulation, which is the most appropriate in the cases of interest to us. We demonstrate that our assumptions are met in a simple case and numerically check our theoretical predictions on this example. In the third chapter we propose a general strategy for constructing control variates for nonequilibrium dynamics. In particular, this method reduces the variance of transport coefficient estimators by ergodic mean. This variance reduction is quantified in a perturbative regime. The control variate is based on the solution of a partial differential equation. In the case of Langevin's equation this equation is hypoelliptic, which motivates the previous chapter. The proposed method is tested numerically on three examples. The fourth chapter is connected to the third since it uses the same idea of a control variate. The aim is to estimate the mobility of a particle in the underdamped regime, where the dynamics are close to being Hamiltonian. This work was done in collaboration with G. Pavliotis during a stay at Imperial College London. The last chapter deals with Piecewise Deterministic Markov Processes, which allow measure sampling in high-dimension. We prove the exponential convergence towards the equilibrium of several dynamics of this type under a general formalism including the Zig-Zag process (ZZP), the Bouncy Particle Sampler (BPS) and the Randomized Hybrid Monte Carlo (RHMC). The dependencies of the bounds on the convergence rate that we demonstrate are explicit with respect to the parameters of the problem. This allows in particular to control the size of the confidence intervals for empirical averages when the size of the underlying phase space is large. This work was done in collaboration with C. Andrieu, A. Durmus and N. Nüsken

La première partie de cette thèse est consacrée a l'estimation par maximum de vraisemblance dans les modèles mixtes a dynamique markovienne. Nous considérons plus précisément des modèles de Markov cachés a effets mixtes et des modèles de diffusion à effets mixtes. Dans le Chapitre 2, nous combinons l'algorithme de Baum-Welch a l'algorithme SAEM pour estimer les paramètres de population dans les modèles de Markov cachés à effets mixtes. Nous proposons également des procédures spéciques pour estimer les paramètres individuels et les séquences d'états cachés. Nous étudions les propriétés de cette nouvelle méthodologie sur des données simulées et l'appliquons sur des données réelles de nombres de crises d'épilepsie. Dans le Chapitre 3, nous proposons d'abord des modèles de diffusion à effets mixtes pour la pharmacocinétique de population. Nous en estimons les paramètres en combinant l'algorithme SAEM a un filtre de Kalman étendu. Nous étudions ensuite les propriétés asymptotiques de l'estimateur du maximum de vraisemblance dans des modèles de diffusion observés sans bruit de mesure continûment sur un intervalle de temps fixé lorsque le nombre de sujets tend vers l'infini. Le Chapitre 4 est consacré à la sélection de covariables dans des modèles mixtes généraux. Nous proposons une version du BIC adaptée au contexte de double asymptotique ou le nombre de sujets et le nombre d'observations par sujet tendent vers l'infini. Nous présentons quelques simulations pour illustrer cette procédure. / The first part of this thesis deals with maximum likelihood estimation in Markovianmixed-effects models. More precisely, we consider mixed-effects hidden Markov models and mixed-effects diffusion models. In Chapter 2, we combine the Baum-Welch algorithm and the SAEM algorithm to estimate the population parameters in mixed-effects hidden Markov models. We also propose some specific procedures to estimate the individual parameters and the sequences of hidden states. We study the properties of the proposed methodologies on simulated datasets and we present an application to real daily seizure count data. In Chapter 3, we first suggest mixed-effects diffusion models for population pharmacokinetics. We estimate the parameters of these models by combining the SAEM algorithm with the extended Kalman filter. Then, we study the asymptotic properties of the maximum likelihood estimatein some mixed-effects diffusion models continuously observed on a fixed time interval when the number of subjects tends to infinity. Chapter 4 is dedicated to variable selection in general mixed-effects models. We propose a BIC adapted to the asymptotic context where both of the number of subjects and the number of observations per subject tend to infinity. We illustrate this procedure with some simulations.

L'objet de cette thèse concerne la proposition de modèles paramétriques des processus de branchement uni et multi-types. Nous mettons en valeur l’intérêt de la théorie des processus de branchement et du développement nécessaire des différents outils et de concepts propres à plusieurs domaines. Pour cela, nous commençons par rappeler quelques définitions et résultats de la théorie des processus de branchement uni et multi-types, et ce en temps discret comme en temps continu. On se consacre par la suite au développement méthodologique de ces modèles.Dans la deuxième partie de ce mémoire, nous étudions seulement l'évolution d’une seule population en temps continu, et présentons quelques familles de lois paramétriques, associées à des processus de branchement homogènes particuliers. Des méthodes récursives de calcul, ainsi que des propriétés pertinentes, concernant ces distributions de probabilité, sont dérivées des fonctions génératrices satisfaisant certaines équations aux dérivées partielles linéaires précisés. Les familles proposées seront utiles à la modélisation de systèmes plus cohérents en dynamique de populations, puisqu'on y montre que les hypothèses usuelles de distributions de Poisson ne peuvent être argumentées.Dans la troisième partie, nous étudions le comportement de l'évolution de plusieurs populations en interactions. Nous y présentons aussi des modèles paramétriques de lois, associés à des processus de branchement multi-types en temps continu et homogènes en temps. Nous considérons ensuite un modèle particulier, où une population ``mère donneuse" autonome alimente en individus K populations filles, qui sont, elles, en interaction. Ce modèle est bien adapté à l'étude des systèmes dynamiques des populations en interaction qui reste à la fois simple, mais riche en variétés de comportement. L'étude du système multi-types se fait via l'évolution des fonctions génératrices de la loi multidimensionnelles des effectifs. Pour cela, utilisant les équations différentielles ordinaires et aux dérivées partielles, nous établissons les équations implicites des distributions temporelles et multidimensionnelles, et discutons des méthodes analytiques ou numériques de leur résolution. Nous développons ensuite des exemples de modèles et en particulier celui concernant 3 et 4 populations.En conclusion, nous argumentons la pertinence de cette approche, et l’interprétation des paramètres, qui sont d'un grand intérêt pour le développement de méthodes d'inférence statistique, pour de nombreux domaines d'applications. / This thesis aims to propose parametric models for single and multi-type branching processes. The importance of the theory of branching processes is pointed out. Hence, developing various tools and specific concepts in several domains is important for applications. For those purpose, we recall some definitions and results of the single-and-multi-type branching processes theory in discrete and continuous case. Afterward, we focus on the methodological development of those models.In the second part, the evolution of a single population in the continuous case has been studied. Then, some parametric distribution families associated to particular branching mechanisms are explored. Recursive computational procedure and relevant properties concerning the associted probability distributions are derived from generating functions that satisfy specified linear partial differential equations. The suggested families are useful for the modeling of systems that are more coherent with population dynamics, contrarily to the usual hypothesis of Poisson distributions, that cannot be argued.In the third part, the evolution of different populations with interaction is explored. Similarly, some parametric models of homogeneous multi-type branching processes in continuous time are proposed. Afterwards, we consider a particular model where an autonomous donor parent population feeds in individuals, K types progeny populations that interacts. This model is well adapted to the study of dynamical systems of populations in interaction. This simple model, but has a rich variety of behaviors.The study of such systems is also done regarding the evolution of generating functions of multidimensional ndividual countrings. To achievea such study, ordinary and partial differential equations are used to establish the implicit equations of temporal and multidimensional distributions. Analytical and numerical methods for equation resolution are then discussed, and examples of particular models are developed.In conclusion, the relevancy of this approach is argumed, censidering parameters interpretation in the development of inference methods for the various applied domains.

La communauté scientifique, en relation avec les gestionnaires forestiers, travaille depuis trente ans à l'élaboration d'outils d'aide à la décision. Cependant, aucune approche de modélisation ne permet une évaluation simultanée de l'impact sur la forêt des changements globaux et de la gestion des services écosystémiques. L'élaboration d'une nouvelle génération de modèles dédiés au couplage sol-plante est indispensable pour aider les gestionnaires forestiers à adapter leurs pratiques sylvicoles face aux changements globaux. Les modèles phénoménologiques issus de la dendrométrie proposent des estimations de croissance et de production tenant compte des pratiques sylvicoles. Leur utilisation est limitée par leur dépendance à un indice de fertilité combinant de manière indifférenciée l'influence du climat et des cycles biogéochimiques. La remise en jeu de la notion d'indice de fertilité nécessite de tirer avantage des concepts issus des modèles à base écophysiologiques et biogéochimiques. La modélisation de la dynamique mensuelle du renouvellement foliaire par le modèle Stand Leaf Canopy Dynamics a permis d'estimer: i) la dynamique de l'indice foliaire, ii) la production de litière foliaire. Ce modèle probabiliste tient compte de l'influence du climat, de la disponibilité en eau de l'écosystème et repose sur des hypothèses écologiques fortes sur les mécanismes impliqués dans la longévité des feuilles. Ce modèle pose les bases de l'intégration du cycle de l'eau aux modèles dendrométriques par couplage avec un modèle écophysiologique. Le modèle de décomposition de la matière organique le long d'un profil de sol vertical proposé par Bosatta et Ågren (1996) a été analysé mathématiquement afin d'étudier le devenir des chutes de litières dans le sol forestier. Nous avons proposé: i) une preuve d'existence et l'unicité de solutions à l'équation de transport intégro-différentielle du modèle, ii) des schémas de différences finies implicites-explicites convergents estimant les solutions du modèle et iii) une discussion autour des hypothèses sous-jacentes à une troncature du modèle proposée par Bosatta et Ågren. Ces deux étapes sont des leviers indispensables à l'étude des cycles biogéochimiques dans une approche à base dendrométrique. L'influence du climat et de la phénologie sur le cycle de l'eau, la décomposition des litières foliaires et la chimie du sol constituera une prochaine étape de modélisation. / The scientific community, in collaboration with forestry managers, developed decision support tools since thirty years. However, current modelling approaches do not allow simultaneous estimations of global changes impact management policies on forests. The design of a new generation of models dedicated to soil-plant coupling is necessary to help forest managers to adapt forestry practices to face global changes. Phenomenological models arising from dendrometry estimate forest growth \& yield and take into account of silvicultural practices. Their use is limited by their dependence to a fertility index, which combines implicitly the influence of climate and biogeochemical cycles. A re-examination of fertility index concept is necessary to take advantage of concepts developed in ecophysiological and biogeochemical models. Modelling the monthly dynamics of foliar renewal with the Stand Leaf Canopy Dynamics model permited: i) to estimate leaf area index dynamics, ii) to model foliar litter production. This probabilistic model takes into account climate influence, ecosystem water availability and lies on strong ecological hypothesis on mechanisms involved in leaf longevity. This model, associated with an ecophysiological model, permits to integrate water cycle to dendrometric models. The soil organic matter decomposition along a vertical soil profile model introduced by Bosatta et Ågren (1996) was analyzed to study the future of litterfall in forest soil. We proposed: i) a proof of existence and uniqueness of solutions to the integro-differential transport equation of the model, ii) convergent implicit-explicit finite differences schemes estimating model solutions and iii) a discussion on hypothesis underlying a model truncation introduced by Bosatta et Ågren. These two modelling steps are essential to consider biogeochemical cycle with a dendrometric approach. Climate and phenology influence on water cycle, litterfall decomposition and soil chemistry are the elements of a future modelling phase.

La simulation numérique de systèmes d’équations différentielles raides ordinaires ou algébriques est devenue partie intégrante dans le processus de conception des systèmes mécaniques à dynamiques complexes. L’objet de ce travail est de développer des méthodes numériques pour réduire les temps de calcul par le parallélisme en suivant deux axes : interne à l’intégrateur numérique, et au niveau de la décomposition de l’intervalle de temps. Nous montrons l’efficacité limitée au nombre d’étapes de la parallélisation à travers les méthodes de Runge-Kutta et DIMSIM. Nous développons alors une méthodologie pour appliquer le complément de Schur sur le système linéarisé intervenant dans les intégrateurs par l’introduction d’un masque de dépendance construit automatiquement lors de la mise en équations du modèle. Finalement, nous étendons le complément de Schur aux méthodes de type "Krylov Matrix Free". La décomposition en temps est d’abord vue par la résolution globale des pas de temps dont nous traitons la parallélisation du solveur non-linéaire (point fixe, Newton-Krylov et accélération de Steffensen). Nous introduisons les méthodes de tirs à deux niveaux, comme Parareal et Pita dont nous redéfinissons les finesses de grilles pour résoudre les problèmes raides pour lesquels leur efficacité parallèle est limitée. Les estimateurs de l’erreur globale, nous permettent de construire une extension parallèle de l’extrapolation de Richardson pour remplacer le premier niveau de calcul. Et nous proposons une parallélisation de la méthode de correction du résidu. / This PhD Thesis deals with the development of parallel numerical methods for solving Ordinary and Algebraic Differential Equations. ODE and DAE are commonly arising when modeling complex dynamical phenomena. We first show that the parallelization across the method is limited by the number of stages of the RK method or DIMSIM. We introduce the Schur complement into the linearised linear system of time integrators. An automatic framework is given to build a mask defining the relationships between the variables. Then the Schur complement is coupled with Jacobian Free Newton-Krylov methods. As time decomposition, global time steps resolutions can be solved by parallel nonlinear solvers (such as fixed point, Newton and Steffensen acceleration). Two steps time decomposition (Parareal, Pita,...) are developed with a new definition of their grids to solved stiff problems. Global error estimates, especially the Richardson extrapolation, are used to compute a good approximation for the second grid. Finally we propose a parallel deferred correction

La performance énergétique est un pilier pour réduire l'utilisation d'énergie non renouvelable, en plus de l'utilisation des énergies renouvelables. En fait, les bâtiments sont au cœur de la politique des performances énergétiques de l'UE puisque 40% de la consommation finale d'énergie et 36% des émissions de gaz à effet de serre provient des bureaux, magasins et autres bâtiments. Les bâtiments peuvent être considérés comme des systèmes dynamiques et le transfert de la chaleur dans les bâtiments peut être représenté en utilisant des modèles dynamiques. De cette façon, le transfert de la chaleur dans les bâtiments peut être décrit par des réseaux thermiques obtenus en utilisant la théorie des graphes et de la thermodynamique, et peuvent être déduits de l'équation de la chaleur classique. Les réseaux thermiques peuvent être exprimés comme un système d'équations différentielles et algébriques (DAE) qui peut être transformé en représentation d'état et obtenir un fonction de transfert à partir de laquelle un modèle autorégressif avec des variables exogènes (ARX) peut être obtenu. Ces différentes structures de modèle peuvent être utilisées pour identifier les paramètres physiques des réseaux thermiques, ce qui implique que la méthode peut être utilisée pour identifier la performance intrinsèque des bâtiments et aider à la réduction de la consommation d'énergie dans les bâtiments.Cela peut faciliter l'évaluation de la performance énergétique des bâtiments dans un cadre reproductible qui permet la comparaison entre différentes solutions constructives.Les principales contributions originales de cette thèse sont: 1) les réseaux thermiques sont présentées à partir de la théorie des graphes et de la thermodynamique, sans considérer l'analogie thermique-électrique; 2) l'équation classique de la chaleur est reliée explicitement avec un système de DAE (réseau thermique) par les éléments finis; 3) différentes transformations pour déduire des modèles du transfert de la chaleur avec signification physique, à partir de l'équation de la chaleur classique, sont présentées toutes ensemble; 4) les transformations entre les modèles sont effectuées à partir des réseaux thermiques jusqu’aux modèles autorégressifs avec des variables exogènes (ARX) et vice-versa; et 5) un critère de sélection de l'ordre du modèle par une analyse de fréquence des mesures est proposé. / Energy efficiency is one of the two pillars to decrease the use of non-renewable energy besides the use of renewables energies. In fact, buildings are central to the EU's energy efficiency policy, as nearly 40% of the final energy consumption and 36% of greenhouse gas emissions take place in houses, offices, shops and other buildings. Buildings may be considered as dynamic systems and heat transfer in buildings may be represented using dynamic models. In this way, heat transfer in buildings may be described by thermal networks which may be stated considering graph theory and thermodynamics, and may be deduced from the classical heat equation. Thermal networks may be expressed as a system of linear differential algebraic equations (DAE) and the system of linear DAE may be transformed into a state-space representation from which an autoregressive model with exogenous (ARX) can be obtained. These different model structures may be used for identifying the physical parameters of thermal networks which implies that this methodology may be useful for identifying the intrinsic performance of buildings and tackling the reduction of non-renewable energy consumption in buildings. This may facilitate the assessment of energy efficiency of buildings within a reproducible framework which allows the comparison between different constructive solutions.The main original contributions of this dissertation are: 1) thermal networks are stated from graph theory and thermodynamics, leaving back the thermal-electrical analogy; 2) classical heat equation is connected explicitly to a system of DAE (thermal network) by using the finite elements; 3) the transformations for deducing heat transfer models with physical meaning from the classical heat equation are put altogether; 4) transformations between models may are done from thermal networks to autoregressive models with exogenous (ARX) and back; and 5) a criterion for selecting the order of the model by frequency analysis of measurements is proposed.

No description available.

Cette thèse porte sur l'étude d'équations différentielles de type frottement, c'est à dire d'équations de type attractive, avec un unique point stable 0, caractérisant la vitesse d'un objet soumis à une force de frottement. La vitesse de cet objet subit des perturbations aléatoires de type Lévy. Dans une première partie, nous nous intéressons aux propriétés fondamentales de ces EDS : existence et unicité de la solution, caractère markovien et ergodique de celle-ci et plus particulièrement le cas des processus de Lévy stable. Dans une deuxième partie, nous étudions la stabilité de la solution de ces EDS lorsque la perturbation est un processus de Lévy stable qui tend vers 0. En effet, nous démontrons l'existence d'un développement limité d'ordre un autour de la solution déterministe pour la vitesse et la position de l'objet. Dans une troisième partie, nous étudions le comportement asymptotique des solutions lorsque la vitesse initiale est nulle et que la perturbation est un processus de Lévy stable symétrique. Nous prouvons dans cette partie que l'accumulation de perturbations entraîne un comportement asymptotique gaussien de la position de l'objet, à condition que l'indice de stabilité du processus de Lévy et la croissance du potentiel soient suffisamment grand. Dans une quatrième partie, nous levons l'hypothèse de symétrie de la perturbation en démontrant le même résultat que dans la troisième partie mais avec une dérive. Pour cela, nous étudions tout d'abord la queue de distribution de la mesure invariante associée à la vitesse de l'objet. Enfin dans une dernière partie, nous nous intéressons au résultat de la troisième partie lorsque la perturbation est la somme d'un mouvement brownien et d'un processus de Lévy purement à sauts. Puis nous commençons l'étude de la dimension deux en traitant le cas où les équations sont découplées mais où les mouvement brownien directeurs sont dépendants. / This thesis deals with the study of friction type differential equations, in other words, attractive equations, with a unique stable point 0, describing the speed of an object submitted to a frictional force. This object's speed is disturbed by Lévy type random perturbations. In a first part, one is interested in fondamental properties of these SDE: existence and unicity of a solution, Markov and ergodic properties, and more particularly the case of stable Lévy processes.In a second part, one study the stability of the solution of these SDE when the perturbation is an stable Lévy process that tends to 0. In fact, one proves the existence of a Taylor expansion of order one around the deterministic solution for the object's speed and position. In a third part, one study the asymptotic behaviour of the solutions when the initial speed is 0 and the perturbation is a symmetric stable Lévy process. One proves that the amount of perturbations, if the stability's index of the Lévy process and the increasing of the potential are big enough, leads to a gaussian asymptotic behaviour for the object's position.In a forth part, one relaxes the assumption of symmetry of the perturbation by proving the same result as in the third part but with a drift. To do so, one first studies the tail of the invariant measure of the object's speed.Finally, in a last part, one is interested in the same result as in the third part when the perturbation is the sum of the Brownian motion and a pure jump stable Lévy process. Then, one begins the study of the dimension two by considering the case where the equations are separated but where the driving Brownian motions are dependent.

Cette thèse aborde deux sujets de recherches, le premier est sur l’existence et l’unicité des solutions des Équations Différentielles Doublement Stochastiques Rétrogrades (EDDSRs) et les Équations aux Dérivées partielles Stochastiques (EDPSs) multidimensionnelles à croissance surlinéaire. Le deuxième établit l’existence d’un contrôle optimal strict pour un système controlé dirigé par des équations différentielles stochastiques progressives rétrogrades (EDSPRs) couplées dans deux cas de diffusions dégénérée et non dégénérée.• Existence et unicité des solutions des EDDSRs multidimensionnels :Nous considérons EDDSR avec un générateur de croissance surlinéaire et une donnée terminale de carré intégrable. Nous introduisons une nouvelle condition locale sur le générateur et nous montrons qu’elle assure l’existence, l’unicité et la stabilité des solutions. Même si notre intérêt porte sur le cas multidimensionnel, notre résultat est également nouveau en dimension un. Comme application, nous établissons l’existence et l’unicité des solutions des EDPS semi-linéaires.• Contrôle des EDSPR couplées :Nous étudions un problème de contrôle avec une fonctionnelle coût non linéaire dont le système contrôlé est dirigé par une EDSPR couplée. L’objective de ce travail est d’établir l’existence d’un contrôle optimal dans la classe des contrôle stricts, donc on montre que ce contrôle vérifie notre équation et qu’il minimise la fonctionnelle coût. La méthode consiste à approcher notre système par une suite de systèmes réguliers et on montre la convergence. En passant à la limite, sous des hypothèses de convexité, on obtient l’existence d’un contrôle optimal strict. on suit cette méthode théorique pour deux cas différents de diffusions dégénérée et non dégénérée. / In this Phd thesis, we considers two parts. The first one establish the existence and the uniquness of the solutions of multidimensional backward doubly stochastic differential equations (BDSDEs in short) and the stochastic partial differential equations (SPDEs in short) in the superlinear growth generators. In the second part, we study the stochastic controls problems driven by a coupled Forward-Backward stochastic differentialequations (FBSDEs in short).• BDSDEs and SPDEs with a superlinear growth generators :We deal with multidimensional BDSDE with a superlinear growth generator and a square integrable terminal datum. We introduce new local conditions on the generator then we show that they ensure the existence and uniqueness as well as the stability of solutions. Our work go beyond the previous results on the subject. Although we are focused on multidimensional case, the uniqueness result we establish is new in one dimensional too. As application, we establish the existence and uniqueness of probabilistic solutions tosome semilinear SPDEs with superlinear growth generator. By probabilistic solution, we mean a solution which is representable throughout a BDSDEs.• Controlled coupled FBSDEs :We establish the existence of an optimal control for a system driven by a coupled FBDSE. The cost functional is defined as the initial value of the backward component of the solution. We construct a sequence of approximating controlled systems, for which we show the existence of a sequence of feedback optimal controls. By passing to the limit, we get the existence of a feedback optimal control. The convexity condition is used to ensure that the optimal control is strict. In this part, we study two cases of diffusions : degenerate and non-degenerate.