Zagadnienia

3. Struktura terminowa stóp procentowych

Liczba godzin 2.
Zakres materiału:
Wartość pieniądza w czasie. Struktura terminowa stóp procentowych. Stopy spotowe i terminowe. Kontrakty FRA i opcje.

3.1. Wprowadzenie

Rozważmy inwestycję o dwu przepływach gotówki C0<0 w momencie T0C1>0 w T1, T1=T0+t (rys. 1.2).

Rys. 3.1. Przepływy gotówki.

Efektywność takiej inwestycji daje się opisać na kilka sposobów. Na przykład możemy pytać o następujące wielkości:
1. Ile kosztuje w chwili T0 1 jednostka monetarna (np. 1 zł), płatna w chwili T1?

B=C0C1

B nazywamy czynnikiem dyskontującym. Zauważmy, że B=1-D, gdzie D stopa dyskonta.
2. Jaka jest (efektywna) stopa zwrotu?

R=C1C01t-1.

3. Jaka jest (efektywna) intensywność oprocentowania?

R¯=1tlnC1C0.

Powyższe wielkości są ze sobą powiązane:

B=e-R¯t=11+Rt;
R¯=-lnBt;
R=eR¯-1=B-1/t-1.

Struktura terminowa to relacja między efektywnością opisanej powyżej inwestycji, a czasem pozostałym do momentu zapadalności. Rt i R¯t,

R,R¯:0,0,,

opisują strukturę terminową stóp zwrotu (stóp procentowych). Ich wykresy nazywa się ,,krzywą dochodowości” lub ,,krzywą zwrotów” (yield curve). Natomiast Bt,

B:0,0,.

opisuje strukturę terminową czynnika dyskontującego.

Mamy następujące ograniczenia (warunki brzegowe) na funkcję B:

B0=1,limtBt=0.

Ponadto Bt jest ściśle malejąca.
Rzeczywiście, gdyby dla pewnych t1<t2, Bt1Bt2, to prościej byłoby otrzymać 1 jednostkę monetarną w momencie t1 i przechować ją do momentu t2.

Rynek dostarcza nam informacje tylko o wartościach B w skończonej liczbie punktów. Dlatego, bez ograniczania ogólności rozważań, w modelowaniu przyjmuje się, że B jest funkcją klasy C1 (być może poza skończoną liczbą punktów) i ma w punkcie 0 pochodną prawostronną. Pozwala to zdefiniować chwilową stopę (intensywność) spot δ0 i chwilowe stopy (intensywności) forward δt, t>0.

δt=-dlnBtdt=-BtBt,δ0=-B0+.

Zauważmy, że skoro B jest ściśle malejąca, to δ jest prawie wszędzie dodatnia. Ponadto

Bt=exp-0tδsds,R¯t=1t0tδsds.

3.2. Interpretacja R¯t i δt

R¯t jest średnią z δs dla 0st. Dlatego też krzywe będące wykresami R¯t i δt mają ,,podobny kształt”. Spośród wszystkich możliwych kształtów tych krzywych wyróżnia się cztery typowe ([13] s. 120, [18] §3.4).

Rys. 3.2. Płaska struktura terminowa.

δt – funkcja stała. Płaska struktura terminowa. Rynek w równowadze.

Rys. 3.3. δt malejąca.

δt – funkcja malejąca. Rynek oczekuje spadku stóp procentowych (oczekiwany jest spadek inflacji lub spadek zapotrzebowania na kapitał związany ze spadkiem aktywności gospodarczej).

Rys. 3.4. δt rosnąca.

δt – funkcja rosnąca. Rynek oczekuje wzrostu stóp procentowych (oczekiwany jest wzrost inflacji lub wzrost zapotrzebowania na kapitał wynikający ze wzrostu aktywności gospodarczej).

Rys. 3.5. Garb.

δt – posiada maksimum tzw. garb (hump). Rynek oczekuje spadku stóp procentowych w dalszej przyszłości, ale równocześnie jest duży popyt na papiery krótkoterminowe.

Powyższe rysunki są ilustracją tzw. faktów stylizowanych, tzn. pewnych uproszczeń sytuacji, z którymi spotykamy się w życiu. Strukturę terminową modelujemy za pomocą rodziny funkcji zależnych od parametrów. Dobry model powinien opisywać fakty stylizowane.

3.3. Przykłady struktur terminowych

Dla większej czytelności zapisu wprowadzimy nową funkcję

Yt=-lnBt=0tδsds.

Niektórzy autorzy nazywają ją funkcją dochodowości (yield function).

3.3.1. Struktura terminowa płaska

δt=δ – funkcja stała. Wówczas:

Yt=tδ,R¯t=δ,Bt=e-tδ,Rt=eδ-1.

Zauważmy, że wartości Bt, dla tt1,t2, są ważonymi średnimi geometrycznymi wartości początkowej i końcowej.

Lemat 3.1

Jeśli 0<t1<t<t2 to

Bt=Bt1t2-tt2-t1Bt2t-t1t2-t1.

Dowód.

Bt1t2-tt2-t1Bt2t-t1t2-t1=
=exp(-t1δt2-tt2-t1-t2δt-t1t2-t1)=exp(-tδ)=B(t).

3.3.2. Struktura terminowa kawałkami płaska

δt – funkcja kawałkami stała.

δt=δi gdy tti-1,ti,

gdzie 0=t0<t1<t2<<tn=+. Wówczas dla tti-1,ti mamy:

Yt=Yti-1+t-ti-1δi,Bt=Bt1exp-t-ti-1δi,R¯t=δi+ti-1tR¯ti-1-δi,R¯0=δ1.
Lemat 3.2

Dla i<n i tti-1,ti zachodzi:

δi=lnBti-lnBti-1ti-ti-1,
Bt=Bti-1ti-tti-ti-1Btit-ti-1ti-ti-1.

Strukturę kawałkami płaską stosujemy gdy znamy tylko skończoną liczbę wartości Bt, a w pozostałych punktach interpolujemy B za pomocą średniej geometrycznej ważonej.

3.3.3. Wzór Stoodleya

δt=p+s1+rest,p,r,s>0.

Wówczas

Yt=0tδsds=p+st-ln1+rest1+r.
R¯t=p+s-1tln1+rest1+r.
Bt=11+re-p+st+r1+re-pt.

Jak widać z powyższego, Bt jest średnią ważoną czynników dyskontujących wyznaczonych przez płaskie struktury terminowe. Wzór Stoodleya stosuje się, gdy inwestor oczekuje, że w przyszłości będzie obwiązywać jedna z dwu płaskich struktur terminowych, ale nie wie która.

Własności funkcji δ.
Funkcja δ jest ściśle malejąca.

δ0=p+s1+r,limtδt=p.
Rys. 3.6. Wykres funkcji Stoodleya.

3.3.4. Wzór Nelsona-Siegela

Wzór Nelsona-Siegela ([24, §15.4]) jest często używany do przybliżonego opisu struktury terminowej.

δt=β0+β1exp-tτ+β2tτexp-tτ,
τ>0,β00,β0+β10,β2<0β2<β1β0+β2expβ1β2-10.

Wówczas

Yt=0tδsds=β0s-τβ1exp-sτ-τβ2sτ+1exp-sτ0t=
=β0t-τβ1exp(-tτ)-τβ2(tτ+1)exp(-tτ)+τβ1+τβ2=
=β0t-τβ1+β2exp-tτ-1-β2texp-tτ,
R¯t=Ytt=β0-β1+β2exp-tτ-1tτ-β2exp-tτ.

Zaletą wzoru Nelsona-Siegela jest liniowa zależność R¯ od parametrów β. Co przy ustalonym parametrze skali τ pozwala je wyznaczyć (na podstawie danych empirycznych) za pomocą metody najmniejszych kwadratów. Zauważmy, że:

R¯0=δ0=β0+β1,limtR¯t=limtδt=β0.

Ponadto

δt=-β11τexp-tτ+β21-tτ1τexp-tτ=
=1τexp-tτ-β1+β2-β2tτ.

Czyli pochodna δ zeruje się i zmienia znak, gdy

tτ=β2-β1β2=1-β1β2.

Zatem intensywność δt posiada ekstremum, gdy

β1β2<1.

Dla parametru β2 dodatniego jest to maksimum, a dla ujemnego minimum. Maksymalna i odpowiednio minimalna intensywność wynoszą

δτ1-β1β2=β0+β2expβ1β2-1.
Rys. 3.7. Wykres funkcji Nelsona-Siegela dla β1=β0 i β2=5β0.

3.3.5. Wzór Vasička

W tym przypadku wzór na funkcję dochodowości Yt jest końcowym wnioskiem ze stochastycznego modelu struktury terminowej, znanego jako model Vasička ([19, §5.2], [24, §7.4.1]). Model Vasička zależy od czterech dodatnich parametrów a, b, σ i r. Ponadto musi być spełniony warunek 2abσ2. Najprostszy opis otrzymamy wprowadzając dwie funkcje pomocnicze A1t i A2t takie, że

Yt=A1t+A2tr i A2A=1-aA2,A20=0,A1A=bA2-12σ2A22,A10=0.

Po scałkowaniu otrzymujemy

A2t=1a1-e-at,
A1t=σ2-2ab2a2A2t-t+σ24aA22t.

Funkcja A2 jest ściśle rosnąca, A20=0 i limtA2t=1a.

Chwilową intensywność otrzymujemy, różniczkując Yt.

δt=YAt=A1At+A2Atr=
=(1-aA2)r+(bA2-12σ2A22)=r+(b-ar)A2-12σ2A22.

Zauważmy, że δ0=r i limtδt=2ab-σ22a2. Co jest uzasadnieniem warunku 2abσ2 (funkcja δt powinna być nieujemna dla wszystkich dodatnich t).

Okazuje się, że dobierając odpowiednio parametry, możemy otrzymać funkcję δ monotoniczną lub z jednym maksimum.

Lemat 3.3

Jeżeli b-ar0, to funkcja δ jest ściśle malejąca.
Jeżeli 0<b-ar<σ2a, to funkcja δ początkowo rośnie, a potem maleje.
Jeżeli σ2ab-ar, to funkcja δ jest ściśle rosnąca.

Dowód.
Funkcję δ możemy zapisać w postaci funkcji złożonej

δt=PA2t,Px=r+b-arx-σ22x2.

Trójmian kwadratowy Px przyjmuje maksimum w punkcie xm=b-arσ2. Badamy monotoniczność Px w przedziale 0,1a, czyli domknięciu obrazu funkcji A2t.

Gdy xm należy do przedziału otwartego 0,1a, to trójmian Px ma ,,garb”, najpierw rośnie, a potem maleje. Gdy xm leży na lewo od tego przedziału (xm0), to trójmian P jest ściśle malejący, a gdy na prawo (xm1a), to ściśle rosnący. Ponieważ funkcja A2t jest ściśle rosnąca, to samo dotyczy funkcji δ.

Interpretacja parametrów modelu:
r – chwilowa stopa procentowa,
σ – wielkość zaburzeń losowych,
ba – położenie równowagi, do którego ,,dąży” model (nie osiągane ze względu na zaburzenia),
a – szybkość reakcji modelu na zaburzenia.

3.3.6. Wzór CIR

Tym razem wzór na funkcję dochodowości Yt jest końcowym wnioskiem ze stochastycznego modelu struktury terminowej, znanego jako model Coxa, Ingersolla, Rosa ([19, §5.2], [24, §7.4.2]). Model ten zależy od czterech dodatnich parametrów a, b, σ i r. Podobnie jak poprzednio, najprościej opisać go, wprowadzając dwie funkcje pomocnicze A1t i A2t takie, że

Yt=A1t+A2tr i A2A=1-aA2-12σ2A22,A20=0,A1A=bA2,A10=0.

Po scałkowaniu otrzymujemy

A2t=sinhγtγcoshγt+a2sinhγt,   2γ=a2+2σ2
A1t=-2bσ2lnγexpat2γcoshγt+a2sinhγt.

Funkcja A2 jest ściśle rosnąca, A20=0 i limtA2t=α=2a+2γ.

Chwilową intensywność otrzymujemy, różniczkując Yt.

δt=YAt=A1At+A2Atr=
=bA2+(1-aA2-12σ2A22)r=r+(b-ar)A2-12σ2rA22.

Zauważmy, że δ0=r i limtδt=bα. Okazuje się, że tak jak w modelu Vasička, dobierając odpowiednio parametry, możemy otrzymać funkcję δ monotoniczną lub z jednym maksimum. Stosując podobne rozumowanie jak w poprzednim podrozdziale, pokazujemy, że:

Lemat 3.4

Jeżeli b-ar0, to funkcja δ jest ściśle malejąca.
Jeżeli 0<b-ar<σ2αr, to funkcja δ początkowo rośnie, a potem maleje.
Jeżeli σ2αrb-ar, to funkcja δ jest ściśle rosnąca.

Interpretacja parametrów jest identyczna jak w modelu Vasička.

3.4. Instrumenty pochodne

3.4.1. Kontrakty FRA

FRA (ang. Forward Rate Agreement) to umowa pomiędzy dwoma kontrahentami, którzy ustalają wysokość stopy procentowej mającej obowiązywać w przyszłości dla określonej kwoty wyrażonej w walucie transakcji dla z góry ustalonego okresu.

W kontrakcie FRA są dwie strony transakcji:

1. Długa pozycja (ang. long position –- a FRA buyer) otrzymuje przepływ determinowany przez stawkę referencyjną, w zamian za stawkę stałą kontraktu (co niekiedy, w celu uniknięcia popełnienia błędu, bywa zapisywane jako: RECEIVE floating and PAY fixed),
2. Krótka pozycja (ang. short position –- a FRA seller) otrzymuje przepływ stały (którego wysokość w rzeczywistości jest ustalona w dniu zawarcia transakcji) w zamian za stawkę referencyjną (RECEIVE fixed and PAY floating).

W kontrakcie FRA występuje:

1. Data zawarcia transakcji (ang. transaction date) Tt,
2. Data ustalenia stawki referencyjnej (ang. fixing date) Tf,
3. Data rozpoczęcia okresu odsetkowego obowiązywania stawki referencyjnej (ang. start date) Ts; przeważnie jest to data spot (tzn. drugi kolejny dzien roboczy) od dnia fixingu stawki,
4. Data zakończenia okresu odsetkowego obowiązywania stawki referencyjnej (ang. end date) Te.

Tt<Tf<Ts<Te

Stopy (stawki) rref i rFRA podawane są w skali rocznej.

Rozliczenie kontraktu następuje w drugim kolejnym dniu roboczym po ustaleniu stawki Tm=Tf+spot, który przeważnie pokrywa się z datą rozpoczęcia okresu odsetkowego obowiązywania stawki referencyjnej. Podczas rozliczania kontraktu nie dochodzi do rzeczywistej wymiany odsetek a jedynie do transferu zdyskontowanej różnicy stóp rref-rFRATe-Ts w odpowiednim kierunku. A zatem, gdy Tm=Ts, to wypłaty (tzw. kwoty kompensacyjne) wynoszą

CFbuyer=Krref-rFRATe-Ts1+rrefTe-Ts=-CFseller,

gdzie K kwota kontraktu (umowna kwota kapitału, na który opiewa kontrakt), a czas jest liczony w latach (por. [4, §8.1.2.3]).

Stawki rFRA są ściśle powiązane ze strukturą terminową stóp procentowych w dniu zawarcia transakcji.

Lemat 3.5
rFRATe-Ts=BTs-TtBTe-Tt-1.

Dowód.
Rozważmy dwie inwestycje:
A. Sprzedaż kontraktu FRA na okres <Ts,Te> na kwotę 1 i kupno 1 obligacji płacącej 1 w chwili Ts.
B. Kupno obligacji płacących 1 w chwili Te za kwotę BTs-Tt i ich sprzedaż w chwili Tf.

Rozliczenie.

W chwili zawarcia transakcji mamy:

CFA,Tt=-BTs-Tt=CFB,Tt.

Rozliczenie obu transakcji nastąpi w chwili Ts=Tf+spot.

CFA,Ts=1+rFRA-rrefTe-Ts1+rrefTe-Ts,
CFB,Ts=BTs-TtBTe-Tt1+rrefTe-Ts.

Po przyrównaniu obydwu wypłat otrzymujemy

1+rrefTe-Ts+rFRA-rrefTe-Ts=BTs-TtBTe-Tt.

Co po uproszczeniu daje tezę lematu.

3.4.2. Kontrakty IRS

Swap stopy procentowej, IRS (ang. interest rate swap), to kontrakt wymiany płatności odsetkowych, jeden z podstawowych instrumentów pochodnych, będący przedmiotem obrotu na rynku międzybankowym. IRS jest umową pomiędzy dwiema stronami, na podstawie której strony wypłacają sobie wzajemnie (w określonych odstępach czasu w trakcie trwania kontraktu) odsetki od umownego nominału kontraktu, naliczane według odmiennie zdefiniowanych stóp procentowych. Transakcja IRS może być traktowana jako seria kontraktów FRA, albo jako wymiana odsetek od dwóch obligacji kuponowych.

Obecnie w coraz większym zakresie transakcje IRS zawierane są przez korporacje. Dokonują tego m.in. w celu zapewnienia sobie stałego kosztu finansowania (receive floating rate and pay fixed), albo pozyskania tańszego finansowania w walucie obcej (receive foreing fixed rate vs. pay domestic fixed rate).

W rodzinie IRS można wyróżnić:

1. Prosty (waniliowy) swap stopy procentowej (ang. plain vanilla IRS) – strony wymieniają się przepływami uzależnionymi od stopy stałej i zmiennej (fixed rate vs. floating rate).
2. Basis swap - obie strony płacą odsetki wg różnej stopy zmiennej, np. WIBOR 3-miesięczny w zamian za WIBOR 6-miesięczny (floating rate vs. floating rate).
3. Walutowy swap stopy procentowej (ang. currency IRS, CIRS) – strony wymieniają się płatnościami denominowanymi w różnych walutach. Nie należy go mylić ze swapem walutowym.

W waniliowej transakcji IRS wyróżnia się dwie pozycje, przy czym to kierunek płatności stawki zmiennej określa jaką pozycję zajmuje kontrahent w transakcji IRS:

1. Kontrahent A zajmuje długą pozycję w stopie zmiennej – referencyjnej, gdy otrzymuje przepływ wyznaczony przez stawkę zmienną w zamian za ustaloną stawkę stałą (RECEIVE floating and PAY fixed).
2. Kontrahent B zajmuje krótką pozycję w stopie zmiennej (ang. short position), gdy płaci odsetki określone przez stawkę zmienną, a w zamian otrzymuje płatności determinowane przez stawkę stałą (RECEIVE fixed and PAY floating).
W charakterze stopy zmiennej występuje zazwyczaj stopa ”rynkowa” (LIBOR, EURIBOR, WIBOR lub inna, zależnie od rynku). Wysokość stopy stałej dla standardowych kontraktów jest kwotowana przez banki i zwana stopą swapową (ang. swap rate). Jest ona dobrana w taki sposób, by początkowa wartość kontraktu była zerowa.

Należy zwrócić uwagę, że wysokość stopy zmiennej płaconej w danym okresie odsetkowym standardowo ustalana jest z góry na początku tego okresu (tak jak dla lokat bankowych). Niekiedy spotykane są kontrakty, w których stopa ta ustalana jest z dołu (tzw. ang. LIA swap lub Libor in arrears swap). Należą one jednak do grupy skomplikowanych w wycenie, tzw. egzotycznych instrumentów pochodnych.

3.4.3. Opcje na stopę procentową

Opcja na górny pułap stopy procentowej (interest rate cap)
jest instrumentem służącym do ochrony różnego rodzaju zobowiązań głównie długoterminowych przed wzrostem stopy procentowej. Przedmiotem zabezpieczenia mogą być na przykład długoterminowe kredyty lub różnego rodzaju instrumenty finansowe, których oprocentowanie jest zmienne i ustalane na bazie danej stawki referencyjnej. Opcja cap jest wielookresowym odpowiednikiem opcji call. W wyniku zawarcia transakcji cap nabywca opcji otrzymuje gwarancję od sprzedawcy, którym jest najczęściej bank, że wzrost stopy procentowej ponad poziom uzgodniony w umowie zostanie zrekompensowany przez sprzedającego

Opcja na dolny pułap stopy procentowej (interest rate floor)
jest instrumentem służącym do ochrony różnego rodzaju należności długoterminowych (choć nie tylko) przed spadkiem stopy procentowej. Przedmiotem tego zabezpieczenia mogą być długoterminowe lokaty lub różnego rodzaju instrumenty finansowe, których oprocentowanie jest zmienne i ustalane na bazie danej stawki referencyjnej. Opcja floor jest wielookresowym odpowiednikiem opcji put. W wyniku zawarcia transakcji typu floor nabywca opcji otrzymuje od sprzedawcy, którym jest najczęściej bank, gwarancję, że spadek stopy procentowej poniżej uzgodnionego w umowie pułapu zostanie zrekompensowany przez sprzedającego. Kompensata polega na przekazaniu przez wystawcę floor kwoty stanowiącej różnicę pomiędzy referencyjną stopą procentową, a stopą procentową ustaloną w transakcji opcyjnej zwaną pułapem. Wystawca płaci ją na początku każdego podokresu wówczas, gdy różnica pomiędzy powyższymi kwotami jest ujemna. W każdym innym przypadku nie dochodzi do żadnej płatności.

3.5. Ćwiczenia

Ćwiczenie 3.1

Rozważamy strukturę terminową stóp procentowych opisaną przez chwilową intensywność

δt=11+t.

Inwestor zainwestował 1 JM na okres n lat. Obliczyć, ile wynosi wypłata i efektywna stopa zwrotu R.

Rozwiązanie. Niech K oznacza wypłatę inwestora po n latach. Wiemy, że
1=KBn. Zatem

K=Bn-1=exp0nδtdt=exp0n11+tdt=expln1+t0n=1+n.

Ponadto

R=Kn-1=n+1n-1.

Odpowiedź. Wypłata jest równa n+1, a efektywna stopa zwrotu n+1n-1.

Ćwiczenie 3.2

Udowodnić lemat 3.4.

Ćwiczenie 3.3

Wyznaczyć stopę rFRA dla półrocznego okresu odsetkowego, który rozpocznie się za dwa lata
a. dla płaskiej struktury terminowej.
b. dla struktury terminowej opisanej wzorem Stoodleya.
c. dla struktury terminowej opisanej wzorem Nelsona-Siegela.
d. dla struktury terminowej opisanej wzorem Vasička.
Zakładamy, że Tm=Ts.

Ćwiczenie 3.4

Struktura terminowa jest wyznaczona przez chwilową stopę δt. Niech rFRAh,t oznacza stopę FRA dla okresu odsetkowego o długości h, który rozpocznie się za t lat. Zakładamy, że Tm=Ts, wyznaczyć

limh0+rFRAh,t.

Treść automatycznie generowana z plików źródłowych LaTeXa za pomocą oprogramowania wykorzystującego LaTeXML.

Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego.

Projekt współfinansowany przez Ministerstwo Nauki i Szkolnictwa Wyższego i przez Uniwersytet Warszawski.