L’être humain est organisé selon une horloge interne d’une période d’environ 24 heures. La pharmacocinétique de certaines classes de médicaments est donc influencée par le rythme circadien. En effet, l’aire sous la courbe de la concentration en médicament en fonction du temps, la concentration maximale en médicament et le temps auquel on obtient la concentration maximale peuvent varier en fonction de l’heure à laquelle a été consommé le médicament. Le but de ce travail est de modéliser la variation de la concentration maximale de ces médicaments selon le moment de la journée auquel ils sont pris. On étudie d’abord un modèle présenté par Godfrey permettant de trouver la concen- tration en médicament en fonction du temps et tenant compte des variations circadiennes. Ce modèle ne permet pas d’illustrer les variations dans la concentration maximale selon le moment de la journée auquel le médicament est pris. Un nouveau modèle à deux com- partiments sera donc développé pour les trois modes d’absorption (orale, intraveineuse, intraveineuse bolus). Les systèmes d’équations différentielles résultants seront étudiés. L’effet de la variation des paramètres de phase sur la concentration maximale sera aussi étudié. La preuve de l’existence des solutions, de leur unicité et de leur positivité sera faite en annexe. / Humans are organised according to an internal clock with a period of approximatively 24 hours. The pharmacokinetic of several classes of drugs are then influenced by circadian rhythms. Indeed, the area under the curve (of the drug concentration as a function of time), the maximal concentration and the time to maximal concentration can change according to the time at which the drug is taken. The objective of this present work is to find a model to represent the variations in the maximal drug concentration according to the absorption’s time. We first study a model presented by Godfrey. It allows to find the drug concentration as a function of time while taking into account circadian rhythms. Unfortunately, this model could not represent the variations in the maximal concentration according to the time at which the drug is taken. We developed a new two-compartmental model for the three ways of absorption (oral, intravenous and intravenous bolus). The resulting systems of ordinary differential equations will be studied. The effect of the phase parameters on the maximal concen- tration will also be studied. Finally, the proof of well-poseness of the model will be developed in the Annex.

La première partie de cette thèse est consacrée aux méthodes numériques pour la simulation de processus aléatoires définis par des équations différentielles stochastiques (EDS). Nous commençons par l’étude de l’algorithme de Beskos et al. [13] qui permet de simuler exactement les trajectoires d’un processus solution d’une EDS en dimension 1. Nous en proposons une extension à des fins de calcul exact d’espérances et nous étudions l’application de ces idées à l’évaluation du prix d’options asiatiques dans le modèle de Black & Scholes. Nous nous intéressons ensuite aux schémas numériques. Dans le deuxième chapitre, nous proposons deux schémas de discrétisation pour une famille de modèles à volatilité stochastique et nous en étudions les propriétés de convergence. Le premier schéma est adapté à l’évaluation du prix d’options path-dependent et le deuxième aux options vanilles. Nous étudions également le cas particulier où le processus qui dirige la volatilité est un processus d’Ornstein-Uhlenbeck et nous exhibons un schéma de discrétisation qui possède de meilleures propriétés de convergence. Enfin, dans le troisième chapitre, il est question de la convergence faible trajectorielle du schéma d’Euler. Nous apportons un début de réponse en contrôlant la distance de Wasserstein entre les marginales du processus solution et du schéma d’Euler, uniformément en temps. La deuxième partie de la thèse porte sur la modélisation de la dépendance en finance et ce à travers deux problématiques distinctes : la modélisation jointe entre un indice boursier et les actions qui le composent et la gestion du risque de défaut dans les portefeuilles de crédit. Dans le quatrième chapitre, nous proposons un cadre de modélisation original dans lequel les volatilités de l’indice et de ses composantes sont reliées. Nous obtenons un modèle simplifié quand la taille de l’indice est grande, dans lequel l’indice suit un modèle à volatilité locale et les actions individuelles suivent un modèle à volatilité stochastique composé d’une partie intrinsèque et d’une partie commune dirigée par l’indice. Nous étudions la calibration de ces modèles et montrons qu’il est possible de se caler sur les prix d’options observés sur le marché, à la fois pour l’indice et pour les actions, ce qui constitue un avantage considérable. Enfin, dans le dernier chapitre de la thèse, nous développons un modèle à intensités permettant de modéliser simultanément, et de manière consistante, toutes les transitions de ratings qui surviennent dans un grand portefeuille de crédit. Afin de générer des niveaux de dépendance plus élevés, nous introduisons le modèle dynamic frailty dans lequel une variable dynamique inobservable agit de manière multiplicative sur les intensités de transitions. Notre approche est purement historique et nous étudions l’estimation par maximum de vraisemblance des paramètres de nos modèles sur la base de données de transitions de ratings passées / The first part of this thesis deals with probabilistic numerical methods for simulating the solution of a stochastic differential equation (SDE). We start with the algorithm of Beskos et al. [13] which allows exact simulation of the solution of a one dimensional SDE. We present an extension for the exact computation of expectations and we study the application of these techniques for the pricing of Asian options in the Black & Scholes model. Then, in the second chapter, we propose and study the convergence of two discretization schemes for a family of stochastic volatility models. The first one is well adapted for the pricing of vanilla options and the second one is efficient for the pricing of path-dependent options. We also study the particular case of an Orstein-Uhlenbeck process driving the volatility and we exhibit a third discretization scheme which has better convergence properties. Finally, in the third chapter, we tackle the trajectorial weak convergence of the Euler scheme by providing a simple proof for the estimation of the Wasserstein distance between the solution and its Euler scheme, uniformly in time. The second part of the thesis is dedicated to the modelling of dependence in finance through two examples : the joint modelling of an index together with its composing stocks and intensity-based credit portfolio models. In the forth chapter, we propose a new modelling framework in which the volatility of an index and the volatilities of its composing stocks are connected. When the number of stocks is large, we obtain a simplified model consisting of a local volatility model for the index and a stochastic volatility model for the stocks composed of an intrinsic part and a systemic part driven by the index. We study the calibration of these models and show that it is possible to fit the market prices of both the index and the stocks. Finally, in the last chapter of the thesis, we define an intensity-based credit portfolio model. In order to obtain stronger dependence levels between rating transitions, we extend it by introducing an unobservable random process (frailty) which acts multiplicatively on the intensities of the firms of the portfolio. Our approach is fully historical and we estimate the parameters of our model to past rating transitions using maximum likelihood techniques

Cette thèse comporte plusieurs procédures d'estimation non-paramétrique de densité de probabilité.Dans chaque cas les variables d'intérêt ne sont pas observées directement, ce qui est une difficulté majeure.La première partie traite un modèle linéaire mixte où des observations répétées sont disponibles.La deuxième partie s'intéresse aux modèles d'équations différentielles stochastiques à effets aléatoires. Plusieurs trajectoires sont observées en temps continu sur un intervalle de temps commun.La troisième partie se place dans un contexte de bruit multiplicatif.Les différentes parties de cette thèse sont reliées par un contexte commun de problème inverse et par une problématique commune: l'estimation de la densité d'une variable cachée. Dans les deux premières parties la densité d'un ou plusieurs effets aléatoires est estimée. Dans la troisième partie il s'agit de reconstruire la densité de la variable d'origine à partir d'observations bruitées.Différentes méthodes d'estimation globale sont utilisées pour construire des estimateurs performants: estimateurs à noyau, estimateurs par projection ou estimateurs construits par déconvolution.La sélection de paramètres mène à des estimateurs adaptatifs et les risques quadratiques intégrés sont majorés grâce à une inégalité de concentration de Talagrand. Une étude sur simulations de chaque estimateur illustre leurs performances. Un jeu de données neuronales est étudié grâce aux procédures mises en place pour les équations différentielles stochastiques. / This thesis contains several nonparametric estimation procedures of a probability density function.In each case, the main difficulty lies in the fact that the variables of interest are not directly observed.The first part deals with a mixed linear model for which repeated observations are available.The second part focuses on stochastic differential equations with random effects. Many trajectories are observed continuously on the same time interval.The third part is in a full multiplicative noise framework.The parts of the thesis are connected by the same context of inverse problems and by a common problematic: the estimation of the density function of a hidden variable.In the first two parts the density of one or two random effects is estimated. In the third part the goal is to rebuild the density of the original variable from the noisy observations.Different global methods are used and lead to well competitive estimators: kernel estimators, projection estimators or estimators built from deconvolution.Parameter selection gives adaptive estimators and the integrated risks are bounded using a Talagrand concentration inequality.A simulation study for each proposed estimator highlights their performances.A neuronal dataset is investigated with the new procedures for stochastic differential equations developed in this work.

Contribution à l’étude du phénomène d’oscillation argumentaire. L’oscillateur argumentaire a un mouvement stable périodique, à une fréquence proche de sa fréquence fondamentale, lorsqu’il est soumis à une excitation provenant d’une source de type harmonique, à une fréquence qui est un multiple de ladite fréquence fondamentale, et agissant de manière telle que son interaction avec le système dépende des coordonnées d’espace du système. La présente thèse étudie quelques systèmes argumentaires et essaie de mettre en évidence des relations symboliques entre les paramètres de ces systèmes et leur comportement observé ou calculé. C’est la représentation de Van der Pol qui a été utilisée la plupart du temps pour représenter l’état du système, car elle est bien adaptée à la méthode de centrage, où l'on cherche une solution sous forme d’un signal de type sinusoïdal, d’amplitude et de phase lentement variables. L’originalité de la présente thèse vis-à-vis des publications antérieures est dans la modélisation, plus proche des systèmes physiques réels, dans les développements symboliques qui donnent des représentations inédites, dans le mode de réalisation des expériences, qui utilisent toutes une visualisation de Van der Pol en temps réel, et dans l’objet de l’expérience de la poutre excitée axialement de manière argumentaire. Au cours de cette thèse, des systèmes simples à un DDL ont été modélisés, construits et expérimentés. Des relations symboliques, notamment concernant les probabilités de capture par des attracteurs, ainsi que des critères de stabilité et une solution symbolique approchée, ont été mis en évidence. Un système continu constitué d’une poutre élancée excitée axialement a ensuite été modélisé à l’aide de deux modèles et expérimenté ; toujours dans le domaine symbolique, des propriétés ont été étudiées, notamment concernant des combinaisons de plages de paramètres permettant au phénomène argumentaire d’exister / Contribution to the study of the argumental oscillation phenomenon. The argumental oscillator has a stable periodic motion at a frequency close to its fundamental frequency when it is subjected to an excitation from a harmonic source at a frequency which is a multiple of said fundamental frequency, and acting in such a way that its interaction with the system depends on the space coordinates of the system. This thesis studies some argumental systems and tries to demonstrate symbolic relations between the parameters of these systems and their observed or calculated behavior. The Van der Pol representation was used most of the time to represent the state of the system, as it is well adapted to the averaging method, where a solution is sought as a signal of sinusoidal type, with slowly varying amplitude and phase. The originality of this thesis with respect to previous publications is in the modeling, closer to real physical systems, in the symbolic developments that give new representations, in the embodiment of the experiments, which all use a real-time Van der Pol visualization, and in the object of the experiment of the beam axially excited in an argumental way. During this thesis, simple systems with one DDL have been modeled, built and tested. Symbolic relationships have been highlighted, in particular with regard to the probabilities of capture by attractors, as well as stability criteria and an approximate symbolic solution. A continuous system consisting of an axially excited slender beam was then modeled using two models, and tested; still in the symbolic domain, properties have been studied, especially concerning combinations of parameter ranges allowing the argumental phenomenon to occur

Les problèmes de gestion des réservoirs sont stochastiques principalement à cause de l’incertitude sur les apports naturels. Ceci entraine des modèles d’optimisation de grande taille pouvant être difficilement traitables numériquement. La première partie de cette thèse réexamine la méthode d’agrégation de scénarios proposée par Rockafellar et Wets (1991). L’objectif consiste à améliorer la vitesse de convergence de l’algorithme du progressive hedgging sur lequel repose la méthode. L’approche traditionnelle consiste à utiliser une valeur fixe pour ce paramètre ou à l’ajuster selon une trajectoire choisie a priori : croissante ou décroissante. Une approche dynamique est proposée pour mettre à jour le paramètre en fonction d’information sur la convergence globale fournie par les solutions à chaque itération. Il s’agit donc d’une approche a posteriori. La thèse aborde aussi la gestion des réservoirs par la programmation dynamique stochastique. Celle-ci se prête bien à ces problèmes de gestion à cause de la nature séquentielle de leurs décisions opérationnelles. Cependant, les applications sont limitées à un nombre restreint de réservoirs. La complexité du problème peut augmenter exponentiellement avec le nombre de variables d’état, particulièrement quand l’approche classique est utilisée, i.e. en discrétisant l’espace des états de « manière uniforme ». La thèse propose une approche d’approximation sur une grille irrégulière basée sur une décomposition simpliciale de l’espace des états. La fonction de valeur est évaluée aux sommets de ces simplexes et interpolée ailleurs. À l’aide de bornes sur la vraie fonction, la grille est raffinée tout en contrôlant l’erreur d’approximation commise. En outre, dans un contexte décision-information spécifique, une hypothèse « uni-bassin », souvent utilisée par les hydrologues, est exploitée pour développer des formes analytiques pour l’espérance de la fonction de valeur. Bien que la méthode proposée ne résolve pas le problème de complexité non polynomiale de la programmation dynamique, les résultats d’une étude de cas industrielle montrent qu’il n’est pas forcément nécessaire d’utiliser une grille très dense pour approximer la fonction de valeur avec une précision acceptable. Une bonne approximation pourrait être obtenue en évaluant cette fonction uniquement en quelques points de grille choisis adéquatement. / Reservoir operation problems are in essence stochastic because of the uncertain nature of natural inflows. This leads to very large optimization models that may be difficult to handle numerically. The first part of this thesis revisits the scenario aggregation method proposed by Rochafellar and Wets (1991). Our objective is to improve the convergence of the progressive hedging algorithm on which the method is based. This algorithm is based on an augmented Lagrangian with a penalty parameter that plays an important role in its convergence. The classical approach consists in using a fixed value for the parameter or in adjusting it according a trajectory chosen a priori: decreasing or increasing. This thesis presents a dynamic approach to update the parameter based on information on the global convergence provided by the solutions at each iteration. Therefore, it is an a posteriori scheme. The thesis also addresses reservoir problems via stochastic dynamic programming. This scheme is widely used for such problems because of the sequential nature of the operational decisions of reservoir management. However, dynamic programing is limited to a small number of reservoirs. The complexity may increase exponentially with the dimension of the state variables, especially when the classical approach is used, i.e. by discretizing the state space into a "regular grid". This thesis proposes an approximation scheme over an irregular grid based on simplicial decomposition of the state space. The value function is evaluated over the vertices of these simplices and interpolated elsewhere. Using bounds on the true function, the grid is refined while controlling the approximation error. Furthermore, in a specific information-decision context, a "uni-bassin" assumption often used by hydrologists is exploited to develop analytical forms for the expectation of the value function. Though the proposed method does not eliminate the non-polynomial complexity of dynamic programming, the results of an industrial case study show that it is not absolutely necessary to use a very dense grid to appropriately approximate the value function. Good approximation may be obtained by evaluating this function at few appropriately selected grid points.

Cette thèse traite d'une part, de questions de gestion, de mesure et de transfert du risque et d'autre part, de problèmes d'analyse stochastique à sauts avec incertitude de modèle. Le premier chapitre est consacré à l'analyse des intégrales de Choquet, comme mesures de risque monétaires non nécessairement invariantes en loi. Nous établissons d'abord un nouveau résultat de représentation des mesures de risque comonotones, puis un résultat de représentation des intégrales de Choquet en introduisant la notion de distorsion locale. Ceci nous permet de donner ensuite une forme explicite à l'inf-convolution de deux intégrales de Choquet, avec des exemples illustrant l'impact de l'absence de la propriété d'invariance en loi. Nous nous intéressons ensuite à un problème de tarification d'un contrat de réassurance non proportionnelle, contenant des clauses de reconstitution. Après avoir défini le prix d'indifférence relatif à la fois à une fonction d'utilité et à une mesure de risque, nous l'encadrons par des valeurs facilement implémentables. Nous passons alors à un cadre dynamique en temps. Pour cela, nous montrons, en adoptant une approche par point fixe, un théorème d'existence de solutions bornées pour une classe d'équations différentielles stochastiques rétrogrades (EDSRs dans la suite) avec sauts et à croissance quadratique. Sous une hypothèse additionnelle classique dans le cadre à sauts, ou sous une hypothèse de convexité du générateur, nous établissons un résultat d'unicité grâce à un principe de comparaison. Nous analysons les propriétés des espérances non linéaires correspondantes. En particulier, nous obtenons une décomposition de Doob-Meyer des surmartingales non-linéaires ainsi que leur régularité en temps. En conséquence, nous en déduisons facilement un principe de comparaison inverse. Nous appliquons ces résultats à l'étude des mesures de risque dynamiques associées, sur une filtration engendrée à la fois par un mouvement brownien et par une mesure aléatoire à valeurs entières, à leur repésentation duale, ainsi qu'à leur inf-convolution, avec des exemples explicites. La seconde partie de cette thèse concerne l'analyse de l'incertitude de modèle, dans le cas particulier des EDSRs du second ordre avec sauts. Nous imposons que ces équations aient lieu au sens presque-sûr, pour toute une famille non dominée de mesures de probabilités qui sont solution d'un problème de martingales sur l'espace de Skorohod. Nous étendons d'abord la définition des EDSRs du second ordre, telles que définies par Soner, Touzi et Zhang, au cas avec sauts. Pour ce faire, nous démontrons un résultat d'agrégation au sens de Soner, Touzi et Zhang sur l'espace des trajectoires càdlàg. Ceci nous permet, entre autres, d'utiliser une version quasi-sûre du compensateur de la mesure des sauts du processus canonique. Nous montrons alors un résultat d'existence et d'unicité pour notre classe d'EDSRs du second ordre. Ces équations sont affectées par l'incertitude portant à la fois sur la volatilité et sur les sauts du processus qui les dirige.

Cette thèse porte sur la théorie de Galois différentielle. Elle est divisée en deux parties. La première concerne la théorie de Galois différentielle paramétrée, et la seconde, les équations aux q-différences. Dans le chapitre 2, nous exposons une généralisation de l'algorithme de Kovacic qui permet de calculer le groupe de Galois paramétré de certaines équations différentielles paramétrées d'ordre 2. Dans le chapitre 3, nous présentons une généralisation du théorème de densité de Ramis qui donne un ensemble de générateurs topologiques du groupe de Galois pour les équations différentielles linéaires paramétrées à coefficients dans un anneau convenable. Nous obtenons une contribution au problème inverse dans cette théorie de Galois, donnons un critère d'isomonodromie, et répondons partiellement à une question posée par Sibuya. Dans le chapitre 4, il est question de confluence et d'équations aux q-différences. Nous prouvons comment la transformée de Borel-Laplace d'une série formelle divergente solution d'une équation différentielle linéaire à coefficients dans C(z) peut être uniformément approchée par un q-analogue de la transformée de Borel-Laplace appliqué à une série formelle solution d'une famille d'équations aux q-différences linéaires qui discrétise l'équation différentielle. Nous faisons directement les calculs dans le cas des séries hypergéométriques basiques, et nous prouvons sous des hypothèses raisonnables, qu'une matrice fondamentale d'une équation différentielle linéaire à coefficients dans C(z) peut être uniformément approchée par une matrice fondamentale d'une famille d'équations aux q-différences linéaires correspondante.

Le sujet de la thèse est l'étude et la construction de méthodes numériques géométriques pour les équations différentielles, qui préservent des propriétés géométriques du flot exact, notamment la symétrie, la symplecticité des systèmes hamiltoniens, la conservation d'intégrales premières, la structure de Poisson, etc.<br />Dans la première partie, on introduit une nouvelle approche de construction d'intégrateurs numériques géométriques d'ordre élevé en s'inspirant de la théorie des équations différentielles modifiées. Le cas des méthodes développables en B-séries est spécifiquement analysé et on introduit une nouvelle loi de composition sur les B-séries. L'efficacité de cette approche est illustrée par la construction d'un nouvel intégrateur géométrique d'ordre élevé pour les équations du mouvement d'un corps rigide. On obtient également une méthode numérique précise pour le calcul de points conjugués pour les géodésiques du corps rigide.<br />Dans la seconde partie, on étudie dans quelle mesure les excellentes performances des méthodes symplectiques, pour l'intégration à long terme en astronomie et en dynamique moléculaire, persistent pour les problèmes de contrôle optimal. On discute également l'extension de la théorie des équations modifiées aux problèmes de contrôle optimal.<br />Dans le même esprit que les équations modifiées, on considère dans la dernière partie des méthodes de pas fractionnaire (splitting) pour les systèmes hamiltoniens perturbés, utilisant des potentiels modifiés. On termine par la construction de méthodes de splitting d'ordre élevé avec temps complexes pour les équations aux dérivées partielles paraboliques, notamment les problèmes de réaction-diffusion en chimie.

Cette thèse est divisée en trois parties. Dans la première partie, nous commençons par décrire les théories de Ziglin, Yoshida et Morales-Ramis et les motiver. Dans la deuxième partie, on étudie l’intégrabilité des équations différentielles de Newton à trois degrés de liberté dont les forces sont des polynômes homogènes de degrés trois. En utilisant une analyse du groupe de Galois différentiel des équations aux variations d’ordre supérieur, nous faisons une classification (presque) complète des forces génériques et intégrables. Dans une dernière partie, nous intéressons à l’intégrabilité d’un système d’équations différentielles homogènes d’ordre un (système A). L’application directe de la théorie de Morales-Ramis ne donne des obstructions à l’intégrabilité. En dérivant le système A par rapport au temps, nous obtenons un système différentiel de Newton homogène d’ordre 2 (système B). L’avantage est que ce dernier possède des solutions particulières algébriquement non triviales et le critère classique de Morales-Ramis nous permet d’établir des conditions nécessaires d’intégrabilité. Nous prouvons qu’il existe des relations explicites entre les intégrales premières des deux systèmes et nous introduisons une nouvelle méthode de recherche d’intégrales premières que l’on appelle « Extension tangente double ». Nous appliquons cette méthode à des systèmes planaires homogènes quadratiques. Comme deuxième application, nous montrons que, sous certaines conditions, les racines newtoniennes d’un système différentiel de Newton avec force centrale sont intégrables par quadratures. Nous présentons plusieurs systèmes intégrables avec deux, trois et quatre degrés de liberté. / This thesis is divided into three parts. In the first part we begin by describing the theories of Ziglin, Yoshida and Morales-Ramis and motivating them. In the second part we study the integrability of three-dimensional differential Newton equations with homogeneous polynomial forces of degree three. Using an analysis of differential Galois group of higher order variational equations, we give an almost complete classification of integrable generic forces. The last part is devoted to a study of the integrability of a system of first order homogeneous differential equations (system A ). The direct application of the Morales-Ramis theory does not lead to obstructions to the integrability. If we differentiate the differential system A with respect to time, we obtain a homogeneous Newtonian system (system B). The advantage is that the system B has a non-trivial particular solution and the classical criterion of Morales-Ramis allows us to establish necessary conditions for integrability. We prove that there are explicit relationships between first integrals of the both systems and we introduce a new method for finding first integrals called ``Double tangent extension method''. We apply the obtained results for a detailed analysis of homogeneous planar differential system. Using the double tangent extension method, we formulate some conditions under which the Newtonian roots of Newton's system with central force are integrable by quadratures. Some new cases of integrability with two, three and four degrees of freedom are found.

Depuis la publication de l'ouvrage de Choquet (1955), la théorie d'espérance non linéaire a attiré avec grand intérêt des chercheurs pour ses applications potentielles dans les problèmes d'incertitude, les mesures de risque et le super-hedging en finance. Shige Peng a construit une sorte d'espérance entièrement non linéaire dynamiquement cohérente par l'approche des EDP. Un cas important d'espérance non linéaire cohérente en temps est la G-espérance, dans laquelle le processus canonique correspondant (B_{t})_{t≥0} est appelé G-mouvement brownien et joue un rôle analogue au processus de Wiener classique. L'objectif de cette thèse est d'étudier, dans le cadre de la G-espérance, certaines équations différentielles stochastiques rétrogrades (G-EDSR) à croissance quadratique avec applications aux problèmes de maximisation d'utilité robuste avec incertitude sur les modèles, certaines équations différentielles stochastiques (G-EDS) réfléchies et équations différentielles stochastiques rétrogrades réfléchies avec générateurs lipschitziens. On considère d'abord des G-EDSRs à croissance quadratique. Dans le Chapitre 2 nous fournissons un resultat d'existence et unicité pour des G-EDSRs à croissance quadratique. D'une part, nous établissons des estimations a priori en appliquant le théorème de type Girsanov, d'où l'on en déduit l'unicité. D'autre part, pour prouver l'existence de solutions, nous avons d'abord construit des solutions pour des G-EDSRs discretes en résolvant des EDPs non-linéaires correspondantes, puis des solutions pour les G-EDSRs quadratiques générales dans les espaces de Banach. Dans le Chapitre 3 nous appliquons les G-EDSRs quadratiques aux problèmes de maximisation d'utilité robuste. Nous donnons une caratérisation de la fonction valeur et une stratégie optimale pour les fonctions d'utilité exponentielle, puissance et logarithmique. Dans le Chapitre 4, nous traitons des G-EDSs réfléchies multidimensionnelles. Nous examinons d'abord la méthode de pénalisation pour résoudre des problèmes de Skorokhod déterministes dans des domaines non convexes et établissons des estimations pour des fonctions α-Hölder continues. A l'aide de ces résultats obtenus pour des problèmes déterministes, nous définissons le G-mouvement Brownien réfléchi et prouvons son existence et son unicité dans un espace de Banach. Ensuite, nous prouvons l'existence et l'unicité de solution pour les G-EDSRs multidimensionnelles réfléchies via un argument de point fixe. Dans le Chapitre 5, nous étudions l'existence et l'unicité pour les équations différentielles stochastiques rétrogrades réfléchies dirigées par un G-mouvement brownien lorsque la barrière S est un processus de G-Itô. / Since the publication of Choquet's (1955) book, the theory of nonlinear expectation has attracted great interest from researchers for its potential applications in uncertainty problems, risk measures and super-hedging in finance. Shige Peng has constructed a kind of fully nonlinear expectation dynamically coherent by the PDE approach. An important case of time-consistent nonlinear expectation is G-expectation, in which the corresponding canonical process (B_{t})_{t≥0} is called G-Brownian motion and plays a similar role to the classical Wiener process. The objective of this thesis is to study, in the framework of the G-expectation, some backward stochastic differential equations (G-BSDE) under a quadratic growth condition on their coefficients with applications to robust utility maximization problems with uncertainty on models, Reflected stochastic differential equations (reflected G-SDE) and reflected backward stochastic differential equations with Lipschitz coefficients (reflected G-BSDE). We first consider G-BSDE with quadratic growth. In Chapter 2 we provide a result of existence and uniqueness for quadratic G-BSDEs. On the one hand, we establish a priori estimates by applying the Girsanov-type theorem, from which we deduce the uniqueness. On the other hand, to prove the existence of solutions, we first constructed solutions for discrete G-BSDEs by solving corresponding nonlinear PDEs, then solutions for the general quadratic G-BSDEs in the spaces of Banach. In Chapter 3 we apply quadratic G-BSDE to robust utility maximization problems. We give a characterization of the value function and an optimal strategy for exponential, power and logarithmic utility functions. In Chapter 4, we discuss multidimensional reflected G-SDE. We first examine the penalization method to solve deterministic Skorokhod problems in non-convex domains and establish estimates for continuous α-Hölder functions. Using these results for deterministic problems, we define the reflected G-Brownian motion and prove its existence and its uniqueness in a Banach space. Then we prove the existence and uniqueness of the solution for the multidimensional reflected G-SDE via a fixed point argument. In Chapter 5, we study the existence and uniqueness of the reflected backward stochastic differential equations driven by a G-Brownian motion when the obstacle S is a G-Itô process.