Я хотел убедиться что SV модель симуляции цен 2х акций соответствует реальности.
Меня беспокоил случай который модель не могла повторить. Очень сильная корреляция в больших убытках, и отсутствие сильной корреляции в больших прибылях (или наоборот). Такую ассиметрию с высокими и меняющимися значениями корреляции модель повторить не может.
Но, оказалось таких случаев нет. Во первых а) даже при очень сильных корреляциях (как золото и серебро или кокакола и пепсикола) таки допускают некую случайность и разброс. б) и они боль менее симметричны.
Итого, что происходит в реальности — корреляция может меняеться в зависимости от волатильности (становиться выше), и это требуется учесть, нужно два параметра корреляции для обычного режима и экстремального (плюс пара параметров для корреляции скрытых процессов волатильности). И по идее этого должно быть достаточно.

Тем не менее, 4 параметра корреляции (2 для прибыли, 2 для волатильности), особенно корреляция хвостов — плохо идентифицируются. Поэтому практически для фиттинга — их подбор лучше сделать в полу ручном режиме. Посмотреть на графиках как корреляция выглядит, и выставить разумные приоры вручную, позволив модели немножко подстроить их.
чтобы она «правильно» жила — её надо отправить в собственный базис и «дать жить» двум собственным значениям и углу поворота
вообще любая положительно-определённая симметричная матрица приводима к С = P^T L P — где P — ортогональная и задаётся углами поворотов в плоскостях, а L — матрица собственных значений
короче, то что вы делали для логарифма сигмы^2 в одномерном случае (там ничто никуда поворачивать не надо) — для двумерного это логарифмы lambda_1 и lambda_2, правда появляется еще угол \phi, на который тоже нужно придумать динамику
с ходу не помню как там точно выразятся [ко]вариации через лямбды и фи, но в итоге у вас получится, что корреляция точно «случайно» не вылезет из [-1;+1]
пока вы пишите динамику для логарифмов лямбд — положительную определенность вы точно не потеряете
про динамику углов — проще всего конечно на неё забить — зафитить и сказать, что хз откуда, но вот так.
стоит понимать, что корреляция зависит не только от углов, но и от того, что у вас есть динамика лямбд — корреляция тоже будет меняться
посложнее — это отпустить углы в брауновское движение, не важно, что они перейдут через 2пи или минус, тут важнее отпустить их в довольно лёгкое движение, не для описани динамики, а просто, чтоб создать фреймворк для учета этого риска
минреверс для углов тоже наверное можно придумать, но как бы это перебор
вообще в n-мерных моделях n дополнительных динамик — это n^2 непонятных коэффициентов в подарок ))))
самое сложное это понять, что по пути не потерялась положительная определенность, а корреляции будут гулять из-за динамики лямбд, даже если вы фиксируете базис