Matemáticas universitarias — Grado 2 · Bachelor Year 2
20Integrales de línea e integrales múltiples
Este capítulo extiende la integración de intervalos a curvas y a dominios del plano y del espacio. Integrales de línea integrar un forma diferencialPdx+Qdy a lo largo de una orientación arco; integrales dobles y triples integran funciones sobre dominios bidimensionales y tridimensionales. Las dos teorías se encuentran en el Teorema de Green-Riemann, el fundamental bidimensional teorema del cálculo, y la principal herramienta computacional es el fórmula de cambio de variables, cuyo factor de distorsión es el jacobiano absoluto determinante.
20.1 Integrales de línea
Definición 20.1(Forma diferencial; integral de linea)
Sea U⊆R2abierto. A diferencial forma de grado 1 y clase C0 en U es una expresión ω=Pdx+Qdy con P,Q:U→Rcontinuo — formalmente, un mapa continuo de U al dual de R2, ω(M)=P(M)e1∗+Q(M)e2∗. por un C1 arco γ:[a,b]→U, γ(t)=(x(t),y(t)), el integral de línea de ω junto γ es
∫γω=∫ab(P(γ(t))x′(t)+Q(γ(t))y′(t))dt.
Las definiciones se extienden textualmente a R3 (formas Pdx+Qdy+Rdz) y a los arcos C1 por partes (suma sobre las piezas).
con signo + si θ es creciente (θ(c)=a) y − si disminuye (los límites se intercambian). ∎
Ejemplo 20.3(Trabajo de una fuerza; circulación)
Si F=(P,Q) es un campo de fuerza, ∫γPdx+Qdy=∫ab⟨F(γ(t)),γ′(t)⟩dt es el trabajar de F junto con γ. Para ω=−ydx+xdya lo largo del círculo unitario γ(t)=(cost,sint) en sentido antihorario:
∫γω=∫02π((−sint)(−sint)+costcost)dt=2π,
dos veces el área adjunto — un primer indicio de Green–Riemann.
Ejemplo 20.4(Una integral, dos parametrizaciones, una trampa de señales)
Calcule ∫γxdy a lo largo de la unidad superior semicírculo de (1,0) a (−1,0). Con γ(t)=(cost,sint), t∈[0,π]:
∫0πcost⋅costdt=2π.
Con la parametrización del gráfico x↦(x,1−x2), xde 1a−1 (¡tenga en cuenta la dirección!):
∫1−1x⋅1−x2−xdx=∫−111−x2x2dx=2π
(x=sinu lo reduce a una integral de Wallis). Mismo valor, como garantiza Proposición 20.2 — pero sólo porque ambas corridas van de (1,0) a (−1,0); revertir el viajar invierte el cartel. Cerrando el camino a lo largo del eje x (donde dy=0) no agrega nada, y el total 2π es el área del medio disco: la primera instancia del área límite fórmulas de Green–Riemann a continuación.
Definición 20.5(Formularios exactos y cerrados)
El formulario ω=Pdx+Qdy de la clase C0 es exacto en U si hay f∈C1(U) (a potencial) con ω=df, es decir P=fx y Q=fy. Un formulario C1 es cerrado si Py=Qx está en U.
Teorema 20.6(Teorema fundamental de la recta integrales)
Si ω=df es exacto y γ es por partes C1 arco en U de A a B, luego
∫γω=f(B)−f(A).
En particular, la integral de un forma exacta a lo largo de cualquier arco cerrado es cero y cada formulario exactoC1 está cerrado.
Demostración.dtdf(γ(t))=fx(γ(t))x′(t)+fy(γ(t))y′(t) por la regla de la cadena (Capítulo 15), por lo que el integrando en Definición 20.1 es la derivada de t↦f(γ(t)), y el teorema fundamental del cálculo da la resultado en cada pieza; El telescopio de valores intermedios. cerrazón de las formas exactoC1 es el teorema de Schwarz: Py=fxy=fyx=Qx. ∎
Ejemplo 20.7(Reconstruyendo un potencial)
Deja ω=yexydx+(xexy+2y)dy en R2. Es cerrado: ambas derivadas cruzadas son iguales. exy(1+xy). Para encontrar un potencial, integre P en x en y fijo:
f(x,y)=∫yexydx=exy+c(y),
luego ajuste c haciendo coincidir fy: xexy+c′(y)=xexy+2yda c(y)=y2. Entonces f(x,y)=exy+y2, y para cualquier arco C1 por partes de (0,0) a (1,1),
∫γω=f(1,1)−f(0,0)=(e+1)−1=e,
independientemente de la ruta — la receta de dos pasos (integrar en x, correcto en y) es lo contrario práctico de Teorema 20.6 en dominios donde está formas cerradasexacto.
Ejemplo 20.8(Cerrado no implica exacto)
En U=R2∖{0}, el forma de ángulo
ω=x2+y2−ydx+xdy
está cerrado (cálculo directo: tanto Py como Qx son iguales (x2+y2)2y2−x2), pero es integral a lo largo de la unidad. el círculo es 2π=0 (el mismo cálculo que Ejemplo 20.3, dividido por 1): ω no es exacto en U. Localmente, ω=dθ para una determinación. θ del ángulo polar; el fracaso es global — el ángulo No se puede definir continuamente alrededor del pinchazo. En dominios sin agujeros la patología desaparece: en un en forma de estrella Conjunto abierto, cada formulario C1 cerrado es exacto (Lema de Poincaré, Ejercicio 20.8).
20.2 Integrales dobles
Damos por sentada la teoría de la integral de Riemann en una variable (volumen del año 1 y Capítulo 9) y boceto su versión de dos variables. Una función fcontinuo en un rectángulo R=[a,b]×[c,d] tiene una integral doble ∬Rf, definido por sumas de Riemann sobre cuadrículas exactamente como en una variable, y calculado por iteración:
Demostración. Establezca F(x)=∫cdf(x,y)dy. Uniforme continuidad de f en el compactoR hace Fcontinuo (estimación dominada: ∣F(x)−F(x′)∣≤(d−c)supy∣f(x,y)−f(x′,y)∣). Ahora subdivida [a,b] y [c,d] en n partes iguales, dando un cuadrícula de celdas Rij de áreaΔxΔy. En cada celda, infRijf⋅ΔxΔy≤∫xi−1xi∫yj−1yjf(x,y)dydx≤supRijf⋅ΔxΔy por monotonicidad del Integral de una variable (aplicada dos veces). Sumando las celdas, el La integral iterada ∫abF se aprieta entre las partes inferior y sumas de Riemann superiores de la cuadrícula; por uniforme continuidad ambas sumas convergen al valor común que define ∬Rf como n→∞. El mismo argumento se aplica con los roles de xy y. intercambiadas, por lo que ambas integrales iteradas son iguales a ∬Rf. ∎
Observación 20.10
Continuidad en un rectángulo compacto está haciendo un trabajo real en La prueba de Fubini: proporciona el uniforme continuidad que exprime las sumas de Riemann. Para integrandos más salvajes, el declaración realmente falla — hay funciones cuyos dos las integrales iteradas existen y difieren. El general honesto El teorema, con la integrabilidad como única hipótesis, es El teorema de Fubini para la integral de Lebesgue, demostrado en el Volumen del año 3; Todo en este capítulo se mantiene dentro del continuo configuración donde se encuentra la prueba elemental anterior completo.
Definición 20.11(Dominios elementales)
Un dominio D⊆R2 es y-elemental si
D={(x,y):a≤x≤b,φ1(x)≤y≤φ2(x)}
con φ1≤φ2continuo en [a,b] (x-elemental: simétricamente). Para fcontinuo en un y-elemental D,
∬Df=∫ab(∫φ1(x)φ2(x)f(x,y)dy)dx,
y uno verifica (extendiendo f mediante un argumento de aproximación, o subdividiendo) que cuando D es elemental en ambas direcciones el dos integrales iteradas concuerdan. Dominios cortados en un número finito Las piezas elementales se manejan por aditividad, y las área de D es Area(D)=∬D1.
Ejemplo 20.12
En el triángulo D={0≤x≤1,0≤y≤x}:
∬Dxydxdy=∫01x(∫0xydy)dx=∫01x⋅2x2dx=81.
Intercambiando el orden (x de y a 1): ∫01y(∫y1xdx)dy=∫01y21−y2dy=81 — mismo valor, diferente cálculo: elegir bien el orden de integración es la mitad del artesanía.
Ejemplo 20.13(Cuando solo funciona un pedido)
Calcule I=∫01∫x1ey2dydx. Como está escrito, la integral interna ∫ey2dy no tiene antiderivada elemental: la el cálculo está atascado. Pero el dominio es el triángulo 0≤x≤y≤1, que es elemental en ambas direcciones; cambiando el orden,
I=∫01∫0yey2dxdy=∫01yey2dy=[21ey2]01=2e−1.
La variable interna x no apareció en ninguna parte del integrando, por lo que integrarlo primero produjo exactamente el factor y que hace que la integral exterior sea inmediata. Moraleja: Fubini no es justo una licencia para iterar — es una licencia para elegir, y el orden correcto puede convertir una integral imposible en una una sola línea. Siempre dibuje el dominio y lea ambos. descripciones antes de comenzar.
Teorema 20.14(Cambio de variables)
Sea Φ:U′→U un difeomorfismo C1 entre conjuntos abierto de R2, sea K⊆U un compacto dominio cortado en pedazos elementales con K′=Φ−1(K), y sea fcontinuo en K. entonces
∬Kf(x,y)dxdy=∬K′f(Φ(u,v))∣detJΦ(u,v)∣dudv.
Demostración.Admitido a este nivel.∎
Observación 20.15
La prueba completo — se aproxima a Φ por su diferencial en una cuadrícula fina y controlar las celdas límite — aunque es largo no profundo; se realiza íntegramente en la teoría de medidas del Año 3, como consecuencia de la teoría de Lebesgue. La heurística es la imagen. ya usado para la superficie área: un pequeño cuadrado de lado du en (u,v) se asigna, en primer orden, al paralelogramo abarcado por Φudu y Φvdv, cuyo área es ∣detJΦ∣dudv (Lema 19.18).
Observación 20.16(Método: elegir el cambio de variables)
Tres reflejos cubren la mayoría de los casos. Simetría de la integrando: x2+y2 pide polar, una estructura de producto Requiere mantener los ejes cartesianos. Forma de la límite: límites u(x,y)=c1, v(x,y)=c2 ruegan por las coordenadas (u,v), como en el hiperbólico ejemplo de región a continuación: el dominio se convierte en un rectángulo, que es toda la victoria. estructura lineal: Las expresiones en x+y y x−y invitan al grado 45. rotación o cizallamiento (Ejemplo 20.18). en total casos, tres casillas a marcar antes de integrar: el mapa es un biyección del nuevo dominio sobre el antiguo; es jacobiano se calcula en la dirección realmente utilizada (invertir en el final si es más fácil); y entra el jacobiano con su valor absoluto.
Ejemplo 20.17(Coordenadas polares)
Φ(ρ,α)=(ρcosα,ρsinα) tiene
JΦ=(cosαsinα−ρsinαρcosα),detJΦ=ρ,
entonces dxdy=ρdρdα. Para el disco DR de radio R:
Comparando con el cuadrado [−R,R]2 (que está apretado entre los discos DR y DR2, todos los integrandos positivos) da (∫−∞∞e−x2dx)2=π:
∫−∞+∞e−x2dx=π
— la integral de Gauss nuevamente, ahora por su prueba más famosa (compare la derivación de una variable en Capítulo 9).
Ejemplo 20.18(Cambios afines de variables)
Para un aplicación afínΦ(u,v)=M(u,v)T+C con M invertible, el jacobiano es la matriz constante M: áreas se multiplican por el factor constante ∣detM∣ — la promesa hecha en Capítulo 17 es ahora un teorema. Dos usos inmediatos. La elipse a2x2+b2y2≤1 es la imagen del disco unitario bajo (u,v)↦(au,bv), por lo que área es ab⋅π — sin cálculo. Y para la integral de f(x+y) sobre el cuadrado K=[0,1]2, la cizalla Φ(u,v)=(u−v,v) (determinante1) lo convierte en una integral de f(u) sobre un paralelogramo, que Fubini corta a u constante: con f=exp,
∬Kex+ydxdy=(∫01exdx)2=(e−1)2,
como lo confirma la estructura del producto. Elegir coordenadas adaptado al integrando — no al dominio — es el la otra mitad de la embarcación.
Ejemplo 20.19(Coordenadas adaptadas a un curvilíneo dominio)
Sea D la región del primer cuadrante acotada por el hipérbolas xy=1 y xy=3 y las líneas y=x y y=3x. En las coordenadas u=xy, v=y/x el dominio se convierte en el cuadrado [1,3]×[1,3]; invertir,
x=u/v,y=uv,detJ=xuyv−xvyu=2v1
(un cálculo de dos líneas con x=u1/2v−1/2, y=u1/2v1/2). Por lo tanto
Area(D)=∫13∫132vdudv=2⋅2ln3=ln3≈1.10.
Intentar cortar D en coordenadas cartesianas significa cortarlo en tres piezas con límites hiperbólicos y lineales — factible, triste y propenso a errores. La moraleja se repite Ejemplo 20.18 con toda su fuerza: lea el ecuaciones de contorno, y deje que a ellos elija la coordenadas; el jacobiano luego convierte el curvilíneo celda de malla área, exactamente como lo hizo ρ para polar coordenadas.
Ejemplo 20.20(Valores medios)
El valor medio de f sobre un dominio D es Area(D)1∬Df. Muestra: el promedio distancia desde el centro para un punto uniformemente elegido de la disco de radio R es
πR21∫02π∫0Rρ⋅ρdρdα=πR22πR3/3=32R,
no R/2: área uniforme pone más masa en radios grandes (el el anillo en el radio ρ tiene un peso proporcional a ρ), por lo que el promedio se sitúa más allá de la mitad del camino. Conseguir esto El factor derecho es exactamente el jacobiano polar en acción, y el misma ponderación explica el centroide zˉ=3R/8 del media bola se calcula más adelante en el capítulo en lugar de R/2.
20.3 El teorema de Green-Riemann
Teorema 20.21(Verde–Riemann)
Sea K⊆R2 un dominio compacto que es elemental en ambas direcciones (o una unión finita de tales segmentos pegados), con límite ∂K a trozos C1 cerrado orientado a curvas sinistrórsum (el dominio permanece en la izquierda). Para P,Q de clase C1 en una vecindad de K:
∮∂KPdx+Qdy=∬K(∂x∂Q−∂y∂P)dxdy.
Demostración. Primero, ambos lados son aditivos bajo el corte K a lo largo de un segmentar en dos piezas K1,K2: las integrales dobles suman por aditividad de ∬; en cuanto a las integrales de frontera, la límites en sentido antihorario de K1 y K2 cada recorrido el interior cortado una vez, en direcciones opuesto, por lo que en la suma
∮∂K1+∮∂K2=∮∂K+(cut, both ways)=∮∂K,
las dos pasadas a lo largo del corte se cancelan (Proposición 20.2) y solo el límite exterior sobrevive. Iterando un número finito de cortes, es suficiente tratar un dominio elemental. Probamos ∮Pdx=−∬KPy en un y-dominio elemental D={a≤x≤b,φ1(x)≤y≤φ2(x)}; la identidad ∮Qdy=∬KQx es simétrico (x-elemental), y el teorema es su suma.
Calcular la integral doble de Fubini y la de una variable teorema fundamental:
∬D∂y∂Pdxdy=∫ab(P(x,φ2(x))−P(x,φ1(x)))dx.
Ahora el límite de D, en sentido antihorario, consta de: la parte inferior gráfico y=φ1(x) recorrido de izquierda a derecha, la derecha segmento vertical x=b (hacia arriba), el gráfico superior y=φ2(x) atravesó derecha a izquierda, el vertical izquierdo segmento x=a (hacia abajo). A lo largo de los segmentos verticales x está constante, por lo que contribuyen 0 a ∮Pdx; los gráficos, parametrizado por x, dé
Demostración. Aplique Verde–Riemann a (P,Q)=(0,x), (−y,0) y 21(−y,x): cada vez Qx−Py=1. ∎
Observación 20.23(Elegir entre las tres fórmulas de área)
Las tres fórmulas de límites son iguales, no intercambiables en práctica. Utilice ∮xdy cuando la parametrización haga dy simple (gráficos sobre el eje y), −∮ydx simétricamente, y la media suma simétrico cuando el la parametrización trata a x y y de manera uniforme — para la elipse produjo un integrando constante, sin trigonometría linealización en absoluto. En límites poligonales la media suma se convierte en la fórmula de cordones de zapatos de Ejercicio 20.12, El algoritmo de los topógrafos. Y cuando se traspasa el límite en el sentido de las agujas del reloj según la parametrización dada, las tres fórmulas devuelve menos o área: un resultado negativo no es un error de cálculo pero un informe de orientación — voltee el signo o la parametrización.
Mismo número, dos cálculos muy diferentes — y eso es el uso práctico: cualquier lado de la identidad de Green es más fácil resulta el cálculo, el otro una verificación. Para circulaciones de campos polinomiales alrededor de curvas cerradas, la La integral doble es casi siempre el lado fácil.
Observación 20.26
Green–Riemann explica Ejemplo 20.8: para un forma cerrada (Qx=Py), la integral alrededor del límite de cualquier El dominio contenido en U desaparece. La forma del ángulo no es exacto sólo porque la punción en el origen impide que el disco se limite por el círculo unitario que se encuentra dentro de U — integrales de línea de formas cerradas detecta los agujeros del dominio. (Empujado más, esta observación se convierte en cohomología de De Rham.)
20.4 Integrales triples
La teoría se extiende a tres variables sin idea nueva: Fubini reduce ∭ a tres integrales de una variable (ya sea por rebanar: ∭Kf=∫(∬Kzf)dzsobre los cortes horizontales Kz, o por apilado: integrando en z primero a lo largo de palos verticales), y la formula de cambio de variables se cumple con el 3×3 Jacobiano.
Ejemplo 20.27(Coordenadas cilíndricas y esféricas)
Cilíndrico(x,y,z)=(ρcosα,ρsinα,z): dxdydz=ρdρdαdz. Esférico(x,y,z)=(rcosθcosφ,rsinθcosφ,rsinφ)(θ longitud, Latitud φ∈[−2π,2π]): expandiendo 3×3determinante a lo largo de la última fila,
descargando por fin la fórmula admitida en los capítulos del volumen de los libros anteriores.
Ejemplo 20.28(El tetraedro, dos veces)
El volumen de T={x,y,z≥0,x+y+z≤1}, por apilamiento: para (x,y) fijo en el triángulo x+y≤1, z pasa por encima de [0,1−x−y], por lo que
V=∫01∫01−x(1−x−y)dydx=∫012(1−x)2dx=61.
Por corte: la sección a la altura z es el triángulo {x,y≥0,x+y≤1−z}, de área2(1−z)2, y V=∫012(1−z)2dz=61 nuevamente — el dos cálculos son las mismas integrales en una forma diferente orden, que es todo lo que afirma Fubini. El valor 61=31⋅21⋅1 es la fórmula del cono. (Ejemplo 20.30) con base triangular, y el n versión tridimensional 1/n! se demuestra exactamente con esto Rebanar el problema del fin de semana.
Ejemplo 20.29(Centroide de media bola)
Para la media bola superior H de radio R (z≥0), la La altura del centroide es zˉ=V1∭Hz, con V=32πR3. En coordenadas esféricas (z=rsinφ, φ∈[0,π/2]):
El punto de equilibrio de un hemisferio sólido se encuentra a tres octavos de el radio por encima de la cara plana — por debajo de la media altura R/2, como debe ser, ya que el sólido es más gordo cerca de la base. Cada cálculo de centroide tiene esta forma: un momento Integral, un volumen, una relación y una verificación de plausibilidad. contra la geometría.
Ejemplo 20.30(Volumen por corte: el cono)
Un cono de base áreaA y altura h (ápice hacia arriba, base en z=0): el corte en altura z es la base escalada por el factor (1−z/h), de áreaA(1−z/h)2. Por lo tanto
V=∫0hA(1−hz)2dz=3Ah:
el tercio de las fórmulas escolares, válidas para la base cualquier forma — cortar la convierte en la integral de un cuadrado.
Ejemplo 20.31(Umbrales de integrabilidad en el avión)
¿Por qué α>0 converge ∬Dρ−αdxdyen el disco unitario perforado D (límite sobre los anillos ε≤ρ≤1)? en polares coordenadas,
∫02π∫ε1ρ−αρdρdα=2π∫ε1ρ1−αdρ,
que converge como ε→0 si 1−α>−1, es decir α<2: en la dimensión 2 la singularidad crítica exponente es la dimensión misma, el ρ extra del Jacobiano suavizando la singularidad por un poder. (Asimismo α<3 para una singularidad puntual en el espacio, vía r2.) Una contabilidad radial de este tipo es la forma de integración decidido de un vistazo en el marco de Lebesgue del Año 3 — y es la razón por la que ∭1/r convergió sin esfuerzo en Ejercicio 20.7.
Observación 20.32(Errores comunes)
(i) Orientación: un integral de línea cambia de signo con el dirección de viaje, y Green-Riemann requiere que el límite en sentido antihorario (dominio de la izquierda); para un dominio con un agujero, se atraviesa el límite interior dextrorso. (ii) El jacobiano entra con valor absoluto: cambio de variables nunca produce área negativo, y olvidar ∣det∣ normalmente voltea las señales exactamente cuando el mapa invierte la orientación. (iii) The polar factor ρ: dxdy=ρdρdα, no dρdα — el error más común en todo capítulo; El análisis dimensional lo capta, ya que dρdα tiene la dimensión de un longitud, no de un área. (iv) Integrales dobles impropias: límites superiores Los discos crecientes y los cuadrados crecientes concuerdan aquí porque el los integrandos son positivos (apretando); para cambiar de signo integrandos el límite puede depender del agotamiento, y no La afirmación se hace sin una convergencia absoluta. (v) Dominios versus integrandos: un integrando de producto en un dominio que no es producto, no factoriza la integral — la factorización necesita ambos, como en Ejemplo 20.17 cuadrado.
Observación 20.33(Perspectivas dentro de este volumen)
La integral de Gauss calculada aquí está silenciosamente en todas partes los capítulos de probabilidad: la constante π dentro La fórmula de Stirling (Teorema 6.13) es la integral de este capítulo y, a través de Stirling, fija la 1/πn asintóticas del retorno del paseo aleatorio probabilidades en Capítulo 21. El Integrales de Wallis de El problema del fin de semana vuelve a aparecer ahí también, conduciendo igual. estimaciones binomiales centrales. En la otra dirección, el área y los elementos de volumen de este capítulo completan la geometría de Capítulo 19, y la fórmula de Green vuelve a calcular el sobreáreas de Capítulo 18 (el astroide, en Ejercicio 20.5). Un capítulo, tres servicios: medida para la geometría, constantes para la probabilidad y la disciplina de cambio de variables utilizada por ambos.
20.5 Ceremonias
Ejercicio 20.1★
Calcule ∫γy2dx+xdy a lo largo de: (a) el segmento de (0,0) a (1,1); (b) el arco de parábola y=x2 de (0,0) a (1,1). ¿El formulario es exacto?
Solución
Solución de Ejercicio 20.1.
(a) Segmento γ(t)=(t,t), t∈[0,1]:
∫γy2dx+xdy=∫01(t2+t)dt=31+21=65.
(b) Parábola γ(t)=(t,t2):
∫01(t4⋅1+t⋅2t)dt=51+32=1513.
Los dos valores difieren, por lo que la integral depende de la trayectoria: la forma es noexacto — consistentemente, Py=2y=1=Qx, entonces ni siquiera está cerrado.
Ejercicio 20.2★
Muestra que ω=(2xy+y3)dx+(x2+3xy2+1)dy está cerrado en R2, busque un potencial y calcule ∫γω a lo largo de cualquier arco desde (0,0) hasta (1,2).
Solución
Solución de Ejercicio 20.2.
P=2xy+y3, Q=x2+3xy2+1: Py=2x+3y2=Qx, cerrado el R2. Busque f con fx=P: f=x2y+xy3+g(y); entonces fy=x2+3xy2+g′(y)=Qfuerza a g′(y)=1, diga g(y)=y. entonces
f(x,y)=x2y+xy3+y
es un potencial (R2 tiene forma de estrella, por lo que un potencial tenía que existir según el lema de Poincaré — pero exhibirlo es más rápido). Por Teorema 20.6, para cualquier arco desde (0,0) hasta (1,2):
∫γω=f(1,2)−f(0,0)=2+8+2=12.
Ejercicio 20.3★
Calcule ∬D(x+y)dxdy donde D es el dominio delimitado por y=x2 y y=x (0≤x≤1), en ambos órdenes de integración.
Usando coordenadas polares, calcule ∬D(1+x2+y2)2dxdy en todo el plano (como límite sobre los discos), y ∬D′xydxdy sobre el cuarto de disco D′={x,y≥0,x2+y2≤1}.
Solución
Solución de Ejercicio 20.4.
Primera integral. En el disco DR, en coordenadas polares:
(El astroide cabe en el disco unitario de áreaπ; tres octavos de π es plausible por su forma de estrella de cuatro cúspides).
Ejercicio 20.6★★
Calcular el volumen del sólido acotado por abajo por el paraboloide z=x2+y2y por arriba por el plano z=1, por ambos métodos: apilamiento (integre 1−x2−y2 sobre el disco unitario, polar coordenadas) y corte (los cortes horizontales son discos de radio z).
Solución
Solución de Ejercicio 20.6.
Apilado: encima de cada (x,y) del disco unitario D, z va desde x2+y2 a 1:
Rebanar: el corte en altura z∈[0,1] es el disco x2+y2≤z, de áreaπz:
V=∫01πzdz=2π.
Ejercicio 20.7★★
(Atracción gravitacional de una pelota — teorema de Newton, especial caso) Demuestre que el volumen de la cáscara esférica a≤r≤bes 34π(b3−a3)y calcule ∭Brdxdydz sobre la bola Bde radio R (r la distancia al origen). (Esférico coordenadas.)
Solución
Solución de Ejercicio 20.7.
En coordenadas esféricas el elemento de volumen es r2cosφdrdθdφ (Ejemplo 20.27), y la parte angular integra a 4π (2π de θ, ∫−π/2π/2cos=2). El volumen de la cáscara es
∫ab4πr2dr=34π(b3−a3).
Para la segunda integral, el integrando 1/r depende únicamente de r:
∭Brdxdydz=∫0R4πr2⋅r1dr=4π2R2=2πR2.
(El integrando explota en el origen, pero sin causar daño: r2/r=r es continuo — la integral sobre los shells ε≤r≤Rconverge como ε→0, que es la precisión sentido del enunciado. Este tipo de cálculo es el primero paso hacia el teorema de Newton de que una bola homogénea se atrae como una masa puntual en su centro.)
Ejercicio 20.8★★★
(Lema de Poincaré, caso en forma de estrella) Sea U en forma de estrella con con respecto a 0 (es decir, M∈U⇒[0,M]⊆U) y ω=Pdx+Qdy un formulario cerrado C1 en U. Definir
f(x,y)=∫01(xP(tx,ty)+yQ(tx,ty))dt.
Usando la diferenciación bajo el signo integral. (Capítulo 9) y Py=Qx, muestran que fx=P y fy=Q: cada forma cerrada en un conjunto abierto en forma de estrella es exacto.
Solución
Solución de Ejercicio 20.8.
El integrando g(t;x,y)=xP(tx,ty)+yQ(tx,ty) es C1 en (x,y), continuo en t, con parcial derivados continuo sobre [0,1]×U; diferenciación bajo el signo integral (Capítulo 9, aplicado en el compactot-intervalo [0,1], la dominación es automática allí) da
entonces el integrando es P(tx,ty)+tdtdP(tx,ty)=dtd[tP(tx,ty)] y
fx(x,y)=[tP(tx,ty)]01=P(x,y).
Simétricamente fy=Q (se utiliza el mismo cálculo con Py=Qx) al revés). Note donde entra la hipótesis: f está definido integrando a lo largo del segmento [0,M], que se encuentra en U Precisamente porque U tiene forma de estrella.
Ejercicio 20.9★★★
(Integral de Dirichlet por doble integración) Justificar y explotar
∫0∞∫0∞e−xysinxdydxvs∫0∞∫0∞e−xysinxdxdy
en [0,A]×[0,∞): muestra ∫0Axsinxdx=2π−∫0∞e−Ay1+y2ysinA+cosAdy y recupera ∫0∞xsinxdx=2π, comparando con la prueba integral de parámetros de Capítulo 9.
Solución
Solución de Ejercicio 20.9.
En la tira [0,A]×[0,∞) la función (x,y)↦e−xysinxno es absolutamente integrable hasta y=∞ uniformemente en un sentido ingenuo, pero cada integral iterada converge y su igualdad se sigue de Fubini en [0,A]×[0,B]más un límite B→∞ (la cola ∫0A∫B∞e−xy∣sinx∣dydx≤∫0Axe−Bx∣sinx∣dx≤∫0Ae−Bxdx→0, utilizando ∣sinx∣≤x).
y primero:∫0∞e−xydy=x1 para x>0, por lo que la primera integral es ∫0Axsinxdx.
x primero: dos integraciones por partes (o tomando la parte imaginaria de ∫0Ae(i−y)xdx) dar
∫0Ae−xysinxdx=1+y21−e−Ay(ysinA+cosA).
Integrando y sobre [0,∞), el término ∫0∞1+y2dy=2π se escinde:
∫0Axsinxdx=2π−∫0∞e−Ay1+y2ysinA+cosAdy.
El resto está limitado por ∫0∞e−Ay1+y2y+1dy≤∫0∞e−Ay⋅1+y21+ydy→0 como A→∞ (convergencia dominada, o la crudo atado 1+y21+y≤23 dando 2A3). Por lo tanto ∫0∞xsinxdx=2π — el mismo valor obtenido en Capítulo 9 diferenciando un integral de parámetros; aquí Fubini hace el trabajo.
Ejercicio 20.10★★★
(Desigualdad isoperimétrica vía Wirtinger) Sea γ un simple curva C1 cerrada de longitud2π, parametrizada por arco longitud en [0,2π], encerrando a áreaA. Usando Corolario 20.22, Parseval y la desigualdad de Wirtinger (ejercicios de Capítulo 14), demostrar A≤π, con igualdad para el círculo. (Normalize ∫02πx(s)ds=0; write 2A=∮xdy−ydxand bound 2A≤∫(x2+y′2)carefully via 2A=∫02π(xy′−yx′)dsand x2+y′2≥2xy′.)
Integrando ∮ydx por partes a lo largo del periodo (límite términos se cancelan por periodicidad), −∫yx′=∫y′x, por lo que de hecho 2A=2∫02πxy′ds. Entonces 2xy′≤x2+y′2 da
Desigualdad de Wirtinger (ejercicios Capítulo 14: para un 2π-función periódica C1 con media cero, ∫x2≤∫x′2) hace que la última integral no sea negativa: A≤π. La igualdad requiere igualdad en Wirtinger (x(s)=acoss+bsins) y en 2xy′≤x2+y′2(y′=xpuntualmente), lo que fuerza a y=asins−bcoss+c: la curva es el círculo unitario (adecuadamente centrado). Desde una curva de longitudL reescala a longitud2π, la declaración general es A≤4πL2: entre todas las curvas cerradas de un perímetro dado, la El círculo encierra el área más grande.
Ejercicio 20.11★★
(Momentos de la bola) Para la bola B de radio R en R3, calcular ∭Bz2dxdydz en coordenadas esféricas, y deducir ∭B(x2+y2+z2)dxdydz por simetría. Verifique el este último contra el cálculo del shell ∫0Rr2⋅4πr2dr.
Solución
Solución de Ejercicio 20.11.
En coordenadas esféricas z=rsinφ y dxdydz=r2cosφdrdθdφ:
Por la simetría de la pelota bajo coordenadas permutantes, ∭Bx2=∭By2=∭Bz2, entonces ∭B(x2+y2+z2)=3⋅154πR5=54πR5. Verificación de shell: ∫0Rr2⋅4πr2dr=54πR5— el integrando r2 es constante en la esfera de radio r, de área4πr2.
Ejercicio 20.12★★
(Fórmula del cordón) Sea K un polígono con vértices (x1,y1),…,(xm,ym) en orden antihorario (índices modelo m). Deduzca de Corolario 20.22 que
Area(K)=21i=1∑m(xiyi+1−xi+1yi),
y verifique la fórmula en el triángulo (0,0), (1,0), (0,1).
Solución
Solución de Ejercicio 20.12.
Parametrizar el borde de (xi,yi) a (xi+1,yi+1) por γ(t)=((1−t)xi+txi+1,(1−t)yi+tyi+1). Su contribución a 21∮(xdy−ydx) es
los términos cruzados se cancelan. Sumar los bordes m da la fórmula del cordón, de Corolario 20.22. Triangulo (0,0),(1,0),(0,1): 21((0⋅0−1⋅0)+(1⋅1−0⋅0)+(0⋅0−0⋅1))=21, el área correcto.
20.6 Problema: el volumen de la pelota en la dimensión n
Problema 20.1
Problema del fin de semana — Vn=πn/2/Γ(2n+1), y la extrañeza de la alta dimensiones
El disco tiene áreaπ, volumen de la bola34π — y entonces? Este problema calcula el volumen de la bola unitaria de Rn por cada n, dos veces (mediante una recursividad de corte impulsada por el Integrales de Wallis, luego a través de la función Γ y el Integral de Gauss de Ejemplo 20.17), y luego lee fuera de la geometría: los volúmenes alcanzan su punto máximo en la dimensión cinco y se precipitan a cero, y casi toda una bola de alta dimensión se esconde en un capa delgada cerca de su límite. Para una función continuo en una bola de Rn, la integral se entiende como la integral iterada n veces (cortando una coordenada a la vez, como en el capítulo n≤3); escribimos Bn(R)para la bola cerrada de radio R centrado en 0, vn(R) para su volumen, y Vn=vn(1), con V0=1 por convención.
Parte I — The slicing recursion.
Sustituyendo xi=Rui en cada uno de los n iterados integrales, muestre vn(R)=VnRn.
Cortando Bn(1) a lo largo de su última coordenada, muestra
Vn=Vn−1∫−11(1−t2)2n−1dt.
Con t=sinθ, identifique la integral como Integral de Wallis: ∫−11(1−t2)2n−1dt=2Wn, donde Wn=∫0π/2cosnθdθ=∫0π/2sinnθdθ.
Demostrar las dos identidades de Wallis (integrar por partes; luego telescopio nWnWn−1):
Wn=nn−1Wn−2(n≥2),WnWn−1=2nπ(n≥1).
Parte II — The recursion solved.
Combine las preguntas 2 a 4 en la recursividad de dos pasos
Tabule V1,…,V7 numéricamente. Usando el relación Vn/Vn−2=2π/n y los valores de 2W5 y 2W6, demuestre que la secuencia (Vn) aumenta hasta su máximo V5=158π2≈5.26 y disminuye a partir de ahí.
Demuestre que Vn→0 es más rápido que cualquier patrón geométrico. secuencia, y que ∑n≥1Vn converge: todos las bolas unitarias juntas tienen un volumen total finito.
Probar la identidad generadora
k≥0∑V2kx2k=eπx2(x∈R),
y deducir ∑k≥0V2k=eπ≈23.14.
Parte III — Second route: Γ and the Gauss integral.
Muestra por Fubini (el integrando es un producto) que
In=∫Rne−∥x∥2dx=(∫−∞+∞e−t2dt)n=πn/2,
la integral de Gauss n-dimensional, entendida como una límite sobre cubos [−R,R]n.
Recuperar Γ(s)=∫0∞ts−1e−tdt (Definición 9.17). De Γ(s+1)=sΓ(s) (Teorema 9.18) y Γ(21)=π (sustituto t=u2 e invocar la integral de Gauss), calcular
Γ(k+1)=k!,Γ(k+23)=2k+11⋅3⋯(2k+1)π.
Demostrar, por inducción mediante la recursividad de la pregunta 5, la fórmula única
Vn=Γ(2n+1)πn/2(n≥1),
y verifique que reproduzca tanto formas cerradas de la pregunta 6.
Muestra ∫0∞e−r2rn−1dr=21Γ(2n) y deduce el identidad
In=nVn∫0∞e−r2rn−1dr.
Interpretarlo: la masa gaussiana de Rn es recogidos a lo largo de conchas esféricas cuyo “(n−1)-dimensional área” en el radio r es nVnrn−1 — ambas partes ahora están probadas de forma independiente, por lo que la interpretación no cuesta nada.
Establezca sn−1=nVn (el área de la unidad de esfera Sn−1, consistente con vn(R)=∫0Rsn−1rn−1dr). Tabular s0,…,s3 y marque s1=2π, s2=4π, s3=2π2.
Parte IV — High dimensions are strange.
De la fórmula de Stirling (Teorema 6.13) aplicado a k!, mostrar incluso n=2k:
Vn∼πn1(n2πe)n/2(n→∞,n even),
y explicar por qué la misma desintegración supergeométrica se limita se extiende al impar n mediante la recursividad.
La bola unitaria se encuentra en el cubo [−1,1]n de volumen 2n. Calcule la relación de llenado Vn/2n para n=2,3,10, y muestra que tiende a 0: en alto dimensión, esencialmente todo el cubo se encuentra en su esquinas.
Demuestre que la fracción de vn(1) que se encuentra dentro la distancia ε de la esfera límite es 1−(1−ε)n→1; numéricamente, ¿qué fracción de una bola de dimensiones 100 se encuentra en el capa exterior de espesor 1%?
Demuestre el Wn∼2nπ(monotonicity of (Wn), the ratio Wn/Wn−2→1, and WnWn−1=2nπ) asintótico de Wallis y el límite inferior Wn≥2(n+1)π para todos los n.
(Concentración en una losa) La fracción de la unidad bola con la primera coordenada más allá de δ es ∫δ1(1−x2)2n−1dx/(2Wn). Usando 1−u≤e−u y la cola atado ∫δ∞e−ax2dx≤2aδe−aδ2, muestra esto la fracción es como máximo
(n−1)δe−(n−1)δ2/2/n+12π
y concluimos: para δ=s/n−1, todos menos un la fracción O(e−s2/2/s) de la pelota se encuentra en el losa ∣x1∣≤s/n−1. una bola de alto dimensión es, estadísticamente, un panqueque delgado en cada dirección a la vez.
Reúna las preguntas 16 a 19 en un párrafo: dónde el volumen de Bn(1) se encuentra (cerca del límite esfera, pero dentro de O(1/n) losas de cada hiperplano que pasa por el centro), y por qué estos dos las declaraciones no se contradicen entre sí.
Part V — Other bodies, and synthesis.
(Simplex) Deje Δn={x∈Rn:xi≥0,∑xi≤1}. Demuestre por corte e inducción que vol(Δn)=n!1.
(Politopo cruzado) Deduzca que Cn={x:∑∣xi∣≤1} tiene volumen n!2n, y verificar el sándwich Cn⊆Bn(1)⊆[−1,1]n a nivel de volúmenes: n!2n≤Vn≤2n.
Calcule V4 de una tercera manera: corte R4=R2×R2, integre el área del disco (z,w) el disco (x,y) en coordenadas polares y recuperar V4=2π2.
(Monte Carlo en problemas) Se sortea un punto uniformemente en el cubo [−1,1]20. Demuestre que el la probabilidad de que caiga en la bola inscrita es V20/220≈2.5⋅10−8, más o menos Se necesitan cuarenta millones de sorteos antes del primer golpe. se espera: estimación de Vn mediante muestreo de rechazo colapsa en alta dimensión (la maldición de dimensionalidad).
Síntesis. Dos derivaciones independientes se encontraron en Vn=πn/2/Γ(2n+1): enumera qué teorema de este capítulo cada uno utilizado (Fubini, cambio de variables, la integral polar de Gauss), y que entradas de una variable (Wallis, Γ, Stirling). ¿Dónde rehace este cálculo el volumen del año 3? con la teoría de Lebesgue, y ¿qué aporta?
Solución
Solución de Problema 20.1.
1. La pelota Bn(R) se describe mediante límites iterados. −R≤xn≤R, luego ∣xn−1∣≤R2−xn2, y así sucesivamente; sustituyendo xi=Rui en cada uno de los n integrales de una variable multiplica cada una por R y asigna los límites a los de Bn(1): vn(R)=Rnvn(1)=VnRn.
2. Cortando a lo largo de xn=t: la rebanada de Bn(1) es la pelota Bn−1(1−t2), entonces, por pregunta 1,
Vn=∫−11vn−1(1−t2)dt=Vn−1∫−11(1−t2)2n−1dt.
3. Con t=sinθ, dt=cosθdθy(1−t2)2n−1=cosn−1θen [−π/2,π/2]:
y θ↦2π−θ intercambia el seno y formas cosenos de Wn.
4. Escribe sinn=sinn−2(1−cos2) y integrar ∫sinn−2cos⋅cos por partes (v=n−1sinn−1):
Wn=Wn−2−n−1Wn⟹Wn=nn−1Wn−2.
De ahí nWnWn−1=(n−1)Wn−1Wn−2: la secuencia (nWnWn−1) es constante, igual a 1⋅W1W0=1⋅2π, por lo que WnWn−1=2nπ.
5. Las preguntas 2–3 dan Vn=2WnVn−1, dos veces:
Vn=2Wn⋅2Wn−1Vn−2=42nπVn−2=n2πVn−2.
6. De V0=1: V2k=2k2πV2k−2=kπV2k−2, entonces V2k=k!πk por inducción. De V1=2: V2k+1=2k+12πV2k−1, entonces
V2k+1=2j=1∏k2j+12π=1⋅3⋅5⋯(2k+1)2k+1πk.
7.V1=2, V2=π≈3.142, V3=34π≈4.189, V4=2π2≈4.935, V5=158π2≈5.264, V6=6π3≈5.168, V7=10516π3≈4.725. La relación de un paso es Vn/Vn−1=2Wn, y (Wn) está disminuyendo (sinn≤sinn−1puntualmente). Ahora 2W5=2⋅54⋅32=1516>1 mientras que 2W6=2⋅65⋅43⋅21⋅2π=165π<1: los ratios superan 1 hasta n=5 y están por debajo de 1 desde n=6 en — (Vn) aumenta a su máximo V5 y luego disminuye.
8. Para n≥13>4π: Vn/Vn−2=2π/n<21, entonces Vn≤C⋅2−n/2 con un fijo constante; mejor, para cualquier q>0, 2π/n<q2 para grandes n, entonces Vn/qn→0: la decadencia supera a toda geometría secuencia. La convergencia de ∑Vn se deriva de la relación Vn/Vn−2→0 (compárese con una serie geométrica de algún rango en).
9.∑k≥0V2kx2k=∑k≥0k!(πx2)k=eπx2, el serie exponencial (Capítulo 11), convergente para cada x. En x=1: ∑kV2k=eπ≈23.14.
10. En el cubo [−R,R]n el integrando es el producto ∏ie−xi2, por lo que la integral iterada factoriza: (∫−RRe−t2dt)n. Dejar R→∞ y usar ∫Re−t2dt=π (Ejemplo 20.17): In=πn/2.
11.t=u2 da Γ(21)=∫0∞t−1/2e−tdt=2∫0∞e−u2du=π. Iterando Γ(s+1)=sΓ(s): Γ(k+1)=k!Γ(1)=k!, y
12. Establecer Fn=πn/2/Γ(2n+1). Desde Γ(2n+1)=2nΓ(2n)=2nΓ(2n−2+1), obtenemos Fn=n2πFn−2: la misma recursividad que Vn (pregunta 5). Bases: F1=π/Γ(23)=π/(2π)=2=V1y F2=π/Γ(2)=π=V2. Por inducción Vn=Fn para todos n; La pregunta 11 convierte esto nuevamente en los dos formas cerradas de pregunta 6.
13. Con r=t: ∫0∞e−r2rn−1dr=21∫0∞t2n−1e−tdt=21Γ(2n). Por lo tanto
Una vez probadas ambas partes, la identidad puede ser leer como la descomposición de capas de la integral gaussiana: la esfera de radio r lleva áreanVnrn−1, y el gaussiano El peso e−r2 está integrado sobre las carcasas.
14.s0=V1=2 (la esfera 0 tiene dos puntos), s1=2V2=2π, s2=3V3=4π, s3=4V4=2π2; y ∫0Rsn−1rn−1dr=VnRn=vn(R): área es la derivada radial del volumen.
15. Para n=2k, Stirling (Teorema 6.13) da k!∼2πk(k/e)k, entonces
V2k=k!πk∼2πk(πe/k)k=πn1(n2πe)n/2(n=2k).
Para n impar: V2k+1=2W2k+1V2k≤2V2k, entonces se mantienen los mismos límites de desintegración supergeométrica (hasta un factor 2 y un desplazamiento de uno en el exponente) — por cada q>0, Vn=o(qn).
16.V2/4=π/4≈0.785; V3/8=π/6≈0.524; V10/210=120⋅1024π5≈0.0025. En general Vn−2/2n−2Vn/2n=4n2π=2nπ→0: la relación tiende a 0 (supergeométricamente). La bola inscrita ocupa un fracción evanescente: el volumen del cubo migra a su esquinas.
17. Por la pregunta 1 la bola interior de radio 1−εtiene volumen Vn(1−ε)n, por lo que la exterior shell lleva la fracción 1−(1−ε)n→1. Para n=100, ε=0.01: (0.99)100=e100ln0.99≈e−1.005≈0.366: aproximadamente 63% de la pelota se encuentra dentro de 1% de su superficie.
18.(Wn) disminuye, por lo que Wn≤Wn−1≤Wn−2=n−1nWn: apretando, Wn−1/Wn→1. Multiplicando por WnWn−1=2nπ: Wn2∼2nπ, es decir Wn∼π/(2n). Límite inferior: Wn2≥WnWn+1=2(n+1)π, entonces Wn≥π/(2(n+1)) por cada n.
19. Numerador: 1−x2≤e−x2 da (1−x2)2n−1≤e−(n−1)x2/2, y con a=2n−1,
Denominador: 2Wn≥2π/(n+1) por la pregunta 18. Al dividir se obtiene el límite mostrado. Para δ=s/n−1 se convierte en 2π(n−1)n+1e−s2/2/s=O(e−s2/2/s), uniformemente en n: fuera de la losa ∣x1∣≤s/n−1 casi no hay volumen, por s moderadamente grande — y por simetría lo mismo vale para cada dirección.
20. Las dos declaraciones coexisten porque describen diferentes coordenadas de un mismo punto. Casi todos los puntos de Bn(1) tiene norma cerca de 1 (pregunta 17: radial concentración cerca de la esfera), sin embargo, cada uno de sus n las coordenadas son pequeñas, de orden 1/n (pregunta 19), lo cual es consistente desde las coordenadas n del tamaño 1/n tiene norma del pedido 1. Concentrados de volumen de alta dimensión donde todas las coordenadas comparten el presupuesto norma por igual — cerca de la esfera, pero lejos de cada polo del eje de coordenadas.
21. Corte Δn en xn=t∈[0,1]: el la porción es {x′∈Rn−1:xi≥0,∑xi≤1−t}=(1−t)Δn−1, de volumen (1−t)n−1vol(Δn−1) por homogeneidad. entonces
22. Los orantes de signo 2n cortan Cn en 2n copias de Δn (con superposiciones insignificantes en las coordenadas hiperplanos): vol(Cn)=n!2n. si ∑∣xi∣≤1 luego ∑xi2≤(∑∣xi∣)2≤1: Cn⊆Bn(1); y Bn(1)⊆[−1,1]n desde ∣xi∣≤∥x∥. Por lo tanto n!2n≤Vn≤2n — consistente con la pregunta 15, que coloca a Vn entre los Escalas factoriales y geométricas.
23. Para (x,y) en el disco unitario, la porción de B4(1) es el disco de radio 1−x2−y2 en el (z,w)-avión, de áreaπ(1−x2−y2). en polares coordenadas:
24. La probabilidad es la relación de volumen. 220V20=10!⋅220π10≈10485760.0258≈2.5⋅10−8. El número de sorteos hasta el primer acierto es de orden inverso, sobre 4⋅107: un muestreador de rechazo que funcionó maravillosamente para el disco (π/4 de visitas) es inútil en la dimensión 20 — la maldición de la dimensionalidad en una línea.
25. Ruta uno (Partes I–II) utilizada: Tipo Fubini corte de la integral iterada, sustitución de una variable en cada coordenada (homogeneidad), y el Integrales de Wallis — cálculo puro de una variable más inducción. Ruta dos (Parte III) utilizado: Fubini para la estructura del producto In, el cambio polar de variables a través de la integral de Gauss de Ejemplo 20.17 y la función Γ. ecuación funcional. Se encuentran en Vn=πn/2/Γ(2n+1), con Stirling (Teorema 6.13) convirtiendo la fórmula en asintóticos. El volumen del año 3 reconstruye todo esto en La integral de Lebesgue: ahí Fubini y el cambio de variables son teoremas para funciones integrables generales, esféricas Las coordenadas existen en todas las dimensiones, y los mismos pelota volúmenes reaparecen como dividendos trabajados del medida del producto y problemas de Stirling — con dominado convergencia reemplazando nuestros apretón hechos a mano.