10. Wykład X, 4.XII.2009
Pod koniec Wykładu IX poznaliśmy już pewne zastosowanie fundamentalnego
Twierdzenia 9.2 (pochodzącego od Karusha i Kuhna-Tuckera, podanego w
Wykładzie IX w w formacie `wteddy' w wersji z obiektami wklęsłymi/wypukłymi),
co prawda zastosowanie dość odległe od analizy portfelowej.
Nasuwa się przy tej okazji pytanie, czy może czytelnik umie rozwiązać
Ćwiczenie 9.5 (Wykład IX) bez Twierdzenia 9.2? Jeśli nie,
to skromniejsze pytanie mogłoby brzmieć tak: czy w rozwiązaniu Ćwiczenia 9.5
podanym w Wykładzie IX koniecznych jest dziesięć warunków nierównościowych?
Np zmniejszenie ilości węzłów z pięciu do trzech dałoby tylko sześć
warunków nierównościowych, lecz (chyba) byłoby to wylanie dziecka
razem z kąpielą.
Pora teraz na zastosowanie bliższe nurtowi (i nazwie) wykładu.
Przez chwilę będzie znowu o modelach ± doskonale skorelowanych.
Czytelnik nabierze wprawy w posługiwaniu się Twierdzeniem 9.2,
rozwiązując samodzielnie (względnie śledząc rozwiązanie ukryte
pod dodatkowym kliknięciem) następujące proste
Ćwiczenie 10.1
W modelu ± doskonale skorelowanym z Przykładu 3.1 w Wykładzie III,
w aspekcie M [z którym związany jest Rysunek 3.2; w Wykładzie III
jeszcze nie znaliśmy aspektu B] znaleźć wszystkie portfele
minimalnego ryzyka, przy czym specjalnie stosując w tym
celu Twierdzenie 9.2.
Wskazówka:
Chodzi o minimalizację wariancji portfela x∈Δ4, która w tym
modelu ± doskonale skorelowanym wyraża się, jeśli pamiętamy to jeszcze
z Wykładu III, wzorem 4x1+x2-x3-4x42. To minimum jest bardzo
łatwo znaleźć tak po prostu, przy czym osiągane jest ono w punktach przecięcia
płaszczyzny krytycznej z sympleksem Δ4 (porównaj Przykład 6.2).
Teraz chodzi nam o uzyskanie tego samego jeszcze raz, wprost z twierdzenia
K-KT.
Rozwiązanie:
Po odpowiednim wyspecyfikowaniu Twierdzenia 9.2, zresztą bardzo podobnym
do tego, które rozwijane jest niżej w tym wykładzie w kontekście relatywnego
minimum ryzyka, szukamy portfeli x=x1,x2,x3,x4T∈Δ4,
dla których istnieje λ=λx∈R takie, że: po pierwsze
(używając nie całej, tylko połowy wariancji)
|
4x1+x2-x3-4x441-1-4+λ1111≥0 |
| (10.1) |
jako wektor w R4 (= nieujemność każdej składowej wektora
po lewej stronie), oraz, po drugie, iloczyn skalarny
wektora stojącego po lewej stronie (10.1) z wektorem
x jest zero:
|
4x1+x2-x3-4x42+λ=0 |
| (10.2) |
(warunek komplementarności przy przyjętej specyfikacji twierdzenia).
Dodajmy teraz stronami na przykład dwie środkowe nierówności w
(10.1).
Dostajemy nieujemność współczynnika λ. Lecz (10.2)
daje niedodatniość λ.
Zatem λ=0, więc też, znowu z (10.2),
4x1+x2-x3-4x4=0. Ryzyko w aspekcie M
minimalizują więc te portfele Markowitza, które … mają
zerowe ryzyko. (Porównaj też jeszcze raz Przykład 6.2.)
Niby nic, lecz ,,wycisnęliśmy tę” informację z Twierdzenia 9.2!
Już w tym rozwiązaniu wychodzą cechy charakterystyczne twierdzenia K-KT.
Po bliższym przyjrzeniu się okazuje się, że jedna część jego tezy, tu
zapisana jako (10.1), ciągnie niejako szukany współczynnik
Lagrange'a (w ćwiczeniu jest on tylko jeden) w jedną stronę.
Natomiast druga część tezy, tu (10.2), ciągnie ten
współczynnik w drugą stronę.
W efekcie wyłania się jakaś równowagowa wartość takiego współczynnika,
która de facto jest zakodowanym opisem ekstremów. Nie przypadkiem
w literaturze francuskojęzycznej, np w [25], twierdzenie
K-KT jest nazywane twierdzeniem o punkcie siodłowym.
Ćwiczenie 10.2
Znaleźć przy pomocy Twierdzenia 9.2 wszystkie portfele
minimalnego ryzyka w aspekcie M w dowolnym modelu ±
doskonale skorelowanym.
Naszym głównym celem w Wykładzie X jest zastosowanie Twierdzenia
Karusha-Kuhna-Tuckera 9.2 (w przyszłości będziemy już pisać
tylko krótko `K-KT') do zautomatyzowanego poszukiwania portfeli
relatywnie minimalnego ryzyka w modelach Markowitza,
które bliższe są warunkom giełdowym niż modele ± doskonale
skorelowane.
Dla (niezdegenerowanych) modeli Blacka w Wykładzie VI polegało to
na dość standardowym zastosowaniu wiedzy z AM II, bo dwa warunki
ograniczające były tylko równościowe. W efekcie gładko dostaliśmy
(w Twierdzeniu 6.1) wszystkie portfele krytyczne Blacka układające
się na prostej krytycznej Blacka.
Teraz jednak dochodzi k ograniczeń nierównościowych związanych
z leżeniem portfeli w sympleksie Δk (nieujemność składowych
portfeli, brak krótkiej sprzedaży), co komplikuje sytuację. Jest
to jakościowo nowy problem. Czym teraz zastąpiona zostanie prosta
krytyczna Blacka?
Zgodnie z tradycją przyjętą w analizie portfelowej, współczynniki
Lagrange'a związane z warunkami równościowymi eTx=1,
μTx=E (wartości E nie będą teraz dowolne, tylko
ograniczone do powłoki wypukłej wartości oczekiwanych stóp zwrotu
ze spółek numer 1, 2,…,k) piszemy jako, odpowiednio, λ
i -λE, zaś nierówności xi≥0 (i=1, 2,…,k)
kodujemy jako ai+2Tx≤bi+2, gdzie
ai+2=0,…,-1↑i,…, 0T,
bi+2=0, i=1, 2,…,k. Jako funkcję f bierzemy
która dla Σ≥0 jest wypukła i różniczkowalna w każdym
punkcie y, zaś jej gradient to ∇fy=Σy. Twierdzenie K-KT
zastosowane w tej sytuacji mówi, że x spełniający wymienione warunki
jest globalnym minimum f przy podanych warunkach (ograniczeniach)
⟺ istnieją λ,λE∈R takie,
że
|
Σx+λe-λEμ≥0 jako wektor w Rkoraz¯xTΣx+λe-λEμ=0warunek komplementarności. |
|
Podamy teraz warunek dostateczny, przy którym takie
portfele relatywnie minimalnego ryzyka istnieją i są jednoznacznie
wyznaczone dla każdej wartości E∈min1≤i≤kμi,max1≤i≤kμi.
Jest on trochę za silny (mógłby być trochę osłabiony), co w praktyce
jednak nie przeszkadza, bo przy estymacji wartości oczekiwanych z
historycznych danych giełdowych praktycznie nigdy nie dostaniemy
pary równych wartości. W zamian zaś pozwala zręcznie opisać
procedurę poszukiwania.
Mianowicie, w dalszym ciągu zakładamy, że Σ>0 oraz
tzn. wszystkie wartości oczekiwane μi są różne.
Te założenia będą obowiązywać do odwołania. (Zauważamy, że
są to silniejsze założenia niż dawna informacja (5.2),
która wystarczała w teorii Blacka. Przy różnych pod-zestawach
zmiennych będziemy w obecnym algorytmie wracać do (5.2).)
Algorytm prowadzący do rozwiązania nie mówi od razu, dla jakich
konkretnie wartości E rozwiązanie znajdziemy na konkretnej ścianie
sympleksu. Całość poszukiwań dzielimy (niestety lub stety) na etapy,
których jest
|
k0+k1+k2+⋯+kk-2= 2k-k-1 . |
| (10.4) |
Etap ∅. Szukanie wśród portfeli
x∈Δk:x1>0,x2>0,…,xk>0.
Warunek komplementarności pociąga wtedy wektorową równość
co łącznie z dwoma ograniczeniami równościowymi w problemie
zapisujemy przy pomocy tzw. macierzy Lagrange'a, tu wymiaru k+2:
Dawne założenie (5.2) oczywiście wynika z (10.3),
więc dostajemy tu jedyne rozwiązanie x,λ,λEE,
na temat którego rozwiązujemy następnie układ nierówności
x1E>0,x2E>0,…,xkE>0.
Wynikiem jest jakiś przedział otwarty wartości E (być może pusty)
i odpowiadające portfele xE. Macierz Lagrange'a, która tu została
przywołana i użyta, to tylko przeformułowanie i pewne uzwarcenie metody
Blacka i noblistów z Wykładu VI. Portfele xE ewentualnie
wyłaniające się w tym etapie są tożsame z portfelami Blacka (6.2)
z Wykładu VI leżącymi wewnątrz sympleksu Δk. Na dalszych etapach
tak już być nie musi (i najczęściej nie będzie); jednak patrz też
Ćwiczenia 11.1 i 11.2 oraz Uwaga 11.2 w Wykładzie XI. Uzwarcenie
metody dawane przez macierz Lagrange'a będzie nam pomocne w
dalszych etapach.
Wreszcie nazwa tego etapu, ∅, oznacza, że pusty jest
w nim zbiór indeksów zmiennych, które przyjmujemy za zero.
Ogólniej, etap OUT, gdzie
OUT⊂1, 2,…,k,
0≤#OUT≤k-2, polega na szukaniu
x∈Δk:xj=0 dla j∈OUT, natomiast xi>0 dla
i∈IN=def1, 2,…,k\OUT,
spełniających warunki w omawianej tu wersji twierdzenia K-KT, tzn.
|
Σx+λe-λEμi=0 dla i∈IN,Σx+λe-λEμi≥0dla i∈OUT |
| (10.5) |
dla jakichś λ,λE∈R. Równości w (10.5)
wynikają z warunku komplementarności, zaś nierówności w (10.5)
wynikają z warunku na gradient minimalizowanej funkcji. Oznacza to
konkretnie, że w takim etapie OUT najpierw rozwiązujemy układ
2+#IN≥2+2=4 równań liniowych
|
ΣINeINμINeINT00μINT00xINλ-λE=01E, |
| (10.6) |
gdzie ΣIN to macierz Σ po wyrzuceniu wierszy i
kolumn o numerach z OUT oraz eIN, μIN
to wektory e,μ po odrzuceniu składowych o numerach z OUT.
Ważne jest, że ΣIN>0 oraz μIN∦eIN, a więc stosuje się teoria Blacka
i współautorów, tyle, że w wymiarze #IN, nie zaś k.
Tu właśnie pracuje zestaw założeń (10.3), Σ>0
i #IN≥2.
Ćwiczenie 10.3
Uzasadnić, że istotnie, założenia Σ>0 i #IN≥2
pociągają ΣIN>0.
Uzasadnić też, że założenia (10.3) i #IN≥2
pociągają μIN∦eIN.
Tym sposobem wiemy, że rozwiązania xIN,λ,λEE
układu (10.6) istnieją i są jedyne. Z pomocą przyszedł nam,
dosyć niespodziewanie, Wykład VI.
Mając już te rozwiązania, rozwiązujemy następnie
ze względu na E układ nierówności
|
xiE>0,i∈IN,ΣxE+λEe-λEEμj≥0,j∈OUT. |
|
W wyniku otrzymujemy przecięcie k półprostych domkniętych i/lub
otwartych, czyli, ogólnie rzecz biorąc, przedział postaci
, lub , lub , lub , albo też ∅.
Definicja 10.1
Przedział ten będziemy w dalszym ciągu nazywać EIN;
uprości nam to w przyszłości przedstawienie kluczowego dla tych
wykładów Algorytmu Prostej Krytycznej (skrót nazwy angielskiej,
nieomalże powszechnie przyjętej, to: CLA).
Tak przebiegamy wszystkie 2k-k-1 etapów, w których
#IN≥2. (Przypominamy, że takie ograniczenie
z dołu na liczebność zbiorów IN było ważne i gwarantowało,
oczywiście przy założeniu (10.3), nierównoległość
zredukowanych wektorów μIN oraz eIN,
a w efekcie możliwość stosowania na każdym etapie OUT
klasycznej teorii Blacka i współautorów.)
Przykład 10.1 (ważny; wracamy do niego już po raz czwarty, za
każdym razem w innym kontekście, po Wykładach: V (Przykład 5.2),
VII (Ćwiczenie 7.3) i IX (Przykład 9.1))
Etap OUT=∅. Nieostre nierówności stają się
(warunek komplementarności!) równościami – w tym etapie szukamy
ewentualnych punktów prostej krytycznej leżących wewnątrz
sympleksu standardowego. Rozwiązujemy zatem równanie (patrz Twierdzenie 5.1)
|
0=9x1+3x2+x3513x1+2x241x1+4x321=10x1+6x3, |
|
oczywiście z zerowym skutkiem. Ten etap nic nie daje, bo
prosta krytyczna nie przechodzi tu przez wnętrze sympleksu:
E1, 2, 3=∅.
Etap OUT=1, w którym x1=0, x2>0, x3>0.
Piszemy odpowiedni układ równań (10.6)
|
2041042111004200x2x3λ-λE=001E, |
|
którego rozwiązania to
|
x2x3λλEE=12E-1-12E+25E-18-5+32E. |
|
Dodatniość x2 i x3 oznacza, że E∈2, 4. Co teraz
z pierwszą nierównością 3x2+x3+λ-5λE≥0,
która też musi być spełniona przez portfele relatywnie minimalnego
ryzyka? Po podstawieniu rozwiązań jest to nierówność 4≥E,
spełniona przez E ze wskazanego przedziału otwartego.
Wniosek: E2, 3=2, 4.
Etap OUT=2, w którym x1>0, x2=0, x3>0.
Odpowiedni układ równań (10.6) to teraz
|
9151142111005200x1x3λ-λE=001E, |
|
którego rozwiązania to z kolei
|
x1x3λλEE=13E-23-13E+53133x1-6x383x1-x3, |
|
gdzie celowo λ i λE zostały uzależnione `tylko'
od x1,x3. Dodatniość x1 i x3 oznacza, że E∈2, 5.
Co tym razem z drugą nierównością 3x1+λ-4λE≥0 ?
Po podstawieniu rozwiązań jej lewa strona to -103x1-2x3,
a więc wielkość ujemna. Wniosek: E1, 3=∅.
Jeśli chodzi o etap OUT=3, to analogiczne rachunki
dają w efekcie E1, 2=4, 5. Zatem, podsumowując, z
całego przedziału Emin,Emax=2, 5 dostajemy tutaj
2, 4∪4, 5. Jest to cały Emin,Emax bez
trzech węzłów wyróżnionych, w terminologii wprowadzanej
w Twierdzeniu 10.1 na końcu tego wykładu.
Ćwiczenie 10.4
Skoro, w tym jednym modelu, dla wszystkich, lub prawie wszystkich
wartości E znaleźliśmy portfele relatywnie minimalnego ryzyka w
aspekcie M, bez trudu możemy już wskazać portfele, które
przechodzą na granicę minimalną Fmin w tym modelu i aspekcie.
Należy to zrobić teraz, po czym porównać odpowiedź
z łamaną efektywną w tym samym modelu, będącą przedmiotem
Ćwiczenia 7.3 w Wykładzie VII. Zbiór portfeli wyłaniający się
z bieżącego ćwiczenia to przykład łamanej portfeli relatywnie
minimalnego ryzyka ilustrujący Twierdzenie 10.1 poniżej.
Zanim sformułujemy zapowiadane twierdzenie ogólne, przyjrzymy się
jeszcze autentycznie przełomowemu przykładowi modelu Markowitza
w wymiarze 4, pochodzącemu od K. Więcha [29], w którym
każdy etap OUT przynosi niepusty zbiór rozwiązań E.
Przykład ten był (tylko) anonsowany w Wykładzie VII. Mianowicie
|
Σ=9716-1618-287831516-1618294978-712-287849788858-1251631516-712-12516238,μ=322305. |
| (10.7) |
Wówczas po wykonaniu wszystkich 11=24-4-1 etapów
OUT uzyskuje się następujące rozwiązania:
Przedział wartości E |
Skrajne wartości parametru λE |
Etap OUT |
3, 5 |
λE3+=-114 |
2,3 |
λE5-=-2532 |
5,387 |
λE5+=-4780 |
1,2 |
λE387=-14232800 |
387,60201051 |
λE60201051=-19514204 |
1 |
60201051,82881287 |
λE82881287=-18015148 |
1,3 |
82881287,731789335 |
λE731789335=-1306774680 |
3 |
731789335,101771129 |
λE101771129=-94118064 |
∅ |
101771129,993101 |
λE933101=59914544 |
2 |
993101,64564 |
λE64564=8659216 |
2,4 |
64564,7499587 |
λE7499587=53419392 |
4 |
7499587, 22 |
λE22-=198 |
3,4 |
22, 30 |
λE22+=16764 |
1,4 |
λE30-=15532 |
Nawet w tej chwili (w roku akademickim 2009/10) przykład ten wygląda
wspaniale, a co dopiero w roku 2001, gdy się pojawił. Proszę np zwrócić
uwagę, jaką formalną symetrię środkową ma ciąg typów przedziałów
wartości E w lewej kolumnie:
Ćwiczenie 10.5 (sprawdzające)
Narysować wykres funkcji
której węzły wykresu są zakodowane w powyższej tabeli. (Pomiędzy
węzłami, jak wynika z opisu etapów algorytmu, funkcja λE
jest liniowa. Ogólnie jest więc ona kawałkami liniowa.)
Rozwiązanie:
Zadanie to jest rozwiązane w ostatnim Wykładzie XV (Rysunek 15.2).
Wszystkie dane do rysunku są jednak dostępne już teraz i każdy
czytelnik może wykonać własną wersję wykresu.
Ogólne twierdzenie na temat zbiorów portfeli relatywnie minimalnego
ryzyka w modelach Markowitza musi … poczekać jeszcze chwilę,
bo właśnie teraz nadarza się najlepsza okazja, by ugruntować
jednocześnie i model Więcha i twierdzenie K-KT.
Portfeli minimalnego ryzyka szukaliśmy dotąd przy pomocy tego twierdzenia
w modelach ± doskonale skorelowanych (Ćwiczenia 10.1 i 10.2 powyżej).
Natomiast w modelach niezdegenerowanych – jak ten koronny przykład
Więcha – jeszcze nie.
Ćwiczenie 10.6
W modelu Więcha znaleźć, oczywiście w aspekcie M, portfel
x~min o najmniejszej wariancji. (Jest to równocześnie,
warto wiedzieć, portfel efektywny w aspekcie M o najmniejszej
wartości oczekiwanej; więcej o tych sprawach będzie w Wykładzie XI.)
Wskazówka:
Zastosować algorytm bazujący na twierdzeniu K-KT, tym razem z jednym
ograniczeniem równościowym eTx=1 zamiast dwóch i z
czterema ograniczeniami nierównościowymi xi≥0,
i=1, 2, 3, 4.
Rozwiązanie:
Inne niż sugerowane we wskazówce rozwiązanie, mniej eleganckie
(zo to wykorzystujące Tabelę powyżej), jest podane w Przykładzie 11.2.
Twierdzenie 10.1 (nt zbioru portfeli relatywnie minimalnego ryzyka w modelach Markowitza)
∙ Przy założeniach Σ>0 oraz (10.3), portfele
relatywnie minimalnego ryzyka w modelu Markowitza tworzą łamaną Ł
o nie więcej niż 2k-k-1 rozłącznych bokach. Każdy bok jest
odcinkiem lub lub lub , i każdy leży
w innej ścianie sympleksu Δk (będącej krawędzią, bądź trójkątem,
bądź czworościanem, bądź …, bądź wnętrzem Δk). Portfele
na różnych bokach Ł mają więc różne składy jakościowe, a na
danym boku niezmienny skład jakościowy.
∙A∙ Nie wszystkie końce boków łamanej Ł należą do
Ł. Takich wyróżnionych przez nienależenie wierzchołków Ł zawsze
jest nie mniej niż dwa i nie więcej niż k – są to niektóre
z wierzchołków sympleksu Δk, przy czym zawsze – wierzchołki
o najmniejszej i największej wartości oczekiwanej. Te wyróżnione
wierzchołki łamanej Ł odpowiadają wartościom parametru
E∈min1≤i≤kμi,max1≤i≤kμi,
które pozostają nieobsłużone po wykonaniu wszystkich 2k-k-1
etapów OUT dla OUT⊂1, 2,…,k,
#OUT≤k-2.
Takie wartości E nazywamy węzłami wyróżnionymi.
∙A∙R∙ Po domknięciu w wierzchołkach wyróżnionych
łamana Ł¯ jest łamaną spójną. Ł, czy też częściej
Ł¯, jest nazywana łamaną wierzchołkową
w danym modelu Markowitza (dużo rzadziej nazywa się ją ”łamaną
portfeli relatywnie minimalnego ryzyka”).