4 votos

Aproximación de la función erf

Estaba tratando de encontrar una solución aproximada a lo siguiente:

$\DeclareMathOperator\erf{erf}$

$$\frac12 \sqrt{\pi} \erf\left (\frac{x-2}{\sqrt{10}}\right) + \frac12 \sqrt{\pi} \erf \left(\frac{x+2}{\sqrt{10}}\right) = \frac25 \sqrt{\pi}$$

Esto equivale naturalmente a

$$\int_0^{\frac{x-2}{\sqrt{10}}} e^{-t^2} dt+ \int_0^{\frac{x+2}{\sqrt{10}}} e^{-t^2} dt = \frac25 \sqrt\pi$$

Lo que intenté entonces fue utilizar la aproximación de taylor y luego resolver las ecuaciones polinómicas.. pero los resultados que obtuve no se acercaban consistentemente a la solución correcta ( $x \sim 1.71$ ), obtenido con wolframio alfa

De hecho, utilizando $\displaystyle \int_0^x e^{-t^2}dt = x - \frac{x^3}{3} + \frac{x^5}{10} - \frac{x^7}{42} + \dots$

He obtenido las siguientes aproximaciones ( $x_i$ indica el resultado obtenido considerando los términos hasta $x^i$ ) $$x_1 \sim 1.12\ \ \ \ x_3 = -4.975 \ \ \ x_5 = 1.59 \ \ \ x_7 = -4.718 | 3.729 | 1.755 \ \ \ \ x_9 = 1.70$$

Preguntas

1) Si tengo en cuenta toda la serie, ¿qué me asegura que sólo se encontrará una solución (con $x_7$ Por ejemplo, encuentro $3$ soluciones reales)

2) Entiendo que mientras esté cerca $0$ la solución será una buena aproximación. Pero a priori no sé qué valores $x$ va a tomar, por lo que el error puede ser tan grande como $x$ y la aproximación de Taylor es inútil.

3) ¿Por qué (eventualmente) el error va a $0$ si uno toma toda la serie?

Creo que no tengo claro el paso entre el polinomio de Taylor (válido sólo cerca de un punto $x_0$ ) y las series (por qué exactamente son convergentes en todas partes (bueno, al menos en el caso de $e^x$ ) si consideramos una cantidad infinita de términos)

1voto

Steven Lu Puntos 866

$\DeclareMathOperator\erf{erf}$ Aproximaciones de $\erf$ según el Manual de funciones matemáticas :

http://people.math.sfu.ca/~cbm/aands/page_299.htm

0voto

bubba Puntos 16773

Como sospechas, las expansiones en serie de Taylor sólo son precisas cerca de un punto. Para obtener una aproximación polinómica que sea buena en todo un intervalo, la técnica estándar es utilizar la aproximación de Chebyshev. La aproximación más sencilla es simplemente hacer una interpolación en algunos nodos cuidadosamente elegidos (el Nodos Chebyshev ). Si tienes acceso a Matlab, hay un complemento llamado Chebfun que hace un muy buen trabajo para construir este tipo de aproximaciones.

0voto

utilice esto $$\text{Solve}\left[\frac{x \left(-2 \left(\sqrt{2 \left(\sqrt{2}+2\right) \pi } \left(\text{erf}\left(\frac{\sqrt{2-\sqrt{2}}-4}{2 \sqrt{10}}\right)+\text{erf}\left(\frac{\sqrt{2-\sqrt{2}}+4}{2 \sqrt{10}}\right)\right)-\sqrt{2 \left(2-\sqrt{2}\right) \pi } \left(\text{erf}\left(\frac{\sqrt{\sqrt{2}+2}-4}{2 \sqrt{10}}\right)+\text{erf}\left(\frac{\sqrt{\sqrt{2}+2}+4}{2 \sqrt{10}}\right)\right)\right) x^2-\left(\sqrt{2}-1\right) \sqrt{\left(2-\sqrt{2}\right) \pi } \left(\text{erf}\left(\frac{\sqrt{\sqrt{2}+2}-4}{2 \sqrt{10}}\right)+\text{erf}\left(\frac{\sqrt{\sqrt{2}+2}+4}{2 \sqrt{10}}\right)\right)+\left(\sqrt{2}+1\right) \sqrt{\left(\sqrt{2}+2\right) \pi } \left(\text{erf}\left(\frac{\sqrt{2-\sqrt{2}}-4}{2 \sqrt{10}}\right)+\text{erf}\left(\frac{\sqrt{2-\sqrt{2}}+4}{2 \sqrt{10}}\right)\right)\right)}{2 \sqrt{2}}=\frac{2 \sqrt{\pi }}{5},x\right]$$ da x aproximadamente 1,7106147726585592

0voto

marty cohen Puntos 33863

$\DeclareMathOperator\erf{erf}$ Desde $\erf'(x) =\dfrac{2}{\sqrt{\pi}}e^{-x^2/2} $ , una vez que tengas una rutina para $\erf(x)$ , se puede utilizar la iteración de Newton $\left(x_{n+1} =x_n-\dfrac{f(x_n)}{f'(x_n)} \right)$ para encontrar la raíz de $f(x) = 0$ .

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X