Распределение Рэлея

редактировать
Рэлей
Функция плотности вероятности График Рэлеевского PDF .
Кумулятивная функция распределения График CDF Рэлея .
Параметрымасштаб: σ>0 {\ displaystyle \ sigma>0}\sigma>0
Поддержка x ∈ [0, ∞) {\ displaystyle x \ in [0, \ infty)}{\ displaystyle x \ in [0, \ infty)}
PDF x σ 2 e - x 2 / (2 σ 2) {\ displaystyle {\ frac {x} {\ sigma ^ {2}}} e ^ {- x ^ {2} / \ left (2 \ sigma ^ {2} \ right)}}{\ displaystyle {\ frac {x} {\ sigma ^ {2}}} e ^ {- x ^ {2} / \ left (2 \ sigma ^ {2} \ right)}}
CDF 1 - e - x 2 / (2 σ 2) {\ displaystyle 1-e ^ {- x ^ {2} / \ left (2 \ sigma ^ {2} \ right)}}{\ displaystyle 1-e ^ {- x ^ {2} / \ left (2 \ sigma ^ {2} \ right)}}
Квантиль Q (F; σ) = σ - 2 ln ⁡ (1 - F) {\ displaystyle Q (F; \ sigma) = \ sigma {\ sqrt {-2 \ ln (1-F)} }}{\ displaystyle Q (F; \ sigma) = \ sigma {\ sqrt {-2 \ ln (1-F)}}}
Среднее σ π 2 {\ displaystyle \ sigma {\ sqrt {\ frac {\ pi} {2}}}}\ sigma {\ sqrt {\ frac {\ pi} {2}}}
Медиана σ 2 ln ⁡ (2) {\ displaystyle \ sigma {\ sqrt {2 \ ln (2)}}}{\ displaystyle \ sigma {\ sqrt {2 \ ln (2)}}}
Режим σ {\ displaystyle \ sigma}\ sigma
Дисперсия 4 - π 2 σ 2 {\ displaystyle {\ frac {4 -\число Пи } {2}} \ sigma ^ {2}}{\ frac {4- \ pi} {2}} \ сигма ^ {2}
Асимметрия 2 π (π - 3) (4 - π) 3/2 {\ displaystyle {\ frac {2 {\ sqrt {\ pi}} (\ pi -3)} {(4- \ pi) ^ {3/2}}}}{\ frac {2 {\ sqrt {\ pi}} (\ pi -3)} {(4- \ pi) ^ {3/2 }}}
Пр. эксцесс - 6 π 2 - 24 π + 16 (4 - π) 2 {\ displaystyle - {\ frac {6 \ pi ^ {2} -24 \ pi +16} {(4- \ pi) ^ { 2}}}}- {\ frac {6 \ pi ^ {2} -24 \ pi +16} {(4- \ pi) ^ {2}}}
Энтропия 1 + ln ⁡ (σ 2) + γ 2 {\ displaystyle 1+ \ ln \ left ({\ frac {\ sigma} {\ sqrt {2}}} \ right) + {\ гидроразрыва {\ gamma} {2}}}1+ \ ln \ left ({\ frac {\ sigma} {\ sqrt {2 }}} \ right) + {\ frac {\ gamma} {2}}
MGF 1 + σ te σ 2 t 2/2 π 2 (erf ⁡ (σ t 2) + 1) {\ displaystyle 1+ \ sigma te ^ {\ sigma ^ {2} t ^ {2} / 2} {\ sqrt {\ frac {\ pi} {2}}} \ left (\ operatorname {erf} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) +1 \ right)}{\ displaystyle 1+ \ sigma te ^ {\ sigma ^ {2} t ^ {2} / 2} {\ sqrt {\ frac {\ pi} {2}}} \ left (\ operatorname {erf} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) +1 \ right)}
CF 1 - σ te - σ 2 t 2/2 π 2 (erfi ⁡ (σ t 2) - i) {\ displaystyle 1 - \ sigma te ^ {- \ sigma ^ {2} t ^ {2} / 2} {\ sqrt {\ frac {\ pi} {2}}} \ left (\ operatorname {erfi} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) -i \ right)}{\ displaystyle 1- \ sigma te ^ {- \ sigma ^ {2} t ^ { 2} / 2} {\ sqrt {\ frac {\ pi} {2}}} \ left (\ operatorname {erfi} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) -i \ right)}

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

. Распределение Рэлея часто наблюдается, когда общая величина вектора связана с его направленными компонентами. Одним из примеров, в котором естественно возникает распределение Рэлея, является анализ скорости ветра в двух измерениях. Если предположить, что каждый компонент некоррелирован, нормально распределен с равной дисперсией и нулевым означает, тогда общая скорость ветра (vector magnitude) будет характеризоваться распределением Рэлея. Второй пример распределения возникает в случае случайных комплексных чисел, действительные и мнимые компоненты которых независимо и одинаково распределены по Гауссу с равной дисперсией и нулевым средним. В этом случае абсолютное значение комплексного числа распределено по Рэлею.

Распределение названо в честь лорда Рэлея ().

Содержание
  • 1 Определение
  • 2 Отношение к длине случайного вектора
  • 3 Свойства
    • 3.1 Дифференциальная энтропия
  • 4 Оценка параметров
    • 4.1 Доверительные интервалы
  • 5 Генерация случайных величин
  • 6 Связанные распределения
  • 7 Приложения
  • 8 См. Также
  • 9 Ссылки
Определение

функция плотности вероятности распределения Рэлея:

f (x; σ) = x σ 2 e - x 2 / (2 σ 2), x ≥ 0, { \ displaystyle f (x; \ sigma) = {\ frac {x} {\ sigma ^ {2}}} e ^ {- x ^ {2} / (2 \ sigma ^ {2})}, \ quad x \ geq 0,}f (x; \ sigma) = {\ frac {x} {\ sigma ^ {2}}} e ^ {-x ^ {2} / (2 \ sigma ^ {2})}, \ quad x \ geq 0,

где σ {\ displaystyle \ sigma}\ sigma - это параметр масштаба распределения. Кумулятивная функция распределения is

F ( Икс; σ) знак равно 1 - е - Икс 2 / (2 σ 2) {\ Displaystyle F (x; \ sigma) = 1-e ^ {- x ^ {2} / (2 \ sigma ^ {2})} }F (x; \ sigma) = 1-e ^ {- х ^ {2} / (2 \ sigma ^ {2})}

для x ∈ [0, ∞). {\ displaystyle x \ in [0, \ infty).}x \ in [0, \ infty).

Отношение к случайной длине вектора

Рассмотрим двумерный вектор Y = (U, V) {\ displaystyle Y = ( U, V)}Y = (U, V) который имеет компоненты, которые нормально распределены, центрированы в нуле и независимы. Тогда U {\ displaystyle U}U и V {\ displaystyle V}V имеют функции плотности

f U (x; σ) = f V (x ; σ) = e - x 2 / (2 σ 2) 2 π σ 2. {\ displaystyle f_ {U} (x; \ sigma) = f_ {V} (x; \ sigma) = {\ frac {e ^ {- x ^ {2} / (2 \ sigma ^ {2})}} {\ sqrt {2 \ pi \ sigma ^ {2}}}}.}{\ displaystyle f_ { U} (x; \ sigma) = f_ {V} (x; \ sigma) = {\ frac {e ^ {- x ^ {2} / (2 \ sigma ^ {2})}} {\ sqrt {2 \ pi \ sigma ^ {2}}}}.}

Пусть X {\ displaystyle X}Х будет длиной Y {\ displaystyle Y}Y . То есть X = U 2 + V 2. {\ displaystyle X = {\ sqrt {U ^ {2} + V ^ {2}}}.}{\ displaystyle X = {\ sqrt {U ^ {2} + V ^ {2}}}.} Тогда X {\ displaystyle X}Х имеет кумулятивную функцию распределения

FX (x; σ) знак равно ∬ D xf U (u; σ) е V (v; σ) d A, {\ displaystyle F_ {X} (x; \ sigma) = \ iint _ {D_ {x }} f_ {U} (u; \ sigma) f_ {V} (v; \ sigma) \, dA,}{\ displaystyle F_ { X} (x; \ sigma) = \ iint _ {D_ {x}} f_ {U} (u; \ sigma) f_ {V} (v; \ sigma) \, dA,}

где D x {\ displaystyle D_ {x}}D_ {x} является диском

D x = {(u, v): u 2 + v 2 < x }. {\displaystyle D_{x}=\left\{(u,v):{\sqrt {u^{2}+v^{2}}}{\ displaystyle D_ {x} = \ left \ {(u, v): {\ sqrt {u ^ {2} + v ^ {2}}} <x \ right \}.}

Записывая двойной интеграл в полярных координатах, он становится

FX (x; σ) = 1 2 π σ 2 ∫ 0 2 π ∫ 0 xre - r 2 / (2 σ 2) drd θ = 1 σ 2 ∫ 0 xre - r 2 / (2 σ 2) dr. {\ displaystyle F_ {X} (x; \ sigma) = {\ frac {1} {2 \ pi \ sigma ^ {2}}} \ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {x} re ^ {- r ^ {2} / (2 \ sigma ^ {2})} \, dr \, d \ theta = {\ frac {1} {\ sigma ^ {2}}} \ int _ {0} ^ {x} re ^ {- r ^ {2} / (2 \ sigma ^ {2})} \, dr.}{\ displaystyle F_ {X} (x; \ sigma) = {\ frac {1} {2 \ pi \ sigma ^ {2}}} \ int _ {0} ^ {2 \ pi} \ int _ {0} ^ {x} re ^ {- r ^ {2} / (2 \ sigma ^ {2})} \, dr \, d \ theta = {\ frac {1} {\ sigma ^ {2}}} \ int _ {0} ^ {x} re ^ {- r ^ {2} / (2 \ sigma ^ {2})} \, dr.}

Наконец, функция плотности вероятности для X {\ displaystyle X}Х - производная его кумулятивной функции распределения, которая согласно фундаментальной теореме исчисления равна

f X (x; σ) = ddx FX (x; σ) = Икс σ 2 е - Икс 2 / (2 σ 2), {\ displaystyle f_ {X} (x; \ sigma) = {\ frac {d} {dx}} F_ {X} (x; \ sigma) = { \ frac {x} {\ sigma ^ {2}}} e ^ {- x ^ {2} / (2 \ sigma ^ {2})},}{\ displaystyle f_ {X} (x; \ sigma) = {\ frac {d} {dx}} F_ {X} (x; \ sigma) = {\ frac {x} {\ sigma ^ {2}}} e ^ {- x ^ {2} / (2 \ sigma ^ {2})},}

которое является распределением Рэлея. Легко обобщить на векторы размерности, отличной от 2. Существуют также обобщения, когда компоненты имеют неравную дисперсию или корреляции, или когда вектор Y следует двумерному t-распределению Стьюдента.

Свойства

исходные моменты задаются следующим образом:

μ j = σ j 2 j / 2 Γ (1 + j 2), {\ displaystyle \ mu _ {j} = \ sigma ^ {j} 2 ^ {j / 2} \, \ Gamma \ left (1 + {\ frac {j} {2}} \ right),}{\ displaystyle \ mu _ {j} = \ sigma ^ {j} 2 ^ {j / 2} \, \ Gamma \ left (1 + {\ frac {j} {2}} \ right),}

где Γ (z) {\ displaystyle \ Gamma (z)}\ Gamma (z) - это гамма-функция.

Таким образом, среднее случайной величины Рэлея:

μ (X) = σ π 2 ≈ 1,253 σ. {\ displaystyle \ mu (X) = \ sigma {\ sqrt {\ frac {\ pi} {2}}} \ \ приблизительно 1,253 \ \ sigma.}{\ displaystyle \ mu (X) = \ sigma {\ sqrt {\ frac {\ pi} {2}}} \ приблизительно 1,253 \ \ sigma.}

дисперсия случайного числа Рэлея переменная:

var ⁡ (X) = μ 2 - μ 1 2 = (2 - π 2) σ 2 ≈ 0,429 σ 2 {\ displaystyle \ operatorname {var} (X) = \ mu _ {2} - \ mu _ {1} ^ {2} = \ left (2 - {\ frac {\ pi} {2}} \ right) \ sigma ^ {2} \ приблизительно 0,429 \ \ sigma ^ {2}}{\ displaystyle \ operatorname {var} (X) = \ mu _ {2} - \ mu _ {1} ^ {2} = \ left (2 - {\ frac {\ pi} {2 }} \ right) \ sigma ^ {2} \ приблизительно 0,429 \ \ sigma ^ {2}}

Режим равен σ, {\ displaystyle \ sigma,}\ sigma, , а максимальный PDF равен

f max = f (σ; σ) = 1 σ e - 1 / 2 ≈ 0.606 σ. {\ displaystyle f _ {\ max} = f (\ sigma; \ sigma) = {\ frac {1} {\ sigma}} e ^ {- 1/2} \ приблизительно {\ frac {0.606} {\ sigma}}.}{\ displaystyle f _ {\ max } = f (\ sigma; \ sigma) = {\ frac {1} {\ sigma}} e ^ {- 1/2} \ приблизительно {\ frac {0.606} {\ sigma}}.}

асимметрия определяется как:

γ 1 = 2 π (π - 3) (4 - π) 3/2 ≈ 0,631 {\ displaystyle \ gamma _ {1} = {\ frac {2 {\ sqrt {\ pi}} (\ pi -3)} {(4- \ pi) ^ {3/2}}} \ приблизительно 0,631}{\ displaystyle \ gamma _ {1} = {\ frac {2 {\ sqrt {\ pi}} (\ pi -3)} {(4- \ pi) ^ {3/2}}} \ приблизительно 0,631}

Избыток эксцесс задается следующим образом:

γ 2 = - 6 π 2 - 24 π + 16 (4 - π) 2 ≈ 0,245 {\ displaystyle \ gamma _ {2} = - {\ frac {6 \ pi ^ {2} -24 \ pi +16} {(4- \ pi) ^ {2}}} \ приблизительно 0,245}\ gamma _ {2} = - {\ frac {6 \ pi ^ {2} -24 \ pi +16} {(4- \ pi) ^ {2}}} \ приблизительно 0,245

Характеристическая функция определяется как:

φ (t) = 1 - σ te - 1 2 σ 2 T 2 π 2 [erfi ⁡ (σ t 2) - я] {\ displaystyle \ varphi (t) = 1- \ sigma te ^ {- {\ frac {1} {2}} \ sigma ^ {2} t ^ {2}} {\ sqrt {\ frac {\ pi} {2}}} \ left [\ operatorname {erfi} \ left ({\ frac {\ sigma t} {\ sqrt {2 }}} \ right) -i \ right]}{\ displaystyle \ varphi (t) = 1- \ sigma te ^ {- {\ frac {1} {2 }} \ sigma ^ {2} t ^ {2}} {\ sqrt {\ frac {\ pi} {2}}} \ left [\ operatorname {erfi} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) -i \ right]}

где erfi ⁡ (z) {\ displaystyle \ operatorname {erfi} (z)}\ operatorname {erfi} (z) - мнимая функция ошибок . порождающая функция момента задается как

M (t) = 1 + σ te 1 2 σ 2 t 2 π 2 [erf ⁡ (σ t 2) + 1] {\ displaystyle M (t) = 1 + \ sigma t \, e ^ {{\ frac {1} {2}} \ sigma ^ {2} t ^ {2}} {\ sqrt {\ frac {\ pi} {2}}} \ left [\ operatorname {erf} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) +1 \ right]}{\ Displaystyle M (t) = 1 + \ sigma t \, e ^ {{\ frac {1} {2}} \ sigma ^ {2} t ^ {2}} {\ sqrt {\ frac { \ pi} {2}}} \ left [\ operatorname {erf} \ left ({\ frac {\ sigma t} {\ sqrt {2}}} \ right) +1 \ right]}

где erf ⁡ (z) {\ displaystyle \ operatorname {erf} (z)}\ operatorname {erf} (z) - это функция ошибок.

Дифференциальная энтропия

дифференциальная энтропия определяется как

H = 1 + пер ⁡ (σ 2) + γ 2 {\ displaystyle H = 1 + \ ln \ left ({\ frac {\ sigma} {\ sqrt {2}}} \ right) + {\ frac {\ gamma} { 2}}}{\ displaystyle H = 1 + \ ln \ left ({\ frac {\ sigma} {\ sqrt {2 }}} \ right) + {\ frac {\ gamma} {2}}}

где γ {\ displaystyle \ gamma}\ гамма - это константа Эйлера – Маскерони.

Оценка параметра

Для выборки N независимые и одинаково распределенные случайные величины Рэлея xi {\ displaystyle x_ {i}}x_ {i} с параметром σ {\ displaystyle \ sigma}\ sigma ,

σ ^ 2 ≈ 1 2 N ∑ я знак равно 1 N xi 2 {\ displaystyle {\ widehat {\ sigma}} ^ {2} \ приблизительно \! \, {\ Frac {1} {2N}} \ s um _ {i = 1} ^ {N} x_ {i} ^ {2}}{\ Displaystyle {\ widehat {\ sig ma}} ^ {2} \ приблизительно \! \, {\ frac {1} {2N}} \ sum _ {i = 1} ^ {N} x_ {i} ^ {2}} - это оценка максимального правдоподобия, а также несмещенная.
σ ^ ≈ 1 2 N ∑ я знак равно 1 N xi 2 {\ displaystyle {\ widehat {\ sigma}} \ приблизительно {\ sqrt {{\ frac {1} {2N}} \ sum _ {i = 1} ^ {N} x_ {i} ^ {2}}}}{\ displaystyle {\ widehat {\ sigma}} \ приблизительно {\ sqrt {{\ frac {1} {2N}} \ sum _ {i = 1} ^ {N} x_ {i} ^ {2}}}} - это смещенная оценка, которую можно исправить с помощью формулы
σ = σ ^ Γ (N) N Γ (N + 1 2) = σ ^ 4 NN! (N - 1)! N (2 N)! π {\ displaystyle \ sigma = {\ widehat {\ sigma}} {\ frac {\ Gamma (N) {\ sqrt {N}}} {\ Gamma (N + {\ frac {1} {2}})}} = {\ widehat {\ sigma}} {\ frac {4 ^ {N} N! (N-1)! {\ sqrt {N}}} {(2N)! {\ sqrt {\ pi}}}}}}{\ displaystyle \ sigma = {\ widehat {\ sigma}} {\ frac {\ Gamma (N) {\ sqrt {N}}} {\ Gamma (N + {\ frac {1} {2}})}} = {\ widehat {\ sigma}} {\ frac {4 ^ {N} N! (N-1)! {\ Sqrt {N} }} {(2N)! {\ Sqrt {\ pi}}}}}

Доверительные интервалы

Чтобы найти доверительный интервал (1 - α), сначала найдите границы [a, b] {\ displaystyle [a, b]}[a, b] где :

п (χ 2 N 2 ≤ a) знак равно α / 2, P (χ 2 N 2 ≤ b) = 1 - α / 2 {\ displaystyle P (\ chi _ {2N} ^ {2} \ leq a) = \ alpha / 2, \ quad P (\ chi _ {2N} ^ {2} \ leq b) = 1- \ alpha / 2}{\ displaystyle P (\ chi _ {2N} ^ {2} \ leq a) = \ alpha / 2, \ quad P (\ chi _ {2N} ^ {2} \ leq b) = 1- \ alpha / 2}

, то параметр масштаба будет находиться в пределах

N x 2 ¯ b ≤ σ ^ 2 ≤ N x 2 ¯ a {\ displaystyle {\ frac {{N} {\ overline {x ^ {2}}}}} {b}} \ leq {\ widehat {\ sigma}} ^ {2} \ leq {\ frac {{N} {\ overline {x ^ {2}}}} {a}}}{\ displaystyle {\ frac {{N} {\ overline { x ^ {2}}}} {b}} \ leq {\ widehat {\ sigma}} ^ {2} \ leq {\ frac {{N} {\ overline {x ^ {2}}}} {a} }}
Генерация случайных величин

Дана случайная величина U, полученная из равномерное распределение в интервале (0, 1), тогда переменная

X = σ - 2 ln ⁡ U {\ displaystyle X = \ sigma {\ sqrt {-2 \ ln U}} \,}{\ displaystyle X = \ sigma { \ sqrt {-2 \ ln U}} \,}

имеет распределение Рэлея с параметром σ {\ displaystyle \ sigma }\ sigma . Это достигается применением метода выборки обратного преобразования.

Связанные распределения
  • R ∼ R ayleigh (σ) {\ displaystyle R \ sim \ mathrm {Rayleigh} (\ sigma)}R \ sim \ mathrm {Rayleigh} (\ sigma) является распределением Рэлея, если R = X 2 + Y 2 {\ displaystyle R = {\ sqrt {X ^ {2} + Y ^ {2}}}}R = {\ sqrt {X ^ {2} + Y ^ {2}}} , где X ∼ N (0, σ 2) {\ displaystyle Икс \ sim N (0, \ sigma ^ {2})}X \ sim N (0, \ sigma ^ {2}) и Y ∼ N (0, σ 2) {\ displaystyle Y \ sim N (0, \ sigma ^ {2})}Y \ sim N (0, \ sigma ^ {2}) - независимые нормальные случайные величины. (Это дает мотивацию к использованию символа «сигма» в приведенной выше параметризации плотности Рэлея.)
  • Если R ∼ R ayleigh (1) {\ displaystyle R \ sim \ mathrm {Rayleigh} (1)}R \ sim \ mathrm {Rayleigh} (1) , то R 2 {\ displaystyle R ^ {2}}R ^ {2} имеет распределение хи-квадрат с параметром N {\ displaystyle N}N , степеней свободы, равное двум (N = 2)
[Q = R 2] ∼ χ 2 (N). {\ Displaystyle [Q = R ^ {2}] \ sim \ chi ^ {2} (N) \.}[Q = R ^ {2}] \ sim \ chi ^ {2} (N) \.
  • Если R ∼ R ayleigh (σ) {\ displaystyle R \ sim \ mathrm {Rayleigh } (\ sigma)}R \ sim \ mathrm {Rayleigh} (\ sigma) , тогда ∑ i = 1 NR i 2 {\ displaystyle \ sum _ {i = 1} ^ {N} R_ {i} ^ {2}}\ sum _ {i = 1} ^ {N} R_ {i} ^ {2} имеет гамма-распределение с параметрами N {\ displaystyle N}N и 2 σ 2 {\ displaystyle 2 \ sigma ^ {2}}2 \ sigma ^ {2}
[Y = ∑ i = 1 NR i 2] ∼ Γ (N, 2 σ 2). {\ displaystyle \ left [Y = \ sum _ {i = 1} ^ {N} R_ {i} ^ {2} \ right] \ sim \ Gamma (N, 2 \ sigma ^ {2}).}\ left [Y = \ sum _ {i = 1} ^ {N} R_ {i} ^ {2} \ right] \ sim \ Gamma (N, 2 \ sigma ^ {2}).
Приложения

Применение оценки σ можно найти в магнитно-резонансной томографии (МРТ). Поскольку изображения МРТ записываются как сложные изображения, но чаще всего рассматриваются как изображения величин, фоновые данные распределены по Рэлею. Следовательно, приведенная выше формула может использоваться для оценки дисперсии шума на МРТ-изображении на основе фоновых данных.

Распределение Рэлея также использовалось в области питания для связи диетического уровни питательных веществ и ответы человека и животного. Таким образом, параметр σ может быть использован для расчета зависимости отклика питательных веществ.

См. Также
Ссылки
Последняя правка сделана 2021-06-03 09:25:24
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте