![]() |
Вначале введем и кратко опишем гамма - функцию, которая встретилась уже и будет встречаться в дальнейшем.Введем математическое определение гамма - функции и применим интегрирование по частям:
,где
,
,
,
.Выполним подстановку:
=
.В результате мы получили рекуррентную формулу для вычисления значений гамма - функции.Если x - целое положительное число, x = n, то в соответствии с этой рекуррентной формулой имеем
.Вычислим G(1) отдельно:
Поэтому
В частности, G(1) = 0×Г(0) = 0! = 1. В дальнейшем нам понадобятся следующие значения гамма - функции от дробных аргументов:
,
,
.Последние два значения гамма-функции получены с использованием только что выведенной рекуррентной формулы.Пусть задана непрерывная дифференцируемая функция от случайной величины x: h = f (x). Известна плотность распределения случайной величины x: j (x). Задача состоит в том, чтобы найти плотность распределения y(y) случайной величины h. Подобная задача возникает в технике, когда случайные процессы или измеряемые величины, возмущенные случайными помехами, претерпевают нелинейные преобразования, и возникает задача прогнозирования характеристик сигнала, который получается в результате этого преобразования.Для вывода необходимого соотношения воспользуемся рис. 20. На этом рисунке представлены: функция преобразования y = f(x) и плотность распределения j(x). Функция преобразования предполагается монотонной, и это свойство функции преобразования практически всегда имеет место в технических устройствах: средствах измерения, измерительныхпреобразователях и регуляторах. В силу взаимной однозначности преобразования случайная величина h принимает значения из интервала Dy, в точности с той же вероятностью, с которой случайная величина x принимает значения из интервала Dx. Поскольку вероятностная мера интервала есть площадь под кривой плотности распределения на этом интервале, это означает, что площади заштрихованных фигур на рис. 20 должны быть равны:
,где
,
- точки, находящиеся внутри выделенных интервалов Dx и
,
- значение искомой плотности распределения в точке
.Из этого выражения следует:
.Заметим здесь, что ширина интервала
, в который преобразуется интервал Dx, не зависит от знака производной функции преобразования, и это обстоятельство мы учтем при выполнении предельного перехода
.В силу инвариантности первого дифференциала производная
выражается через производную от обратной функции:
Кроме того в выражении для y(y) необходимо выразить аргумент x плотности распределения j(x) через y с помощью обратной функции:
. В итоге окончательно получим:
. Пример: Случайная величина h образуется в результате линейного преобразования случайной величины
:
. В данном случае реализуется функциональное преобразование
, обратная функция
, производная от нее по y равна
. В результате подстановки в общую формулу получим:
.Это означает, что любое линейное преобразование не изменяет вид плотности распределения случайной величины. Изменяется лишь масштаб и смещение от начала координат.
![]() |
17.Вывод формулы для плотности распределения случайной величины и для ее характеристической функции,гдеx–случайная величина с нормальной плотностью распределения N(0;1) Случайная величина x распределена нормально:
. Функция преобразования
. Эта ситуация представлена на рис. 21. Из рисунка видно, что в силу двузначности обратной функции случайная величина h принимает значения из интервала
, когда случайная величина x принимает значения в одном из двух выделенных интервалов Dx. Поэтому для данного примера исходное выражение должно быть изменено следующим образом:
. Из этого следует соответствующее изменение общей формулы:
. Для данного примера
,
,
,
В конечном итоге, после подстановки в общую формулу получим искомую плотность распределения: .Найдем характеристическую функцию этого распределения.
.Сделаем замену переменной интегрирования:
.В результате этой замены получим выражение с участием гамма-функции:
Поскольку
, окончательно получим:
18.Вывод формулы для плотности распределения ,где
–произвольная функция распределения непрерывной случайной величины Пусть
- интегральная функция распределения случайной величины x. Образуем случайную величину h, как функцию от случайной величины x:
. Задача состоит в том, чтобы найти плотность распределения y (y) случайной величины h. Воспользуемся полученным ранее выражением
.В нашем случае в качестве функции y = f(x) выступает функция
, производная от которой по x есть плотность распределения j(x). Поэтому
.Таким образом, оказывается, что случайная величина h, полученная в результате функционального преобразования любой непрерывной случайной величины x путем ее подстановки в ее же интегральную функцию распределения, распределена равномерно в интервале [0,1] вне зависимости от вида функции распределения величины x. Полученный результат имеет два полезных применения. Первое. Машинное моделирование случайных чисел с заданной интегральной функцией распределения. Технология моделирования такова:- задается функция распределения F(x),- по стандартным программам генерируются случайные числа
, распределенные равномерно в интервале [0, 1], - случайные числа
, распределенные в соответствии с заданной функцией распределения F(x) получаются, как решения уравнений
. Второе. Статистическое оценивание параметров и характеристик случайных величин по результатам экспериментов, вне зависимости от вида распределения исследуемой случайной величины. Это применение будет изложено ниже в разделе 2. Математическая статистика.