Cycle space: Difference between revisions
en>Beroal m →The binary cycle space: free vector space |
→The binary cycle space: Changed wikilink |
||
Line 1: | Line 1: | ||
{{Distinguish|Rayleigh mixture distribution}} | |||
{{Probability distribution| | |||
name =Rayleigh| | |||
type =density| | |||
pdf_image =[[Image:Rayleigh distributionPDF.svg|325px|Plot of the Rayleigh PDF]]<br /><small></small>| | |||
cdf_image =[[Image:Rayleigh distributionCDF.svg|325px|Plot of the Rayleigh CDF]]<br /><small></small>| | |||
parameters =<math>\sigma>0\,</math>| | |||
support =<math>x\in [0,+\infty)</math>| | |||
pdf =<math>\frac{x}{\sigma^2} e^{-x^2/2\sigma^2}</math>| | |||
cdf =<math>1 - e^{-x^2/2\sigma^2}</math>| | |||
mean =<math>\sigma \sqrt{\frac{\pi}{2}}</math>| | |||
median =<math>\sigma\sqrt{\ln(4)}\,</math>| | |||
mode =<math>\sigma\,</math>| | |||
variance =<math>\frac{4 - \pi}{2} \sigma^2</math>| | |||
skewness =<math>\frac{2\sqrt{\pi}(\pi - 3)}{(4-\pi)^{3/2}}</math>| | |||
kurtosis =<math>-\frac{6\pi^2 - 24\pi +16}{(4-\pi)^2}</math>| | |||
entropy =<math>1+\ln\left(\frac{\sigma}{\sqrt{2}}\right)+\frac{\gamma}{2}</math>| | |||
mgf =<math>1+\sigma t\,e^{\sigma^2t^2/2}\sqrt{\frac{\pi}{2}} | |||
\left(\textrm{erf}\left(\frac{\sigma t}{\sqrt{2}}\right)\!+\!1\right)</math>| | |||
char =<math>1\!-\!\sigma te^{-\sigma^2t^2/2}\sqrt{\frac{\pi}{2}}\!\left(\textrm{erfi}\!\left(\frac{\sigma t}{\sqrt{2}}\right)\!-\!i\right)</math>| | |||
}} | |||
In [[probability theory]] and [[statistics]], the '''Rayleigh distribution''' {{IPAc-en|ˈ|r|eɪ|l|i}} is a [[continuous probability distribution]] for positive-valued [[random variable]]s. | |||
A Rayleigh distribution is often observed when the overall magnitude of a vector is related to its directional [[Euclidean_vector#Vector_components|components]]. One example where the Rayleigh distribution naturally arises is when wind velocity is analyzed into its orthogonal 2-dimensional vector components. Assuming that the magnitudes of each component are [[uncorrelated]], [[Normal distribution|normally distributed]] with equal [[variance]], and zero [[mean]], then the overall wind speed ([[Euclidean vector|vector]] magnitude) will be characterized by a Rayleigh distribution. A second example of the distribution arises in the case of random complex numbers whose real and imaginary components are i.i.d. (independently and identically distributed) [[normal distribution|Gaussian]] with equal variance and zero mean. In that case, the absolute value of the complex number is Rayleigh-distributed. | |||
The distribution is named after [[John Strutt, 3rd Baron Rayleigh|Lord Rayleigh]].{{Citation needed|date=April 2013}} | |||
==Definition== | |||
The [[probability density function]] of the Rayleigh distribution is<ref Name=PP>Papoulis, Athanasios; Pillai, S. (2001) ''Probability, Random Variables and Stochastic Processe''. ISBN 0073660116, ISBN 9780073660110 {{Page needed|date=April 2013}}</ref> | |||
:<math>f(x;\sigma) = \frac{x}{\sigma^2} e^{-x^2/2\sigma^2}, \quad x \geq 0,</math> | |||
where <math>\sigma >0,</math> is the [[scale parameter]] of the distribution. The [[cumulative distribution function]] is<ref Name=PP/> | |||
:<math>F(x) = 1 - e^{-x^2/2\sigma^2}</math> | |||
for <math>x \in [0,\infty).</math> | |||
==Relation to random vector lengths== | |||
Consider the two-dimensional vector <math> Y = (U,V) </math> which has components that are Gaussian-distributed and independent. Then <math> f_U(u; \sigma) = \frac{e^{-u^2/2\sigma^2}}{\sqrt{2\pi\sigma^2}} </math>, and similarly for <math> f_V(v; \sigma) </math>. | |||
Let <math> x </math> be the length of <math> Y </math>. It is distributed as | |||
: <math>f(x; \sigma) = \frac{1}{2\pi\sigma^2} \int_{-\infty}^\infty du \, \int_{-\infty}^\infty dv \, e^{-u^2/2\sigma^2} e^{-v^2/2\sigma^2} \delta(x-\sqrt{u^2+v^2}).</math> | |||
By transforming to the [[polar coordinate system]] one has | |||
: <math> f(x; \sigma) = \frac{1}{2\pi\sigma^2} \int_0^{2\pi} \, d\phi \int_0^\infty dr \, \delta(r-x) r e^{-r^2/2\sigma^2}= \frac{x}{\sigma^2} e^{-x^2/2\sigma^2}, | |||
</math> | |||
which is the Rayleigh distribution. It is straightforward to generalize to vectors of dimension other than 2. | |||
There are also generalizations when the components have unequal variance or correlations. | |||
==Properties== | |||
The raw [[moment (mathematics)|moments]] are given by: | |||
:<math>\mu_k = \sigma^k2^\frac{k}{2}\,\Gamma\left(1 + \frac{k}{2}\right)</math> | |||
where <math>\Gamma(z)</math> is the [[Gamma function]]. | |||
The [[mean]] and [[variance]] of a Rayleigh [[random variable]] may be expressed as: | |||
:<math>\mu(X) = \sigma \sqrt{\frac{\pi}{2}}\ \approx 1.253 \sigma</math> | |||
and | |||
:<math>\textrm{var}(X) = \frac{4 - \pi}{2} \sigma^2 \approx 0.429 \sigma^2</math> | |||
The mode is <math>\sigma </math> and the maximum pdf is | |||
:<math> f_\text{max} = f(\sigma;\sigma) = \frac{1}{\sigma} e^{-\frac{1}{2}} \approx \frac{1}{\sigma} 0.606</math> | |||
The [[skewness]] is given by: | |||
:<math>\gamma_1 = \frac{2\sqrt{\pi}(\pi - 3)}{(4 - \pi)^\frac{3}{2}} \approx 0.631</math> | |||
The excess [[kurtosis]] is given by: | |||
:<math>\gamma_2 = -\frac{6\pi^2 - 24\pi + 16}{(4 - \pi)^2} \approx 0.245</math> | |||
The [[characteristic function (probability theory)|characteristic function]] is given by: | |||
:<math>\varphi(t) = 1 - \sigma te^{-\frac{1}{2}\sigma^2t^2}\sqrt{\frac{\pi}{2}} \left[\textrm{erfi} \left(\frac{\sigma t}{\sqrt{2}}\right) - i\right]</math> | |||
where <math>\operatorname{erfi}(z)</math> is the imaginary [[error function]]. The [[moment generating function]] is given by | |||
:<math> | |||
M(t) = 1 + \sigma t\,e^{\frac{1}{2}\sigma^2t^2}\sqrt{\frac{\pi}{2}} | |||
\left[\textrm{erf}\left(\frac{\sigma t}{\sqrt{2}}\right) + 1\right]</math> | |||
where <math>\operatorname{erf}(z)</math> is the [[error function]]. | |||
===Differential entropy=== | |||
The [[differential entropy]] is given by{{Citation needed|date=April 2013}} | |||
:<math>H = 1 + \ln\left(\frac{\sigma}{\sqrt{2}}\right) + \frac{\gamma}{2}</math> | |||
where <math>\gamma</math> is the [[Euler–Mascheroni constant]]. | |||
== Parameter estimation == | |||
Given a sample of ''N'' [[independent and identically distributed]] Rayleigh random variables <math>x_i</math> with parameter <math>\sigma</math>, | |||
:<math>\widehat{\sigma^2}\approx \!\,\frac{1}{2N}\sum_{i=1}^N x_i^2</math> is an unbiased [[maximum likelihood]] estimate. | |||
:<math>\hat{\sigma}\approx \!\,\sqrt{\frac{1}{2N}\sum_{i=1}^N x_i^2}</math> is a biased estimator that can be corrected via the formula | |||
:<math>\sigma = \hat{\sigma} \frac {\Gamma(N)\sqrt{N}} {\Gamma(N + \frac {1} {2})} = \hat{\sigma} \frac {4^N N!(N-1)!\sqrt{N}} {(2N)!\sqrt{\pi}}</math><ref>[https://archive.org/details/jresv68Dn9p1005 Siddiqui, M. M. (1964) "Statistical inference for Rayleigh distributions", ''The Journal of Research of the National Bureau of Standards, Sec. D: Radio Science'', Vol. 68D, No. 9, p. 1007]</ref> | |||
=== Confidence intervals === | |||
To find the (1 − ''α'') confidence interval, first find <math>\chi_1^2, \ \chi_2^2</math> where: | |||
: <math>Pr(\chi^2(2n) \leq \chi_1^2) = \alpha/2, \quad Pr(\chi^2(2n) \leq \chi_2^2) = 1 - \alpha/2</math> | |||
then | |||
: <math>\frac{2n\overline{x^2}}{\chi_2^2} \leq \widehat{\sigma}^2 \leq \frac{2n\overline{x^2}}{\chi_1^2}</math><ref>[http://nvlpubs.nist.gov/nistpubs/jres/66D/jresv66Dn2p167_A1b.pdf Siddiqui, M. M. (1961) "Some Problems Connected With Rayleigh Distributions", ''The Journal of Research of the National Bureau of Standards, Sec. D: Radio Propagation'', Vol. 66D, No. 2, p. 169]</ref> | |||
== Generating random variates == | |||
Given a random variate ''U'' drawn from the [[uniform distribution (continuous)|uniform distribution]] in the interval <nowiki>(0, 1)</nowiki>, then the variate | |||
:<math>X=\sigma\sqrt{-2 \ln(U)}\,</math> | |||
has a Rayleigh distribution with parameter <math>\sigma</math>. This is obtained by applying the [[inverse transform sampling]]-method. | |||
==Related distributions== | |||
*<math>R \sim \mathrm{Rayleigh}(\sigma)</math> is Rayleigh distributed if <math>R = \sqrt{X^2 + Y^2}</math>, where <math>X \sim N(0, \sigma^2)</math> and <math>Y \sim N(0, \sigma^2)</math> are independent [[Normal_distribution|normal random variables]].<ref>[http://home.kpn.nl/jhhogema1966/skeetn/ballist/sgs/sgs.htm#_Toc96439743 Hogema, Jeroen (2005) "Shot group statistics"]</ref> (This gives motivation to the use of the symbol "sigma" in the above parameterization of the Rayleigh density.) | |||
*The [[chi distribution]] with ''v'' = 2 is equivalent to Rayleigh Distribution with ''σ'' = 1. I.e., if <math>R \sim \mathrm{Rayleigh} (1)</math>, then <math>R^2</math> has a [[chi-squared distribution]] with parameter <math>N</math>, degrees of freedom, equal to two (''N'' = 2) | |||
:: <math>[Q=R^2] \sim \chi^2(N)\ .</math> | |||
*If <math>R \sim \mathrm{Rayleigh}(\sigma)</math>, then <math>\sum_{i=1}^N R_i^2</math> has a [[gamma distribution]] with parameters <math>N</math> and <math>\frac{1}{2\sigma^2}</math> | |||
:: <math>\left[Y=\sum_{i=1}^N R_i^2\right] \sim \Gamma(N,\frac{1}{2\sigma^2}) .</math> | |||
*The [[Rice distribution]] is a generalization of the Rayleigh distribution. | |||
*The [[Weibull distribution]] is a generalization of the Rayleigh distribution. In this instance, parameter <math>\sigma</math> is related to the Weibull scale parameter <math>\lambda</math>: <math>\lambda = \sigma \sqrt{2} .</math> | |||
*The [[Maxwell–Boltzmann distribution]] describes the magnitude of a normal vector in three dimensions. | |||
*If <math>X</math> has an [[exponential distribution]] <math>X \sim \mathrm{Exponential}(\lambda)</math>, then <math>Y=\sqrt{2X\sigma^2\lambda} \sim \mathrm{Rayleigh}(\sigma) .</math> | |||
== Applications == | |||
An application of the estimation of σ can be found in [[magnetic resonance imaging]] (MRI). As MRI images are recorded as [[complex numbers|complex]] images but most often viewed as magnitude images, the background data is Rayleigh distributed. Hence, the above formula can be used to estimate the noise variance in an MRI image from background data.<ref>Sijbers J., den Dekker A. J., Raman E. and Van Dyck D. (1999) "Parameter estimation from magnitude MR images", ''International Journal of Imaging Systems and Technology'', 10(2), 109–114</ref> | |||
==See also== | |||
*[[Rayleigh fading]] | |||
*[[Rayleigh mixture distribution]] | |||
{{More footnotes|date=April 2013}} | |||
== References == | |||
{{reflist}} | |||
{{ProbDistributions|continuous-semi-infinite}} | |||
{{DEFAULTSORT:Rayleigh Distribution}} | |||
[[Category:Continuous distributions]] | |||
[[Category:Exponential family distributions]] | |||
[[Category:Probability distributions]] |
Revision as of 13:20, 20 June 2013
Template:Probability distribution
In probability theory and statistics, the Rayleigh distribution Template:IPAc-en is a continuous probability distribution for positive-valued random variables.
A Rayleigh distribution is often observed when the overall magnitude of a vector is related to its directional components. One example where the Rayleigh distribution naturally arises is when wind velocity is analyzed into its orthogonal 2-dimensional vector components. Assuming that the magnitudes of each component are uncorrelated, normally distributed with equal variance, and zero mean, then the overall wind speed (vector magnitude) will be characterized by a Rayleigh distribution. A second example of the distribution arises in the case of random complex numbers whose real and imaginary components are i.i.d. (independently and identically distributed) Gaussian with equal variance and zero mean. In that case, the absolute value of the complex number is Rayleigh-distributed.
The distribution is named after Lord Rayleigh.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
Definition
The probability density function of the Rayleigh distribution is[1]
where is the scale parameter of the distribution. The cumulative distribution function is[1]
Relation to random vector lengths
Consider the two-dimensional vector which has components that are Gaussian-distributed and independent. Then , and similarly for .
Let be the length of . It is distributed as
By transforming to the polar coordinate system one has
which is the Rayleigh distribution. It is straightforward to generalize to vectors of dimension other than 2. There are also generalizations when the components have unequal variance or correlations.
Properties
The raw moments are given by:
where is the Gamma function.
The mean and variance of a Rayleigh random variable may be expressed as:
and
The mode is and the maximum pdf is
The skewness is given by:
The excess kurtosis is given by:
The characteristic function is given by:
where is the imaginary error function. The moment generating function is given by
where is the error function.
Differential entropy
The differential entropy is given byPotter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
where is the Euler–Mascheroni constant.
Parameter estimation
Given a sample of N independent and identically distributed Rayleigh random variables with parameter ,
- is an unbiased maximum likelihood estimate.
Confidence intervals
To find the (1 − α) confidence interval, first find where:
then
Generating random variates
Given a random variate U drawn from the uniform distribution in the interval (0, 1), then the variate
has a Rayleigh distribution with parameter . This is obtained by applying the inverse transform sampling-method.
Related distributions
- is Rayleigh distributed if , where and are independent normal random variables.[4] (This gives motivation to the use of the symbol "sigma" in the above parameterization of the Rayleigh density.)
- The chi distribution with v = 2 is equivalent to Rayleigh Distribution with σ = 1. I.e., if , then has a chi-squared distribution with parameter , degrees of freedom, equal to two (N = 2)
- If , then has a gamma distribution with parameters and
- The Rice distribution is a generalization of the Rayleigh distribution.
- The Weibull distribution is a generalization of the Rayleigh distribution. In this instance, parameter is related to the Weibull scale parameter :
- The Maxwell–Boltzmann distribution describes the magnitude of a normal vector in three dimensions.
- If has an exponential distribution , then
Applications
An application of the estimation of σ can be found in magnetic resonance imaging (MRI). As MRI images are recorded as complex images but most often viewed as magnitude images, the background data is Rayleigh distributed. Hence, the above formula can be used to estimate the noise variance in an MRI image from background data.[5]
See also
References
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
55 yrs old Metal Polisher Records from Gypsumville, has interests which include owning an antique car, summoners war hack and spelunkering. Gets immense motivation from life by going to places such as Villa Adriana (Tivoli).
my web site - summoners war hack no survey ios
- ↑ 1.0 1.1 Papoulis, Athanasios; Pillai, S. (2001) Probability, Random Variables and Stochastic Processe. ISBN 0073660116, ISBN 9780073660110 Template:Page needed
- ↑ Siddiqui, M. M. (1964) "Statistical inference for Rayleigh distributions", The Journal of Research of the National Bureau of Standards, Sec. D: Radio Science, Vol. 68D, No. 9, p. 1007
- ↑ Siddiqui, M. M. (1961) "Some Problems Connected With Rayleigh Distributions", The Journal of Research of the National Bureau of Standards, Sec. D: Radio Propagation, Vol. 66D, No. 2, p. 169
- ↑ Hogema, Jeroen (2005) "Shot group statistics"
- ↑ Sijbers J., den Dekker A. J., Raman E. and Van Dyck D. (1999) "Parameter estimation from magnitude MR images", International Journal of Imaging Systems and Technology, 10(2), 109–114