Главная » Просмотр файлов » Наброски лекций. 2017

Наброски лекций. 2017 (1134111), страница 9

Файл №1134111 Наброски лекций. 2017 (Наброски лекций. 2017) 9 страницаНаброски лекций. 2017 (1134111) страница 92019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 9)

. . ,а через произвольные промежутки ∆t > 0. Доопределим матричнозначную функцию π( · ) в точке t = 0, наложив совершенно очевидное условие π(0) = I, где I –единичная (s × s)-матрица, т. е.(1, i = j,πij (0) = δij =(4.8)0, i 6= j.Это условие говорит о том, что за нулевое время система не сможет совершитьникакого перехода и останется в своём состоянии.В дальнейшем нас будут интересовать аналитические свойства матричнозначнойфункции π( · ) : R+ 7→ Ms : непрерывность, дифференцируемость и пр. Условимся,что матричные равенства, связанные с предельными переходами, мы будем понимать как поэлементные, напримерlim π(h) = Ah→+0⇐⇒lim πij (h) = Aijh→+0а равенстваZdπ= D,π̇(t) ≡dtZi, j = 1, .

. . , s,t2π(t) dt = Jt1эквивалентныdπijπ̇ij (t) ≡= Dij ,dtдля всехt2πij (t) dt = Vijдля всех i, j = 1, . . . , s.t1Здесь A, D, J ∈ Mat(s).4.1. Система уравнений Колмогорова для матриц перехода. В нашихматематических выкладках очень важную роль будет играть условие s < ∞. Онопозволит нам менять местами предельные переходы и любые суммы или произведения по состояниям:limh→+0sXi=1(·) =sXi=1lim ( · ),limh→+0h→+040sYi=1(·) =sYi=1lim ( · ).h→+0Для бесконечного (счётного) числа состояний эти равенства требуют обоснования,но многие утверждения тем не менее останутся верными, возможно, при некоторыхдополнительных условиях.Докажем несколько лемм.Лемма 4.2. При любом t > 0 матрица π(t) невырождена.Доказательство. Заметим, что в силу (4.8) определитель матрицы π(h) какнекоторая вполне конкретная (конечная) сумма произведений её элементов стремится при h → +0 к определителю матрицы I, т.

е. к единице. Поэтому найдётся такое h0 > 0, что det π(h) 6= 0 при h < h0 . Это означает, что матрица π(h)невырождена при достаточно малых h. Далее, для произвольного фиксированногоt > 0 всегда можно найти такое натуральное n, что t/n < h0 . Тогда по уравнениюЧепмена–Колмогороваndet π(h) = det π n (t/n) = det π(t/n) 6= 0,t/n < h0 ,а это и означает, что матрица det π(h) невырождена при любом h > 0.Теперь покажем непрерывность матрицы перехода.

Для t = 0 непрерывностьбудем понимать, конечно же, как непрерывность справа и постулируем эту непрерывность без доказательства: потребуемπij (h) → δijприh → +0(4.9)или, в матричной записи,π(t) → π(0) = Iприt → +0.(4.10)Это условие согласовано с (4.8) и означает, что, попав в любое состояние xi , случайный процесс с вероятностью единица остаётся в этом состоянии некоторое положительное время. Точнее, если случайная величина τi равна времени пребыванияв i-м состоянии, то P (τi > 0) = 1.Лемма 4.3. Матричнозначная функция π( · ) : R+ 7→ Ms непрерывна при всехt > 0, т. е.

для любого фиксированного t > 0lim π(t + h) − π(t)) = lim π(t) − π(t − h)) = 0,h→+0h→+0где в правой части равенства 0 означает матрицу, все элементы которой равнынулю2) .Доказательство. Непрерывность справа непосредственно вытекает из уравнения Чепмена–Колмогорова (4.7) и условия (4.10):πij (t + h) − πij (t) =sXk=1πik (t)πkj (h) − πij (t) =sXk=1πik (t) πkj (h) − δkj → 0,h → +0,при любых i, j = 1, . . .

, s для любого фиксированного t > 0.2) Мы не вводим специальных обозначений, отличающих 0 как нулевую матрицу от числа 0;предполагается, что природа равенства должна сама подсказать, о каком нуле идёт речь.41Непрерывность слева также следует из уравнения Чепмена–Колмогорова, но потребуются чуть более сложные рассуждения. Для t > 0 и 0 < h < t имеемπ(t) − π(t − h) = π(t − h) π(h) − I ,отсюда в силу 0 6 πik (t − h) 6 1 получаем оценку sX|πij (t) − πij (t − h)| = πik (t − h) πkj (h) − δij 6k=16sXk=1πik (t − h)|πkj (h) − δkj | 6sXk=1|πkj (h) − δkj | → 0,h → +0,при любых i, j = 1, . . . , s для любого фиксированного t > 0.

Лемма доказана.Теперь исследуем матрицу π̇(t) с элементамиπ̇ij (t) =π(t + h) − π(t)dπij= lim,h→0dthi, j = 1, . . . , s.Лемма 4.4. Зададим матричнозначную функциюZ t2π(u) du,0 < t1 < t2 .J(t1 , t2 ) =t1Тогда при достаточно малой разности t2 − t1 существует обратная матрицаR(t1 , t2 ) = J −1 (t1 , t2 ).Доказательство. Пользуясь непрерывностью π( · ) как матричнозначной функции на R+ , запишем формулу среднего для каждого элемента матрицы: для любых0 < t1 < t2 при каждых i, j ∈ {1, . .

. , s}Z t21πij (u) du = πij (t1 + hij ),0 6 hij 6 t2 − t1 .t2 − t1 t1Опять же в силу непрерывности π( · ) для любого ε > 0 найдётся δij > 0 такая, что|πij (t1 + hij ) − πij (t1 )| < ε, если hij < δij . Выберем минимальную по i, j = 1, . . . , sвеличину δij и обозначим её как δ. Тогда для любого t1 > 0 c учётом введённого вусловии леммы обозначения Jij (t1 , t2 ) = |πij (t1 + hij ) − πij (t1 )| < ε, если 0 < t2 − t1 < δ. (4.11)−π(t)ij1 t2 − t1Итак, для любого ε > 0 найдётся δ = δ(t1 , ε) такое, что det J(t1 , t2 ) < ε,−detπ(t)0 < t2 − t1 < δ.ij1 t2 − t1Поскольку det π(t1 ) 6= 0 при любом t1 > 0, мы можем выбрать какое-либо положительное ε < | det π(t1 )| и зафиксировать его.

Теперь δ определяется только значением t1 и1det J(t1 , t2 ) > det πij (t1 ) − ε 6= 0,t2 − t1если0 < t2 − t1 < δ.Следовательно, при 0 < t2 − t1 < δ существует R(t1 , t2 ) = J −1 (t1 , t2 ).42Лемма 4.5. Матричнозначная функция π(·) имеет правую производную в точке t = 0, т. е. существует матрица Λ размера s × s такая, чтоπ(h) − π(0)= Λ.h→+0h(4.12)limДоказательство. Воспользуемся уравнением (4.7) Чепмена–Колмогорова и запишем цепочку матричных равенств для h > 0 и любых 0 < t1 < t2π(h) − IZt2π(u) du =t1Z=t2t1Zt2t1π(h)π(u) du −π(u + h) du −Zt2Zπ(u) du =t1π(u) du =t1t2Zt2 +ht1 +hπ(u) du −Zt2π(u) du.t1Теперь будем считать, что 0 < h < t2 − t1 , тогда t1 < t1 + h < t2 < t2 + h иZt2 +hπ(u) du −t1 +hZt2π(u) du =t1Zt2 +hπ(u) du −t2Zt1 +hπ(u) du.t1В обозначениях леммы 4.4 мы можем записать результат нашей цепочки преобразований какZ t2 +hZ t1 +hπ(u) du −π(u) du,(4.13)π(h) − I J(t1 , t2 ) =t2t1и это равенство верно для любых 0 < t1 < t2 и 0 < h < t2 − t1 .Пусть теперь t2 − t1 < δ, где δ = δ(t1 ) взято из леммы 4.4.

Тогдаπ(h) − I=h Z t2 +hZ1 t1 +h1π(u) du −π(u) du R(t1 , t2 ).h t2h t1При h → +0 с учетом формул (4.11) для средних значений интегралов в правойчасти равенства существуетlimh→+0Zt2 +ht2π(u) du −Следовательно,limh→+0Zt1 +ht1π(u) du = π(t2 ) − π(t1 ) R(t1 , t2 ).π(h) − I= π(t2 ) − π(t1 ) R(t1 , t2 ),h(4.14)причём это равенство верно при любых достаточно близких друг другу положительных временах t1 < t2 . Выберем какие-либо удовлетворяющие всем необходимымусловиям t1 и t2 , зафиксируем их и для этих фиксированных времён обозначимматрицу в правой части равенства (4.14) через Λ.

Тогда, учитывая, что I = π(0),получаемπ(h) − π(0)= Λ.limh→+0hЛемма доказана.43Теорема 4.1. Матрица π( · ) непрерывно дифференцируема при t > 0 и удовлетворяет уравнениям Колмогороваπ̇(t) = Λπ(t),(4.15)π̇(t) = π(t)Λ.(4.16)Доказательство. Воспользуемся уравнением (4.7). Для произвольного t > 0имеемπ(h) − π(0)π(h) − π(0)π(t + h) − π(t)=π(t) = π(t),hhhоткудаπ(t + h) − π(t)lim= Λπ(t) = π(t)Λ,(4.17)h→+0hт. е. матрица π(t) имеет правую производную при всех t > 0 и для правой производной верны уравнения (4.15) и (4.16).Теперь найдём левую производную.

Зафиксируем t > 0 и выберем 0 < h < t,тогдаπ(h) − π(0)π(h) − π(0)π(t) − π(t − h)=π(t − h) = π(t − h),hhhперейдём к пределу при h → +0, учитывая непрерывность π(t − h) → π(t), получимπ(t) − π(t − h)= Λπ(t) = π(t)Λ,h→+0h(4.18)limи эти равенства вместе с (4.17) влекут уравнения (4.15), (4.16).Непрерывность матричнозначной функции π̇(t), t > 0, непосредственно следуетиз непрерывности матричнозначной функции π( · ) и любого из уравнений (4.15),(4.16). Теорема доказана.Замечание 4.6.

В случае счётного числа состояний система уравнений Колмогорова также верна, но следует наложить на матрицы π(t) некоторые дополнительные условия.4.2. Система уравнений Колмогорова для вероятностей состояний. Положим pi (t) = P (ξ(t) = xi ), i = 1, . .

. , s, и выведем из уравнения (4.16) системудифференциальных уравнений для функций pi ( · ). Имеемdefpj (t) == P (ξ(t) = xj ) =sXi=1P (ξ(t) = xj | ξ(0) = xi )P (ξ(0) = xi ) =где ai – начальная вероятность. Отсюда в силу уравнения (4.16)ṗj (t) ==sXi=1sXk=1π̇ij (t)ai =ΛkjXsss XXi=144πik (t)Λkj ai =k=1πik (t)aii=1=sXk=1pk (t)Λkj ,sXi=1πij (t)ai ,Итак, мы имеемṗj (t) =sXpk (t)Λkj ,(4.19)j = 1, . . .

, s.k=1В ряде случаев уравнение (4.19) удобно записать в другом виде. Заметим, чтодля любого j = 1, . . . , s и любого h > 0sXk=1πjk (h) − δjk =sXk=1поскольку выполнено условие стохастичностиsXk=1sΛjkπjk (h) − 1 = 0,Psk=1πjk (h) = 1.

Отсюда1X= limπjk (h) − δjk = 0,h→+0 hi=1Λjj = −XΛjk .k : k6=jТаким образом, получаемṗj (t) =sXk : k6=jpk (t)Λkj − pj (t)XΛjk ,j = 1, . . . , s.(4.20)k : k6=jЭтот вид уравнения полностью отвечает физическому смыслу переходов в марковском процессе. Перепишем (4.20) как следующее приближённое равенство прималом h > 0, учтя, что Λ = π̇(0):pj (t + h) − pj (t) ≈=sXk : k6=jsXk : k6=jXpk (t)(πkj (h) − δkj −pj (t)(πjk (h) − δjk =k : k6=jpk (t)πkj (h) −Xpj (t)πjk (h).(4.21)k : k6=jРассмотрим систему из большого числа N частиц, в которой каждая из частицможет находиться в одном из состояний x1 , .

. . , xs , а динамика распределений по состояниям осуществляется по законам марковского процесса. Заменим вероятностьна долю частиц в системе, находящихся в момент времени t в состоянии xj , т. е. положим pj (t) ≈ Nj (t)/N . Тогда (4.21) описывает, как изменяется количество частиц,находящихся в момент t в состоянии xj , за малое время h (уравнение баланса):Nj (t + h) − Nj (t) ≈sXk : k6=jNk (t)Pk→j (h) −XNj (t)Pj→k (h),(4.22)k : k6=jгде Pk→j (h) – вероятность перехода из xk в xj за время h. При этом первая суммав правой части (4.22) даёт число частиц, пришедших в состояние xj из всех другихсостояний, а вторая – число частиц, ушедших из состояния xj в какое-либо другоесостояние. Мы полагаем, что за малое время происходит не более одного перехода,что согласуется с разложением π(h) = I + Λh + o(h).Система уравнений (4.19) или (4.20) дополняется начальным условием pj (0) = aj ,Psj = 1, .

Характеристики

Тип файла
PDF-файл
Размер
468,8 Kb
Тип материала
Высшее учебное заведение

Список файлов лекций

Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6553
Авторов
на СтудИзбе
299
Средний доход
с одного платного файла
Обучение Подробнее