Характеризация отношения к риску
Так, например, решением а может служить структура инвестиционного портфеля, состоянием среды и - доходности ценных бумаг, входящих в портфель, а результатом ? - доходность портфеля в целом.
Отношение к риску
Пусть X - совокупность всех рисков (для определенности - случайных величин); рассмотрим произвольный риск X G X с функцией распределения F(z) = P{X z} и математическим ожиданием рх = EX; в качестве меры полезности риска будем использовать u(X) - среднее значение некоторой функции полезности U на этом риске:
u(X) = EU (X). (43)
Исследуем связь формы функции полезности с отношением ее обладателя к риску. Значения X здесь будем трактовать, как доход (чем больше, тем лучше), а функцию полезности U считать возрастающей, так что индивидуум с данной функцией полезности стремится максимизировать значение u(-). В зависимости от соотношения u(X) и рх будем различать следующие варианты отношения индивидуума к риску:
- нейтральное отношение: u(X) = рх, X G X;
- склонность к риску: u(X) рх, X G X;
- неприятие риска: u(X) рх, X G X.
Будем обозначать классы функций полезности, описывающих нейтральное отношение, склонность к риску и неприятие риска Ur, Ur Цд, соответственно:
Un = {U : EU(X) = U(EX), X G X},
Ur = {U : EU(X) U(EX), X G X},
Ua = {U : EU(X) U(EX), X G X}.
5.1.1 Нейтралитет
Рассмотрим сначала линейную функцию полезности U (z) = az + b. В этом случае
u(X) = EU (X) = E(aX + b) = aEX + b = U (EX) = U (р)
ввиду линейности математического ожидания. Это означает, что линейная функция полезности описывает нейтральное отношение к риску и любая линейная функция U является элементом Ur.
5.1.2 Склонность к риску
Пусть теперь функция полезности U выпукла, т.е. удовлетворяет условию U(ay + (1 a)z) aU(y) + (1 a)U(z), a G [0,1],
и, более общо,
n n n
U(X az*) X aiU(z*), a* 0, Xa* = 1.
i= 1 i=1 i= 1
Тогда, полагая распределение X для простоты дискретным со значениями zi и соответствующими вероятностями pi; i = 1, получаем
EU(X) = / U(z) dF(z) = ? U(zi)pi U(E = U(EX) = U(д), (45)
i=1 i=1
т.е. выпуклая функция полезности описывает склонность к риску. Отметим здесь, что вывод формулы (45) справедлив для произвольного, а не только дискретного распределения X.
5.1.3 Неприятие риска
Рассмотрим теперь случай вогнутой функции полезности, удовлетворяющей условию
Здесь аналогично получаем неравенство
означающее, что вогнутые функции полезности описывают неприятие риска (risk aversion). Верно и обратное утверждение: неприятие риска описывается вогнутой функцией полезности; сформулируем его в виде теоремы.
Теорема 5.1 Если для произвольного риска X выполняется неравенство
EU(X) U(EX), X Е X
то функция U является вогнутой, т.е. удовлетворяет условию (46).
Замечание 5.1 Поскольку реальные участники рынков с рисками не приемлют риск, т.е. предпочитают детерминированный актив дх риску X с EX = дх, всюду в дальнейшем, будем рассматривать только строго вогнутые функции полезности. Без существенного ущерба для общности можно считать U достаточно гладкой функцией.
Таким образом, можно дать следующее определение:
Определение 5.1 Функцией полезности называется дважды непрерывно дифференцируем,ая функция, обладающая свойствами
U(0) = 0, U'(x) 0, U(x) 0, x 0. (48)
Упражнение 5.1 Доказать теорему 5.1.
Упражнение 5.2 Доказать аналогичную теорему:
EU (X) U (EX), X ЕХ^4 U eUr
Упражнение 5.3 Доказать:
EU (X) = U (EX), X ЕХ^4 U eUn
Количественное выражение неприятия риска
5.2.1 Цена риска
Пусть X - произвольный риск, w - начальный капитал индивидуума. Очевидно, w + X также является риском. Ясно, что если DX 0, то EU(X) U(EX) и, следовательно, для произвольного w 0 имеем
EU(w + X) U(w + EX), (49)
так что уравнение
EU(w + X) = U (w + EX - п) (50)
имеет единственное решение п 0, зависящее от риска X, начального капитала w и функции полезности U.
Определение 5.2 Решение п = n(X) = nwU(X) уравнения (50) называется ценой риска X.
Явное выражение для цены риска имеет вид
п = w + EX - U-1[EU(w + X)], (51)
откуда, с учетом (49) очевидна его положительность.
Представляет интерес зависимость цены риска от формы функции полезности. Введем сначала количественное понятие неприятия риска.
5.2.2 Неприятие риска
Определение 5.3 Неприятием риска (для заданной строго вогнутой функции полезности U) называется отношение
a(z) = -U(z)/U'(z).
Интересующая нас связь описывается теоремой Пратта [7]. Рассмотрим две функции полезности U1, U2 и обозначим a1,a2 соответствующие величины неприятия риска, п1, п2 - цены риска.
5.2.3 Теорема Пратта
Теорема 5.2 Для двух произвольных функций полезности U1, U2 следующие утверждения эквивалентны:
а) a1(z) a2(z), z 0;
б) n1(X) П2 (X), VX;
в) Существует функция T такая, что T'(z) 0, T(z) 0, z 0 м
U1(z)= T(U2(z)), z 0.
Доказательство проведем по схеме а) ^ в), б) ^ в), следуя изложению в [б]. 1. а) ^ в). Имеем а\(?) a2(z), z 0. Зададим функцию T выражением
T (v) = U1(U2-1())
(52)
и покажем, что она обладает требуемыми свойствами. Подстановкой v = U2(z) легко проверяется, что U(z) = T(U2(z)), поэтому остается проверить вогнутость T. Проделаем проверку прямым вычислением производных и определением их знака. Имеем:
U1 (U2-1(v)) U2 (u-1(v))
T' (v) = U1 (U2-1(v))(U2~1y(v)
U2 (U-1(v)) U1' (U-) (U-1(v))' - U1 (U-1(v)) U (U--1(v)) (U-\v))'
T (v)
U2 (U-1(v)
U2 (U2-1(W) U1' fa) U (U2-1(0 U2' (Ц!) 1 , rfU' („-1fv)1V
U2 (U2-1 (v)) - U2 (U2-1(v))A1U2A (v)4
U1 (z) fU' (z) U2'(z ) 1 U1 (z)
{a2(z) - a1(z)}
(U2(z))H U'(z) U2(z)J (U2(z))‘
что строго меньше 0 по предположению.
2. в) ^ а). Пусть U1 = T(U2) и T' 0, T'' 0. Отсюда U1 = T'U2 и U1' =
T'U2' + T''(U2)2, поэтому
T'U2' + T''(U2')2
U1'
T
- u2
a2
2)
U1
T' U2
так что, очевидно, a1 a2.
3. в) ^ б). Пусть снова U1 = T(U2) и T' 0, T'' 0. Поскольку
П = w + EX - U-1(EUi(w + X)), i = 1, 2,
имеем:
П1 - П2 = U2-1(EU2(w + X)) - Uf1 (EU1(w + X)) = U2-1(EY) - U-1(ET(Y)),
где Y = U2(w + X). По неравенству Йенсена имеем ET(Y) T(EY), а так как функция U-1 является возрастающей, то и U1-1(ET(Y)) U1-1(T(EY)). Таким образом,
П1 - П2 U-1 (EY) - U-1 (T(EY)) = U- 1 (EY) - U- 1 (U (U2-1 (EY))) = 0,
что и требовалось.
4. б) ^ в). Зададим снова T формулой (52), из п1 п2 получим с учстом U2 = T-1(U1), U2-1 = U-1 (T): U2-1(EY) U-1(ET(y)), откуда U-1(T(EY)) U-1(ET(Y)) и, следовательно, ввиду монотонности U-\ T(EY) ET(Y). Выполнение последнего соотношения для произвольной случайной величины Y эквивалентно вогнутости T.
Упражнение 5.4 Для функции полезности U(z) = 1 exp(az) (a 0) выписать меру неприятия риска и явное выражение для цены риска (рисковой премии) п. Какие функции данного класса лежат в Ua с-
Упражнение 5.5 То же для функции полезности U(z) = z" (a 0). Какие функции данного класса лежат в Ua с-
Упражнение 5.6 В теорем,е 5.2 утверждение в) можно прочитать так: если Ui = T (U2) и преобразов анис T таково, что T' 0,, T 0, то обладатель функции полезности Ui менее рискован. Однако, ввиду условия, T' 0 функция, T взаимно однозначна, следовательно, им,ест обратную T-1, так ч,то U2 = T-i (Ui) и функция, U2 должна бы описывать менее рискованного индивидуума.
В чем заключается, асимметрия свойств T, не позволяющая, сделать такой, вывод?
Упражнение 5.Т Теорема 5.2 доказана по схеме, требующей, четырех логических цепочек доказательств. Попробуйте найти, способ доказательства, использующий минимальное возможное (3) количество цепочек (например, а) ^ б) ^ в) ^ а)).
Упражнение 5.8 Проверить, является, ли, цена риска nwU(X) монотонной функцией, w при, фиксированных X, U. Каков характер монотонности?
Простейший процесс риска
Рассмотрим простейшую модель процесса риска в терминах следующей игры ([3], гл. XIV): в игре участвуют два игрока, в каждой партии первый игрок выигрывает единицу с вероятностью p G (0,1) и проигрывает единицу с вероятностью q = 1 p. Суммарный начальный капитал обоих игроков равен а, начальный капитал первого игрока равен z; здесь а, z - целые числа, 0 z а. Таким образом, процесс описывается уравнением
X (t) = z + ? Zi, (53)
i=i
где Z1, Z2,... последовательность независимых одинаково распределенных случайных величин с дискретным распределением,
1, с вероятностью р,
1, с вероятностью q =1 р.
(54)
Zi
Игра заканчивается, когда обнуляется капитал одного из игроков (капитал первого игрока становится равным 0 или а), что трактуется, как разорение соответствующего игрока. Ясно, что при р 1/2 игра невыгодна доя первого игр ока ввиду EZi 0, при р 1/2 - выгоднa (EZi 0), а при р =1/2 является нейтральной, справедливой: EZi = 0. Обозначим tz - момент разорения первого игрока: tz = min{t : X(t) = 0} (несобственная случайная величина), qz вероятность разорения первого игрока при начальном капитале z: qz = P{tz то}.
Отметим, что ввиду известной симметрии игры вероятность разорения qz второго игрока при начальном капитале первого, равном z, может быть вычислена формальной заменой а на z а и перестановкой р и q в выражении для qz.
6.2 Уравнение для вероятности разорения
Если начальный капитал первого игрока равен 0 или а, то игра не проводится и соответствующие значения вероятностей разорения равны
Если же 0 z а, то после первой партии капитал первого игрока принимает значение z + 1 с вероятностью p или значение z 1 с вероятностью q. Поэтому, по формуле полной вероятности,
(56) представляет собой разностное уравнение второго порядка с характеристическим уравнением
p p
корни которого равны 1 и q/p, соответственно. Обозначим второй корень
(58)
Q p'
6.3 Вычисление вероятностей разорения
Если p = 1/2, то у = 1, корни различны, и общее решение уравнения (56) имеет вид
qz = Ci + C2yz. (59)
Используя краевые условия (55), находим значения постоянных C1, C2 и соответ-
ствующее частное решение:
(60)
1
qz
1 ya 1 1 ' 2
Вычислим вероятность разорения второго игрока упомянутым формальным приемом. Перестановка p и q означает замену у на у-1, так что
y-(a-z) _ y-a yz 1 1 yz
qz 1 _ y-a ya _ 1 1 _ ya
Нетрудно заметить, что qz + qz = 1, так что в игре без ограничения времени с вероятностью 1 происходит разорение одного из игроков.
При p = 1/2 имеем y = 1, так что корни характеристического уравнения вещественны и совпадают. Второе линейно независимое решение уравнения (56) имеет в этом случае вид z, общее решение - qz = C1 + C2z, частное решение -
К тому же результату можно прийти, раскрывая неопределенность в (60) при y 1 по правилу Лопиталя. Вероятность разорения второго игрока имеет вид
а z z
Qz
а а
так что снова с вероятностью 1 происходит разорение одного из игроков.
Функция qz в (60) и (61) монотонно убывает от 1 до 0 при изменении z от 0 до а. Если p 1/2, то y 1, и (60) является выпуклой; если же p 1/2, то y 1, и (60) явля?тся вогнутой
Игра с бесконечно богатым противником
Изучим поведение вероятности разорения при а ^ го; для случая p 1/2 при произвольном фиксированном z имеем
lim qz = 1,
а^ж
что означает бесперспективность игры с бесконечно богатым противником в случае невыгодности (р 1/2) или нейтральности (р = 1/2) игры. При p 1/2 получаем
lim qz = yz = e-az,
а^-ж
т.е. зависимость вероятности разорения от начального капитала является эксиопеп-циальной с показателем а = ln у = ln(q/p) 0. Отметим, что условие положительности рисковой надбавки EZ 0 эквивалентно р 1/2.
Классический процесс риска
Классический процесс риска изучался на протяжении всего 20 века, начиная с работы Лундберга [4]. Уравнением этого процесса описывается динамический портфель страховой компании, банка, других финансовых организаций, являющихся перерас-пределителями финансовых потоков в окружении рискованной среды.
Среди других приложений можно упомянуть описание уровня воды в водохранилище.
Определение
Рассмотрим определение процесса риска на примере работы страховой компании. Пусть страховые премии поступают равномерным потоком с интенсивностью с, а в случайные моменты времени 0 T1 Т2 ... наступают страховые события, наносящие ущерб случайного размера Z\, Z2,..., соответственно. Тогда размер капитала компании в момент времени t при условии, что начальный капитал (в момент времени T = 0) равен ж, описывается выражением
_ N(t) _
i= 1
где
количество страховых событий, наступивших в интервале времени [0,t]. Поскольку моменты времени Ti, i = 1, 2,... случайны, случайными оказываются и промежутки времени между последовательными страховыми событиями
Случайный процесс вида (62) называется классическим процессом риска, если случайные величины ?і, i = 1, 2,... являются независимыми, одинаково распределенными и имеют показательное распределение с параметром Л 0:
случайные величины Zi, i = 1, 2,... также являются независимыми и одинаково распределенными и имеют функцию распределения
При этом, как известно [3], количество страховых событий N(t) имеет распределение Пуассона с параметром Лк
а накопленный размер страховых убытков
N(t) _
i=1
на интервале времени [0, t] является случайной величиной с так называемым составным распределением Пуассона, функция распределения которого имеет вид
где F*' означает k - кратную свертку функции распределения F с собой, т.е. функцию распределения суммы k независимых одинаково распределенных случайных величин с функцией распределения F.
В тех случаях, когда необходимо подчеркнуть зависимость значения процесса X(t) от случайного аргумента ш G Q, будем использовать обозначение X(w,t), в частности, отдельную траекторию процесса при фиксированном ш будем обозначать
Как видим, классический процесс риска вполне определяется значениями четырех параметров (ж, с, Л, Fgr), удовлетворяющих условиям
Произвольный классический процесс риска с фиксированными значениями параметров, удовлетворяющих условиям (70), будем обозначать
а совокупность всех классических процессов риска с такими параметрами -
(72)
X = {X(ж, с, Л, F7) : ж 0, с 0, Л 0, F~(0) = 0}.
Разорение процесса
Под разорением процесса (62) понимается достижение уровня 0, то есть событие
R(x) = R(x, Л, С, F~) = {(л Е О : 31 0, X(л, t) 0}, (73)
при этом моментом разорения называется случайная величина
т(х) = т(х, Л, С, F~) = min{t : X(t) 0}. (74)
Эта случайная величина зависит от параметров процесса (62) и может оказаться несобственной, с положительной вероятностью принимая значение то; такая ситуация соответствует траекториям, которые не разоряются на всей временной полуоси [, то).
Вероятностью разорения процесса (62) называется величина
Р{т(х) то} = P(R(x)),
т.е. вероятностная мера множества тех траекторий, которые разоряются за конечное время. Эта величина также является, очевидно, функцией параметров процесса, что будем подчеркивать обозначением
R(х) = R(х, Л, С, F-) = Р{т(х) то}. (75)
В некоторых случаях более удобной характеристикой процесса риска оказывается вероятность выживания процесса
5(х) = 5(х, Л, С, F~) = 1 RR(s). (76)
Зависимость вероятности разорения процесса от параметров
Отметим следующие свойства монотонности вероятности разорения, как функции параметров процесса.
- R является невозрастающей функцией х;
- R является невозрастающей функцией С;
- R является неубывающей функцией Л;
- R является невозрастающей функцией F~, если порядок на множестве функций
распределения задан отношением F2 ^\(х) ^2(х) ?х.
Проверка перечисленных свойств производится путем сравнения событий разорения (73) при соответствующих значениях параметров.
Агрегированный процесс риска
В данном разделе рассматривается процесс риска в дискретном времени, который оказывается тесно связанным с классическим процессом риска посредством операции агрегирования.
Операция агрегирования
Рассмотрим классический процесс риска (62), зафиксируем число 5 0 и разобьем положительную полуось R+ на интервалы Д* = [(i 1)5, i5) длины 5. Далее, сгруппируем все премиальные поступления и страховые убытки, произошедшие в этих интервалах времени. Тогда размер премиальных поступлений за любой период Д* равен, очевидно,
c = c5. (77)
Размер Z* накопленных страховых убытков на интервале Д* вычисляется следующим образом. Для каждого i размер Z* является случайной величиной, причем для i = 1 ее распределение уже вычислено в (68), следует лишь подставить значение длины интервала времени t = 5:
~ (\5)k _
Fz(v) = P{Zi v} = P{^[o;5) v} = ? e-xs Ц-)-Ff (v), v 0. (78)
k=0 k
Далее, ввиду стационарности потока страховых событий и независимости и одинаковой распределенности убытков классического процесса риска {Zk, к = 1, 2,...} размеры убытков {Z*, i = 1, 2,...} на интервалах Д* также являются независимыми одинаково распределенными случайными величинами с функцией распределения
(78).
Таким образом, значение классического процесса риска (62) в момент времени n5 при целых n равно
X(n) = x + cn ^2 Z%, n = 0,1, 2,.... (79)
i=\
Процесс (79) называется агрегированным процессом риска и является аппроксимирующей моделью для классического процесса риска при 5 ^ 0.
Ясно, что процесс (79) зависит от трех параметров (x, c, Fz) вычисляемых по параметрам исходного классического процесса (x, c, Л, Ff) и значению параметра агрегирования 5 по формулам (77), (78). Будем обозначать
X = X(x, c, Fz) (80)
агрегированный процесс риска, определяемый параметрами x, c, Fz а X - совокупность всех аргегпрованных процессов риска при всевозможных допустимых значениях параметров:
(81)
X = {X(x, c, Fz); x 0, c 0, Fz(0) = 0}.
Разорение
Обозначимсобытие разорения агрегированного процесса,момент разорения ивероятность разорения и выживания процесса, соответственно. Рассматривая отношение включения событий разорения (выживания) при различных значениях параметров процесса риска, нетрудно заключить, что вероятность выживания является неубывающей функцией начального капитала x, интенсивности премиального потока c и функции распределения убытков FZ, если на множестве функций распределения рассмотреть естественный частичный порядок
Fi F2 Fi(v) F2(v), v 0.
8.3 Случайное блуждание
Агрегированный процесс риска заменой переменных Y = c Zj, i = 1, 2,... можно представить в виде случайного блуждания [3]
П
X (n) = x + ^ Y (86)
i=1
с независимыми одинаково распределенными шагами1^, i = 1, 2,.... Отсюда, используя теорему п. 2.XII.2 из [3], выводим
Теорема 8.1 Агрегированный процесс риска (79) может принадлежатъ к одному и только к одному из трех типов в зависимости от соотношения между его параметрами:
1. c = EZ1; осциллирующий тип; процессы этого типа с вероятностью единица достигают любого наперед заданного уровня; точнее:
| Р{ inf X(n) = to} = 1 n? N |
(87) |
| и | |
| P{sup X (n) = to} = 1. | (88) |
| neN | |
| 2. c EZ^- разоряющийся тип; для процессов этого вероятностью 1 существует конечный максимум | типа выполняется (87) и с |
| M(x) = max X (n). neN |
(89) |
m(x) = min X (n). (90)
n6 N
Из приведенной теоремы ясно, что процессы первого типа ввиду (87) разоряются с вероятностью 1. Поскольку для процессов второго типа (87) также имеет место, они тоже являются разоряющимися с вероятностью 1. Для процессов третьего типа вероятность выживания есть
S(x) = P{m(x) 0} (91)
и может быть, вообще говоря, положительной. Более того, поскольку для произвольных x, y 0 имеет место m(x) = m(y) + (x y), для Gx - функции распределения m(x) - получаем
Gx(v) = P{m(x) v} = P{m(y) y v x} = Gy(v + y x) ^ 0
при x ^ то и фиксированных v, y, поэтому справедлива
Теорема 8.2 Пусть выполнено
(92)
(93)
c EZi.
Тогда
lim S(x) = 1.
Замечание 8.1 Условие (92) обычно называется условием, положительности рисковой надбавки c EZ1.
8.4 Уравнение для вероятности разорения
Используя формулу полной вероятности, нетрудно вывести интегральное уравнение для вероятности выживания процесса (79), как функции начального капитала. Зафиксируем параметры c, Fz, и определим событие выживание при условии, что начальный капитал равен x: а для интервала I - аналогичное событие выживание при условии, что начальный капитал принадлежал интервалу I: Ясно, что S(x) С S1 = {Z1 x + c}. Для произвольного целого m 0 разобьем отрезок [0, x + c) на m равных частей Ik длины y = (x + c)/m: Ik = [(k 1)д, кд), k = 1,..., m. Тогда по формуле полной вероятности имеем P{S (x)} = Y) S (x + c Ik )P{Zi Е Ik} = X) S (x + c Ik) (Fz (кд) Fz (к 1)д)). Правая часть последнего выражения представляет собой, очевидно, интегральную сумму для интеграла
/ S(x + c v) dFZ(v), J 0
и сходится к нему при m ^ то, а левая часть не зависит от m и равна S(x), так что
S(x) = / S(x + c v) dFZ(v). (96)
J 0
Это интегральное уравнение позволяет изучать многие свойства агрегированного процесса риска, выраженные в терминах вероятностей выживания или разорения.
Пример: простейший процесс риска
Положим c = 1, p Е (0, 1) и
( 0, v 0,
Fz(v) = p, 0 v 2, ll, v 2.
Нетрудно заметить, что при этом агрегированный процесс риска (79) превращается в простейший процесс риска (53), а уравнение для вероятности выживания приобретает вид
S (x) = S (x + 1) + (1 p)S (x 1),
что согласуется с (56).
Время жизни процессов риска
В параграфах 6-8 получены уравнения для вероятности разорения. Представляет интерес также и момент разорения.
В настоящем параграфе изучим среднее время жизни процессов риска.
Простейший процесс риска
Рассмотрим простейший процесс риска (53):
П
Xn = z + ? Zi, (97)
i= 1
где Zi, i = 1, 2,... - независимые одинаково распределенные случайные величины, принимающие значение 1 с вероятностью p и значение -1с вероятностью q =1 p, и обозначим tz момент окончания игры (разорения одного из игроков) при условии, что первый игрок имел начальный капитал z, а второй: a z:
tz = min{n 0 : Xn = 0 или Xn = a}. (98)
Как мы помним, на бесконечном горизонте времени разорение одного из игроков наступает с вероятностью 1, так что tz является собственной случайной величиной.
Обозначим m(z) = Etz ее среднее значение, и поставим задачу вычисления этой величины. Из формулы полной вероятности для математических ожиданий нетрудно получить уравнение для m(z) вида m(z) = 1 + pm(z + 1) + qm(z 1), откуда
Характеристическое уравнение для (99) имеет корни Аі = 1; Л2 = q/p = y. Рассмотрим сначала случай некратных корней у = 1 (p =1/2). Тогда частное решение неоднородного уравнения (99) имеет вид m0(z) = p(yz_1), так что его общее решение запишется в форме
z
m(z) = C1 + C2yZ + 7-77 - (100)
p(y 1)
Произвольные постоянные C1, C2 можно определить из краевых условий m(0) = m(a) = 0, которые дают
C2 = Cl
a
p(1 У)(1 Уa) ’
так что решение уравнения (99) есть
1
p(1 У)
1 yz
a--z
1 ya
(101)
m(z)
В случае кратных корней y = 1 (p =1/2) общее решение уравнения (99) имеет вид C1 + C2z + m0(z), где частное решение m0(z) легко находится в классе квадратичных функций и равно m0(z) = z2, так что искомое решение после определения постоянных из краевых условий записывается в форме
m(z) = z(a z). (102)
Интересно отметить, что решение (102) может быть получено предельным переходом из (101) с помощью двукратного применения правила Лопиталя.
Рассмотрим, наконец, игру с бесконечно богатым противником (a ^ то).При y 1 из (101), видно, что среднее время жизни неограниченно увеличивается. Этот эффект связан с тем, что при y 1 имеем p 1/2, и даже в игре с бесконечно богатым противником первый игрок разоряется с вероятностью, меньшей 1, то есть с положительной вероятностью время жизни процесса бесконечно.
При y = 1 го (102) также вытекает m(z) ^ то при a ^ то. Отсюда вытекает, что хотя в этом случае (p = 1/2) первый игрок и разоряется с вероятностью 1 в игре с бесконечно богатым противником, однако среднее время ожидания разорения бесконечно велико.
При y 1 (p 1/2) из (101) следует, что ожидаемое время до разорения в игре с бесконечно богатым противником есть
m(z) = ( z lV p(y 1)
и линейно возрастает вместе с z.
Игра в кошки мышки
Рассмотрим еще одну задачу на вычисление ожидаемой продолжительности процесса. Пусть в двух комнатах обитают кот и мышь, и перемещаются из комнаты в комнату независимо друг от друга в соответствии с матрицами переходных вероятностей
0.3 0.7 0.6 0.4
0.2 0.8
0.8 0.2
(103)
где элемент с- матрицы С есть вероятность перемещения кота из комнаты i в комнату j; i,j = 1, 2, а элементы матрицы M описывают аналогичные вероятности для мыши. Если в какой-либо из моментов времени кот и мышь оказываются в одной комнате, то кот съедает мышь.
Определить среднее время жизни мыши при различных начальных расположениях.
Опишем задачу в виде марковской цепи с четырьмя состояниями E = (i, j); i, j = 1, 2: кот находится в комнате i, а мышь - в комнате j, и обозначим t- среднее время оставшейся жизни мыши в состоянии (i, j). Обозначим En состояние цепи в момент времени п, тогда, очевидно, tn = t22 = 0, а для t12, t2i по формуле полной вероятности запишем уравнения
t12 = 1 + t12P{E1 = (1, 2) |E0 = (1, 2)} + t21P{E1 = (2, 1) |E0 = (1, 2)} t21 = 1 + t12P{E1 = (1, 2)|E0 = (2, 1)} + t21P{E1 = (2, 1)|E0 = (2, 1)}.
Вычисляя переходные вероятности с учетом условия независимости, получаем систему уравнений
0.94t12 - 0.56t21 = 1,
-0.48t12 + 0.92t21 = 1,
которая имеет решение t12 = 2.483, t21 = 2.383.
