автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Автономная немарковская система массового обслуживания и ее применение в задачах демографии

На правах рукописи

Носова Мария Геннадьевна

АВТОНОМНАЯ НЕМАРКОВСКАЯ СИСТЕМА МАССОВОГО ОБСЛУЖИВАНИЯ И ЕЁ ПРИМЕНЕНИЕ В ЗАДАЧАХ ДЕМОГРАФИИ

05.13.18 - математическое моделирование, численные методы и комплексы программ

Автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

0034Э2383

Томск-2010

003492383

Работа выполнена на кафедре теории вероятностей и математической статистики факультета прикладной математики и кибернетики Томского государственного университета.

Научный руководитель: доктор технических наук,

профессор Назаров Анатолий Андреевич.

Официальные оппоненты: доктор физико-математических наук,

профессор Змеев Олег Алексеевич.

кандидат физико-математических наук, доцент Гарайшииа Ирина Рашитовна.

Ведущая организация: Сибирский федеральный университет, г.Красноярск.

Защита состоится 18 марта 2010 года в 10^ на заседании диссертационного совета Д 212.267.08 при Томском государственном университете по адресу: 634050, г.Томск, пр.Ленина, 36, корп. 2, ауд. 102.

Отзывы на автореферат (в двух экземплярах), заверенные гербовой печатью организации, просим направлять по адресу: 634050, г.Томск, пр.Ленина, 36, ученому секретарю ТГУ Буровой Н.Ю.

С диссертацией можно ознакомиться в Научной библиотеке Томского государственного университета по адресу: 634050, г.Томск, пр.Ленина, 34а.

Автореферат разослан 12 февраля 2010 года.

Ученый секретарь

диссертационного совета Д 212.267.08 доктор технических наук, профессор

А.В. Скворцов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

При анализе и прогнозировании демографических процессов наиболее эффективным инструментом является использование методов математического моделирования.

К первым математическим моделям демографических процессов относятся детерминированные модели роста человечества, а именно модели линейного и экспоненциального роста, описание которых встречается в работах И.Г. Венецкого, Д. Кендалла, X. Касвелла и других авторов.

Наиболее известная детерминированная модель — модель стабильного населения (непрерывные и дискретные модели). Ее создание связывают с такими именами как Л. Эйлер, Г. Кнапп, В. Лсксис, Дж. Лотка, В. Борткевич, П. Лссли. Большой вклад в разработку методов практического применения стабильного населения внесли советские демографы С.А. Новосельский, В.В. Паевский, А.Я. Боярский, И.Г. Венецкий и др.

Современный этап развития теории стабильного населения связывают с обобщением основных выводов на случай демографических процессов с переменными интенсивностями. В этом направлении работали А. Коул и А. Лопес. А также Р. Лии, Д. Ахлбург, 3. Сайкес, М. Алхо, Б. Спенсер, Дж. Поллард, Дж. Кохен, Н. Кейфиц и другие.

В конце XIX века учеными была предпринята попытка использовать в качестве модели, описывающей рост человечества, логистическую кривую. Этой модели придерживались Р. Пирль и Л. Рид.

Еще одной детерминированной моделью воспроизводства населения является модель гиперболического изменения численности населения, предложенная в 1960 г. X. Ферстером, П. Мором и Л. Амио-том. В настоящее время в этом направлении работает С.П. Капица.

Значительный вклад в математическую демографию сделан О.В. Староверовым, который рассматривал демографические процессы в виде марковских моделей в форме цепей Маркова. Староверовым О.В. предложена стохастическая модель развития населения с дискретным временем, учитывающая случайность в рождаемости и смертности.

Метод компонент или метод передвижки разработан П.К. Уэл-птоном. В России перспективными исчислениями численности насе-

ления методом возрастных передвижек занимались С.Г. Струмилин, А.Я. Боярский, П.П. Шушерин, М.С. Бедный, а в последнее время Государственный комитет Российской Федерации по статистике, Центр демографии и экологии человека, Отдел населения ООН.

Приведенный анализ множества моделей, показал, что в моделировании демографических процессов наиболее распространены детерминированные модели (дискретные и непрерывные) и стохастические дискретные. Преимущество стохастических моделей перед детерминированными моделями заключается в том, что они учитывают отклонение частот демографических событий от их вероятностей.

Однако демографические процессы протекают в непрерывном времени и являются стохастическими. Методы исследования таких процессов в математической демографии недостаточно развиты. Именно поэтому актуальной является задача существенного расширения математических моделей процесса изменения демографической ситуации, а также развития методов их исследования.

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

В рамках указанной цели были поставлены следующие задачи:

1. Разработка математического метода исследования автономных стохастических систем с непрерывным и дискретным временем.

2. Разработка стохастических моделей демографии в виде упрощенных автономных систем массового обслуживания с РН-распределением времени обслуживания и развитие методов моментов и асимптотического анализа для их исследования.

3. Разработка новой математической модели процесса изменения демографической ситуации в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания.

4. Разработка метода виртуальных фаз и модификация метода асимптотического анализа для исследования автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания.

5. Применение разработанных математических моделей к исследованию процесса изменения демографической ситуации.

Научная новизна состоит в следующем:

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

2. Разработан метод виртуальных фаз для исследования автономной немарковской системы массового обслуживания с неограниченным числом приборов, позволяющий находить многомерное распределение вероятностей численностей обслуживаемых заявок различных возрастов и основные характеристики, определяющие эти распределения.

3. Для исследования упрощенных моделей (трехфазной и пяти-фазной автономных систем, автономной системы с двумя типами заявок, автономной системы с РН-распределением времени обслуживания) предложены модификации метода моментов и метода асимптотического анализа. Найдены распределения вероятностей числа обслуживаемых заявок в таких системах и их основные характеристики.

4. На основе метода прямой передвижки возрастных групп предложен метод обратной передвижки, позволяющий выполнять интерполяцию демографических данных между датами переписей населения, а также восстанавливать функцию интенсивности процесса рождаемости. Метод прямой и обратной передвижки возрастных групп применены для оценки величины людских потерь РФ в годы ВОВ. Значение оценки составляет 17 млн. человек.

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

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

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

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

Апробация работы. Основные положения работы и отдельные ее результаты докладывались и обсуждались:

1. III Всероссийская научно-практическая конференция «Информационные технологии и математическое моделирование (ИТММ-2004)», г.Анжеро-Судженск, 2004 г.

2. Международная научная конференция «Математические методы повышения эффективности функционирования телекоммуникационных сетей», г.Минск, 2005 г.

3. IX Всероссийская научно-практическая конференция «Научное творчество молодежи», г.Анжеро-Судженск, 2005 г.

4. V Международная научно-практическая конференция «Информационные технологии и математическое моделирование», г.Анжеро-Судженск, 2006 г.

5. XI Всероссийская научно-практическая конференция «Научное творчество молодежи», г.Анжеро-Судженск, 2007 г.

6. Международная научная конференция «Математические методы повышения эффективности функционирования телекоммуникационных сетей», г.Минск, 2007 г.

7. VII Российская конференция с международным участием «Новые информационные технологии в исследовании сложных структур», г.Томск, 2008 г.

Работа выполнена при поддержке АВЦП «Развитие научного потенциала высшей школы (2009-2010 гг.)» Федерального агентства по образованию РФ по проекту «Разработка методов исследования немарковских систем массового обслуживания и их применения к сложным экономическим системам и компьютерным сетям связи».

Публикации. По результатам выполненных исследований автором опубликованы 13 печатных работ, в том числе б статей, из них 3 в изданиях, рекомендованных списком ВАК.

Личное участие автора в получении результатов, изложенных в диссертации.

Постановка изложенных в диссертации задач была сделана научным руководителем аспиранта, доктором технических наук, профессором, Назаровым A.A., который указал основные направления исследования и принимал участие в обсуждении результатов. Доказательство и обоснование полученных в диссертации результатов, математические выкладки и численные расчеты выполнены лично автором работы. В совместных публикациях научному руководителю Назарову A.A. принадлежат постановки задач и указания основных направлений исследований, а основные результаты, выкладки и численные расчеты выполнены диссертантом.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения, списка литературы из 158 наименований и приложений. Общий объем составляет 204 страницы, в том числе основной текст 160 страница.

СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована важность и актуальность работы, сформулирована цель и задачи диссертационного исследования, показана его научная новизна, выявлена теоретическая и практическая ценность полученных результатов, кратко излагается содержание диссертационной работы.

В главе 1 предлагаются решения некоторых задач математической демографии с помощью метода математического моделирования. В разделе 1.1 дано математическое обоснование метода прямой передвижки возрастных групп. Метод прямой передвижки применяется для определения оценок значений численности N(x + \,t + ч) группы лиц возраста х + т в году t + x при условии, что известна численность N(x,t), т -шагпрогнозирования.

На основе метода прямой передвижки возрастных групп предложен метод обратной передвижки, который позволяет находить оценку N(x-x,t-x) значения численности демографической группы лиц возраста х-х ъ году t - т.

В разделе 1.2 найдены условные математическое ожидание и дисперсия величины N(x-x,t-x) при условии, что выполняется равенство N(x, t) = п.

В разделе 1.3 на основе полученных формул применены методы прямой и обратной передвижки возрастов к решению задачи оценки величины людских потерь РФ в годы ВОВ. Показано, что суммарные значения оценки людских потерь при методах прямой и обратной передвижки возрастных групп совпадают и составляют S = 17 млн. чел.

В разделе 1.4 представлено математическое обоснование нецелесообразности аппроксимации в математической демографии процесса рождаемости потоком Пуассона при долгосрочном прогнозировании. Для этого на первом этапе рассматривается наиболее простая для исследования математическая модель случайного потока, управляемого цепью Маркова. Такая модель определяется как поток заявок в автономной системе массового обслуживания с неограниченным числом приборов и средним значением m(t) числа занятых приборов. Для системы определено распределение вероятностей P{i,i) и составлена система дифференциальных уравнений Колмогорова, найден аналитический вид характеристической функции H(u,t)

лЯ

ъ ¿V ь -1)ехр{(& -т

Ъ -Vм - •1)ехр{(Ь- И)'}

Далее исследуется открытая система массового обслуживания с неограниченным числом приборов, на вход которой поступает нестационарный пуассоновский поток с интенсивностью , выбираемой из условия совпадения средних значений числа приборов, занятых в автономной системе и в открытой системе с пуассоновским входящим потоком. Для системы записано и решено дифференциальное уравнение для характеристический функции числа занятых приборов С(и,1). Характеристичная функция (?(«,/) имеет вид

С(и,0 =4-г"'" +е->ие-*}Ы ехр|> -1)}

Через обратное преобразование Фурье найдены распределения вероятностей /'(г,/) и РА{1,1), определяемые характеристическими функциями С/(г<,0 и Н(и, ?).

Получены значения расстояния р(й,г) между распределениями вероятностей Р(/,/) и РА{для различных значений параметров Ь и

t, показано, что расстояние p(b,í) возрастает с ростом значений параметров but. Поскольку, характеристики таких систем массового обслуживания существенно различаются, следовательно, делается вывод о нецелесообразности аппроксимации процесса рождаемости потоком Пуассона при долгосрочном прогнозировании.

В главе 2 предлагается математическая модель женского населения в виде автономной системы массового обслуживания с РН-распределением времени обслуживания и неограниченным числом приборов. Особенностью такой системы является отсутствие внешнего источника заявок.

Продолжительность обслуживания т каждой заявки в такой системе складывается из продолжительностей случайного числа фаз т(. Величины т, являются независимыми экспоненциально распределенными случайными величинами с параметрами ц,, i = 1, v.

В терминах демографии обслуживаемая заявка интерпретируется как женщина, время обслуживания заявки - продолжительность жизни этой женщины, величина bj(t) - интенсивность рождения девочек у женщины í-й возрастной группы в году t. Входящим потоком заявок является процесс рождения девочек.

В этой главе рассматриваются некоторые упрощенные варианты автономной системы с РН-распределением времени обслуживания и методы их исследования, а именно, трехфазная система (раздел 2.2), пятифазная система (раздел 2.3) и система с двумя типами заявок (раздел 2.4).

Состояние упрощенной автономной системы в момент времени t определяется многомерным вектором n(t). Для указанных моделей определены многомерные распределения вероятностей P(n,t), где п = {п\,П2,щ,...\, для которых составлены системы дифференциальных уравнений Колмогорова и записаны дифференциальные уравнения для многомерных характеристических функций H(u,t) числа заявок, обслуживаемых в системах, здесь и = {щ,и2,щ,...}.

Для нахождения распределения вероятностей, а также основных его характеристик применены методы моментов и асимптотического анализа, модифицированные к рассматриваемой задаче.

В данном случаем при методе асимптотического анализа уравнение для характеристической функции Я (и, О рассматривалось в асимптотическом условии большой численности групп, пропорциональных некоторой бесконечно большой величине N, что вполне

обосновано в демографических исследованиях, где численности стандартных возрастных групп измеряются миллионами.

В первой асимптотике были предложены замены u = zw, H(u,i)-Fi(w,t,е), где 8 = 1 IN, £ — малый положительный параметр.

Показано, что

lim (w,r,6) = ^[(w,/) = exp{./'/«(Оw},

E—>0

где аналитическое среднее Nm(t) является решением системы обыкновенных дифференциальных уравнений первого порядка.

Для нахождения асимптотики второго порядка в задаче для F\(w,t,z) выполнены замены

Я (и, 0 = cxp{jNm(t)u}H2 (u,t), тогда для функций H2(u,t) получены задачи, в которых, обозначая 2

е -MN, заменяем

и-zw, H2(u,t) - Fjiw,!,!:).

Показано, что

lim FjiwJ,?-) - Fi(\v,t) - exp-j - —wäiv7 [, E—>0 12 J

где элементы матрицы ковариаций NR-[NRtj] являются решением

системы обыкновенных дифференциальных уравнений первого порядка.

Определено, что многомерная величина n{t) имеет асимптотически нормальное распределение с вектором математических ожиданий Nm(t) и матрицей ковариаций NR = [NRy ]. Величины i и j определяются типом системы массового обслуживания.

Отметим, что асимптотические средние значения Nm,-(t) числа обслуживаемых заявок в рассматриваемых системах в момент времени f, найденные в предельном асимптотическом условии N <я, совпадают с допредельными, полученными методом моментов. Заметим, что такое совпадение справедливо и для вторых моментов в предельной и допредельной моделях.

Предложенная упрощенная стохастическая модель в виде автономной системы массового обслуживания с PH-распределением времени обслуживания и частные ее виды могут применяться для моделирования широкого класса технических систем и социально-экономических явлений. В частности в диссертации предложено применение для прогнозирования процесса изменения демографической ситуации, как это продемонстрировано в разделах 2.2.3, 2.3.3, 2.4.3.

В главе 3 предложена стохастическая модель в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания.

Определим процесс функционирования данной системы массового обслуживания в терминах теории массового обслуживания.

Каждая заявка входящего потока в момент своего поступления занимает свободный прибор и находится на нем в течение всего времени обслуживания, продолжительность которого случайна. Продолжительности обслуживания различных требований стохастически независимы, имеют одинаковое распределение, определяемое функцией 5(х) = 1 -В(х), где В(х) - функция распределения времени обслуживания заявки. Завершив обслуживание, заявка покидает систему.

Каждая заявка возраста х в момент времени I с интенсивностью Ь(х, 0 генерирует новое требование. Новая заявка в момент своего появления занимает свободный прибор и начинает процесс своего обслуживания, генерируя требования нового поколения.

В терминах демографии обслуживаемая заявка интерпретируется как женщина,. время обслуживания заявки - продолжительность жизни этой женщины, Б(х) — функция дожития для женщин, возраст* заявки - возраст женщины в рассматриваемый момент времени, функция ¿(х,?) ~ интенсивность рождения девочек у женщины возраста х в году ? (функция фертильности).

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

В разделе 3.2 исследована общая система массового обслуживания с РН-распределением времени обслуживания. В этой модели время обслуживания аппроксимировано суммой случайного числа независимых одинаково экспоненциально распределенных случайных величин. Для распределения вероятностей -Р(и1>«2,...,/) числа обслуживаемых заявок записано дифференциальное уравнение Колмогорова.

Наряду с исследованием допредельной модели методом моментов, проведено ее асимптотическое исследование. Асимптотическое условие, а также замены асимптотик первого и второго порядков аналогичны заменам главы 2.

Показано, что в условии предельного уменьшения длительностей фаз, то есть при ц оо, рассматриваемое РН-распределение сходится к распределению, определяемому функцией S(x). Поэтому, сделав замены

щ = u(i / ц), щ =

\

H(u,t) = H(u,t,\i)-M< exp

/2 « . i=l

( ■ ^ i

m

Ц

и полагая, что при (.i —>оо, i —> со отношение —» х и существуют следующие пределы

lim н(//ц) = и(х),

i/\i->x

f ( 00

lim Н(и,(,ц.) = М<ехр\ j\u{x)^{x,t)dx > = F(u,t),

Ivo J

i!\i-wr

а также существует предел по распределению

lim |1и(//ц, /) = £(*, О,

г->оо,ц->да

в разделе 3.4 перешли к рассмотрению автономной немарковской системы массового обслуживания с произвольным временем обслуживания. Случайная функция \{x,t) назовем стохастической плотностью численности заявок возраста х в момент времени t.

Для характеристического функционала

(Г со V

F(u,t) = Mlcxp jju(x)E1(x,t)dx Ivo записано основное уравнение

а о Г ад J

Уравнение (1) будет решено методом асимптотического анализа в предположении, что случайная функция i;(x,i) принимает достаточно большие значения, пропорциональные некоторой бесконечно большой величине N, то есть 7V —> со. Для уравнения найдены асимптотики первого и второго порядков его решения.

Для нахождения асимптотики первого порядка в уравнении (1) выполняются замены

и = ш, F(u,t) = Fi(w,i,s). где £ = 1 / Л', а £ - малый положительный параметр.

Из (1) получим равенство

& ) й» Теорема 3.3. При £ 0 решение задачи (2) имеет вид

^ (и-, /) = ехр|у | м>(х)т(х, (3)

где является решением уравнения

(4)

з/ & ^о)

с краевым условием

00

т(0,1)= \т(х,ф(х,1)с1х. О

Для того чтобы найти асимптотику второго порядка в уравнении (1) выполним замену

Г(и,0 = Н2 («,Оехр| №\и{х)т{х,1)ск

о

а затем в полученном равенстве, обозначив б = 1/7У и введя замены и = ш, Н.1 (и, = Б) >

получим уравнение

а -1 '0 Ых) Г ; ; 5(х)

- ]т'(х) } А- ^(м-,*, Е)| | ^(х)^^Л + (5)

I о %

о ]

Теорема 3.4. Решение задачи (5) при £ 0 имеет вид Г 1 0000 1

(н-, 0 = ехр - -11 Му)Мг)Я(у, 2,1)ОусЬ к (6)

I ^ 00 ]

где слагаемые кросскорреляционной функции вида

ДОМ,/) = а2(г)5(у)5(2) + г(у,ф(г) + г(г,1Щу) + ф,1)ф,ф(у - г) определяются начальными и краевыми условиями

о2(х,0) = 0, г(ж,0) = 0, ст2(0) = 0, ст2(0,/) = 0, г(0,0 = 0 и являются решением системы дифференциальных уравнений в частных производных

dt дх < JS(x)

^ - r(x,f)is'(0) + + a2 (x,t)b(x,t), (7)

dt дх [ S(x)

2 оо

90 W = 25"(0)ст2(0 + 2\r(x,t)b(x,t)dx+ m(Q,t).

dt о

Полученные аналитические решения полностью определяют параметры гауссовского распределения, которому удовлетворяет распределение вероятностей стохастической плотности численности заявок E,(x,t) возраста х в году t.

В главе 4 предложено применение исследованной в главе 3 автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания к изучению процесса изменения демографической ситуации.

Функция дожития S(x) выбрана в виде модели Гомперца-Мейкема

S(x) = expj- J (А + ВеСи = ехр{- Ах - В(еСх - 1)/с} где параметр А учитывает риски, связанные с несчастными случаями,

Сх

а слагаемое Be учитывает влияние возраста на смертность.

Интенсивность процесса рождаемости b(x,t) определим в виде

А(х,0 = 0,488т\(0ч'(х,0.

где r|(i) - суммарный коэффициент рождаемости в момент времени I, а функция имеет смысл плотности распределения вероятностей

значений репродуктивного возраста женщины и задана как плотность смещенного двупарамерического у-распределения.

На основе сделанных предположений относительно динамики суммарного коэффициента рождаемости r|(i) в разделе 4.3 построены пессимистический и оптимистический сценарии развития демографической ситуации в РФ на долгосрочную перспективу. Найдены дисперсия численности возрастных групп, а также ковариация между численностью возрастной группы нулевого возраста и численностью возрастной группы возраста л.

Заметим, что автономная немарковская система массового обслуживания с неограниченным числом приборов и ее упрощенные варианты, а также методы их исследования, предложенные в диссерта-

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

В Заключении диссертации приведены основные результаты.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

В журналах перечня ВАК (редакция от апреля 2008 года):

1. Назаров A.A. Многофазная автономная система массового обслуживания и ее применение к задачам демографии / A.A. Назаров, М.Г. Носова // Известия Томского политехнического университета. Серия Управление, вычислительная техника и информатика. - 2009. -Том 315, №5. —С. 183-186.

2. Назаров A.A. Исследование математической модели демографических процессов в виде пятифазной системы массового обслуживания / A.A. Назаров, М.Г. Носова // Вестник Сибирского государственного аэрокосмического университета им. академика Решетнева. -2010.-Том 1.-С. 53-58.

3. Назаров A.A. Математическая модель процесса изменения демографической ситуации и ее исследование / A.A. Назаров, М.Г. Носова // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2009. - Т. 2 (20). - С. 100-105.

В других изданиях:

4. Назаров A.A. Стохастическая модель демографических процессов как автономная система массового обслуживания / A.A. Назаров, М.Г. Носова // Обозрение прикладной и промышленной математики. - 2009. - Т. 16, Вып. 6. - С, 1098-1099.

5. Назаров A.A. О нецелесообразности аппроксимации процесса рождаемости потоками Пуассона при долгосрочном прогнозировании / A.A. Назаров, М.Г. Носова // Вестник Томского государственного университета. Серия Управление, вычислительная техника и информатика. - 2009. - №3(8). - С. 75-80.

6. Назаров A.A. Метод передвижки возрастных групп в демографии и его приложения / A.A. Назаров, М.Г. Носова // Вестник Томского государственного университета. Серия Управление, вычислительная техника и информатика. - 2009. — №3(8). - С. 67-75.

7. Носова М.Г. Статистическая инвариантность во времени распределения вероятностей значений репродуктивного возраста // Информационные технологии и математическое моделирование (ИТММ-2004): Материалы III Всероссийской научно-практической конференции. - Томск: ТГУ, 2004. - Ч. 2. - С. 32-33.

8. Назаров A.A. Исследование демографических процессов методами теории массового обслуживания / A.A. Назаров, М.Г. Носова // ММПЭИТС: Материалы международной научной конференции. -Минск: БГУ, 2005. - С. 156-161.

9. Морозова A.C. Исследование основных характеристик математической модели сценария развития демографической ситуации / A.C. Морозова, М.Г. Носова // Научное творчество молодежи: Материалы IX Всероссийской научно-практической конференции. - Томск: ТГУ, 2005. - Ч. 1. - С. 42-45.

10. Назаров A.A. Долгосрочный прогноз значений процесса изменения демографической ситуации / A.A. Назаров, М.Г. Носова // Научное творчество молодежи: Материалы IX Всероссийской научно-практической конференции. - Томск: ТГУ, 2005. - Ч. 1. - С. 37-40.

11. Назаров A.A. Упрощенная стохастическая модель демографии в виде автономной системы обслуживания с неограниченным числом приборов / A.A. Назаров, М.Г. Носова // Информационные технологии и математическое моделирование: Материалы V Международной научно-практической конференции. - Томск: ТГУ, 2006. -С. 134-137.

12. Назаров A.A. Исследование системы с РН-распределением времени обслуживания / A.A. Назаров, М.Г. Носова // Научное творчество молодежи: Материалы XI Всероссийской научно-практической конференции. - Томск: ТГУ, 2007. - Ч. 1. - С. 41-44.

13. Назаров A.A. Автономная немарковская система обслуживания с неограниченным числом приборов / A.A. Назаров, М.Г. Носова // ММПЭИТС: Материалы международной научной конференции. -Минск: РИВШ, 2007. - С. 175-180.

Отпечатано в ООО «НИП» г. Томск, ул. Советская, 47, «Штемпельная Мастерская» тел.-. 53-14-70, 52-83-10 Заказ № 476-34, Тираж 125 экз.

Оглавление.

Введение.

Глава 1 Некоторые математические аспекты применения метода математического моделирования в демографии.

1.1 Метод прямой передвижки возрастных групп.

1.2 Метод обратной передвижки возрастных групп.

1.3 Применение методов прямой и обратной передвижки возрастных групп для оценки величины людских потерь РФ в годы ВОВ.

1.4 О нецелесообразности аппроксимации процесса рождаемости потоком Пуассона при долгосрочном прогнозировании.

Резюме.

Глава 2 Упрощенные автономные системы массового обслуживания в демографии.

2.1 Автономная упрощенная стохастическая система массового обслуживания.

2.2 Исследование стохастической модели демографии в виде трехфазной автономной системы массового обслуживания.

2.2.1 Метод моментов.

2.2.2 Асимптотики первого и второго порядков.

2.2.3 Применение трехфазной системы массового обслуживания к задачам демографии.

2.3 Исследование стохастической модели демографии в виде пятифазной автономной системы массового обслуживания.

2.3.1 Метод моментов.

2.3.2 Асимптотики первого и второго порядков.

2.3.3 Применение пятифазной системы массового обслуживания к задачам демографии.

2.4 Исследование стохастической модели демографии в виде автономной системы массового обслуживания с двумя типами заявок.

2.4.1 Метод моментов.

2.4.2 Асимптотики первого и второго порядков.

2.4.3 Применение системы массового обслуживания с двумя типами заявок к задачам демографии.

Резюме.

Глава 3 Стохастическая модель демографии в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов.

3.1 Автономная немарковская система массового обслуживания с неограниченным числом приборов.

3.2 Автономная система массового обслуживания с РН-распределением времени обслуживания.

3.3 Исследование автономной системы массового обслуживания с РН-распределением времени обслуживания.

3.3.1 Метод моментов.

3.3.2 Асимптотика первого порядка.

3.3.3 Асимптотика второго порядка.

3.4 Автономная система массового обслуживания с произвольным временем облуживания.

3.4.1 Асимптотика первого порядка.

3.4.2 Асимптотика второго порядка.

Резюме.

Глава 4 Применение модели и методов теории массового обслуживания к исследованию процесса изменения демографической ситуации.

4.1 Функция дожития. Аналитический закон смертности в виде модели Гомперца-Мейкема (Gompertz-Makeham).

4.2 Определение интенсивности процесса рождаемости.

4.3 Сценарий развития демографической ситуации. Численные результаты.

4.3.1 Пессимистический сценарий.

4.3.2 Оптимистический сценарий.

Резюме.

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

Среди методов демографического анализа выделяют следующие основные группы [16]: статистические методы изучения населения, социологические методы, графоаналитические и картографические методы, использование математических моделей в изучении населения. Использование математических моделей относится к современному методу демографического анализа и заключается в определении некоторых исходных предположений относительно динамики демографических процессов.

Для демографического прогнозирования численности и состава населения полезными являются статистические методы и методы математического моделирования (рисунок 1). Однако наиболее эффективным инструментом является математическое моделирование.

Рисунок 1 — Классификация методов, используемых для демографического прогнозирования

Остановимся более подробно на известных методах математического моделирования демографических процессов.

В зависимости от того, учитывает ли демографическая модель отклонение частот демографических событий от их вероятностей, демографические модели принято делить на стохастические (или вероятностные) и детерминированные. Также демографические модели подразделяют на непрерывные и дискретные.

К первым математическим моделям воспроизводства относятся модели воспроизводства биологических популяций. Простейшей моделью является задача о размножении одной пары кроликов, которая была описана в 1202 г. Leonardo Pisano (около 1170 - около 1250 гг.), известным как Фибоначчи

Fibonacci) [94]. Согласно его теории рост популяции происходит за счет рождения каждой парой кроликов новой пары.

Следующей по времени моделью общей численности популяции в демографии считается модель английского ученого Т. Мальтуса (Thomas Robert Malthus) (1766-1834 гг.) [15, 17, 31, 54, 57-58, 94, 103]. В модели Мальтуса население имеет экспоненциальный рост. Заметим, что модели популяции Мальтуса и популяции Фибоначчи являются детерминированными непрерывными моделями, в частном случае дискретными моделями.

Что касается простых детерминированных моделей роста человечества, то к ним относятся, прежде всего, модели линейного и экспоненциального роста [1, 17, 24, 36, 57-58, 100-101, 126, 156]. В модели линейного роста изменение численности населения происходит при постоянстве абсолютного прироста. Следует отметить, что модель линейного роста дает удовлетворительные результаты только на короткий период, продление же на более длительный срок ни в прошлое, ни в будущее не дает адекватных результатов.

В литературе также встречается стационарная модель воспроизводства населения [12, 17-18, 52, 99, 103, 120]. В такой модели предполагается, что рождаемость равна смертности, то есть при постоянстве численности прирост равен нулю. Таким образом, в стационарном населении остаются постоянными общая численность и численность в возрастно-половых группах.

Наиболее известный подход, позволяющий описать модель воспроизводства населения в целом, — теория стабильного населения [13, 17—19, 24, 26, 28, 30, 36, 43, 52, 57-58, 60, 94, 96, 101, 103, 120, 125, 139-140, 145]. Стабильное население характеризуется постоянными во времени возрастными интенсивно-стями рождаемости, смертности и возрастной структурой населения. Население, в котором не изменяются функция дожития, возрастная структура населения и коэффициент прироста называется экспоненциальным населением. Ему соответствуют все основные формулы теории стабильного населения, в которые не входит функция рождаемости. Очевидно, что стационарное население -частный случай стабильного населения.

Известны непрерывные и дискретные аналоги модели стабильного населения. В основе непрерывных моделей лежит интегральное уравнение воспроизводства населения (уравнение Лотки), в основе дискретных — матричная модель (матрица Лесли).

Идея экспоненциального роста населений как основы стабильного населения впервые представлена английским актуарием М. Хейлом (М. Hale). Разработку теории стабильного населения связывают с такими именами как Л. Эйлер (L. Euler) (1707-1783 гг.), Г. Кнапп (G. Кпарр) (1842-1926 гг.), В. Лексис (V. Lexis) (1837-1914 гг.), Дж. Лотка (J. Lotka) (1880-1949 гг.), В. Борткевич (1868—1931 гг.). Заметим, что Дж. Лотка разработал теорию стабильного населения в непрерывной форме для одного пола. Дискретную теорию стабильного населения, используя дискретный подход У. Феллера (W. Feller) (1906—1970 гг.), построил П. Лесли (P. Leslie) [137]. Обобщением теории стабильного населения занимались П. Венсент (P. Vincent), Дж. Хаджнал (J. Hajnal) (1924-2008 гг.), П. Кармел (P. Carmel), К. Джини (К. Gini) (1884-1965 гг.), Ж. Буржуа-Пиша (J.Bourgeois-Pichat) (1912-1990 гг.) [111], Л. Анри (L. Anry), Н. Кейфиц (N. Keyfltz) [127-133], А. Я. Боярский (1906-1985 гг.), а также X. Хюрениус (Н. Hyrenius) (1914-1979 гг.), А. Роджерс (A. Rogers) [148-149] и В. Ф. Шукайло (1932-1983 гг.). Большой вклад в разработку методов практического применения стабильного населения привнесли советские демографы С. А. Новосельский (1872-1953 гг.), В. В. Паевский (1893-1934 гг.), А. Я. Боярский (1906-1985 гг.), И. Г. Венецкий (1914-1981 гг.) и др. [5-9, 17, 26, 56].

Современный этап развития теории стабильного населения связывают с обобщением основных выводов на случай демографических процессов с переменными интенсивностями (процессы рождаемости и смертности в виде временных рядов, случайных процессов) - А. Коул (A. Coale) (1917-2002 гг.) [114, 144] и А. Лопес (A. Lopes). А также Р. Лии (R. Lee) [135-136], Д. Ахлбург (D. Ahlburg) [107], 3. Сайкес (Z. Sykes) [153-154], Картер (Carter), М. Алхо (М. Alho) [108-109], Б. Спенсер (В. Spencer) [108], Дж. Поллард (J. Pollard) [142143], Дж. Кохен (J. Cohen) [116-118], Н. Кейфица (N. Keyfltz) [127-133], X.

Касвелл (Н. Caswell) [112-113], М. Тревисан (М. Trevisan) [113], Л. Гудман (L. Goodman) [122-124] и другие [26, 110, 115, 119]. Применение теории стабильного населения встречается в работах [3, 18, 23, 32].

В конце ХЕК века учеными была предпринята попытка использовать в качестве модели, описывающей рост человечества, логистическую кривую [17, 24, 29, 57-58, 101, 103, 113, 120, 126, 139, 146]. Этой модели придерживались Р. Пирль (R. Pearl) (1879-1940 гг.) и JL Рид (L. Reed) [146]. Уравнение этой кривой было введено впервые в 1835 г. бельгийским математиком Ф. Ферхюльстом (F. Verhulst) (1804—1849 гг.), общая идея которого состояла в наложении на экспоненциальный рост численности населения некоторого фактора, замедляющего этот рост и увеличивающего свое действие по мере возрастания численности населения [27, 30]. Надо иметь в виду, что использование Р. Пирлем и другими демографами логистической кривой для абстрактных расчетов численности населения означало игнорирование влияния социально-экономических факторов на численность населения. В связи с этим попытки использовать логистическую кривую для прогнозирования численности населения на достаточно отдаленное будущее привели к ненадежным и неправильным выводам. Логистическая кривая может быть полезной лишь при краткосрочном прогнозе численности населения, а также как приближение динамики некоторых демографических показателей, например рождаемости, смертности, функции дожития и т. д. [134].

Еще одной детерминированной моделью воспроизводства населения является модель гиперболического изменения численности населения. Впервые предположение, что скорость роста численности населения пропорциональна квадрату численности, было замечено в 1960 г. Хайнцом фон Ферстером (V. Foerster) (1911-2002 гг.), П. Мором (P. Mora) и Л. Амиотом (L. Amiot). По этой же формуле И. С. Шкловский прогнозировал численность населения земного шара в период с 1600 по 1960 г. [17]. В настоящее время в этом направлении работает С. П. Капица [34—41]. Модель Ферстера встречается также в [46].

Значительный вклад в математическую демографию сделан Олегом Васильевичем Староверовым (1933-2006 гг.) [94-96], который рассматривал демографические процессы в виде марковских моделей в форме цепей Маркова. Староверов О. В., применяя теорию пуассоновских потоков, предполагал, что все демографические процессы, включая миграцию, обладают вероятностными свойствами: стационарностью, отсутствием последействия, ординарностью. Из простейших предположений Староверовым О. В. получены модели демографии, миграции, межотраслевого и социального движения населения. Модели естественного движения населения исследованы в дискретном и непрерывном виде, для них в [94, 96] получены основные уравнения.

В [95] Староверовым О. В. предлагается стохастическая модель развития населения с дискретным временем, учитывающая! случайность, как в рождаемости, так и в смертности. Суть исследования с помощью такой модели состоит в следующем: вводится понятие состояние индивидуума' в момент времени, определяемое возрастной группой, в которой он находится, и определяются возможные- переходы и вероятности этих переходов. Кроме этого, применение марковских цепей встречается в [42, 83-85, 121].

В отдельную группу можно > выделить модели, которые рассматривают демографический процесс изменения численности как ветвящийся процесс, в частном случае как процесс гибели и размножения [122-124, 142-143, 147].

Особое место в демографическом моделировании занимает метод передвижки возрастов [10-12, 17, 57-58, 105, 138, 145, 151, 157-158]. Метод компонент или метод передвижки разработан П. К. Уэлптоном (Р.К. Whelpton) (1893— 1964 гг.) [157—158]. В этом методе за основу принимают распределение населения по возрастам и постепенно передвигают численности отдельных возрастных групп в соответствии с показателями таблиц смертности.

В России перспективными расчетами численности населения методом возрастной передвижки занимались С. Г. Струмили» (1877-1974 гг.), А. Я. Боярский (1906-1985 гг.), П. П. Шушерин, М. С. Бедный [12, 48, 103].

В последнее время расчетами половозрастной структуры населения методом возрастной передвижки занимаются Государственный комитет Российской Федерации по статистике, Центр демографии и экологии человека, Отдел населения ООН, Международный Институт прикладного системного анализа (IISA) в г. Лаксенбург (Австрия), В. Лутц (W. Lutz), У. Сандерсон (Sanderson W.) и С. Щербов [19, 25, 57—58, 77-80]. Применение метода передвижки встречается в работах [2, 13, 21, 44, 51, 55, 86, 88-90, 92, 99], а также [17, 23, 100, 107, 109, 152].

Следует заметить, что соотношения для модели естественного движения населения в дискретном времени (матричная модель - матрица Лесли), полученные также О.В. Староверовым [94, 96], являются, так называемой передвижкой возрастов, то есть применение подхода Лесли приводит к тем же результатам, что и передвижка по возрастам. А модель Лотки есть не что иное, как аналог метода передвижки в непрерывном времени.

Общая теория стабильного населения на случай демографических процессов с переменными интенсивностями, применяемая А. Коулом, А. Лопесом, Р. Лии, Д. Ахлбургом, 3. Сайкесом, Картером, М. Алхо, Б. Спенсером, Дж. Поллардом, Дж. Кохеном, Н. Кейфицем, X. Касвеллом, М. Тревисаном, Л. Гуд-маном и другими, выделяются здесь в отдельный класс методов - стохастический вариант метода передвижки по возрастам [107-109, 112—114, 116-118, 122-124, 127-133, 135-136, 142-144, 153-154].

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

Подведем итог, проведенный анализ множества моделей, показывает, что в моделировании демографических процессов наиболее распространены детерминированные модели (дискретные и непрерывные) и стохастические дискретные. Преимущество стохастических моделей перед детерминированными моделями заключается в том, что они учитывают отклонение частот демографических событий от их вероятностей.

Однако демографические процессы протекают в непрерывном времени и являются, стохастическими. Методы исследования таких процессов в математической демографии недостаточно развиты. Именно поэтому актуальной является задача существенного расширения^ математических моделей процесса изменения демографической ситуации, а также развития методов их исследования.

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

Таким образом, поставлены следующие задачи:

1. Разработка математического метода исследования* автономных стохастических систем с непрерывным и дискретным временем.

2. Разработка стохастических моделей демографии в- виде упрощенных автономных систем массового обслуживания1 с РН-распределением. времени обслуживания и развитие методов моментов и асимптотического анализа для их исследования.

3. Разработка новой математической модели процесса- изменения демографической ситуации в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем, обслуживания.

4. Разработка метода виртуальных фаз и модификация метода асимптотического анализа для исследования автономной немарковской системы массового обслуживания, с неограниченным числом приборов и произвольным временем обслуживания.

5. Применение разработанных математических моделей к исследованию процесса изменения демографической ситуации.

Научная новизна и результаты, выносимые на защиту, состоят в следующем:

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

2. Разработан метод виртуальных фаз. для-исследования автономной немарковской системы массового обслуживания1 с неограниченным числом приборов, позволяющий находить многомерное распределение вероятностей чис-ленностей обслуживаемых заявок различных возрастов и основные характеристики, определяющие эти распределения.

3. Для исследования упрощенных моделей (трехфазной"и пятифазной автономных систем, автономной системы с двумя типами заявок,, автономной системы с РН-распределением времени обслуживания) предложены модификации метода моментов и метода асимптотического г анализа; Найдены распределения вероятностей числа обслуживаемых заявок в таких системах и> их основные характеристики. 4. На основе метода прямой передвижки возрастных групп предложен метод обратной передвижки, позволяющий выполнять интерполяцию демографических данных между датами переписей населения, а также восстанавливать функцию интенсивности процесса рождаемости. Метод прямой и обратной передвижки возрастных групп применены для оценки величины людских потерь РФ в годы ВОВ. Значение оценки составляет 17 млн. человек.

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

Результаты, представленные в работе, имеют как теоретическое, так и практическое значение.

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

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

Публикации. По материалам диссертации опубликованы 13 работ, из них 3 статьи [65, 69, 71] в журналах перечня ВАК:

В журналах перечня ВАК:

1. Назаров А. А. Многофазная автономная система массового обслуживания и ее применение к задачам демографии / А. А. Назаров, М. Г. Носова // Известия Томского политехнического университета. Серия Управление, вычислительная техника и информатика. - 2009. — Том 315, №5. - С. 183-186.

2. Назаров А. А. Исследование математической модели демографических процессов в виде пятифазной системы массового обслуживания / А. А. Назаров, М. Г. Носова // Вестник Сибирского государственного аэрокосмического университета им. академика Решетнева. — 2010. — Том 1. - С. 53-58.

3. Назаров А. А. Математическая модель процесса изменения демографической ситуации и ее исследование / А. А. Назаров, М. Г. Носова // Доклады Томского государственного университета систем управления и радиоэлектроники. - 2009. - Т. 2 (20). - С. 100-105.

В других изданиях:

4. Назаров А. А. Стохастическая модель демографических процессов как автономная система массового обслуживания / А. А. Назаров, М. Г. Носова // Обозрение прикладной и промышленной математики. — 2009. — Т. 16, Вып. 6. — С. 1098-1099.(

5. Назаров А. А. О нецелесообразности аппроксимации процесса рождаемости потоками Пуассона при долгосрочном прогнозировании / А. А. Назаров, М. Г. Носова // Вестник Томского государственного университета. Серия Управление, вычислительная техника и информатика. — 2009. — №3(8). — С. 75— 80.

6. Назаров А. А. Метод передвижки возрастных групп в демографии и его приложения / А. А. Назаров, М. Г. Носова // Вестник Томского государственного университета. Серия Управление, вычислительная техника и информатика. — 2009.-№3(8).-С. 67-75.

7. Носова М. Г. Статистическая инвариантность во времени распределения вероятностей значений репродуктивного возраста // Информационные технологии и математическое моделирование (ИТММ-2004): Материалы Ш Всероссийской научно-практической конференции. - Томск: ТТУ, 2004. — Ч. 2. — С. 32-33.

8. Назаров А. А. Исследование демографических процессов методами теории массового обслуживания / А. А. Назаров, М. Г. Носова // ММПЭИТС: Материалы международной научной конференции. — Минск: БГУ, 2005. - С. 156-161.

9. Морозова А. С. Исследование основных характеристик математической, модели сценария развития демографической ситуации / А. С. Морозова, М. Г.

Носова // Научное творчество молодежи: Материалы IX Всероссийской научно-практической конференции. - Томск: ТГУ, 2005. — Ч. 1. — С. 42-^45.

10. Назаров А. А. Долгосрочный прогноз значений процесса изменения демографической ситуации / А. А. Назаров, М. Г. Носова // Научное творчество молодежи: Материалы IX Всероссийской научно-практической конференции. -Томск: ТГУ, 2005. - Ч. 1. - С. 37^0.

11. Назаров А. А. Упрощенная стохастическая модель демографии в виде автономной системы обслуживания с неограниченным числом приборов / А. А. Назаров, М. Г. Носова // Информационные технологии и математическое моделирование: Материалы V Международной научно-практической конференции. - Томск: ТГУ, 2006. - С. 134-137.

12. Назаров А. А. Исследование системы с РН-распределением времени обслуживания / А. А. Назаров, М. Г. Носова // Научное творчество молодежи: Материалы XI Всероссийской научно-практической конференции. — Томск: ТГУ, 2007. - Ч. 1. - С. 41-44.

13. Назаров А. А. Автономная немарковская система обслуживания с неограниченным числом приборов / А. А. Назаров, М. Г. Носова // ММПЭИТС: Материалы международной научной конференции. — Минск: РИВШ, 2007. — С. 175-180.

Апробация работы. Основные положения работы и отдельные ее результаты докладывались и обсуждались:

1. III Всероссийская научно-практическая конференция «Информационные технологии и математическое моделирование (ИТММ-2004)», г. Анжеро-Судженск, 2004 г.

2. Международная научная конференция «Математические методы повышения эффективности функционирования телекоммуникационных сетей», г. Минск, 2005 г.

3. IX Всероссийская научно-практическая конференция «Научное творчество молодежи», г. Анжеро-Судженск, 2005 г.

4. V Международная научно-практическая конференция «Информационные технологии и математическое моделирование», г. Анжеро-Судженск, 2006 г.

5. XI Всероссийская научно-практическая конференция «Научное творчество молодежи», г. Анжеро-Судженск, 2007 г.

6. Международная научная конференция «Математические методы повышения эффективности функционирования телекоммуникационных сетей», г. Минск, 2007 г.

7. VII Российская конференция с международным участием «Новые информационные технологии в исследовании сложных структур», г. Томск, 2008 г.

Работа выполнена при поддержке АВЦП «Развитие научного потенциала высшей школы (2009—2010 гг.)» Федерального агентства по образованию РФ по проекту «Разработка методов исследования немарковских систем массового обслуживания и их применения к сложным экономическим системам и компьютерным сетям связи».

Структура работы. Работа состоит из введения, четырех глав, заключения, списка использованной литературы и приложений.

В первой главе предлагаются решения некоторых задач математической демографии с помощью математических методов.

В разделе 1.1 дано математическое обоснование метода прямой передвижки возрастных групп [157-158, 70]. Метод прямой передвижки применяется для определения оценок значений численности 7V(x + x,/ + x) группы лиц возраста х-Н; в году t+z при условии, что известна численность N(pc,t), т - шаг прогнозирования. Получены основные равенства для применения метода передвижки возрастных групп. В частности, для метода прямой передвижки основное равенство записывается в виде

N(x + x,t + x) = N(x, t) + б2 > где S(x) — функция дожития [102], N(x,t) задано, a N(x + x,t + т) является оценкой значения численности демографической группы лиц возраста х+т в году t+т, а 82 - случайная ошибка с математическим ожиданием Ме2=0.

Показано, что оценка значения численности возрастных групп по методу прямой передвижки обладает достаточно высокой точностью, поэтому ошибкой е2 здесь можно пренебречь.

Введено понятие обратной передвижки возрастов. Для метода обратной передвижки основное равенство записывается в виде

Щх- z,t - т) = Щх,t) + е4>

S(x) где N(x,t) задано, a N(x - x,t - т) является оценкой значения численности демографической группы лиц возраста х-т в году t—т, е4 - случайная ошибка с математическим ожиданием Ме4=0. Данное равенство позволяет определить оценку значения численности демографической группы в прошлые моменты времени.

В разделе 1.2 проведено дополнительное исследование оценки значения численности N(x-i,t -т), полученной методом обратной передвижки возрастных групп. Найдены условные математическое ожидание и дисперсия величины N{x — x,t — x) при условии, что выполняется равенство N(x,f)=n.

Заметим, что в смысле значений коэффициентов вариации, показано, что оценки, полученные методом обратной передвижки, на порядок (в 10 раз) точнее оценок, полученных методом прямой передвижки, при одинаковом горизонте прогнозирования т. Поэтому случайной ошибкой е4 здесь также можно пренебречь.

В разделе 1.3 на основе полученных формул применены методы прямой и обратной передвижки возрастов к решению задачи оценки величины людских потерь Российской Федерации в годы Великой Отечественной Войны. Показано, что суммарные значения людских потерь РФ в годы ВОВ при методах прямой и обратной передвижки возрастных групп совпадают и составляют 5=17 млн. чел.

В разделе 1.4 представлено математическое обоснование нецелесообразности аппроксимации в математической демографии процесса рождаемости потоком Пуассона при долгосрочном прогнозировании [72, 75, 96]. Для этого на первом этапе рассматривается наиболее простая для исследования математическая модель случайного потока, управляемого цепью Маркова. Такая модель определяется как поток заявок в автономной системе массового обслуживания с неограниченным числом приборов и средним значением m(f) числа занятых приборов. Для системы определено распределение вероятностей P(i,t) и составлена система дифференциальных уравнений Колмогорова, найден аналитический вид характеристической функции H(u,t)

•л N i±i£ 1)ехр{(Ь-[i)t]

H(u,t) b b eju - H (eju -1) exp{(£ - ц)*} b

Далее исследуется, открытая система массового обслуживания с неограниченным числом приборов, на вход которой поступает нестационарный пуас-соновский поток с интенсивностью выбираемой из условия совпадения средних значений числа приборов, занятых в автономной системе и в открытой системе с пуассоновским входящим потоком. Аналогично предыдущей системе массового обслуживания записано и решено дифференциальное уравнение для характеристической функции числа занятых приборов G(u;t). Характеристичная функция G(u,t) имеет вид

G{u,t) = {l-e^ +ejue-v)N ехр{(<> -l)}

Через обратное преобразование-Фурье [14] найдены распределения вероятностей P\(i,t) и Pi(i,t); определяемые характеристическими функциями G(u,t) и H(u,t). Определено расстояние между распределениями вероятностей P\(i,t) и Р2(л0 как p(b,t) - max

UPx(n,t)-P2{n,t)) п=О

Получены значения расстояния p{b,t) для различных значений параметров but, показано, что расстояние р (b,t) возрастает с ростом значений параметров but. Поскольку, характеристики таких систем массового обслуживания существенно различаются, следовательно, делается вывод о нецелесообразности аппроксимации процесса рождаемости потоком Пуассона при долгосрочном прогнозировании.

В главе 2 [65, 71, 74] предлагается математическая модель женского населения в виде автономной немарковской системы массового обслуживания с РН-распределением времени обслуживания и неограниченным числом приборов. Особенностью такой системы является отсутствие внешнего источника1 заявок.

Продолжительность обслуживания т каждой заявки складывается из про-должительностей случайного числа фаз

Т = Т] Т2 где т,- — продолжительность /-й фазы обслуживания. Величины т, являются независимыми экспоненциально распределенными случайными величинами с параметрами ц,/, характеризующими продолжительности фаз обслуживания, i - U v .

В терминах демографии обслуживаемая заявка интерпретируется как женщина, время обслуживания заявки — продолжительность жизни этой женщины, величина b,(t) — интенсивность рождения девочек у женщины i-й возрастной группы в году t. Входящим потоком заявок является* процесс рождения девочек, то есть последовательность моментов рождения девочек от всей совокупности женщин.

В этой главе рассматриваются некоторые упрощенные варианты автономной системы с РН-распределением времени обслуживания и методы их исследования, а именно, трехфазная система массового обслуживания (раздел 2.2), пятифазная система массового облуживания (раздел 2.3) и система массового обслуживания с двумя типами заявок (раздел 2.4).

Различия между моделями следующие: в трехфазной модели репродукция новой-заявки, возможна только на 2-й фазе с интенсивностью b(t)=b, в пятифазной - на 2-й, 3-й и 4-й фазах с интенсивностью b\{t)=b\, b2(t)=b2 и Ь3(г)=6з соответственно. В системе массового обслуживания с двумя типами заявок рассматриваются заявки первого типа, продолжительность обслуживания, которых в системе состоит из трех фаз, и заявки второго типа, их время обслуживания состоит из единственной фазы. В терминах демографии под заявкой первого типа подразумевается женщина, а под заявкой второго типа - мужчина, репродуктивная 2-я фаза - репродуктивный возраст женщины, а время обслуживания в системе - продолжительность жизни человека. В перечисленных системах массового обслуживания рассматриваются фазы разной продолжительности.

Состояние упрощенной автономной системы в момент времени t определяется вектором n{t) соответствующей размерности. Для указанных моделей определены многомерные распределения! вероятностей P(n,t), где п={п\,п2, и3,.}, для которых составлены системы дифференциальных уравнений Колмогорова [76] и записаны дифференциальные уравнения для многомерных характеристических функций H{u,t) числа заявок, обслуживаемых в системах, здесь и—{щ, г/2, из,.}.

Для нахождения распределения вероятностей, а также основных его характеристик применены методы моментов и асимптотического анализа [62].

Метод моментов позволил найти числовые характеристики (моменты) многомерного случайного*процесса n(t). Для исследования распределения вероятностей его значений использован метод асимптотического анализа, модифицированный к рассматриваемой задаче.

При методе асимптотического анализа уравнение для характеристической функции H{u,t) рассматривалось в асимптотическом условии большой численности групп, пропорциональных некоторой бесконечно большой величине N , что- вполне обосновано в демографических исследованиях, где численности стандартных возрастных групп измеряются миллионами.

В1 первой асимптотике были предложены замены u = ew, H(u,t) = F\(w,t,z), где £=1 IN, е — малый положительный параметр [53], которые позволили перейти к сингулярно возмущенным по малому параметру 8 задачам для Fi(w,t,E).

Считаем, что получаемые в результате замен сингулярно возмущенные задачи обладают некоторыми особенностями, которые приводят к тому, что для решений Fi(w,t,е) сингулярно возмущенных задач существуют пределы F\(w,f), которые определяются решением предельных невозмущенных задач. Данное предположение относится ко всем сингулярно возмущенным задачам, рассматриваемым в диссертации.

Показано, что где аналитическое среднее Nm(t) является решением системы обыкновенных дифференциальных уравнений первого порядка.

Для нахождения асимптотики второго порядка в задаче для Fx{w,t,£) выполнены замены где элементы матрицы ковариаций NR=[NRjj] являются также решением системы обыкновенных дифференциальных уравнений первого порядка.

Определено, что многомерная величина n(t) имеет асимптотически нормальное распределение с вектором математических ожиданий Nm(t) и матрицей ковариаций NR=[NRij\, где Rij=M{}, очевидно, что Rf=Rji. Величины i и j определяются типом системы массового обслуживания.

На примере упрощенных стохастических моделей сформулированы основные идеи метода асимптотического анализа в условии большой численности lim Fy (w,t, e) = Fx (w,t) = exp {jm(t)w}, групп, которые изложены в доказательствах Теорем 2.1 и 2.2 раздела 2.2, Теорем 2.3 и 2.4 раздела 2.3, Теорем 2.5 и 2.6 раздела 2.4. Все рассмотренные асимптотики данной диссертации реализованы в условии большой численности групп.

Отметим, что асимптотические средние значения Nm\{t) числа заявок, обслуживаемых в рассмотренных системах в момент времени t, найденные в предельном асимптотическом условии JV—>оо, совпадают с допредельными, полученными методом моментов. Заметим, что такое совпадение справедливо и для вторых моментов в предельной и допредельной моделях.

Предложенная упрощенная стохастическая модель в виде автономной системы массового- обслуживания с РН-распределением времени обслуживания, а также частные ее виды, применяется для анализа демографической ситуации. В таком случае, аналитические результаты метода моментов и метода асимптотического анализа позволяют построить долгосрочный сценарий развития демографической ситуации, как это продемонстрировано в разделах 2.2.3, 2.3.3, 2.4.3. t

В главе 3 [63, 66-69, 73] предложена стохастическая модель в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания. Принципиальное отличие от модели главы .2 заключается в том, что здесь исследуется стохастическая плотность численности заявок. Предложен метод ее исследования, названный методом виртуальных фаз, заключающийся в аппроксимации времени обслуживания суммой случайного числа независимых одинаково экспоненциально распределенных случайных величин и предельным переходом при неограниченном возрастании числа фаз, а также пропорциональном уменьшении длительности каждой фазы.

В разделе 3.2 исследована общая система массового обслуживания с РН-распределением-времени обслуживания. В этой модели время обслуживания аппроксимировано суммой случайного числа независимых одинаково экспоненциально распределенных случайных величин. Для распределения вероятностей P(ni,n2,.,t) числа обслуживаемых заявок записана система дифференциальных уравнений Колмогорова [76] P(nl-l,n2,.,t){(nl-l)bl{t) + 1 ntbt{[)} + 2

00 ЦХ (и,- + \){P{nx,n2,.,ni +1/2/+1,.,0(1-//) +

1=1 Р(щ,П2,.,П; + 1,И/+1 -1 ,Л/+2.,*>/}, где

Г Л (i-i\

S

1 м- J вероятность перехода заявки на следующую фазу, a .S'(x)=l—fi(x), где 5(х) — функция распределения времени обслуживания заявки.

Наряду с исследованием допредельной модели, предложено ее асимптотическое исследование. Асимптотическое условие, а также замены асимптотик первого и второго порядков аналогичны главе 2. Основные идеи метода асимптотического анализа изложены в доказательствах Теорем 3.1 и 3.2. Найдена асимптотика первого порядка оо

Hx(u,t) = exp{/iV£>v/wv(/)}>

V=1 где компоненты m^t) являются решением системы обыкновенных дифференциальных уравнений первого порядка + t),

I '=1 т; (0 = -\xnii (0 + (0, /' > 1.

Асимптотика второго порядка определена в виде

00 J 00 00

Н2 (и, 0 = exp {jN^uvmv (t) --Nj^ Z uiRy (*)uj }» v=i 2 i=i j=i где ^ = 1,°°, у = l,co, определяются из неоднородной системы обыкновенных дифференциальных уравнений первого порядка

ЯМО = -2\iRn(t)+Y,bv(t){Rlv(t)+Rvm+qu(t),

V=1

R'n (0 = -2^1 (0+Д/-11 (0 ■+ X К (0 ** (0 + 4i\ (0,

V—1

00

R[j Сf) = -2цЯ1у (0 + (0 + E К (t)RVJ (0 + <7U (0,

V=1

0 = -2[iRjj (t) + м^-Л-и (О + K/-1 Vi W + )» а функции имеют вид

00 n(0 = №(0+X>v(0&v(0>

V = 1

0 = Vv+lv (0 = OK > V = 1, 00 , tfw (0 = \unv (0 + (O'V-I» V = 2, GO,

0 = 0, N>2, /,y=2,oo.

Заметим, что равенства, определяющие средние характеристики в допредельной и предельной моделях совпадают.

Показано, что в условии предельного уменьшения длительностей фаз, то есть при ц—>оо, рассматриваемое РН-распределение сходится к распределению, определяемому функцией S(x), где S(x) — функция дожития, вид которой можно выбирать достаточно произвольно. Поэтому, сделав замены щ=и{И\х), nt{t) = n(i/\x,t),

Н(и, t) = Н(и, t, ц) = М< ехр

00 ( I \ (i \ 1 т

I 1 -,t vH- ) И и полагая, что при (i—>оо, г—>оо отношение И\ь—>х и существуют следующие пределы lim u(i/\x) = и(х),

7ц->дс lim Н(и

1->оо i/\i->x ехр V jju(x)^(x,t)dx \ = F(u,t), V о / а также существует предел по распределению lim \m (i I \x., = x, t), в разделе 3.4 перешли к рассмотрению автономной системы массового обслуживания с произвольным временем обслуживания. Случайная функция ££x,f) названа стохастической плотностью численности заявок возраста х в момент времени t. Записано основное уравнение для характеристического функционала

Г ( 00 Л

F(u,t) = M<exр j\u(x)ll(x,t)dx I V о вок стохастической плотности численности заяdF(u,t) j*jdF(u,t)

1 - б>(0) }Ь(х, 0 + (е~Мх) -1) - ju'(x)\clx.

S(x) I dt 0 ди(х)

Для уравнения найдены асимптотики первого и второго порядков его решения. Асимптотическое условие, а также замены асимптотик первого и второго порядков аналогичны заменам главе 2. Получены аналитические решения, определяющие средние характеристики Nm(x,t) и вторые моменты стохастической плотности численности заявок x,t).

Полученные аналитические решения полностью определяют параметры гауссовского распределения, которому удовлетворяет распределение вероятностей стохастической плотности численности заявок £,(х,/) возраста х в году t, то есть полностью решают поставленную задачу исследования автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания.

В главе 4 [61, 64, 82] предложено применение исследованной в главе 3 автономной немарковской системы массового обслуживания с неограниченным числом приборов и произвольным временем обслуживания к изучению процесса изменения демографической ситуации. Функционирование математической модели определяется двумя функциональными параметрами S(x) и b(x,t), которые имеют очевидную демографическую интерпретацию функции дожития женщин и функции фертильности [26]. Для применения полученных результатов к исследованию демографических процессов функции S(pc) и b(x,t) выбраны в следующем виде.

Функция дожития S(x) в виде модели Гомперца-Мейкема [102, 141, 155]

S{x) = ехр{-)(А + BeCu)du} = exp{-Jjt-В(еСх -1)/С}, о где параметр А учитывает риски, связанные с несчастными случаями, а слагае

Сх мое Be учитывает влияние возраста на смертность.

Интенсивность процесса рождаемости b(x,t) определим в виде b(x, t) = 0.488г|(?)\|/(х, 0, где т](У) — суммарный коэффициент рождаемости в момент времени t [26], а функция имеет смысл плотности распределения вероятностей значений репродуктивного возраста женщины. Функция y(x,t) задана как плотность смещенного двупарамерического у-распределения [97].

На основе сделанных предположений относительно динамики суммарного коэффициента рождаемости г|(?) в разделе 4.3 построены сценарии развития демографической ситуации в Российской Федерации на долгосрочную перспективу. А также найдены дисперсия численности возрастных групп, а также кова-риация между численностью возрастной группы нулевого возраста и численностью возрастной группы возраста х на период с 2005 по 2105 год.

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

Заключение

В представленной диссертационной работе предлагается математическая модель демографии в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов, а также ее частные случаи, а именно автономная система с РН-распределением времени обслуживания, трехфазная система массового обслуживания, пятифазная система массового обслуживания и система массового обслуживания с двумя типами заявок. Для исследования таких систем предложены метод моментов и метод асимптотического анализа и представлена возможность их применения к исследованию демографической ситуации.

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

В главе 2 предлагается математическая модель динамики численности женского населения в виде автономной системы массового обслуживания с неограниченным числом приборов и РН-распределением времени обслуживания, которая формирует новый класс систем массового обслуживания. Особенностью такой системы является, то, что в системе отсутствует внешний источник заявок, а продолжительность обслуживания т каждой заявки складывается из продолжительностей конечного числа фаз. В главе рассматриваются некоторые упрощенные варианты автономных систем массового обслуживания с РН-распределением времени обслуживания и методы их исследования, а именно, трехфазная система массового обслуживания (раздел 2.2), пятифазная система массового облуживания (раздел 2.3) и система массового обслуживания с двумя типами заявок (раздел 2.3).

Для указанных моделей составлены системы дифференциальных уравнений Колмогорова для распределения вероятностей числа обслуживаемых заявок и записаны дифференциальные уравнения для характеристических функций. Для решения дифференциальных уравнений для характеристических функций применены методы моментов и асимптотического анализа. Показано, что асимптотическое среднее значение Nm\(t) числа заявок, обслуживаемых в рассмотренных системах в момент времени t, найденные в предельном асимптотическом условии iV—>оо, совпадает с допредельным, полученным методом моментов. Заметим, как показали численные результаты, вторые моменты в предельной и допредельной моделях также совпадают.

В главе 3 предложена стохастическая демографическая модель с произвольным временем обслуживания в виде автономной немарковской системы массового обслуживания с неограниченным числом приборов. Принципиальное отличие от модели главы 2 заключается в том, что здесь исследуется стохастическая плотность численности заявок £(x,t) возраста х.

Предложен метод ее исследования, названный методом виртуальных фаз. Согласно которому, в разделе 3.2 исследована вспомогательная общая система массового обслуживания с РН-распределением времени обслуживания. Для распределения вероятностей числа обслуживаемых заявок записана система дифференциальных уравнений Колмогорова. В изучении характеристик вспомогательной системы с РН-распределением времени обслуживания применены методы асимптотического анализа и моментов.

Показано, что при ji—>оо рассматриваемое РН-распределение сходится к распределению, определяемому функцией S(x), где S(x) — функция дожития. Поэтому, сделав соответствующие замены, в разделе 3.4 перешли к рассмотрению автономной системы с произвольным временем обслуживания.

Для основного уравнения относительно характеристического функционала найдены асимптотики первого и второго порядков его решения. Получены аналитические решения, определяющие средние характеристики Nm(pc,t) и вторые моменты стохастической плотности численности заявок £(x,t) возраста jc в году t.

В главе 4 предложено применение исследованной автономной системы массового обслуживания с произвольным временем обслуживания к изучению процесса изменения демографической ситуации. Для применения полученных результатов к исследованию демографических процессов функции S(x) и b(x,t) выбраны в аналитическом виде. Функция дожития — в виде модели Гомперца-Мейкема, а плотность распределения вероятностей значений репродуктивного возраста женщины была задана как плотность смещенного двупарамерического у-распределения.

На основе сделанных предположений относительно динамики суммарного коэффициента рождаемости в разделе 4.3 построены сценарии изменения демографической ситуации на долгосрочную перспективу.

Таким образом, результаты работы имеют как теоретическое, так и практическое значение. Теоретическое значение заключается в дальнейшем развитии теории массового обслуживания, математической демографии и их аналитических методов, а именно

1. Сформирован новый класс систем массового обслуживания, включающий в себя автономную немарковскую систему массового обслуживания с произвольным временем обслуживания заявок и ее упрощенные варианты.

2. Разработаны методы исследования автономных немарковских систем массового обслуживания с неограниченным числом приборов, позволяющие находить распределение вероятностей численности заявок, обслуживаемых в системе, и основные характеристики, определяющие это распределение.

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

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

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

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

1. Аева Н. В. Статистический анализ демографического развития Республики Хакасия: дис. . канд. экон. наук: 08.00.12 / Н. В. Аева; Московский гос. унт экон., статистики и информатики. — М., 2001. 169 с.

2. Башарин Г. П. Математические модели прогнозирования демографических характеристик женского населения / Г. П. Башарин, Н. В. Новожилова // Вестн. Рос. ун-та дружбы народов. Сер.: Прикладная математика и информатика.- 1998.-№ l.-c. 48-64.

3. Белавин В. А. Математическая модель глобальных демографических процессов с учетом пространственного распределения / В. А. Белавин, С. П. Капица, С. П. Курдюмов // Журн. вычислительной математики и мат. физ. — 1998. — № 6. С. 885-902.

4. Библиографический бюллетень 1994 год Электронный ресурс. // Центр по изучению проблем народонаселения экон. фак. МГУ. Электрон, дан. — М., 1994. - URL: http://dmo.econ.msu.ru/demografia/Biblio/ CEDUM/b94.htm (дата обращения 24.05.2009).

5. Библиографический бюллетень 1995 год Электронный ресурс. // Центр по изучению проблем народонаселения экон. фак. МГУ. Электрон, дан. - М., 1994. - URL: http://dmo.econ.msu.ru/demografia/Biblio/ CEDUM/b95.htm (дата обращения 24.05.2009).

6. Библиографический бюллетень 1996 год Электронный ресурс. // Центр по изучению проблем народонаселения экон. фак. МГУ. Электрон, дан. - М., 1994. - URL: http://dmo.econ.msu.ru/demografia/Biblio/ CEDUM/b96.htm (дата обращения 24.05.2009).

7. Библиографический бюллетень 1998 год Электронный ресурс. // Центр по изучению проблем народонаселения экон. фак. МГУ. Электрон, дан. - М., 1994. - URL: http://dmo.econ.msu.ru/demografia/Biblio/ CEDUM/b98.htm (дата обращения 24.05.2009).

8. Библиографический бюллетень 1999 год Электронный ресурс. // Центр по изучению проблем народонаселения экон. фак. МГУ. — Электрон, дан. М., 1994. - URL: http://dmo.econ.msu.ru/demografia/Biblio/ CEDUM/b99.htm (дата обращения 24.05.2009).

9. Борисов В. А. Демография / В. А. Борисов. М.: Издат. дом NOTA BENE, 1999.-272 с.

10. Борисов В. А. Демография: учеб. для вузов / В. А. Борисов. М.: NOTA BENE Медиа Трейд Компания, 2005. - 344 с.

11. Боярский А. Я. Основы демографии / А. Я. Боярский, Д. И. Валентей, А. Я. Кваша; под ред. А. Я. Боярского. М.: Статистика, 1980. - 295 с.

12. Бреев Б. Д. Межотраслевой баланс движения населения и трудовых ресурсов (методологические вопр.) / Б. Д. Бреев, В. П. Крюков. — М.: Наука, 1974.- 183 с.

13. Бронштейн И. Н. Справочник по математике для инженеров и учащихся втузов / И. Н. Бронштейн, К. А. Семендяев. 13-е изд., испр. - М.: Наука, 1986.- 544 с.

14. Бутов В. И. Демография: учеб. пособие / под ред. В. Р. Игнатова. 2-е изд., перераб. и доп. - Ростов-н/Д.: МарТ, 2005. — 576 с.

15. Валентей Д. И. Основы демографии / Д. И. Валентей, А. Я. Кваша. М.: Мысль, 1989.-286 с.

16. Венецкий И. Г. Статистические методы в демографии / И. Г. Венецкий. -М.: Статистика, 1977. 207 с.

17. Винокурова Т. 3. Возрастная структура населения как фактор формирования трудового потенциала региона (на примере Республики Саха (Якутия)): дис. . канд. экон. наук: 08.00.05 / Т. 3. Винокурова; Ин-т регион, экон. Ан РС(Я) Якутск, 2001.-170 с.

18. Гнеденко Б. В. Введение в теорию массового обслуживания / Б. В. Гне-денко, И. Н. Коваленко. 3-е изд., испр. и доп. - М.: КомКнига, 2005. - 397 с.

19. Гольдин В. Я. Нестационарная математическая модель роста населения Земли / В. Я. Гольдин, В. Я. Пестрякова // Мат. моделирование. 1998. - Том 10, №3.-С. 39-47.

20. Гришагин В. А. Математическое моделирование стабильного населения: учеб. пособие. / В. А. Гришагин, Н. Р. Стронгина. Н. Новгород: Изд-во Ниже-город. гос. ун-та, 2002. - 49 с.

21. Демографические модели: сб. статей / под ред. Е. М. Андреева, А. Г. Волкова. -М.: Статистика, 1977. 182 с.

22. Демографические перспективы России: статистический сб. / ред. О. В. Осипова; Ин-т стат. и экономических исслед. (ИСЭИ) Госкомстата России. Отделение демографии. М.: Республиканский информ.-издат. центр , 1993. — 64 с.

23. Демографический энциклопедический словарь / под ред. Д. И. Валентея. -М.: Советская энциклопедия, 1985. 608 с.

24. Демография и статистика населения: учеб. / И. И. Елисеева, Э. К. Васильева, М. А. Клупт и др.; под ред. И. И. Елисеевой. — М.: Финансы и статистика, 2006. 688 с.

25. Демография: учеб. для втузов / под ред. Н. А. Волгина, JI. JI. Рыбаковско-го. М.: Логос, 2005. - 280 с.

26. Денисенко М. Б. Демография: учеб. пособие / М. Б. Денисенко, Н. М. Калмыкова. М.: Инфра-М, 2007. - 424 с.

27. Динамическая теория биологических популяций / под ред. Р. А. Полуэк-това. М.: Наука, 1974. - 455 с.

28. Джансеитов К. К. К 200-летию демографической теоремы Т. Р. Мальтуса // Вестн. КрасГУ, Сер. Физико-мат. науки. 2002. - Вып. 1. - С. 121-125.

29. Жафяров А. Ж. Об одной нелинейной модели движения стабильного населения // Оптимизация. — 1985. — Вып. 36. — С. 69-86.

30. Змеев О. А. Математический анализ. Часть I: учеб. пособие / О. А. Змеев, А. Ф. Терпугов, Р. Т. Якупов. Томск: Изд-во HTJI, 2006. - 176 с.

31. Капица С. П. Модель роста населения Земли и предвидимое будущее цивилизации Электронный ресурс. // Мир России. — 2002. — Т. 11, № 3. — Электрон. версия печат. публ. URL:http://www.ecsocman.edu.ru/mirros/msg/293878.html (дата обращения: 14.04.2009).

32. Капица С. П. Общая теория роста человечества: сколько людей жило, живет и будет жить на земле / С. П. Капица. М.: Наука, 1999. - 190 с.

33. Капица С.П. Рост населения Земли как главная глобальная проблема человечества // Глобальные проблемы биосферы. — М., 2001. С. 40-61.

34. Капица С. П. Рост населения Земли и предвидимое будущее цивилизации // Журн. теоретической и практической мед. 2003. — №1. - С. 73-86.

35. Капица С. П. Синергетика и прогнозы будущего / С. П. Капица, С. П. Курдюмов, Г. Г. Малинецкий. М.: Едиториал УРСС , 2003. - 283 с.

36. Капица С. П. Сколько людей жило, живет и будет жить на земле. Очерк теории роста человечества Электронный ресурс. / С. П. Капица. М., 1999.

37. Электрон, версия печат. публ. — URL:http://rogov.zwz.ru/Macroevolution/kapica.pdf (дата обращения: 30.04.2009).

38. Капица С. П. Феноменологическая теория роста населения Земли Электронный ресурс. // Успехи физ. наук. — 1996. № 1. - Электрон, версия печат. публ. - URL: http://ufii.ru/ru/articles/1996/1/с/ (дата обращения: 14.04.2009).

39. Карелова О. JI. Применение Марковских цепей для прогнозирования демографической ситуации в мире / О. JI. Карелова, М. А. Банько // Математическое моделирование. 2006. — Том 18, № 2. — С. 43-50.

40. Кваша А. Я. Что такое демография? / А. Я. Кваша. М.: Мысль, 1985. — 126 с.

41. Кетова К.В. Оптимальное распределение капиталовложений с учетом демографического прогноза: дис. . канд. физ.-мат. наук: 08.00.13 / К. В. Кетова; ИжГТУ. Ижевск, 2004 - 151 с.

42. Кильдишев Г. С. Статистика населения с основами демографии / Г. С. Кильдишев, JI. JI. Козлова, С. П. Ананьева. М.: Финансы и статистика, 1990. — 312 с.

43. Коровкин А. Г. Балансы движения населения и трудовых ресурсов и моделирование на их основе // Экономика и мат. методы. 1989. - Т. 25, Вып. 2. -С. 260-268.

44. Корчак-Чепурковский Ю. А. Избранные демографические исследования / Ю. А. Корчак-Чепурковский. М.: Статистика, 1970. - 386 с.

45. Крамер Г. Математические методы статистики: пер. с англ. / Г. Крамер; под ред. А. Н. Колмогорова. — 2-е изд., стер. М.: Мир, 1975. - 648 с.

46. Краснов М. JI. Интегральные уравнения / М. JI. Краснов, А. И. Киселев, Г. И. Макаренко. М.: Наука, 1976. - 216 с.

47. Курс демографии: учеб. пособие / А. Я. Боярский, Д. И. Валентей, А. Г. Вишневский и др.; под ред. проф. А. Я. Боярского. — 3-е изд., перераб. и доп. -М.: Финансы и статистика, 1985. 391 с.

48. Ломов С. А. Введение в общую теорию сингулярных возмущений / С. А. Ломов. М.: Наука, 1981.-400 с.

49. Мальтус Т. Р. Опыт о законе народонаселения: пер. с англ. / Т. Р. Мальтус. -Петрозаводск: Петроком, 1993. 139 с.

50. Математические методы в социологии: сб. науч. тр. / ред. Ф. М. Бород-кин. Новосибирск, 1974. - 215 с.

51. Модели демографических связей: сб. статей / под ред. А. Я. Боярского. -М.: Статистика, 1972. 128 с.

52. Медков В. М. Введение в демографию: учебное пособие для студентов вузов, обучающихся по экон. и социальным специальностям / В. М. Медков. -М.: Академический проспект, Фонд Мир, 2006. 428 с.

53. Медков В. М. Демография: учеб. / В. М. Медков. 2-е изд. - М.: ИНФРА-М, 2007. - 683 с.

54. Методические указания к решению задач по курсу "Дифференциальные уравнения" интегрирование линейных систем дифференциальных уравнений с постоянными коэффициентами / под ред. Л. А. Нежельской. — Томск: РИО ТГУ, 1996.-26 с.

55. Методы исследования: сб. статей / под ред. Л. Л. Рыбаковского. М.: Мысль, 1986. - 184 с.

56. Морозова А. С. Исследование основных характеристик математической модели сценария развития демографической ситуации / А. С. Морозова, М. Г.

57. Носова // Науч. творчество молодежи: Материалы IX Всероссийской науч.-практической конф. Томск: ТГУ, 2005. - Ч. 1. - С. 42-45.

58. Назаров А. А. Метод асимптотического анализа в теории массового обслуживания / А. А. Назаров, С. П. Моисеева. Томск: Изд-во HTJI, 2006. - 112 с.

59. Назаров А. А. Автономная немарковская система обслуживания с неограниченным числом приборов / А. А. Назаров, М. Г. Носова // ММПЭИТС: Материалы международной науч. конф. Минск: РИВШ, 2007. - С. 175-180.

60. Назаров А. А. Долгосрочный прогноз значений процесса изменения демографической ситуации / А. А. Назаров, М. Г. Носова // Науч. творчество молодежи: Материалы IX Всероссийской науч.-практической конф. Томск: ТГУ, 2005.-Ч. 1.-С. 37-40.

61. Назаров А. А. Исследование демографических процессов методами теории массового обслуживания / А. А. Назаров, М. Г. Носова // ММПЭИТС: Материалы международной науч. конф. Минск: БГУ, 2005. - С. 156-161.

62. Назаров А. А. Исследование системы с РН-распределением времени обслуживания / А. А. Назаров, М. Г. Носова // Науч. творчество молодежи: Материалы XI Всероссийской науч.-практической конф. Томск: ТТУ, 2007. — Ч. 1. С. 41-44.

63. Назаров А. А. Математическая модель процесса изменения демографической ситуации и ее исследование / А. А. Назаров, М. Г. Носова // Доклады Томского государственного университета систем управления и радиоэлектроники. — 2009. Т. 2 (20). - С. 100-105.

64. Назаров А. А. Метод передвижки возрастных групп в демографии и его приложения / А. А. Назаров, М. Г. Носова // Вестн. Томского гос. ун-та. Сер. Упр., вычислительная техника и информатика. — 2009. №3(8). — С. 67—75.

65. Назаров А. А. Стохастическая модель демографических процессов как автономная система массового обслуживания / А. А. Назаров, М. Г. Носова // Обозрение прикладной и промышленной математики. — 2009. — Т. 16, Вып. 6. — С. 1098-1099.

66. Назаров А. А. Теория вероятностей и случайных процессов: учеб. пособие / А. А. Назаров, А. Ф. Терпугов. Томск: Изд-во HTJI, 2006. - 204 с.

67. Назаров А. А. Теория массового обслуживания / А. А. Назаров, А. Ф. Терпугов. Томск: Изд-во НТЛ, 2004. - 228 с.

68. Население России 1999. Седьмой ежегодный демографический доклад / под ред. А. Г. Вишневского. М.: Книжный дом Ун-т, 2000. - 176 с.

69. Население России 2000. Восьмой ежегодный демографический доклад / под ред. А. Г. Вишневского. М.: Книжный дом Ун-т, 2001. — 176 с.

70. Население России 2001. Девятый ежегодный демографический доклад / под ред. А. Г. Вишневского. М.: Книжный дом Ун-т, 2002. — 216 с.

71. Население России 2002. Десятый ежегодный демографический доклад / под ред. А. Г. Вишневского. М.: Книжный дом Ун-т, 2004. — 224 с.

72. Новое в синергетике: Взгляд в третье тысячелетие: сб. статей / отв. ред. Г. Г. Малинецкий, С. П. Курдюмов. М.: Наука, 1999. - 479 с.

73. Орлов Ю. Н. Кинетические уравнения для нестационарных демографических моделей: препр. Электронный ресурс. // Ин-т прикладной математики РАН. Электрон, дан. - 2001. - URL:http://library.keldysh.ra/prepvw.asp?pid=1802 (дата обращения: 11.01.2009).

74. Орлов Ю. Н. Кинетические уравнения для некоторых моделей демографии / Ю. Н. Орлов Ю. Н., В. М. Суслин // Мат. моделирование. 2003. - Т. 15, № 3. - С. 43-54.

75. Перцев Н. В. Об одной модели динамики численности населения с учетом формирования и распада семейных пар // Мат. структуры и моделирование. — 2000.-№6.-С. 101-106.

76. Российский статистический ежегодник: стат. сб.: офиц. изд. // Госкомитет РФ по статистике / Под ред. B.JI. Соколин (пред.) и др. М.: 2001. - 642 с.

77. Русяк И. Г. К вопросу о выводе уравнения динамики возрастного состава / И. Г. Русяк, К. В. Кетова // Вестн. ИжГТУ. 2004. - № 2. - С. 49-52.

78. Русяк И. Г. К вопросу о математическом моделировании демографических показателей / И. Г. Русяк, К. В. Кетова, В. Б. Петряев // Высокие и информ. технологии в механике: Науч.—практическая конф., Ижевск, 2003. — Ижевск: Изд-во ИжГТУ, 2003. С. 163-169.

79. Русяк И. Г. Математическое моделирование демографических показателей / И. Г. Русяк, К. В. Кетова. Ижевск: Изд-во ИжГТУ, 2002. - С. 163-169.

80. Саати Т. JI. Элементы теории массового обслуживания и ее приложения / Т. JI. Саати. М.: Советское радио, 1971. - 519 с.

81. Соболева С. В. Демографические процессы в региональном социально-экономическом развитии / С. В. Соболева. Новосибирск: Наука, Сибирское отделение, 1988.-205 с.

82. Современная демография / под ред. А. Я. Кваши, В. А. Ионцевой. М: Изд-во МГУ, 1995. - 272 с.

83. Староверов О. В. Азы математической демографии / О. В. Староверов. -М.: Наука, 1997.-158 с.

84. Староверов О. В. Марковские модели мобильности населения // Экономика и мат. методы. 1983. - Т. 19. - Вып. №4. - С. 618-628.

85. Староверов О. В. Модели движения населения / О. В. Староверов. М.: Наука, 1979.-230 с.

86. Справочник по прикладной статистики / под ред. Э. Ллойда, У. Ледерма-на. М.: Финансы и статистика, 1990. - 508 с.

87. Терпугов А. Ф. Математическая статистика (конспект лекций) / А. Ф. Терпугов. Томск: Изд-во Томского университета, 1974. — 136 с.

88. Тихомиров Н. П. Демография. Методы анализа и прогнозирования: учеб. для вузов / Н. П. Тихомиров. М.: Экзамен, 2005. - 256 с.

89. Тихонова Л. Е. Государственное регулирование демографических процессов в Республике Беларусь: дис. . докт. экон. наук: 08.00.05 / Л. Е. Тихонова; М-во экон. Республики Беларусь, Научно-исследовательский экономический ин-т. М., 2002. - 296 с.

90. Уильямсон М. Анализ биологических популяций / М. Уильямсон; пер. с англ. А. Д. Базыкина; под ред. Ю. М. Свирежева. -М.: Мир, 1975. -271 с.

91. Фалин Г. И. Введение в актуарную математику / Г. И. Фалин, А. И. Фалин. М.: Изд-во МГУ, 1994. - 86 с.

92. Харченко JI. П. Демография: учеб. пособие для студентов вузов, обучающихся по специальности 080601 "Статистика" и др. экон. специальностям / JI. П. Харченко. - 2-е изд. - М.: Омега-JI, 2007. - 350 с.

93. Четыркин Е. М. Пенсионные фонды. Зарубежный опыт для отечественных предприятий, актуарные расчеты / Е. М. Четыркин. М.: Арго, 1993. - 100 с.

94. Шахотько JL П. Компьютерное решение задачи построения демографических прогнозов / JI. П. Шахотько, С. М. Терещенко // Вопросы статистики. — 1999. -№10. -С. 57-66.

95. Эльсгольц JI. Э. Дифференциальные уравнения и вариационное исчисление / JI. Э. Эльсгольц. М.: Наука, 1969. - 424 с.

96. Ahlburg D., Vaupel J. Alternative projections of the U.S. population Электронный ресурс. // Demography. — 1990. — Vol. 27, N 4. — Электрон, версия печат. публ. Доступ из базы данных "Jstor". — URL: http://www.jstor.org (дата обращения: 11.01.2009).

97. Alho J. Statistical demography and forecasting New York Электронный ресурс. / J. Alho, B. Spencer. Springer, 2005. - Электрон, версия печат. публ. -Доступ из базы данных "Springerlink". - URL: http://www.springerlink.com (дата обращения: 24.05.2009).

98. Caswell Н. Matrix population models: construction, analysis, and interpretation / H. Caswell. Sunderland, Massachusetts, Sinauer Associates, 1989. - 722 p.

99. Caswell H., Trevisan M. Sensitivity analysis of periodic matrix models Электронный ресурс. // Ecology. 1994. - Vol. 75, N 5. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". — URL: http://www.jstor.org (дата обращения: 11.01.2009).

100. Donald Т. R. Demographic methods and concepts / Т. R. Donald. Oxford: Oxford University Press, 2006. - 523 p.

101. Goodman L. Population growth of the sexes Электронный ресурс. // Biome-trika. 1953. - Vol. 9, N 2. - Электрон, версия печат. публ. — Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 23.04.2009).

102. Hinde A. Demographic methods / A. Hinde. Arnold, 1998. - 305 p.

103. Keyfitz N. Applied mathematical demography Электронный ресурс. / N. Keyfitz, Н. Caswell. — third edition. Электрон, версия печат. публ. - Springer

104. New York, 2005. Доступ из базы данных "Springerlink". - URL: http://www.springerlink.com (дата обращения: 28.04.2009).

105. Keyfitz N. Introduction to the mathematics of population with revisions / N. Keyfitz. Wesley Pub. Co, 1977. - 490 p.

106. Keyfitz N. Models Электронный ресурс. // Demography. — 1971. — Vol. 8, N 4. Электрон, версия печат. публ. — Доступ из базы данных "Jstor". — URL: http://www.jstor.org (дата обращения: 11.01.2009).

107. Keyfitz N. Murphy matrix and multiple decrement in population analysis Электронный ресурс. // Biometrics. 1967. - Vol. 23, N 3. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 01.05.2009).

108. Keyfitz N. On future population Электронный ресурс. // Journal of the American Statistical Association. 1972. - Vol. 67, N 338. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 11.01.2009).

109. Keyfitz N. On the momentum of population growth Электронный ресурс. // Demography. 1971. - Vol. 8, N 1. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 11.01.2009).

110. Namboodiri К. Demographic analysis: a stochastic approach / К. Namboodiri.- Academic Press, 1990. 390 p.

111. Newell C. Methods and models in demography / C. Newell. Belhaven, 1988. -217 p.

112. Olshansky S., Carnes B. Ever Since Gompertz Электронный ресурс. // Demography. 1997. - Vol. 35, № 1. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 02.03.2009).

113. Preston S., Coale A. Age structure, growth, attrition, and accession: A new synthesis Электронный ресурс. // Population Index. 1982. - Vol. 48, N 2.

114. Электрон, версия печат. публ. — Доступ из базы данных "Jstor". URL: http://www.jstor.org (дата обращения: 11.01.2009).

115. Preston S. Demography: measuring and modeling population processes / S. Preston, P. Heuveline, M. Guillot. Wiley-Blackwell, 2001. - 291 p.

116. Shryock Н. S. Forecasts of population in the United States Электронный ресурс. // Population Studies. 1950. - Vol. 3, № 4. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". - URL: http://www.jstor.org (дата обращения: 24.05.2009).

117. Sweeney S. Н., Konty К. J. Population forecasting with nonstationary multiregional growth matrices // Geogr.Anal. 2002. - N 4. - P. 289-312.

118. Sykes Z., Kim J. Dynamics of some special populations with NRR = 1 Электронный ресурс. // Demography. 1978. - Vol. 15, N 4. - Электрон, версия печат. публ. - Доступ из базы данных "Jstor". — URL: http://www.jstor.org (дата обращения: 11.01.2009).

119. Эмпирические значения функции дожития