The main topic of this dissertation is the estimation of the Hurst index H of the solutions of stochastic differential equations (SDEs) driven by the fractional Brownian motion (fBm). Firstly, the limit behavior of the first and second order quadratic variations of the solutions of SDEs driven by the fBm is analyzed. This yields several strongly consistent estimators of the Hurst index H. Secondly, it is proved that in case the solution of the SDE is replaced by its Milstein approximation, the estimators remain strongly consistent. Additionally, the possibilities of applying the increment ratios (IR) statistic based estimator of H originally obtained by J. M. Bardet and D. Surgailis in 2010 to the fractional geometric Brownian motion are examined. Furthermore, this dissertation derives the convergence rate of the modified Gladyshev’s estimator of the Hurst index to its real value. The estimators obtained in the dissertation were compared with several other known estimators of the Hurst index H, namely the naive and ordinary least squares Gladyshev and eta-summing oscillation estimators, the variogram estimator and the IR estimator. The models chosen for comparison of these estimators were the fractional Ornstein-Uhlenbeck (O-U) process and the fractional geometric Brownian motion (gBm). The initial inference about the behavior of these estimators was drawn for the O-U process which is Gaussian, while the gBm process was used to check how the estimators behave in a... [to full text] / Pagrindinė šios disertacijos tema – stochastinių diferencialinių lygčių (SDL), valdomų trupmeninio Brauno judesio (tBj), sprendinių Hursto indekso H vertinimas. Pirmiausia disertacijoje išnagrinėta SDL, valdomų tBj, sprendinių pirmos ir antros eilės kvadratinių variacijų ribinė elgsena. Iš šių rezultatų seka keli stipriai pagrįsti Hursto indekso H įvertiniai. Įrodyta, kad šie įvertiniai išlieka stipriai pagrįsti, jei tikra sprendinio trajektorija keičiama jos Milšteino aproksimacija. Taip pat išnagrinėtos pokyčių santykio (increment ratios) statistikos H įvertinio, gauto J. M. Bardeto ir D. Surgailio 2010 m., taikymo trupmeninio geometrinio Brauno judesio Hursto indekso vertinimui galimybės bei nustatytas modifikuoto Gladyševo H įvertinio konvergavimo i tikrąją parametro reikšme greitis. Gauti įvertiniai palyginti su kai kuriais kitais žinomais Hursto indekso H įvertiniais: naiviais bei mažiausių kvadratų Gladyševo ir eta-sumavimo osciliacijos įvertiniais, variogramos įvertiniu ir pokyčių santykio statistikos įvertiniu. Įvertinių elgsena buvo palyginta trupmeniniam Ornšteino-Ulenbeko (OU) procesui bei trupmeniniam geometriniam Brauno judesiui (gBj). Pradinės išvados buvo padarytos O-U procesui, kuris yra Gauso, o gBj procesas buvo naudojamas patikrinti, kaip šie įvertiniai elgiasi, kai procesas yra ne Gauso. Disertaciją sudaro įvadas, 3 pagrindiniai skyriai, išvados, literatūros sąrašas, autoriaus publikacijų disertacijos tema sąrašas ir du priedai.

Cette thèse traite des équations différentielles stochastiques rétrogrades réfléchies du second ordre dans une filtration générale . Nous avons traité tout d'abord la réflexion à une barrière inférieure puis nous avons étendu le résultat dans le cas d'une barrière supérieure. Notre contribution consiste à démontrer l'existence et l'unicité de la solution de ces équations dans le cadre d'une filtration générale sous des hypothèses faibles. Nous remplaçons la régularité uniforme par la régularité de type Borel. Le principe de programmation dynamique pour le problème de contrôle stochastique robuste est donc démontré sous les hypothèses faibles c'est à dire sans régularité sur le générateur, la condition terminal et la barrière. Dans le cadre des Équations Différentielles Stochastiques Rétrogrades (EDSRs ) standard, les problèmes de réflexions à barrières inférieures et supérieures sont symétriques. Par contre dans le cadre des EDSRs de second ordre, cette symétrie n'est plus valable à cause des la non linéarité de l'espérance sous laquelle est définie notre problème de contrôle stochastique robuste non dominé. Ensuite nous un schéma d'approximation numérique d'une classe d'EDSR de second ordre réfléchies. En particulier nous montrons la convergence de schéma et nous testons numériquement les résultats obtenus. / This thesis deals with the second-order reflected backward stochastic differential equations (2RBSDEs) in general filtration. In the first part , we consider the reflection with a lower obstacle and then extended the result in the case of an upper obstacle . Our main contribution consists in demonstrating the existence and the uniqueness of the solution of these equations defined in the general filtration under weak assumptions. We replace the uniform regularity by the Borel regularity(through analytic measurability). The dynamic programming principle for the robust stochastic control problem is thus demonstrated under weak assumptions, that is to say without regularity on the generator, the terminal condition and the obstacle. In the standard Backward Stochastic Differential Equations (BSDEs) framework, there is a symmetry between lower and upper obstacles reflection problem. On the contrary, in the context of second order BSDEs, this symmetry is no longer satisfy because of the nonlinearity of the expectation under which our robust stochastic non-dominated stochastic control problem is defined. In the second part , we get a numerical approximation scheme of a class of second-order reflected BSDEs. In particular we show the convergence of our scheme and we test numerically the results.

>Magister Scientiae - MSc / We present various methods of pricing Asian options. The methods include Monte Carlo simulations designed using control and antithetic variates, numerical solution of partial differential equation and using lower bounds.The price of the Asian option is known to be a certain risk-neutral expectation. Using the Feynman-Kac theorem, we deduce that the problem of determining the expectation implies solving a linear parabolic partial differential equation. This partial differential equation does not admit explicit solutions due to the fact that the distribution of a sum of lognormal variables is not explicit. We then solve the partial differential equation numerically using finite difference and Monte Carlo methods.Our Monte Carlo approach is based on the pseudo random numbers and not deterministic sequence of numbers on which Quasi-Monte Carlo methods are designed. To make the Monte Carlo method more effective, two variance reduction techniques are discussed.Under the finite difference method, we consider explicit and the Crank-Nicholson’s schemes. We demonstrate that the explicit method gives rise to extraneous solutions because the stability conditions are difficult to satisfy. On the other hand, the Crank-Nicholson method is unconditionally stable and provides correct solutions. Finally, we apply the pricing methods to a similar problem of determining the price of a European-style arithmetic basket option under the Black-Scholes framework. We find the optimal lower bound, calculate it numerically and compare this with those obtained by the Monte Carlo and Moment Matching methods.Our presentation here includes some of the most recent advances on Asian options, and we contribute in particular by adding detail to the proofs and explanations. We also contribute some novel numerical methods. Most significantly, we include an original contribution on the use of very sharp lower bounds towards pricing European basket options.

On étudie une équation différentielle stochastique animée par un processus ponctuel de Poisson, qui modélise un changement continu de lénvironnement d'une population et la fixation stochastique de mutations bénéfiques pour compenser ce changement. La probabilité de fixation d'une mutation augmente dès que le retard phénotypique $X_t$ entre la population et l'optimum augmente. On suppose que les mutations favorables se fixent instantanément induisant un saut adaptatif. En premier lieu, on a étudié le comportement à long terme de la solution de cette équation sachant qu'on ne considère qu'un seul trait phénotypique de la population et on a trouvé les conditions sous lesquelles $X_t$ est récurrent (possibilité de survie) ou transient (extinction inévitable). Ensuite, on a généralisé nos résultats en considérant un vecteur de traits phénotypiques de la population, essentiellement dans $mathbb R^2$. A la fin, on introduit une limite des petits sauts pour caractériser et comprendre le cas récurrent. / We study a stochastic differential equation driven by a Poisson point process, which models continuous changes in a population's environment, as well as the stochastic fixation of beneficial mutations that might compensate for this change. The fixation probability of a given mutation increases as the phenotypic lag $X_t$ between the population and the optimum grows larger, and successful mutations are assumed to fix instantaneously (leading to an adaptive jump). First, we study the large time behavior of the solution of this SDE taking into consideration one phenotypic trait of the population and we find the conditions under which $X_t$ is recurrent (possibility of survival) or transient (doomed to exctinction).Then we generalize our results to the case of a phenotypic traits vector, essentially in $R^2$. Finally, we introduce a small jumps limit to characterize and understand the recurrent case.

Reatores tubulares de polimerização podem apresentar um perfil de velocidade bastante distorcido. Partindo desta observação, um modelo estocástico baseado no modelo de dispersão axial foi proposto para a representação matemática da fluidodinâmica de um reator tubular para produção de poliestireno. A equação diferencial foi obtida inserindo a aleatoriedade no parâmetro de dispersão, resultando na adição de um termo estocástico ao modelo capaz de simular as oscilações observadas experimentalmente. A equação diferencial estocástica foi discretizada e resolvida pelo método Euler-Maruyama de forma satisfatória. Uma função estimadora foi desenvolvida para a obtenção do parâmetro do termo estocástico e o parâmetro do termo determinístico foi calculado pelo método dos mínimos quadrados. Uma análise de convergência foi conduzida para determinar o número de elementos da discretização e o modelo foi validado através da comparação de trajetórias e de intervalos de confiança computacionais com dados experimentais. O resultado obtido foi satisfatório, o que auxilia na compreensão do comportamento fluidodinâmico complexo do reator estudado. / The velocity profile of polymerization tubular reactors may be very distorted. Based on this observation, a stochastic model based on the axial dispersion model was proposed for the mathematical representation of the fluid dynamics of a tubular reactor for polystyrene production. The differential equation was built by inserting randomness in the dipersion coefficient, which added a stochastic term to the model. This term was capable of simulating the experimentally observed fluctuations. The stochastic differential equation was discretized and solved by the Euler-Maruyama method adequately. An estimator function has been developed to calculate the parameter of the stochastic term, while the parameter of the deterministic term was estimated by a least squares method. A convergence analysis was carried out in order to determine the number of elements needed for the time discretization. The model was validated through comparisons of sample paths and computational confidence intervals with experimental data. The result was considered satisfactory, allowing a better understanding of the complex fluid dynamic behaviour of the analised reactor. Key-words: modelling, simulation, stochastic differential equation, polymerization tubular reactor, time residence distribution.

Cette thèse traite les jeux différentiels stochastiques de somme non nulle (JDSNN) dans le cadre de Markovien et de leurs liens avec les équations différentielles stochastiques rétrogrades (EDSR) multidimensionnelles. Nous étudions trois problèmes différents. Tout d'abord, nous considérons un JDSNN où le coefficient de dérive n'est pas borné, mais supposé uniquement à croissance linéaire. Ensuite certains cas particuliers de coefficients de diffusion non bornés sont aussi considérés. Nous montrons que le jeu admet un point d'équilibre de Nash via la preuve de l'existence de la solution de l'EDSR associée et lorsque la condition d'Isaacs généralisée est satisfaite. La nouveauté est que le générateur de l'EDSR, qui est multidimensionnelle, est de croissance linéaire stochastique par rapport au processus de volatilité. Le deuxième problème est aussi relatif au JDSNN mais les payoffs ont des fonctions d'utilité exponentielles. Les EDSRs associées à ce jeu sont de type multidimensionnelles et quadratiques en la volatilité. Nous montrons de nouveau l'existence d’un équilibre de Nash. Le dernier problème que nous traitons, est un jeu bang-bang qui conduit à des hamiltoniens discontinus. Dans ce cas, nous reformulons le théorème de vérification et nous montrons l’existence d’un équilibre de Nash qui est du type bang-bang, i.e., prenant ses valeurs sur le bord du domaine en fonction du signe de la dérivée de la fonction valeur ou du processus de volatilité. L'EDSR dans ce cas est un système multidimensionnel couplé, dont le générateur est discontinu par rapport au processus de volatilité. / This dissertation studies the multiple players nonzero-sum stochastic differential games (NZSDG) in the Markovian framework and their connections with multiple dimensional backward stochastic differential equations (BSDEs). There are three problems that we are focused on. Firstly, we consider a NZSDG where the drift coefficient is not bound but is of linear growth. Some particular cases of unbounded diffusion coefficient of the diffusion process are also considered. The existence of Nash equilibrium point is proved under the generalized Isaacs condition via the existence of the solution of the associated BSDE. The novelty is that the generator of the BSDE is multiple dimensional, continuous and of stochastic linear growth with respect to the volatility process. The second problem is of risk-sensitive type, i.e. the payoffs integrate utility exponential functions, and the drift of the diffusion is unbounded. The associated BSDE is of multi-dimension whose generator is quadratic on the volatility. Once again we show the existence of Nash equilibria via the solution of the BSDE. The last problem that we treat is a bang-bang game which leads to discontinuous Hamiltonians. We reformulate the verification theorem and we show the existence of a Nash point for the game which is of bang-bang type, i.e., it takes its values in the border of the domain according to the sign of the derivatives of the value function. The BSDE in this case is a coupled multi-dimensional system, whose generator is discontinuous on the volatility process.

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.

No description available.

This thesis investigated Stochastic Runge–Kutta Lawson (SRKL) schemes and their application to the Heston model. Two distinct SRKL discretization methods were used to simulate a single asset’s dynamics under the Heston model, notably the Euler–Maruyama and Midpoint schemes. Additionally, standard Monte Carlo and variance reduction techniques were implemented. European and Asian option prices were estimated and compared with a benchmark value regarding accuracy, effectiveness, and computational complexity. Findings showed that the SRKL Euler–Maruyama schemes exhibited promise in enhancing the price for simple and path-dependent options. Consequently, integrating SRKL numerical methods into option valuation provides notable advantages by addressing challenges posed by the Heston model’s SDEs. Given the limited scope of this research topic, it is imperative to conduct further studies to understand the use of SRKL schemes within other models.

Ce travail, à la frontière de l’analyse et des probabilités, s’intéresse à l’étude de systèmes différentiels a priori mal posés. Nous cherchons, grâce à des techniques issues de la théorie des chemins rugueux et de l’étude trajectorielle des processus stochastiques, à donner un sens à de tels systèmes puis à les résoudre, tout en montrant que les notions proposées ici étendent bien les notions classiques de solutions. Cette thèse se décompose en trois chapitres. Le premier traite des systèmes différentiels ordinaires perturbés additivement par des processus irréguliers éventuellement stochastiques ainsi que des effets de régularisation de tels processus. Le deuxième chapitre concerne l’équation de transport linéaire perturbée multiplicativement par des chemins rugueux ; enfin, le dernier chapitre s’intéresse à une équation de la chaleur non linéaire perturbée par un bruit blanc espace-temps, l’équation de quantisation stochastique phi4 en dimension 3. / In this work we investigate a priori ill-posed differential systems from an analytic and probabilistic point of view. Thanks to technics inspired by the rough path theory and pathwise study of stochastic processes, we want to define those ill-posed systems and then study them. The first chapter of this thesis is related to ordinary differential equations perturbed by some irregular (stochastic) processes and the effects induced by the regularization of such processes. The second chapter deals with the linear transport equation multiplicatively perturbed by a rough path. Finally, in the last chapter we investigate the stochastic quantization equation Phi4 in three dimensions.