simBySolution

Симулируйте приближенное решение диагонального дрейфа процессы HWV

свернуть все на странице

Синтаксис

[Paths,Times,Z] = simBySolution(MDL,NPeriods)
[Paths,Times,Z] = simBySolution(___,Name,Value)

Описание

пример

[Paths,Times,Z] = simBySolution(MDL,NPeriods) симулирует приближенное решение диагонального дрейфа для Гауссовой Диффузии Hull-White/Vasicek (HWV) процессы.

пример

[Paths,Times,Z] = simBySolution(___,Name,Value) добавляют дополнительные аргументы пары "имя-значение".

Примеры

свернуть все

Создайте hwv объект представлять модель:

dXt=0.2(0.1−Xt)dt+0.05dWt.

hwv = hwv(0.2, 0.1, 0.05)  % (Speed, Level, Sigma)
hwv = 

   Class HWV: Hull-White/Vasicek
   ----------------------------------------
     Dimensions: State = 1, Brownian = 1
   ----------------------------------------
      StartTime: 0
     StartState: 1
    Correlation: 1
          Drift: drift rate function F(t,X(t)) 
      Diffusion: diffusion rate function G(t,X(t)) 
     Simulation: simulation method/function simByEuler
          Sigma: 0.05
          Level: 0.1
          Speed: 0.2

simBySolution функция симулирует вектор состояния Xt с помощью приближения решения закрытой формы диагонального дрейфа HWV модели. Каждый элемент вектора состояния Xt описывается как сумма NBrowns коррелированые Гауссовы случайные ничьи добавляются к детерминированному переменному временем дрейфу. 

NPeriods = 100
[Paths,Times,Z] = simBySolution(hwv, NPeriods,'NTrials', 10);

Входные параметры

свернуть все

Режим Hull-White/Vasicek (HWV) в виде hwv объект, который создается с помощью hwv.

Типы данных: object

Количество периодов симуляции в виде положительного скалярного целого числа. Значение NPeriods определяет количество строк симулированного выходного ряда.

Типы данных: double

Аргументы name-value

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: [Paths,Times,Z] = simBySolution(HWV,NPeriods,'DeltaTime',dt,'NTrials',10)

Симулированные испытания (демонстрационные пути) NPeriods наблюдения каждый в виде разделенной запятой пары, состоящей из 'NTrials' и положительное скалярное целое число.

Типы данных: double

Положительное время постепенно увеличивается между наблюдениями в виде разделенной запятой пары, состоящей из 'DeltaTimes' и скаляр или NPeriods- 1 вектор-столбец.

DeltaTime представляет знакомый dt, найденный в стохастических дифференциальных уравнениях, и определяет времена, в которые сообщают о симулированных путях переменных состояния вывода.

Типы данных: double

Количество промежуточных временных шагов в течение каждого раза постепенно увеличивает dt (заданный как DeltaTime) в виде разделенной запятой пары, состоящей из 'NSteps' и положительное скалярное целое число.

simBySolution функциональные разделы каждый раз постепенно увеличивают dt в NSteps подынтервалы длины dt/NSteps, и совершенствовал симуляцию путем оценки симулированного вектора состояния в NSteps − 1 промежуточные точки. Несмотря на то, что simBySolution не сообщает вектор состояния вывода в этих промежуточных точках, улучшение улучшает точность, позволяя симуляции более тесно аппроксимировать базовый процесс непрерывного времени.

Типы данных: double

Отметьте, чтобы указать ли simBySolution использует прямо противоположную выборку, чтобы сгенерировать Гауссовы случайные варьируемые величины, которые управляют вектором Броуновского движения (винеровские процессы) в виде разделенной запятой пары, состоящей из 'Antithetic' и скалярный логический флаг со значением True или False.

Когда вы задаете True, simBySolution выполняет выборку таким образом, что все первичные и прямо противоположные пути симулированы и сохранены в последовательных парах соответствия:

  • Нечетные испытания (1,3,5,...) соответствуйте первичным Гауссовым путям.

  • Даже испытания (2,4,6,...) соответствующие прямо противоположные пути каждой пары, выведенной путем отрицания Гауссовых ничьих соответствующего первичного (нечетного) испытания.

Примечание

Если вы задаете входной процесс шума (см. Z), simBySolution игнорирует значение Antithetic.

Типы данных: логический

Прямая спецификация зависимого случайного шумового процесса раньше генерировала вектор Броуновского движения (винеровский процесс), который управляет симуляцией в виде разделенной запятой пары, состоящей из 'Z' и функция или как (NPeriods * NSteps)- NBrowns- NTrials 3D массив зависимых случайных варьируемых величин.

Входной параметр Z позволяет вам непосредственно задавать шумовой процесс генерации. Этот процесс более приоритетен по сравнению с Correlation параметр входа gbm возразите и значение Antithetic введите флаг.

Примечание

Если вы задаете Z как функция, это должно возвратить NBrowns- 1 вектор-столбец, и необходимо вызвать его с двумя входными параметрами:

  • Скалярное время наблюдения с действительным знаком t.

  • NVars- 1 вектор состояния Xt.

Типы данных: double | function

Отметьте, который указывает как выходной массив Paths хранится и возвратился в виде разделенной запятой пары, состоящей из 'StorePaths' и скалярный логический флаг со значением True или False.

Если StorePaths True (значение по умолчанию), или не задано, simBySolution возвращает Paths как 3D массив временных рядов.

Если StorePaths False (логический 0), simBySolution возвращает Paths выходной массив как пустая матрица.

Типы данных: логический

Последовательность процессов конца периода или корректировок вектора состояния формы в виде разделенной запятой пары, состоящей из 'Processes' и функциональный или массив ячеек функций формы

Xt=P(t,Xt)

simBySolution функционируйте запуски, обрабатывающие функции в каждый раз интерполяции. Они должны принять текущее время интерполяции t и вектор текущего состояния Xt, и возвратить вектор состояния, который может быть корректировкой состояния ввода.

simBySolution применяет функции обработки в конце каждого периода наблюдения. Эти функции должны принять текущее время наблюдения t и вектор текущего состояния X t, и возвратить вектор состояния, который может быть корректировкой состояния ввода.

Конец периода Processes аргумент позволяет вам отключать данное испытание рано. В конце каждого временного шага, simBySolution тестирует вектор состояния Xt на все-NaN условие. Таким образом, чтобы сигнализировать о раннем завершении данного испытания, всех элементах вектора состояния Xt должен быть NaN. Этот тест включает пользовательский Processes функционируйте, чтобы сигнализировать о раннем завершении испытания и предложениях значительные выигрыши в производительности в некоторых ситуациях (например, оценивая разоренные барьерные опционы).

Если вы задаете больше чем одну функцию обработки, simBySolution вызывает функции в порядке, в котором они появляются в массиве ячеек. Можно использовать этот аргумент, чтобы задать граничные условия, предотвратить отрицательные цены, накопить статистику, графики графика и т.д.

Типы данных: cell | function

Выходные аргументы

свернуть все

Симулированные пути коррелированых переменных состояния, возвращенных как (NPeriods + 1)- NVars- NTrials 3D массив временных рядов.

Для данного испытания, каждой строки Paths транспонирование вектора состояния X t во время t. Когда входной флаг StorePaths = False, simBySolution возвращает Paths как пустая матрица.

Времена наблюдения сопоставлены с симулированными путями, возвращенными как (NPeriods + 1)- 1 вектор-столбец. Каждый элемент Times сопоставлен с соответствующей строкой Paths.

Зависимые случайные варьируемые величины раньше генерировали вектор Броуновского движения (винеровские процессы), которые управляют симуляцией, возвращенной как (NPeriods * NSteps)- NBrowns- NTrials 3D массив временных рядов.

Больше о

свернуть все

Прямо противоположная выборка

Методы симуляции позволяют вам указывать, что популярный метод сокращения отклонения вызвал прямо противоположную выборку.

Этот метод пытается заменить одну последовательность случайных наблюдений с другим из того же ожидаемого значения, но меньшее отклонение. В типичной симуляции Монте-Карло каждый демонстрационный путь независим и представляет независимое испытание. Однако прямо противоположная выборка генерирует демонстрационные пути в парах. Первый путь пары упоминается как первичный путь и второе как прямо противоположный путь. Любая данная пара независима от любой другой пары, но эти два пути в каждой паре высоко коррелируются. Прямо противоположная литература выборки часто рекомендует составить в среднем обесцененные выплаты каждой пары, эффективно деля на два количество испытаний Монте-Карло.

Этот метод пытается уменьшать отклонение путем стимулирования отрицательной зависимости между парными входными выборками, идеально приведения к отрицательной зависимости между парными выходными выборками. Чем больше степень отрицательной зависимости, тем более эффективная прямо противоположная выборка.

Алгоритмы

simBySolution метод симулирует NTrials демонстрационные пути NVars коррелированые переменные состояния, управляемые NBrowns Источники броуновского движения риска по NPeriods последовательные периоды наблюдения, аппроксимируя непрерывное время Hull-White/Vasicek (HWV) приближением решения закрытой формы.

Рассмотрите отделимую, модель HWV с векторным знаком формы:

dXt=S(t)[L(t)−Xt]dt+V(t)dWt

где:

  • X является NVars-by-1 вектор состояния переменных процесса.

  • S является NVars-by-NVars матрица скоростей возвращения к среднему уровню (уровень возвращения к среднему уровню).

  • L является NVars-by-1 вектор из уровней возвращения к среднему уровню (отдаленное среднее значение или уровень).

  • V является NVars-by-NBrowns мгновенная матрица уровня энергозависимости.

  • W является NBrowns-by-1 Вектор броуновского движения.

simBySolution метод симулирует вектор состояния Xt с помощью приближения решения закрытой формы моделей диагонального дрейфа.

При выполнении выражений, simBySolution принимает, что все параметры модели являются кусочно-постоянными за каждый период симуляции.

В общем случае это не точное решение моделей, потому что вероятностные распределения симулированных и истинных векторов состояния идентичны только для кусочно-постоянных параметров.

Когда параметры являются кусочно-постоянными за каждый период наблюдения, симулированный процесс точен в течение времен наблюдения, в которые производится Xt.

Гауссовы модели диффузии, такие как hwv, позвольте отрицательные состояния. По умолчанию, simBySolution не делает ничего, чтобы предотвратить отрицательные состояния, и при этом это не гарантирует, что модель строго возвращается среднее значение. Таким образом модель может показать ошибочный или взрывной рост.

Ссылки

[1] Aït-Sahalia, Yacine. “Тестируя Модели Непрерывного времени Точечной Процентной ставки”. Анализ Финансовых Исследований 9, № 2 (апрель 1996): 385–426.

[2] Aït-Sahalia, Yacine. “Плотность перехода для Процентной ставки и Другой Нелинейной Диффузии”. Журнал Финансов 54, № 4 (август 1999): 1361–95.

[3] Глассермен, Пол. Методы Монте-Карло в финансовой разработке. Нью-Йорк: Springer-Verlag, 2004.

[4] Оболочка, Джон К. Опции, фьючерсы и Другие Производные. 7-й редактор, Prentice Hall, 2009.

[5] Джонсон, Норман Ллойд, Сэмюэль Коц и Нэраянэсвами Бэлэкришнэн. Непрерывные Одномерные распределения. 2-й редактор Вайли Серис в Вероятности и Математической Статистике. Нью-Йорк: Вайли, 1995.

[6] Shreve, Стивен Э. Стохастическое исчисление для финансов. Нью-Йорк: Springer-Verlag, 2004.

Смотрите также

| | |

Темы

Введенный в R2008a

Документация Financial Toolbox

Поддержка

Для просмотра документации необходимо авторизоваться на сайте