M.C. Veraar
Please Note
23 records found
1
In many cases, population data for all species is unavailable. To make accurate estimates for the unknown or hidden data, the extended Kalman filter can be applied. Which, through a combination of the data and the mathematical model, creates an estimate for the population. An exponential bound for the error of this estimation is derived in expectation ...
In many cases, population data for all species is unavailable. To make accurate estimates for the unknown or hidden data, the extended Kalman filter can be applied. Which, through a combination of the data and the mathematical model, creates an estimate for the population. An exponential bound for the error of this estimation is derived in expectation
The main challenge in the numerical approximation of such equations is the computation of conditional expectations over potentially high-dimensional spaces. In classical settings, where the dimensionality of the underlying randomness is moderate, many approaches have been proposed in the literature. However, for high-dimensional problems, one has to resort toMonte Carlo methods. In recent years, a new class of regression Monte Carlo methods has arisen in the literature, so called deep BSDE methods, which practically approximate the solution of BSDEs in a neural network regression Monte Carlo framework, after forming a suitable loss function motivated either by stochastic optimal control or the martingale representation theorem. These classes of methods can roughly be divided into two main categories. Forward methods, where the solution of the associated backward SDE is simultaneously optimized in a global optimization, minimizing a loss function stemming from a stochastic target problem reformulation. Alternatively, backward methods have been investigated, where the numerical resolution of the equation is decomposed into smaller sub-optimizations corresponding to a discrete set of points in a suitable time discretization. These methods enabled the numerical treatment of longstanding open challenges, such as the pricing and deltahedging of multi-asset financial options up to d = 100 risk factors and beyond.
The goal of this thesis is to analyze such modern machine learning based numerical methods, and apply them in the financial mathematical context. We propose numerical extensions of these methods in high-dimensional frameworks, analyze their convergence properties in discrete time, and investigate their robustness and accuracy in practical applications such as hedging and stochastic optimal control. Ourmain contributions in each chapter can be summarized as follows…
...
The main challenge in the numerical approximation of such equations is the computation of conditional expectations over potentially high-dimensional spaces. In classical settings, where the dimensionality of the underlying randomness is moderate, many approaches have been proposed in the literature. However, for high-dimensional problems, one has to resort toMonte Carlo methods. In recent years, a new class of regression Monte Carlo methods has arisen in the literature, so called deep BSDE methods, which practically approximate the solution of BSDEs in a neural network regression Monte Carlo framework, after forming a suitable loss function motivated either by stochastic optimal control or the martingale representation theorem. These classes of methods can roughly be divided into two main categories. Forward methods, where the solution of the associated backward SDE is simultaneously optimized in a global optimization, minimizing a loss function stemming from a stochastic target problem reformulation. Alternatively, backward methods have been investigated, where the numerical resolution of the equation is decomposed into smaller sub-optimizations corresponding to a discrete set of points in a suitable time discretization. These methods enabled the numerical treatment of longstanding open challenges, such as the pricing and deltahedging of multi-asset financial options up to d = 100 risk factors and beyond.
The goal of this thesis is to analyze such modern machine learning based numerical methods, and apply them in the financial mathematical context. We propose numerical extensions of these methods in high-dimensional frameworks, analyze their convergence properties in discrete time, and investigate their robustness and accuracy in practical applications such as hedging and stochastic optimal control. Ourmain contributions in each chapter can be summarized as follows…
The preliminary concepts necessary to interpret such equations are collected in the first chapter.
The second chapter is devoted to analyzing a class of equations which generalizes the spatial Whittle–Matérn stochastic partial differential equations (SPDEs) to space–time. We define solution concepts, establish their well-posedness and equivalence, and study their spatiotemporal regularity and covariance structure. The abstract spatial operator in these equations is assumed to be the negative generator of a strongly continuous semigroup on a separable Hilbert space; for certain spatial regularity results, it is moreover assumed to be analytic.
In the third chapter, we define various higher-order Markov properties for Hilbert-space-valued stochastic processes and investigate the relations between them. For solutions to linear abstract SPDEs, we identify two sets of additional conditions under which locality of the precision operator is either necessary or sufficient to satisfy the weakest Markov property. We show that the mild solutions from Chapter 2 satisfy a higher-order Markov property if the orders of the equations are non-fractional (i.e., integer) and that, conversely, a necessary condition for the weakest Markov property is, in general, not satisfied if this parameter is fractional. We moreover establish that an infinite-dimensional analog of the fractional Brownian motion can be obtained as a limiting case of this class of equations.
The fourth chapter concerns the deterministic natural Dirichlet problem for nonlocal abstract space–time operators posed on an arbitrary Banach space. We impose that the whole "past" of the solution equals a given function, i.e., its values are prescribed at all times preceding some (non-trivial) initial time. If the semigroup associated to the spatial operator is exponentially stable, we show that the problem is well-posed in an L^p-sense with p in [1, infinity]. Whenever such solutions are continuous, they satisfy a mild solution formula which expresses them in terms of the initial data and the semigroup, thus generalizing the well-known variation of constants formula for the first-order abstract Cauchy problem. A comparison to analogous solution concepts arising from Riemann–Liouville and Caputo type initial value problems is included.
Finally, in the fifth chapter, we study the convergence of a sequence of semilinear parabolic stochastic evolution equations posed on a sequence of Banach spaces approximating a limiting space. The abstract "discrete-to-continuum approximation" setting is encoded using projection and lifting operators. These allow us to define and compare the discretized equations, and to formulate conditions under which their solutions are well posed and convergent when lifted to a common state space. Our framework is applied to the case where the limiting problem is an SPDE whose linear part is a generalized Whittle–Matérn operator on a manifold, discretized by a sequence of geometric graphs constructed from a (random) point cloud. ...
The preliminary concepts necessary to interpret such equations are collected in the first chapter.
The second chapter is devoted to analyzing a class of equations which generalizes the spatial Whittle–Matérn stochastic partial differential equations (SPDEs) to space–time. We define solution concepts, establish their well-posedness and equivalence, and study their spatiotemporal regularity and covariance structure. The abstract spatial operator in these equations is assumed to be the negative generator of a strongly continuous semigroup on a separable Hilbert space; for certain spatial regularity results, it is moreover assumed to be analytic.
In the third chapter, we define various higher-order Markov properties for Hilbert-space-valued stochastic processes and investigate the relations between them. For solutions to linear abstract SPDEs, we identify two sets of additional conditions under which locality of the precision operator is either necessary or sufficient to satisfy the weakest Markov property. We show that the mild solutions from Chapter 2 satisfy a higher-order Markov property if the orders of the equations are non-fractional (i.e., integer) and that, conversely, a necessary condition for the weakest Markov property is, in general, not satisfied if this parameter is fractional. We moreover establish that an infinite-dimensional analog of the fractional Brownian motion can be obtained as a limiting case of this class of equations.
The fourth chapter concerns the deterministic natural Dirichlet problem for nonlocal abstract space–time operators posed on an arbitrary Banach space. We impose that the whole "past" of the solution equals a given function, i.e., its values are prescribed at all times preceding some (non-trivial) initial time. If the semigroup associated to the spatial operator is exponentially stable, we show that the problem is well-posed in an L^p-sense with p in [1, infinity]. Whenever such solutions are continuous, they satisfy a mild solution formula which expresses them in terms of the initial data and the semigroup, thus generalizing the well-known variation of constants formula for the first-order abstract Cauchy problem. A comparison to analogous solution concepts arising from Riemann–Liouville and Caputo type initial value problems is included.
Finally, in the fifth chapter, we study the convergence of a sequence of semilinear parabolic stochastic evolution equations posed on a sequence of Banach spaces approximating a limiting space. The abstract "discrete-to-continuum approximation" setting is encoded using projection and lifting operators. These allow us to define and compare the discretized equations, and to formulate conditions under which their solutions are well posed and convergent when lifted to a common state space. Our framework is applied to the case where the limiting problem is an SPDE whose linear part is a generalized Whittle–Matérn operator on a manifold, discretized by a sequence of geometric graphs constructed from a (random) point cloud.
The local well-posedness of quasi-linear problems in a critical setting by the maximal $L^p$-regularity theory from Chapter 18 of Hytönen, van Neerven, Veraar, and Weis (2024) are investigated, where we use the non-autonomous setting with non-constant domains from Di Giorgio, Lunardi, and Schnaubelt (2005). Dominant examples in the literature of such problems are problems with multidimensional, non-constant Neumann boundary conditions influencing the domain of the operator. In this thesis, we look for ways to ensure the short-timescale existence and uniqueness of solutions to these problems and research the possibility of applying them to the model problem with Neumann boundary conditions. By applying non-autonomous linear theory from Chapter 3 Part II of Yagi (2010), we find a result that allows us to determine the existence and uniqueness of short timescale mild solutions. When applied, however, we see that, unlike the work of Yagi (2010), we can only guarantee the local well-posedness of the model Neumann problem by using averaging functions because of our more strict critical setting. In the future, results can be based on different regularity types, or the given result can be applied on spaces with negative smoothness. ...
The local well-posedness of quasi-linear problems in a critical setting by the maximal $L^p$-regularity theory from Chapter 18 of Hytönen, van Neerven, Veraar, and Weis (2024) are investigated, where we use the non-autonomous setting with non-constant domains from Di Giorgio, Lunardi, and Schnaubelt (2005). Dominant examples in the literature of such problems are problems with multidimensional, non-constant Neumann boundary conditions influencing the domain of the operator. In this thesis, we look for ways to ensure the short-timescale existence and uniqueness of solutions to these problems and research the possibility of applying them to the model problem with Neumann boundary conditions. By applying non-autonomous linear theory from Chapter 3 Part II of Yagi (2010), we find a result that allows us to determine the existence and uniqueness of short timescale mild solutions. When applied, however, we see that, unlike the work of Yagi (2010), we can only guarantee the local well-posedness of the model Neumann problem by using averaging functions because of our more strict critical setting. In the future, results can be based on different regularity types, or the given result can be applied on spaces with negative smoothness.
Using stochastic compactness arguments the existence of weak martingale solutions is established for linear noise in effective dimension two (the physical dimension) and for nonlinear noise in dimension one. A key ingredient is the derivation of a-priori estimates on solutions to the equation as well as finding consistent approximations converging to a non-negative (and hence physically meaningful) limit. In the nonlinear noise case, the situation of an almost everywhere positive and non-negative initial value are treated separately. While in the former case the energy of the system can be estimated uniformly in time, the latter case allows only for a control of lower order functionals called α-entropies along the dynamics. While this forces us to rely on a weaker notion of solutions for non-negative initial values, it applies to a larger class of noises driving the equation.
Subsequently, based on stochastic maximal regularity techniques, stochastic thin-film equations are shown to be well-posed for strictly positive initial values in any spatial dimension until the profile touches down or blows up in suitable function spaces. In dimension one, the latter possibility is excluded by establishing a-priori estimates on the solution under the additional consideration of repulsive interaction forces between the molecules of the fluid and the substrate. Consequently, the equation admits unique solutions globally in time in this case for linear and nonlinear gradient noise terms. We also show that these solutions become as smooth as the noise permits. ...
Using stochastic compactness arguments the existence of weak martingale solutions is established for linear noise in effective dimension two (the physical dimension) and for nonlinear noise in dimension one. A key ingredient is the derivation of a-priori estimates on solutions to the equation as well as finding consistent approximations converging to a non-negative (and hence physically meaningful) limit. In the nonlinear noise case, the situation of an almost everywhere positive and non-negative initial value are treated separately. While in the former case the energy of the system can be estimated uniformly in time, the latter case allows only for a control of lower order functionals called α-entropies along the dynamics. While this forces us to rely on a weaker notion of solutions for non-negative initial values, it applies to a larger class of noises driving the equation.
Subsequently, based on stochastic maximal regularity techniques, stochastic thin-film equations are shown to be well-posed for strictly positive initial values in any spatial dimension until the profile touches down or blows up in suitable function spaces. In dimension one, the latter possibility is excluded by establishing a-priori estimates on the solution under the additional consideration of repulsive interaction forces between the molecules of the fluid and the substrate. Consequently, the equation admits unique solutions globally in time in this case for linear and nonlinear gradient noise terms. We also show that these solutions become as smooth as the noise permits.
gence in probability have been widely applied. However, in that case we restrict
ourselves to only a part of all the measurable functions and to an underlying prob-
ability space. One of the main aims of this thesis is to generalize this concept to
the set of all measurable functions with the usual a.e. equivalence classes (which
we call L0) and (possibly) non-finite measure spaces. The other main aim is to
establish an ordered structure on this L0 space ...
gence in probability have been widely applied. However, in that case we restrict
ourselves to only a part of all the measurable functions and to an underlying prob-
ability space. One of the main aims of this thesis is to generalize this concept to
the set of all measurable functions with the usual a.e. equivalence classes (which
we call L0) and (possibly) non-finite measure spaces. The other main aim is to
establish an ordered structure on this L0 space
For such a symbol, we define its Newton polygon $N(P)$ as a certain convex hull of points on $[0,\infty)^2$, which serves as a geometric description of the order of $P$. We define the weight function $W_P$ to be a positive polynomial with the orders found on the vertices of the Newton polygon $N(P)$. We define a notion of order functions, and show the order function of $P$ is $d_\gamma(P)$.
We define a notion of parameter-ellipticity and parabolicity for symbols in $S(L_t\times L_x)$ based on the Newton polygon $N(P)$, namely N-parameter-ellipticity and N-parabolicity. Using various results from the work of R. Denk and M. Kaip \cite{Denk Kaip} and results I introduced myself, we then prove the equivalence between N-parameter-ellipticity and having non-vanishing $\gamma$-principal parts. ...
For such a symbol, we define its Newton polygon $N(P)$ as a certain convex hull of points on $[0,\infty)^2$, which serves as a geometric description of the order of $P$. We define the weight function $W_P$ to be a positive polynomial with the orders found on the vertices of the Newton polygon $N(P)$. We define a notion of order functions, and show the order function of $P$ is $d_\gamma(P)$.
We define a notion of parameter-ellipticity and parabolicity for symbols in $S(L_t\times L_x)$ based on the Newton polygon $N(P)$, namely N-parameter-ellipticity and N-parabolicity. Using various results from the work of R. Denk and M. Kaip \cite{Denk Kaip} and results I introduced myself, we then prove the equivalence between N-parameter-ellipticity and having non-vanishing $\gamma$-principal parts.
In [8] Kwapien proved that every mean zero
function f ∈ L∞[0, 1] we can
write as f = g ◦ T − g for some g ∈ L∞[0, 1] and some measure preserving
transformation T of [0, 1]. However, as was
discovered in [4] there is a gap
in the proof for the case that f is not
continuous. The aim of this bachelor
thesis is filling in that gap in the proof. We
first extend Kwapien’s proof for continuous functions to certain other measure
spaces. Thereafter, we use the method of proof suggested by Kwapien, to proof the theorem for mean zero
function f ∈ L∞[0, 1] for which λ(f−1({x})) = 0 for all x ∈ R.
Using this result we then proof that every mean zero function f ∈ L∞[0, 1] can be written as a sum f =(g1 ◦ T1 − g1) +
(g2 ◦ T2 − g2) where g1, g2
∈ L∞[0, 1]
and where T1, T2 are
measure preserving transformations of [0, 1]. We
finish this thesis with an
application of Kwapien’s theorem in the study to
singular traces
...
In [8] Kwapien proved that every mean zero
function f ∈ L∞[0, 1] we can
write as f = g ◦ T − g for some g ∈ L∞[0, 1] and some measure preserving
transformation T of [0, 1]. However, as was
discovered in [4] there is a gap
in the proof for the case that f is not
continuous. The aim of this bachelor
thesis is filling in that gap in the proof. We
first extend Kwapien’s proof for continuous functions to certain other measure
spaces. Thereafter, we use the method of proof suggested by Kwapien, to proof the theorem for mean zero
function f ∈ L∞[0, 1] for which λ(f−1({x})) = 0 for all x ∈ R.
Using this result we then proof that every mean zero function f ∈ L∞[0, 1] can be written as a sum f =(g1 ◦ T1 − g1) +
(g2 ◦ T2 − g2) where g1, g2
∈ L∞[0, 1]
and where T1, T2 are
measure preserving transformations of [0, 1]. We
finish this thesis with an
application of Kwapien’s theorem in the study to
singular traces
Het Plank Probleem
Het toegankelijk maken van een open probleem