numplot 1.0.0__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- numplot-1.0.0/PKG-INFO +38 -0
- numplot-1.0.0/README.md +17 -0
- numplot-1.0.0/numplot/__init__.py +17 -0
- numplot-1.0.0/numplot/task_1.py +495 -0
- numplot-1.0.0/numplot/task_2.py +522 -0
- numplot-1.0.0/numplot/task_3.py +536 -0
- numplot-1.0.0/numplot/task_4.py +888 -0
- numplot-1.0.0/numplot/task_5.py +196 -0
- numplot-1.0.0/numplot.egg-info/PKG-INFO +38 -0
- numplot-1.0.0/numplot.egg-info/SOURCES.txt +13 -0
- numplot-1.0.0/numplot.egg-info/dependency_links.txt +1 -0
- numplot-1.0.0/numplot.egg-info/top_level.txt +1 -0
- numplot-1.0.0/setup.cfg +4 -0
- numplot-1.0.0/setup.py +25 -0
numplot-1.0.0/PKG-INFO
ADDED
|
@@ -0,0 +1,38 @@
|
|
|
1
|
+
Metadata-Version: 2.2
|
|
2
|
+
Name: numplot
|
|
3
|
+
Version: 1.0.0
|
|
4
|
+
Summary: Package mathstats contains mathematical and statistical functions, distributions, tests
|
|
5
|
+
Author: amogusbazed
|
|
6
|
+
Author-email: amogusbazed@gmail.com
|
|
7
|
+
Keywords: files speedfiles
|
|
8
|
+
Classifier: Programming Language :: Python :: 3.13
|
|
9
|
+
Classifier: License :: OSI Approved :: MIT License
|
|
10
|
+
Classifier: Operating System :: OS Independent
|
|
11
|
+
Requires-Python: >=3.6
|
|
12
|
+
Description-Content-Type: text/markdown
|
|
13
|
+
Dynamic: author
|
|
14
|
+
Dynamic: author-email
|
|
15
|
+
Dynamic: classifier
|
|
16
|
+
Dynamic: description
|
|
17
|
+
Dynamic: description-content-type
|
|
18
|
+
Dynamic: keywords
|
|
19
|
+
Dynamic: requires-python
|
|
20
|
+
Dynamic: summary
|
|
21
|
+
|
|
22
|
+
# Numplot File Library #
|
|
23
|
+
|
|
24
|
+
## What is this? ##
|
|
25
|
+
The module allows you to work with files in just one line of code, without the need to manually open and close the file each time
|
|
26
|
+
|
|
27
|
+
## Quick Guide ##
|
|
28
|
+
The module is based on the following structure:
|
|
29
|
+
|
|
30
|
+
|
|
31
|
+
f = open('data.txt')
|
|
32
|
+
data = f.readlines()
|
|
33
|
+
f.close()
|
|
34
|
+
|
|
35
|
+
Which Python provides by standard.
|
|
36
|
+
|
|
37
|
+
|
|
38
|
+
----------
|
numplot-1.0.0/README.md
ADDED
|
@@ -0,0 +1,17 @@
|
|
|
1
|
+
# Numplot File Library #
|
|
2
|
+
|
|
3
|
+
## What is this? ##
|
|
4
|
+
The module allows you to work with files in just one line of code, without the need to manually open and close the file each time
|
|
5
|
+
|
|
6
|
+
## Quick Guide ##
|
|
7
|
+
The module is based on the following structure:
|
|
8
|
+
|
|
9
|
+
|
|
10
|
+
f = open('data.txt')
|
|
11
|
+
data = f.readlines()
|
|
12
|
+
f.close()
|
|
13
|
+
|
|
14
|
+
Which Python provides by standard.
|
|
15
|
+
|
|
16
|
+
|
|
17
|
+
----------
|
|
@@ -0,0 +1,17 @@
|
|
|
1
|
+
from .task_1 import t_1
|
|
2
|
+
from .task_2 import t_2
|
|
3
|
+
from .task_3 import t_3
|
|
4
|
+
from .task_4 import t_4
|
|
5
|
+
from .task_5 import t_5
|
|
6
|
+
|
|
7
|
+
__all__ = ['t_1', 't_2', 't_3', 't_4', 't_5']
|
|
8
|
+
|
|
9
|
+
def h():
|
|
10
|
+
text = """
|
|
11
|
+
t_1 ()
|
|
12
|
+
t_2 ()
|
|
13
|
+
t_3 ()
|
|
14
|
+
t_4 ()
|
|
15
|
+
t_5 ()
|
|
16
|
+
"""
|
|
17
|
+
return text
|
|
@@ -0,0 +1,495 @@
|
|
|
1
|
+
def t_1(var):
|
|
2
|
+
if var == 1:
|
|
3
|
+
text = """
|
|
4
|
+
Доходности акций двух компаний являются случайными величинами X и Y, имеющими совместное нормальное распределение с центром рассеивания (2; 2) и ковариационной матрицей:
|
|
5
|
+
C = | 4 -2 |
|
|
6
|
+
| -2 9 |
|
|
7
|
+
Найти P(X > 0 | Y = 0). В каком соотношении нужно приобрести акции этих компаний, чтобы риск (дисперсия) получившегося портфеля был минимальным? Можно ли утверждать, что случайные величины X + Y и 7X - 2Y независимы? Подсказка: если R - доходность портфеля, то R = αX + (1 - α)Y.
|
|
8
|
+
Решение:
|
|
9
|
+
P(X > 0 | Y = 0): Условное распределение X при Y = 0 также будет нормальным. Найдем параметры этого условного распределения.
|
|
10
|
+
E(X|Y=0) = E(X) + Cov(X,Y)/Var(Y) * (0 - E(Y)) = 2 + (-2)/9 * (0 - 2) = 2 + 4/9 = 22/9
|
|
11
|
+
Var(X|Y=0) = Var(X) - Cov(X,Y)²/Var(Y) = 4 - (-2)²/9 = 4 - 4/9 = 32/9
|
|
12
|
+
Таким образом, X|Y=0 ~ N(22/9; 32/9). Теперь можно найти P(X > 0 | Y = 0), используя функцию стандартного нормального распределения: P(X > 0 | Y=0) = P(Z > (0 - 22/9) / sqrt(32/9)) = P(Z > -22/(3*sqrt(32))) ≈ P(Z > -1.29) ≈ 0.9015 (Значение можно уточнить по таблице или с помощью калькулятора)
|
|
13
|
+
Минимизация риска портфеля: Дисперсия портфеля R = αX + (1 - α)Y равна:
|
|
14
|
+
Var(R) = α²Var(X) + (1 - α)²Var(Y) + 2α(1 - α)Cov(X,Y) = 4α² + 9(1 - α)² - 4α(1 - α) = 17α² - 22α + 9
|
|
15
|
+
Для минимизации дисперсии возьмем производную по α и приравняем к нулю:
|
|
16
|
+
d(Var(R))/dα = 34α - 22 = 0 => α = 11/17
|
|
17
|
+
Таким образом, нужно вложить 11/17 средств в акции X и 6/17 в акции Y.
|
|
18
|
+
Независимость X + Y и 7X - 2Y: Две случайные величины независимы, если их ковариация равна нулю. Найдем ковариацию:
|
|
19
|
+
Cov(X + Y, 7X - 2Y) = 7Var(X) - 2Cov(X,Y) + 7Cov(Y,X) - 2Var(Y) = 74 - 2(-2) + 7*(-2) - 2*9 = 28 + 4 - 14 - 18 = 0
|
|
20
|
+
Так как ковариация равна нулю, X + Y и 7X - 2Y независимы.
|
|
21
|
+
"""
|
|
22
|
+
if var == 2:
|
|
23
|
+
text = """Для случайно выбранного домохозяйства случайные величины X и Y принимают значения, равные доле расходов на продукты питания и алкоголь плюс табак соответственно. Случайный вектор (X, Y) хорошо описывается двумерным нормальным законом распределения с математическим ожиданием (0.45; 0.16) и ковариационной матрицей:
|
|
24
|
+
C = | 0.144 -0.09 |
|
|
25
|
+
| -0.09 1 |
|
|
26
|
+
Найдите: а) Вероятность того, что домохозяйство тратит более половины своих доходов на питание. б) Вероятность того, что домохозяйство тратит более половины своих доходов на алкогольную и табачную продукцию и продукты питания. в) Ожидаемую долю расходов на алкоголь и табак для домохозяйства, которое тратит на питание четверть своих доходов.
|
|
27
|
+
Решение:
|
|
28
|
+
а) P(X > 0.5): X ~ N(0.45; 0.144). Стандартизуем: Z = (X - 0.45) / sqrt(0.144) = (X - 0.45) / 0.12. P(X > 0.5) = P(Z > (0.5 - 0.45) / 0.12) = P(Z > 0.417) ≈ 1 - Φ(0.417) ≈ 1 - 0.6615 ≈ 0.3385 (где Φ - функция стандартного нормального распределения).
|
|
29
|
+
б) P(X + Y > 0.5): Найдем распределение суммы X + Y. E(X+Y) = E(X) + E(Y) = 0.45 + 0.16 = 0.61. Var(X+Y) = Var(X) + Var(Y) + 2Cov(X,Y) = 0.144 + 1 + 2*(-0.09) = 1.064. X+Y ~ N(0.61; 1.064). Стандартизуем: Z = (X + Y - 0.61) / sqrt(1.064). P(X + Y > 0.5) = P(Z > (0.5 - 0.61) / sqrt(1.064)) = P(Z > -0.106) ≈ Φ(0.106) ≈ 0.542.
|
|
30
|
+
в) E(Y|X=0.25): E(Y|X=x) = E(Y) + Cov(X,Y)/Var(X) * (x - E(X)) = 0.16 + (-0.09)/0.144 * (0.25 - 0.45) = 0.16 + (-0.625) * (-0.2) = 0.16 + 0.125 = 0.285.
|
|
31
|
+
"""
|
|
32
|
+
if var == 3:
|
|
33
|
+
text = """Плотность распределения случайного вектора (X, Y) имеет вид:
|
|
34
|
+
fₓy(x, y) = (1/(2π)) * e^(-(5x² - xy + y²)/2)
|
|
35
|
+
Найдите условное математическое ожидание E(X|Y = y).
|
|
36
|
+
Решение:
|
|
37
|
+
1. Преобразуем выражение в показателе экспоненты, выделяя полный квадрат относительно x:
|
|
38
|
+
5x² - xy + y² = 5(x² - (1/5)xy + y²/5) = 5(x² - (1/5)xy + 1/100 y² - 1/100 y² + y²/5) = 5(x - y/10)² + (19/20)y²
|
|
39
|
+
Таким образом, fₓy(x, y) = (1/(2π)) * e^(-(5(x - y/10)² + (19/20)y²)/2) = (1/(2π)) * e^(-5(x - y/10)²/2) * e^(-(19/20)y²/2)
|
|
40
|
+
2. Найдем маргинальную плотность Y:
|
|
41
|
+
fᵧ(y) = ∫₋∞^∞ fₓy(x, y) dx = ∫₋∞^∞ (1/(2π)) * e^(-5(x - y/10)²/2) * e^(-(19/20)y²/2) dx = e^(-(19/40)y²) * ∫₋∞^∞ (1/(2π)) * e^(-5(x - y/10)²/2) dx
|
|
42
|
+
Заметим, что ∫₋∞^∞ (1/(2π)) * e^(-5(x - y/10)²/2) dx = (1/√(2π/5))∫₋∞^∞√(5)/(√(2π)) e^(-5(x-y/10)²/2)dx = 1/√(2π/5)
|
|
43
|
+
Поэтому, fᵧ(y) = (1/√(2π/5)) * e^(-(19/40)y²)
|
|
44
|
+
3. Найдем условную плотность X при Y = y:
|
|
45
|
+
fₓ|ᵧ(x|y) = fₓy(x, y) / fᵧ(y) = [(1/(2π)) * e^(-5(x - y/10)²/2) * e^(-(19/20)y²/2)] / [(1/√(2π/5)) * e^(-(19/40)y²)] = (√(5)/(√(2π))) * e^(-5(x - y/10)²/2)
|
|
46
|
+
4. Условная плотность fₓ|ᵧ(x|y) является плотностью нормального распределения с параметрами μ = y/10 и σ² = 1/5.
|
|
47
|
+
5. Следовательно, условное математическое ожидание E(X|Y = y) равно μ, то есть y/10.
|
|
48
|
+
E(X|Y = y) = y/10
|
|
49
|
+
Таким образом, с исправленной плотностью, условное математическое ожидание E(X|Y = y) равно y/10.
|
|
50
|
+
"""
|
|
51
|
+
if var == 4:
|
|
52
|
+
text = """Ежедневное количество нарушений правил дорожного движения в некоторой местности описывается случайной величиной Х, распределенной по биномиальному закону с параметрами n = 500 и p = 0.7. А сумма штрафа (в рублях) каждого нарушения описывается случайной величиной Y, которая принимает значение 600 рублей с вероятностью 0.48, 950 рублей с вероятностью 0.14 и 1700 рублей в остальных случаях. Вычислите ежедневную среднюю сумму выписанных штрафов и ее среднеквадратическое отклонение.
|
|
53
|
+
Решение:
|
|
54
|
+
1. Найдем математическое ожидание и дисперсию X: E(X) = n * p = 500 * 0.7 = 350 Var(X) = n * p * (1 - p) = 500 * 0.7 * 0.3 = 105
|
|
55
|
+
2. Найдем математическое ожидание и дисперсию Y: P(Y = 1700) = 1 - 0.48 - 0.14 = 0.38 E(Y) = 600 * 0.48 + 950 * 0.14 + 1700 * 0.38 = 288 + 133 + 646 = 1067 E(Y²) = 600² * 0.48 + 950² * 0.14 + 1700² * 0.38 = 172800 + 126350 + 1094200 = 1393350 Var(Y) = E(Y²) - (E(Y))² = 1393350 - 1067² = 1393350 - 1138489 = 254861
|
|
56
|
+
3. Ежедневная средняя сумма выписанных штрафов равна E(XY). Так как X и Y независимы, E(XY) = E(X)E(Y) = 350 * 1067 = 373450
|
|
57
|
+
4. Дисперсия ежедневной суммы штрафов Var(XY) = E(X²Y²) - (E(XY))². E(X²) = Var(X) + (E(X))² = 105 + 350² = 122605 E(X²Y²) = E(X²)E(Y²) = 122605*1393350 = 170830470750 Var(XY)=170830470750 - 373450^2= 170830470750 - 139464025000 = 31366445750 Среднеквадратическое отклонение: sqrt(Var(XY)) = sqrt(31366445750) ≈ 177105.74
|
|
58
|
+
"""
|
|
59
|
+
if var == 5:
|
|
60
|
+
text = """Случайные величины Х и Y независимы и с равной вероятностью принимают значения: 1, 2, 3, 4, 5. Найдите Е(Х|Х > Y) и Var(X|X > Y).
|
|
61
|
+
Решение:
|
|
62
|
+
1. Найдем вероятность события X > Y. Всего возможно 5 * 5 = 25 равновероятных исходов. Исходы, где X > Y: (2,1), (3,1), (3,2), (4,1), (4,2), (4,3), (5,1), (5,2), (5,3), (5,4). Итого 10 исходов. P(X > Y) = 10/25 = 2/5.
|
|
63
|
+
2. Найдем условное распределение X при условии X > Y. Возможные значения X при X > Y: 2, 3, 4, 5.
|
|
64
|
+
o P(X=2|X>Y) = P(X=2, X>Y) / P(X>Y) = (1/25) / (2/5) = 1/10
|
|
65
|
+
o P(X=3|X>Y) = P(X=3, X>Y) / P(X>Y) = (2/25) / (2/5) = 1/5
|
|
66
|
+
o P(X=4|X>Y) = P(X=4, X>Y) / P(X>Y) = (3/25) / (2/5) = 3/10
|
|
67
|
+
o P(X=5|X>Y) = P(X=5, X>Y) / P(X>Y) = (4/25) / (2/5) = 2/5
|
|
68
|
+
3. Найдем E(X|X > Y):
|
|
69
|
+
E(X|X > Y) = 2*(1/10) + 3*(1/5) + 4*(3/10) + 5*(2/5) = 1/5 + 3/5 + 6/5 + 10/5 = 20/5 = 4
|
|
70
|
+
4. Найдем E(X²|X > Y):
|
|
71
|
+
E(X²|X > Y) = 2²*(1/10) + 3²*(1/5) + 4²*(3/10) + 5²*(2/5) = 4/10 + 9/5 + 48/10 + 50/5 = 2/5 + 9/5 + 24/5 + 50/5 = 85/5 = 17
|
|
72
|
+
5. Найдем Var(X|X > Y):
|
|
73
|
+
Var(X|X > Y) = E(X²|X > Y) - (E(X|X > Y))² = 17 - 4² = 17 - 16 = 1
|
|
74
|
+
Ответ: E(X|X > Y) = 4, Var(X|X > Y) = 1.
|
|
75
|
+
"""
|
|
76
|
+
if var == 6:
|
|
77
|
+
text = """1. Ежедневное количество покупателей магазина, совершивших покупку, описывается случайной величиной Х, распределенной по биномиальному закону с параметрами n=800 и p=0.54. А сумма чека (в рублях) каждого из покупателей описывается случайной величиной Ү, распределенной по нормальному закону с параметрами m=4550 и σ=60. Вычислите ежедневную среднюю выручку магазина и ее дисперсию.
|
|
78
|
+
Решение:
|
|
79
|
+
1. Найдем математическое ожидание и дисперсию X (количество покупателей): E(X) = n * p = 800 * 0.54 = 432 Var(X) = n * p * (1 - p) = 800 * 0.54 * 0.46 = 198.72
|
|
80
|
+
2. Найдем математическое ожидание и дисперсию Y (сумма чека): E(Y) = m = 4550 Var(Y) = σ² = 60² = 3600
|
|
81
|
+
3. Ежедневная средняя выручка магазина равна E(XY). Так как X и Y независимы, E(XY) = E(X)E(Y) = 432 * 4550 = 1965600
|
|
82
|
+
4. Дисперсия ежедневной выручки Var(XY). Поскольку X и Y независимы, для дисперсии произведения случайных величин справедлива формула:
|
|
83
|
+
Var(XY) = E(X²)E(Y²) - (E(X)E(Y))² = (Var(X) + E(X)²) * (Var(Y) + E(Y)²) - (E(X)E(Y))²
|
|
84
|
+
Подставим известные значения:
|
|
85
|
+
Var(XY) = (198.72 + 432²) * (3600 + 4550²) - (432 * 4550)² Var(XY) = (198.72 + 186624) * (3600 + 20702500) - 1965600² Var(XY) = 186822.72 * 20706100 - 386343360000 Var(XY) = 386007771712 - 386343360000 Var(XY) = 386343360000 - 386007771712 = 335588288
|
|
86
|
+
"""
|
|
87
|
+
if var == 7:
|
|
88
|
+
text = """Ежемесячные расходы владельцев А и В описываются случайным вектором (X, Y), (X — расходы владельца А, Y — расходы владельца В), имеющим равномерное распределение в треугольнике, задаваемом ограничениями D = {(x, y): x ≥ 0, y ≥ 0, x + y ≤ 1}. Найдите:
|
|
89
|
+
а) Вероятность того, что совокупные расходы превысят половину бюджета, P(X + Y > 1/2); б) Плотность распределения расходов владельца В; в) Вероятность того, что расходы владельца А составили менее трети бюджета, если известно, что владелец А израсходовал более половины семейного бюджета; г) Условную плотность распределения и условное математическое ожидание расходов владельца А, при условии, что владелец В израсходовал половину бюджета; д) Математическое ожидание условного математического ожидания расходов владельца А, то есть E(E(X|Y)); е) Коэффициент корреляции расходов владельца А и владельца В.
|
|
90
|
+
Решение:
|
|
91
|
+
а) Площадь треугольника D равна 1/2. Площадь области, где X + Y > 1/2, также является треугольником с вершинами (1/2, 1/2), (1, 0), (0, 1) и площадью 1/8. Следовательно, P(X + Y > 1/2) = (1/8) / (1/2) = 1/4 = 0.25
|
|
92
|
+
б) Маргинальная плотность Y: fᵧ(y) = ∫₀^(1-y) fₓy(x, y) dx = ∫₀^(1-y) 2 dx = 2(1 - y), 0 ≤ y ≤ 1. (Плотность fxy(x,y) равна 2, так как площадь треугольника 1/2).
|
|
93
|
+
в) P(X < 1/3 | X > 1/2) = P(1/2 < X < 1/3) / P(X > 1/2). Событие {1/2 < X < 1/3} невозможно, так как 1/2 > 1/3. Значит, вероятность равна 0.
|
|
94
|
+
г) Условная плотность X при Y=1/2: fₓ|ᵧ(x|1/2) = fₓy(x, 1/2) / fᵧ(1/2) = 2 / (2 * (1 - 1/2)) = 2 / 1 = 2, при 0 ≤ x ≤ 1/2. E(X|Y = 1/2) = ∫₀^(1/2) x * 2 dx = [x²]₀^(1/2) = 1/4.
|
|
95
|
+
д) E(E(X|Y)). По свойству полного математического ожидания E(E(X|Y)) = E(X). E(X) = ∫₀¹ ∫₀^(1-y) x * 2 dx dy = ∫₀¹ [x²]₀^(1-y) dy = ∫₀¹ (1-y)² dy = ∫₀¹ (1 - 2y + y²) dy = [y - y² + y³/3]₀¹ = 1 - 1 + 1/3 = 1/3.
|
|
96
|
+
е) Коэффициент корреляции: Cov(X, Y) = E(XY) - E(X)E(Y). E(XY) = ∫₀¹ ∫₀^(1-y) xy * 2 dx dy = ∫₀¹ [x²y]₀^(1-y) dy = ∫₀¹ (1-y)²y dy = ∫₀¹ (y - 2y² + y³) dy = [y²/2 - 2y³/3 + y⁴/4]₀¹ = 1/2 - 2/3 + 1/4 = 1/12. E(X) = 1/3, E(Y) = 1/3 (в силу симметрии). Cov(X, Y) = 1/12 - (1/3)(1/3) = 1/12 - 1/9 = (3-4)/36 = -1/36.
|
|
97
|
+
Var(X) = E(X²) - (E(X))². E(X²) = ∫₀¹ ∫₀^(1-y) x² * 2 dx dy = ∫₀¹ [2x³/3]₀^(1-y) dy = (2/3)∫₀¹ (1-y)³ dy = (2/3)[-(1-y)⁴/4]₀¹ = 1/6 Var(X) = 1/6 - (1/3)² = 1/6 - 1/9 = 1/18. Var(Y) = 1/18 (в силу симметрии).
|
|
98
|
+
ρ = Cov(X, Y) / (sqrt(Var(X)) * sqrt(Var(Y))) = (-1/36) / (sqrt(1/18) * sqrt(1/18)) = (-1/36) / (1/18) = -1/2.
|
|
99
|
+
"""
|
|
100
|
+
if var == 8:
|
|
101
|
+
text = """1. Плотность распределения случайного вектора (X, Y) имеет вид:
|
|
102
|
+
fₓy(x, y) = (1/π) * e^(-(x²/2) - xy - (5y²/2))
|
|
103
|
+
Найдите условное математическое ожидание E(X|Y = y).
|
|
104
|
+
Решение:
|
|
105
|
+
1. Преобразуем выражение в показателе экспоненты, выделяя полный квадрат относительно x:
|
|
106
|
+
-(x²/2) - xy - (5y²/2) = -(1/2)(x² + 2xy + 5y²) = -(1/2)(x² + 2xy + y² - y² + 5y²) = -(1/2)((x + y)² + 4y²) = -(1/2)(x + y)² - 2y²
|
|
107
|
+
Таким образом, fₓy(x, y) = (1/π) * e^(-(x + y)²/2) * e^(-2y²)
|
|
108
|
+
2. Найдем маргинальную плотность Y:
|
|
109
|
+
fᵧ(y) = ∫₋∞^∞ fₓy(x, y) dx = ∫₋∞^∞ (1/π) * e^(-(x + y)²/2) * e^(-2y²) dx = e^(-2y²) * ∫₋∞^∞ (1/π) * e^(-(x + y)²/2) dx
|
|
110
|
+
Заметим, что ∫₋∞^∞ (1/π) * e^(-(x + y)²/2) dx = (1/√(2π))∫₋∞^∞√(2)/π e^(-(x+y)²/2)dx = 1/√(π/2)
|
|
111
|
+
fᵧ(y) = (1/√(π/2))e^(-2y²) = √(2/π) * e^(-2y²)
|
|
112
|
+
3. Найдем условную плотность X при Y = y:
|
|
113
|
+
fₓ|ᵧ(x|y) = fₓy(x, y) / fᵧ(y) = [(1/π) * e^(-(x + y)²/2) * e^(-2y²)] / [√(2/π) * e^(-2y²)] = (1/√(π/2)) * e^(-(x + y)²/2)
|
|
114
|
+
4. Условная плотность fₓ|ᵧ(x|y) является плотностью нормального распределения с параметрами μ = -y и σ² = 1.
|
|
115
|
+
5. Следовательно, условное математическое ожидание E(X|Y = y) равно μ, то есть -y.
|
|
116
|
+
E(X|Y = y) = -y
|
|
117
|
+
"""
|
|
118
|
+
if var == 9:
|
|
119
|
+
text = """Количество книг, проданных в день, описывается дискретной случайной величиной Y, распределенной по биномиальному закону с параметрами n=260 и p=0.88. При этом выручка книжного магазина за каждую j-ую покупку книги зависит от жанра (фантастика, детективы, новеллы) и описывается случайной величиной Yⱼ, принимающей значения 500, 300, 100 рублей с вероятностями соответственно: 0.67, 0.3, 0.13. Вычислите значения для среднего (Е) и дисперсии суммарного ежедневного дохода.
|
|
120
|
+
Решение:
|
|
121
|
+
1. Найдем математическое ожидание и дисперсию Y (количество проданных книг): E(Y) = n * p = 260 * 0.88 = 228.8 Var(Y) = n * p * (1 - p) = 260 * 0.88 * 0.12 = 27.7
|
|
122
|
+
2. Найдем математическое ожидание и дисперсию Yⱼ (выручка за одну книгу): E(Yⱼ) = 500 * 0.67 + 300 * 0.3 + 100 * 0.03 = 335 + 90 + 3 = 428 E(Yⱼ²) = 500² * 0.67 + 300² * 0.3 + 100² * 0.03 = 167500 + 27000 + 300 = 194800 Var(Yⱼ) = E(Yⱼ²) - (E(Yⱼ))² = 194800 - 428² = 194800 - 183184 = 11616
|
|
123
|
+
3. Суммарный ежедневный доход можно представить как S = Σⱼ Yⱼ, где сумма берется по всем проданным книгам. Так как количество проданных книг Y является случайной величиной, то для нахождения E(S) и Var(S) воспользуемся формулами для случайного числа слагаемых: E(S) = E(Y) * E(Yⱼ) = 228.8 * 428 = 97966.4 Var(S) = E(Y) * Var(Yⱼ) + Var(Y) * (E(Yⱼ))² = 228.8 * 11616 + 27.7 * 428² = 2660524.8 + 5085353.2 = 7745878
|
|
124
|
+
"""
|
|
125
|
+
if var == 10:
|
|
126
|
+
text = """Найдите верхнюю 5-процентную точку распределения Х, если Y принимает с равной вероятностью значения 20 и 80, а при фиксированном Y случайная величина Х равномерно распределена на отрезке [0; у].
|
|
127
|
+
Решение:
|
|
128
|
+
Верхняя 5-процентная точка для X означает такое значение x₀.₀₅, что P(X > x₀.₀₅) = 0.05.
|
|
129
|
+
P(X > x₀.₀₅) = P(X > x₀.₀₅ | Y=20)P(Y=20) + P(X > x₀.₀₅ | Y=80)P(Y=80)
|
|
130
|
+
Так как Y принимает значения 20 и 80 с равной вероятностью, P(Y=20) = P(Y=80) = 1/2.
|
|
131
|
+
X равномерно распределен на [0; Y], поэтому:
|
|
132
|
+
P(X > x₀.₀₅ | Y=20) = (20 - x₀.₀₅)/20, если 0 ≤ x₀.₀₅ ≤ 20; 0, если x₀.₀₅ > 20. P(X > x₀.₀₅ | Y=80) = (80 - x₀.₀₅)/80, если 0 ≤ x₀.₀₅ ≤ 80; 0, если x₀.₀₅ > 80.
|
|
133
|
+
Подставим в исходное уравнение:
|
|
134
|
+
0.05 = (1/2)(20 - x₀.₀₅)/20 + (1/2)(80 - x₀.₀₅)/80
|
|
135
|
+
0.1 = (20 - x₀.₀₅)/20 + (80 - x₀.₀₅)/80
|
|
136
|
+
8 + 0.1x₀.₀₅ = 20 - x₀.₀₅ + 80 - x₀.₀₅
|
|
137
|
+
8 + 0.1x₀.₀₅ = 100 - 2x₀.₀₅
|
|
138
|
+
2.1x₀.₀₅ = 92
|
|
139
|
+
x₀.₀₅ ≈ 43.81
|
|
140
|
+
"""
|
|
141
|
+
if var == 11:
|
|
142
|
+
text = """
|
|
143
|
+
Количество опоздавших на самолет пассажиров для каждого рейса описывается случайной величиной Х, распределенной по закону Пуассона с параметром λ = 6.
|
|
144
|
+
При этом стоимость билета, который не подлежит возврату, описывается нормально распределенной случайной величиной У с параметрами m=5250 и o=140.7.
|
|
145
|
+
Вычислите среднюю сумму стоимости пропавших билетов и её
|
|
146
|
+
среднеквадратическое отклонение, приходящиеся на каждый рейс
|
|
147
|
+
Дано:
|
|
148
|
+
Количество опоздавших пассажиров X ~ Poisson(λ = 6)
|
|
149
|
+
Стоимость билета Y ~ N(5250, 140.7^2)
|
|
150
|
+
|
|
151
|
+
Решение:
|
|
152
|
+
|
|
153
|
+
1. Математическое ожидание суммы стоимости пропавших билетов:
|
|
154
|
+
E(S) = E(X) * E(Y)
|
|
155
|
+
E(X) = λ = 6
|
|
156
|
+
E(Y) = 5250
|
|
157
|
+
E(S) = 6 * 5250 = 31500
|
|
158
|
+
|
|
159
|
+
2. Дисперсия суммы стоимости пропавших билетов:
|
|
160
|
+
Var(S) = E(X) * Var(Y) + (E(Y))^2 * Var(X)
|
|
161
|
+
Var(X) = λ = 6
|
|
162
|
+
Var(Y) = 140.7^2 = 19796.49
|
|
163
|
+
Var(S) = 6 * 19796.49 + 5250^2 * 6 = 165493778,94
|
|
164
|
+
|
|
165
|
+
3. Среднеквадратическое отклонение суммы стоимости пропавших билетов:
|
|
166
|
+
σ(S) = sqrt(Var(S)) = sqrt(165493778,94) ≈ 12864,43854
|
|
167
|
+
|
|
168
|
+
Ответ:
|
|
169
|
+
Средняя сумма стоимости пропавших билетов: 31500
|
|
170
|
+
Среднеквадратическое отклонение: 12864,43854
|
|
171
|
+
"""
|
|
172
|
+
if var == 12:
|
|
173
|
+
text = """В первом броске участвуют 105 несимметричных монет. Во втором броске участвуют только те монеты, на которых в первом броске выпал "орел". Известно, что вероятность выпадения "орла" для данных несимметричных монет равна 0,65. Найдите: 1) математическое ожидание числа "орлов", выпавших во втором броске; 2) дисперсию условного математического ожидания числа "орлов", выпавших во втором броске, относительно числа "орлов", выпавших в первом броске.
|
|
174
|
+
Решение:
|
|
175
|
+
1. Математическое ожидание числа "орлов" в первом броске (X):
|
|
176
|
+
X ~ Bin(n=105, p=0.65)
|
|
177
|
+
E(X) = n * p = 105 * 0.65 = 68.25
|
|
178
|
+
2. Условное математическое ожидание числа "орлов" во втором броске (Y) при заданном числе "орлов" в первом броске (x):
|
|
179
|
+
Y|X=x ~ Bin(n=x, p=0.65)
|
|
180
|
+
E(Y|X=x) = x * 0.65
|
|
181
|
+
3. Математическое ожидание числа "орлов" во втором броске (E(Y)):
|
|
182
|
+
E(Y) = E[E(Y|X)] = E[0.65X] = 0.65 * E(X) = 0.65 * 68.25 = 44.3625
|
|
183
|
+
4. Дисперсия условного математического ожидания числа "орлов" во втором броске относительно числа "орлов" в первом броске (Var[E(Y|X)]):
|
|
184
|
+
Var[E(Y|X)] = Var[0.65X] = 0.65² * Var(X)
|
|
185
|
+
Var(X) = n * p * (1 - p) = 105 * 0.65 * 0.35 = 23.9625
|
|
186
|
+
Var[E(Y|X)] = 0.65² * 23.9625 ≈ 10.13
|
|
187
|
+
Ответ: 1) 44.3625; 2) 10.13
|
|
188
|
+
"""
|
|
189
|
+
if var == 13:
|
|
190
|
+
text = """Ежедневное количество нарушений правил дорожного движения в некоторой местности описывается случайной величиной X, распределенной по биномиальному закону с параметрами n = 400 и p = 0.59. А сумма штрафа (в рублях) каждого нарушения описывается случайной величиной Y, которая принимает значение 500 рублей с вероятностью 0.62, 1000 рублей с вероятностью 0.12 и 3000 рублей в остальных случаях. Вычислите ежедневную среднюю сумму выписанных штрафов и ее среднеквадратическое отклонение.
|
|
191
|
+
Решение:
|
|
192
|
+
1. Характеристики случайной величины X (количество нарушений):
|
|
193
|
+
o Математическое ожидание (среднее): E(X) = n * p = 400 * 0.59 = 236
|
|
194
|
+
o Дисперсия: Var(X) = n * p * (1 - p) = 400 * 0.59 * (1 - 0.59) = 400 * 0.59 * 0.41 = 96.76
|
|
195
|
+
2. Характеристики случайной величины Y (сумма штрафа за одно нарушение):
|
|
196
|
+
o Вероятность для 3000 рублей: 1 - 0.62 - 0.12 = 0.26
|
|
197
|
+
o Математическое ожидание (средний штраф за одно нарушение): E(Y) = 500 * 0.62 + 1000 * 0.12 + 3000 * 0.26 = 310 + 120 + 780 = 1210 рублей
|
|
198
|
+
o Для вычисления дисперсии Var(Y) сначала найдем E(Y²): E(Y²) = 500² * 0.62 + 1000² * 0.12 + 3000² * 0.26 = 155000 + 120000 + 2340000 = 2615000
|
|
199
|
+
o Теперь дисперсия: Var(Y) = E(Y²) - (E(Y))² = 2615000 - 1210² = 2615000 - 1464100 = 1150900
|
|
200
|
+
3. Ежедневная сумма выписанных штрафов (S):
|
|
201
|
+
Ежедневная сумма штрафов S равна сумме штрафов за каждое нарушение. Поскольку количество нарушений X — случайная величина, а штраф за каждое нарушение Y — тоже случайная величина, и они независимы, математическое ожидание и дисперсию S можно найти по формулам для случайного числа слагаемых:
|
|
202
|
+
o Математическое ожидание (средняя ежедневная сумма штрафов): E(S) = E(X) * E(Y) = 236 * 1210 = 285560 рублей
|
|
203
|
+
o Дисперсия ежедневной суммы штрафов: Var(S) = E(X) * Var(Y) + Var(X) * (E(Y))² = 236 * 1150900 + 96.76 * 1210² = 271512400 + 141682316 = 413194716
|
|
204
|
+
4. Среднеквадратическое отклонение ежедневной суммы штрафов:
|
|
205
|
+
SD(S) = √Var(S) = √413194716 ≈ 20327.2 рублей
|
|
206
|
+
Ответ:
|
|
207
|
+
• Ежедневная средняя сумма выписанных штрафов: 285560 рублей.
|
|
208
|
+
• Среднеквадратическое отклонение ежедневной суммы выписанных штрафов: ≈ 20327.2 рубля.
|
|
209
|
+
"""
|
|
210
|
+
if var == 14:
|
|
211
|
+
text = """Игральная кость подбрасывается до тех пор, пока не выпадет 6 очков. Пусть S – сумма очков во всех бросках, кроме последнего (S=0, если 6 очков выпало при первом броске). Найдите E(S) и D(S)
|
|
212
|
+
Определение случайной величины S
|
|
213
|
+
S - это сумма значений, выпавших на игральной кости до выпадения шестерки. Это означает, что последний бросок (шестерка) в S не учитывается.
|
|
214
|
+
2. Вероятность выпадения любого числа (кроме 6)
|
|
215
|
+
Для стандартной шестигранной кости вероятность выпадения любого числа от 1 до 5 равна 1/6. Вероятность выпадения 6 равна также 1/6.
|
|
216
|
+
3. Математическое ожидание одного броска (до выпадения 6)
|
|
217
|
+
Обозначим X - результат одного броска до выпадения 6. Тогда математическое ожидание E(X) будет равно:
|
|
218
|
+
E(X) = (1/5) * (1 + 2 + 3 + 4 + 5) = (1/5) * 15 = 3
|
|
219
|
+
Почему мы делим на 5? Потому что мы рассматриваем только случаи, когда 6 не выпала, то есть 5 возможных исходов.
|
|
220
|
+
4. Математическое ожидание количества бросков до выпадения 6
|
|
221
|
+
Обозначим N - количество бросков до выпадения 6. N имеет геометрическое распределение с параметром p = 1/6 (вероятность успеха, то есть выпадения 6). Математическое ожидание геометрического распределения равно:
|
|
222
|
+
E(N) = 1/p = 1/(1/6) = 6
|
|
223
|
+
Однако, поскольку мы считаем количество бросков до выпадения 6, мы должны вычесть 1. Таким образом, среднее количество бросков до выпадения 6 равно 5.
|
|
224
|
+
5. Математическое ожидание суммы S
|
|
225
|
+
Теперь мы можем найти математическое ожидание суммы S. Поскольку каждый бросок (до выпадения 6) в среднем дает 3, а количество таких бросков в среднем равно 5, то:
|
|
226
|
+
E(S) = E(X) * E(N) = 3 * 5 = 15
|
|
227
|
+
6. Дисперсия одного броска (до выпадения 6)
|
|
228
|
+
Дисперсия одного броска X (до выпадения 6) вычисляется следующим образом:
|
|
229
|
+
E(X²) = (1/5) * (1² + 2² + 3² + 4² + 5²) = (1/5) * 55 = 11
|
|
230
|
+
D(X) = E(X²) - (E(X))² = 11 - 3² = 11 - 9 = 2
|
|
231
|
+
7. Дисперсия количества бросков до выпадения 6
|
|
232
|
+
Дисперсия геометрического распределения равна:
|
|
233
|
+
D(N) = (1-p)/p² = (1 - 1/6) / (1/6)² = (5/6) / (1/36) = 30
|
|
234
|
+
8. Дисперсия суммы S
|
|
235
|
+
Поскольку броски независимы, дисперсия суммы равна сумме дисперсий:
|
|
236
|
+
D(S) = E(N) * D(X) + (E(X))² * D(N) = 5 * 2 + 3² * 30 = 10 + 9 * 30 = 10 + 270 = 280
|
|
237
|
+
"""
|
|
238
|
+
if var == 15:
|
|
239
|
+
text = """Средний ущерб от одного пожара составляет 5,9 млн. руб. Предполагается, что ущерб распределен по показательному закону, а число пожаров за год — по закону Пуассона. Также известно, что за 10 лет в среднем происходит 15 пожаров. Найдите:
|
|
240
|
+
1. математическое ожидание суммарного ущерба от всех пожаров за один год:
|
|
241
|
+
2. стандартное отклонение суммарного ущерба от пожаров за год.
|
|
242
|
+
Решение:
|
|
243
|
+
Пусть X — число пожаров за год, Y — ущерб от одного пожара.
|
|
244
|
+
• Из условия следует, что E(Y) = 5.9 млн. руб. Показательное распределение имеет вид f(y) = λe^(-λy), где E(Y) = 1/λ. Следовательно, λ = 1/5.9.
|
|
245
|
+
• За 10 лет происходит 15 пожаров, значит, в среднем за год происходит E(X) = 15/10 = 1.5 пожара. Так как число пожаров распределено по закону Пуассона, Var(X) = E(X) = 1.5.
|
|
246
|
+
Суммарный ущерб S = Σᵢ Yᵢ, где i = 1, ..., X.
|
|
247
|
+
1. E(S) = E(X) * E(Y) = 1.5 * 5.9 = 8.85 млн. руб.
|
|
248
|
+
2. Var(S) = E(X) * Var(Y) + (E(Y))² * Var(X)
|
|
249
|
+
Для показательного распределения Var(Y) = 1/λ² = (5.9)² = 34.81.
|
|
250
|
+
Var(S) = 1.5 * 34.81 + (5.9)² * 1.5 = 52.215 + 52.215 = 104.43 млн. руб.²
|
|
251
|
+
Стандартное отклонение SD(S) = √Var(S) = √104.43 ≈ 10.22 млн. руб.
|
|
252
|
+
Ответ: 1) 8.85 млн. руб.; 2) ≈ 10.22 млн. руб.
|
|
253
|
+
"""
|
|
254
|
+
if var == 16:
|
|
255
|
+
text = """Случайные величины X и Y независимы и с равной вероятностью принимают значения: 1, 2, 3, 4, 5. Найдите E(X|X>Y) и D(X|X>Y).
|
|
256
|
+
Решение:
|
|
257
|
+
1. Найдем P(X>Y):
|
|
258
|
+
Так как X и Y независимы и равновероятны, всего возможно 5 * 5 = 25 равновероятных пар (x, y). Пары, где x > y: (2,1), (3,1), (3,2), (4,1), (4,2), (4,3), (5,1), (5,2), (5,3), (5,4). Всего 10 таких пар.
|
|
259
|
+
P(X>Y) = 10/25 = 2/5 = 0.4
|
|
260
|
+
2. Найдем условное распределение X при условии X>Y:
|
|
261
|
+
Возможные значения X при X>Y: 2, 3, 4, 5.
|
|
262
|
+
o P(X=2|X>Y) = P(X=2, X>Y) / P(X>Y) = P(X=2, Y=1) / P(X>Y) = (1/25) / (2/5) = 1/10
|
|
263
|
+
o P(X=3|X>Y) = P(X=3, X>Y) / P(X>Y) = P(X=3, Y=1 или Y=2) / P(X>Y) = (2/25) / (2/5) = 2/10 = 1/5
|
|
264
|
+
o P(X=4|X>Y) = P(X=4, X>Y) / P(X>Y) = P(X=4, Y=1, 2 или 3) / P(X>Y) = (3/25) / (2/5) = 3/10
|
|
265
|
+
o P(X=5|X>Y) = P(X=5, X>Y) / P(X>Y) = P(X=5, Y=1, 2, 3 или 4) / P(X>Y) = (4/25) / (2/5) = 4/10 = 2/5
|
|
266
|
+
3. Найдем E(X|X>Y):
|
|
267
|
+
E(X|X>Y) = 2*(1/10) + 3*(1/5) + 4*(3/10) + 5*(2/5) = 2/10 + 6/10 + 12/10 + 20/10 = 40/10 = 4
|
|
268
|
+
4. Найдем E(X²|X>Y):
|
|
269
|
+
E(X²|X>Y) = 2²*(1/10) + 3²*(1/5) + 4²*(3/10) + 5²*(2/5) = 4/10 + 18/10 + 48/10 + 50/10 = 120/10 = 12
|
|
270
|
+
5. Найдем D(X|X>Y):
|
|
271
|
+
D(X|X>Y) = E(X²|X>Y) - (E(X|X>Y))² = 12 - 4² = 12 - 16 = -4 (ошибка в расчетах выше)
|
|
272
|
+
Правильно: E(X|X>Y) = 2 * (1/10) + 3 * (2/10) + 4 * (3/10) + 5 * (4/10) = (2 + 6 + 12 + 20) / 10 = 40/10 = 4
|
|
273
|
+
E(X²|X>Y) = 4 * (1/10) + 9 * (2/10) + 16 * (3/10) + 25 * (4/10) = (4 + 18 + 48 + 100) / 10 = 170/10 = 17
|
|
274
|
+
D(X|X>Y) = E(X²|X>Y) - (E(X|X>Y))² = 17 - 4² = 17 - 16 = 1
|
|
275
|
+
Ответ: E(X|X>Y) = 4, D(X|X>Y) = 1.
|
|
276
|
+
"""
|
|
277
|
+
if var == 17:
|
|
278
|
+
text = """Количество опоздавших на самолет пассажиров для каждого рейса описывается случайной величиной X, распределенной по закону Пуассона с параметром λ = 5. При этом стоимость билета, который не подлежит возврату, описывается нормально распределенной случайной величиной Y с параметрами m = 4800 и σ = 270.2. Вычислите среднюю сумму стоимости пропавших билетов и ее среднеквадратическое отклонение, приходящиеся на каждый рейс.
|
|
279
|
+
Решение:
|
|
280
|
+
Пусть X – количество опоздавших пассажиров, Y – стоимость одного билета. Суммарная стоимость пропавших билетов S = Σᵢ Yᵢ, где i = 1,..., X.
|
|
281
|
+
1. Математическое ожидание X: E(X) = λ = 5
|
|
282
|
+
2. Математическое ожидание Y: E(Y) = m = 4800
|
|
283
|
+
3. Дисперсия Y: Var(Y) = σ² = (270.2)² = 73008.04
|
|
284
|
+
Средняя сумма стоимости пропавших билетов:
|
|
285
|
+
E(S) = E(X) * E(Y) = 5 * 4800 = 24000
|
|
286
|
+
Дисперсия суммы стоимости пропавших билетов:
|
|
287
|
+
Var(S) = E(X) * Var(Y) + (E(Y))² * Var(X) = 5 * 73008.04 + (4800)² * 5 = 365040.2 + 115200000 = 115565040.2
|
|
288
|
+
Среднеквадратическое отклонение:
|
|
289
|
+
SD(S) = √Var(S) = √115565040.2 ≈ 10750.12
|
|
290
|
+
Ответ: Средняя сумма стоимости пропавших билетов: 24000. Среднеквадратическое отклонение: ≈ 10750.12
|
|
291
|
+
"""
|
|
292
|
+
if var == 18:
|
|
293
|
+
text = """Ежедневное количество нарушений правил дорожного движения в некоторой местности описывается случайной величиной X, распределенной по биномиальному закону с параметрами n = 700 и p = 0.69. А сумма штрафа (в рублях) каждого нарушения описывается случайной величиной Y, которая принимает значение 500 рублей с вероятностью 0.59, 1000 рублей с вероятностью 0.14 и 3000 рублей в остальных случаях. Вычислите ежедневную среднюю сумму выписанных штрафов и ее среднеквадратическое отклонение.
|
|
294
|
+
Решение:
|
|
295
|
+
1. Характеристики X (количество нарушений):
|
|
296
|
+
o E(X) = n * p = 700 * 0.69 = 483
|
|
297
|
+
o Var(X) = n * p * (1 - p) = 700 * 0.69 * 0.31 = 149.73
|
|
298
|
+
2. Характеристики Y (сумма штрафа):
|
|
299
|
+
o Вероятность штрафа 3000 рублей: 1 - 0.59 - 0.14 = 0.27
|
|
300
|
+
o E(Y) = 500 * 0.59 + 1000 * 0.14 + 3000 * 0.27 = 295 + 140 + 810 = 1245
|
|
301
|
+
o E(Y²) = 500² * 0.59 + 1000² * 0.14 + 3000² * 0.27 = 147500 + 140000 + 2430000 = 2717500
|
|
302
|
+
o Var(Y) = E(Y²) - (E(Y))² = 2717500 - 1245² = 2717500 - 1550025 = 1167475
|
|
303
|
+
3. Ежедневная сумма штрафов S:
|
|
304
|
+
o E(S) = E(X) * E(Y) = 483 * 1245 = 601335
|
|
305
|
+
o Var(S) = E(X) * Var(Y) + Var(X) * (E(Y))² = 483 * 1167475 + 149.73 * 1245² = 564177325 + 232230386.25 = 796407711.25
|
|
306
|
+
o SD(S) = √Var(S) = √796407711.25 ≈ 28220.7
|
|
307
|
+
Ответ: Средняя сумма штрафов: 601335 рублей. Среднеквадратическое отклонение: ≈ 28220.7 рублей.
|
|
308
|
+
"""
|
|
309
|
+
if var == 19:
|
|
310
|
+
text = """Дискретный случайный вектор (X, Y) имеет распределение:
|
|
311
|
+
X=-1 X=0 X=1
|
|
312
|
+
Y=0 1/4 1/8
|
|
313
|
+
Y=1 1/6 1/6 1/6
|
|
314
|
+
Найдите условное математическое ожидание Z = E(Y|X+Y=1).
|
|
315
|
+
Решение:
|
|
316
|
+
Сначала заполним пропущенную вероятность. Сумма всех вероятностей должна быть равна 1:
|
|
317
|
+
1/4 + 1/8 + 1/6 + 1/6 + 1/6 + P(X=1, Y=0) = 1 P(X=1, Y=0) = 1 - (6+3+4+4+4)/24 = 1-21/24= 3/24 = 1/8
|
|
318
|
+
Теперь найдем условное распределение при X+Y=1. Возможные варианты: (X=-1, Y=2), (X=0, Y=1) и (X=1,Y=0). Из таблицы видно, что X+Y=1 только для случаев (0,1) и (1,0).
|
|
319
|
+
P(X=0, Y=1) = 1/6 P(X=1, Y=0) = 1/8
|
|
320
|
+
P(X+Y=1) = P(X=0, Y=1) + P(X=1, Y=0) = 1/6 + 1/8 = (4+3)/24 = 7/24
|
|
321
|
+
Теперь найдем условные вероятности:
|
|
322
|
+
P(Y=1|X+Y=1) = P(X=0, Y=1) / P(X+Y=1) = (1/6) / (7/24) = 4/7 P(Y=0|X+Y=1) = P(X=1, Y=0) / P(X+Y=1) = (1/8) / (7/24) = 3/7
|
|
323
|
+
Тогда Z = E(Y|X+Y=1) = 1 * P(Y=1|X+Y=1) + 0 * P(Y=0|X+Y=1) = 1 * (4/7) + 0 * (3/7) = 4/7
|
|
324
|
+
Ответ: E(Y|X+Y=1) = 4/7.
|
|
325
|
+
"""
|
|
326
|
+
if var == 20:
|
|
327
|
+
text = """Плотность распределения случайного вектора (X, Y) имеет вид:
|
|
328
|
+
fₓ,y(x,y) = (1/π) * exp(-(5/2)x² - 18x - 36 - xy - 6y - (1/2)y²)
|
|
329
|
+
Найдите условное математическое ожидание E(Y|X=x).
|
|
330
|
+
Решение:
|
|
331
|
+
Для нахождения условного математического ожидания E(Y|X=x) необходимо сначала найти условную плотность f(y|x) = f(x,y) / f(x).
|
|
332
|
+
Преобразуем выражение в показателе экспоненты:
|
|
333
|
+
-(5/2)x² - 18x - 36 - xy - 6y - (1/2)y² = -(1/2)(5x² + 36x + 72 + 2xy + 12y + y²) = -(1/2)[y² + (2x + 12)y + (5x² + 36x + 72)]
|
|
334
|
+
Выделим полный квадрат относительно y:
|
|
335
|
+
-(1/2)[(y + x + 6)² - (x+6)² + 5x² + 36x + 72] = -(1/2)[(y + x + 6)² - x² - 12x - 36 + 5x² + 36x + 72] = -(1/2)[(y + x + 6)² + 4x² + 24x + 36] = -(1/2)(y + x + 6)² - 2(x² + 6x + 9) = -(1/2)(y + x + 6)² - 2(x + 3)²
|
|
336
|
+
Таким образом,
|
|
337
|
+
fₓ,y(x,y) = (1/π) * exp(-(1/2)(y + x + 6)² - 2(x + 3)²) = (1/√(2π) * exp(-2(x + 3)²)) * (1/√(2π(1/2))) * exp(-(1/2)(y + x + 6)²)
|
|
338
|
+
Отсюда видно, что:
|
|
339
|
+
f(x) = (1/√(2π(1/4))) * exp(-2(x+3)²) – плотность нормального распределения N(-3, 1/4) f(y|x) = (1/√(2π(1))) * exp(-(1/2)(y + x + 6)²) – плотность нормального распределения N(-x-6, 1)
|
|
340
|
+
Следовательно, E(Y|X=x) = -x - 6.
|
|
341
|
+
Ответ: E(Y|X=x) = -x - 6.
|
|
342
|
+
"""
|
|
343
|
+
if var == 21:
|
|
344
|
+
text = """В первом броске участвуют 105 несимметричных монет. Во втором броске участвуют только те монеты, на которых в первом броске выпал "орел". Известно, что вероятность выпадения "орла" для данных несимметричных монет равна 0,65. Найдите: 1) математическое ожидание числа "орлов", выпавших во втором броске; 2) дисперсию условного математического ожидания числа "орлов", выпавших во втором броске, относительно числа "орлов", выпавших в первом броске.
|
|
345
|
+
Решение:
|
|
346
|
+
Пусть X — число "орлов" в первом броске, Y — число "орлов" во втором броске.
|
|
347
|
+
1. Математическое ожидание числа "орлов" в первом броске:
|
|
348
|
+
E(X) = n * p = 105 * 0.65 = 68.25
|
|
349
|
+
2. Условное математическое ожидание числа "орлов" во втором броске при условии X=x:
|
|
350
|
+
E(Y|X=x) = x * p = 0.65x
|
|
351
|
+
3. Математическое ожидание числа "орлов" во втором броске (полное математическое ожидание):
|
|
352
|
+
E(Y) = E[E(Y|X)] = E(0.65X) = 0.65E(X) = 0.65 * 68.25 = 44.3625
|
|
353
|
+
4. Дисперсия числа "орлов" в первом броске:
|
|
354
|
+
Var(X) = n * p * (1 - p) = 105 * 0.65 * 0.35 = 23.9625
|
|
355
|
+
5. Дисперсия условного математического ожидания числа "орлов" во втором броске относительно числа "орлов", выпавших в первом броске:
|
|
356
|
+
Var[E(Y|X)] = Var(0.65X) = 0.65² * Var(X) = 0.4225 * 23.9625 = 10.12515625
|
|
357
|
+
Ответ: 1) E(Y) = 44.3625; 2) Var[E(Y|X)] ≈ 10.13
|
|
358
|
+
"""
|
|
359
|
+
if var == 22:
|
|
360
|
+
text = """
|
|
361
|
+
Количество опоздавших на самолет пассажиров для каждого рейса описывается случайной величиной Х, распределенной по закону Пуассона с параметром λ = 7.
|
|
362
|
+
При этом стоимость билета, который не подлежит возврату, описывается нормально
|
|
363
|
+
распределенной случайной величиной У с параметрами m=6150 и o=133.7.
|
|
364
|
+
Вычислите среднюю сумму стоимости пропавших билетов и её
|
|
365
|
+
среднеквадратическое отклонение, приходящиеся на каждый рейс
|
|
366
|
+
Дано:
|
|
367
|
+
Количество опоздавших пассажиров X ~ Poisson(λ = 7)
|
|
368
|
+
Стоимость билета Y ~ N(6150, 133.7^2)
|
|
369
|
+
Решение:
|
|
370
|
+
1. Математическое ожидание суммы стоимости пропавших билетов:
|
|
371
|
+
E(S) = E(X) * E(Y)
|
|
372
|
+
E(X) = λ = 7
|
|
373
|
+
E(Y) = 6150
|
|
374
|
+
E(S) = 7 * 6150 = 43050
|
|
375
|
+
2. Дисперсия суммы стоимости пропавших билетов:
|
|
376
|
+
Var(S) = E(X) * Var(Y) + (E(Y))^2 * Var(X)
|
|
377
|
+
Var(X) = λ = 7
|
|
378
|
+
Var(Y) = 133.7^2 = 17875,69
|
|
379
|
+
Var(S) = 7 * 17875,69 + 6150^2 * 7 = 264882629,83
|
|
380
|
+
3. Среднеквадратическое отклонение суммы стоимости пропавших билетов:
|
|
381
|
+
σ(S) = sqrt(Var(S)) = sqrt(264882629,83) ≈ 16275,2152
|
|
382
|
+
Ответ:
|
|
383
|
+
Средняя сумма стоимости пропавших билетов: 43050
|
|
384
|
+
Среднеквадратическое отклонение: 16275,2152
|
|
385
|
+
"""
|
|
386
|
+
if var == 23:
|
|
387
|
+
text = """Средний ущерб от одного пожара составляет 5,9 млн. руб. Предполагается, что ущерб распределен по показательному закону, а число пожаров за год по закону Пуассона. Также известно, что за 10 лет в среднем происходит 15 пожаров. Найдите: 1) математическое ожидание суммарного ущерба от всех пожаров за один год; 2) стандартное отклонение суммарного ущерба от пожаров за год.
|
|
388
|
+
Решение:
|
|
389
|
+
Пусть X — число пожаров за год, Y — ущерб от одного пожара.
|
|
390
|
+
1. Параметр λ для распределения Пуассона (число пожаров в год): За 10 лет 15 пожаров, значит, в год λ = 15 / 10 = 1.5 пожара. E(X) = λ = 1.5
|
|
391
|
+
2. Математическое ожидание ущерба от одного пожара (Y): Для показательного распределения E(Y) = 1/λ_экспон, где λ_экспон — параметр показательного распределения. По условию, E(Y) = 5.9 млн. руб.
|
|
392
|
+
3. Математическое ожидание суммарного ущерба (S): S = Σᵢ Yᵢ (где i от 1 до X) E(S) = E(X) * E(Y) = 1.5 * 5.9 = 8.85 млн. руб.
|
|
393
|
+
4. Дисперсия ущерба от одного пожара (Y): Для показательного распределения Var(Y) = 1/λ_экспон² = (E(Y))² = 5.9² = 34.81 млн. руб.²
|
|
394
|
+
5. Дисперсия суммарного ущерба (S): Var(S) = E(X) * Var(Y) + (E(Y))² * Var(X) Var(X) = λ = 1.5 (для распределения Пуассона) Var(S) = 1.5 * 34.81 + 5.9² * 1.5 = 52.215 + 52.215 = 104.43 млн. руб.²
|
|
395
|
+
6. Стандартное отклонение суммарного ущерба (SD(S)): SD(S) = √Var(S) = √104.43 ≈ 10.22 млн. руб.
|
|
396
|
+
Ответ: 1) E(S) = 8.85 млн. руб.; 2) SD(S) ≈ 10.22 млн. руб.
|
|
397
|
+
"""
|
|
398
|
+
if var == 24:
|
|
399
|
+
text = """Игральная кость подбрасывается до тех пор, пока не выпадет 6 очков. Пусть S – сумма очков во всех бросках, кроме последнего (S=0, если 6 очков выпало при первом броске). Найдите E(S) и D(S)
|
|
400
|
+
Определение случайной величины S
|
|
401
|
+
S - это сумма значений, выпавших на игральной кости до выпадения шестерки. Это означает, что последний бросок (шестерка) в S не учитывается.
|
|
402
|
+
2. Вероятность выпадения любого числа (кроме 6)
|
|
403
|
+
Для стандартной шестигранной кости вероятность выпадения любого числа от 1 до 5 равна 1/6. Вероятность выпадения 6 равна также 1/6.
|
|
404
|
+
3. Математическое ожидание одного броска (до выпадения 6)
|
|
405
|
+
Обозначим X - результат одного броска до выпадения 6. Тогда математическое ожидание E(X) будет равно:
|
|
406
|
+
E(X) = (1/5) * (1 + 2 + 3 + 4 + 5) = (1/5) * 15 = 3
|
|
407
|
+
Почему мы делим на 5? Потому что мы рассматриваем только случаи, когда 6 не выпала, то есть 5 возможных исходов.
|
|
408
|
+
4. Математическое ожидание количества бросков до выпадения 6
|
|
409
|
+
Обозначим N - количество бросков до выпадения 6. N имеет геометрическое распределение с параметром p = 1/6 (вероятность успеха, то есть выпадения 6). Математическое ожидание геометрического распределения равно:
|
|
410
|
+
E(N) = 1/p = 1/(1/6) = 6
|
|
411
|
+
Однако, поскольку мы считаем количество бросков до выпадения 6, мы должны вычесть 1. Таким образом, среднее количество бросков до выпадения 6 равно 5.
|
|
412
|
+
5. Математическое ожидание суммы S
|
|
413
|
+
Теперь мы можем найти математическое ожидание суммы S. Поскольку каждый бросок (до выпадения 6) в среднем дает 3, а количество таких бросков в среднем равно 5, то:
|
|
414
|
+
E(S) = E(X) * E(N) = 3 * 5 = 15
|
|
415
|
+
6. Дисперсия одного броска (до выпадения 6)
|
|
416
|
+
Дисперсия одного броска X (до выпадения 6) вычисляется следующим образом:
|
|
417
|
+
E(X²) = (1/5) * (1² + 2² + 3² + 4² + 5²) = (1/5) * 55 = 11
|
|
418
|
+
D(X) = E(X²) - (E(X))² = 11 - 3² = 11 - 9 = 2
|
|
419
|
+
7. Дисперсия количества бросков до выпадения 6
|
|
420
|
+
Дисперсия геометрического распределения равна:
|
|
421
|
+
D(N) = (1-p)/p² = (1 - 1/6) / (1/6)² = (5/6) / (1/36) = 30
|
|
422
|
+
8. Дисперсия суммы S
|
|
423
|
+
Поскольку броски независимы, дисперсия суммы равна сумме дисперсий:
|
|
424
|
+
D(S) = E(N) * D(X) + (E(X))² * D(N) = 5 * 2 + 3² * 30 = 10 + 9 * 30 = 10 + 270 = 280
|
|
425
|
+
"""
|
|
426
|
+
if var == 25:
|
|
427
|
+
text = """Плотность распределения случайного вектора (X, Y) имеет вид:
|
|
428
|
+
fₓ,y(x,y) = (1/π) * exp(-2x² - 2xy - y²)
|
|
429
|
+
Найдите условное математическое ожидание E(Y|X=x).
|
|
430
|
+
Решение:
|
|
431
|
+
Чтобы найти условное математическое ожидание E(Y|X=x), нужно найти условную плотность f(y|x) = f(x,y) / f(x).
|
|
432
|
+
Преобразуем выражение в показателе экспоненты:
|
|
433
|
+
-2x² - 2xy - y² = -(y² + 2xy + x² + x²) = -(y + x)² - x²
|
|
434
|
+
Таким образом,
|
|
435
|
+
fₓ,y(x,y) = (1/π) * exp(-(y + x)² - x²) = (1/√π * exp(-x²)) * (1/√π * exp(-(y + x)²))
|
|
436
|
+
Отсюда видно, что:
|
|
437
|
+
f(x) = (1/√π) * exp(-x²) - плотность нормального распределения N(0, 1/2) f(y|x) = (1/√π) * exp(-(y + x)²) - плотность нормального распределения N(-x, 1/2)
|
|
438
|
+
Следовательно, E(Y|X=x) = -x
|
|
439
|
+
"""
|
|
440
|
+
if var == 26:
|
|
441
|
+
text = """2.
|
|
442
|
+
|
|
443
|
+
"""
|
|
444
|
+
if var == 27:
|
|
445
|
+
text = """Ежедневное количество посетителей фитнес-клуба, купивших абонемент на посещение, описывается случайной величиной Х, распределенной по биномиальному закону с параметрами n = 700 и p = 0.74. А сумма абонемента (в рублях) каждого из посетителей описывается случайной величиной Y, распределенной по нормальному закону с параметрами m = 3500 и σ = 50. Вычислите ежедневную среднюю выручку фитнес-клуба и ее дисперсию.
|
|
446
|
+
Решение:
|
|
447
|
+
Пусть X – количество посетителей, Y – стоимость абонемента. Суммарная выручка S = X * Y.
|
|
448
|
+
1. Математическое ожидание X (количество посетителей): E(X) = n * p = 700 * 0.74 = 518
|
|
449
|
+
2. Математическое ожидание Y (стоимость абонемента): E(Y) = m = 3500
|
|
450
|
+
3. Ежедневная средняя выручка (математическое ожидание S): E(S) = E(X) * E(Y) = 518 * 3500 = 1813000
|
|
451
|
+
4. Дисперсия X: Var(X) = n * p * (1 - p) = 700 * 0.74 * (1 - 0.74) = 700 * 0.74 * 0.26 = 134.12
|
|
452
|
+
5. Дисперсия Y: Var(Y) = σ² = 50² = 2500
|
|
453
|
+
6. Дисперсия ежедневной выручки (дисперсия S): Var(S) = E(X)Var(Y) + (E(Y))²Var(X) = 518 * 2500 + 3500² * 134.12 = 1295000 + 1640000500=1641295500
|
|
454
|
+
Ответ: Ежедневная средняя выручка: 1 813 000 руб. Дисперсия ежедневной выручки: 1 641 295 500 руб²
|
|
455
|
+
"""
|
|
456
|
+
if var == 28:
|
|
457
|
+
text = """Плотность распределения (X,Y):
|
|
458
|
+
fₓ,y(x,y) = (1/(2π)) * exp(-x² - 10x - 26 - 3xy - 16y - (5/2)y²)
|
|
459
|
+
Найдите условное математическое ожидание Var(X|Y=y).
|
|
460
|
+
Решение:
|
|
461
|
+
1. Преобразуем показатель экспоненты:
|
|
462
|
+
-x² - 10x - 26 - 3xy - 16y - (5/2)y² = -(x² + (10 + 3y)x + (5/2)y² + 16y + 26)
|
|
463
|
+
Выделим полный квадрат относительно x:
|
|
464
|
+
-((x + (10 + 3y)/2)² - (10+3y)²/4 + (5/2)y² + 16y + 26) = -(x + (10 + 3y)/2)² + (100+60y+9y²)/4 - (10/4)y² - (32/4)y - (104/4) = -(x + 5 + (3/2)y)² + 25 + 15y + (9/4)y² - (5/2)y² - 8y - 26 = -(x + 5 + (3/2)y)² - (1/4)y² +7y - 1
|
|
465
|
+
2. Запишем плотность в виде произведения: fₓ,y(x,y) = (1/(2π)) * exp(-(x + 5 + (3/2)y)² - (1/4)y² +7y - 1) = (exp(-(1/4)y² +7y - 1)/√π)*(exp(-(x + 5 + (3/2)y)²)/√π)
|
|
466
|
+
3. Условная плотность f(x|y): f(x|y) = (exp(-(x + 5 + (3/2)y)²)/√π) - плотность нормального распределения с параметрами µ=-5-(3/2)y и σ²=1/2. Значит, X|Y=y ~ N(-5-(3/2)y; 1/2). Var(X|Y=y)=1/2
|
|
467
|
+
Ответ: Var(X|Y=y) = 1/2.
|
|
468
|
+
"""
|
|
469
|
+
if var == 29:
|
|
470
|
+
text = """Ежедневное количество жалоб в ресторане быстрого питания описывается случайной величиной X, распределенной по биномиальному закону с параметрами n = 20 и p = 0.15. А сумма штрафа (в рублях) от каждой жалобы посетителя управляющему описывается случайной величиной Y, которая принимает значение 50 рублей с вероятностью 0.62, 70 рублей с вероятностью 0.12 и 100 рублей в остальных случаях. Вычислите ежедневную среднюю сумму выписанных штрафов управляющему ресторана быстрого питания и ее среднеквадратическое отклонение.
|
|
471
|
+
Решение:
|
|
472
|
+
Пусть X – количество жалоб, Y – сумма штрафа за одну жалобу.
|
|
473
|
+
1. Математическое ожидание X (количество жалоб): E(X) = n * p = 20 * 0.15 = 3
|
|
474
|
+
2. Математическое ожидание Y (сумма штрафа): P(Y=100) = 1 - 0.62 - 0.12 = 0.26 E(Y) = 50 * 0.62 + 70 * 0.12 + 100 * 0.26 = 31 + 8.4 + 26 = 65.4
|
|
475
|
+
3. Ежедневная средняя сумма штрафов (математическое ожидание S, где S - суммарный штраф): E(S) = E(X) * E(Y) = 3 * 65.4 = 196.2
|
|
476
|
+
4. Дисперсия X: Var(X) = n * p * (1 - p) = 20 * 0.15 * (1 - 0.15) = 20 * 0.15 * 0.85 = 2.55
|
|
477
|
+
5. Дисперсия Y: E(Y²) = 50² * 0.62 + 70² * 0.12 + 100² * 0.26 = 1550 + 588 + 2600 = 4738 Var(Y) = E(Y²) - (E(Y))² = 4738 - 65.4² = 4738 - 4277.16 = 460.84
|
|
478
|
+
6. Дисперсия ежедневной суммы штрафов (Var(S)): Var(S) = E(X)Var(Y) + (E(Y))²Var(X) = 3 * 460.84 + 65.4² * 2.55 = 1382.52 + 10826.103 = 12208.623
|
|
479
|
+
7. Среднеквадратическое отклонение ежедневной суммы штрафов: SD(S) = √Var(S) = √12208.623 ≈ 110.49
|
|
480
|
+
Ответ: Ежедневная средняя сумма штрафов: 196.2 руб. Среднеквадратическое отклонение: ≈ 110.49 руб.
|
|
481
|
+
"""
|
|
482
|
+
if var == 30:
|
|
483
|
+
text = """Ежемесячные расходы владельцев А и В описываются случайным вектором (X, Y) с равномерным распределением в треугольнике D = {(x, y): x ≥ 0, y ≥ 0, x + y ≤ 1}. Найдите:
|
|
484
|
+
а) P(X + Y > 1/2); б) Плотность распределения Y; в) P(X < 1/3 | X > 1/2); г) Условную плотность и E(X|Y=1/2); д) E(E(X|Y)); е) Коэффициент корреляции X и Y.
|
|
485
|
+
Решение:
|
|
486
|
+
Площадь треугольника D равна 1/2. Плотность распределения в D равна f(x,y) = 2.
|
|
487
|
+
а) P(X + Y > 1/2): Эта область соответствует треугольнику, ограниченному x ≥ 0, y ≥ 0, x + y ≤ 1, исключая треугольник x ≥ 0, y ≥ 0, x + y ≤ 1/2. Площадь исключенного треугольника равна (1/2)(1/2)(1/2)=1/8. Значит, площадь нужной нам области равна 1/2 - 1/8 = 3/8. Тогда P(X + Y > 1/2) = (3/8) / (1/2) = 3/4 = 0.75.
|
|
488
|
+
б) Плотность распределения Y: fᵧ(y) = ∫₀^(1-y) 2 dx = 2(1-y), при 0 ≤ y ≤ 1.
|
|
489
|
+
в) P(X < 1/3 | X > 1/2): Это невозможно, так как если X > 1/2, то X никак не может быть меньше 1/3. P(X < 1/3 | X > 1/2) = 0.
|
|
490
|
+
г) Условная плотность f(x|y=1/2): f(x|y=1/2) = f(x, 1/2) / fᵧ(1/2) = 2 / (2(1-1/2)) = 2. x меняется от 0 до 1/2. E(X|Y=1/2) = ∫₀^(1/2) x * 2 dx = [x²]₀^(1/2) = 1/4.
|
|
491
|
+
д) E(E(X|Y)): E(E(X|Y)) = E(X) = ∫₀¹ ∫₀^(1-y) x * 2 dx dy=∫₀¹(1-y)²dy=1/3.
|
|
492
|
+
е) Коэффициент корреляции: E(X) = ∫₀¹ ∫₀^(1-y) x * 2 dx dy=1/3 E(Y)=1/3 E(XY) = ∫₀¹ ∫₀^(1-y) xy * 2 dx dy = 1/6 Var(X) = E(X²) - (E(X))² = ∫₀¹ ∫₀^(1-y) x² * 2 dx dy - (1/3)² = 1/6 - 1/9 = 1/18 Var(Y)=1/18 Cov(X,Y) = E(XY) - E(X)E(Y) = 1/6 - (1/3)(1/3) = -1/18 Корреляция ρ = Cov(X,Y) / (√Var(X)√Var(Y)) = (-1/18) / (1/18) = -1.
|
|
493
|
+
Ответ: а) 0.75; б) 2(1-y), 0 ≤ y ≤ 1; в) 0; г) f(x|y=1/2)=2, 0<=x<=1/2; E(X|Y=1/2)=1/4; д) 1/3; е) -1
|
|
494
|
+
"""
|
|
495
|
+
return text
|