Ортогональные полиномы и кривые распределения вероятностей

полиномы и кривые распределения вероятностей" width="117" height="68" />,

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

Выражения для  будет иметь вид:

.

Выражения для коэффициентов  будут следующими:

.

Вводя для сокращения обозначение

через , запишем выражение для  в таком виде:

.

Для  выражение будет иметь вид

.

Что касается величин  и , то они равны соответственно

 и .

Теперь перейдем к определению коэффициентов  в выражении

.

Для получим выражение

.

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

.

Также упростятся выражения для

 и .

Функция  станет равной , функции  определяются путем последовательных подстановок выражений  в формулы

.

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

.

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

Обозначим сумму квадратов отклонений через . Тогда можно написать

.

будет равняться

,

а выражать рекуррентно через по формуле

.

Итак,

, , ,

, , , ,

, , , , .

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

§ 2. Обобщение Грамма - Шарлье.

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

 (4)

где  - есть к–ая производная функции . Здесь полагаем, что

.

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

Производные функции  мы можем представить в виде [3]

,

тогда можем записать

где функции  должны удовлетворять следующему свойству:

       если     (5)

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

=

=

Отсюда следует, что

.

На практике в этом разложении мы используем только четыре первых члена, и коэффициенты перед ними есть:

          

          

Коэффициенты  имеют четкий статистический смысл, а именно: коэффициент , выраженный через , отвечает за асимметрию закона распределения, и коэффициент  выраженный через  - за эксцесс или дефект кривой распределения.

Свойство (5) есть свойство ортогональности полиномов, т. е.  по определению является системой ортогональных полиномов, которая получена по способу Чебышева в предыдущем параграфе [3], [5].

§ 3. Весовые функции и системы ортогональных полиномов.

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

Полиномы Чебышева - Эрмита.

Пусть многочлен (2) не имеет корней, тогда уравнение Пирсона (1) после переноса начала координат запишется в виде

,

тогда решение этого уравнения запишется в виде

    (6).

Линейным преобразованием независимого переменного

эта функция приводится с точностью до постоянного множителя  к весовой функции многочленов Чебышева – Эрмита, которая имеет вид

.

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

.

В этом случае условие ортогональности запишется в виде:

если

Полиномы Чебышева - Лагерра.

Пусть теперь многочлен (2) имеет один корень. Тогда уравнение (1) представимо в виде

.

Тогда его решение запишется в виде

.

Многочлены, ортогональные с таким весом, можно рассматривать как обобщение многочленов Чебышева – Лагерра, ортогональных с весом

.

Причем и здесь можно выразить эти многочлены через многочлены Чебышева – Лагерра , а условие ортогональности будет:

           если

Полиномы Якоби.

Предположим, что многочлен (2) имеет два различных действительных нуля. Тогда , и уравнение Пирсона (1) представимо в виде

,

где  и  - некоторые постоянные и . Тогда решение уравнения (1)

представимо в виде

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

.

И соответственно условие ортогональности будет иметь вид:

       если

Многочлены Чебышева I рода являются частным случаем многочленов Якоби, так как весовая функция, относительно которой ортогональны эти многочлены, имеет вид:

и получается при подстановке в весовую функцию многочленов Якоби параметров .

Многочлены Чебышева II рода так же являются частным случаем многочленов Якоби, так как весовая функция многочленов Чебышева II рода имеет вид

и получается при подстановке в весовую функцию многочленов Якоби параметров .

Следует так же отметить, что многочлены Лежандра являются частным случаем многочленов Якоби, так как весовая функция многочленов Лежандра

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

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

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

§ 1. Примеры нахождения кривых распределения вероятностей.

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

.

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

Ниже рассмотрено несколько примеров наиболее часто встретившихся форм кривой распределения I типа.

Пример 1.

Рассмотрим выборку: