11 svar
113 visningar
Moffen 1875
Postad: 8 dec 2020 15:20

Stokastisk analys/options prissättning

Hej!

Jag sitter fast på följande uppgift:

Jag känner mig väldigt osäker så fort vi blandar in någon typ av "max av {...}".

Jag vet inte riktigt vart jag ska börja, någonting med risk-neutral valuation?

Tack!

Albiki 5096 – Fd. Medlem
Postad: 8 dec 2020 16:47 Redigerad: 8 dec 2020 17:31

Hej,

Du kan börja med att skriva systemet på matrisform.

    dSt=AStdt+ΣStdWt\displaystyle dS_t = AS_t\,dt + \Sigma S_t \,dW_t

där

    St=St1St2S_t = \begin{pmatrix}S^1_t\\S^2_t\end{pmatrix}

och

    A=α100α2A = \begin{pmatrix}\alpha_1&0\\0&\alpha_2\end{pmatrix}

och

    Σ=σ10ρσ21-ρ2σ2\Sigma=\begin{pmatrix}\sigma_1&0\\\rho\sigma_2&\sqrt{1-\rho^2}\sigma_2\end{pmatrix}

samt

    Wt=Wt1Wt2.W_t=\begin{pmatrix}W^1_t\\W^2_t\end{pmatrix}.

Korrigering: Matristyperna passar inte ihop i multiplikationen ΣStdWt\Sigma S_t dW_t så denna term berhöver rättas till.

Moffen 1875
Postad: 8 dec 2020 17:48 Redigerad: 8 dec 2020 18:03

Just det, tack. Jag tror att jag är någonting på spåret genom att använda "state-space reduction", dvs. skriva om funktionen Φs1,s2=s2Φs1s2,1\Phi\left(s_1,s_2\right)=s_{2}\Phi\left(\frac{s_1}{s_2},1\right), och sen använda z=s1s2z=\frac{s_1}{s_2}, låter det rimligt?

Så vår prissättningsfunktion är Ft,s=s2Gt,zF\left(t,s\right)=s_{2}G\left(t,z\right)

Vi behöver alltså lösa följande PDE:

Gt+12z2(C11-2C12+C22)Gzz=0G(T,z)=ϕ(z,1)

där Φz,1=maxS1tS2t,1\Phi\left(z,1\right)=max\left(\frac{S_{1}\left(t\right)}{S_{2}\left(t\right)},1\right) och C11=σ1,C12=C21=ρσ2,C22=1-ρ2σ2C_{11}=\sigma_{1}, C_{12}=C_{21}=\rho \sigma_{2}, C_{22}=\sqrt{1-\rho^2}\sigma_2.

Har jag tänkt fel någonstans? Hur kommer jag vidare? Känns som ett hav av variabler och matriser att ha koll på...

EDIT: Black-Scholes ekvation kanske? Det passar in i PDE'n ovan.

EDIT 2: Jag menar också att Ft,s=Ft,s1,s2F\left(t,s\right)=F\left(t,s_{1},s_{2}\right).

Albiki 5096 – Fd. Medlem
Postad: 9 dec 2020 00:38
Moffen skrev:

Just det, tack. Jag tror att jag är någonting på spåret genom att använda "state-space reduction", dvs. skriva om funktionen Φs1,s2=s2Φs1s2,1\Phi\left(s_1,s_2\right)=s_{2}\Phi\left(\frac{s_1}{s_2},1\right), och sen använda z=s1s2z=\frac{s_1}{s_2}, låter det rimligt?

Så vår prissättningsfunktion är Ft,s=s2Gt,zF\left(t,s\right)=s_{2}G\left(t,z\right)

Vi behöver alltså lösa följande PDE:

Gt+12z2(C11-2C12+C22)Gzz=0G(T,z)=ϕ(z,1)

där Φz,1=maxS1tS2t,1\Phi\left(z,1\right)=max\left(\frac{S_{1}\left(t\right)}{S_{2}\left(t\right)},1\right) och C11=σ1,C12=C21=ρσ2,C22=1-ρ2σ2C_{11}=\sigma_{1}, C_{12}=C_{21}=\rho \sigma_{2}, C_{22}=\sqrt{1-\rho^2}\sigma_2.

Har jag tänkt fel någonstans? Hur kommer jag vidare? Känns som ett hav av variabler och matriser att ha koll på...

EDIT: Black-Scholes ekvation kanske? Det passar in i PDE'n ovan.

EDIT 2: Jag menar också att Ft,s=Ft,s1,s2F\left(t,s\right)=F\left(t,s_{1},s_{2}\right).

Ingenting av det du skrivit verkar vettigt för mig.

  • Vad är Φ\Phi och varför skulle den kunna reduceras?
  • Varför är prissättningen så som du skriver
  • och varför styrs den av den PDE som du angivit?
  • Är Φ\Phi och ϕ\phi samma funktion? Är ϕ\phi differentierbar?
Moffen 1875
Postad: 11 dec 2020 14:19 Redigerad: 11 dec 2020 14:23

Hej!
Jag ber om ursäkt, härledningen ska finnas i boken "Arbitrage theory in continuous time" av Tomas Björk, vet tyvärr inte riktigt vilken utgåva.

Ja, phi och phi ska vara samma funktion, fick inte till falluppdelnings klammern i LaTeX här på PA. Jag tror inte att dom behöver vara differentierbara.

Hur som helst, om vi tar det som givet så är jag fast vid (om jag nu gjort rätt, inte helt säker men jag bifogar min lösning iallafall) att beräkna 𝔼maxZt,1\mathbb{E}\left[\text{max}\left\{Z\left(t\right),1\right\}\right] där Zt=S1(t)S2(t)Z\left(t\right)=\dfrac{S_{1}(t)}{S_{2}(t)}.

Visa spoiler

Ser du något uppenbart fel?

Jag har ingen aning hur jag kommer vidare.

EDIT: Jag inser nu också att det var ett dåligt val att använda σ\sigma som matris, element i matris, och därmed "element" i matris.

Albiki 5096 – Fd. Medlem
Postad: 11 dec 2020 17:31 Redigerad: 11 dec 2020 17:34

Hej,

Det första hindret som måste överbryggas är korrelationen mellan bruset som styr S2S^2 och bruset som styr S1S^1; vanligtvis antas dessa brus vara oberoende processer.

Dynamiken skrivs på matrisform

    dSt=diag(St)αdt+Σρdiag(St)dWt\displaystyle dS_t = \text{diag}(S_t)\alpha \,dt + \Sigma_\rho \text{diag}(S_t)\,dW_t

där volatilitetsmatrisen Σρ\Sigma_\rho hanterar korrelationen mellan de två brusprocesserna. Vektor WtW_t består som vanligt av två oberoende Wienerprocesser.

Beteckningen diag(St)\text{diag}(S_t) är diagonalmatris bildad av vektor StS_t.

    diag(St)=St100St2\text{diag}(S_t) = \begin{pmatrix}S^1_t&0\\0&S^2_t\end{pmatrix} och α=α1α2\alpha = \begin{pmatrix}\alpha_1\\\alpha_2\end{pmatrix}

och för varje fixerad korrelation ρ\rho definieras

    Σρ=σ10ρσ2σ21-ρ2\Sigma_\rho=\begin{pmatrix}\sigma_1&0\\\rho \sigma_2&\sigma_2 \sqrt{1-\rho^2}\end{pmatrix}.

Nu är dynamiken uttryckt på standardform och existerande resultat kan tillämpas för att bestämma den sökta prisprocessen för det aktuella kontraktet χ\chi.

Moffen 1875
Postad: 11 dec 2020 17:44
Albiki skrev:

Hej,

Det första hindret som måste överbryggas är korrelationen mellan bruset som styr S2S^2 och bruset som styr S1S^1; vanligtvis antas dessa brus vara oberoende processer.

Dynamiken skrivs på matrisform

    dSt=diag(St)αdt+Σρdiag(St)dWt\displaystyle dS_t = \text{diag}(S_t)\alpha \,dt + \Sigma_\rho \text{diag}(S_t)\,dW_t

där volatilitetsmatrisen Σρ\Sigma_\rho hanterar korrelationen mellan de två brusprocesserna. Vektor WtW_t består som vanligt av två oberoende Wienerprocesser.

Beteckningen diag(St)\text{diag}(S_t) är diagonalmatris bildad av vektor StS_t.

    diag(St)=St100St2\text{diag}(S_t) = \begin{pmatrix}S^1_t&0\\0&S^2_t\end{pmatrix} och α=α1α2\alpha = \begin{pmatrix}\alpha_1\\\alpha_2\end{pmatrix}

och för varje fixerad korrelation ρ\rho definieras

    Σρ=σ10ρσ2σ21-ρ2\Sigma_\rho=\begin{pmatrix}\sigma_1&0\\\rho \sigma_2&\sigma_2 \sqrt{1-\rho^2}\end{pmatrix}.

Nu är dynamiken uttryckt på standardform och existerande resultat kan tillämpas för att bestämma den sökta prisprocessen för det aktuella kontraktet χ\chi.

Hej!

Tack för svar.

Jag tror att detta (tillvägagångssättet i min början till en lösning) är den enda metod vi gått igenom för "multidimensional cases". Men du har nog rätt i att det kanske blir lite knas då S1tS_{1}\left(t\right) och S2tS_{2}\left(t\right) inte är helt okorrolerade.

Har du någon annan idé hur man kan lösa detta? Eller en annan idé genom att använda "state-space reduction" (metoden i min lösning)?

Albiki 5096 – Fd. Medlem
Postad: 11 dec 2020 17:44 Redigerad: 11 dec 2020 17:45

För att du ska kunna använda existerande resultat om kontraktsfunktionen χ\chi behöver du se till att det aktuella kontraktet kan skrivas på formen

    χ:Φ(x,y)=x·Φ(z,1)\chi: \Phi(x,y) = x\cdot \Phi(z,1) där z=yxz=\frac{y}{x}.

Här är Φ(x,y)=min(x,y)\Phi(x,y) = \min(x,y) och du behöver undersöka om man kan skriva

    min(x,y)=x·min(z,1)\min(x,y) = x \cdot \min(z,1) där z=y/xz = y/x.

Albiki 5096 – Fd. Medlem
Postad: 11 dec 2020 17:50
Moffen skrev:
Albiki skrev:
Visa spoiler

Hej,

Det första hindret som måste överbryggas är korrelationen mellan bruset som styr S2S^2 och bruset som styr S1S^1; vanligtvis antas dessa brus vara oberoende processer.

Dynamiken skrivs på matrisform

    dSt=diag(St)αdt+Σρdiag(St)dWt\displaystyle dS_t = \text{diag}(S_t)\alpha \,dt + \Sigma_\rho \text{diag}(S_t)\,dW_t

där volatilitetsmatrisen Σρ\Sigma_\rho hanterar korrelationen mellan de två brusprocesserna. Vektor WtW_t består som vanligt av två oberoende Wienerprocesser.

Beteckningen diag(St)\text{diag}(S_t) är diagonalmatris bildad av vektor StS_t.

    diag(St)=St100St2\text{diag}(S_t) = \begin{pmatrix}S^1_t&0\\0&S^2_t\end{pmatrix} och α=α1α2\alpha = \begin{pmatrix}\alpha_1\\\alpha_2\end{pmatrix}

och för varje fixerad korrelation ρ\rho definieras

    Σρ=σ10ρσ2σ21-ρ2\Sigma_\rho=\begin{pmatrix}\sigma_1&0\\\rho \sigma_2&\sigma_2 \sqrt{1-\rho^2}\end{pmatrix}.

Nu är dynamiken uttryckt på standardform och existerande resultat kan tillämpas för att bestämma den sökta prisprocessen för det aktuella kontraktet χ\chi.

 

Hej!

Tack för svar.

Jag tror att detta (tillvägagångssättet i min början till en lösning) är den enda metod vi gått igenom för "multidimensional cases". Men du har nog rätt i att det kanske blir lite knas då S1tS_{1}\left(t\right) och S2tS_{2}\left(t\right) inte är helt okorrolerade.

Har du någon annan idé hur man kan lösa detta? Eller en annan idé genom att använda "state-space reduction" (metoden i min lösning)?

Jag har givit dig lösning till problemet; allt är upplagt för att du ska tillämpa gängse metodik för oberoende Wienerprocesser på den aktuella situationen. Den enda utmaningen i uppgiften handlade om hur korrelationen skulle hanteras.

Moffen 1875
Postad: 11 dec 2020 18:02 Redigerad: 11 dec 2020 18:31

Jag har givit dig lösning till problemet; allt är upplagt för att du ska tillämpa gängse metodik för oberoende Wienerprocesser på den aktuella situationen. Den enda utmaningen i uppgiften handlade om hur korrelationen skulle hanteras.

Ok. Så min nuvarande lösning som jag la i spoilern är OK (förutsatt att dynamiken jag skrivit upp för ZtZ\left(t\right) är korrekt).

Då återstår det fortfarande att beräkna 𝔼maxZt,1\mathbb{E}\left[\text{max}\left\{Z\left(t\right),1\right\}\right], och det är väl där jag sitter fast just nu.

Tips?

EDIT: Dynamiken för ZtZ\left(t\right) hittade jag genom att jämföra dess PDE med Ito's formel.

Albiki 5096 – Fd. Medlem
Postad: 11 dec 2020 18:55

Var kommer max in i bilden? Pratar vi inte om samma problem längre? Uppgiften handlar om min-kontrakt, inte om max-kontrakt. Bokens exempel handlar om max-kontrakt... 

Moffen 1875
Postad: 11 dec 2020 19:13 Redigerad: 11 dec 2020 19:14
Albiki skrev:

Var kommer max in i bilden? Pratar vi inte om samma problem längre? Uppgiften handlar om min-kontrakt, inte om max-kontrakt. Bokens exempel handlar om max-kontrakt... 

Ursäkta!

Jag menar självklart min-kontrakt...

Svara
Close