Построение адаптивных моделей прогнозов динамики логистических потоков на основе фильтров Брауна
Важную роль в процессе проектирования и управления функционированием РЛС ПП играют прогнозные оценки функционирования РЛС ПП, которые можно получать с помощью разных методов, для выбора которых важно соотнести имеющееся количество и вид данных с требованиями к минимальному количеству исходных данных, на основе которых можно сделать прогноз тем или иным методом [109, 132].
Наиболее часто для прогнозирования используются методы, основанные на анализе временных рядов. Эти методы имеют разную сложность.
К простым методам сглаживания данных можно отнести метод экспоненциального сглаживания с одним параметром и метод арифметического сглаживания. Важнейшая предпосылка любого метода сглаживания состоит в использовании последних данных ряда, поскольку информация имеет свойство устаревания.
Часто используются модели экспоненциального сглаживания с одним, двумя параметрами. В большинстве случаев они линейны. Исходя из этого предлагается полиномиальная модель второго порядка, параметры которой рассчитываются по методике Брауна [63].
Данная методика позволяет рассчитать параметры модели по достаточно простым рекуррентным формулам, через функцию сглаживания второго порядка, по-разному учесть влияние наблюдений временного ряда. Варьируя параметрами сглаживания, можно адаптировать прогнозные модели к различным условиям применения (различным видам временных рядов).
Рассмотрим модель построения прогноза работы логистической системы пассажирских перевозок на основе фильтров Брауна[51,55]. Для анализа используем полиномиальную модель.
- неизвестные параметры модели;
q- степень полинома.
В расчетах воспользуемся квадратичной модлеью, основным достинством которой явялется простота:
C- параметры модели
t- текущее время.
Оценки
параметров модели (3.2.2) могут быть найдены по
МНК. Этот метод наиболее разработан в теоретическом отношении и позволяет получить несмещенные, состоятельные, а в случае независимого остатка и эффективные оценки. Однако данный метод не единственно возможный.
Обычно наибольшее воздействие на будущий период оказывают более «поздние» наблюдения, и поэтому им следует придавать больший вес по сравнению с «ранними». Именно этим принципам отвечает разработанный Брауном [2, 3] метод экспоненциального сглаживания.
В методе Брауна взвешенная скользящая средняя с экспоненциально распределенными весами характеризует значение процесса на конце интервала сглаживания. Это позволяет для исходного ряда получить оценки параметров, определяющих неслучайную функцию (тренд процесса) и характеризующих не средний уровень процесса, а тенденцию, которая сложилась к моменту последнего наблюдения.
Учитывая вышеизложенное, в работе предлагается методика прогноза и оценки параметров фильтров Брауна для краткосрочного прогнозирования.
По Брауну, предсказуемое значение функции P(t)может быть выражено рядом Тейлора
члены которого связаны рекуррентной формулой
определяющей экспоненциальную среднюю n-го порядка для ряда P(t). Параметр сглаживания
учитывает влияние исходного ряда
наблюдений на результаты прогнозирования.
Для фильтра Брауна второго порядка имеют место следующие соотношения:
В качестве начального используется условие
где
- значение, соответствующее точке, с которой начинается прогноз.
Качество прогноза, оцениваемое среднеквадратической и максимальной относительной ошибками, существенно зависит от выбора параметра сглаживания а. В зависимости от величины этого параметра прогнозные оценки по-разному учитывают влияние исходного ряда наблюдений: чем больше а, тем существеннее вклад последних наблюдений в результат прогнозирования (влияние начальных условий быстро убывает) [5].
Помимо параметра а, на качество прогноза оказывает влияние и выбор величины предыстории к, используемой для прогноза.
Величина ошибки
выражена в процентах от
фактических значений нагрузки для различных значений параметра a∈[0,1;0,9].
Таким образом, рассматриваемая модель прогноза является адаптируемой по параметрам к различным интервалам упреждения.
Разброс прогнозных значений
процесса P(t)для фильтра Брауна характеризуется величиной дисперсии
где
- разность между процессом изменения пассажиропотока и
средней составляющей.
Дисперсия прогнозного значения пассажиропотока позволяет определить доверительную верхнюю Рв и нижнюю Рн границы прогноза
где tβ - критическая точка распределения Стьюдента для k-2степеней свободы и уровня значимости β =0,05.
В условиях поставленной задачи для того, чтобы свести к минимуму выход пассажиропотока за установленное ограничение, достаточно воспользоваться верхней границей PB(t)прогнозного значения пассажиропотока, определяемой соотношением (9).
Процедура коррекции выборок, используемых для прогнозирования пассажиропотока.
Прогнозируемый процесс (перевозка пассажиров) на отрезке коррекции
- объем предыстории, по которой выполняется прогноз (строится модель)) представляется в виде соотношения
где
откорректированный процесс,
- суммарная объем перевозок, величина которого изменилась на отрезке коррекции;
- нагрузка /-го маршрута и j-го уровня
транспортной сети, изменившего режим работы;
- характер
изменения режима работы /-го маршрута и j-го уровня сети: -1 - означает, что маршрут не работает; +1 - функционирует в полном объеме.
После реализации процедуры коррекции выполняются расчет параметров
математической модели, вычисление прогнозного значения объема перевозок и его верхней доверительной границы.
Обеспечение качественного анализа функционирования региональных логистических систем основано на среднесрочных и долгосрочных прогнозах объемов перевозок пассажиров, грузов и платного багажа по видам транспорта, включенных в логистическую систему. Рассмотрим данное обстоятельство на примере развития региональной логистической системы Саратовской области. В качестве основы выступают перспективные планы развития, сформированные министерством транспорта области, а также показатели статистической отчетности развития экономики на 20082014 годы и тенденций до 2038 г.
Проведение прогнозных расчетов осуществлялось на основе углубленного анализа перечисленных документов. Учитывая, что Саратовская область до настоящего времени является дотационным регионом Российской Федерации, в основу прогноза объемов авиационных перевозок заложена мультипликативная регрессионная модель вида
где Yt- объемы пассажирских, либо грузовых перевозок из аэропорта в t - м году;
- параметры модели;
J1t- индекс соотношения средней заработной платы к бюджетному прожиточному минимуму;
J 2t- индекс авиационных тарифов в t-м году по отношению к среднему уровню тарифов в базовом году.
Приведенная регрессионная модель соответствует так называемому пессимистическому сценарию развития экономики Саратовской области.
Для построения оптимистического сценария прогноза используется регрессионное уравнение вида
где Y1- перевозки на воздушном транспорте в расчете на одного жителя;
Y2- перевозки на железнодорожном транспорте в межгородском сообщении в расчете на одного жителя;
C1- средний тариф на воздушном транспорте;
C2- средний тариф на железнодорожном транспорте;
B- расходы населения на транспортные услуги воздушного и железнодорожного транспорта в расчете на одного жителя;
a- параметры модели.
Прогнозные расчеты пассажирских авиаперевозок, приведенные в соответствии с моделью (3.2.9), показывают, что в 2038 году перевозки пассажиров составят 1 035 тыс. человек, индекс роста по сравнению с базовым 2007 годом составит 5,1, среднегодовой индекс роста составит 1 053, что ниже прогнозируемого среднегодового индекса роста ВРП.
Прогнозные расчеты пассажирских перевозок, приведенные в соответствии с моделью (3.2.10), показывают, что в 2038 году перевозки составят 2 930 тыс. человек, индекс роста по сравнению с 2007 г. - 14,45, среднегодовой индекс роста - 1,09, что выше прогнозируемого
среднегодового индекса роста ВРП.
Из-за невозможности принимать средне- и дальнемагистральные воздушные суда, в качестве исходных условий также принимается, что в аэропорту г. Саратова существует неудовлетворенный спрос на перевозки.
Для расчета неудовлетворенного спроса в качестве опорного эквивалента используется аэропорт города Самары. При расчете неудовлетворенного спроса были применены понижающие коэффициенты, учитывающие разницу в численности населения в зонах тяготения аэропортов и уровне экономического развития регионов.
94
Таким образом, данные варианты расчетов позволят разрабатывать рациональные варианты обслуживания формирующихся пассажиропотоков с учетом перспективных вариантов развития всей региональной логистической системы пассажирского транспорта.
Таблица 3.2.1 - Прогноз перевозок пассажиров, грузов и почты через аэропорт г. Саратова за 2008-2038 годы (по данным министерства транспорта Саратовской области)
| Показатели | Прогноз (годы) | ||||||
| 2015 | 2016 | 2017 | 2018 | 2019 | 2020 | 2021 | |
| Пассажирские перевозки, всего (чел.), в т. ч. | 564 753 | 704 371 | 802 366 | 872 974 | 939 320 | 1 004 039 | 1 070 306 |
| Отправлено (чел.) | 278 986 | 347 957 | 396 366 | 431 246 | 464 021 | 495 992 | 528 727 |
| в том числе | |||||||
| по внутренним ВЛ | 244 174 | 304 539 | 346 907 | 377 435 | 406 120 | 434 102 | 462 752 |
| по международным ВЛ | 34 812 | 43 418 | 49 459 | 53 811 | 57 901 | 61 890 | 65 975 |
| Принято (чел.) | 281 739 | 351 390 | 400 277 | 435 501 | 468 600 | 500 886 | 533 945 |
| в том числе | |||||||
| по внутренним ВЛ | 247 548 | 308 746 | 351 700 | 382 649 | 411 732 | 440 099 | 469 146 |
| по международным ВЛ | 34 191 | 42 644 | 48 577 | 52 852 | 56 868 | 60 787 | 64 799 |
| Прямой транзит (чел.) | 4 028 | 5 024 | 5 723 | 6 227 | 6 699 | 7 161 | 7 634 |
| Грузовые перевозки, всего (тонн), в т. ч. | 1 632,6 | 2 036,8 | 2 320,7 | 2 525,3 | 2 717,6 | 2 905,3 | 3 097,5 |
| Отгружено (тонн) | 308 | 386 | 442 | 483 | 522 | 560 | 600 |
| в том числе | |||||||
| по внутренним ВЛ | 281,8 | 353,2 | 404,5 | 442,0 | 477,7 | 512,4 | 549,0 |
| по международным ВЛ | 26,2 | 32,8 | 37,5 | 41,0 | 44,3 | 47,6 | 51,0 |
| Разгружено (тонн) | 1 324,6 | 1 650,8 | 1 878,7 | 2 042,3 | 2 195,6 | 2 345,3 | 2 497,5 |
| в том числе | |||||||
| по внутренним ВЛ | 1 212,1 | 1 510,6 | 1 719,1 | 1 868,8 | 2 009,1 | 2 146,2 | 2 285,4 |
| по международным ВЛ | 112,5 | 140,2 | 159,6 | 173,5 | 186,5 | 199,1 | 212,1 |
Таблица 3.2.2 - Прогноз вылетов воздушных судов по аэропорту г.
Саратова в период 2016-2023 годов (по данным министерства транспорта
Саратовской области)
| № п/п | Тип ВС | 2016 | 2017 | 2018 | 2019 | 2020 | 2021 | 2022 | 2023 |
| Ближнемагистральные ВС | |||||||||
| 1 | Як-42 | 1892 | 1785 | 1612 | 1499 | 1187 | 1020 | 875 | 764 |
| 2 | Як-40 | ||||||||
| 3 | Ан-24 | ||||||||
| 4 | Ан-74 (и модификации) | 10 | 11 | 12 | 13 | 13 | 14 | 15 | 15 |
| 5 | Ан140 (и модификации) | 216 | 364 | 468 | 555 | 628 | 733 | 813 | 890 |
| 6 | SRJ-100-75 | 757 | 862 | 934 | 1000 | 1063 | 1125 | 1189 | 1249 |
| 7 | SRJ-100-95 | 2643 | 3402 | 4043 | 4576 | 5327 | 5906 | 6468 | 6977 |
| Среднемагистральные ВС | |||||||||
| 8 | Ту-154М | 106 | 128 | 111 | 103 | 88 | 76 | 61 | 50 |
| 9 | Ту-204 (и модификации) | 181 | 237 | 259 | 279 | 302 | 325 | 348 | 372 |
| 10 | A319 (и модификации) | 285 | 372 | 404 | 437 | 471 | 508 | 544 | 582 |
| 11 | B737-300 (700) | 271 | 355 | 388 | 418 | 453 | 487 | 522 | 558 |
| 12 | B737-500 (600) | 317 | 415 | 452 | 488 | 528 | 568 | 610 | 651 |
| Дальнемагистральные ВС | |||||||||
| 13 | самолеты класса A-310-300 | ||||||||
| Грузовые ВС | |||||||||
| 14 | Ан-26 | ||||||||
| 15 | Ан-72 | 27 | 25 | 22 | 20 | 18 | 16 | 14 | 12 |
| 16 | Ан-74Т (и модификации) | 102 | 118 | 130 | 141 | 150 | 160 | 170 | 180 |
| Вертолеты | |||||||||
| 17 | Ми-8 (и замена Ми-8) | 373 | 391 | 408 | 417 | 428 | 440 | 452 | 467 |
| Малая и бизнес-авиация | |||||||||
| 18 | Н25В, F-900, FA-50 и т.п. | 163 | 192 | 207 | 225 | 240 | 255 | 273 | 288 |
Представленные перспективные показатели работы аэропорта г. Саратова показывают уровень загрузки терминалов. Примененные модели прогнозирования учитывают многофакторный подход при оценке показателей, однако необходимо отметить, что полученные значения рассчитаны и на ближайшую перспективу, но сложность оценки полученных пессимистических и оптимистических прогнозов заключается в том, что аэропорт не введен в эксплуатацию. Данное обстоятельство подталкивает нас
к размышлениям о необходимости корректировки полученных прогнозов и учету тех факторов, которые не введены в данную модель, особенно сложно формализуемых и оцененных количественно. Именно влияние таких факторов как политическая обстановка, экономический кризис необходимо учитывать в первую очередь, так как они в значительной степени оказали влияние на развитие РЛС 1П1, а именно в Саратовской области. Указанные направления исследования в дополнение к уже изложенным лишь подтверждают правильность выбранного курса и необходимость разработки интеллектуальных моделей управления РЛС ПП.
Точкой отсчета должны стать данные о пассажиропотоках на функционирующем аэровокзале, и уже опираясь на них, необходимо выстраивать новые прогнозы. В дополнение к этому в расчетах нами предложено использовать полиномиальные модели с подбором коэффициентов модели с использованием фильтров Брауна.
Для обеспечения качественного прогноза была разработана программа «Прогнозирование объемов перевозок в системах транспорта и логистики» [79], которая позволяет осуществлять:
- загрузку ретроспективных данных.
- визуализацию исходных данных в виде графика и таблицы;
- установку и изменение длины предыстории;
- установку и изменение параметров прогнозирующей модели:
• порядка прогнозной модели;
• параметра сглаживания, для адаптации модели к различным ретроспективным данным;
• интервала упреждения, на который может быть осуществлен прогноз;
- визуализацию результатов прогнозирования в виде графика и таблиц;
- расчет и отображение относительных ошибок прогнозирования на различных интервалах прогноза;
- выделение максимальной относительной ошибки прогнозирования и временного интервала ее возникновения;
- подстройку параметров прогнозной модели.
Создание программного продукта для построения адаптивных моделей прогнозов работы систем транспорта и логистики с помощью методики Брауна.
Изложены теоретические положения построения прогнозных моделей на основе полиномиальной модели второго порядка с адаптацией параметров с использованием фильтров Брауна. Предполагается процедура коррекции выборок, используемых для прогнозирования. Практическая реализация предложенной модели прогнозирования для систем логистики и транспорта позволяет повысить качество принятия управленческих решений.
Таблица 3.2.3 - Объемы перевозок пассажиров в г. Саратове (по данным МКУ «Транспортное управление» администрации г. Саратова)
| Год | Объем перевозок пассажиров, |
| тыс. пасс. | |
| 2003 | 66564 |
| 2004 | 76318 |
| 2005 | 93867,3 |
| 2006 | 100610 |
| 2007 | 104038 |
| 2008 | 106916,7 |
| 2009 | 107163 |
| 2010 | 109662 |
| 2011 | 109477 |
| 2012 | 111536 |
| 2013 | 114723 |
| 2014 | 112385 |
Результаты прогнозирования при помощи данной программы
Рисунок 3.2.2 - Расчет величины относительной ошибки в программе «Прогнозирование объемов перевозок в системах транспорта и логистики»
Рисунок 3.2.3 - Фрагмент интерфейса программы
Процедура построения прогностической модели выполнена и проверена на показателях области о работе пассажирского транспорта за предыдущие годы по Саратовской области. При подборе статистической информации мы столкнулись с проблемой возможности использования ее данных для построения прогноза, так как до 1999 года была иная система учета статистических данных. В связи с этим выполним построение прогноза начиная с 2003 года. Длина предыстории составила 10 лет [72].
Величина полученного прогноза объема перевозок пассажиров зависит от выбора параметра сглаживания а. В зависимости от величины этого параметра прогнозные оценки по-разному учитывают влияние исходного ряда наблюдений: чем больше а, тем существеннее вклад последних наблюдений в результат прогнозирования. В таблице приведены зависимости изменения параметра сглаживания и относительной ошибки прогноза.
По данным таблицы видно, что минимальная ошибка прогноза получена при α, равном 0,2. Изменения величины параметра сглаживания не оправданны большими значениями относительной ошибки прогноза и получением некачественного прогноза.
3.3