5 votos

De Euler-Maclaurin suma de $e^{-x^2}$

Quiero aproximado de la suma $$\sum_{k=0}^\infty e^{-k^2}$$ el uso de Euler-Maclaurin fórmula $$\sum_{k=0}^\infty f(k) = \int_0^\infty f(x) \, dx + \frac{1}{2}(f(0) + f(\infty)) + \frac{1}{12}(f'(\infty) - f'(0)) - \frac{1}{720}(f'''(\infty) - f'''(0)) + \ldots$$

donde$f(\infty)$$\lim_{x \to \infty} f(x)$. Pero poner $f(x) = e^{-x^2}$, no es demasiado difícil ver que todos los derivados de orden impar son una extraña polinomio veces $f$, lo que implica que se desvanecen en $x = 0$. Pero esto significa que todos los términos de Euler-Maclaurin fórmula desaparecer, excepto para los dos primeros! Pero yo que esto está mal porque he evaluado la suma numérica en Mathematica.

Lo que está mal con mis cálculos?

8voto

Chris Benard Puntos 1430

No hay nada malo con los detalles de la computación, pero que están tratando de Euler-Maclaurian fórmula como la igualdad, cuando no lo es.

Antes de continuar, voy a cambiar la suma de un poco: se está trabajando con $$S:=\sum_{k=0}^{\infty} e^{-k^2}.$$ Las fórmulas son un poco más limpio en términos de $$T:=\sum_{k=-\infty}^{\infty} e^{-k^2}.$$ Estos están relacionados por $T=2S-1$, por lo que es fácil cambiar entre ellos.

El de Euler-Macluarin fórmula con el resto término, combinado con su correcto cálculo que $f^{(k)}(x)$$0$$x \to \pm \infty$, nos dice que tenemos $$\sum_{k=-\infty}^{\infty} e^{-k^2} = \int_{x=-\infty}^{\infty} e^{-x^2} dx + \int_{-\infty}^{\infty} \frac{B_N(x-\lfloor x \rfloor)}{N!} \frac{d^N e^{-x^2}}{(dx)^N} dx$$ para cualquier entero positivo $N$. Aquí $B_N$ $N$- th Bernouli polinomio y $\lfloor x \rfloor$ $x$ redondea a un número entero. Para algunas funciones, el resto de los $\int \frac{B_N(x-\lfloor x \rfloor)}{N!} \frac{d^N f}{(dx)^N} dx$$0$$N \to \infty$, pero a menudo no lo hace, y no en su caso.

Para dar algunos otros ejemplos, si calcular el asymptotics de $\sum_{k \leq M} \frac{1}{k}$ el uso de Euler-Maclaurin, el de Euler-Mascheroni constante se produce como límite a la del resto término. Si usted derivar la fórmula de Stirling por Euler-Maclaurin suma de $\sum \log k$, $\log (2 \pi)$ se produce como límite a la del resto término.

Hay una muy buena explicación de esto en el Capítulo 9 de Concreto de las Matemáticas, por Graham, Knuth y Patashnik. Usted puede disfrutar particularmente el cuarto ejemplo, en el Capítulo 9.6, donde el uso de Euler de Maclaurin suma de demostrar que, para cualquier $N$, tenemos $$\sum_{k=-\infty}^{\infty} e^{-k^2/t} = \sqrt{\pi t} + O(t^{-N}) \ \mbox{as}\ t \to \infty.$$

4voto

Marko Riedel Puntos 19255

Muchos ejemplos de Euler-Maclaurin de totalización son en realidad armónica sumas y puede ser tratada por Mellin transformar los métodos.

En el presente caso poner $$S(x) = \sum_{k\ge 1} e^{-x^2 k^2}$$ con lo que estamos interesados en $S(1/\sqrt{t})$ como $t\rightarrow\infty.$

Esta suma puede ser evaluado por la inversión de sus Mellin transformar.

Recordar que el armónico suma de identidad $$\mathfrak{M}\left(\sum_{k\ge 1} \lambda_k g(\mu_k x);s\right) = \left(\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} \right) g^*(s)$$ donde $g^*(s)$ es la transformada de Mellin $g(x).$

En el presente caso tenemos $$\lambda_k = 1, \quad \mu_k = k \quad \text{y} \quad g(x) = e^{-x^2}.$$

Necesitamos la Mellin transformar $g^*(s)$ $g(x)$ que es $$\int_0^\infty e^{-x^2} x^{s-1} dx.$$ El uso de la sustitución de $x^2 = u$, de modo que $2x \; dx = du$ para obtener $$\int_0^\infty e^{-u} u^{1/2s-1/2} \frac{1/2 \; du}{\sqrt{u}} = \frac{1}{2} \int_0^\infty e^{-u} u^{1/2s-1} du = \frac{1}{2} \Gamma(s/2).$$

Los fundamentales de la tira de esta transformación es Mellin $\langle 0,\infty\rangle.$

De ello se desprende que la Mellin transformar $Q(s)$ de la suma de armónicos $S(x)$ está dado por

$$Q(s) = \frac{1}{2} \Gamma(s/2) \zeta(s) \quad\text{porque}\quad \sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} = \sum_{k\ge 1} \frac{1}{k^s} = \zeta(s)$$ para $\Re(s) > 1.$

El Mellin de inversión integral para esta transformación es $$\frac{1}{2\pi i} \int_{3/2-i\infty}^{3/2+i\infty} P(s)/x^s ds$$ que nos evaluar por desplácelo hacia la izquierda para una expansión de alrededor de cero (hay que recordar que como $t\rightarrow\infty$ tenemos $1/\sqrt{t}\rightarrow 0.$)

Observar que los polos están en $s=1$ a partir de la función zeta y a término el no positivo números enteros a partir de la función gamma plazo. Sin embargo todos ellos, excepto el uno a cero se cancela por el trivial los ceros de la función zeta plazo, dejando sólo los polos $s=0$ y $s=1.$

Para estos dos que hemos $$\mathrm{Res}\left(Q(s)/x^s; s=1\right) = \frac{1}{2} \Gamma(1/2)\frac{1}{x} = \frac{\sqrt{\pi}}{2x}$$ y $$\mathrm{Res}\left(Q(s)/x^s; s=0\right) = \frac{1}{2} \times 2 \times -\frac{1}{2} = -\frac{1}{2}.$$

De ello se sigue que como $t\rightarrow\infty$ hemos $$S(1/\sqrt{t}) \sim \frac{1}{2} \sqrt{\pi t} - \frac{1}{2}$$ y, en particular, $$2S(1/\sqrt{t})+1 = \sum_{k=-\infty}^\infty e^{-k^2/t} \sim \sqrt{\pi t}.$$

Como para el término de error si hemos cambiado la integral a la línea de $\Re(s) = -q/2$ $q>1$ $q$ impar tenemos para la norma de la zeta función plazo en la línea de $-q/2+iv$ de los dependientes de $|v|^{1/2+q/2}$ y el la función gamma plazo decae exponencialmente en $v$ a lo largo de las líneas verticales y en $q$ a los valores en $-q/2$, por lo que el término de error decae exponencialmente también.

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