Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,593

ИССЛЕДОВАНИЕ АППРОКСИМАЦИИ ФУНКЦИЙ РАПРЕДЕЛЕНИЯ С ПРИМЕНЕНИЕМ АППАРАТА «СШИТЫХ» ФУНКЦИЙ

Киреева Н.В. 1 Чупахина Л.Р. 1
1 ФГОБУ ВПО «Поволжский государственный университет телекоммуникаций и информатики»
Произведен анализ функций распределения Вейбулла при аппроксимации ее плотности: суммой затухающих экспонент и аппроксимации на двух участках – участок около нуля аппроксимации с помощью полиномов, и на участке «тяжелого» хвоста распределения – с помощью суммы затухающих экспонент. Данное решение показало, что наилучшая аппроксимация достигается при втором случаи. Для объединения двух функций использован алгоритм «сшивания» функций, и далее построен график одной их характеристик сети среднее время ожидания пакетов в очереди после решения спектральным способом интегрального уравнения Линдли.
аппроксимация
экспонента
уравнение Линдли
распределение с «тяжелым» хвостом
«сшитые» функции
1. Мищенко В.А. Метод селектирующих функций в нелинейных задачах контроля и управления. – Москва: Сов. радио, 1973. – 184 с.
2. Метод аппроксимации произвольной плотности распределения суммами экспонент / Блатов И.А., Карташевский В.Г., Киреева Н.В., Чупахина Л.Р. // Вестник ВГУ. – 2013. – № 2. – С. 53–57.
3. Решение уравнения Линдли спектральным методом для систем массового обслуживания общего вида / Блатов И.А., Карташевский В.Г., Киреева Н.В., Чупахина Л.Р. // Электросвязь. – 2014. – № 11. – С. 48–50.

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

В работе представлено исследование разработанного алгоритма решения интегрального уравнения (ИУ) Линдли спектральным методом для системы массового облуживания (СМО) типа G/G/1, на основе аппроксимаций с помощью суммы затухающих экспонент. Данная аппроксимация позволяет получить выражение для функции распределения времени ожидания и для среднего времени ожидания пакета в очереди.

Однако при проведении анализа аппроксимации суммой затухающих экспонент возникла проблема в минимизации и сведения к нулю погрешности аппроксимации. Полученные результаты исследования, на примере распределений «тяжелого» хвоста (РТХ), доказывают, что аппроксимация произвольной функции плотности распределения вероятностей, суммой затухающих экспонент целесообразна, так как полученные теоретические априорные оценки погрешности минимальны. Но, на примере функции распределения Вейбулла видно, что осталась не решенной проблема на начальном участке аппроксимации, так как аппроксимация суммой затухающих экспонент плохо описывает исследуемую функцию на участке около нуля.

Задача

Поэтому была предпринята попытка исследования двух участков функции распределения, так называемых «сшитых» функций распределения.

Есть функция Вейбулла f(x), которую аппроксимируем по участкам от (0, x0) и до (x0, ∞), для более точного решения ИУ Линдли, чтобы свести погрешность аппроксимации R(x) → 0.

Плотность функции распределения Вейбулла имеет следующий вид (рис. 1):

kir01.wmf (1)

kireev1.tif

Рис. 1. Функция плотности распределения Вейбулла при α = 1,6 и β = 0,5

Как видно из рис. 2 около нуля возникла проблема с плохой аппроксимацией суммой затухающих экспонент, поэтому будем аппроксимировать функции отдельно и затем «сшивать» их [1].

kireev2.tif

Рис. 2. Сравнение двух плотностей распределения Вейбулла

Представляем аппроксимацию f(x) на участке (0, x0) полиномами:

kir02.wmf

Рассматриваем аппроксимацию f(x) на участке (x0, ∞) суммой затухающих экспонент [2, 3]:

kir03.wmf

Обязательно необходимо учесть условие нормировки для обеих частей «сшитой» функции.

Рассмотрим «сшитие» функций в одну [1]:

kir05.wmf

kir06.wmf (2)

kir07.wmf

kir08.wmf

Учитываем условие нормировки

kir09.wmf

kireev3.tif

Рис. 3. Аппроксимация функции распределения Вейбулла f(x) на участке (0, x0) полиномами

Необходимо «сшить» две функции φ1(x) и φ2(x) в одну φc(x) а затем найти преобразование Лапласа для решения ИУ Линдли спектральным методом [2, 3].

kir10.wmf (3)

где A(s) и B(s) – преобразование Лапласа плотности распределения промежутков времени между поступлениями пакетов и плотности распределения времени обслуживания, соответственно.

Аналитический вид «сшитой» функции:

kir11.wmf

В нашем случае получаем:

kir12.wmf

kir13.wmf

так как функция Вейбулла f(x) существует на отрезке от (0, x0).

В нашем случае примем a = 0, b = ∞ и получим интеграл вида (4)

kir14.wmf kir15.wmf (4)

Для функции Вейбулла f(x), аппроксимированной полиномами на участке (0, x0), имеем (5):

kir16.wmf, (5)

где kir17.wmf.

Для функции Вейбулла f(x), аппроксимированной суммой затухающих экспонент на участке (x0, ∞), имеем (6):

kir18.wmf (6)

где kir19.wmf.

Пример «сшитых» функций

Рассмотрим случай, когда промежутки времени между поступлениями пакетов имеют распределение Вейбулла с параметрами α = 1,6 и β = 0,5, а промежутки времени обслуживания имеют распределение Парето с параметрами α = 1,5 и β = 1. Экспериментальным путем определяем точку x0, при которой достигается максимальной значение функции распределения Вейбулла при данных параметрах, в этом случае x0 = 0,3.

В соответствии с алгоритмом аппроксимации суммой затухающих экспонент, после получение преобразования Лапласа I «сшитой» функции, необходимо решить ИУ Линдли спектральным методом.

A(s) – преобразование Лапласа плотности распределения промежутков времени между поступлениями пакетов от функции Вейбулла. Преобразование Лапласа от сшитой функции это I = A(s), B(s) – преобразование Лапласа плотности распределения времени обслуживания от функции Парето:

kir20.wmf

Решаем ИУ Линдли спектральным методом [2] получаем функцию времени ожидания пакетов в очереди W(t):

kir21.wmf

kir22.wmf

kireev4.tif

Рис. 4. Аппроксимация f(x) на участке (x0, ∞) суммой затухающих экспонент

kireev5.tif

Рис. 5. Функция времени ожидания пакетов в очереди W(t)

Выводы

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


Библиографическая ссылка

Киреева Н.В., Чупахина Л.Р. ИССЛЕДОВАНИЕ АППРОКСИМАЦИИ ФУНКЦИЙ РАПРЕДЕЛЕНИЯ С ПРИМЕНЕНИЕМ АППАРАТА «СШИТЫХ» ФУНКЦИЙ // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 5-1. – С. 17-20;
URL: https://applied-research.ru/ru/article/view?id=9173 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674