Matlab  |  Mathcad  |  Maple  |  Mathematica  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Банк задач  |  Консультации & Форум  |  Download  |  Ссылки  |  Конкурсы
Научно-практический журнал "Exponenta Pro. Математика в приложениях". Вышел 1/2004 номер журнала
 
Лабораторные работы по курсу "Теория вероятностей"

Вернуться на страницу <Методические разработки>
В начало | Л.р.1 | Л.р.2 | Л.р.3 | Л.р.4 | Л.р.5 | Л.р.6 | Л.р.7

Л.р.9 | Л.р.10 | Л.р.11 | Л.р.12 | Л.р.13 | Л.р.14 | Л.р.15

Лабораторная работа 8
Функции от случайных величин

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

matrix([[xi, x[1], x[2], %?], [p, p[1], p[2], %?]])...

и пусть y=g(x) - монотонная функция от действительного аргумента x. Тогда дискретная случайная величина h =g(x), которая является функцией от x , имеет ряд распределения

matrix([[eta, g(x[1]), g(x[2]), %?], [p, p[1], p[2]...

Если y=g(x) - немонотонная функция, то среди ее значений могут быть равные. В этом случае столбцы с равными значениями функции объединяют в один столбец, а соответствующие вероятности складывают.

Пусть теперь  x - непрерывная случайная величина с функцией распределения F[xi] (x) и плотностью вероятности f[xi] (x) и пусть функция y=g(x) монотонно возрастает, x= g^(-1) (y) - обратная функция. Тогда

F[eta] (y)=P{ eta <y}=P{g(  x )<y}=P{   x < g^(-1) (y)}= F[xi] ( g^(-1) (y)).

Дифференцируя последнее равенство по y, имеем (если g(x) дифференцируема)

diff(F[eta],y) (y)= diff(F[xi],x) ( g^(-1) (y)) diff(g^(-1),y) (y),

откуда получаем соотношение между плотностями

f[eta] (y)= f[xi] ( g^(-1) (y)) diff(g^(-1),y) (y).

Если y=g(x) - монотонно убывающая функция, то аналогично получаются следующие соотношения:

F[eta] (y)=1- F[xi] ( g^(-1) (y)), f[eta] (y)=- f[xi] ( g^(-1) (y)) diff(g^(-1),y) (y).

Если случайные величины xi[1] и xi[2] абсолютно непрерывны и независимы, то случайная величина xi[1]+xi[2] тоже абсолютно непрерывна и ее плотность распределения определяется по формуле

f[xi[1]+x[2]](x) = int(f[xi[1]](u)*f[xi[2]](x-u),u ... .

Задача 1 . Дискретная случайная величина  x имеет ряд распределения matrix([[xi, Pi/4, Pi/2, 3*Pi/4], [p, .2, .7, .1]])...

Построить ряд распределения случайной величины eta = sin(xi) .

Решение. Построим ряд распределения случайной величины h =sin(   x )

> matrix([[eta,sin(Pi/4), sin(Pi/2), sin(3*Pi/4)], [p,0.2, 0.7, 0.1]]);

matrix([[eta, 1/2*sqrt(2), 1, 1/2*sqrt(2)], [p, .2,...

или

> matrix([[eta,sqrt(2)/2,1],[p,0.3,0.7]]);

matrix([[eta, 1/2*sqrt(2), 1], [p, .3, .7]])

Задача 2 . Дискретная случайная величина  x имеет ряд распределения matrix([[xi, -2, -1, 0, 1, 2], [p, .1, .2, .3, .3, ... . Построить ряд распределения случайных величин hx2 +1; z =|  x |.

Задача 3 . Непрерывная случайная величина   x имеет плотность вероятности f[xi] (x). Найти плотность вероятности случайной величины h =3  x .

Воспользуемся формулой f[eta] (y)= f[xi] ( g^(-1) (y)) diff(g^(-1),y) (y) в силу того, что функция g(x)=3x монотонно возрастающая. Обратная функция g^(-1) (y)=y/3. Найдем ее производную

> diff(y/3,y);

1/3

Очевидно, что

> f[eta](y):=1/3*f[xi](y/3);

f[eta](y) := 1/3*f[xi](1/3*y)

Задача 4 . Случайная величина x  распределена по нормальному закону с параметрами ( 0, sigma^2 ). Найти закон распределения обратной ей величины eta = 1/xi .

Задача 5 . Случайная величина x  имеет показательное распределение с плотностью вероятности f[xi] = exp(-x) , x>0. Найти закон распределения случайной величины eta = exp(-xi) .

Задача 6 . Случайная величина x  распределена равномерно на отрезке [-1, 2]. Найти закон распределения случайной величины eta = xi^2 .

Задача 7 . Докажите, что если случайная величина x  имеет нормальное распределение с параметрами ( a, sigma^2 ), то случайная величина eta = A*xi+B , A <> 0 , имеет нормальное распределение с параметрами ( A*a, A^2*sigma^2 ).

Задача 8 . Величины x 1 и x 2 независимы, одинаково распределены и имеют показательное распределение: f[xi[1]](x) = alpha*exp(-alpha*x) , f[xi[2]](x) = alpha*exp(-alpha*x) , x>0. Найти закон распределения их суммы.

В начало | Л.р.1 | Л.р.2 | Л.р.3 | Л.р.4 | Л.р.5 | Л.р.6 | Л.р.7

| Л.р.9 | Л.р.10 | Л.р.11 | Л.р.12 | Л.р.13 | Л.р.14 | Л.р.15
Вернуться на страницу <Методические разработки>

Карта сайта | На первую страницу | Поиск |О проекте |Сотрудничество |
Exponenta Pro | Matlab.ru

Наши баннеры


Copyright © 2000-2003. Компания SoftLine. Все права защищены.

Дата последнего обновления информации на сайте: 11.05.04
Сайт начал работу 1.09.00

Программное обеспечение Microsoft, Macromedia, VERITAS, Novell, Borland, Symantec, Oracle и др.