Matemáticas universitarias — Grado 2 · Bachelor Year 2
6Comparación de funciones
Análisis asintótico: el arte de reemplazar una cantidad complicada por uno simple más un error controlado — se inició en el Año 1 volumen con expansiones de Taylor. Este capítulo la convierte en una disciplina de propio: expansiones a lo largo de escalas generales, la serie–integral comparación con su fuerza asintótica total, la fórmula de Stirling (probada completamente), y el estudio sistemático de los conceptos implícitamente definidos. secuencias. Estas técnicas son el pan de cada día de los asintóticos. análisis, y cada capítulo posterior que estima algo — series, integrales, probabilidades — come de esta tabla.
6.1 Relaciones de comparación y escalas.
Definición 6.1
Cerca de un punto a (a∈R o ±∞), para funciones (o secuencias, con n→∞): f=o(g), f=O(g), f∼g como en el volumen del Año 1. Un comparación escala en a es una familia de positivos funciones, comparables por pares, totalmente ordenadas por o(⋅) — la escala estándar en +∞ es
xα(lnx)β(α,β∈R),
ordenado lexicográficamente en (α,β), refinado cuando sea necesario por exponenciales eγx.
Definición 6.2(expansión asintótica)
f admite el expansión asintótica
f=c1φ1+c2φ2+⋯+ckφk+o(φk)(φi+1=o(φi) in the scale)
cuando los restos sucesivos satisfacen las estimaciones mostradas. el los coeficientes son entonces únicos: c1=limf/φ1, y inductivamente ci+1=lim(f−∑j≤icjφj)/φi+1.
Ejemplo 6.3
Las expansiones de Taylor son expansiones asintóticas a lo largo de la escala (x−a)ken a. Pero la noción es estrictamente más amplia: en+∞,
x−lnx1=x1⋅1−xlnx1=x1+x2lnx+o(x2lnx),
una expansión a lo largo de la escala mixta — no se aplica ningún teorema de Taylor, sólo la expansión geométrica y el cálculo de o.
Ejemplo 6.4(La báscula estándar está realmente ordenada)
El reclamo lexicográfico de Definición 6.1 necesita una línea de prueba por caso. Comparar xα(lnx)β y xα′(lnx)β′ en +∞. Si α<α′: la relación es xα−α′(lnx)β−β′→0, porque un La potencia negativa de x aplasta cualquier potencia de lnx (establecer x=et: e(α−α′)ttβ−β′→0 por el límite polinómico de latidos exponenciales del año 1 volumen). Si α=α′ y β<β′: la relación es (lnx)β−β′→0 directamente. Entonces las parejas (α,β), ordenado lexicográficamente, ordena la escala por o(⋅) — y la sustitución x=et es la truco único para comparaciones de registros de potencia mixtos.
Ejemplo 6.5(Clasificación de una colección de animales.)
Las escalas deben ser ordenado; Aquí está el taladro estándar. en +∞, comparar n10, elnn⋅n, 2n y nlnn tomando logaritmos:
10lnn≪(lnn)2≪nlnn≪nln2,
donde an≪bn significa an=o(bn); la segunda entrada es ln(nlnn). Los exponenciales preservan estas brechas estrictas (si lnun−lnvn→−∞ luego un/vn→0), entonces
n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).
La moraleja, doble: comparar siempre mediante logaritmos (diferencias de registros, no proporciones de registros), y nunca concluir un∼vn de lnun∼lnvn — el par n10 y nlnn tiene una relación ln que tiende a ∞, pero 2n y 4n tienen una relación ln exactamente 2 y son tremendamente desigual.
Si ∫1∞f diverge, las sumas parciales satisfacen ∑k=1nf(k)=∫1nf+C+o(1) para algunos constante C: la diferencia ∑k≤nf(k)−∫1nfconverge.
Demostración. El bracketing f(k+1)≤∫kk+1f≤f(k) (disminución) fue el dispositivo del Año 1; sumando k≥n+1, resp. k≥n, da (1). Para (2), configure uk=f(k)−∫kk+1f: por el entre paréntesis, 0≤uk≤f(k)−f(k+1), por lo que las sumas parciales de ∑uk están delimitados por el f(1)−f(n+1)≤f(1)telescópico: la serie converge. Además, la secuencia (∫nn+1f)nno es creciente (f disminuye) y no es negativa, por lo tanto convergente. escribiendo
k=1∑nf(k)−∫1nf=k=1∑nuk+∫nn+1f,
el lado derecho converge como n→∞: la diferencia converge a una constante C, que es la declaración (2). ∎
Ejemplo 6.7(La expansión armónica)
Para f(t)=t1: Hn=lnn+γ+o(1), recuperando Constante de Euler (volumen del año 1) con una prueba más limpia. Empujando uno ordenar más (Ejercicio 6.3):
Hn=lnn+γ+2n1+o(n1).
Los números hacen visible la ganancia en n=10: H10=2.928968…yln10=2.302585…, por lo que la estimación bruta de γ es H10−ln10=0.626383, desviado por 0.049; restando la corrección 201 se obtiene 0.576383, fuera de γ=0.577216 por solo 8.3⋅10−4 — que es a su vez el siguiente término 12⋅1001 del expansión, como lo demuestra el problema del fin de semana (pregunta 8).
Ejemplo 6.8(Un ln(n!) tosco sin Stirling)
El dispositivo de fijación por sí solo ya localiza ln(n!). desde ln aumenta,
∫k−1klntdt≤lnk≤∫kk+1lntdt,
y sumando k=2,…,n (con ∫1nln=nlnn−n+1):
nlnn−n+1≤ln(n!)≤(n+1)ln(n+1)−n.
Ambas vallas son nlnn−n+O(lnn): de ahí ln(n!)=nlnn−n+O(lnn), y en particular ln(n!)∼nlnn. Lo que Stirling agrega son los siguientes dos peldaños: el 21lnn y la constante ln2π — que cuesta más fino telescópico de Teorema 6.13. sabiendo La precisión que adquiere cada herramienta es la mitad del arte asintótico.
Ejemplo 6.9(Cancelación exige ampliaciones)
Calcule el límite de n2+n−n. Ambos términos son ∼n y “∼n−n” no tienen sentido: equivalentes no se puede restar. Expandir en su lugar:
límite 21, con la velocidad de aproximación 8n1 como bonificación. El mecanismo merece un nombre: una diferencia de dos grandes cantidades equivalentes vive enteramente en su próximo términos, por lo que uno debe expandirse al primer orden en el que los dos lados difieren — y lleve el resto para certificar que nada más sobrevive en ese orden.
Para f(t)=tlnt1 en [2,+∞) (continuo, positivo, decreciente): ∫2xf=lnlnx−lnln2→∞, entonces por Teorema 6.6 (2),
k=2∑nklnk1=lnlnn+C+o(1)
para alguna constante C. Dos lecciones. Primero, la divergencia es real pero glacial: la suma parcial supera por primera vez 4 alrededor de n≈ee4−C, astronómicamente grande. En segundo lugar, el formalnlnn fue entregado por una antiderivada, no adivinado: para términos monótonos, la integral es la canónica dispositivo sumador, y la constante C — como γ de Euler — es la memoria de los términos iniciales.
6.3 La fórmula de Stirling.
Lema 6.11(Integrales de Wallis, revisadas)
Deje Wn=∫0π/2sinntdt. Luego nWnWn−1=2π para n≥1, (Wn)disminuye y Wn∼2nπ.
Demostración. La integración por partes da nWn=(n−1)Wn−2 (n≥2), por lo que nWnWn−1 es constante en n, igual a 1⋅W1W0=2π. Disminución: sinn+1≤sinn en [0,2π]. El apretón, en detalle: monotonicidad da Wn+1≤Wn≤Wn−1, y dividiendo por Wn−1>0,
n+1n=Wn−1Wn+1≤Wn−1Wn≤1,
la identidad izquierda de la recurrencia en el índice n+1. ambos los límites tienden a 1: Wn∼Wn−1, de donde
nWn2∼nWnWn−1=2π⟹Wn∼2nπ.
∎
Ejemplo 6.12(Las primeras integrales de Wallis)
De W0=2π, W1=1 y la recurrencia nWn=(n−1)Wn−2:
W2=4π,W3=32,W4=163π,W5=158,W6=325π.
Los índices pares llevan un π, los impares son racionales — los dos productos intercalados de las formas cerradas. Numéricamente W6≈0.4909 contra lo asintótico π/12≈0.5116: en n=6 el equivalente es ya dentro de 5%, y la identidad del producto es exacta en cada n: 6W6W5=6⋅325π⋅158=2π. Las mesas pequeñas como esta son la forma más económica de atrapar un error de álgebra antes de que infecte un argumento asintótico.
Teorema 6.13(Stirling)
n!∼2πn(en)n.
Demostración.Step 1: n!∼Cn(n/e)n for some constant C>0. Conjunto
por la expansión de Taylor de ln(1+n1). Por tanto, la serie ∑(dn−dn+1)converge absolutamente (comparación con ∑n−2), por lo que (dn)converge, digamos a d; exponenciando, n!∼Cn(n/e)ncon C=ed.
Step 2: C=2π via Wallis. El formulario cerrado W2p=4p(p!)2(2p)!⋅2π (de la recurrencia, Cálculo del año 1 rehecho en Lema 6.11 ajuste) se combina con el Paso 1:
la probabilidad de que un paseo aleatorio simétrico regrese a 0 en el momento 2n es ∼πn1 — un anuncio de Capítulo 22.
Observación 6.15(Perspectivas dentro de este volumen)
Cada capítulo cuantitativo que sigue habla de la esencia de este capítulo. idioma. Capítulo 7 clasifica series comparando términos a la escala n−α(lnn)−β — su El problema del fin de semana mapea esa frontera por completo. Capítulo 9 hace lo mismo para integrales impropias, con idéntica escala en la variable continuo. Capítulo 11 calcula radios de convergencia a partir de limsup∣an∣1/n, un equivalente de n-ésimas raíces ejercicio donde Stirling es la clave estándar (nn!∼en, Ejercicio 6.4). y el capítulos de probabilidad en efectivo Stirling directamente: el local estimaciones de Capítulo 22 para coeficientes binomiales son Ejemplo 6.14 y Ejemplo 6.21 palabra por palabra. La asintótica no es un capítulo aquí; es el acento del volumen.
Método 6.16(La lista de verificación de arranque)
Antes de confiar en una expansión iniciada, audite cuatro puntos. (1) La existencia primero: la raíz o secuencia debe estar fijada hacia abajo (monotonicidad, valores intermedios) antes de cualquier expansión — los símbolos sin referentes se expanden bellamente y significan nada. (2) Un pedido por pase: cada sustitución puede sólo se debe confiar en el orden de la estimación ingresada; extrayendo dos nuevos términos de una sola pasada es la fuente clásica de errores coeficientes. (3) Los restos siguen adelante: lleva el o(⋅) a través de cada paso algebraico y dejar que la absorción (términos más pequeños absorbidos por restos más grandes) ocurren en el final, explícitamente. (4) Auditoría numérica: evaluar a la vez valor honesto de n; un error de coeficiente sobrevive al algebraico La rederivación es sorprendentemente frecuente y casi nunca sobrevive. aritmética.
Observación 6.17(Errores comunes)
(i) Equivalentes agregar mal: de un∼n+lnn y vn∼−n se puede no concluir un+vn∼lnn; las cancelaciones exigen ampliaciones con restos explícitos, nunca desnudos equivalentes. (ii) Nunca exponenciar una equivalencia: n+1∼n pero en+1∼en; la caja fuerte dirección es tomar logaritmos de equivalentes que tienden a +∞ (problema de fin de semana de este capítulo, pregunta 24). (iii) Un expansión asintótica está conectado a un escala: escritura f=x1+o(x21) reclama más de f=x1+o(x1) y mezclar los dos invalida álgebra posterior. (iv) En bootstraps, sustituir el entero expansión actual, resto incluido — cayendo un pase medio o(⋅) produce resultados plausibles pero incorrectos. coeficientes. (v) La serie: la comparación integral necesita monotonicidad: para términos oscilantes falla rotundamente (compárese ∑ksink, Capítulo 7).
Ejemplo 6.18(Stirling en números)
En n=10: la fórmula da 20π(10/e)10≈3598696 contra 10!=3628800: error relativo 8.3⋅10−3, notable por una declaración "asintótica" en n=10. El error tiene una estructura — el refinamiento exacto n!=2πn(n/e)n(1+12n1+O(n−2))— cuya primera corrección 1201≈8.3⋅10−3 explica casi exactamente la brecha observada. el La maquinaria Euler-Maclaurin del problema del fin de semana es precisamente la fuente sistemática de tales términos de corrección.
Observación 6.19(Dónde se utiliza este capítulo)
La comparación asintótica es la gramática de todo lo cuantitativo. downstream: las pruebas de convergencia y el panorama Bertrand de Capítulo 7, los criterios de integrabilidad de Capítulo 9, los cálculos del radio de convergencia de Capítulo 11, y los teoremas de límite de Capítulo 22 (donde Stirling dirige el de Moivre–estimaciones de Laplace). El volumen del Año 3 industrializa el una idea que probamos a mano aquí: extraiga el término principal, cota el resto — al método de Laplace y dominó la convergencia.
Ejemplo 6.20(Una integral comparada consigo misma: ∫2xlntdt)
La caja de herramientas de comparación también funciona con integrales. Sea F(x)=∫2xlntdt (el integrando es continuo en [2,∞)). Integrar por partes:
Por lo tanto F(x)∼lnxx. Lectores que conocieron al mejor teorema de números en el problema de fin de semana de este capítulo reconocerá F: es la integral logarítmica, el mejor estimador de π(x), y el cálculo muestra que concuerda con lnxx al primer pedido.
Ejemplo 6.21(Stirling en un binomio torcido)
La misma rutina trifactorial que para Ejemplo 6.14 da, para (n3n)=n!(2n)!(3n)!:
La tasa exponencial 427=2233 es e3nH(1/3) en la notación entrópica de información. teoría: los binomios asimétricos crecen estrictamente más lento que el central 4n por dos pasos — aquí (27/4)1/3≈1.89<2 por paso. Cada binomio asintótico en combinatoria y probabilidad. (Capítulo 22) es este cálculo con diferentes pesos.
6.4 Secuencias definidas implícitamente
Método 6.22
Para encontrar las asintóticas de las soluciones xn de una ecuación F(x,n)=0:
Localizar: demostrar la existencia y unicidad de xn en un intervalo definido (monotonicidad, valor intermedio teorema), y encontrar su comportamiento crudo (límite, orden de crecimiento).
Oreja: sustituye la forma cruda xn=(main term)(1+εn) en la ecuación y resuelve el siguiente orden de εn; repetir, cada pase refinando un pedido.
Ejemplo 6.23
Para n≥1, la ecuación tanx=x tiene exactamente una solución xn en (nπ−2π,nπ+2π) (la función tanx−x aumenta de −∞ a +∞ allí, es siendo el derivado tan2x≥0). Crudo:xn=nπ+2π−yncon yn∈(0,π); ya que xn→∞ytanxn=xn→+∞, xn se acerca al asíntota desde la izquierda: yn→0. Oreja:tanxn=cotyn=tanyn1∼yn1, y la ecuación cotyn=xn∼nπ da yn∼nπ1. Por lo tanto
xn=nπ+2π−nπ1+o(n1),
y el proceso continúa a cualquier orden (Ejercicio 6.6).
Ejemplo 6.24(Una segunda ejecución del método)
Resuelva x+lnx=n asintóticamente. Localizar:x↦x+lnxaumenta de−∞a+∞ en (0,+∞): una raíz única xn y xn→∞. Crudo:lnxn=o(xn) da xn∼n. Oreja: de xn=n−lnxn y lnxn=lnn+o(1)(logaritmos de equivalentes, ambos lados →∞):
xn=n−lnn+o(1);
una pasada más, con lnxn=ln(n−lnn+o(1))=lnn−nlnn+o(nlnn):
xn=n−lnn+nlnn+o(nlnn).
(Verifique en n=100: la raíz es x≈95.4415; la La fórmula de tres términos da 100−4.6052+0.0461=95.4409, el dos términos uno 95.3948 — cada pase gana el pronóstico orden.) Mismo bucle, tercer paisaje: el método de Método 6.22 no le importa lo que el parece la ecuación, sólo que cada pasada aísla la dominante desconocido.
6.5 Ceremonias
Ejercicio 6.1★
Expandir en +∞, dos términos más allá del principal:
x2+x+1,ln(x2+x)−2lnx,x−lnxx+sinx.
Solución
Solución de Ejercicio 6.1.
x2+x+1=x1+x1+x21=x+21+83⋅x1+o(x1) (binomial expansión: 21u−81u2 con u=x1+x21 da 2x1+2x21−8x21=2x1+8x23, luego multiplica por x).
Ordene las aportaciones en la escala +∞: xlnx≫x1≥xsinx≫x2(lnx)2. Por lo tanto, los dos términos después del 1 inicial son xlnx, luego el término de oscilación acotada xsinx:
x−lnxx+sinx=1+xlnx+xsinx+O(x2(lnx)2).
Ejercicio 6.2★
Dé la naturaleza (convergencia/divergencia) y, cuando sea divergente, la asintóticas principales de ∑k≤nkα para α>−1, α=−1, α<−1, vía Teorema 6.6.
Solución
Solución de Ejercicio 6.2.
f(t)=tα (t≥1).
α>−1: divergencia, y por Teorema 6.6 (2), ∑k≤nkα=α+1nα+1+C+o(1)si α<0(donde fdisminuye); para α≥0(f en aumento) el El mismo corchetes con desigualdades invertidas da ∑k≤nkα∼α+1nα+1.
Para n≥2, demuestre que xn+x=1 tiene una solución única xn∈(0,1), que xn→1, y establezca
xn=1−nlnn+o(nlnn).
(From xnn=1−xn: take logarithms and bootstrap with xn=1−εn.)
Solución
Solución de Ejercicio 6.5.
g(x)=xn+x−1 aumenta estrictamente en [0,1] desde −1 a 1: raíz única xn. Desde xnn=1−xn∈(0,1): si xn≤c<1a lo largo de una subsecuencia, entonces xnn≤cn→0, entonces 1−xn→0: contradicción con xn≤c. Por lo tanto xn→1.
Escriba xn=1−εn, εn→0+. el la ecuación dice (1−εn)n=εn, es decir
nln(1−εn)=lnεn⟹−nεn(1+o(1))=lnεn.
Entonces nεn=−lnεn(1+o(1))→+∞, y tomando logaritmos nuevamente: lnn+lnεn=ln(−lnεn)+o(1). Desde ln(−lnεn)=o(ln(1/εn)), esto da lnεn∼−lnn, de donde εn=n−lnεn(1+o(1))∼nlnn:
Determinar limn→∞n!1∑k=0nk!(limita la suma de todos los términos excepto los dos últimos), y deducir el expansión asintótica∑k≤nk!=n!(1+n1+O(n−2)).
Entonces 1≤n!1∑k!≤1+n2: el límite es 1. Refinación: n!(n−1)!=n1 y el crudo consolidado ∑k≤n−2k!≤(n−1)! se puede afilar de la misma manera: ∑k≤n−2k!=(n−2)!(1+O(n1))=O(n2n!). Por lo tanto
k=0∑nk!=n!(1+n1+O(n21)).
Ejercicio 6.8★★★
Sean u0>0 y un+1=un+un1. Demuestre que un→∞, luego que un∼2n(study un2: its increments are 2+un−2; sum) y refine:
un=2n(1+8nlnn+o(nlnn)).
(From un2=2n+∑k<nuk−2+u02 and uk2∼2k: the sum is ∼21lnn by Teorema 6.6.)
Solución
Solución de Ejercicio 6.8.
(un) aumenta; si está acotado, convergería a ℓ con ℓ=ℓ+ℓ1: absurdo. Entonces un→∞.
Cuadrados: un+12=un2+2+un−2, entonces
un2=u02+2n+k=0∑n−1uk21.
La suma es o(n) (los términos tienden a 0, Cesàro), por lo que un2∼2n y un∼2n.
Refinamiento: uk21∼2k1, así que en comparación (Teorema 6.6, o equivalentes de parcial sumas de series positivas) ∑k<nuk−2∼21lnn. Por lo tanto
(Una suma de Riemann con un giro) Determinar el comportamiento asintótico de
Sn=k=1∑nn+klnn1.
(Factor n: Sn=n1∑k(1+nklnn)−1; recognize a Riemann-type sum with a slowly varying parameter t=lnn, compute ∫011+tudu=tln(1+t), and conclude Sn∼lnnlnlnn.)
Solución
Solución de Ejercicio 6.9.
Factorice n y configure t=lnn:
Sn=n1k=1∑n1+tnk1.
Para t fijo, la suma es una suma de Riemann de u↦1+tu1 sobre [0,1]; la función es monótona en u, por lo que la La suma de Riemann está entre paréntesis por la integral desplazada una malla:
∫011+tudu−n1≤Sn≤∫011+tudu+n1
(comparación de las sumas de Riemann de una función monótona con su integral, válido para cada n con su propio t=lnn). ahora ∫011+tudu=tln(1+t), y n1=o(tlnt): por lo tanto
Sn=lnnln(1+lnn)+O(n1)∼lnnlnlnn.
Ejercicio 6.10★
Demuestre la identidad (lnn)lnn=nlnlnn y luego clasifique el siguiendo en orden creciente o(⋅) al infinito, con pruebas: n2, (lnn)lnn, 2n, n!, nn.
Solución
Solución de Ejercicio 6.10.
Identidad: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn. Clasificación: comparar logaritmos. ln(n2)=2lnn; ln((lnn)lnn)=lnnlnlnn; ln(2n)=nln2; ln(n!)=nlnn−n+O(lnn) (Stirling, o el más crudo entre corchetes lnn!∼nlnn); ln(nn)=nlnn. desde 2lnn=o(lnnlnlnn), lnnlnlnn=o(n), nln2=o(nlnn−n)y nlnn−n∼nlnnpero n!/nn→0(la diferencia de registros es−n+O(lnn)→−∞):
n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).
(Para cada paso: la diferencia de logaritmos tiende a +∞, por lo que la relación tiende a 0.)
Ejercicio 6.11★★
(Cola de ∑1/k2, dos términos) Usando el telescopio exacto ∑k>nk(k+1)1=n+11 y el descomposición k21=k(k+1)1+k2(k+1)1, probar
k>n∑k21=n1−2n21+O(n31).
Solución
Solución de Ejercicio 6.11.
Descomponer k21=k(k+1)1+k2(k+1)1 y suma para k>n:
k>n∑k21=n+11+k>n∑k2(k+1)1,
la primera suma se telescopa exactamente (k(k+1)1=k1−k+11). Por el segundo: k2(k+1)1=k31+O(k41)(desde k2(k+1)1−k31=k3(k+1)−1), y por el comparación integral ∑k>nk31=2n21+O(n31), ∑k>nk41=O(n31). Por lo tanto
Sumar vk+1−vk=1+O(1) primero da vn=n+O(n), por lo tanto vn≥cn eventualmente; resumiendo con 2vk1=O(k1) da vn=n+O(lnn). uno mas pasar: 2vk1=2k1(1+O(klnk)), entonces
6.6 Problema: Bootstrapping, desde Euler-Maclaurin hasta el primos
Una cantidad implícita o acumulada rara vez entrega su asintóticos a la vez; uno los extrae en pasadas, cada pasada introducir la estimación anterior en la relación definitoria. Este fin de semana se entrena un problema que gira en torno a nuevas ecuaciones, lo demuestra el Fórmula de Euler-Maclaurin de primer orden (el trapezoide Actualización de la serie: comparación integral, con error riguroso. barras), invierte xlnx=n y cobra el método más famoso. comprobar: del teorema de los números primos admitido, la asintótica ley pn∼nlnn del n-ésimo primo.
Problema 6.1
Problema de fin de semana — el Euler–Maclaurin corrección y las asintóticas del n-ésimo primo
Parte I — The bootstrap loop on a fresh equation.
Demostrar la afirmación de unicidad de Definición 6.2: si f=∑i≤kciφi+o(φk)=∑i≤kci′φi+o(φk) está en la misma escala, entonces ci=ci′ para todos los i. Luego empuja el curso mixto ejemplo un peldaño más allá:
x−lnx1=x1+x2lnx+x3(lnx)2+o(x3(lnx)2)(x→+∞),
y explique por qué no aparece ningún término x2c.
Demuestre que para cada n≥1 la ecuación ex+x=ntiene exactamente una solución real xn, y que xn→+∞con xn∼lnn.
Arranca dos veces:
xn=lnn−nlnn−2n2(lnn)2+o(n2(lnn)2).
Verifique numéricamente en n=1000: compare x1000≈6.90083 con los de uno, dos y tres términos valores de la pregunta 3, a cinco decimales.
Parte II — Euler–Maclaurin, order one.
Demuestre la identidad del núcleo trapezoidal: para g de la clase C2 en [0,1],
∫01g(t)dt=2g(0)+g(1)−21∫01t(1−t)g′′(t)dt
(integrate 21t(1−t)g′′ by parts twice).
Sea fC2 en [1,+∞) con ∫1∞∣f′′∣<∞. mostrar eso
En=k=1∑nf(k)−∫1nf−2f(1)+f(n)
converge a una constante E, con la cola unida ∣E−En∣≤81∫n∞∣f′′∣: el Fórmula de Euler-Maclaurin al primer pedido.
Extrae el siguiente coeficiente: show εn=−12n21+o(n21)(the increments of En are 21∫01t(1−t)f′′(n+t)dt=121f′′(n)+o(f′′(n)); sum the tail with Teorema 6.6).
Aplicar la pregunta 6 a f=ln: volver a derivar en tres líneas la convergencia de dn=lnn!−(n+21)lnn+n (Paso 1 de Teorema 6.13), con el error de bonificación tasa dn=d+O(n1).
Aplicar la pregunta 6 a f(t)=t1: mostrar
k=1∑nk1=2n+c+2n1+O(n3/21)
para alguna constante c y evalúe todos los términos en n=104(la constante es c≈−1.4604).
Parte III — Inversion: the equation xlnx=n.
Demuestre que xlnx=n tiene exactamente una solución xn∈[1,+∞) para n≥1, que xn→∞y que lnxn∼lnn.
Deduce la inversión de un término xn∼lnnn, luego arranca una vez más:
Prueba en n=106: la raíz verdadera es x≈87848; comparar con el de un término (≈72382) y valores de dos términos (≈86140), y explique los ganancia lenta (el parámetro de expansión es lnnlnlnn, solo ≈0.19 en n=106).
Ahora admitir teorema de los números primos: el número π(x) de números primos ≤x satisface a π(x)∼lnxx como x→∞ (probado honestamente en el volumen del año 3). Escribiendo pn para el n-ésimo primo, justifique π(pn)=n y ejecute la inversión de preguntas 11–12 para probar
pn∼nlnn.
Dividendos: (a) mostrar ∑k≤npk∼2n2lnn(compare ∑klnk with ∫tlntdt); (b) calcule el valor aproximado probabilidad de que un entero uniformemente aleatorio con 100 dígitos es primo (ln10100≈230.26: aproximadamente uno en 230).
Parte IV — The method exported: xtanx=1.
Demuestre que para cada n≥1 la ecuación tanx=x1tiene exactamente una solución xn en (nπ,nπ+2π), y que zn=xn−nπ→0+.
Un término: zn∼nπ1.
Demuestre que la expansión de zn tiene non2c término: zn=nπ1+O(n31).
Tres términos: usando arctanu=u−3u3+O(u5)y xn1=nπ1−(nπ)2zn+O(n−3⋅zn2), pruebe
xn=nπ+nπ1−3π3n34+o(n31).
Verifique en n=3: raíz verdadera x3≈9.5293344; compare los valores de uno y tres términos, y contraste en una frase con el tanx=x del curso (Ejemplo 6.23): donde se ubica cada secuencia en su ventana y por qué.
Part V — A dynamical bootstrap, rules of the game, synthesis.
Dejemos u0∈(0,π) y un+1=sinun. Muestre un→0 de forma decreciente y calcule el límite de un+121−un21(expand sin−2 via sinu=u−6u3+o(u3)).
Deduce, vía Cesàro significa (volumen Año 1), el clásico
un∼n3.
(Números certificados) Utilizando el límite riguroso de la pregunta 7, muestra que evaluar lnn+γ+2n1 en n=106 produce H106 con error como máximo 1.25⋅10−13 — una suma millonaria calculada para trece dígitos por tres términos.
(Reglas del juego) Probar o refutar, con pruebas o contraejemplos: (a) si un∼vn→+∞ entonces lnun∼lnvn; (b) si un∼vn entonces eun∼evn; (c) si f∼g en +∞ (f,g diferenciable) y luego f′∼g′.
(Síntesis) En una frase cada uno: el bucle bootstrap de Método 6.22 como se usa en las Partes I, III, IV; ¿Qué aporta la corrección trapezoidal? Teorema 6.6; por qué inversión de xlnx es exactamente el puente de π(x) a pn; y cuál de las reglas de la pregunta 24 protegía cuál paso. Nombra las dos cumbres: la fórmula Euler-Maclaurin (primer orden), y la ley asintótica del n-ésimo prima.
Solución
Solución de Problema 6.1.
1. Restando las dos expansiones: ∑i(ci−ci′)φi=o(φk). Si algún coeficiente difiere, sea i0 sea el primero: dividir por φi0 y usar φj=o(φi0) para j>i0 da ci0−ci0′=o(1): cero, contradicción. Para la ampliación: con u=xlnx→0,
No aparece ningún término x2c porque la expansión es una serie geométrica en u=xlnx: cada término lleva como muchos poderes de lnx a partir de x1 más allá del primero; el El peldaño de escala x21 (coeficiente de (lnx)0) es simplemente ausente, con coeficiente 0.
2.f(x)=ex+x es continuo, estrictamente creciente, con límites −∞ y +∞: una biyección R→R, por lo que xn=f−1(n) existe y es único, y xn→+∞ (f−1 aumenta a +∞). De exn=n−xn: xn=ln(n−xn)≤lnn, entonces xn/n→0 y xn=lnn+ln(1−xn/n)=lnn+o(1)∼lnn.
3. Escribe un=xn/n. Segundo pase: un=nlnn+o(1), entonces
4. En n=1000: ln1000≈6.90776 (error 7⋅10−3); dos términos: 6.90085 (error 2⋅10−5); tres términos: 6.90082 (error a continuación 10−5), contra x1000≈6.90083. Cada pase compra aproximadamente el factor previsto nlnn.
5. Dos integraciones por partes, empezando por la derecha: con dtd[21t(1−t)]=21−tyt(1−t) desapareciendo en ambos extremos,
Desde 0≤t(1−t)≤41: ∣En+1−En∣≤81∫nn+1∣f′′∣, cuya suma sobre n converge por hipótesis: (En) converge (incrementos absolutamente sumables) a algunos E, con
con c=E−23. En n=104: 2n=200, c≈−1.46035, 2n1=0.005: previsto 198.54465 y, de hecho, ∑k≤104k−1/2=198.544645… — tres términos, siete dígitos.
11.t↦tlnt es continuo y estrictamente aumentando en [1,∞) (derivada lnt+1≥1), de 0 a +∞: existe un xn único y un xn→∞(de lo contrario, xnlnxn permanecería limitado). tomando logaritmos en xnlnxn=n: lnxn+lnlnxn=lnn; desde lnlnxn=o(lnxn), al dividir por lnxn se obtiene lnxnlnn→1: lnxn∼lnn.
12. De xn=lnxnn y lnxn∼lnn: xn∼lnnn. Siguiente paso: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1), entonces lnxn=lnn−lnlnn+o(1) y
13. En n=106: lnnn≈72382 (desactivado por 18%), dos términos dan ≈86140 (desactivado por 1.9%), contra el verdadero x≈87848. La ganancia por pasa solo el factor lnnlnlnn≈13.82.63≈0.19: las escalas logarítmicas convergen con una lentitud enloquecedora: una realidad dondequiera que se encuentren los números primos. involucrados.
14. Hay exactamente n números primos ≤pn (es decir, p1,…,pn): π(pn)=n. El teorema de los números primos (admitido; volumen del año 3) da n=π(pn)∼lnpnpn, es decir pn∼nlnpn: este es el ecuación xlnx≈n leída al revés. Tomando logaritmos: lnpn=lnn+lnlnpn+o(1) y lnlnpn=o(lnpn)fuerzan a lnpn∼lnn como en la pregunta 11. Sustituyendo de nuevo:
pn∼nlnpn=nlnnlnnlnpn∼nlnn.
15. (a) Reparar ε>0; para k grande, (1−ε)klnk≤pk≤(1+ε)klnk. Por comparación con el creciente tlnt (horquillado tipo Teorema 6.6), ∑k≤nklnk=∫1ntlntdt+O(nlnn)=2n2lnn−4n2+O(nlnn)∼2n2lnn. Por lo tanto ∑k≤npk=2n2lnn(1+O(ε)+o(1))para cada ε: ∑k≤npk∼2n2lnn. (b) Por la prima teorema de los números, entre los números enteros hasta 10100 a proporción ∼ln101001=230.26…1 son primos: un entero uniformemente aleatorio de dígitos 100 es primo con probabilidad aproximadamente 2301.
16. En (nπ,nπ+2π), g(x)=tanx−x1es continuo y aumenta estrictamente (g′=1+tan2x+x21>0), con g→−nπ1<0 en el extremo izquierdo y g→+∞ en el derecho: exactamente uno raíz xn. Desde tanzn=tanxn=xn1→0 con zn∈(0,2π): zn→0+.
17.tanzn∼zn y xn1∼nπ1: zn∼nπ1.
18.zn=arctanxn1 y arctanu=u+O(u3). Con zn=O(n1):
entonces zn=nπ1+O(n31): el El renglón n2c lleva el coeficiente 0, porque el La primera corrección a xn1 es en sí misma de tamaño. n2zn=O(n−3).
19. Inserte zn=nπ1+O(n−3) en el visualización anterior:
xn1=nπ1−n3π31+O(n51),
entonces zn=arctanxn1=xn1−31(xn1)3+O(n51)=nπ1−n3π31−3n3π31+O(n51):
xn=nπ+nπ1−3π3n34+O(n51).
20. En n=3: un término 9.53088, tres términos 9.52929, raíz verdadera 9.52933: errores 1.5⋅10−3 y 5⋅10−5. Contraste: para tanx=x la raíz debe hacer tan enorme, por lo que abraza el extremo biennπ+2π de la ventana, a distancia ∼nπ1 antes de la asíntota; para xtanx=1 la raíz debe hacer tan pequeño, por lo que se encuentra justo después del extremo izquierdanπ, a la distancia ∼nπ1 después del cero. Mismo método, Geografía espejo.
21.sinu<u en los mapas (0,π) y sin(0,π) en (0,1]⊆(0,π): después del uno paso u1∈(0,1], luego (un) disminuye y es acotado abajo por 0: converge a un punto fijo de sin, es decir, a 0. Expansión: sinu=u(1−6u2+o(u2)), entonces
entonces un2∼n3 y, siendo todos los términos positivos, un∼3/n.
23. Por la pregunta 7, Hn−lnn−γ−2n1≤8n21. En n=106 este límite es 8⋅10121=1.25⋅10−13: tres computados términos entregan la suma armónica de un millón de términos a trece dígitos, con un certificado de error totalmente riguroso — el conjunto punto de una fórmula asintótica con resto explícito.
24. (a) Verdadero: lnun−lnvn=lnvnun→0 mientras que lnvn→+∞, por lo que la relación de logaritmos tiende a 1. (b) Falso: un=n+1∼vn=n, pero eun/evn=e=1. La equivalencia tolera errores aditivos o(1) en el exponente, no O(1). (c) Falso: f(x)=x+sin(x2)∼g(x)=x en +∞, pero f′(x)=1+2xcos(x2)oscila ilimitadamente mientras que g′=1: Las derivadas de funciones equivalentes no necesitan ser comparables en todos.
25. Se ejecutó el bucle de Método 6.22. idénticamente tres veces: localizar la raíz, extraer un crudo plazo, retroalimentarlo para el siguiente pedido — en ex+x=n (Parte I), en xlnx=n (Parte III), en xtanx=1 (Parte IV). La corrección trapezoidal actualiza la serie: integral comparación de "la diferencia converge" a una definición explícita Término 2f(1)+f(n) con un resto certificado O(∫n∞∣f′′∣) — constantes y barras de error en lugar de mera convergencia. El puente hacia los números primos es pura inversión: el El teorema de los números primos dice π(x)lnx≈x, por lo que pn, definido por π(pn)=n, resuelve una ecuación xlnx=n — y hereda sus asintóticas. Regla a) de la pregunta 24 legitimó cada pasaje de un∼vn a lnun∼lnvn (preguntas 11, 14); la falsedad de (b) es la razón por la que nunca exponenciar equivalencias. Cumbres: Euler-Maclaurin fórmula de primer orden (pregunta 6), y la ley asintótica pn∼nlnn del n-ésimo primo (pregunta 14).