Für zufällige dynamische Systeme mit stetiger Zeit existieren zwei wichtige Klassen von Generatoren: Zum einen stationäre zufällige ifferentialgleichungen, i.e. gewöhnliche Differentialgleichungen, die von einem stationärer zufälligen Vektorfeld getrieben werden, und zum anderen stochastische Stratonovichdifferentialgleichungen mit weißem Rauschen. Während die erste Klasse sich gut in den ergodentheoretischen Rahmen der Theorie der zufälligen dynamischen Systeme einfügt, widersetzte sich die zweite Klasse lange Zeit der dynamischen Untersuchung aufgrund des "Konflikts zwischen Ergodentheorie und stochastischer Analysis". In dieser Arbeit wird gezeigt, daß beide Klassen von zufälligen dynamischen Systemen nicht wesentlich verschieden sind, genauer: Zu jeder stochastischen Stratonovichdifferentialgleichung mit weißem Rauschen (unter den üblichen Regularitätsforderungen an die Vektorfelder, die die Existenz von Flüssen garantieren) existiert eine stationäre zufällige Differentialgleichung derart, daß die erzeugten zufälligen dynamischen Systeme konjugiert sind. Als Anwendung wird eine Version des lokalen Linearisierungssatzes von Hartman/Grobman für stochastische Stratonovichdifferentialgleichungen bewiesen. / For continuous time random dynamical systems there exist two important classes of generators: on the one hand stationary random differential quations, i.e. ordinary differential equations driven by a stationary random vector field, and on the other hand stochastic Stratonovich differential equations with white noise. While the first class fits well into the framework of the theory of random dynamical systems, the second class resisted for a long time the dynamical investigation due to the "conflict between ergodic theory and stochastic analysis". The main result of this thesis is that both classes of random dynamical systems are not essentially distinct, more precisely: For each stochastic Stratonovich differential equation with white noise (under usual regularity assumptions) there exists a stationary random differential equation such that the corresponding random dynamical systems are conjugate. As an application a version of the local Hartman/Grobman theorem for stochastic differential equations is proved.

We apply some recent numerical solutions to diffusion bridges written in Iacus (2008). One is an approximate scheme from Bladt and S{\o}rensen (2007), another one, from  Beskos et al (2006), is an algorithm which is exact: no numerical error at given grid points!

In this thesis, we investigate the quantum mechanical interaction of light with matter in the form of a gas of ultracold atoms: the atomic ensemble. We present a theoretical analysis of two problems, which involve the interaction of quantized electromagnetic fields (called signal and idler) with the atomic ensemble (i) cascade two-photon emission in an atomic ladder configuration, and (ii) photon frequency conversion in an atomic diamond configuration. The motivation of these studies comes from potential applications in long-distance quantum communication where it is desirable to generate quantum correlations between telecommunication wavelength light fields and ground level atomic coherences. In the two systems of interest, the light field produced in the upper arm of an atomic Rb level scheme is chosen to lie in the telecom window. The other field, resonant on a ground level transition, is in the near-infrared region of the spectrum. Telecom light is useful as it minimizes losses in the optical fiber transmission links of any two long-distance quantum communication device. We develop a theory of correlated signal-idler pair correlation. The analysis is complicated by the possible generation of multiple excitations in the atomic ensemble. An analytical treatment is given in the limit of a single excitation assuming adiabatic laser excitations. The analysis predicts superradiant timescales in the idler emission in agreement with experimental observation. To relax the restriction of a single excitation, we develop a different theory of cascade emission, which is solved by numerical simulation of classical stochastic differential equation using the theory of open quantum systems. The simulations are in good qualitative agreement with the analytical theory of superradiant timescales. We further analyze the feasibility of this two-photn source to realize the DLCZ protocol of the quantum repeater communication system. We provide a quantum theory of near-infrared to telecom wavelength conversion in the diamond configuration. The system provides a crucial part of a quantum-repeater memory element, which enables a "stored" near-infrared photon to be converted to a telecom wavelength for transmission without the destruction of light-atom quantum correlation. We calculate the theoretical conversion efficiency, analyzing the role of optical depth of the ensemble, pulse length, and quantum fluctuations on the process.

<p>We apply some recent numerical solutions to diffusion bridges written in Iacus (2008). One is an approximate scheme from Bladt and S{\o}rensen (2007), another one, from  Beskos et al (2006), is an algorithm which is exact: no numerical error at given grid points!</p>

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 į tikrąją parametro reikšmę 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. Įvertiniu 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. / 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]

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.