On multivariate Multi-Resolution Analysis, using generalized (non homogeneous) polyharmonic splines. or: A way for deriving RBF and associated MRA
|
|
- Hilary Rose
- 5 years ago
- Views:
Transcription
1 MAIA conference Erice (Italy), September 6, 3 On multivariate Multi-Resolution Analysis, using generalized (non homogeneous) polyharmonic splines or: A way for deriving RBF and associated MRA Christophe Rabut, Mira Bozzini and Milvia Rossini University of Toulouse (INSA ; IMT, IREM, MAIAA), France. Some known tools... using Fourier Transform. Extension of polyharmonic splines, and associated MRA
2 Part : Some known tools... using Fourier Transform
3 A word on (odd degree) polynomial splines Definition: σ m = Argmin i [:n], f(xi )=y i IR (f (m) (x)) dx = σ m (x) = λ i x x i m + p m (x) i=:n with q IP m, λ iq(x i ) = and p m IP m i=:n Radial basis functions : writing u m (x) = (m!) x m σ m (x) = µ i u m (x x i ) + q m (x) i=:n (same constraints on µ and q m ) Derivative and Fourier Transform: E(u m ) := ( ) m u (m) m = Dirac û m (ω) = ω m 3
4 ... and on associated functions Basic spline, or B-spline : ϕ m = ( ) m δ m u m ϕm (ω) = sin ω m ω m Lagrangian spline, or L-spline : L m () = ; j ZZ\{}, L m (j) = L m = j Z L m (ω) = L m (j/) L m ( j) ; ϕ m = l Z ϕ m (ω) ϕ m (ω πl) = j Z ϕm (ω) j Z ϕ m (j) exp( i j ω) = ϕ m (j) L m ( j) l Z ω m (ω πl) m Semi-orthogonal wavelet: ψ m = ( ) m D m L m (to be normalized) ψ m (ω) = ω m L m (ω) = l Z ω m (ω πl) 4m Orthogonal wavelet: ψ (ω) = ( l Z ψ(ω) ψ( ω πl )) = ω m (ω πl) 4m l Z 4
5 Linear case B-spline L-spline psi-spline function B ; norm(b) =.48 function L ; norm(l) =.4 function psi ; norm(psi) = κ = ( ) ; γ = ( ) ; κ = ( ) ; γ = ( ) ; κ = ( ) ; γ = ( ) ; psi-ortho psi-hat psi-ortho-hat 3.5 function psi ortho ; norm(psi ortho ) = function psi hat function psi ortho hat κ = ( ) ; γ = ( ) ; 3 3 κ = ( ) ; γ = ( ) ; 3 3 κ = ( ) ; γ = ( ) ; 5
6 6
7 Cubic case B-spline L-spline psi-spline function B ; norm(b) =.346 function L ; norm(l) =.467 function psi ; norm(psi) = κ = ( ) ; γ = ( ) ; κ = ( ) ; γ = ( ) ; κ = ( ) ; γ = ( ) ; psi-ortho psi-hat psi-ortho-hat 3 function psi ortho ; norm(psi ortho ) = 6 function psi hat function psi ortho hat κ = ( ) ; γ = ( ) ; 3 3 κ = ( ) ; γ = ( ) ; 3 3 κ = ( ) ; γ = ( ) ; 7
8 8
9 Use of the associated functions (translation invariant spaces) ŷ is defined by ŷ(ω) = j Z y j exp( i j ω). B-spline approximation of vector y : (or of points P j for B-spline curve) σ m = j Z y j ϕ m ( j) σ = ŷ ϕ m. Interpolating spline of vector y : (or P j instead of y j for interpolating spline curve) σ m = y j L m ( j) j Z σ = ŷ L. Wavelet decomposition of some f L (IR): f = j Z l Z (f, ψ ( l j)) ψ ( l j) (f, ψ ( l j)) = l exp(i l j) ( f, ψ ( l ) 9
10 Spline under tension Definition and Fourier transform σ m,k = Argmin i Z, f(xi )=y i IR (f (m) (x)) dx+ρ IR (f (k) (x)) dx (say k < m) E ρ (u) := ( ) m D m u + ( ) k ρ D k u = Dirac ; û ρ (ω) = ω m + ρ ω k ϕ ρ (ω) = sinm ω + ρ sin k ω ω m + ρ ω k L ρ (ω) = ψ ρ (ω) = ω m L ρ(ω) + ρ ω k L ρ(ω) = l Z (ω m + ρ ω k ) ( (ω πl) m + ρ (ω πl) k) ( ω m + ρ ω k) ( (ω πl) m + ρ (ω πl) k) ψρ (ω) = l Z l Z ( ω m + ρ ω k) ((ω πl) m + ρ (ω πl) k )
11 function B ; norm(b) = κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; Tension splines (linear-cubic) B-spline L-spline psi-spline function L ; norm(l) = κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; function psi ; norm(psi) = κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; psi-ortho psi-hat psi-ortho-hat function psi ortho ; norm(psi ortho ) = 8 function psi hat function psi ortho hat κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ;
12
13 Tension splines (cubic-quintic) B-spline L-spline psi-spline function B ; norm(b) = function L ; norm(l) = function psi ; norm(psi) = κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; psi-ortho psi-hat psi-ortho-hat 3 function psi ortho ; norm(psi ortho ) = 7 x 4 function psi hat.4 function psi ortho hat κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( 3 ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 3
14 4
15 Definition and Fourier transform Fractional splines Let s [..) and m IN such that α := m + s > d/ σ α = Argmin i [:n], f(xi )=y i IR (f (α) (x)) dx = Argmin i [:n], f(xi )=y i IR (ω ) s ( F(f (m) (ω) ) dω E α (u α ) := ( ) α D α u α = Dirac ; û α (ω) = (ω ) α u α (x) = c α x α u α (x) = c α x α ln x if α is not an even integer number. if α is an even integer number. (c α is some known real valued constant) They too are in some place between order (linear) and order (cubic) splines, but are different from splines under tension. B-spline : ϕα (ω) = sin ω ω α 5
16 Fractional splines (linear-cubic) B-spline L-spline psi-spline function B ; norm(b) = function L ; norm(l) = function psi ; norm(psi) = κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; psi-ortho psi-hat psi-ortho-hat 3.5 function psi ortho ; norm(psi ortho ) = 6 function psi hat function psi ortho hat κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; 3 3 κ = ( ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; (.6 ) ; ( ) ; 6
17 7
18 comparison tension B-splines versus fractional B-splines.8 function B ; norm(b) = κ = (.84 ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; ( ) ; (.6 ) ; 8
19 comparison tension L-splines versus fractional L-splines function L ; norm(l) = κ = (.84 ) ; ( ) ; ( ) ; ( ) ; γ = ( ) ; (.3 ) ; ( ) ; (.6 ) ; 9
20 (d-dimensional, real order) polyharmonic splines Definition and expression Let s [..) and m IN such that α := m + s > / σ α = Argmin i [:n], f(xi )=y i IR d D α f(x) dx = Argmin i [:n], f(xi )=y i IR d ω s F(D m f)(ω) dx = σ α (x) = λ i u α (x x i ) + p m (x) i=:n with q IP m, where E α (u α ) := ( ) m α u α = Dirac i=:n λ iq(x i ) = and p m IP m u α (x) = c α x α d if α d is not an even integer number. = c α x α d ln x if α d is an even integer number. Radial basis functions (no other known writing)
21 (d-dimensional, real order) polyharmonic splines Derivative and Fourier Transform α u α = Dirac û α (ω) = ω α Polyharmonic B-spline: ϕ α = α u α ϕα (ω) = sin ω ω α Lagrangian spline, or L-spline : L α (ω) = ω α ( ω πl ) α l Z d Wavelets: ψ α = α L α ψ α (ω) = ( ω ) α L α (ω) = ( ω ) α ( ω πl ) α ψ α (ω) = ( ω ) α (( ω πl ) α l Z d l Z d
22 Biharmonic ( thin plate ) B-spline Elementary Thin Plate B Spline 3 Elementary Thin Plate B Spline
23 Biharmonic ( thin plate ) L and ψ-spline Lagrangean function polyharmonic wavelet
24 Biharmonic ( thin plate ) ψ and ψ -spline 4
25 (d-dimensional) scaled Whittle-Matérn-Sobolev kernel Idea Regularize the operator m of polyharmonic splines via a regularization of the differential operator: instead of ( ) m, use the operator E = ( + κ I) m, κ being some real number( scaled is for κ ). Its Fourier transform is clearly E(ω) = ( ω + κ ) m. Associated kernel (u function) u m,κ (x) = m (m )! κd/ m x m d/ K m d/ (κ x ) where K m d/ is the order m d/ Bessel function of the third kind (note the mythic m d/ constraint, actually here for ensuring u m being a continuous function and the differential operator having a meaning). Same type of functions generated. 5
26 A generalisation of Whittle-Matérn kernel (/) Mira Bozzini, Milvia Rossini, Robert Schaback, Idea: Mix different values of κ for each scaled iterated Laplacian operator: use E = m j= ( + κ ji) instead of ( + κ I) m, We get: E κ (ω) = m j= ( ω + κ j) The fundamental solution meets ûκ(ω) = j= ( ω + κ j) and u κ is a continuous radial kernel, and is a convolution product of the Kernels u,κj whose Fourier transform is ( ω + κ j ) (no computable by this way!). However we have the following explicit value u κ (x) = m ( ) m [κ /,..., κ x m/] z z m d/ To prove it we use the m-th divided difference relation : ( ) m m j= (s + t j) = [t,..., t m ] z (s + z) 6 K d/ ( x z).
27 Form of some u-functions 7 κ = ( ), (3 7), ( 3 4)
28 A generalisation of Whittle-Matérn kernel (/) Associated B-spline Lagrangian function ϕ κ (ω) = L κ (ω) = m j= l Z d sin ω + κ j ω + κ j m j= ( ω + κ j) m j= ( ω πl + κ j) Wavelets: ψ κ = ψ κ (ω) = ψ κ (ω) = m m j= ( + κ j I) j= ( ω + κ j) m m ( l Z d j= ( ω + κ j) L [κ κ] L [κ κ] (ω) = j= ( ω πl + κ j ) l Z d m j= ( ω + κ j) m j= ( ω πl + κ j 8
29 Proof of the expression of the pre-wavelet Theorem Let E = m ( + j= κ ji). Let L be the Lagrangian E -function, and let ψ = E(L ) (so ψ(ω) = Ê(ω) L (ω)). Let ψ e = ψ( +e) = (E(L ))( +e) where e {, } d \ d. Then ψ e is orthogonal to any cardinal E-function. So ψ e is a pre-wavelet (semi-orthogonal wavelet). Proof Let σ = j Z d λ j u( j) be a cardinal E-function (u is such that E(u) = Dirac). Then : (σ, ψ e ) = (σ, E(L ( +e) ) = (E(σ), L ( +e) ) = ( j Z d λ j E(u( j)), L ( +e) ) = ( j Z d λ j Dirac j, L ( +e) ) = j Z d λ j L (j + e) =. Actually an even more general property (quite general E, and also for scattered data Lagrangean function L and pre-wavelet ψ = E(L )). 9
30 Part : A proposal for a global extension of polyharmonic splines: functions generated by the differential operator E = m i= ( + κ ji) α j 3
31 The proposed generalization of polyharmonic splines Idea: We want both tension and continuous choice of the order. So we choose polyharmonic splines - regularized by the (extended) Whittle-Matérn coefficients (for tension), - generalized by real exponents, (for continuous choice of the order and of the shape). So : E κ,α = m i= ( + κ ji) α j Since the differential operator is a mix of the ( + κ ji) α j, the obtained functions will be a kind of mix of the fractional regularized polyharmonic splines. We get (with the condition m α j > d/) j= E κ,α (ω) = m i= ( ω + κ j) α j The fundamental u solution is a radial function which meets û κ,α (ω) = m j= ( ω + κ j) α j 3
32 Associated B-function Lagrangian function Wavelets: ψ κ,α = ψ κ,α (ω) = ψκ,α(ω) = m m Associated functions ϕ κ,α (ω) = L κ,α (ω) = E κ,α (sin ω) E κ,α (ω) l Z d j= ( + κ j I) α j j= ( ω + κ j) α j l Z (Ê(ω)) d Ê(ω πl) = m = m j= ( sin ω + κ j) α j m j= ( ω + κ j) α j m i= ( ω + κ j) α j j= ( ω πl + κ j) α j L [κ κ;α α] L [κ κ α α] (ω) = l Z d m l Z d m j= ( ω + κ j) α j m j= ( ω πl + κ j) α j j= ( ω + κ j) α j m j= ( ω πl + κ j ) α j 3
33 Examples of B- and L- functions B-function L-function κ = ( ) κ = (3 7) 33
34 34
35 Examples of ψ- and ψ functions ψ-function ψ function κ = ( ) κ = (3 7) 35
36 Minimization property, when all α j are integer Let γ be such that m γ i x i = m (x + i= j= κ j). γ i is positive and is the sum of all product of m i different κ j. In particular : γ = m j= κ j, γ = m i= j i κ j, γ m = m j= κ j, γ m = ; γ j >. Let H κ : (f, g) κ = IR d m j= γ j (D j f D j g) and f κ = IR d m j= γ j D j f Theorem Let A be a set of elements in IR d, and (y a ) a A some associated real numbers. If A is an infinite set, suppose furthermore that there exists f H κ such that a A, f(a) = y a. Then the set of all functions f in H κ meeting a A, f(a) = y a has a unique element, denoted σ A,κ with minimal (semi-)norm f κ. Besides, there exist real numbers (λ a ) a A such that σ A,y meets the relation E(σ A,y ) = λ aδ a, so σ A,y = λ au κ (x a) + p l, a A a A where l is the number of κ j being and p l is some polynomial in IP l. 36
37 Role and influence of κ : tension coefficients Theorem (direct consequence of previous theorem) Let σ A be in the form σ A = a A λ au κ ( a) + p l (), and let y = (y a ) a A such that y a = σ A (a). Then σ A minimizes f κ on all the functions f H κ such that a A, f(a) = y a. Discussion m Since f κ = IR d γ j D j f, we see that any function in the form () j= minimizes a quantity which is a linear combination of various polyharmonic seminorms, the weights of the linear combination being the γ j, which are connected to the κ j by the relation m γ i x i = m (x + i= j= κ j). Since it is known that a polyharmonic spline is all the more oscillating than its order is high, we can interpret γ j as tension coefficients, giving all the more tension than j is small. Intuition of the impact of a vector κ is quite easy when all but one κ j are zero, but is not so easy otherwise. 37
38 Role and influence of real valued α j s Very easy when only one α j is non integer and when the associated κ j is zero : we have real order polyharmonic splines under tension, with tension parameters as explained above for integer α j s. As said before, we know that real order polyharmonic splines are somewhere between the polyharmonic spline of the integer part of the order and of the integer part plus. Otherwise, we need work directly on the Fourier transform of E, and the minimized quantity is less intuitive. There is still work to do to improve intuition about that. 38
39 Two remarks On B-functions Note that for (one dimensional) polynomial splines, the B-spline is non-negative. This is NOT true for many other B-functions, as some tension B-splines, polyharmonic B-splines and all their d-dimensional extensions. We claim that the important property is NOT non-negativity, but positive definite functions, as they all are since ϕ is non-negative (in all extensions we showed here). On the strategy Starting on the differential operator and/or its Fourier transform is an interesting way to work on, and gives other possibilities, depending on what we want to obtain. An operator involving iterated Laplacian operator will give radial basis functions. Of course appropriate properties for having a MRA have to be checked. For example if we want something between quintic and linear we can use κ = ( ρ ρ ) (which gives E = D 6 + (ρ + ρ )D 4 + ρ ρ D in one dimension) and so get Ê(ω) = ω6 + (ρ + ρ )ω 4 ρ ρ ω (or equivalent with and ω in many dimensions), but we also can try to minimize IR (f (3) (x)) +ρ (f (x)) dx, which gives E = D 6 + ρ D and Ê(ω) = ω6 ρ ω (or equivalent with and ω in many dimensions), which is not included in the presented extension (not done!). 39
40 Computational considerations Extensive use of the fast Fourier transform All the functions are given by their Fourier transform. The whole computation is done via FFT and IFFT. For the basis functions (ϕ, ψ, ψ, L: sampling the Fourier transform of the function, and then IFFT of the so-obtained (d-dimensional) vector. the longest computation is for computing the sample of the Fourier transform of the functions. For spline functions in the form y j B( j h), we compute ŷ(ω) (via the FFT), i Z d and then we compute IFFT(ŷ B). Note that this is no more complicated for any type of extension presented here. As usual, there are some side effects for boundly supported vectors or functions. 4
41 Extension to scattered data analysis Surprisingly this is possible, still by using this type of kernel. However everything is not yet done. The main difficulty is to choose the appropriate points close to the point where we want to discretize the derivatives. But we can define the basis functions (one for each centre!) via the Fourier transform, and go on pretty well. However this need still work theoretical as well as numerical... Besides there are some true normalization problems 4
42 A classical test : Lena We use the associated filters on one Lena s eye : 4
43 and obtain : κ = ( ) κ = (3 7) 43
44 Enjoy your reseach, Enjoy your life,...and take care! 44
45 Some curves for κ = [;.5; ; ] function B ; norm(b) = function L ; norm(l) = κ = ( ) ; (.5 ) ; ( ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; κ = ( ) ; (.5 ) ; ( ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 4 function psi ; norm(psi) = 4 function psi ortho ; norm(psi ortho ) = κ = ( ) ; (.5 ) ; ( ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; κ = ( ) ; (.5 ) ; ( ) ; ( ) ; γ = ( ) ; ( ) ; ( ) ; ( ) ; 45
2D penalized spline (continuous-by-continuous interaction)
2D penalized spline (continuous-by-continuous interaction) Two examples (RWC, Section 13.1): Number of scallops caught off Long Island Counts are made at specific coordinates. Incidence of AIDS in Italian
More informationModeling of Price. Ximing Wu Texas A&M University
Modeling of Price Ximing Wu Texas A&M University As revenue is given by price times yield, farmers income risk comes from risk in yield and output price. Their net profit also depends on input price, but
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationA Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options
A Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options Luis Ortiz-Gracia Centre de Recerca Matemàtica (joint work with Cornelis W. Oosterlee, CWI) Models and Numerics
More informationarxiv: v1 [math.st] 18 Sep 2018
Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit
More informationOrthogonality to the value group is the same as generic stability in C-minimal expansions of ACVF
Orthogonality to the value group is the same as generic stability in C-minimal expansions of ACVF Will Johnson February 18, 2014 1 Introduction Let T be some C-minimal expansion of ACVF. Let U be the monster
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 informationsymmys.com 3.2 Projection of the invariants to the investment horizon
122 3 Modeling the market In the swaption world the underlying rate (3.57) has a bounded range and thus it does not display the explosive pattern typical of a stock price. Therefore the swaption prices
More informationMonetary Economics Final Exam
316-466 Monetary Economics Final Exam 1. Flexible-price monetary economics (90 marks). Consider a stochastic flexibleprice money in the utility function model. Time is discrete and denoted t =0, 1,...
More informationPh.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program August 2017
Ph.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program August 2017 The time limit for this exam is four hours. The exam has four sections. Each section includes two questions.
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Risk Measures Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Reference: Chapter 8
More informationThe rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx
1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that
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 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 informationRiemannian Geometry, Key to Homework #1
Riemannian Geometry Key to Homework # Let σu v sin u cos v sin u sin v cos u < u < π < v < π be a parametrization of the unit sphere S {x y z R 3 x + y + z } Fix an angle < θ < π and consider the parallel
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 informationEstimating Macroeconomic Models of Financial Crises: An Endogenous Regime-Switching Approach
Estimating Macroeconomic Models of Financial Crises: An Endogenous Regime-Switching Approach Gianluca Benigno 1 Andrew Foerster 2 Christopher Otrok 3 Alessandro Rebucci 4 1 London School of Economics and
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationFourier Space Time-stepping Method for Option Pricing with Lévy Processes
FST method Extensions Indifference pricing Fourier Space Time-stepping Method for Option Pricing with Lévy Processes Vladimir Surkov University of Toronto Computational Methods in Finance Conference University
More informationEstimating Term Structure of U.S. Treasury Securities: An Interpolation Approach
Estimating Term Structure of U.S. Treasury Securities: An Interpolation Approach Feng Guo J. Huston McCulloch Our Task Empirical TS are unobservable. Without a continuous spectrum of zero-coupon securities;
More information1 Dynamic programming
1 Dynamic programming A country has just discovered a natural resource which yields an income per period R measured in terms of traded goods. The cost of exploitation is negligible. The government wants
More informationCHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION
CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION Szabolcs Sebestyén szabolcs.sebestyen@iscte.pt Master in Finance INVESTMENTS Sebestyén (ISCTE-IUL) Choice Theory Investments 1 / 65 Outline 1 An Introduction
More informationMethod of Characteristics
The Ryan C. Trinity University Partial Differential Equations January 22, 2015 Linear and Quasi-Linear (first order) PDEs A PDE of the form A(x,y) u x +B(x,y) u y +C 1(x,y)u = C 0 (x,y) is called a (first
More informationIntegrating rational functions (Sect. 8.4)
Integrating rational functions (Sect. 8.4) Integrating rational functions, p m(x) q n (x). Polynomial division: p m(x) The method of partial fractions. p (x) (x r )(x r 2 ) p (n )(x). (Repeated roots).
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 informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
More informationDetermining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series
Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series M.B. de Kock a H.C. Eggers a J. Schmiegel b a University of Stellenbosch, South Africa b Aarhus University, Denmark VI
More informationECON 6022B Problem Set 2 Suggested Solutions Fall 2011
ECON 60B Problem Set Suggested Solutions Fall 0 September 7, 0 Optimal Consumption with A Linear Utility Function (Optional) Similar to the example in Lecture 3, the household lives for two periods and
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationStability in geometric & functional inequalities
Stability in geometric & functional inequalities A. Figalli The University of Texas at Austin www.ma.utexas.edu/users/figalli/ Alessio Figalli (UT Austin) Stability in geom. & funct. ineq. Krakow, July
More informationMathematics in Finance
Mathematics in Finance Steven E. Shreve Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 USA shreve@andrew.cmu.edu A Talk in the Series Probability in Science and Industry
More informationDecision theoretic estimation of the ratio of variances in a bivariate normal distribution 1
Decision theoretic estimation of the ratio of variances in a bivariate normal distribution 1 George Iliopoulos Department of Mathematics University of Patras 26500 Rio, Patras, Greece Abstract In this
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
More informationChapter 8. Markowitz Portfolio Theory. 8.1 Expected Returns and Covariance
Chapter 8 Markowitz Portfolio Theory 8.1 Expected Returns and Covariance The main question in portfolio theory is the following: Given an initial capital V (0), and opportunities (buy or sell) in N securities
More informationCorrelation Structures Corresponding to Forward Rates
Chapter 6 Correlation Structures Corresponding to Forward Rates Ilona Kletskin 1, Seung Youn Lee 2, Hua Li 3, Mingfei Li 4, Rongsong Liu 5, Carlos Tolmasky 6, Yujun Wu 7 Report prepared by Seung Youn Lee
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 informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationA Trust Region Algorithm for Heterogeneous Multiobjective Optimization
A Trust Region Algorithm for Heterogeneous Multiobjective Optimization Jana Thomann and Gabriele Eichfelder 8.0.018 Abstract This paper presents a new trust region method for multiobjective heterogeneous
More informationExercises on the New-Keynesian Model
Advanced Macroeconomics II Professor Lorenza Rossi/Jordi Gali T.A. Daniël van Schoot, daniel.vanschoot@upf.edu Exercises on the New-Keynesian Model Schedule: 28th of May (seminar 4): Exercises 1, 2 and
More informationStochastic Volatility (SV) Models
1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to
More informationSYLLABUS AND SAMPLE QUESTIONS FOR MSQE (Program Code: MQEK and MQED) Syllabus for PEA (Mathematics), 2013
SYLLABUS AND SAMPLE QUESTIONS FOR MSQE (Program Code: MQEK and MQED) 2013 Syllabus for PEA (Mathematics), 2013 Algebra: Binomial Theorem, AP, GP, HP, Exponential, Logarithmic Series, Sequence, Permutations
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More informationinduced by the Solvency II project
Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université
More informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationModeling dynamic diurnal patterns in high frequency financial data
Modeling dynamic diurnal patterns in high frequency financial data Ryoko Ito 1 Faculty of Economics, Cambridge University Email: ri239@cam.ac.uk Website: www.itoryoko.com This paper: Cambridge Working
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
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 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 informationMulti-period Portfolio Choice and Bayesian Dynamic Models
Multi-period Portfolio Choice and Bayesian Dynamic Models Petter Kolm and Gordon Ritter Courant Institute, NYU Paper appeared in Risk Magazine, Feb. 25 (2015) issue Working paper version: papers.ssrn.com/sol3/papers.cfm?abstract_id=2472768
More informationMath-Stat-491-Fall2014-Notes-V
Math-Stat-491-Fall2014-Notes-V Hariharan Narayanan December 7, 2014 Martingales 1 Introduction Martingales were originally introduced into probability theory as a model for fair betting games. Essentially
More information3 The Model Existence Theorem
3 The Model Existence Theorem Although we don t have compactness or a useful Completeness Theorem, Henkinstyle arguments can still be used in some contexts to build models. In this section we describe
More informationNon replication of options
Non replication of options Christos Kountzakis, Ioannis A Polyrakis and Foivos Xanthos June 30, 2008 Abstract In this paper we study the scarcity of replication of options in the two period model of financial
More informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
More informationAsset pricing in the frequency domain: theory and empirics
Asset pricing in the frequency domain: theory and empirics Ian Dew-Becker and Stefano Giglio Duke Fuqua and Chicago Booth 11/27/13 Dew-Becker and Giglio (Duke and Chicago) Frequency-domain asset pricing
More informationApplications of Good s Generalized Diversity Index. A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK
Applications of Good s Generalized Diversity Index A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK Internal Report STAT 98/11 September 1998 Applications of Good s Generalized
More informationSYLLABUS AND SAMPLE QUESTIONS FOR MS(QE) Syllabus for ME I (Mathematics), 2012
SYLLABUS AND SAMPLE QUESTIONS FOR MS(QE) 2012 Syllabus for ME I (Mathematics), 2012 Algebra: Binomial Theorem, AP, GP, HP, Exponential, Logarithmic Series, Sequence, Permutations and Combinations, Theory
More informationOn the White Noise of the Price of Stocks related to the Option Prices from the Black-Scholes Equation
IAENG International Journal of Applied Mathematics, 48:, IJAM_48 4 On the White Noise of the Price of Stocks related to the Option Prices from the Black-Scholes Equation A Kananthai, Kraiwiradechachai
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationWe follow Agarwal, Driscoll, and Laibson (2012; henceforth, ADL) to estimate the optimal, (X2)
Online appendix: Optimal refinancing rate We follow Agarwal, Driscoll, and Laibson (2012; henceforth, ADL) to estimate the optimal refinance rate or, equivalently, the optimal refi rate differential. In
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 informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationPh.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program June 2017
Ph.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program June 2017 The time limit for this exam is four hours. The exam has four sections. Each section includes two questions.
More information4: SINGLE-PERIOD MARKET MODELS
4: SINGLE-PERIOD MARKET MODELS Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 4: Single-Period Market Models 1 / 87 General Single-Period
More informationGPD-POT and GEV block maxima
Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,
More informationLecture 7. The consumer s problem(s) Randall Romero Aguilar, PhD I Semestre 2018 Last updated: April 28, 2018
Lecture 7 The consumer s problem(s) Randall Romero Aguilar, PhD I Semestre 2018 Last updated: April 28, 2018 Universidad de Costa Rica EC3201 - Teoría Macroeconómica 2 Table of contents 1. Introducing
More informationAnalytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model
Analytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model Advances in Computational Economics and Finance Univerity of Zürich, Switzerland Matthias Thul 1 Ally Quan
More informationTHE NUMBER OF UNARY CLONES CONTAINING THE PERMUTATIONS ON AN INFINITE SET
THE NUMBER OF UNARY CLONES CONTAINING THE PERMUTATIONS ON AN INFINITE SET MICHAEL PINSKER Abstract. We calculate the number of unary clones (submonoids of the full transformation monoid) containing the
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
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 information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationA Robust Option Pricing Problem
IMA 2003 Workshop, March 12-19, 2003 A Robust Option Pricing Problem Laurent El Ghaoui Department of EECS, UC Berkeley 3 Robust optimization standard form: min x sup u U f 0 (x, u) : u U, f i (x, u) 0,
More informationRho and Delta. Paul Hollingsworth January 29, Introduction 1. 2 Zero coupon bond 1. 3 FX forward 2. 5 Rho (ρ) 4. 7 Time bucketing 6
Rho and Delta Paul Hollingsworth January 29, 2012 Contents 1 Introduction 1 2 Zero coupon bond 1 3 FX forward 2 4 European Call under Black Scholes 3 5 Rho (ρ) 4 6 Relationship between Rho and Delta 5
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More informationProspect Theory, Partial Liquidation and the Disposition Effect
Prospect Theory, Partial Liquidation and the Disposition Effect Vicky Henderson Oxford-Man Institute of Quantitative Finance University of Oxford vicky.henderson@oxford-man.ox.ac.uk 6th Bachelier Congress,
More informationMartingales. by D. Cox December 2, 2009
Martingales by D. Cox December 2, 2009 1 Stochastic Processes. Definition 1.1 Let T be an arbitrary index set. A stochastic process indexed by T is a family of random variables (X t : t T) defined on a
More informationProbability distributions relevant to radiowave propagation modelling
Rec. ITU-R P.57 RECOMMENDATION ITU-R P.57 PROBABILITY DISTRIBUTIONS RELEVANT TO RADIOWAVE PROPAGATION MODELLING (994) Rec. ITU-R P.57 The ITU Radiocommunication Assembly, considering a) that the propagation
More informationInvariant Variational Problems & Integrable Curve Flows. Peter J. Olver University of Minnesota olver
Invariant Variational Problems & Integrable Curve Flows Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Cocoyoc, November, 2005 1 Variational Problems x = (x 1,..., x p ) u = (u 1,...,
More informationPartial Fractions. A rational function is a fraction in which both the numerator and denominator are polynomials. For example, f ( x) = 4, g( x) =
Partial Fractions A rational function is a fraction in which both the numerator and denominator are polynomials. For example, f ( x) = 4, g( x) = 3 x 2 x + 5, and h( x) = x + 26 x 2 are rational functions.
More informationFinancial Modelling Using Discrete Stochastic Calculus
Preprint typeset in JHEP style - HYPER VERSION Financial Modelling Using Discrete Stochastic Calculus Eric A. Forgy, Ph.D. E-mail: eforgy@yahoo.com Abstract: In the present report, a review of discrete
More informationInterpolation of κ-compactness and PCF
Comment.Math.Univ.Carolin. 50,2(2009) 315 320 315 Interpolation of κ-compactness and PCF István Juhász, Zoltán Szentmiklóssy Abstract. We call a topological space κ-compact if every subset of size κ has
More informationStatistics for Business and Economics
Statistics for Business and Economics Chapter 5 Continuous Random Variables and Probability Distributions Ch. 5-1 Probability Distributions Probability Distributions Ch. 4 Discrete Continuous Ch. 5 Probability
More informationEstimating Maximum Smoothness and Maximum. Flatness Forward Rate Curve
Estimating Maximum Smoothness and Maximum Flatness Forward Rate Curve Lim Kian Guan & Qin Xiao 1 January 21, 22 1 Both authors are from the National University of Singapore, Centre for Financial Engineering.
More informationarxiv: v2 [q-fin.cp] 22 Jun 2017
arxiv:1701.01429v2 [q-fin.cp] 22 Jun 2017 Chebyshev Reduced Basis Function applied to Option Valuation. Javier de Frutos and Víctor Gatón June 23, 2017 Abstract We present a numerical method for the frequent
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 informationA Non-Parametric Technique of Option Pricing
1 A Non-Parametric Technique of Option Pricing In our quest for a proper option-pricing model, we have so far relied on making assumptions regarding the dynamics of the underlying asset (more or less realistic)
More informationMarshall and Hicks Understanding the Ordinary and Compensated Demand
Marshall and Hicks Understanding the Ordinary and Compensated Demand K.J. Wainwright March 3, 213 UTILITY MAXIMIZATION AND THE DEMAND FUNCTIONS Consider a consumer with the utility function =, who faces
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationEconomics 101. Lecture 3 - Consumer Demand
Economics 101 Lecture 3 - Consumer Demand 1 Intro First, a note on wealth and endowment. Varian generally uses wealth (m) instead of endowment. Ultimately, these two are equivalent. Given prices p, if
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 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 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 informationNew Business Start-ups and the Business Cycle
New Business Start-ups and the Business Cycle Ali Moghaddasi Kelishomi (Joint with Melvyn Coles, University of Essex) The 22nd Annual Conference on Monetary and Exchange Rate Policies Banking Supervision
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
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 informationBROWNIAN MOTION II. D.Majumdar
BROWNIAN MOTION II D.Majumdar DEFINITION Let (Ω, F, P) be a probability space. For each ω Ω, suppose there is a continuous function W(t) of t 0 that satisfies W(0) = 0 and that depends on ω. Then W(t),
More informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationOptimal Investment with Deferred Capital Gains Taxes
Optimal Investment with Deferred Capital Gains Taxes A Simple Martingale Method Approach Frank Thomas Seifried University of Kaiserslautern March 20, 2009 F. Seifried (Kaiserslautern) Deferred Capital
More information