A model reduction approach to numerical inversion for parabolic partial differential equations
|
|
- Charles Jenkins
- 5 years ago
- Views:
Transcription
1 A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann Arbor 2 University of Houston 3 Schlumberger Doll Research, Boston Liliana Borcea (University of Michigan) Model Reduction for Inversion / 29
2 Problem formulation Controlled Source Electromagnetics Inverse problem: determine underground resistivity r(x) from measurements of the magnetic field H(t, x) at the surface. Model: diffusion Maxwell system (negligible displacement current) [ ] r(x) H(x, t) = t H(x, t) Setup: x = (x, z), r := r(x) for x Ω R n, simply connected with boundary B = B A B I. Measurements are on accessible B A. Polarization H(x, t) = u(x, t)e z and the boundary conditions n u(, t) = 0 and BA u(, t) = 0. BI Excitation is in initial conditions and measurements are u(, t) Liliana Borcea (University of Michigan) Model Reduction for Inversion 2 / 29 BA
3 Problem formulation Generic parabolic problem t u(x, t) = (r(x) u(x, t)), x Ω, t > 0 Boundary conditions u(, t) BI = 0, n u(, t) = 0. BA Initial conditions u(x, 0) = u o (x) supported on B A. The sources/receivers are on B A and give knowledge of the measurement operator M(u o ) = u(x, t) for all u o supported in B A and t > 0. BA Inverse problem is to find r(x) given M. Liliana Borcea (University of Michigan) Model Reduction for Inversion 3 / 29
4 Problem formulation Model reduction The measurement operator M is nonlinear in the unknown r(x), but it depends linearly on the initial conditions u o. Suffices to consider integral kernel of mapping from u o to M(u o ). In -D we have Ω = [0, ], with B A = {0} and B I = {} and u(0, t) = u o M(δ(x)) with scalar response y(t) := M(δ). In 2-D B A is a line segment and y ij (t) := M(δ(x j )), x i, x j B A. xi Laplace transform of y(t) is the transfer function of the system. Reduced model rational interpolation of transfer function. Liliana Borcea (University of Michigan) Model Reduction for Inversion 4 / 29
5 Problem formulation Model reduction motivation Benefit for forward problem: Reduced model is small (fast computations) and yet approximates well the transfer function. We may relate the model to a finite difference spatial discretization of the PDE on a special grid, with few points. Benefit for inversion: The nonlinear map R from the function space of the unknown resistivity r(x) to the low-dimensional space of the discrete resistivities may be used to precondition the inversion. We can adapt the size of the reduced model to the noise level (regularization) so we get a natural stability-resolution trade-off. Liliana Borcea (University of Michigan) Model Reduction for Inversion 5 / 29
6 Model order reduction and inversion -D case. Semi-discrete case (for simplicity). Discretize Ω = [0, ] on a fine grid with spacing h = N+ t u(t) = A(r)u(t), u(0) = h e Differential operator (r ) is replaced by the matrix A(r) = D T diag(r)d, r R N + Inverse problem: given y(t; r) = e T u(t) find r RN +. Forward model reduction: we know r and seek a reduced model that approximates y. Inverse model reduction: obtain a reduced model from y, then find r which has this reduced model. Liliana Borcea (University of Michigan) Model Reduction for Inversion 6 / 29
7 Model reduction Model order reduction and inversion Transfer function of full model A(r) R N N, b R N G(s; r) = 0 dt y(t; r)e st = b T (si A(r)) b, s > 0, b = e h Reduced model (A m, b m ) with A m R m m, b m R m and m N m G m (s) = b T m(si m A m ) c j b m = s + θ j with θ j = eigenvalues of A m, and c j = (b T mz j ) 2, where z j = normalized eigenvectors. j= Rational interpolation s k G m (σ j ) = s k G(σ j ) at nodes σ j [0, + ), l for j =,..., l, k = 0,..., 2M j and m = M j. j= Liliana Borcea (University of Michigan) Model Reduction for Inversion 7 / 29
8 Model order reduction and inversion Projection-based model reduction Grimme 997 showed that interpolation is equivalent to projection-based model reduction: A m = V T AV R m m, b m = V T b R m, V R N m has orthonormal columns spanning rational Krylov space } K m (σ) = span {(σ j I A) k b j =,..., l; k =,..., M j The choice of the points σ j matter as we shall see! Liliana Borcea (University of Michigan) Model Reduction for Inversion 8 / 29
9 Model order reduction and inversion Connection to finite differences The matrix A m is dense in general, so how do we connect the reduced model to a finite difference scheme? We have the (Stieltjes) continued fraction representation G m (s) = κ s + κ κ m s + κ m G m (s) = W (s) the response of second-order difference scheme ( Wj+ W j W ) j W j sw j = 0, j =,..., m κ j κ j κ j where r(x) is encoded in {κ j, κ j }. Liliana Borcea (University of Michigan) Model Reduction for Inversion 9 / 29
10 Model order reduction and inversion Optimization formulation Data is d(t) = y(t; r true ) + N (t) where N (t) = noise and discretization error. Reconstruction is r = arg min r R N + 2 Q(d(t)) Q(y(t; r)) 2 2 Instead of minimizing misfit between d(t) and y(t; r) = [F(r)] (t), where F is the forward map, we work with R(r) = Q F(r) where Q is used as a preconditioner. We seek r that gives the reduced model calculated for the measurements d(t). Liliana Borcea (University of Michigan) Model Reduction for Inversion 0 / 29
11 Model order reduction and inversion The map R(r) = Q(y( ; r)) = Q F(r). R : r (a) A(r) (b) V (c) (d) A m {(c j, θ j )} m (e) j= {(log κ j, log κ j )} m j= (a) is definition A(r) = D T diag(r)d. (b) is calculation of orthonormal basis of rational Krylov subspace projection matrix V. (c) is definition A m = V T A(r)V and b m = V T b. (d) here θ m are the eigenvalues of A m and c j = (b T mz j ) 2, where z j are the eigenvectors. (e) Conversion of rational interpolant to continued fraction (Lanczos) and then logarithm of the coefficients. Here all the steps are stable! Liliana Borcea (University of Michigan) Model Reduction for Inversion / 29
12 Model order reduction and inversion The chain of mappings in Q(d(t) Q : d(t) (a) G(s) (b) G m (s) (c) {(c j, θ j )} m (d) j= {(log κ j, log κ j )} m j=. (a) is Laplace transform of measured data. (b) is rational (Padé) interpolation. This is the only ill-conditioned step and we regularize by reducing m until we get positive continued fraction coefficients. (c) Calculate the poles θ j and residues c j of the rational G m (s). (d) Calculate the continued fraction coefficients (Lanczos iteration). Liliana Borcea (University of Michigan) Model Reduction for Inversion 2 / 29
13 Model order reduction and inversion The map R(r) = Q F(r) and discrete resistivities Change coordinates x ξ so that r /2 x = ξ and w(ξ, s) = Laplace transform of u(x, t) satisfies Reduced model corresponds to ( Wj+ W j κ j κ j r /2 ξ (r /2 ξ w) sw = 0 W j W j κ j ) sw j = 0...., m The discrete resistivity is defined by κ 2 j and κ 2 j up to scaling by mesh steps, which we do not need to know. In the optimization we work with the log of κ j and κ j, so grid factors cancel out. R takes r to the (log of ) discrete resistivities. If we interpolated these grid values in Ω, we would expect that R(r) would be close to the identity for smooth r i.e., Q is like inverse of F. Liliana Borcea (University of Michigan) Model Reduction for Inversion 3 / 29
14 Model order reduction and inversion Optimal grids calculated from κ o j and κ o j for r o m = m = Compare three choices of matching conditions: (, ) Moment matching at zero, K m (0) = span { A b, A 2 b,..., A m b } (, ) Interpolation K m (σ) = span { (σ j I A) b j =,..., m } at geometrically spaced nodes σ j = σ ( + 2/m) j (, ) Interpolation at faster growing nodes σ Liliana Borcea (University of Michigan) Model Reduction for Inversion 4 / 29
15 Model order reduction and inversion Conditioning of the Jacobian DR cond(dr) ( ) N ( ) N κ κ Rows and r k k= r k k= are almost collinear when grid points get too close poor conditioning m Condition number growth for moment matching at zero ( ), interpolation at σ ( ), interpolation at σ ( ) Liliana Borcea (University of Michigan) Model Reduction for Inversion 5 / 29
16 Model order reduction and inversion Sensitivity functions κ k r (blue) and κ k r (red) k =,... m = 6 from left to right, top to bottom Liliana Borcea (University of Michigan) Model Reduction for Inversion 6 / 29
17 Reconstruction Regularization and numerical results We solve min r Q(d( )) R(r) 2 with Gauss-Newton iteration. Jacobian DR R 2m N has a large null space 2m N. At each iteration the Gauss-Newton update is in the 2m dimensional range of the pseudoinverse DR. We can improve results with prior info, by adding correction δr that minimizes regularization term P(r + δr) over δr null(dr). E.g. weighted discrete H seminorm P(r) = 2 W/2 Dr 2 2. Liliana Borcea (University of Michigan) Model Reduction for Inversion 7 / 29
18 Regularization and numerical results Numerical experiments: D setup Synthetic noisy data d j = y(t j ; r true ) + N (t j ), j =,..., N T Noise model N = ɛ diag(χ,..., χ NT )y, with independent χ k Gaussian mean zero, variance. Reduced model size m chosen based on noise level ɛ m 3 4 ɛ Initial guess r, five Gauss-Newton iterations n GN = 5. Liliana Borcea (University of Michigan) Model Reduction for Inversion 8 / 29
19 Regularization and numerical results Numerical results: smooth resistivity m = 3, ɛ = m = 4, ɛ = e e True resistivity r true (quadratic) Reconstruction r (2) after one Gauss-Newton iteration Reconstruction r (6) after five Gauss-Newton iterations Liliana Borcea (University of Michigan) Model Reduction for Inversion 9 / 29
20 Regularization and numerical results Numerical results: smooth resistivity m = 3, ɛ = m = 4, ɛ = e e True resistivity r true (linear + Gaussian) Reconstruction r (2) after one Gauss-Newton iteration Reconstruction r (6) after five Gauss-Newton iterations Liliana Borcea (University of Michigan) Model Reduction for Inversion 20 / 29
21 Regularization and numerical results Numerical results: piecewise constant resistivity m = 4, ɛ = m = 5, ɛ = e e True resistivity r true (jump of contrast 2) Reconstruction r (2) after one Gauss-Newton iteration Reconstruction r (6) after five Gauss-Newton iterations Liliana Borcea (University of Michigan) Model Reduction for Inversion 2 / 29
22 2-D case Philosophy -D problem is formally determined: D unknown r(x) and D data y(t). 2-D problem is overdetermined: 2D unknown r(x) but 3D data y ij (t), where (i, j) are source-detector pairs. We proceed as in -D and construct separately reduced models (i.e., maps R ij ) for data subsets indexed by i, j =,..., N d. The maps R ij are coupled by their dependence on the same r(x) and are all used in inversion. Liliana Borcea (University of Michigan) Model Reduction for Inversion 22 / 29
23 2-D case Semi-discretization on fine grid with N points For the j source in B A (a line on the surface) we have t u (j) (t) = A(r)u (j) (t), u (j) (0) = b (j). The response matrix y ij (t) is the restriction of the solution u (j) (t) = e A(r)t b (j) to the support of the i th receiver y ij (t) = b (i)t e A(r)t b (j). Its Laplace transform is the transfer function G (ij) (s). The reduced model is (A (ij) m, b (ij) m ) and its transfer function is ( ) G (ij) m (s) = b (i)t m si m A (ij) (j) m b m = m q= c (ij) q s + θ (ij) q Liliana Borcea (University of Michigan) Model Reduction for Inversion 23 / 29
24 2-D case Model reduction We can do this for any (i, j), but we are interested in using the continued fraction representation of G (ij) m (s). {( Coefficients κ (ij) l )}, κ (ij) m l c (ij) q = l= (b (i)t m are guaranteed to be positive if ) ) (b (j)t > 0 z (ij) q m z (ij) q Thus, we work only with the diagonal G (jj) m (s). We have as before R j (r) = Q ( y jj ( ; r) ) ( = log κ (jj) l, log κ (jj) l We solve optimization r = arg min r R N + 2 N d Q(d j (t)) R j (r)) 2 2. j= ) m l= Liliana Borcea (University of Michigan) Model Reduction for Inversion 24 / 29
25 2-D case Sensitivity functions κjj k r (left) and κjj k r (right) k =,... m = 5 (top to bottom) for a single source/receiver. Liliana Borcea (University of Michigan) Model Reduction for Inversion 25 / 29
26 2-D case Reconstructions: two checkerboard inclusions Top: true r(x, y). Bottom: reconstruction after a single Gauss-Newton iteration. Constant initial guess r 0 (x, y). Sensors:. Liliana Borcea (University of Michigan) Model Reduction for Inversion 26 / 29
27 2-D case Reconstructions: two layered inclusions Top: true r(x, y). Bottom: reconstruction after a single Gauss-Newton iteration. Constant initial guess r 0 (x, y). Sensors:. Liliana Borcea (University of Michigan) Model Reduction for Inversion 27 / 29
28 2-D case Reconstructions: single skewed inclusion Top: true r(x, y). Bottom: reconstruction after a single Gauss-Newton iteration. Constant initial guess r 0 (x, y). Sensors:. Liliana Borcea (University of Michigan) Model Reduction for Inversion 28 / 29
29 Conclusion End notes All the details are in: A model reduction approach to numerical inversion for a parabolic partial differential equation. L. Borcea, V. Druskin, A.V. Mamonov and M. Zaslavsky, Inverse Problems 30 (204), 250 (33pp). Still plenty left to be done. The higher dimensional problem treatment is not optimal. Dealing with one source and multiple receivers means doing matrix rational approximation... 3-D problems? Liliana Borcea (University of Michigan) Model Reduction for Inversion 29 / 29
A model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 2, Mikhail Zaslavsky 2 University of Michigan, Ann
More informationStatistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography
Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography Aku Seppänen Inverse Problems Group Department of Applied Physics University of Eastern Finland
More informationReduced Basis Methods for MREIT
Reduced Basis Methods for MREIT Dominik Garmatter garmatter@math.uni-frankfurt.de Group for Numerics of Partial Differential Equations, Goethe University Frankfurt, Germany Joint work with Bastian Harrach
More informationMagnet Resonance Electrical Impedance Tomography (MREIT): convergence of the Harmonic B z Algorithm
Magnet Resonance Electrical Impedance Tomography (MREIT): convergence of the Harmonic B z Algorithm Dominik Garmatter garmatter@math.uni-frankfurt.de Group for Numerics of PDEs, Goethe University Frankfurt,
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationUsing radial basis functions for option pricing
Using radial basis functions for option pricing Elisabeth Larsson Division of Scientific Computing Department of Information Technology Uppsala University Actuarial Mathematics Workshop, March 19, 2013,
More informationExact shape-reconstruction by one-step linearization in EIT
Exact shape-reconstruction by one-step linearization in EIT Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München, Germany Joint work with Jin Keun Seo, Yonsei
More informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationNumerical Solution of Two Asset Jump Diffusion Models for Option Valuation
Numerical Solution of Two Asset Jump Diffusion Models for Option Valuation Simon S. Clift and Peter A. Forsyth Original: December 5, 2005 Revised: January 31, 2007 Abstract Under the assumption that two
More informationHeinz W. Engl. Industrial Mathematics Institute Johannes Kepler Universität Linz, Austria
Some Identification Problems in Finance Heinz W. Engl Industrial Mathematics Institute Johannes Kepler Universität Linz, Austria www.indmath.uni-linz.ac.at Johann Radon Institute for Computational and
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationOn the Ross recovery under the single-factor spot rate model
.... On the Ross recovery under the single-factor spot rate model M. Kijima Tokyo Metropolitan University 11/08/2016 Kijima (TMU) Ross Recovery SMU @ August 11, 2016 1 / 35 Plan of My Talk..1 Introduction:
More informationFinite Element Method
In Finite Difference Methods: the solution domain is divided into a grid of discrete points or nodes the PDE is then written for each node and its derivatives replaced by finite-divided differences In
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Program Advisor: Howard Elman Department of Computer Science May 5, 2016 Tengfei
More information9.1 Principal Component Analysis for Portfolios
Chapter 9 Alpha Trading By the name of the strategies, an alpha trading strategy is to select and trade portfolios so the alpha is maximized. Two important mathematical objects are factor analysis and
More informationPartitioned Analysis of Coupled Systems
Partitioned Analysis of Coupled Systems Hermann G. Matthies, Rainer Niekamp, Jan Steindorf Technische Universität Braunschweig Brunswick, Germany wire@tu-bs.de http://www.wire.tu-bs.de Coupled Problems
More informationHIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR 1D PARABOLIC EQUATIONS. Ahmet İzmirlioğlu. BS, University of Pittsburgh, 2004
HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR D PARABOLIC EQUATIONS by Ahmet İzmirlioğlu BS, University of Pittsburgh, 24 Submitted to the Graduate Faculty of Art and Sciences in partial fulfillment of
More informationTrust Region Methods for Unconstrained Optimisation
Trust Region Methods for Unconstrained Optimisation Lecture 9, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Trust
More informationDirect Methods for linear systems Ax = b basic point: easy to solve triangular systems
NLA p.1/13 Direct Methods for linear systems Ax = b basic point: easy to solve triangular systems... 0 0 0 etc. a n 1,n 1 x n 1 = b n 1 a n 1,n x n solve a n,n x n = b n then back substitution: takes n
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationModelling, Estimation and Hedging of Longevity Risk
IA BE Summer School 2016, K. Antonio, UvA 1 / 50 Modelling, Estimation and Hedging of Longevity Risk Katrien Antonio KU Leuven and University of Amsterdam IA BE Summer School 2016, Leuven Module II: Fitting
More informationSparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models
Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationIntroduction to Numerical PDEs
Introduction to Numerical PDEs Varun Shankar February 16, 2016 1 Introduction In this chapter, we will introduce a general classification scheme for linear second-order PDEs, and discuss when they have
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
More informationExact shape-reconstruction by one-step linearization in EIT
Exact shape-reconstruction by one-step linearization in EIT Bastian von Harrach harrach@math.uni-mainz.de Zentrum Mathematik, M1, Technische Universität München, Germany Joint work with Jin Keun Seo, Yonsei
More informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More informationPICOF, Palaiseau, April 2-4, 2012
The Sobolev gradient regularization strategy for optical tomography coupled with a finite element formulation of the radiative transfer equation Fabien Dubot, Olivier Balima, Yann Favennec, Daniel Rousse
More informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More informationParameter estimation in SDE:s
Lund University Faculty of Engineering Statistics in Finance Centre for Mathematical Sciences, Mathematical Statistics HT 2011 Parameter estimation in SDE:s This computer exercise concerns some estimation
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationPDE Project Course 1. Adaptive finite element methods
PDE Project Course 1. Adaptive finite element methods Anders Logg logg@math.chalmers.se Department of Computational Mathematics PDE Project Course 03/04 p. 1 Lecture plan Introduction to FEM FEM for Poisson
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
More informationResearch Article Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation
Applied Mathematics Volume 1, Article ID 796814, 1 pages doi:11155/1/796814 Research Article Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation Zhongdi
More informationInfinite Reload Options: Pricing and Analysis
Infinite Reload Options: Pricing and Analysis A. C. Bélanger P. A. Forsyth April 27, 2006 Abstract Infinite reload options allow the user to exercise his reload right as often as he chooses during the
More informationIs Greedy Coordinate Descent a Terrible Algorithm?
Is Greedy Coordinate Descent a Terrible Algorithm? Julie Nutini, Mark Schmidt, Issam Laradji, Michael Friedlander, Hoyt Koepke University of British Columbia Optimization and Big Data, 2015 Context: Random
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationAdvanced Numerical Methods for Financial Problems
Advanced Numerical Methods for Financial Problems Pricing of Derivatives Krasimir Milanov krasimir.milanov@finanalytica.com Department of Research and Development FinAnalytica Ltd. Seminar: Signal Analysis
More informationUnit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform
The Laplace Transform Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform Engineering 5821: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationAmerican Options; an American delayed- Exercise model and the free boundary. Business Analytics Paper. Nadra Abdalla
American Options; an American delayed- Exercise model and the free boundary Business Analytics Paper Nadra Abdalla [Geef tekst op] Pagina 1 Business Analytics Paper VU University Amsterdam Faculty of Sciences
More informationOptimal investments under dynamic performance critria. Lecture IV
Optimal investments under dynamic performance critria Lecture IV 1 Utility-based measurement of performance 2 Deterministic environment Utility traits u(x, t) : x wealth and t time Monotonicity u x (x,
More informationFX Smile Modelling. 9 September September 9, 2008
FX Smile Modelling 9 September 008 September 9, 008 Contents 1 FX Implied Volatility 1 Interpolation.1 Parametrisation............................. Pure Interpolation.......................... Abstract
More informationNeuro-Dynamic Programming for Fractionated Radiotherapy Planning
Neuro-Dynamic Programming for Fractionated Radiotherapy Planning Geng Deng Michael C. Ferris University of Wisconsin at Madison Conference on Optimization and Health Care, Feb, 2006 Background Optimal
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationInterest Rate Curves Calibration with Monte-Carlo Simulatio
Interest Rate Curves Calibration with Monte-Carlo Simulation 24 june 2008 Participants A. Baena (UCM) Y. Borhani (Univ. of Oxford) E. Leoncini (Univ. of Florence) R. Minguez (UCM) J.M. Nkhaso (UCM) A.
More informationAn adaptive cubic regularization algorithm for nonconvex optimization with convex constraints and its function-evaluation complexity
An adaptive cubic regularization algorithm for nonconvex optimization with convex constraints and its function-evaluation complexity Coralia Cartis, Nick Gould and Philippe Toint Department of Mathematics,
More informationCalibration of local volatility surfaces under PDE constraints
Calibration of local volatility surfaces under PDE constraints Love Lindholm Abstract he calibration of a local volatility surface to option market prices is an inverse problem that is ill-posed as a result
More informationA local RBF method based on a finite collocation approach
Boundary Elements and Other Mesh Reduction Methods XXXVIII 73 A local RBF method based on a finite collocation approach D. Stevens & H. Power Department of Mechanical Materials and Manufacturing Engineering,
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationIMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS
IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More information1 Explicit Euler Scheme (or Euler Forward Scheme )
Numerical methods for PDE in Finance - M2MO - Paris Diderot American options January 2018 Files: https://ljll.math.upmc.fr/bokanowski/enseignement/2017/m2mo/m2mo.html We look for a numerical approximation
More informationEE/AA 578 Univ. of Washington, Fall Homework 8
EE/AA 578 Univ. of Washington, Fall 2016 Homework 8 1. Multi-label SVM. The basic Support Vector Machine (SVM) described in the lecture (and textbook) is used for classification of data with two labels.
More informationCalibration Lecture 1: Background and Parametric Models
Calibration Lecture 1: Background and Parametric Models March 2016 Motivation What is calibration? Derivative pricing models depend on parameters: Black-Scholes σ, interest rate r, Heston reversion speed
More informationSlides for DN2281, KTH 1
Slides for DN2281, KTH 1 January 28, 2014 1 Based on the lecture notes Stochastic and Partial Differential Equations with Adapted Numerics, by J. Carlsson, K.-S. Moon, A. Szepessy, R. Tempone, G. Zouraris.
More informationDiscussion Paper No. DP 07/05
SCHOOL OF ACCOUNTING, FINANCE AND MANAGEMENT Essex Finance Centre A Stochastic Variance Factor Model for Large Datasets and an Application to S&P data A. Cipollini University of Essex G. Kapetanios Queen
More information1.1 Forms for fractions px + q An expression of the form (x + r) (x + s) quadratic expression which factorises) may be written as
1 Partial Fractions x 2 + 1 ny rational expression e.g. x (x 2 1) or x 4 x may be written () (x 3) as a sum of simpler fractions. This has uses in many areas e.g. integration or Laplace Transforms. The
More informationReduced models for sparse grid discretizations of the multi-asset Black-Scholes equation
Reduced models for sparse grid discretizations of the multi-asset Black-Scholes equation The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationBACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin
BACHELIER FINANCE SOCIETY 4 th World Congress Tokyo, 26 Investments and forward utilities Thaleia Zariphopoulou The University of Texas at Austin 1 Topics Utility-based measurement of performance Utilities
More informationRisk minimizing strategies for tracking a stochastic target
Risk minimizing strategies for tracking a stochastic target Andrzej Palczewski Abstract We consider a stochastic control problem of beating a stochastic benchmark. The problem is considered in an incomplete
More informationSome mathematical results for Black-Scholes-type equations for financial derivatives
1 Some mathematical results for Black-Scholes-type equations for financial derivatives Ansgar Jüngel Vienna University of Technology www.jungel.at.vu (joint work with Bertram Düring, University of Mainz)
More informationOptimal order execution
Optimal order execution Jim Gatheral (including joint work with Alexander Schied and Alla Slynko) Thalesian Seminar, New York, June 14, 211 References [Almgren] Robert Almgren, Equity market impact, Risk
More informationThe Forward Kolmogorov Equation for Two Dimensional Options
The Forward Kolmogorov Equation for Two Dimensional Options Antoine Conze (Nexgenfs bank), Nicolas Lantos (Nexgenfs bank and UPMC), Olivier Pironneau (LJLL, University of Paris VI) March, 04 Abstract Pricing
More informationPackage multiassetoptions
Package multiassetoptions February 20, 2015 Type Package Title Finite Difference Method for Multi-Asset Option Valuation Version 0.1-1 Date 2015-01-31 Author Maintainer Michael Eichenberger
More informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
More informationStable Local Volatility Function Calibration Using Spline Kernel
Stable Local Volatility Function Calibration Using Spline Kernel Thomas F. Coleman Yuying Li Cheng Wang January 25, 213 Abstract We propose an optimization formulation using the l 1 norm to ensure accuracy
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationTo apply SP models we need to generate scenarios which represent the uncertainty IN A SENSIBLE WAY, taking into account
Scenario Generation To apply SP models we need to generate scenarios which represent the uncertainty IN A SENSIBLE WAY, taking into account the goal of the model and its structure, the available information,
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationWeak Reflection Principle and Static Hedging of Barrier Options
Weak Reflection Principle and Static Hedging of Barrier Options Sergey Nadtochiy Department of Mathematics University of Michigan Apr 2013 Fields Quantitative Finance Seminar Fields Institute, Toronto
More informationGlobal Currency Hedging
Global Currency Hedging JOHN Y. CAMPBELL, KARINE SERFATY-DE MEDEIROS, and LUIS M. VICEIRA ABSTRACT Over the period 1975 to 2005, the U.S. dollar (particularly in relation to the Canadian dollar), the euro,
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has
More informationJUMPS WITHOUT TEARS: A NEW SPLITTING TECHNOLOGY FOR BARRIER OPTIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages c Institute for Scientific Computing and Information JUMPS WITHOUT TEARS: A NEW SPLITTING TECHNOLOGY FOR BARRIER OPTIONS ANDREY
More informationInt. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p approach
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS001) p.5901 What drives short rate dynamics? approach A functional gradient descent Audrino, Francesco University
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationA way to improve incremental 2-norm condition estimation
A way to improve incremental 2-norm condition estimation Jurjen Duintjer Tebbens Institute of Computer Science Academy of Sciences of the Czech Republic duintjertebbens@cs.cas.cz Miroslav Tůma Institute
More informationConvex-Cardinality Problems Part II
l 1 -norm Methods for Convex-Cardinality Problems Part II total variation iterated weighted l 1 heuristic matrix rank constraints Prof. S. Boyd, EE364b, Stanford University Total variation reconstruction
More informationGraph signal processing for clustering
Graph signal processing for clustering Nicolas Tremblay PANAMA Team, INRIA Rennes with Rémi Gribonval, Signal Processing Laboratory 2, EPFL, Lausanne with Pierre Vandergheynst. What s clustering? N. Tremblay
More informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
More informationPricing Algorithms for financial derivatives on baskets modeled by Lévy copulas
Pricing Algorithms for financial derivatives on baskets modeled by Lévy copulas Christoph Winter, ETH Zurich, Seminar for Applied Mathematics École Polytechnique, Paris, September 6 8, 26 Introduction
More informationarxiv: v1 [math.st] 6 Jun 2014
Strong noise estimation in cubic splines A. Dermoune a, A. El Kaabouchi b arxiv:1406.1629v1 [math.st] 6 Jun 2014 a Laboratoire Paul Painlevé, USTL-UMR-CNRS 8524. UFR de Mathématiques, Bât. M2, 59655 Villeneuve
More informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationarxiv: v1 [q-fin.cp] 1 Nov 2016
Essentially high-order compact schemes with application to stochastic volatility models on non-uniform grids arxiv:1611.00316v1 [q-fin.cp] 1 Nov 016 Bertram Düring Christof Heuer November, 016 Abstract
More informationA Monte Carlo Based Analysis of Optimal Design Criteria
A Monte Carlo Based Analysis of Optimal Design Criteria H. T. Banks, Kathleen J. Holm and Franz Kappel Center for Quantitative Sciences in Biomedicine Center for Research in Scientific Computation North
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationPricing and Hedging of Credit Derivatives via Nonlinear Filtering
Pricing and Hedging of Credit Derivatives via Nonlinear Filtering Rüdiger Frey Universität Leipzig May 2008 ruediger.frey@math.uni-leipzig.de www.math.uni-leipzig.de/~frey based on work with T. Schmidt,
More informationPricing with a Smile. Bruno Dupire. Bloomberg
CP-Bruno Dupire.qxd 10/08/04 6:38 PM Page 1 11 Pricing with a Smile Bruno Dupire Bloomberg The Black Scholes model (see Black and Scholes, 1973) gives options prices as a function of volatility. If an
More informationGalerkin Least Square FEM for the European option price with CEV model
Galerkin Least Square FEM for the European option price with CEV model A Major Qualifying Project Submitted to the Faculty of Worcester Polytechnic Institute In partial fulfillment of requirements for
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationCS227-Scientific Computing. Lecture 6: Nonlinear Equations
CS227-Scientific Computing Lecture 6: Nonlinear Equations A Financial Problem You invest $100 a month in an interest-bearing account. You make 60 deposits, and one month after the last deposit (5 years
More informationVladimir Spokoiny (joint with J.Polzehl) Varying coefficient GARCH versus local constant volatility modeling.
W e ie rstra ß -In stitu t fü r A n g e w a n d te A n a ly sis u n d S to c h a stik STATDEP 2005 Vladimir Spokoiny (joint with J.Polzehl) Varying coefficient GARCH versus local constant volatility modeling.
More information