10 results
Search Results
Now showing 1 - 10 of 10
- ItemSolving linear parabolic rough partial differential equations(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Bayer, Christian; Belomestny, Denis; Redmann, Martin; Riedel, Sebastian; Schoenmakers, JohnWe study linear rough partial differential equations in the setting of [Friz and Hairer, Springer, 2014, Chapter 12]. More precisely, we consider a linear parabolic partial differential equation driven by a deterministic rough path W of Hölder regularity with 1=3 < 1=2. Based on a stochastic representation of the solution of the rough partial differential equation, we propose a regression Monte Carlo algorithm for spatio-temporal approximation of the solution. We provide a full convergence analysis of the proposed approximation method which essentially relies on the new bounds for the higher order derivatives of the solution in space. Finally, a comprehensive simulation study showing the applicability of the proposed algorithm is presented.
- ItemAn H2-type error bound for time-limited balanced truncation(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Redmann, Martin; Kürschner, PatrickWhen solving partial differential equations numerically, usually a high order spatial discretization is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatially-discretized systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is balanced truncation (BT). However, if one aims at finding a good ROM on a certain finite time interval only, time-limited BT (TLBT) can be a more accurate alternative. So far, no error bound on TLBT has been proved. In this paper, we close this gap in the theory by providing an H2 error bound for TLBT with two different representations. The performance of the error bound is then shown in several numerical experiments.
- ItemBalanced truncation and singular perturbation approximation model order reduction for stochastically controlled linear systems(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Redmann, Martin; Freitag, Melina A.When solving linear stochastic differential equations numerically, usually a high order spatial discretisation is used. Balanced truncation (BT) and singular perturbation approximation (SPA) are well-known projection techniques in the deterministic framework which reduce the order of a control system and hence reduce computational complexity. This work considers both methods when the control is replaced by a noise term. We provide theoretical tools such as stochastic concepts for reachability and observability, which are necessary for balancing related model order reduction of linear stochastic differential equations with additive Lévy noise. Moreover, we derive error bounds for both BT and SPA and provide numerical results for a specific example which support the theory.
- ItemDynamic programming for optimal stopping via pseudo-regression(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Bayer, Christian; Redmann, Martin; Schoenmakers, John G.M.We introduce new variants of classical regression-based algorithms for optimal stopping problems based on computation of regression coefficients by Monte Carlo approximation of the corresponding L2 inner products instead of the least-squares error functional. Coupled with new proposals for simulation of the underlying samples, we call the approach pseudo regression. We show that the approach leads to asymptotically smaller errors, as well as less computational cost. The analysis is justified by numerical examples.
- ItemType II singular perturbation approximation for linear systems with Lévy noise(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Redmann, MartinWhen solving linear stochastic partial differential equations numerically, usually a high order spatial discretisation is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatially-discretised systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is singular perturbation approximation (SPA), a method which has been extensively studied for deterministic systems. As so-called type I SPA it has already been extended to stochastic equations. We provide an alternative generalisation of the deterministic setting to linear systems with Lévy noise which is called type II SPA. It turns out that the ROM from applying type II SPA has better properties than the one of using type I SPA. In this paper, we provide new energy interpretations for stochastic reachability Gramians, show the preservation of mean square stability in the ROM by type II SPA and prove two different error bounds for type II SPA when applied to Lévy driven systems.
- ItemTowards time-limited H2-optimal model order reduction(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Goyal, Pawan; Redmann, MartinIn order to solve partial differential equations numerically and accurately, a high order spatial discretization is usually needed. Model order reduction (MOR) techniques are often used to reduce the order of spatially-discretized systems and hence reduce computational complexity. A particular class of MOR techniques are H2-optimal methods such as the iterative rational Krylov subspace algorithm (IRKA) and related schemes. However, these methods are used to obtain good approximations on a infinite time-horizon. Thus, in this work, our main goal is to discuss MOR schemes for time-limited linear systems. For this, we propose an alternative time-limited H2-norm and show its connection with the time-limited Gramians. We then provide first-order optimality conditions for an optimal reduced order model (ROM) with respect to the time-limited H2-norm. Based on these optimality conditions, we propose an iterative scheme which upon convergences aims at satisfying these conditions. Then, we analyze how far away the obtained ROM is from satisfying the optimality conditions.We test the efficiency of the proposed iterative scheme using various numerical examples and illustrate that the newly proposed iterative method can lead to a better reduced-order compared to unrestricted IRKA in the time interval of interest.
- ItemLow-dimensional approximations of high-dimensional asset price models(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Redmann, Martin; Bayer, Christian; Goyal, PawanWe consider high-dimensional asset price models that are reduced in their dimension in order to reduce the complexity of the problem or the effect of the curse of dimensionality in the context of option pricing. We apply model order reduction (MOR) to obtain a reduced system. MOR has been previously studied for asymptotically stable controlled stochastic systems with zero initial conditions. However, stochastic differential equations modeling price processes are uncontrolled, have non-zero initial states and are often unstable. Therefore, we extend MOR schemes and combine ideas of techniques known for deterministic systems. This leads to a method providing a good pathwise approximation. After explaining the reduction procedure, the error of the approximation is analyzed and the performance of the algorithm is shown conducting several numerical experiments. Within the numerics section, the benefit of the algorithm in the context of option pricing is pointed out.
- ItemEnergy estimates and model order reduction for stochastic bilinear systems(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Redmann, MartinIn this paper, we investigate a large-scale stochastic system with bilinear drift and linear diffusion term. Such high dimensional systems appear for example when discretizing a stochastic partial differential equations in space. We study a particular model order reduction technique called balanced truncation (BT) to reduce the order of spatially-discretized systems and hence reduce computational complexity. We introduce suitable Gramians to the system and prove energy estimates that can be used to identify states which contribute only very little to the system dynamics. When BT is applied the reduced system is obtained by removing these states from the original system. The main contribution of this paper is an L2-error bound for BT for stochastic bilinear systems. This result is new even for deterministic bilinear equations. In order to achieve it, we develop a new technique which is not available in the literature so far.
- ItemRunge--Kutta methods for rough differential equations(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Redmann, Martin; Riedel, SebastianWe study Runge-Kutta methods for rough differential equations which can be used to calculate solutions to stochastic differential equations driven by processes that are rougher than a Brownian motion. We use a Taylor series representation (B-series) for both the numerical scheme and the solution of the rough differential equation in order to determine conditions that guarantee the desired order of the local error for the underlying Runge-Kutta method. Subsequently, we prove the order of the global error given the local rate. In addition, we simplify the numerical approximation by introducing a Runge-Kutta scheme that is based on the increments of the driver of the rough differential equation. This simplified method can be easily implemented and is computational cheap since it is derivative-free. We provide a full characterization of this implementable Runge-Kutta method meaning that we provide necessary and sufficient algebraic conditions for an optimal order of convergence in case that the driver, e.g., is a fractional Brownian motion with Hurst index 1/4 < H ≤ 1/2. We conclude this paper by conducting numerical experiments verifying the theoretical rate of convergence.
- ItemType II balanced truncation for deterministic bilinear control systems(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Redmann, MartinWhen solving partial differential equations numerically, usually a high order spatial discretisation is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatially-discretised systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is balanced truncation (BT), a method which has been extensively studied for deterministic linear systems. As so-called type I BT it has already been extended to bilinear equations, an important subclass of nonlinear systems. We provide an alternative generalisation of the linear setting to bilinear systems which is called type II BT. The Gramians that we propose in this context contain information about the control. It turns out that the new approach delivers energy bounds which are not just valid in a small neighbourhood of zero. Furthermore, we provide an H1-error bound which so far is not known when applying type I BT to bilinear systems.