stringtranslate.com

Función de error

En matemáticas, la función de error (también llamada función de error de Gauss ), a menudo denotada por erf , es una función definida como: [1]

La integral aquí es una integral de contorno compleja que no depende de la trayectoria porque es holomorfa en todo el plano complejo . En muchas aplicaciones, el argumento de la función es un número real, en cuyo caso el valor de la función también es real.

En algunos textos antiguos, [2] la función de error se define sin el factor de . Esta integral no elemental es una función sigmoidea que aparece a menudo en probabilidad , estadística y ecuaciones diferenciales parciales .

En estadística, para valores reales no negativos de x , la función de error tiene la siguiente interpretación: para una variable aleatoria real Y que se distribuye normalmente con media 0 y desviación estándar , erf x es la probabilidad de que Y caiga en el rango [− x , x ] .

Dos funciones estrechamente relacionadas son la función de error complementaria se define como

y la función de error imaginaria se define como

donde i es la unidad imaginaria .

Nombre

El nombre "función de error" y su abreviatura erf fueron propuestos por JWL Glaisher en 1871 debido a su conexión con "la teoría de la probabilidad, y en particular la teoría de los errores ". [3] El complemento de la función de error también fue discutido por Glaisher en una publicación separada en el mismo año. [4] Para la "ley de facilidad" de errores cuya densidad está dada por (la distribución normal ), Glaisher calcula la probabilidad de un error que se encuentra entre p y q como:

Gráfico de la función de error Erf(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1
Gráfico de la función de error Erf(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1

Aplicaciones

Cuando los resultados de una serie de mediciones se describen mediante una distribución normal con desviación estándar σ y valor esperado 0, entonces erf ( a/σ 2 ) ​​es la probabilidad de que el error de una sola medición se encuentre entrea y + a , para a positivo. Esto es útil, por ejemplo, para determinar la tasa de error de bits de un sistema de comunicación digital.

Las funciones de error y error complementaria ocurren, por ejemplo, en soluciones de la ecuación de calor cuando las condiciones de contorno están dadas por la función escalón de Heaviside .

La función de error y sus aproximaciones se pueden utilizar para estimar resultados que se cumplen con alta o baja probabilidad. Dada una variable aleatoria X ~ Norm[ μ , σ ] (una distribución normal con media μ y desviación estándar σ ) y una constante L > μ , se puede demostrar mediante integración por sustitución:

donde A y B son ciertas constantes numéricas. Si L está suficientemente lejos de la media, específicamente μLσ ln k , entonces:

Por lo tanto, la probabilidad tiende a 0 cuando k → ∞ .

La probabilidad de que X esté en el intervalo [ L a , L b ] se puede derivar como

Propiedades

Parcelas en el plano complejo

La propiedad erf (− z ) = −erf z significa que la función error es una función impar . Esto resulta directamente del hecho de que el integrando e t 2 es una función par (la antiderivada de una función par que es cero en el origen es una función impar y viceversa).

Dado que la función de error es una función entera que convierte números reales en números reales, para cualquier número complejo z : donde z es el conjugado complejo de z .

El integrando f = exp(− z 2 ) y f = erf z se muestran en el plano z complejo en las figuras de la derecha con coloración de dominio .

La función de error en +∞ es exactamente 1 (ver integral de Gauss ). En el eje real, erf z tiende a la unidad en z → +∞ y a −1 en z → −∞ . En el eje imaginario, tiende a ± i .

Serie de Taylor

La función de error es una función entera ; no tiene singularidades (excepto la del infinito) y su expansión de Taylor siempre converge. Sin embargo, para x >> 1 , la cancelación de los términos principales hace que la expansión de Taylor no sea práctica.

La integral definitoria no se puede evaluar en forma cerrada en términos de funciones elementales (véase el teorema de Liouville ), pero al expandir el integrando e z 2 en su serie de Maclaurin e integrar término por término, se obtiene la serie de Maclaurin de la función de error como: que se cumple para cada número complejo z . Los términos del denominador son la secuencia A007680 en la OEIS . 

Para el cálculo iterativo de la serie anterior, puede ser útil la siguiente formulación alternativa : porque−(2 k − 1) z 2/k (2k + 1) expresa el multiplicador para convertir el k -ésimo término en el ( k  + 1) -ésimo término (considerando z como el primer término).

La función de error imaginaria tiene una serie de Maclaurin muy similar, que es: que es válida para cada número complejo z . 

Derivada e integral

La derivada de la función de error se sigue inmediatamente de su definición: De esto, la derivada de la función de error imaginaria también es inmediata: Una antiderivada de la función de error, obtenible por integración por partes , es Una antiderivada de la función de error imaginaria, también obtenible por integración por partes, es Las derivadas de orden superior están dadas por donde H son los polinomios de Hermite de los físicos . [5]

Serie Bürmann

Una expansión, [6] que converge más rápidamente para todos los valores reales de x que una expansión de Taylor, se obtiene utilizando el teorema de Hans Heinrich Bürmann : [7] donde sgn es la función de signo . Manteniendo solo los dos primeros coeficientes y eligiendo c 1 = 31/200 y c 2 = − 341/8000 , la aproximación resultante muestra su mayor error relativo en x = ±1,3796 , donde es menor que 0,0036127:

Funciones inversas

Función de error inverso

Dado un número complejo z , no existe un único número complejo w que satisfaga erf w = z , por lo que una función inversa verdadera sería multivaluada. Sin embargo, para −1 < x < 1 , existe un único número real denotado erf −1 x que satisface

La función de error inversa se define habitualmente con el dominio (−1,1) , y está restringida a este dominio en muchos sistemas de álgebra computacional. Sin embargo, se puede extender al disco | z | < 1 del plano complejo, utilizando la serie de Maclaurin [8] donde c 0 = 1 y

Así tenemos la expansión de la serie (los factores comunes han sido cancelados de los numeradores y denominadores): (Después de la cancelación, los valores del numerador y denominador en OEIS : A092676 y OEIS : A092677 respectivamente; sin cancelación, los términos del numerador son valores en OEIS : A002067 ). El valor de la función de error en  ±∞ es igual a  ±1 .

Para | z | < 1 , tenemos erf(erf −1 z ) = z .

La función de error complementaria inversa se define como Para x real , existe un único número real erfi −1 x que satisface erfi(erfi −1 x ) = x . La función de error imaginaria inversa se define como erfi −1 x . [9]

Para cualquier x real , se puede utilizar el método de Newton para calcular erfi −1 x , y para −1 ≤ x ≤ 1 , la siguiente serie de Maclaurin converge: donde c k se define como arriba.

Expansión asintótica

Una expansión asintótica útil de la función de error complementaria (y por lo tanto también de la función de error) para x reales grandes es donde (2 n − 1)!! es el factorial doble de (2 n − 1) , que es el producto de todos los números impares hasta (2 n − 1) . Esta serie diverge para cada x finito , y su significado como expansión asintótica es que para cualquier entero N ≥ 1 se tiene donde el resto es que se sigue fácilmente por inducción, escribiendo e integrando por partes.

El comportamiento asintótico del término restante, en notación de Landau , es cuando x → ∞ . Esto se puede encontrar mediante Para valores suficientemente grandes de x , solo se necesitan los primeros términos de esta expansión asintótica para obtener una buena aproximación de erfc x (mientras que para valores no demasiado grandes de x , la expansión de Taylor anterior en 0 proporciona una convergencia muy rápida).

Expansión de fracciones continuas

Laplace encontró una expansión de fracción continua de la función de error complementaria : [10] [11]

Serie factorial

La serie factorial inversa: converge para Re( z 2 ) > 0 . Aquí z n denota el factorial ascendente y s ( n , k ) denota un número de Stirling con signo de primera especie . [12] [13] También existe una representación mediante una suma infinita que contiene el factorial doble :

Aproximaciones numéricas

Aproximación con funciones elementales

Tabla de valores

Funciones relacionadas

Función de error complementaria

La función de error complementaria , denotada erfc , se define como

Gráfico de la función de error complementaria Erfc(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1
Gráfico de la función de error complementaria Erfc(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1

que también define erfcx , la función de error complementaria escalada [25] (que se puede utilizar en lugar de erfc para evitar el desbordamiento aritmético [25] [26] ). Otra forma de erfc x para x ≥ 0 se conoce como fórmula de Craig, en honor a su descubridor: [27] Esta expresión es válida solo para valores positivos de x , pero se puede utilizar junto con erfc x = 2 − erfc(− x ) para obtener erfc( x ) para valores negativos. Esta forma es ventajosa porque el rango de integración es fijo y finito. Una extensión de esta expresión para el erfc de la suma de dos variables no negativas es la siguiente: [28]

Función de error imaginario

La función de error imaginaria , denotada erfi , se define como

Gráfico de la función de error imaginaria Erfi(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1
Gráfico de la función de error imaginaria Erfi(z) en el plano complejo desde -2-2i hasta 2+2i con colores creado con la función ComplexPlot3D de Mathematica 13.1

donde D ( x ) es la función de Dawson (que puede utilizarse en lugar de erfi para evitar el desbordamiento aritmético [25] ).

A pesar del nombre "función de error imaginaria", erfi x es real cuando x es real.

Cuando se evalúa la función de error para argumentos complejos arbitrarios z , la función de error compleja resultante generalmente se analiza en forma escalada como la función de Faddeeva :

Función de distribución acumulativa

La función de error es esencialmente idéntica a la función de distribución acumulativa normal estándar , denotada Φ , también llamada norm( x ) por algunos lenguajes de software [ cita requerida ] , ya que difieren solo en la escala y la traducción. De hecho,

La función de distribución acumulativa normal representada en el plano complejo.
La función de distribución acumulativa normal representada en el plano complejo.

o reorganizado para erf y erfc :

En consecuencia, la función de error también está estrechamente relacionada con la función Q , que es la probabilidad de cola de la distribución normal estándar. La función Q se puede expresar en términos de la función de error como

La inversa de Φ se conoce como función cuantil normal o función probit y puede expresarse en términos de la función de error inversa como

La función de distribución normal estándar se utiliza con más frecuencia en probabilidad y estadística, y la función de error se utiliza con más frecuencia en otras ramas de las matemáticas.

La función de error es un caso especial de la función Mittag-Leffler , y también puede expresarse como una función hipergeométrica confluente (función de Kummer):

Tiene una expresión simple en términos de la integral de Fresnel . [ se necesita más explicación ]

En términos de la función gamma regularizada P y la función gamma incompleta , sgn x es la función de signo .

Integrales iteradas de la función de error complementaria

Las integrales iteradas de la función de error complementaria están definidas por [29]

La fórmula general de recurrencia es

Tienen la serie de potencias de las cuales se siguen las propiedades de simetría y

Implementaciones

Como función real de un argumento real

Como función compleja de un argumento complejo

Referencias

  1. ^ Andrews, Larry C. (1998). Funciones especiales de las matemáticas para ingenieros. SPIE Press. p. 110. ISBN 9780819426161.
  2. ^ Whittaker, Edmund Taylor ; Watson, George Neville (2021). Moll, Victor Hugo (ed.). Un curso de análisis moderno (quinta edición revisada). Cambridge University Press . p. 358. ISBN 978-1-316-51893-9.
  3. ^ Glaisher, James Whitbread Lee (julio de 1871). "Sobre una clase de integrales definidas". Londres, Edimburgo y Dublín Philosophical Magazine and Journal of Science . 4. 42 (277): 294–302. doi :10.1080/14786447108640568 . Consultado el 6 de diciembre de 2017 .
  4. ^ Glaisher, James Whitbread Lee (septiembre de 1871). "Sobre una clase de integrales definidas. Parte II". Londres, Edimburgo y Dublín Philosophical Magazine and Journal of Science . 4. 42 (279): 421–436. doi :10.1080/14786447108640600 . Consultado el 6 de diciembre de 2017 .
  5. ^ Weisstein, Eric W. "Erf". MundoMatemático .
  6. ^ Schöpf, HM; Supancic, PH (2014). "Sobre el teorema de Bürmann y su aplicación a problemas de transferencia y difusión de calor lineal y no lineal". The Mathematica Journal . 16 . doi : 10.3888/tmj.16-11 .
  7. ^ Weisstein, Eric W. "Teorema de Bürmann". MathWorld .
  8. ^ Dominici, Diego (2006). "Análisis asintótico de las derivadas de la función de error inversa". arXiv : math/0607230 .
  9. ^ Bergsma, Wicher (2006). "Sobre un nuevo coeficiente de correlación, su descomposición ortogonal y pruebas de independencia asociadas". arXiv : math/0604627 .
  10. Pierre-Simon Laplace , Traité de mécanique céleste , tomo 4 (1805), libro X, página 255.
  11. ^ Cuyt, Annie AM ; Petersen, Vigdis B.; Verdonk, Brigitte; Waadeland, Haakon; Jones, William B. (2008). Manual de fracciones continuas para funciones especiales . Springer-Verlag. ISBN 978-1-4020-6948-2.
  12. ^ Schlömilch, Oskar Xavier (1859). "Ueber facultätenreihen". Zeitschrift für Mathematik und Physik (en alemán). 4 : 390–415.
  13. ^ Nielson, Niels (1906). Handbuch der Theorie der Gammafunktion (en alemán). Leipzig: BG Teubner. pag. 283 ecuación. 3 . Consultado el 4 de diciembre de 2017 .
  14. ^ Chiani, M.; Dardari, D.; Simon, MK (2003). "Nuevos límites exponenciales y aproximaciones para el cálculo de la probabilidad de error en canales con desvanecimiento" (PDF) . IEEE Transactions on Wireless Communications . 2 (4): 840–845. CiteSeerX 10.1.1.190.6761 . doi :10.1109/TWC.2003.814350. 
  15. ^ Tanash, IM; Riihonen, T. (2020). "Aproximaciones y límites minimax globales para la función Q gaussiana mediante sumas de exponenciales". IEEE Transactions on Communications . 68 (10): 6514–6524. arXiv : 2007.06939 . doi :10.1109/TCOMM.2020.3006902. S2CID  220514754.
  16. ^ Tanash, IM; Riihonen, T. (2020). "Coeficientes para aproximaciones minimax globales y límites para la función Q gaussiana mediante sumas de exponenciales [conjunto de datos]". Zenodo . doi :10.5281/zenodo.4112978.
  17. ^ Karagiannidis, GK; Lioumpas, AS (2007). "Una aproximación mejorada para la función Q gaussiana" (PDF) . IEEE Communications Letters . 11 (8): 644–646. doi :10.1109/LCOMM.2007.070470. S2CID  4043576.
  18. ^ Tanash, IM; Riihonen, T. (2021). "Coeficientes mejorados para las aproximaciones de Karagiannidis–Lioumpas y límites para la función Q gaussiana". IEEE Communications Letters . 25 (5): 1468–1471. arXiv : 2101.07631 . doi :10.1109/LCOMM.2021.3052257. S2CID  231639206.
  19. ^ Chang, Seok-Ho; Cosman, Pamela C .; Milstein, Laurence B. (noviembre de 2011). "Límites de tipo Chernoff para la función de error gaussiana". IEEE Transactions on Communications . 59 (11): 2939–2944. doi :10.1109/TCOMM.2011.072011.100049. S2CID  13636638.
  20. ^ Winitzki, Sergei (2003). "Aproximaciones uniformes para funciones trascendentales" . Computational Science and Its Applications – ICCSA 2003. Lecture Notes in Computer Science. Vol. 2667. Springer, Berlín. págs. 780–789. doi :10.1007/3-540-44839-X_82. ISBN. 978-3-540-40155-1.
  21. ^ Zeng, Caibin; Chen, Yang Cuan (2015). "Aproximaciones globales de Padé de la función generalizada de Mittag-Leffler y su inversa". Cálculo fraccionario y análisis aplicado . 18 (6): 1492–1506. arXiv : 1310.5592 . doi :10.1515/fca-2015-0086. S2CID  118148950. De hecho, Winitzki [32] proporcionó la denominada aproximación global de Padé.
  22. ^ Winitzki, Sergei (6 de febrero de 2008). "Una aproximación práctica para la función de error y su inversa".
  23. ^ Press, William H. (1992). Recetas numéricas en Fortran 77: el arte de la computación científica . Cambridge University Press. pág. 214. ISBN 0-521-43064-X.
  24. ^ Dia, Yaya D. (2023). "Integrales incompletas aproximadas, aplicación a la función de error complementaria". Revista electrónica SSRN . doi :10.2139/ssrn.4487559. ISSN  1556-5068.
  25. ^ abc Cody, WJ (marzo de 1993), "Algoritmo 715: SPECFUN: un paquete FORTRAN portátil de rutinas de funciones especiales y controladores de pruebas" (PDF) , ACM Trans. Math. Softw. , 19 (1): 22–32, CiteSeerX 10.1.1.643.4394 , doi :10.1145/151271.151273, S2CID  5621105 
  26. ^ Zaghloul, MR (1 de marzo de 2007), "Sobre el cálculo del perfil de línea de Voigt: una única integral propia con un integrando de seno amortiguado", Monthly Notices of the Royal Astronomical Society , 375 (3): 1043–1048, Bibcode :2007MNRAS.375.1043Z, doi : 10.1111/j.1365-2966.2006.11377.x
  27. ^ John W. Craig, Un resultado nuevo, simple y exacto para calcular la probabilidad de error para constelaciones de señales bidimensionales Archivado el 3 de abril de 2012 en Wayback Machine , Actas de la Conferencia de Comunicación Militar IEEE de 1991, vol. 2, págs. 571–575.
  28. ^ Behnad, Aydin (2020). "Una nueva extensión de la fórmula de función Q de Craig y su aplicación en el análisis del rendimiento de EGC de doble rama". IEEE Transactions on Communications . 68 (7): 4117–4125. doi :10.1109/TCOMM.2020.2986209. S2CID  216500014.
  29. ^ Carslaw, HS ; Jaeger, JC (1959). Conducción de calor en sólidos (2.ª ed.). Oxford University Press. pág. 484. ISBN 978-0-19-853368-9.
  30. ^ "math.h - declaraciones matemáticas". opengroup.org . 2018 . Consultado el 21 de abril de 2023 .
  31. ^ "Funciones especiales – Documentación de GSL 2.7".

Lectura adicional

Enlaces externos