Calcular el área bajo la campana de Gauss
Introducción
La función definida de la siguiente forma:
se denomina
función gaussiana y su gráfica tiene forma de campana. Tomando ciertos valores de
y
obtenemos que esta función es la función de densidad de una variable aleatoria normal:
Si
es una variable aleatoria que sigue una distribución normal de media
y desviación típica
(suele escribirse también de la forma
), entonces
tiene como función de densidad a:
Para comprobar que
es una función de densidad debemos comprobar estas dos condiciones:
1.-
2.-
La primera condición es evidente, al ser
(por definición de desviación típica) y por ser la exponencial siempre positiva. La comprobación de la segunda condición consiste simplemente en el cálculo de esa integral impropia…
¿Simplemente?Quienes no conozcan esta integral pueden intentar calcularla encontrando una primitiva y evaluando después en
y en
. Tiempo perdido. No podemos encontrar dicha primitiva. Pero, por suerte, sí podemos calcular esa integral, aunque no con cálculo integral de una variable sino con
cálculo integral de dos variables.
Explicación previa
Para empezar es evidente que
es simétrica respecto de la recta
(sustituid
por
y
y veréis que se obtiene el mismo resultado). Por simplicidad para los cálculos tomaremos
y
, es decir, la función de densidad de la conocida
. Esto es, calcularemos la integral siguiente:
Con ello, por tanto, obtenemos una
función par (simétrica respecto del eje Y), por lo que bastará con calcular la integral entre
e
y después multiplicar por dos el resultado.
Hemos dicho que vamos a utilizar cálculo integral de dos variables. Concretamente calcularemos una integral doble (en dos variables) y deduciremos de ese cálculo cuánto vale la integral buscada. Para el cálculo de dicha integral doble serán esenciales las
coordenadas polares.
Cálculo de la integral
La integral doble que vamos a calcular es la siguiente:
Para ello realizamos un cambio a coordenadas polares:
Al realizar este cambio los intervalos en los que toman valores
y
son los siguientes:
Por otra parte se tiene que en polares
, por lo que el exponente de la exponencial quedará
. Y por realizar un cambio de variable en una integral doble debemos multiplicar por el valor absoluto del jacobiano del cambio de variable, que en este caso es
.
Ya estamos preparados para plantear la integral:
Al quedarnos una función independiente de
obtenemos:
Calculamos dicha integral:
![=\frac{\pi}{2} \cdot (-1) \cdot (e^{\frac{-r^2}{2}}) \Bigg ]_0^{\infty}=\frac{\pi}{2}](http://s0.wp.com/latex.php?latex=%3D%5Cfrac%7B%5Cpi%7D%7B2%7D+%5Ccdot+%28-1%29+%5Ccdot+%28e%5E%7B%5Cfrac%7B-r%5E2%7D%7B2%7D%7D%29+%5CBigg+%5D_0%5E%7B%5Cinfty%7D%3D%5Cfrac%7B%5Cpi%7D%7B2%7D&bg=ffffff&fg=000000&s=0 "=\frac{\pi}{2} \cdot (-1) \cdot (e^{\frac{-r^2}{2}}) \Bigg ]_0^{\infty}=\frac{\pi}{2}")
Obtenemos entonces la primera parte de nuestro cálculo:
Por otro lado:
Sacamos la función independiente de
(la segunda exponencial) de la integral respecto de
:
![\displaystyle{=\int_0^{\infty} e^{\frac{-y^2}{2}} \left [\int_0^{\infty} e^{\frac{-x^2}{2}} \, dx \right ] dy=}](http://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%3D%5Cint_0%5E%7B%5Cinfty%7D+e%5E%7B%5Cfrac%7B-y%5E2%7D%7B2%7D%7D+%5Cleft+%5B%5Cint_0%5E%7B%5Cinfty%7D+e%5E%7B%5Cfrac%7B-x%5E2%7D%7B2%7D%7D+%5C%2C+dx+%5Cright+%5D+dy%3D%7D&bg=ffffff&fg=000000&s=0 "\displaystyle{=\int_0^{\infty} e^{\frac{-y^2}{2}} \left [\int_0^{\infty} e^{\frac{-x^2}{2}} \, dx \right ] dy=}")
Como la integral que queda dentro es independiente de
puede sacarse completamente de la integral respecto de
, quedando:
Obtenemos entonces un producto de dos integrales iguales simplemente renombrando la
de la segunda integral a
(esto no supone ningún problema, es simplemente un cambio de nombre de la variable). Tenemos entonces la integral al cuadrado. Tomando el principio y el final de este paso llegamos a:
Uniendo los resultados obtenidos de estos dos pasos llegamos a:
Tomando raíces cuadradas se tiene lo siguiente:
Recordemos que la función a integrar era par, por lo que el valor de la integral completa será el resultado de multiplicar lo obtenido por dos, es decir:
Ahora tomamos la definición inicial de la función y llegamos al resultado buscado:
