Размещено на http://www.allbest.ru/

1

Научный журнал КубГАУ, №92(08), 2013 года

http://ej.kubagro.ru/2013/08/pdf/76.pdf

ИСПОЛЬЗОВАНИЕ АППРОКСИМАЦИИ РОЗЕНБЛАТТА-ПАРЗЕНА ДЛЯ ВОССТАНОВЛЕНИЯ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ НЕПРЕРЫВНОЙ СЛУЧАЙНОЙ ВЕЛИЧИНЫ С ОГРАНИЧЕННЫМ ОДНОМОДАЛЬНЫМ ЗАКОНОМ РАСПРЕДЕЛЕНИЯ

Поршнев Сергей Владимирович, д.т.н., профессор

Копосов Александр Сергеевич, аспирант

Уральский федеральный университет имени первого Президента России Б.Н. Ельцина, Екатеринбург

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

Ключевые слова: НЕПАРАМЕТРИЧЕСКАЯ СТАТИСТИКА, ВОССТАНОВЛЕНИЕ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ, АППРОКСИМАЦИЯ РОЗЕНБЛАТТА-ПАРЗЕНА, ПАРАМЕТР РАЗМЫТОСТИ, МЕТОД МНИМЫХ ИСТОЧНИКОВ

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

(1)

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

(2)

где функция Хэвисайда

в соответствие с выбранной мерой близости, зависящей, вообще говоря, от вида распределения [3].

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

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

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

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

(3)

где - ядерная функция, удовлетворяющая следующим условиям:

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

б) функция, симметричная относительно 0;

в) при

- параметр «размытости», определяющий гладкость получаемой оценки.

Соответственно, плотности вероятности вычисляется по формуле

(4)

Оценки плотности вероятности (4) предложены Розенблаттом [7] и исследованы Парзеном [8]. На практике наиболее часто в качестве ядерных функций используются функции, представленные в таблице 1 [2].

Таблица 1. Ядерные функции, наиболее часто используемые на практике

аппроксимация розенблатт парзен распределение

Оптимальные значения ядерной функции и параметра h находятся из условия достижения функционалом

(5)

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

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

(6)

и далее выбирать ту ядерную функцию , для которой значение функции

 (7)

в точке будет наибольшим.

При этом качество аппроксимации функции распределения в рассматриваемом методе оказывается зависящим от вида ядерной функции и значения параметра размытости [6], что иллюстрирует рис. 1.

Рисунок 1 Плотность распределения случайной последовательности , сгенерированной в соответствие с нормальным законом распределения :1 гистограмма случайной последовательности; 2 нормальное ядро, , где пока не определено; 3 нормальное ядро, , где пока не определено; 4 нормальное ядро, , где пока не определено.

Отметим, что формально метод Розенблатта-Парзена позволяет, построить аппроксимацию функции распределения любой конечной случайной последовательности, которая при условии правильного выбора параметра h, оказывается достаточно гладкой (см. [6]). При оценка (3) является непараметрической и асимптотически оптимальной оценкой функции распределения в классе непрерывных функций [6]. Однако на практике случайные последовательности для которых ищется функция распределения, имеют конечную длину. В этой ситуации метод Розенблатта-Парзена, вообще говоря, не гарантирует близости аппроксимирующей функции, найденной с его помощью, и истинной функции распределения. В этой связи представляет практический интерес провести анализ точности восстановления методом Розенблатта-Парзена функции распределения случайных последовательностей конечной длины, теоретический закон распределения которых известен.

В статье изложены результаты анализа точности восстановления функции распределения случайной величины с ограниченной областью рассеяния, основные свойства которой описаны в [3], методом Розенблатта-Парзена [7,8] и методом мнимых источников, в котором параметры аппроксимирующей функции находятся с помощью генетических алгоритмов [9].

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

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