<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Моделирование потоков платежей с помощью марковских случайных процессов</article-title>
      </title-group>
      <fpage>111</fpage>
      <lpage>117</lpage>
      <abstract>
        <p>1 - Уральский государственный университет путей сообщения (Екатеринбург) 2 - Уральский федеральный университет (Екатеринбург) Рассматриваются вероятностные модели потока платежей для инвестиционного проекта и для кредитного портфеля. Описаны подходы к получению оценок ожидаемой доходности и риска проекта и кредитного портфеля на основе расчета статистических моментов многошаговой системы, зависящей от марковской цепи. Предложены модели, учитывающие влияние макроэкономических факторов. Ключевые слова: поток платежей; оценивание; марковский процесс; чистая приведенная стоимость.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>что этот доход заранее не известен и зависит от марковского случайного процесса с дискретным временем
ξt = {ξ(1), . . . , ξ(t)}.</p>
      <p>При сравнении инвестиционных проектов могут быть использованы различные финансовые показатели:
чистая приведенная стоимость проекта – Net Present Value (NPV), внутренняя норма доходности (IRR),
срок окупаемости (PBP), эффективный индекс доходности (DPI) и другие показатели [4]. В данной статье
доходность проекта будем оценивать на основе его чистой приведенной стоимости (NPV), где [5]:
здесь r(t) – коэффициент дисконтирования, то есть доходность безрисковых вложений за тот же период.</p>
      <p>
        Будем рассматривать модель динамики N P V (t) инвестиционного проекта в форме многошаговой
системы. Обозначим через B(t, ξt) значение N P V (t, ξt). Из определения (1) чистой приведенной стоимости
проекта получаем, что многошаговая система [6]
(1)
(2)
(
        <xref ref-type="bibr" rid="ref1">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">4</xref>
        )
B(t + 1, ξt+1) = B(t, ξt) + a(t + 1)A(t + 1, ξt+1), B(0) = A(0),
      </p>
      <p>a(t + 1) = (1 + r(t))−1a(t), a(0) = 1.
описывает изменение чистой приведенной стоимости и может быть основой для оценки статистических
параметров (среднего значения, дисперсии, квантили) B(T, ξT ) = N P V (T, ξT ) в случае, если получено
вероятностное описание дохода за период A(t) = A(t, ξt).</p>
      <p>Для сравнения чистой приведенной стоимости различных проектов или различных стратегий
управления проектом в случае моделирования доходов как случайных величин используются следующие критерии
[7]:
• Максимум средней ожидаемой прибыли max N P V (T ) за заданный промежуток времени T при
фиксированном значении начальных инвестиций A(0).
• Минимум риска. Под риском проекта в финансовом анализе понимается вероятность получения
убытков, то есть вероятность отрицательного значения N P V (T )</p>
      <p>α(T ) = P{N P V (T, ξT )} &lt; 0
• Сочетание этих двух подходов, как двух критериев. Выбирается эффективный проект, то есть
Паретооптимальный по двум критериям.
3</p>
      <p>Вероятностное моделирование дохода по проекту за период
Обычно доходы и расходы по инвестиционному проекту за каждый период описываются, как разность
доходов и расходов за фиксированный период. Например, в работах [7, 8] используется следующая формула:</p>
      <p>A(t) = b1(t)b2(t) − b3(t)u(t) − c(t),
где b1(t) – количество проданных товаров, b2(t) – цена продажи единицы продукции, b3(t) – затраты на
закупку единицы сырья, u(t) – объем закупаемого сырья, c(t) – условно-постоянные расходы.</p>
      <p>Будем рассматривать модель, в которой параметры b(t) = {b1(t), b2(t), b3(t)} являются случайными и
зависят от состояния экономики. Изменение состояния экономики будем описывать, как марковскую цепь
ξ(t) с 3-мя возможными состояниями:
• S1 – кризисное состояние,
• S2 – равновесное состояние,
• S3 – экономический подъем.
Обозначим через e1, e2, e3 базисные вектора, пусть случайный процесс ζ(t) является индикатором
состояния экономики, т.е. ζ(t) = ej, если экономика находится в j-ом состоянии.</p>
      <p>Будем предполагать, что каждый из параметров проекта можно представить в виде двух слагаемых:
bij среднего значения i-го параметра (например, цены продажи) в j-ом состоянии экономики и случайного
отклонения от среднего ожидаемого значения ciηi(t). Например, количество проданного товара (b1(t))
моделируется как b1(t) = b12 + c1η1(t), если экономика находится в равновесном состоянии, т.е. ζ(t) =
{0, 1, 0}T = e2. Вектор параметров b(t) при условии ζ(t) = ej имеет вид</p>
      <p>
        bj(t) = Bjξ(t) + Cη(t),
где η(t), t = 1, 2, ... – дискретный случайный процесс с независимыми значениями и стандартным
нормальным распределением компонент. Сокращенно математическую модель изменения параметров проекта
можно записать в форме
b(t, ξ) = Bζ(t) + Cη(t),
(
        <xref ref-type="bibr" rid="ref3">5</xref>
        )
где B – заданная матрица размера 3 × 3.
      </p>
      <p>Несмотря на то, что параметры b1 и b2 входят выражение для A(t, ζ(t), η(t)) мультипликативно, в
рассматриваемом случае независимых марковских случайных процессов ζt, ηt можно найти условное
математическое ожидание и дисперсию Aj(t) = A(t|ζ(t) = ej) при условии, что экономика в момент t находится
в j-ом состоянии. Отсюда, используя соотношения для моментов многошаговой системы, зависящей от
простой марковской цепи [9], получаем прогноз для статистических моментов состояния системы (2).
Таким образом, находим прогноз статистических моментов N P V (T, uT ). Обозначим через mT = m(T, uT )
математическое ожидание чистой приведенной стоимости проекта при заданной стратегии закупки сырья
uT = {u(1), . . . , u(T )}, через sT = s(T, uT ) – среднее квадратичное отклонение этого параметра.</p>
      <p>Имитационное моделирование с помощью программы MathCad показало, что распределение чистой
приведенной стоимости проекта в рассматриваемых условиях близко к нормальному, что позволяет оценивать
риск проекта на основе функции нормального распределения
α(T, uT ) = P{N P V (T, ξT )} &lt; 0 = 0.5 − Φ
mT
sT
где Φ(z) – интегральная функция Лапласа, P(B) – вероятность события B.</p>
      <p>
        В случае заданной матрицы переходных вероятностей марковской цепи ζt для выбранного управления
uT соотношения (
        <xref ref-type="bibr" rid="ref2">4</xref>
        )–(
        <xref ref-type="bibr" rid="ref3">5</xref>
        ) вместе с уравнением потока платежей (2) и описанием динамики марковской цепи
определяют ожидаемую доходность и риск проекта. Однако, матрица переходных вероятностей известна
в данном случае не точно, кроме того, стоит задача выбора оптимального управления (стратегии закупки
сырья) uT , которая пока не решена.
      </p>
      <p>Модель динамики структуры кредитного портфеля
В работе [6] предложено описание потоков платежей по кредитному портфелю, в форме случайного
процесса зависящего от марковской цепи ξt = {ξ(0), . . . , ξ(t)} с k состояниями и дискретным временем. При
моделирование динамики структуры портфеля потребительских кредитов, весь портфель разбивается на
группы по признаку наличия или отсутствия задолженности и ее срокам. Под состоянием кредита
подразумевается принадлежность той или иной группе. Приведем пример простейшего разбиения кредитов на
группы для модели марковской цепи с шагом по времени равным 1 кварталу.</p>
      <p>1. кредиты без задержек платежей, в том числе новые (S1);
2. кредиты с задержкой по выплате менее 90 дней (S2);
3. кредиты с задержкой по выплате от 91 до 180 дней (S3);
4. кредиты с задержкой более 180 дней – проблемные (S4);
5. погашенные кредиты (S5).
Обозначим вероятность того, что система находится в i-ом состоянии в момент t через xi(t):
Рис. 2: Граф состояний системы
xi(t) = P{ξ(t) = ei}, i = 1, . . . , k,
{e1, . . . , ek} базис Rk.</p>
      <p>Обозначим через pij вероятность перехода из состояния Si в состояние Sj за один шаг</p>
      <p>P{ξ(t + 1) = ej|ξ(t) = ei} = pij.</p>
      <p>Обозначим через P = {pij} k × k матрицу переходных вероятностей и x(t) = {x1(t), . . . , xk(t)}&gt; –
вектор вероятностей состояний (долей портфеля). В случае постоянных переходных вероятностей динамика
вектора x(t) описывается уравнением
x(t + 1) = P &gt;x(t),</p>
      <p>t = 0, 1, . . . , T.</p>
      <p>x(t + 1) = P &gt;(Y (t), U (t))x(t),
Следуя подходу, предложенному [10], в работе [11] для описания динамики структуры кредитного
портфеля используется модель c переменной матрицей переходных вероятностей
где матрица переходных вероятностей P зависит от внешних экономических факторов Y (t) и действий
менеджмента банка U (t).</p>
      <p>
        Можно использовать другой вариант описания зависимости структуры кредитного портфеля от
состояния экономики. Будем полагать, что в зависимости от состояния экономики Sj матрица переходных
вероятностей принимает одно из 3-х возможных значений: P (1), P (2) или P (
        <xref ref-type="bibr" rid="ref1">3</xref>
        ). Таким образом, изменение
структуры портфеля будет описываться системой со случайной структурой
      </p>
      <p>
        x(t + 1) = P &gt;(ζ(t), U )x(t),
где ζ(t) – марковская цепь, описывающая изменение состояния экономики.
(
        <xref ref-type="bibr" rid="ref4">6</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">7</xref>
        )
(
        <xref ref-type="bibr" rid="ref6">8</xref>
        )
      </p>
      <p>Поток платежей по кредитному портфелю
Рассмотрим поток платежей по кредитному договору с постоянными выплатами (аннуитетными
платежами) с начальной суммой D(0). Обозначим ежеквартальный платеж по договору через d = dbD(0), где
коэффициент db определяется по известной формуле [5] исходя из процентной ставки по кредиту b.</p>
      <p>Через D(t, ξt) обозначим сумму основного долга в период t. Платеж A(t, ξt) в период t + 1 является
случайной величиной и связан с переходом договора из одной группы в другую. Если номер группы не
меняется в период t+1 (за исключением кредитов из проблемной группы), то это означает, что аннуитетный
платеж d оплачен:</p>
      <p>A(t + 1, ξt+1) = d, D(t + 1, ξt+1) = (1 + b)D(t, ξt) − d.
Если платеж не оплачен, то номер группы увеличивается на 1 или не изменяется для группы проблемных
кредитов:</p>
      <p>A(t + 1, ξt+1) = 0, D(t + 1, ξt+1) = (1 + b)D(t, ξt),
и т.д. В случае погашения кредита в период t + 1 кредит переходит в последнюю группу (погашенных
кредитов)</p>
      <p>A(t + 1, ξt+1) = (1 + b)D(t, ξt), D(t + 1, ξt+1) = 0.</p>
      <p>Получаем многошаговую систему, описывающую платежи по кредиту A(t, ξt) и изменение основного
долга D(t, ξt), вида [6]:</p>
      <p>
        A(t + 1) = ξ(t)&gt;(d · C1 + (1 + b)D(t)C2)ξ(t + 1), t = 0, . . . , T,
D(t + 1) = (1 + b)D(t) − ξ(t)&gt;(d · C1 + (1 + b)D(t)C2)ξ(t + 1), D(0) = A(0),
(
        <xref ref-type="bibr" rid="ref7">9</xref>
        )
(
        <xref ref-type="bibr" rid="ref8">10</xref>
        )
где C1 и C2 – заданные k × k матрицы, учитывающие, в том числе, выплату штрафов за просроченные
платежи.
      </p>
      <p>
        Уравнения (
        <xref ref-type="bibr" rid="ref7">9</xref>
        )–(
        <xref ref-type="bibr" rid="ref8">10</xref>
        ) вместе с уравнением потока платежей (2) и описанием динамики долей кредитного
портфеля (
        <xref ref-type="bibr" rid="ref4">6</xref>
        ) определяют распределение случайной величины B(T ) = N P V (T ). Если матрица переходных
вероятностей P задана, то мы можем найти среднее и дисперсию чистой приведенной стоимости портфеля
(N P V (T )) и её , используя соотношения из статьи [9] для моментов линейной многошаговой стохастической
системы, зависящей от простой марковской цепи. На основе этих параметров оценивается риск дефолта
по портфелю, т.е. вероятность α(T ) события N P V (T ) &lt; 0.
      </p>
      <p>
        Если при анализе динамики структуры портфеля выявлены зависимости от внешнеэкономических
факторов в форме (
        <xref ref-type="bibr" rid="ref5">7</xref>
        ), то поток платежей по кредитному портфелю, а также оценка ожидаемой доходности и
риска портфеля зависят от прогноза макроэкономических показателей YˆT = Yˆ (t), . . . , Yˆ (T ). В этом случае
для прогнозирования чистой приведенной стоимости портфеля удобнее применять имитационное
моделирование, с помощью которого получаются оценки статистических моментов случайной величины B(t, YˆT ),
удовлетворяющей соотношениям (1), (
        <xref ref-type="bibr" rid="ref7">9</xref>
        )–(
        <xref ref-type="bibr" rid="ref8">10</xref>
        ), (
        <xref ref-type="bibr" rid="ref5">7</xref>
        ).
Заключение
В статье предложена модель прогнозирования чистой приведенной стоимости проекта для случая, когда
поток платежей по проекту описывается марковским случайным процессом, зависящим от дискретной
марковской цепи связанной с состоянием экономики.
      </p>
      <p>Описан подход к получению оценок ожидаемой доходности и риска проекта и кредитного портфеля на
основе расчета статистических моментов стохастической многошаговой системы, учитывающей прогноз
макроэкономических факторов.
Список литературы
[1] G. Timofeeva, N. Timofeev. Evaluation of Payment Flows Based on Markov Chain Model with</p>
      <p>
        Incomplete Information. AIP Conference Proceedings, 631:17–22, 2014.
[2] V. I. Klokov, S. I. Kichko. Probabilistic assessment of the effectiveness of investment projects with
use of discounting in risk. Ekonomika i upravlenie, (
        <xref ref-type="bibr" rid="ref1">3</xref>
        ):82–87, 2009. (in Russian) = В. И. Клоков,
С. И. Кичко. Вероятностная оценка эффективности инвестиционных проектов с использованием
дисконтирования в условиях риска. Экономика и управление, (
        <xref ref-type="bibr" rid="ref1">3</xref>
        ):82–87, 2009.
[11] G. Timofeeva, Ya. Bozhalkina. Markov model of the loan portfolio dynamics considering influence of
management and external economic factors. AiP Conference Proceedings, 1789(020023), 2016.
      </p>
      <p>Modelling of cash flows by means of Markov processes</p>
    </sec>
    <sec id="sec-2">
      <title>Galina A. Timofeeva</title>
      <p>Ural State University of Railway Transport (Yekaterinburg, Russia)
Ural Federal University (Yekaterinburg, Russia)</p>
    </sec>
    <sec id="sec-3">
      <title>Yana A. Bozhalkina</title>
      <p>Ural State University of Railway Transport (Yekaterinburg, Russia)</p>
      <p>Abstract. Probabilistic models of cash flow for an investment project and for a loan portfolio are considered.
The approaches to the estimation of profitability and risk of the project and the loan portfolio are obtained on
the basis of calculation of statistical moments multistep system, which depends on the Markov chain. Models,
taking into account the changes in macroeconomic factors, are proposed.</p>
      <p>Keywords: payment flow, evaluation, Markov process, net present value.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>D.</given-names>
            <surname>Fantazzini</surname>
          </string-name>
          .
          <source>Credit Risk Management. Applied Econometrics</source>
          ,
          <volume>12</volume>
          (
          <issue>4</issue>
          ):
          <fpage>105</fpage>
          -
          <lpage>138</lpage>
          ,
          <year>2008</year>
          . = Д.
          <article-title>Фантацци- ни. Управление кредитным риском</article-title>
          .
          <source>Прикладная эконометрика</source>
          ,
          <volume>12</volume>
          (
          <issue>4</issue>
          ):
          <fpage>105</fpage>
          -
          <lpage>138</lpage>
          ,
          <year>2008</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>E. V.</given-names>
            <surname>Bucenko</surname>
          </string-name>
          .
          <article-title>Development of an expert system for investment design</article-title>
          .
          <source>Ekonomicheskie issledovanija</source>
          , (
          <volume>3</volume>
          ):
          <fpage>3</fpage>
          -
          <lpage>15</lpage>
          ,
          <year>2012</year>
          .
          <article-title>(in Russian) = Е</article-title>
          . В. Буценко.
          <article-title>Разработка экспертной системы инвестиционного проектирования</article-title>
          .
          <source>Экономические исследования</source>
          , (
          <volume>3</volume>
          ):
          <fpage>3</fpage>
          -
          <lpage>15</lpage>
          ,
          <year>2012</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [5]
          <string-name>
            <surname>V. I. Shiryaev.</surname>
          </string-name>
          <article-title>Finansovaja matematika: Potoki platezhej</article-title>
          ,
          <source>proizvodnye finansovye instrumenty [Financial Mathematics: Payment Flows</source>
          , Derivative FinancialInstruments]. Moscow, USSR.RU,
          <year>2011</year>
          .
          <article-title>(in Russian) = В</article-title>
          . И. Ширяев.
          <article-title>Финансовая математика: Потоки платежей, производные финан- совые инструменты</article-title>
          . Москва, USSR.ru,
          <year>2011</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>G. A.</given-names>
            <surname>Timofeeva</surname>
          </string-name>
          .
          <article-title>Estimation of the dynamics of the process depending on the Markov chain with incomplete information. 12-oe Vserossijskoe soveshhanie po problemam upravlenija</article-title>
          .
          <source>Trudy [ The 12th All-Russian Conference on Control. Proceedings]</source>
          . Moscow, ISS RAS:
          <fpage>1225</fpage>
          -
          <lpage>1230</lpage>
          ,
          <year>2014</year>
          .
          <article-title>(in Russian) = Г</article-title>
          . А. Тимофеева.
          <article-title>Оценивание динамики процесса, зависящего от марковской цепи с неполной информацией. 12-ое Всероссийское совещание по проблемам управления</article-title>
          .
          <source>Труды. Москва</source>
          , ИПУ РАН:
          <fpage>1225</fpage>
          -
          <lpage>1230</lpage>
          ,
          <year>2014</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>G. M.</given-names>
            <surname>Vakulina</surname>
          </string-name>
          <article-title>Probabilistic models in the investment projects assessment</article-title>
          .
          <source>Herald of USURT</source>
          ,
          <volume>9</volume>
          (
          <issue>1</issue>
          ):
          <fpage>93</fpage>
          -
          <lpage>100</lpage>
          ,
          <year>2011</year>
          .
          <article-title>(in Russian) = Г</article-title>
          . М. Вакулина.
          <article-title>Вероятностные модели в оценке инвестиционных про- ектов</article-title>
          .
          <source>Вестник УрГУПС</source>
          ,
          <volume>9</volume>
          (
          <issue>1</issue>
          ):
          <fpage>93</fpage>
          -
          <lpage>100</lpage>
          ,
          <year>2011</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>D. S.</given-names>
            <surname>Zavalishchhin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T. V.</given-names>
            <surname>Zav</surname>
          </string-name>
          <article-title>'yalova. Application of stochastic modeling to evaluate the effectiveness of an investment project</article-title>
          .
          <source>Informacionnye tehnologii modelirovanija i upravlenija</source>
          ,
          <volume>101</volume>
          (
          <issue>5</issue>
          ):
          <fpage>371</fpage>
          -
          <lpage>381</lpage>
          ,
          <year>2016</year>
          .
          <article-title>(in Russian) = Д</article-title>
          . С. Завалищин, Т. В. Завьялова.
          <article-title>Применение стохастического моделирования для оценки эффективности инвестиционного проекта</article-title>
          .
          <source>Информационные технологии моделирования и управления</source>
          ,
          <volume>101</volume>
          (
          <issue>5</issue>
          ):
          <fpage>371</fpage>
          -
          <lpage>381</lpage>
          ,
          <year>2016</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>T. V.</given-names>
            <surname>Zav</surname>
          </string-name>
          <article-title>'yalova, I. Ya</article-title>
          . Kats,
          <string-name>
            <given-names>G. A.</given-names>
            <surname>Timofeeva</surname>
          </string-name>
          .
          <article-title>Stochastic System with a Ramdom Jump in Phase Trajectoty: Stability of Their Motion</article-title>
          .
          <source>Automation and Remote Control</source>
          ,
          <volume>63</volume>
          (
          <issue>7</issue>
          ):
          <fpage>1070</fpage>
          -
          <lpage>1079</lpage>
          ,
          <year>2002</year>
          . = Т. В. Завьялова, И. Я. Кац, Г. А. Тимофеева.
          <article-title>Об устойчивости систем со случайным условием скачка фазовой трактории</article-title>
          .
          <source>Автоматика и телемеханика</source>
          , (
          <volume>7</volume>
          ):
          <fpage>33</fpage>
          -
          <lpage>43</lpage>
          ,
          <year>2002</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>A.</given-names>
            <surname>Caja</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Q.</given-names>
            <surname>Guibert</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Planchet</surname>
          </string-name>
          .
          <article-title>Influence of Economic Factors on the Credit Rating Transitions and Defaults of Credit Insurance Business</article-title>
          . Working Papers,
          <fpage>hal</fpage>
          -
          <lpage>01178812</lpage>
          , НAL,
          <year>2015</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>