Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
262 Y. Achdou • ω i are positive weights, • J R is a suitable regularizing functional, • the prices P i (0,S ◦ ) are computed by solving problem (6)–(8), with t = t i and P ◦ (S) =(x i − S) + . Evaluating the functional requires solving #I variational inequalities. This approach was chosen in [Ach05, AP05b] for calibrating models of local volatility (i.e. the volatility is a function of t and S) with American options. In the present case, it is possible to choose a better approach: we call (τ,S) ↦→ P (τ, S, t, x) the pricing function for the vanilla American put with maturity t and strike x. Hereafter, we use the notation P ◦ (x) =(x − S) + . (9) It can be seen that the solution of (6)–(8) is of the form P (τ, S, t, x) =xg(ξ,y), y = S x ∈ R +, ξ = t − τ ∈ (0,τ), where g is the solution of a complementarity problem independent of x, easily deduced from (6)–(8). For brevity, we do not write this problem. From this observation, easy calculations show that, as a function of t and x, P (0,S,t,x) satisfies the following forward problem: P (t =0)=P ◦ and for t ∈ (0,T]andx>0, ( ∂P ∂t − σ2 x 2 ∂ 2 ) P 2 ∂x 2 + rx∂P ∂x + BP ≥ 0, (10) P (t, x) ≥ P ◦ (x), (11) ( ∂P ∂t − σ2 x 2 ∂ 2 ) P 2 ∂x 2 + rx∂P ∂x + BP (P − P ◦ )=0, (12) where the integral operator B is defined by ∫ ( (Bu)(x) =− k(z) x(e z − 1) ∂u ) ∂x (x)+ez (u(xe −z ) − u(x)) dz. (13) z∈R The problem (10)–(12) can also be obtained by probabilistic arguments. The new approach for calibrating the Lévy process is to minimize the functional (σ, k) ↦→ ∑ i∈I ω i| ¯P i − P (t i ,x i )| 2 + J R (σ, k) forσ and k in a suitable set, where the prices P (t i ,x i ) are computed by solving (10)–(12), with P ◦ (x) =(x − S ◦ ) + . In contrast with the previous approach, evaluating the functional requires solving one variational inequality only. Such a forward problem is reminiscent of the forward equation which is often used for the calibration of the local volatility with vanilla European options. This equation is known as Dupire’s equation in the finance community, see [Dup97, AP05a]. Note that the arguments used to obtain (10)–(12) are easier than those used for getting Dupire’s equation, because the operator in (6)–(8) is invariant by any change of variable S ↦→ λS, λ>0, which is not the case with local volatility. Note also that finding a forward linear complementarity problem in the variables t and x is not possible in the case of American options with local volatility.
Calibration of Lévy Processes with American Options 263 Calibration of σ and k is an inverse problem for finding the coefficients of a variational inequality involving a partial integro-differential operator. The main goal of the paper is to study the last least square optimization problem theoretically, for a special parameterization of k, see (25) below, with σ bounded away from 0, and to give necessary optimality conditions. The results presented here have their discrete counterparts when the variational inequalities are discretized with finite elements of finite differences. Numerical results will be presented in a forthcoming paper. 2 Preliminary Results 2.1 Change of Unknown Function in the Forward Problem It is helpful to change the unknown function: we set u ◦ (x) =(S − x) + , u(t, x) =P (t, x) − x + S. (14) The function u satisfies: for t ∈ (0,T]andx>0, ∂u ∂t − σ2 x 2 ∂ 2 u 2 ∂x 2 + rx∂u + Bu ≥−rx, (15) ∂x u(t, x) ≥ u ◦ (x), (16) ( ∂u ∂t − σ2 x 2 ∂ 2 ) u 2 ∂x 2 + rx∂u ∂x + Bu + rx (u − u ◦ )=0. (17) The initial condition for u is u(t =0,x)=u ◦ (x), x > 0. (18) For writing the variational inequalities stemming from (15)–(18), we need to introduce suitable weighted Sobolev spaces. In particular, fractional order weighted Sobolev spaces will be useful for studying the non-local part of the operator. 2.2 Functional Setting Sobolev Spaces on R For a real number s, let the Sobolev space H s (R) be defined as follows: the distribution w defined on R belongs to H s (R) if and only if its Fourier transform ŵ satisfies ∫ (1 + ξ 2 ) s |ŵ(ξ)| 2 dξ < +∞. R
- Page 207 and 208: 206 A. Bonito et al. References [AM
- Page 209 and 210: 208 A. Bonito et al. [Set96] J. A.
- Page 211 and 212: 210 J. Hao et al. due to shear flow
- Page 213 and 214: 212 J. Hao et al. The backward reac
- Page 215 and 216: 214 J. Hao et al. where u and p den
- Page 217 and 218: 216 J. Hao et al. * * * * * * * * *
- Page 219 and 220: 218 J. Hao et al. and solve for V n
- Page 221 and 222: 220 J. Hao et al. Table 2. The calc
- Page 223 and 224: 222 J. Hao et al. References [ASS80
- Page 225 and 226: Computing the Eigenvalues of the La
- Page 227 and 228: Eigenvalues of the Laplace-Beltrami
- Page 229 and 230: Eigenvalues of the Laplace-Beltrami
- Page 231 and 232: Eigenvalues of the Laplace-Beltrami
- Page 233 and 234: A Fixed Domain Approach in Shape Op
- Page 235 and 236: Shape Optimization Problems with Ne
- Page 237 and 238: Shape Optimization Problems with Ne
- Page 239 and 240: Shape Optimization Problems with Ne
- Page 241 and 242: Shape Optimization Problems with Ne
- Page 243 and 244: Reduced-Order Modelling of Dispersi
- Page 245 and 246: Reduced-Order Modelling of Dispersi
- Page 247 and 248: Reduced-Order Modelling of Dispersi
- Page 249 and 250: Reduced-Order Modelling of Dispersi
- Page 251 and 252: Reduced-Order Modelling of Dispersi
- Page 253 and 254: Reduced-Order Modelling of Dispersi
- Page 255 and 256: Calibration of Lévy Processes with
- Page 257: Calibration of Lévy Processes with
- Page 261 and 262: Calibration of Lévy Processes with
- Page 263 and 264: We have proved Calibration of Lévy
- Page 265 and 266: Calibration of Lévy Processes with
- Page 267 and 268: Calibration of Lévy Processes with
- Page 269 and 270: Calibration of Lévy Processes with
- Page 271 and 272: Note that p ∗ satisfies Calibrati
- Page 273 and 274: Calibration of Lévy Processes with
- Page 275 and 276: 280 S. Ikonen and J. Toivanen the p
- Page 277 and 278: 282 S. Ikonen and J. Toivanen Merto
- Page 279 and 280: 284 S. Ikonen and J. Toivanen For H
- Page 281 and 282: 286 S. Ikonen and J. Toivanen and {
- Page 283 and 284: 288 S. Ikonen and J. Toivanen 1.6 1
- Page 285 and 286: 290 S. Ikonen and J. Toivanen 8 Con
- Page 287: 292 S. Ikonen and J. Toivanen [Mer7
Calibration of Lévy Processes with American Options 263<br />
Calibration of σ <strong>and</strong> k is an inverse problem for finding the coefficients<br />
of a variational inequality involving a partial integro-differential operator.<br />
The main goal of the paper is to study the last least square optimization<br />
problem theoretically, for a special parameterization of k, see (25) below,<br />
with σ bounded away from 0, <strong>and</strong> to give necessary optimality conditions. The<br />
results presented here have their discrete counterparts when the variational<br />
inequalities are discretized with finite elements of finite differences. Numerical<br />
results will be presented in a forthcoming paper.<br />
2 Preliminary Results<br />
2.1 Change of Unknown Function in the Forward Problem<br />
It is helpful to change the unknown function: we set<br />
u ◦ (x) =(S − x) + , u(t, x) =P (t, x) − x + S. (14)<br />
The function u satisfies: for t ∈ (0,T]<strong>and</strong>x>0,<br />
∂u<br />
∂t − σ2 x 2 ∂ 2 u<br />
2 ∂x 2 + rx∂u + Bu ≥−rx, (15)<br />
∂x<br />
u(t, x) ≥ u ◦ (x), (16)<br />
( ∂u<br />
∂t − σ2 x 2 ∂ 2 )<br />
u<br />
2 ∂x 2 + rx∂u ∂x + Bu + rx (u − u ◦ )=0. (17)<br />
The initial condition for u is<br />
u(t =0,x)=u ◦ (x), x > 0. (18)<br />
For writing the variational inequalities stemming from (15)–(18), we need<br />
to introduce suitable weighted Sobolev spaces. In particular, fractional order<br />
weighted Sobolev spaces will be useful for studying the non-local part of the<br />
operator.<br />
2.2 Functional Setting<br />
Sobolev Spaces on R<br />
For a real number s, let the Sobolev space H s (R) be defined as follows: the distribution<br />
w defined on R belongs to H s (R) if <strong>and</strong> only if its Fourier transform<br />
ŵ satisfies<br />
∫<br />
(1 + ξ 2 ) s |ŵ(ξ)| 2 dξ < +∞.<br />
R