Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MathCadTV 5semestr

.pdf
Скачиваний:
35
Добавлен:
10.05.2015
Размер:
1.93 Mб
Скачать

Искомое значение n= 6764.

Задание 4 Непрерывные случайные величины

Наиболее распространенные распределения непрерывных случайных величин

Равномерное распределение. Непрерывная случайная величина , принимающая значения на отрезке [a,b] распределена равномерно на [a,b], если плотность распределения p (x) и функция распределения случайной величины имеют соответственно вид:

0, x [a, b]

0, x a

 

a

 

 

 

 

x

 

p (x)

1

 

F (x)

 

 

, a x b

, x [a, b]

 

a

 

 

 

b

 

 

a

 

b

 

1, x b

 

 

 

 

 

 

 

Ниже приведены построенные в MathCAD графики плотности вероятностей и функции распределения случайной величины , принимающей значения на отрезке [0,1] и имеющей равномерное распределение.

1

 

1

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dun if(x 0 1)

0.5

 

 

 

 

 

 

punif(x 0 1)

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0.5

 

2

1

0.5

2

 

1

 

 

 

x

2

 

1

 

 

 

x

2

В MathCAD значения в точке х плотности распределения и функции распределения случайной величины, имеющей равномерное распределение на отрезке [a,b], вычисляются встроенными функциями соответственно dunif(x,a,b) и punif(x,a,b).

Экспоненциальное (показательное) распределение. Непрерывная случайная величина имеет показательное распределение с параметром 0 , если плотность распределения p (x) имеет вид

0, x 0,

p (x) e x , x 0.

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

0, x 0,

F (x) 1 e x , x 0.

Ниже приведены графики плотности вероятностей и функций распределения случайных величин, имеющих показательное распределение с параметрами =1 и =2, построенные в

MathCAD.

dexp(x 1)

1

 

 

pexp(x 1)

1

 

 

0.5

 

 

0.5

 

 

 

 

 

 

 

 

 

 

0

2

5

0

5

 

 

x

 

x

 

 

2

 

 

1

 

dexp(x 2)

1

 

pexp(x 2)

0.5

 

 

 

 

 

 

 

 

0

2

5

0

5

 

x

 

 

x

 

В MathCAD значения в точке x плотности распределения и функции распределения случайной величины, имеющей экспоненциальное распределение с параметром , вычисляются встроенными функциями соответственно dexp(x, ) и рехр(х, ).

Нормальное распределение. Это распределение играет исключительно важную роль в теории вероятностей и математической статистике. Случайная величина нормально распределена с параметрами a и , 0 , если ее плотность распределения имеет вид

 

 

1

 

 

 

(x a)

2

 

 

 

 

 

 

 

 

p (x)

 

 

 

 

2

 

2

exp

2

 

.

 

 

 

 

 

 

 

Если случайная величина имеет нормальное распределение с параметрами а и , то будем записывать это в виде ~ N(a, ) .

Случайная величина имеет стандартное нормальное распределение, если

a 0 и

1,

~ N(0,1) . Плотность стандартного нормального распределения имеет вид

 

 

 

 

 

 

1

 

 

 

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

p (x)

 

 

 

 

 

 

 

 

 

 

 

2

 

exp

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а его функция распределения — F (x) (x) , где (x)

 

 

- функция Лапласа:

 

 

1

 

 

x

 

 

 

 

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

dz .

 

 

2

 

 

 

 

2

 

 

 

Функция распределения нормальной величины ~ N(a, ) также выражается через функцию Лапласа:

x a

F (x)

 

.

 

 

 

Ниже приведены построенные в MathCAD графики плотности вероятностей и функций распределения для ~ N(0,1) , ~ N(1,2) .

dnorm(x 0 1)

1

 

 

 

 

 

 

 

 

 

 

 

 

pnormx( 0 1)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

1

 

 

 

 

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

3

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dnorm(x 1 2)

0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

pnormx( 1 2)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

0.33

2.33

5

3

2

1

0

1

2

3

4

5

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

В MathCAD значения в точке x плотности распределения и функции распределения нормальной случайной величины с параметрами а, а вычисляются встроенными функциями соответственно dnorm(x,a, ) и pnorm(x,a, ).

Распределение xи-квадрат ( 2-распределение). Пусть 1, 2, … n - независимые случайные величины, каждая из которых имеет стандартное нормальное распределение N(0, 1). Составим случайную величину

2 12 22 ... n2 .

Ее распределение называется 2-распределением с n степенями свободы. Для справочных целей приведем здесь выражение плотности распределения этой случайной величины:

Pn (z)

где Г(x) — гамма-функция Эйлера:

0, z 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2

 

 

x

 

 

1

 

 

 

 

 

 

 

n

 

z

 

e

 

, z 0 ,

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

Г

 

 

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Г(x) xz 1e z dz .

0

Ниже приведены графики плотности вероятностей и функций распределения для 2- распределения с двумя, четырьмя и восемью степенями свободы, построенные в MathCAD. Для сравнения приведены графики для ~ N(0,1) .

0.6

chisq(x 2)

chisq(x 4) 0.4 chisq(x 8)

norm(x 0 1)

0.2

 

 

 

 

 

 

 

 

 

 

4

2

0

2

4

6

8

10

 

 

 

 

x

 

 

 

 

1

chisq(x 2) chisq(x 4)

chisq(x 8) 0.5

normx( 0 1)

 

 

 

 

 

5

0

5

10

x

В MathCAD значения в точке x плотности распределения и функции 2-распределения с n степенями свободы вычисляется встроенными функциями соответственно dchisq(x,n) и pchisq(x, n).

ЗАДАНИЕ 5.6

Постройте графики плотности распределения и функции распределения 2 c указанным числом степеней свободы, равным k = №. Здесь № номер варианта

Порядок выполнения задания

1.Постройте графики плотности распределения 2 суказанным числом степеней свободы.

Пример выполнения задания

Постройте графики плотности распределения и функции распределения 2 с n = 2, 4, 8 степенями свободы. Примерный вариант выполнения задания приведен выше.

Распределение Стьюдента. Пусть случайная величина имеет стандартное нормальное распределение, а случайная величина n2 2-распределение с n степенями свободы. Если

и 2

независимы, то про случайную величину

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

2

n

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

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

 

 

 

 

 

n 1

 

 

 

 

 

 

n 1

 

 

 

1

 

 

Г

 

 

 

 

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

p (x)

 

 

 

 

1

 

 

 

 

 

, x R

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

 

n

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

При больших н распределение Стьюдента практически не отличается от N(0,1).

Ниже приведены графики плотности вероятностей и функций распределения для n =2, 5,10, построенные в MathCAD. Для сравнения приведены графики для ~ N(0,1) .

dt(x 2)

0.4

 

dt(x 5)

dt(x 10)

dnorm(x 0 1)

0.2

 

 

 

 

 

2

1

0

1

2

 

 

x

 

 

1

pt(x 2)

pt(x 5)

pt(x 10)

0.5

pnormx( 0 1)

 

 

 

 

 

 

 

3

1.5

0

1.5

3

 

 

 

x

 

 

 

В MathCAD значения в точке x плотности распределения и функции Стьюдента с n степенями свободы вычисляются встроенными функциями соответственно dt(x,n) и pt(x,n).

ЗАДАНИЕ 5.7

Постройте графики плотности распределения и функции распределения Стьюдента с указанным числом степеней свободы, равным k = №. Здесь № номер варианта.

Порядок выполнения задания

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

Пример выполнения задания

Постройте графики плотности распределения и функции распределения 2 с n = 2, 5,10 степенями свободы. Примерный вариант выполнения задания приведен выше.

F-распределение Фишера. Пусть случайные величины n2 и m2 независимы и имеют распределение 2 с n и m степенями свободы соответственно. Тогда случайная величина

Fn,m

2

n

n

 

имеет F-распределение с плотностью вероятности

2

m

 

m

 

 

 

n m

 

 

 

 

 

 

 

 

n 2

 

Г

 

 

 

 

 

 

n

n 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2

 

 

 

 

 

pF (x)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, x > 0.

n

 

m

 

 

 

 

 

 

 

 

 

n m

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nx 2

 

Г

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ниже приведены построенные в MathCAD графики плотности вероятностей и функций распределения для n = 2, 5 и m = 5, 2.

1

dF(x 2 5)

dF(x 5 2) 0.5

 

 

 

 

 

 

 

 

0

0.5

1

1.5

2

 

 

 

x

 

 

 

1

pF(x 2 5)

 

pF(x 5 2)

0.5

 

 

 

 

 

 

1

0

1

2

3

x

Указание. В MathCAD значения в точке x плотности распределения и функции Фишера вычисляются встроенными функциями соответственно dF(x, n,m) и pF(x, n,m).

ЗАДАНИЕ 5.8

Постройте графики плотности распределения и функции распределения Фишера для указанных значений n и m..

N

п

т

 

 

N

п т

 

N

п т

 

N

п

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

6

4

 

 

 

6

9

5

 

11

4

2

 

16

9

5

2

7

4

 

 

 

7

8

5

 

12

5

2

 

17

10

5

3

8

4

 

 

 

8

7

5

 

13

6

3

 

18

9

5

4

9

4

 

 

 

9

6

3

 

14

7

4

 

19

8

5

5

10

4

 

 

 

10

5

3

 

15

8

4

 

20

7

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Порядок выполнения задания

1.Постройте графики плотности распределения Фишера для указанных значений n и m.

2.Постройте графики функции распределения Фишера для указанных значений n и m.

Пример выполнения задания

Пример выполнения задания для распределения Фишера со значениями n = 2, 5 и m= 5, 2 приведен выше.

Распределение Парето. Распределение Парето часто применяется в экономических исследованиях. Плотность вероятностей для случайной величины, распределенной по Парето, имеет вид

0, x a

p (x) a x ( 1) , x a

где a 0, p 0 .

Как видно из этой формулы, случайная величина, распределенная по Парето, принимает значения только в области x > а.

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

соглашения об именах для имеющихся статистических функций, как функции переменной x и параметров a и p: dP(x, p, a) и рР(x, p, a) — соответственно для плотности вероятностей и функции распределения Парето (параметр обозначен как p).

dP(x p a)

 

p ap x ( p 1)

if x a

 

 

 

0

if

x a

 

pP(x p a)

 

if

x a

 

 

0

 

 

 

x

 

 

 

 

 

 

 

 

p ap t ( p 1) dt if x a

 

 

a

 

 

 

1

dP(x 1 1) dP(x 1 2) 0.5

 

 

 

 

 

 

 

0

1

2

3

4

5

 

 

 

x

 

 

1

pP(x 1 1) pP(x 1 2) 0.5

 

 

 

 

 

 

 

 

0

5

10

15

20

 

 

 

x

 

 

 

 

ЗАДАНИЕ 5.9

Постройте графики плотности распределения и функции распределения Парето для указанных значений a и .

N

а

 

N

а

 

N

а

 

N

а

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

1

6

1.9

1

11

2.1

2

16

2.9

2

2

1.5

1

7

2.0

1.5

12

2.5

2

17

3.0

2.5

3

1.6

1

8

2.1

1.5

13

2.6

2

18

3.1

2.5

4

1.7

1

9

2.2

1.5

14

2.7

2

19

3.2

2.5

5

1.8

1

10

2.3

1.5

15

2.8

2

20

3.3

2.5

 

 

 

 

 

 

 

 

 

 

 

 

Порядок выполнения задания

1.Определите плотность распределения Парето как функцию x и параметров a и .

2.2. Определите функцию распределения Парето как функцию х и параметров a и .

3.Постройте графики плотности распределения Парето для заданных значений a и .

4.Постройте графики функции распределения Парето для заданных значений a и ..

Пример выполнения задания

Постройте графики плотности распределения и функции распределения Парето для = 1 и a= 1, 2. Примерный вариант выполнения задания приведен выше.

Логистическое распределение. Это еще одно распределение, широко применяемое в экономических исследованиях. Оно имеет следующую функцию распределения:

F (x)

где и - параметры распределения.

1

 

 

, x R

 

 

 

 

 

x

1 exp

 

 

 

 

 

 

 

 

 

 

 

Плотность распределения вероятностей для логистического распределения вычисляется по формуле:

p (x) 1

 

 

 

x

 

 

exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, x R

 

 

 

 

x

2

1 exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По своим свойствам логистическое распределение очень похоже на нормальное.

Ниже приведены графики плотности вероятностей и функций распределения для 0 ,1,2 , построенные в MathCAD. Для сравнения приведены графики для нормального распределения N(0,2).

dlogis(x 0 1)

0.2

dlogis(x 0 2)

dnorm(x 0 2)

0.1

 

 

 

 

 

2

0

2

x

1

plogis(x 0 1)

plogis(x 0 2) 0.5

pnormx( 0 2)

 

 

 

 

 

 

10

5

0

5

10

 

 

x

 

 

 

В MathCAD значения в точке x плотности распределения и функции распределения для логистического распределения вычисляются встроенными функции соответственно dlogis(x,

, ) и plogis(x, , ).

ЗАДАНИЕ 5.10

Постройте графики плотности распределения и функции распределения для логистического распрения при значениях параметров = а и = для значений а и из задания 5.9.

Порядок выполнения задания

1.Постройте графики плотности распределения для логистического распределения с указанными значениями параметров.

2.Постройте графики функции распределения для логистического распределения с указанными значениями параметров.

Пример выполнения задания

Постройте графики плотности распределения и функции распределения для логистического распределения с параметрами = 0 и = 1,2. Примерный вариант выполнения задания приведен выше.

Квантили

При решении практических задач часто требуется найти значение я, при котором функция распределения случайной величины принимает заданное значение, т.е. требуется решить уравнение Fξ (x) = p. Решения такого уравнения в теории вероятностей называются квантилями.

Квантилъю xp (p-квантилью, квантилью уровня p) случайной величины ξ, имеющей функцию распределения Fξ(x), называют решение xp уравнения Fξ (x) = p, p (0,1).

Для некоторых p уравнение Fξ (x) = p может иметь несколько решений, для некоторых — ни одного. Это означает, что для соответствующей случайной величины некоторые квантили определены неоднозначно, а некоторые квантили не существуют.

Квантили, наиболее часто встречающиеся в практических задачах, имеют свои названия: медиана — квантиль уровня 0.5;

нижняя квартиль — квантиль уровня 0.25; верхняя квартиль — квантиль уровня 0.75;

децили — квантили уровней 0.1, 0.2, ..., 0.9; процентили — квантили уровней 0.01, 0.02, ..., 0.99.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]