• 検索結果がありません。

ITÔ-SKOROHOD STOCHASTIC EQUATIONS AND APPLICATIONS TO FINANCE

N/A
N/A
Protected

Academic year: 2022

シェア "ITÔ-SKOROHOD STOCHASTIC EQUATIONS AND APPLICATIONS TO FINANCE"

Copied!
11
0
0

読み込み中.... (全文を見る)

全文

(1)

AND APPLICATIONS TO FINANCE

CIPRIAN A. TUDOR

Received 6 November 2003 and in revised form 8 June 2004

We prove an existence and uniqueness theorem for a class of It ˆo-Skorohod stochastic equations. As an application, we introduce a Black-Scholes market model where the price of the risky asset follows a nonadapted equation.

1. Introduction

The introduction of the anticipating (or Skorohod) integral in [8] and of the anticipat- ing stochastic calculus in [7] has opened the question of solving anticipating stochastic differential equations. In general, the existence and uniqueness of the solution for these equations is not known. The difficulty of solving such equations is due to the fact that the classical method of Picard iterations cannot be applied because the mean square formula for the Skorohod integral involves the Malliavin derivation in a such way that we cannot find “closed” formulas. Only in few particular cases do some results exist; see, for exam- ple, [1,2,3]. We have recently proved in [9] that the set of Skorohod integrals coincides with a set of integrals of It ˆo type. In the present work, using this correspondence between Skorohod integrals and It ˆo-Skorohod integrals, we introduce a class of anticipating equa- tions (calledItˆo-Skorohod equations) that can be solved using standard techniques. As an application, we introduce a market model where the price of the risky asset follows such an equation with a random initial condition (the price at the transaction time). We prove that our model is complete and has no arbitrage opportunities and we derive a Black- Scholes formula when the initial price of the risky asset is given by a standard normal random variable.

We organized the paper as follows.Section 2contains some preliminaries on the an- ticipating stochastic calculus. InSection 3, we define the class of It ˆo-Skorohod equations and we prove the existence and uniqueness of the solution. InSection 4, we introduce a market model with price dynamics following an It ˆo-Skorohod equation and we obtain a Black-Scholes option valuation formula and the expression of the replicant portfolio.

Copyright©2004 Hindawi Publishing Corporation

Journal of Applied Mathematics and Stochastic Analysis 2004:4 (2004) 359–369 2000 Mathematics Subject Classification: 60H05, 60H07

URL:http://dx.doi.org/10.1155/S1048953304311044

(2)

2. Preliminaries

We start with some elements of the Malliavin calculus. We refer to [6] for a complete presentation of this topic. Let (W(t))t[0,1]be a standard Wiener process on the canonical Wiener space (Ω,F,P) and let (Ft)t[0,1]be the filtration generated byW. A functional of the Brownian motion of the form

F=fWt1

,. . .,Wtn, (2.1)

witht1,. . .,tn[0, 1] and f Cb(Rn), is called a smooth random variable and this class is denoted by᏿. The Malliavin derivative is defined on᏿as

DtF= n i=1

∂ f

∂xi

Wt1

,. . .,Wtn

1[0,ti](t), t[0, 1], (2.2)

ifFhas the form (2.1). The operatorDis closable and can be extended to the closure of

᏿with respect to the seminorm

Fk,pp =E|F|p+ k j=1

ED(j)FLp2([0,1]), (2.3)

whereD(i) denotes theith iterated derivative. The adjoint ofDis denoted byδ and is called the Skorohod integral. That is,δis defined on its domain

Dom(δ)=

uL2[0, 1]×/E1

0usDsF dsCFL2(Ω)

(2.4) and is given by the duality relationship

EFδ(u)=E1

0usDsF ds, uDom(δ),F. (2.5) Recall that the variance of the Skorohod integral is

Eδ2(u)=E1

0u2α+E1

0

1

0DβuαDαuβdα dβ. (2.6) ByLk,p, we denote the setL2([0, 1];Dk,p), fork1 and p2, and we note thatLk,pis a subset of the domain ofδ. The following version of the Ocone-Clark formula was given in [7]:

F=EF/F[s,t]c

+ t

sEDαF/F[α,t]c

dW(α), forFD1,2. (2.7)

(3)

We will need the integration-by-parts formula Fδ(u)=δ(Fu) +

[0,1]DsFusds (2.8)

if all above terms are defined. Recall also that ifFis a random variable, Malliavin differ- entiable, and measurable with respect to aσ-algebraFA,AᏮ(R), then

DF=0, onAc×. (2.9)

We define, fork1 andp2, the sets of processes ᏹk,p=

X=

Xtt[0,1],Xt= t

0usdWs,uLk,p, ᏺk,p=

Y=

Yt

t[0,1],Yt= t

0E vs/F[s,t]c

dWs,vLk,p.

(2.10)

We will refer to the elements ofᏺk,p asItˆo-Skorohod integral processes and to the ele- ments ofᏹk,pasSkorohod integral processes.It has been proved in [9] that for sufficiently regular integrands, the two classes coincide. As a consequence, to study Skorohod inte- gral processes, it suffices to study It ˆo-Skorohod integral processes, which have two inter- esting properties. Firstly, note that the integralYt=t

0E[uα/F[α,t]c]dWα exists even for uL2([0, 1]×Ω) and has similarities to a classical It ˆo integral. Observe, by (2.6), that this integral is an “isometry”:

Et

0E uα/F[α,t]c

dWα 2

=Et

0

E uα/F[α,t]c

2

dα. (2.11)

Secondly, if we define, for everyλt,Ytλ=λ

0E[uα/F[α,t]c]dWα, then the process (Ytλ)λt

is anF(λ,t]c-martingale and we have

λlimt,λtYtλ=Yt a.s. and inL2. (2.12) We will now define the stochastic integral with respect to Itˆo-Skorohod integral processes.

Definition 2.1. Letu,vL2([0, 1]×Ω) be adapted processes and consider moreYt=Y0+ t

0E[uα/F[α,t]c]dWα+0tE[vα/F[α,t]c]dα. By definition, for any adapted square integrable processX,

t

0XsdYs:= t

0XsdsYts, (2.13)

where

Ytλ=Y0+ λ

0E uα/F[α,t]c

dWα+ λ

0E vα/F[α,t]c

(2.14)

and the integral on the right-hand side of (2.13) is understood in the semimartingale sense.

(4)

3. Itˆo-Skorohod stochastic equations

In this section, we state and prove an existence and uniqueness theorem for a class of anticipating stochastic differential equations using the method of Picard iterations. It is known that in the anticipating stochastic calculus, this method cannot be applied because the formula of the mean square of the Skorohod integral involves the Malliavin deriva- tive and one cannot find “closed” formulas. We define here a new class of anticipating equations, located “between” It ˆo and Skorohod equations, that can be solved by classical techniques. Consider the following stochastic differential equation:

Xt=Z+ t

0σs,EXs/F[s,t]c

dWs+ t

0bs,Xsds. (3.1) Note that the stochastic integral above is a Skorohod integral since the integrand is not adapted and the initial condition is anticipating. The solution will also be anticipating.

In what follows, the coefficientsσ(t,x),b(t,x) : [0, 1]×RRare given and satisfy the following standard conditions.

(H1) (Measurability):σandbare jointly measurable in (t,x)[0, 1]×R.

(H2) (Lipschitz condition): there exists aD >0 such that for allt[0, 1] andxR, σ(t,x)σ(t,y)+b(t,x)b(t,y)D|xy|. (3.2) (H3) (Linear growth condition): there exists a C >0 such that for all t[0, 1] and

xR,

σ(t,x)2+b(t,x)2C21 +|x|2

. (3.3)

We also make a hypothesis concerning the initial valueZ.

(H4)Zis a random variable withE|Z|2<.

A square integrable process that satisfies a.s. (3.1) is called a solution of (3.1). For given coefficientsσandb, any solutionXwill depend on the initial valueZ. We will say that the solution is unique if, for everyt[0, 1],P(Xt1=Xt2)=1 for any two solutionsX1andX2 with the same initial condition.

We start by proving the existence and the uniqueness of the solution of (3.1).

Theorem3.1. Under assumptions (H1), (H2), (H3), and (H4), stochastic equation (3.1) has a unique solutionXton[0, 1]with

sup

0t1

EXt2<. (3.4)

Proof. Throughout this proof,Kwill denote a generic constant depending only onDand E|Z|2. We consider the usual Picard iterationsXt(0)=Zand

Xt(n+1)=Z+ t

0σs,EXs(n)/F[s,t]c

dWs+ t

0bs,Xs(n)ds. (3.5)

(5)

We first prove the existence of the solution. We have, from (2.6), (H3), and H¨older’s inequalities, that

EXt(1)Xt(0)22Et

0σs,EZ/F[s,t]c

dWs

2+ 2Et

0b(s,Z)ds

2

2Et

0

σs,EZ/F[s,t]c2ds+ 2tEt

0

b(s,Z)2ds

2C2Et

0

1 +EZ/F[s,t]c2

ds+ 2tC2Et

0

1 +|Z|2 ds

Kt.

(3.6)

Using the same arguments and condition (H4), we obtain

EXt(n+1)Xt(n)22E t

0

σs,EXs(n)/F[s,t]c

σs,EXs(n1)/F[s,t]c

dWs

2 + 2Et

0

bs,Xs(n)bs,Xs(n1)ds

2

2Et

0

σs,EXs(n)/F[s,t]c

σs,EXs(n1)/F[s,t]c2

ds + 2tEt

0

bs,Xs(n)bs,Xs(n1)2ds

2D2(1 +t) t

0EXs(n)Xs(n1)2ds.

(3.7)

By induction, one can show that there existsK >0 such that for allt[0, 1] andn1, EXt(n+1)Xt(n)2(Kt)n+1

(n+ 1)!. (3.8)

Relation (3.8) and standard arguments imply the convergence, inL2(Ω), of the successive approximationsXt(n)to a limitXtdefined byXt=Z+n=0(Xt(n+1)Xt(n)).

To prove thatXis a solution, we take theL2(Ω)-limit in (3.5) asn→ ∞. Obviously,

2Et

0

σs,EXs(n)/F[s,t]c

σs,EXs/F[s,t]c

dWs

2

K t

0EXs(n)Xs2ds−→n→∞0, Et

0

bs,Xs(n)bs,Xsds

2

K t

0EXs(n)Xs2ds−→n→∞0.

(3.9)

(6)

The uniqueness of the solution is given by Gronwall’s lemma since for any two solutions X,Y with the same initial condition and for everyt[0, 1], we have

EXtYt2K t

0EXsYs2ds. (3.10) Concerning bound (3.4), we will only note that standard techniques apply (see, e.g., [4]).

Remark 3.2. We define the following stochastic differential equation:

Xt=EZ/Ftc

+ t

0σs,EXs/F[s,t]c

dWs+ t

0bs,EXs/F[s,t]c

ds. (3.11)

Following the lines of the proof ofTheorem 3.1, one can show that (3.11) admits a unique solutionXwith sup0t1E|Xt|2<.

In the particular case of linear coefficients, one can explicitly obtain the solution of (3.11).

Corollary3.3. Letσ,bRandX0L2(Ω). Consider the equation Xt=EX0/Ftc

+ t

0σEXs/F[s,t]c

dWs+ t

0bEXs/F[s,t]c

ds. (3.12)

Then the unique solution of (3.12) is given by Xt=EX0/Ftc

eσWt+(bσ2/2)t. (3.13) Proof. DenoteMt=eσWt+(bσ2/2)t. ThenMtsatisfies the equation

Mt=1 + t

0σMsdWs+ t

0bMsds, (3.14)

and using (2.8) and (2.9), we obtain Xt=EX0/Ftc

Mt

=EX0/Ftc +

t

0σEX0/Ftc

MsdWs+ t

0bEX0/Ftc Msds

=EX0/F[s,t]c

+ t

0σEEX0/Fsc

Ms/F[s,t]c

dWs+ t

0bEEX0/Fsc

Ms/F[s,t]c

ds

=EX0/Ftc +

t

0σEXs/F[s,t]c dWs+

t

0bEXs/F[s,t]c ds.

(3.15)

4. Black-Scholes model driven by Itˆo-Skorohod stochastic differential equations We introduce, in this section, a market model with price dynamics following an It ˆo- Skorohod stochastic equation. As usual, we will consider two assets on the probability

(7)

space (Ω,F,P, (Ft)t[0,1]): the safe investmentA=(At)t[0,1]satisfyingAt=1 +r0tAsds and the risky assetS=(St)t[0,1]with price dynamics following the stochastic differential equation

St=ES0/Ftc +

t

0σESs/F[s,t]c dWs+

t

0bESs/F[s,t]c

ds. (4.1)

Clearly,At=ertandCorollary 3.3implies that St=ES0/Ftc

eσWt+(bσ2/2)t. (4.2) The value of the portfolio at the instanttis defined by

Vt=htAt+HtSt, (4.3)

where the componentsh,HL2([0, 1]×Ω) are adapted to the Brownian filtration and represent the quantities of the safe asset and the risky asset at the instantt.

We say that the portfolio (ht,Ht)t[0,1]isself-financingif Vt=EV0/Ftc

+ t

0hsdAs+ t

0HsdSs, (4.4)

where the differentialdSis understood in the sense ofDefinition 2.1.

Remark 4.1. Note thatDefinition 2.1can be used although the initial value depends ont because, by the Ocone-Clark formula (2.7), we can write

ES0/Ftc

=S0 t

0EDsS0/F[s,t]c

dWs. (4.5)

In other words, the self-financing condition (4.4) can be written as Vt=EV0/Ftc

+ t

0hsrersds t

0HsEDsS0/F[s,t]c dWs

+ t

0HsσESs/F[s,t]c

dWs+ t

0bHsESs/F[s,t]c

ds.

(4.6)

In the following, we will denote by ˜St=ertStthe discounted risky asset price. A nec- essary and sufficient condition for the portfolio, to be self-financing, is given in the next result.

Proposition4.2. Assume thath,HL2([0, 1]×Ω)and let the processVbe given by (4.3).

DenoteV˜t=ertVt.Then the portfolio is self-financing if and only if V˜t=EV0/Ftc

+ t

0HsdS˜s for everyt[0, 1]. (4.7)

(8)

Proof. Suppose thatV satisfies (4.4). Define, for everyλ[0,t], Vλ,t=EV0/Ftc

+ λ

0hsrersds λ

0HsEDsS0/F[s,t]c

dWs

+ λ

0HsσESs/F[s,t]c

dWs+ λ

0bHsESs/F[s,t]c

ds.

(4.8)

It is not difficult to check thatVλ,t=E(Vλ/F[λ,t]c).

We can write It ˆo’s formula foreVλ,t since, for fixedt, the process (Vλ,t)λ[0,t]is a F[λ,t]c-semimartingale. It holds, taking the limit (a.s. or inL2) asλt, that

V˜t=EV0/Ftc

+ t

0ershsrersds t

0ersHsEDsS0/F[s,t]c

dWs

+ t

0ersσHsESs/F[s,t]c

dWs+ t

0ersbHsESs/F[s,t]c

ds+ t

0Vs,trersds

=EV0/Ftc

t

0ersHsEDsS0/F[s,t]c

dWs

+ t

0ersσHsESs/F[s,t]c

dWs+ t

0ers(br)HsESs/F[s,t]c

ds.

(4.9)

On the other hand, writing It ˆo’s formula foreSλ,twith Sλ,t=ES0/Ftc

+ λ

0σESs/F[s,t]c dWs+

λ

0bESs/F[s,t]c

ds=ESλ/F[λ,t]c

, (4.10) we get

S˜t=ES0/Ftc

+ t

0ersσESs/F[s,t]c

dWs+ t

0ers(br)ESs/F[s,t]c

ds. (4.11)

Identity (4.7) follows from (4.9) and the above equation usingDefinition 2.1. The proof

of the necessary part is not more difficult.

LetTbe the exercise time. In the classical Black-Scholes settings, to prove the nonexis- tence of arbitrage, it suffices to exhibit a probability measure equivalent toPunder which the discounted price ˜Sis a martingale. In our case, we have the following.

Proposition4.3. The unique probability measureP˜equivalent toPunder which the pro- cessS˜t/E(S0/Ftc)is a martingale is given by the Radon-Nikodym derivative

dP˜

dP=exprµ σ WT1

2

(rµ)2

σ2 T P-a.s. (4.12)

Under probabilityP, the process˜ W˜t=Wt+ ((br)/σ)tis a standard Brownian motion and the discounted priceS˜satisfies the equation

S˜t=ES˜0/Ftc

+ t

0σES˜s/F[s,t]c

dW˜s. (4.13)

(9)

Proof. DenoteZt=S˜t/E(S0/Ftc)=erteσWt+(bσ2/2)t. It is well known that there exist the unique probability ˜P and the Brownian motion ˜W as above and it holds thatZt=1 + t

0σZsdW˜s. Taking into account that the natural filtrations ofWand ˜Wcoincide, we get S˜t=ES˜0/Ftc

Zt=ES˜0/Ftc

+ t

0σES˜0/Ftc

ZsdW˜s=ES˜0/Ftc

+ t

0σES˜s/F[s,t]c

dW˜s. (4.14)

Remark 4.4. Note that, byCorollary 3.3, we have ˜St=E( ˜S0/Ftc)eσW˜t2/2)t. Also, an im- mediate consequence of Propositions4.2and4.3is the fact that the market is complete and has no arbitrage opportunities.

ConsiderVT=(STK)+the payofffunction of the European call option with exercise timeT and strike priceK. Denote by ˜Ethe expectation with respect to ˜Pand by ˜Dthe Malliavin derivative with respect to ˜W. By formulas (4.7) and (4.13), we have

V˜t=E˜V0/Ftc

t

0HsE˜D˜sS0/F[s,t]c

dW˜s+ t

0σHsE˜S˜s/F[s,t]c

dW˜s. (4.15)

Taking the conditional expectation with respect to theσ-algebraFt, we obtain E˜V˜t/Ft

=E˜V0

t

0HsE˜D˜sS0/Fs

dW˜s+ t

0σHsE˜S˜s/Fs

dW˜s. (4.16)

Therefore, the process ( ˜E( ˜Vt/Ft))t[0,1]is a martingale and, for everytT, it holds that E˜( ˜Vt/Ft)=E˜( ˜VT/Ft) or

Vt/Ft

=E˜er(Tt)VT/Ft

. (4.17)

We have the following option valuation Black-Scholes formula.

Proposition 4.5. Assume that the terminal value is given by VT = f(ST)with f(x)= (xK)+and the initial price of the risky asset isS0=W1+c, wherecis a positive constant.

Then

Vt/Ft

=G

t, St

W1Wt+c

, (4.18)

where

G(t,x)= 1 2π(1T)

×

x

Reu2/(1T)(u+c)Nd1

du

Ker(Tt)

Reu2/(1T)Nd2

du

(4.19)

(10)

with

d1=d1(x,u)=lnK/x(u+c)+r+σ2/2(Tt)

σTt ,

d2=d2(x,u)=d1σTt,

(4.20)

andN(d)=(1/2π)−∞d ex2/2dx.

Proof. Using the fact that the increments of the Wiener process are independent on dis- joint intervals, the Markov property and (4.17) imply that

Vt/Ft

=E˜ er(Tt)fST

/Ft

=E˜ er(Tt)feσWt+(rσ2/2)tE˜S0/FTc

eσ(WTWt)+(rσ2/2)(Tt)

=G

t, St

S0/Ftc

,

(4.21)

where

G(t,x)=er(Tt)E˜ fxES0/FTc

eσ(WTWt)+(rσ2/2)(Tt). (4.22) SinceE(S0/FTc)=W1WT+c, by the joint normal distribution of (W1WT,WTWt),

G(t,x)= 1 2π(1T)

×

Reu2/2(1T)

er(Tt) 2π(Tt)

Rfx(u+c)eσv+(rσ2/2)(Tt)ev2/2(Tt)dv

du.

(4.23) We refer to classical arguments (see [5]) to get

er(Tt) 2π(Tt)

Rfx(u+c)eσv+(rσ2/2)(Tt)ev2/2(Tt)dv

=x(u+c)Nd1(x,u)Ker(Tt)Nd2(x,u),

(4.24)

whered1,d2are given by (4.20) and the conclusion follows.

Since the market is complete, every bounded contingent claim is attainable. There- fore, it is of importance to find the expression of the replicant portfolio. This is given in Proposition 4.6.

Proposition 4.6. Under the hypothesis of Proposition 4.5 and denoting g(t,x)= ertG(t,ertx), the replicant portfolio is given by

Ht=

S˜t

S0/Ftc1 1

σ∂g

∂x

t, S˜t

S0/Ftc

, ht=G

t, St W1Wt+c

cHt

S˜tS0/Ftc

.

(4.25)

(11)

Proof. DenoteMt=S˜t/E˜(S0/Ftc). We utilize the classical procedure to determine the un- known quantitieshandH. We have that

V˜t/Ft

=ertE˜Vt/Ft

=ertG

t, St

S0/Ftc

=ertG

t,ert S˜t

S0/Ftc

=ertGt,ertMt

(4.26)

withGaCfunction on [0,T)×R. Writing It ˆo’s formula forg(t,Mt), we obtain gt,Mt

=g0,M0

+ t

0σ∂g

∂x u,Mu

MudW˜u

+ t

0

∂g

∂t

u,Mudu+1 2

t

0

2g

∂x2

u,Muσ2Mu2du.

(4.27)

Note first that the bounded variation part is zero. On the other hand, by (4.7), E˜V˜t/Ft

=E˜V0

t

0HsED˜sS0/Fs +

t

0σHsE˜S˜s/Fs

dW˜s. (4.28) By (4.27) and (4.28), the natural candidate for H satisfies σ(∂g/∂x)(s,Ms)Ms= σHsE˜( ˜Ss/Fs) and since ˜E( ˜Ss/Fs)=E˜(S0)Ms, we obtain relation (4.25).

References

[1] R. Buckdahn,Linear Skorohod stochastic differential equations, Probab. Theory Related Fields 90(1991), no. 2, 223–240.

[2] , Skorohod stochastic differential equations of diffusion type, Probab. Theory Related Fields93(1992), no. 3, 297–323.

[3] R. Buckdahn and D. Nualart,Linear stochastic differential equations and Wick products, Probab.

Theory Related Fields99(1994), no. 4, 501–526.

[4] P. E. Kloeden and E. Platen,Numerical Solution of Stochastic Differential Equations, Applica- tions of Mathematics, vol. 23, Springer-Verlag, Berlin, 1992.

[5] D. Lamberton and B. Lapeyre,Introduction au Calcul Stochastique Appliqu´e `a la Finance[Intro- duction to Stochastic Calculus Applied to Finance], Edition Ellipses, Paris, 1991.

[6] D. Nualart,The Malliavin Calculus and Related Topics, Probability and Its Applications (New York), Springer-Verlag, New York, 1995.

[7] D. Nualart and ´E. Pardoux,Stochastic calculus with anticipating integrands, Probab. Theory Related Fields78(1988), no. 4, 535–581.

[8] A. V. Skorohod,On a generalization of a stochastic integral, Teor. Verojatnost. i Primenen20 (1975), 219–233.

[9] C. A. Tudor,Martingale-type stochastic calculus for anticipating integral processes, Bernoulli10 (2004), no. 2, 313–325.

Ciprian A. Tudor: Laboratoire de Probabilit´es et Mod`eles Al´eatoires, Universit´e de Paris 6, 4 Place Jussieu, 75252 Paris Cedex 5, France

E-mail address:[email protected]

参照

関連したドキュメント

In this paper we study optimal control problems governed by fractional stochastic partial neutral functional integro-differential equations with infinite delay in Hilbert spaces..

Construction of Lyapunov functionals is usually used for investigation of stability of hereditary systems which are described by functional differential equations or Volterra

In this state space model, the stochastic system model is represented by the stochastic Equations (4) and (5) and the probability distributions given in Section (2.3); the

have considered the exponential stability of neutral stochastic delay partial differential equations by the Lyapunov functional approach; in 8, Dauer and Mahmudov have analyzed

This article demonstrates a systematic derivation of stochastic Taylor methods for solving stochastic delay differential equations (SDDEs) with a constant time lag, r &gt; 0..

A gener- alization of this idea to stochastic differential equations (SDE) and parabolic equations (PE) allows us to prove similar necessary and sufficient conditions for global in

[7] Kunita, H., Stochastic differential equation with jumps and stochastic flows of diffeomorphisms, Itˆ o’s stochastic calculus and probability theory,

As we see in Chapter 2, Malliavin calculus is well known as a method to know regularity properties of distributions of so- lutions of stochastic differential equations driven by