QuantLib: a free/open-source library for quantitative finance
fully annotated source code - version 1.34
Loading...
Searching...
No Matches
hestonprocess.hpp
Go to the documentation of this file.
1/* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
2
3/*
4 Copyright (C) 2005, 2007, 2009, 2014 Klaus Spanderen
5
6 This file is part of QuantLib, a free-software/open-source library
7 for financial quantitative analysts and developers - http://quantlib.org/
8
9 QuantLib is free software: you can redistribute it and/or modify it
10 under the terms of the QuantLib license. You should have received a
11 copy of the license along with this program; if not, please email
12 <quantlib-dev@lists.sf.net>. The license is also available online at
13 <http://quantlib.org/license.shtml>.
14
15 This program is distributed in the hope that it will be useful, but WITHOUT
16 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
17 FOR A PARTICULAR PURPOSE. See the license for more details.
18*/
19
20/*! \file hestonprocess.hpp
21 \brief Heston stochastic process
22*/
23
24#ifndef quantlib_heston_process_hpp
25#define quantlib_heston_process_hpp
26
29#include <ql/quote.hpp>
30
31namespace QuantLib {
32
33 //! Square-root stochastic-volatility Heston process
34 /*! This class describes the square root stochastic volatility
35 process governed by
36 \f[
37 \begin{array}{rcl}
38 dS(t, S) &=& \mu S dt + \sqrt{v} S dW_1 \\
39 dv(t, S) &=& \kappa (\theta - v) dt + \sigma \sqrt{v} dW_2 \\
40 dW_1 dW_2 &=& \rho dt
41 \end{array}
42 \f]
43
44 \ingroup processes
45 */
47 public:
57
61 Real v0,
62 Real kappa,
63 Real theta,
64 Real sigma,
65 Real rho,
67
68 Size size() const override;
69 Size factors() const override;
70
71 Array initialValues() const override;
72 Array drift(Time t, const Array& x) const override;
73 Matrix diffusion(Time t, const Array& x) const override;
74 Array apply(const Array& x0, const Array& dx) const override;
75 Array evolve(Time t0, const Array& x0, Time dt, const Array& dw) const override;
76
77 Real v0() const { return v0_; }
78 Real rho() const { return rho_; }
79 Real kappa() const { return kappa_; }
80 Real theta() const { return theta_; }
81 Real sigma() const { return sigma_; }
82
83 const Handle<Quote>& s0() const;
86
87 Time time(const Date&) const override;
88
89 // probability densitiy function,
90 // semi-analytical solution of the Fokker-Planck equation in x=ln(s)
91 Real pdf(Real x, Real v, Time t, Real eps=1e-3) const;
92
93 private:
94 Real varianceDistribution(Real v, Real dw, Time dt) const;
95
100 };
101}
102#endif
1-D array used in linear algebra.
Definition: array.hpp:52
Concrete date class.
Definition: date.hpp:125
Shared handle to an observable.
Definition: handle.hpp:41
Square-root stochastic-volatility Heston process.
Array drift(Time t, const Array &x) const override
returns the drift part of the equation, i.e.,
Size size() const override
returns the number of dimensions of the stochastic process
Array evolve(Time t0, const Array &x0, Time dt, const Array &dw) const override
Real pdf(Real x, Real v, Time t, Real eps=1e-3) const
Real varianceDistribution(Real v, Real dw, Time dt) const
Matrix diffusion(Time t, const Array &x) const override
returns the diffusion part of the equation, i.e.
const Handle< YieldTermStructure > & dividendYield() const
Time time(const Date &) const override
Size factors() const override
returns the number of independent factors of the process
Discretization discretization_
Array initialValues() const override
returns the initial values of the state variables
Handle< YieldTermStructure > dividendYield_
Array apply(const Array &x0, const Array &dx) const override
const Handle< Quote > & s0() const
const Handle< YieldTermStructure > & riskFreeRate() const
Handle< YieldTermStructure > riskFreeRate_
Matrix used in linear algebra.
Definition: matrix.hpp:41
multi-dimensional stochastic process class.
const DefaultType & t
Date d
Real Time
continuous quantity with 1-year units
Definition: types.hpp:62
QL_REAL Real
real number
Definition: types.hpp:50
std::size_t Size
size of a container
Definition: types.hpp:58
Definition: any.hpp:35
ext::shared_ptr< BlackVolTermStructure > v
purely virtual base class for market observables
stochastic processes
Interest-rate term structure.