Mathematics · Book 4 · Bachelor Year 2

Matemáticas universitarias — Grado 2

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 diferencial P ⁣dx+Q ⁣dyP\,\dd x + Q\,\dd y 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 UR2U \subseteq \R^2 abierto. A diferencial forma de grado 11 y clase C0\mathcal{C}^0 en UU es una expresión ω=P ⁣dx+Q ⁣dy\omega = P\,\dd x + Q\,\dd y con P,Q ⁣:URP, Q \colon U \to \R continuo — formalmente, un mapa continuo de UU al dual de R2\R^2, ω(M)=P(M)e1+Q(M)e2\omega(M) = P(M)\,e_1^* + Q(M)\,e_2^*. por un C1\mathcal{C}^1 arco γ ⁣:[a,b]U\gamma \colon [a, b] \to U, γ(t)=(x(t),y(t))\gamma(t) = (x(t), y(t)), el integral de línea de ω\omega junto γ\gamma es

γω=ab(P(γ(t))x(t)+Q(γ(t))y(t)) ⁣dt.\int_\gamma \omega = \int_a^b \Bigl(P(\gamma(t))\,x'(t) + Q(\gamma(t))\,y'(t)\Bigr)\,\dd t .

Las definiciones se extienden textualmente a R3\R^3 (formas P ⁣dx+Q ⁣dy+R ⁣dzP\,\dd x + Q\,\dd y + R\,\dd z) y a los arcos C1\mathcal{C}^1 por partes (suma sobre las piezas).

Proposición 20.2 (Invariancia y orientación)

El integral de línea no cambia bajo un C1\mathcal{C}^1 creciente. cambio de parametro, y cambia de signo a uno decreciente. eso por lo tanto solo depende del orientado arco geométrico.

Demostración. Si θ ⁣:[c,d][a,b]\theta \colon [c, d] \to [a, b] es un cambio de parametro y γ~=γθ\tilde\gamma = \gamma \circ \theta, luego por la regla de la cadena y el cambio de una variable de variables t=θ(u)t = \theta(u),

γ~ω=cd(P(γ(θ(u)))x(θ(u))+Q(γ(θ(u)))y(θ(u)))θ(u) ⁣du=±ab(Px+Qy)(t) ⁣dt,\int_{\tilde\gamma}\omega = \int_c^d \bigl(P(\gamma(\theta(u)))\,x'(\theta(u)) + Q(\gamma(\theta(u)))\,y'(\theta(u))\bigr)\,\theta'(u)\,\dd u = \pm\int_a^b \bigl(Px' + Qy'\bigr)(t)\,\dd t ,

con signo ++ si θ\theta es creciente (θ(c)=a\theta(c) = a) y - si disminuye (los límites se intercambian).

Ejemplo 20.3 (Trabajo de una fuerza; circulación)

Si F=(P,Q)F = (P, Q) es un campo de fuerza, γP ⁣dx+Q ⁣dy=abF(γ(t)),γ(t) ⁣dt\int_\gamma P\dd x + Q\dd y = \int_a^b \langle F(\gamma(t)), \gamma'(t)\rangle\,\dd t es el trabajar de FF junto con γ\gamma. Para ω=y ⁣dx+x ⁣dy\omega = -y\,\dd x + x\,\dd ya lo largo del círculo unitario γ(t)=(cost,sint)\gamma(t) = (\cos t, \sin t) en sentido antihorario:

γω=02π((sint)(sint)+costcost) ⁣dt=2π,\int_\gamma \omega = \int_0^{2\pi}\bigl((-\sin t)(-\sin t) + \cos t\cos t\bigr)\,\dd t = 2\pi ,

dos veces el área adjunto — un primer indicio de Green–Riemann.

Ejemplo 20.4 (Una integral, dos parametrizaciones, una trampa de señales)

Calcule γx ⁣dy\int_\gamma x\,\dd y a lo largo de la unidad superior semicírculo de (1,0)(1, 0) a (1,0)(-1, 0). Con γ(t)=(cost,sint)\gamma(t) = (\cos t, \sin t), t[0,π]t \in \intcc0\pi:

0πcostcost ⁣dt=π2.\int_0^\pi\cos t\cdot\cos t\,\dd t = \frac\pi2 .

Con la parametrización del gráfico x(x,1x2)x \mapsto (x, \sqrt{1 - x^2}), xxde 11a1-1 (¡tenga en cuenta la dirección!):

11xx1x2 ⁣dx=11x21x2 ⁣dx=π2\int_1^{-1}x\cdot\frac{-x}{\sqrt{1 - x^2}}\,\dd x = \int_{-1}^{1}\frac{x^2}{\sqrt{1 - x^2}}\,\dd x = \frac\pi2

(x=sinux = \sin u 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)(1,0) a (1,0)(-1,0); revertir el viajar invierte el cartel. Cerrando el camino a lo largo del eje xx (donde  ⁣dy=0\dd y = 0) no agrega nada, y el total π2\frac\pi2 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 ω=P ⁣dx+Q ⁣dy\omega = P\,\dd x + Q\,\dd y de la clase C0\mathcal{C}^0 es exacto en UU si hay fC1(U)f \in \mathcal{C}^1(U) (a potencial) con ω= ⁣df\omega = \dd f, es decir P=fxP = f_x y Q=fyQ = f_y. Un formulario C1\mathcal{C}^1 es cerrado si Py=QxP_y = Q_x está en UU.

Teorema 20.6 (Teorema fundamental de la recta integrales)

Si ω= ⁣df\omega = \dd f es exacto y γ\gamma es por partes C1\mathcal{C}^1 arco en UU de AA a BB, luego

γω=f(B)f(A).\int_\gamma \omega = f(B) - f(A) .

En particular, la integral de un forma exacta a lo largo de cualquier arco cerrado es cero y cada formulario exacto C1\mathcal{C}^1 está cerrado.

Demostración.  ⁣d ⁣dtf(γ(t))=fx(γ(t))x(t)+fy(γ(t))y(t)\frac{\dd}{\dd t}f(\gamma(t)) = f_x(\gamma(t))x'(t) + f_y(\gamma(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 tf(γ(t))t \mapsto f(\gamma(t)), y el teorema fundamental del cálculo da la resultado en cada pieza; El telescopio de valores intermedios. cerrazón de las formas exacto C1\mathcal{C}^1 es el teorema de Schwarz: Py=fxy=fyx=QxP_y = f_{xy} = f_{yx} = Q_x.

Ejemplo 20.7 (Reconstruyendo un potencial)

Deja ω=yexy ⁣dx+(xexy+2y) ⁣dy\omega = y\,\eu^{xy}\,\dd x + (x\,\eu^{xy} + 2y)\,\dd y en R2\R^2. Es cerrado: ambas derivadas cruzadas son iguales. exy(1+xy)\eu^{xy}(1 + xy). Para encontrar un potencial, integre PP en xx en yy fijo:

f(x,y)=yexy ⁣dx=exy+c(y),f(x, y) = \int y\,\eu^{xy}\,\dd x = \eu^{xy} + c(y),

luego ajuste cc haciendo coincidir fyf_y: xexy+c(y)=xexy+2yx\,\eu^{xy} + c'(y) = x\,\eu^{xy} + 2yda c(y)=y2c(y) = y^2. Entonces f(x,y)=exy+y2f(x,y) = \eu^{xy} + y^2, y para cualquier arco C1\mathcal C^1 por partes de (0,0)(0,0) a (1,1)(1,1),

γω=f(1,1)f(0,0)=(e+1)1=e,\int_\gamma\omega = f(1,1) - f(0,0) = (\eu + 1) - 1 = \eu ,

independientemente de la ruta — la receta de dos pasos (integrar en xx, correcto en yy) es lo contrario práctico de Teorema 20.6 en dominios donde está formas cerradas exacto.

Ejemplo 20.8 (Cerrado no implica exacto)

En U=R2{0}U = \R^2 \setminus \{0\}, el forma de ángulo

ω=y ⁣dx+x ⁣dyx2+y2\omega = \frac{-y\,\dd x + x\,\dd y}{x^2 + y^2}

está cerrado (cálculo directo: tanto PyP_y como QxQ_x son iguales y2x2(x2+y2)2\frac{y^2 - x^2}{(x^2+y^2)^2}), pero es integral a lo largo de la unidad. el círculo es 2π02\pi \neq 0 (el mismo cálculo que Ejemplo 20.3, dividido por 11): ω\omega no es exacto en UU. Localmente, ω= ⁣dθ\omega = \dd\theta para una determinación. θ\theta 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\mathcal{C}^1 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 ff continuo en un rectángulo R=[a,b]×[c,d]R = [a, b] \times [c, d] tiene una integral doble Rf\iint_R f, definido por sumas de Riemann sobre cuadrículas exactamente como en una variable, y calculado por iteración:

Teorema 20.9 (Fubini sobre un rectángulo)

Para ff continuo en R=[a,b]×[c,d]R = [a,b] \times [c,d],

Rf=ab(cdf(x,y) ⁣dy) ⁣dx=cd(abf(x,y) ⁣dx) ⁣dy.\iint_R f = \int_a^b \Bigl(\int_c^d f(x, y)\,\dd y\Bigr)\dd x = \int_c^d \Bigl(\int_a^b f(x, y)\,\dd x\Bigr)\dd y .

Demostración. Establezca F(x)=cdf(x,y) ⁣dyF(x) = \int_c^d f(x, y)\,\dd y. Uniforme continuidad de ff en el compacto RR hace FF continuo (estimación dominada: F(x)F(x)(dc)supyf(x,y)f(x,y)\abs{F(x) - F(x')} \leq (d - c)\sup_y\abs{f(x,y) - f(x',y)}). Ahora subdivida [a,b][a,b] y [c,d][c,d] en nn partes iguales, dando un cuadrícula de celdas RijR_{ij} de área ΔxΔy\Delta x\,\Delta y. En cada celda, infRijfΔxΔyxi1xiyj1yjf(x,y) ⁣dy ⁣dxsupRijfΔxΔy\inf_{R_{ij}} f \cdot \Delta x \Delta y \leq \int_{x_{i-1}}^{x_i}\int_{y_{j-1}}^{y_j} f(x,y)\,\dd y\,\dd x \leq \sup_{R_{ij}} f \cdot \Delta x \Delta y por monotonicidad del Integral de una variable (aplicada dos veces). Sumando las celdas, el La integral iterada abF\int_a^b F 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\iint_R f como nn \to \infty. El mismo argumento se aplica con los roles de xxy yy. intercambiadas, por lo que ambas integrales iteradas son iguales a Rf\iint_R f.

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 DR2D \subseteq \R^2 es yy-elemental si

D={(x,y):axb, φ1(x)yφ2(x)}D = \{(x, y) : a \leq x \leq b,\ \varphi_1(x) \leq y \leq \varphi_2(x)\}

con φ1φ2\varphi_1 \leq \varphi_2 continuo en [a,b][a,b] (xx-elemental: simétricamente). Para ff continuo en un yy-elemental DD,

Df=ab(φ1(x)φ2(x)f(x,y) ⁣dy) ⁣dx,\iint_D f = \int_a^b\Bigl( \int_{\varphi_1(x)}^{\varphi_2(x)} f(x,y)\,\dd y\Bigr)\dd x ,

y uno verifica (extendiendo ff mediante un argumento de aproximación, o subdividiendo) que cuando DD 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 DD es Area(D)=D1\operatorname {Area}(D) = \iint_D 1.

Ejemplo 20.12

En el triángulo D={0x1, 0yx}D = \{0 \leq x \leq 1,\ 0 \leq y \leq x\}:

Dxy ⁣dx ⁣dy=01x(0xy ⁣dy) ⁣dx=01xx22 ⁣dx=18.\iint_D xy \,\dd x\,\dd y = \int_0^1 x\Bigl(\int_0^x y\,\dd y\Bigr)\dd x = \int_0^1 x\cdot\frac{x^2}{2}\,\dd x = \frac18 .

Intercambiando el orden (xx de yy a 11): 01y(y1x ⁣dx) ⁣dy=01y1y22 ⁣dy=18\int_0^1 y\bigl(\int_y^1 x\,\dd x\bigr)\dd y = \int_0^1 y\,\frac{1 - y^2}{2}\,\dd y = \frac18 — 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 ⁣ ⁣x1ey2 ⁣dy ⁣dxI = \displaystyle\int_0^1\!\!\int_x^1 \eu^{y^2}\,\dd y\,\dd x. Como está escrito, la integral interna ey2 ⁣dy\int\eu^{y^2}\dd y no tiene antiderivada elemental: la el cálculo está atascado. Pero el dominio es el triángulo 0xy10 \leq x \leq y \leq 1, que es elemental en ambas direcciones; cambiando el orden,

I=01 ⁣ ⁣0yey2 ⁣dx ⁣dy=01yey2 ⁣dy=[12ey2]01=e12.I = \int_0^1\!\!\int_0^y \eu^{y^2}\,\dd x\,\dd y = \int_0^1 y\,\eu^{y^2}\,\dd y = \Bigl[\tfrac12\eu^{y^2}\Bigr]_0^1 = \frac{\eu - 1}{2} .

La variable interna xx no apareció en ninguna parte del integrando, por lo que integrarlo primero produjo exactamente el factor yy 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 Φ ⁣:UU\Phi \colon U' \to U un difeomorfismo C1\mathcal{C}^1 entre conjuntos abierto de R2\R^2, sea KUK \subseteq U un compacto dominio cortado en pedazos elementales con K=Φ1(K)K' = \Phi^{-1}(K), y sea ff continuo en KK. entonces

Kf(x,y) ⁣dx ⁣dy=Kf(Φ(u,v))detJΦ(u,v) ⁣du ⁣dv.\iint_K f(x, y)\,\dd x\,\dd y = \iint_{K'} f\bigl(\Phi(u, v)\bigr)\, \abs{\det J_\Phi(u, v)}\,\dd u\,\dd v .

Demostración. Admitido a este nivel.

Observación 20.15

La prueba completo — se aproxima a Φ\Phi 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\dd u en (u,v)(u, v) se asigna, en primer orden, al paralelogramo abarcado por Φu ⁣du\Phi_u\,\dd u y Φv ⁣dv\Phi_v\,\dd v, cuyo área es detJΦ ⁣du ⁣dv\abs{\det J_\Phi}\,\dd u\,\dd v (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+y2x^2 + y^2 pide polar, una estructura de producto Requiere mantener los ejes cartesianos. Forma de la límite: límites u(x,y)=c1u(x,y) = c_1, v(x,y)=c2v(x,y) = c_2 ruegan por las coordenadas (u,v)(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+yx + y y xyx - y invitan al grado 4545. 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α)\Phi(\rho, \alpha) = (\rho\cos\alpha,\ \rho\sin\alpha) tiene

JΦ=(cosαρsinαsinαρcosα),detJΦ=ρ,J_\Phi = \begin{pmatrix} \cos\alpha & -\rho\sin\alpha\\ \sin\alpha & \rho\cos\alpha \end{pmatrix}, \qquad \det J_\Phi = \rho ,

entonces  ⁣dx ⁣dy=ρ ⁣dρ ⁣dα\dd x\,\dd y = \rho\,\dd\rho\,\dd\alpha. Para el disco DRD_R de radio RR:

DRe(x2+y2) ⁣dx ⁣dy=02π ⁣ ⁣0Reρ2ρ ⁣dρ ⁣dα=π(1eR2)Rπ.\iint_{D_R} e^{-(x^2 + y^2)}\,\dd x\,\dd y = \int_0^{2\pi}\!\!\int_0^R e^{-\rho^2}\rho\,\dd\rho\,\dd\alpha = \pi\bigl(1 - e^{-R^2}\bigr) \xrightarrow[R\to\infty]{} \pi .

Comparando con el cuadrado [R,R]2[-R, R]^2 (que está apretado entre los discos DRD_R y DR2D_{R\sqrt2}, todos los integrandos positivos) da (ex2 ⁣dx)2=π\bigl(\int_{-\infty}^\infty e^{-x^2}\dd x\bigr)^2 = \pi:

 +ex2 ⁣dx=π \boxed{\ \int_{-\infty}^{+\infty} e^{-x^2}\,\dd x = \sqrt{\pi}\ }

— 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\Phi(u, v) = M(u, v)^{\mathsf T} + C con MM invertible, el jacobiano es la matriz constante MM: áreas se multiplican por el factor constante detM\abs{\det M} — la promesa hecha en Capítulo 17 es ahora un teorema. Dos usos inmediatos. La elipse x2a2+y2b21\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1 es la imagen del disco unitario bajo (u,v)(au,bv)(u, v) \mapsto (au, bv), por lo que área es abπab \cdot \pi — sin cálculo. Y para la integral de f(x+y)f(x + y) sobre el cuadrado K=[0,1]2K = \intcc01^2, la cizalla Φ(u,v)=(uv,v)\Phi(u, v) = (u - v, v) (determinante 11) lo convierte en una integral de f(u)f(u) sobre un paralelogramo, que Fubini corta a uu constante: con f=expf = \exp,

Kex+y ⁣dx ⁣dy=(01ex ⁣dx)2=(e1)2,\iint_K \eu^{x+y}\,\dd x\,\dd y = \Bigl(\int_0^1 \eu^x\,\dd x\Bigr)^2 = (\eu - 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 DD la región del primer cuadrante acotada por el hipérbolas xy=1xy = 1 y xy=3xy = 3 y las líneas y=xy = x y y=3xy = 3x. En las coordenadas u=xyu = xy, v=y/xv = y/x el dominio se convierte en el cuadrado [1,3]×[1,3]\intcc13 \times \intcc13; invertir,

x=u/v,y=uv,detJ=xuyvxvyu=12vx = \sqrt{u/v}, \qquad y = \sqrt{uv}, \qquad \det J = x_uy_v - x_vy_u = \frac{1}{2v}

(un cálculo de dos líneas con x=u1/2v1/2x = u^{1/2}v^{-1/2}, y=u1/2v1/2y = u^{1/2}v^{1/2}). Por lo tanto

Area(D)=13 ⁣ ⁣13 ⁣du ⁣dv2v=2ln32=ln31.10.\operatorname{Area}(D) = \int_1^3\!\!\int_1^3\frac{\dd u\,\dd v}{2v} = 2\cdot\frac{\ln 3}{2} = \ln 3 \approx 1.10 .

Intentar cortar DD 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 ρ\rho para polar coordenadas.

Ejemplo 20.20 (Valores medios)

El valor medio de ff sobre un dominio DD es 1Area(D)Df\frac1{\operatorname{Area}(D)}\iint_Df. Muestra: el promedio distancia desde el centro para un punto uniformemente elegido de la disco de radio RR es

1πR202π ⁣ ⁣0Rρρ ⁣dρ ⁣dα=2πR3/3πR2=2R3,\frac{1}{\pi R^2}\int_0^{2\pi}\!\!\int_0^R \rho\cdot\rho\,\dd\rho\,\dd\alpha = \frac{2\pi R^3/3}{\pi R^2} = \frac{2R}3 ,

no R/2R/2: área uniforme pone más masa en radios grandes (el el anillo en el radio ρ\rho tiene un peso proporcional a ρ\rho), 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\bar z = 3R/8 del media bola se calcula más adelante en el capítulo en lugar de R/2R/2.

20.3 El teorema de Green-Riemann

Teorema 20.21 (Verde–Riemann)

Sea KR2K \subseteq \R^2 un dominio compacto que es elemental en ambas direcciones (o una unión finita de tales segmentos pegados), con límite K\partial K a trozos C1\mathcal{C}^1 cerrado orientado a curvas sinistrórsum (el dominio permanece en la izquierda). Para P,QP, Q de clase C1\mathcal{C}^1 en una vecindad de KK:

KP ⁣dx+Q ⁣dy=K(QxPy) ⁣dx ⁣dy.\oint_{\partial K} P\,\dd x + Q\,\dd y = \iint_K \Bigl(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\Bigr)\,\dd x\,\dd y .

Demostración. Primero, ambos lados son aditivos bajo el corte KK a lo largo de un segmentar en dos piezas K1,K2K_1, K_2: las integrales dobles suman por aditividad de \iint; en cuanto a las integrales de frontera, la límites en sentido antihorario de K1K_1 y K2K_2 cada recorrido el interior cortado una vez, en direcciones opuesto, por lo que en la suma

K1+K2=K+(cut, both ways)=K,\oint_{\partial K_1} + \oint_{\partial K_2} = \oint_{\partial K} + (\text{cut, both ways}) = \oint_{\partial 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 P ⁣dx=KPy\oint P\,\dd x = -\iint_K P_y en un yy-dominio elemental D={axb, φ1(x)yφ2(x)}D = \{a \leq x \leq b,\ \varphi_1(x) \leq y \leq \varphi_2(x)\}; la identidad Q ⁣dy=KQx\oint Q\,\dd y = \iint_K Q_x es simétrico (xx-elemental), y el teorema es su suma.

Calcular la integral doble de Fubini y la de una variable teorema fundamental:

DPy ⁣dx ⁣dy=ab(P(x,φ2(x))P(x,φ1(x))) ⁣dx.\iint_D \frac{\partial P}{\partial y}\,\dd x\,\dd y = \int_a^b \bigl(P(x, \varphi_2(x)) - P(x, \varphi_1(x))\bigr) \,\dd x .

Ahora el límite de DD, en sentido antihorario, consta de: la parte inferior gráfico y=φ1(x)y = \varphi_1(x) recorrido de izquierda a derecha, la derecha segmento vertical x=bx = b (hacia arriba), el gráfico superior y=φ2(x)y = \varphi_2(x) atravesó derecha a izquierda, el vertical izquierdo segmento x=ax = a (hacia abajo). A lo largo de los segmentos verticales xx está constante, por lo que contribuyen 00 a P ⁣dx\oint P\,\dd x; los gráficos, parametrizado por xx, dé

DP ⁣dx=abP(x,φ1(x)) ⁣dxabP(x,φ2(x)) ⁣dx=DPy ⁣dx ⁣dy.\oint_{\partial D} P\,\dd x = \int_a^b P(x, \varphi_1(x))\,\dd x - \int_a^b P(x, \varphi_2(x))\,\dd x = -\iint_D \frac{\partial P}{\partial y}\,\dd x\,\dd y . \qedhere

Corolario 20.22 (Área por el límite)

Bajo las hipótesis de Teorema 20.21,

Area(K)=Kx ⁣dy=Ky ⁣dx=12Kx ⁣dyy ⁣dx.\operatorname{Area}(K) = \oint_{\partial K} x\,\dd y = -\oint_{\partial K} y\,\dd x = \frac12\oint_{\partial K} x\,\dd y - y\,\dd x .

Demostración. Aplique Verde–Riemann a (P,Q)=(0,x)(P, Q) = (0, x), (y,0)(-y, 0) y 12(y,x)\frac12(-y, x): cada vez QxPy=1Q_x - P_y = 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 x ⁣dy\oint x\,\dd y cuando la parametrización haga  ⁣dy\dd y simple (gráficos sobre el eje yy), y ⁣dx-\oint y\,\dd x simétricamente, y la media suma simétrico cuando el la parametrización trata a xx y yy 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.

Ejemplo 20.24 (Área de la elipse)

Para x=acostx = a\cos t, y=bsinty = b\sin t, t[0,2π]t \in [0, 2\pi]:

Area=1202π(acostbcostbsint(asint)) ⁣dt=ab202π ⁣dt=πab.\operatorname{Area} = \frac12\int_0^{2\pi}\bigl(a\cos t \cdot b\cos t - b\sin t\cdot(-a\sin t)\bigr)\,\dd t = \frac{ab}{2}\int_0^{2\pi}\dd t = \pi ab .

Ejemplo 20.25 (Verde–Riemann como verificación cruzada)

Tome P=y3P = -y^3, Q=x3Q = x^3 en el disco unitario cerrado DD. Lado límite, con γ(t)=(cost,sint)\gamma(t) = (\cos t, \sin t):

DP ⁣dx+Q ⁣dy=02π(sin4t+cos4t) ⁣dt=2π(38+38)=3π2,\oint_{\partial D}P\,\dd x + Q\,\dd y = \int_0^{2\pi}\bigl(\sin^4 t + \cos^4 t\bigr)\dd t = 2\pi\cdot\Bigl(\frac38 + \frac38\Bigr) = \frac{3\pi}2 ,

por linealización (sin4+cos4=34+14cos4t\sin^4 + \cos^4 = \tfrac34 + \tfrac14\cos4t). Lado interior:

D(QxPy) ⁣dx ⁣dy=D3(x2+y2) ⁣dx ⁣dy=302π ⁣ ⁣01ρ3 ⁣dρ ⁣dα=3π2.\iint_D(Q_x - P_y)\,\dd x\,\dd y = \iint_D 3(x^2 + y^2)\,\dd x\,\dd y = 3\int_0^{2\pi}\!\!\int_0^1\rho^3\,\dd\rho\,\dd\alpha = \frac{3\pi}2 .

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=PyQ_x = P_y), la integral alrededor del límite de cualquier El dominio contenido en UU 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 UUintegrales 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 \iiint a tres integrales de una variable (ya sea por rebanar: Kf=(Kzf) ⁣dz\iiint_K f = \int\bigl(\iint_{K_z} f\bigr)\dd zsobre los cortes horizontales KzK_z, o por apilado: integrando en zz primero a lo largo de palos verticales), y la formula de cambio de variables se cumple con el 3×33 \times 3 Jacobiano.

Ejemplo 20.27 (Coordenadas cilíndricas y esféricas)

Cilíndrico (x,y,z)=(ρcosα,ρsinα,z)(x, y, z) = (\rho\cos\alpha, \rho\sin\alpha, z):  ⁣dx ⁣dy ⁣dz=ρ ⁣dρ ⁣dα ⁣dz\dd x\,\dd y\,\dd z = \rho\,\dd\rho\,\dd\alpha\,\dd z. Esférico (x,y,z)=(rcosθcosφ, rsinθcosφ, rsinφ)(x, y, z) = (r\cos\theta\cos\varphi,\ r\sin\theta\cos\varphi,\ r\sin\varphi)(θ\theta longitud, Latitud φ[π2,π2]\varphi \in [-\frac\pi2, \frac\pi2]): expandiendo 3×33 \times 3 determinante a lo largo de la última fila,

detJ=r2cosφ, ⁣dx ⁣dy ⁣dz=r2cosφ   ⁣dr ⁣dθ ⁣dφ.\det J = r^2\cos\varphi , \qquad \dd x\,\dd y\,\dd z = r^2\cos\varphi\;\dd r\,\dd\theta\,\dd\varphi .

Volumen de la bola de radio RR:

V=0R ⁣ ⁣02π ⁣ ⁣π/2π/2r2cosφ   ⁣dφ ⁣dθ ⁣dr=R332π2=43πR3,V = \int_0^R\!\!\int_0^{2\pi}\!\!\int_{-\pi/2}^{\pi/2} r^2\cos\varphi\;\dd\varphi\,\dd\theta\,\dd r = \frac{R^3}{3}\cdot 2\pi \cdot 2 = \boxed{\frac43\pi R^3} ,

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,z0, x+y+z1}T = \{x, y, z \geq 0,\ x + y + z \leq 1\}, por apilamiento: para (x,y)(x, y) fijo en el triángulo x+y1x + y \leq 1, zz pasa por encima de [0,1xy]\intcc0{1 - x - y}, por lo que

V=01 ⁣ ⁣01x(1xy) ⁣dy ⁣dx=01(1x)22 ⁣dx=16.V = \int_0^1\!\!\int_0^{1-x}(1 - x - y)\,\dd y\,\dd x = \int_0^1\frac{(1 - x)^2}{2}\,\dd x = \frac16 .

Por corte: la sección a la altura zz es el triángulo {x,y0, x+y1z}\{x, y \geq 0,\ x + y \leq 1 - z\}, de área(1z)22\frac{(1-z)^2}2, y V=01(1z)22 ⁣dz=16V = \int_0^1\frac{(1-z)^2}2\,\dd z = \frac16 nuevamente — el dos cálculos son las mismas integrales en una forma diferente orden, que es todo lo que afirma Fubini. El valor 16=13121\frac16 = \frac13\cdot\frac12\cdot1 es la fórmula del cono. (Ejemplo 20.30) con base triangular, y el nn versión tridimensional 1/n!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 HH de radio RR (z0z \geq 0), la La altura del centroide es zˉ=1VHz\bar z = \frac1{V}\iiint_H z, con V=23πR3V = \frac23\pi R^3. En coordenadas esféricas (z=rsinφz = r\sin\varphi, φ[0,π/2]\varphi \in \intcc0{\pi/2}):

Hz=0Rr3 ⁣dr02π ⁣dθ0π/2sinφcosφ ⁣dφ=R442π12=πR44,\iiint_H z = \int_0^R r^3\,\dd r\int_0^{2\pi}\dd\theta \int_0^{\pi/2}\sin\varphi\cos\varphi\,\dd\varphi = \frac{R^4}4\cdot2\pi\cdot\frac12 = \frac{\pi R^4}4 ,

entonces

zˉ=πR4/42πR3/3=3R8:\bar z = \frac{\pi R^4/4}{2\pi R^3/3} = \frac{3R}8 :

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/2R/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 área AA y altura hh (ápice hacia arriba, base en z=0z = 0): el corte en altura zz es la base escalada por el factor (1z/h)(1 - z/h), de áreaA(1z/h)2A(1 - z/h)^2. Por lo tanto

V=0hA(1zh)2 ⁣dz=Ah3:V = \int_0^h A\Bigl(1 - \frac zh\Bigr)^2\dd z = \frac{Ah}{3} :

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\alpha > 0 converge Dρα ⁣dx ⁣dy\iint_{D}\rho^{-\alpha}\,\dd x\,\dd yen el disco unitario perforado DD (límite sobre los anillos ερ1\varepsilon \leq \rho \leq 1)? en polares coordenadas,

02π ⁣ ⁣ε1ραρ ⁣dρ ⁣dα=2πε1ρ1α ⁣dρ,\int_0^{2\pi}\!\!\int_\varepsilon^1\rho^{-\alpha}\, \rho\,\dd\rho\,\dd\alpha = 2\pi\int_\varepsilon^1\rho^{1-\alpha}\,\dd\rho ,

que converge como ε0\varepsilon \to 0 si 1α>11 - \alpha > -1, es decir α<2\alpha < 2: en la dimensión 22 la singularidad crítica exponente es la dimensión misma, el ρ\rho extra del Jacobiano suavizando la singularidad por un poder. (Asimismo α<3\alpha < 3 para una singularidad puntual en el espacio, vía r2r^2.) 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\iiint 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\abs{\det} normalmente voltea las señales exactamente cuando el mapa invierte la orientación. (iii) The polar factor ρ\rho:  ⁣dx ⁣dy=ρ ⁣dρ ⁣dα\dd x\,\dd y = \rho\,\dd\rho\,\dd\alpha, no  ⁣dρ ⁣dα\dd\rho\,\dd\alpha — el error más común en todo capítulo; El análisis dimensional lo capta, ya que  ⁣dρ ⁣dα\dd\rho\,\dd\alpha 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 π\sqrt\pi 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/πn1/\sqrt{\pi 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 γy2 ⁣dx+x ⁣dy\int_\gamma y^2\,\dd x + x\,\dd y a lo largo de: (a) el segmento de (0,0)(0,0) a (1,1)(1,1); (b) el arco de parábola y=x2y = x^2 de (0,0)(0,0) a (1,1)(1,1). ¿El formulario es exacto?

Solución

Solución de Ejercicio 20.1.

(a) Segmento γ(t)=(t,t)\gamma(t) = (t, t), t[0,1]t \in [0,1]:

γy2 ⁣dx+x ⁣dy=01(t2+t) ⁣dt=13+12=56.\int_\gamma y^2\dd x + x\dd y = \int_0^1 (t^2 + t)\,\dd t = \frac13 + \frac12 = \frac56 .

(b) Parábola γ(t)=(t,t2)\gamma(t) = (t, t^2):

01(t41+t2t) ⁣dt=15+23=1315.\int_0^1 \bigl(t^4\cdot 1 + t\cdot 2t\bigr)\dd t = \frac15 + \frac23 = \frac{13}{15} .

Los dos valores difieren, por lo que la integral depende de la trayectoria: la forma es no exacto — consistentemente, Py=2y1=QxP_y = 2y \neq 1 = Q_x, entonces ni siquiera está cerrado.

Ejercicio 20.2

Muestra que ω=(2xy+y3) ⁣dx+(x2+3xy2+1) ⁣dy\omega = (2xy + y^3)\,\dd x + (x^2 + 3xy^2 + 1)\,\dd y está cerrado en R2\R^2, busque un potencial y calcule γω\int_\gamma\omega a lo largo de cualquier arco desde (0,0)(0, 0) hasta (1,2)(1, 2).

Solución

Solución de Ejercicio 20.2.

P=2xy+y3P = 2xy + y^3, Q=x2+3xy2+1Q = x^2 + 3xy^2 + 1: Py=2x+3y2=QxP_y = 2x + 3y^2 = Q_x, cerrado el R2\R^2. Busque ff con fx=Pf_x = P: f=x2y+xy3+g(y)f = x^2y + xy^3 + g(y); entonces fy=x2+3xy2+g(y)=Qf_y = x^2 + 3xy^2 + g'(y) = Qfuerza a g(y)=1g'(y) = 1, diga g(y)=yg(y) = y. entonces

f(x,y)=x2y+xy3+yf(x, y) = x^2y + xy^3 + y

es un potencial (R2\R^2 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)(0,0) hasta (1,2)(1,2):

γω=f(1,2)f(0,0)=2+8+2=12.\int_\gamma\omega = f(1, 2) - f(0, 0) = 2 + 8 + 2 = 12 .

Ejercicio 20.3

Calcule D(x+y) ⁣dx ⁣dy\iint_D (x + y)\,\dd x\,\dd y donde DD es el dominio delimitado por y=x2y = x^2 y y=xy = x (0x10 \leq x \leq 1), en ambos órdenes de integración.

Solución

Solución de Ejercicio 20.3.

El dominio es {0x1, x2yx}\{0 \leq x \leq 1,\ x^2 \leq y \leq x\}. yy-primero:

01 ⁣x2x(x+y) ⁣dy ⁣dx=01(x(xx2)+x2x42) ⁣dx=01(3x22x3x42) ⁣dx=1214110=320.\int_0^1\!\int_{x^2}^{x}(x + y)\,\dd y\,\dd x = \int_0^1\Bigl(x(x - x^2) + \frac{x^2 - x^4}{2}\Bigr)\dd x = \int_0^1\Bigl(\frac{3x^2}{2} - x^3 - \frac{x^4}{2}\Bigr)\dd x = \frac12 - \frac14 - \frac1{10} = \frac{3}{20} .

xx-primero: el corte a la altura y[0,1]y \in [0, 1] es yxyy \leq x \leq \sqrt y, por lo que

01 ⁣yy(x+y) ⁣dx ⁣dy=01(yy22+y(yy)) ⁣dy=1416+2513=320.\int_0^1\!\int_{y}^{\sqrt y}(x + y)\,\dd x\,\dd y = \int_0^1\Bigl(\frac{y - y^2}{2} + y(\sqrt y - y)\Bigr)\dd y = \frac14 - \frac16 + \frac25 - \frac13 = \frac{3}{20} .

Ejercicio 20.4 ★★

Usando coordenadas polares, calcule D ⁣dx ⁣dy(1+x2+y2)2\iint_D \frac{\dd x\,\dd y}{(1 + x^2 + y^2)^2} en todo el plano (como límite sobre los discos), y Dxy ⁣dx ⁣dy\iint_{D'} xy\,\dd x\,\dd y sobre el cuarto de disco D={x,y0, x2+y21}D' = \{x, y \geq 0,\ x^2 + y^2 \leq 1\}.

Solución

Solución de Ejercicio 20.4.

Primera integral. En el disco DRD_R, en coordenadas polares:

DR ⁣dx ⁣dy(1+x2+y2)2=02π ⁣ ⁣0Rρ ⁣dρ ⁣dα(1+ρ2)2=2π[12(1+ρ2)]0R=π(111+R2)Rπ.\iint_{D_R}\frac{\dd x\,\dd y}{(1 + x^2 + y^2)^2} = \int_0^{2\pi}\!\!\int_0^R \frac{\rho\,\dd\rho\,\dd\alpha}{(1 + \rho^2)^2} = 2\pi\Bigl[-\frac{1}{2(1 + \rho^2)}\Bigr]_0^R = \pi\Bigl(1 - \frac{1}{1 + R^2}\Bigr) \xrightarrow[R \to \infty]{} \pi .

Segunda integral. El cuarto de disco es 0απ20 \leq \alpha \leq \frac\pi2, 0ρ10 \leq \rho \leq 1y xy=ρ2cosαsinαxy = \rho^2\cos\alpha\sin\alpha:

Dxy ⁣dx ⁣dy=0π/2 ⁣ ⁣cosαsinα ⁣dα01ρ3 ⁣dρ=1214=18.\iint_{D'}xy\,\dd x\,\dd y = \int_0^{\pi/2}\!\!\cos\alpha\sin\alpha\,\dd\alpha \int_0^1 \rho^3\,\dd\rho = \frac12\cdot\frac14 = \frac18 .

Ejercicio 20.5 ★★

Calcule el área encerrado por el astroide x=cos3tx = \cos^3 t, y=sin3ty = \sin^3 t, t[0,2π]t \in [0, 2\pi], usando Corolario 20.22. (Linearize sin2tcos2t\sin^2 t\cos^2 t.)

Solución

Solución de Ejercicio 20.5.

Por Corolario 20.22 con x=cos3tx = \cos^3 t, y=sin3ty = \sin^3 t: x=3cos2tsintx' = -3\cos^2 t\sin t, y=3sin2tcosty' = 3\sin^2 t\cos t, entonces

xyyx=3cos4tsin2t+3sin4tcos2t=3sin2tcos2t=34sin22t=38(1cos4t).xy' - yx' = 3\cos^4 t\sin^2 t + 3\sin^4 t\cos^2 t = 3\sin^2 t\cos^2 t = \frac{3}{4}\sin^2 2t = \frac{3}{8}(1 - \cos 4t) .

Por lo tanto

A=1202π38(1cos4t) ⁣dt=3162π=3π8.A = \frac12\int_0^{2\pi}\frac38(1 - \cos 4t)\,\dd t = \frac{3}{16}\cdot 2\pi = \frac{3\pi}{8} .

(El astroide cabe en el disco unitario de área π\pi; tres octavos de π\pi 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+y2z = x^2 + y^2y por arriba por el plano z=1z = 1, por ambos métodos: apilamiento (integre 1x2y21 - x^2 - y^2 sobre el disco unitario, polar coordenadas) y corte (los cortes horizontales son discos de radio z\sqrt z).

Solución

Solución de Ejercicio 20.6.

Apilado: encima de cada (x,y)(x, y) del disco unitario DD, zz va desde x2+y2x^2 + y^2 a 11:

V=D(1x2y2) ⁣dx ⁣dy=02π ⁣ ⁣01(1ρ2)ρ ⁣dρ ⁣dα=2π(1214)=π2.V = \iint_D \bigl(1 - x^2 - y^2\bigr)\dd x\,\dd y = \int_0^{2\pi}\!\!\int_0^1 (1 - \rho^2)\rho\,\dd\rho\,\dd\alpha = 2\pi\Bigl(\frac12 - \frac14\Bigr) = \frac\pi2 .

Rebanar: el corte en altura z[0,1]z \in [0, 1] es el disco x2+y2zx^2 + y^2 \leq z, de área πz\pi z:

V=01πz ⁣dz=π2.V = \int_0^1 \pi z\,\dd z = \frac\pi2 .

Ejercicio 20.7 ★★

(Atracción gravitacional de una pelota — teorema de Newton, especial caso) Demuestre que el volumen de la cáscara esférica arba \leq r \leq bes 43π(b3a3)\frac43\pi(b^3 - a^3)y calcule B ⁣dx ⁣dy ⁣dzr\iiint_{B} \frac{\dd x\,\dd y\,\dd z}{r} sobre la bola BBde radio RR (rr 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φ ⁣dr ⁣dθ ⁣dφr^2\cos\varphi\,\dd r\,\dd\theta\,\dd\varphi (Ejemplo 20.27), y la parte angular integra a 4π4\pi (2π2\pi de θ\theta, π/2π/2cos=2\int_{-\pi/2}^{\pi/2}\cos = 2). El volumen de la cáscara es

ab4πr2 ⁣dr=43π(b3a3).\int_a^b 4\pi r^2\,\dd r = \frac43\pi\bigl(b^3 - a^3\bigr).

Para la segunda integral, el integrando 1/r1/r depende únicamente de rr:

B ⁣dx ⁣dy ⁣dzr=0R4πr21r ⁣dr=4πR22=2πR2.\iiint_B \frac{\dd x\,\dd y\,\dd z}{r} = \int_0^R 4\pi r^2\cdot\frac1r\,\dd r = 4\pi\,\frac{R^2}{2} = 2\pi R^2 .

(El integrando explota en el origen, pero sin causar daño: r2/r=rr^2/r = r es continuo — la integral sobre los shells εrR\varepsilon \leq r \leq Rconverge como ε0\varepsilon \to 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 UU en forma de estrella con con respecto a 00 (es decir, MU[0,M]UM \in U \Rightarrow [0, M] \subseteq U) y ω=P ⁣dx+Q ⁣dy\omega = P\dd x + Q\dd y un formulario cerrado C1\mathcal{C}^1 en UU. Definir

f(x,y)=01(xP(tx,ty)+yQ(tx,ty)) ⁣dt.f(x, y) = \int_0^1 \bigl(x\,P(tx, ty) + y\,Q(tx, ty)\bigr)\dd t .

Usando la diferenciación bajo el signo integral. (Capítulo 9) y Py=QxP_y = Q_x, muestran que fx=Pf_x = P y fy=Qf_y = 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)g(t; x, y) = xP(tx, ty) + yQ(tx, ty) es C1\mathcal{C}^1 en (x,y)(x, y), continuo en tt, con parcial derivados continuo sobre [0,1]×U[0,1] \times U; diferenciación bajo el signo integral (Capítulo 9, aplicado en el compacto tt-intervalo [0,1][0,1], la dominación es automática allí) da

fx(x,y)=01(P(tx,ty)+txPx(tx,ty)+tyQx(tx,ty)) ⁣dt.f_x(x, y) = \int_0^1 \bigl(P(tx, ty) + tx\,P_x(tx, ty) + ty\,Q_x(tx, ty)\bigr)\dd t .

Usando el cierre Qx=PyQ_x = P_y:

txPx(tx,ty)+tyPy(tx,ty)=t ⁣d ⁣dt[P(tx,ty)],tx\,P_x(tx, ty) + ty\,P_y(tx, ty) = t\,\frac{\dd}{\dd t}\bigl[P(tx, ty)\bigr] ,

entonces el integrando es P(tx,ty)+t ⁣d ⁣dtP(tx,ty)= ⁣d ⁣dt[tP(tx,ty)]P(tx, ty) + t\frac{\dd}{\dd t}P(tx, ty) = \frac{\dd}{\dd t}\bigl[t\,P(tx, ty)\bigr] y

fx(x,y)=[tP(tx,ty)]01=P(x,y).f_x(x, y) = \Bigl[t\,P(tx, ty)\Bigr]_0^1 = P(x, y) .

Simétricamente fy=Qf_y = Q (se utiliza el mismo cálculo con Py=QxP_y = Q_x) al revés). Note donde entra la hipótesis: ff está definido integrando a lo largo del segmento [0,M][0, M], que se encuentra en UU Precisamente porque UU tiene forma de estrella.

Ejercicio 20.9 ★★★

(Integral de Dirichlet por doble integración) Justificar y explotar

0 ⁣ ⁣0exysinx   ⁣dy ⁣dxvs0 ⁣ ⁣0exysinx   ⁣dx ⁣dy\int_0^\infty\!\!\int_0^\infty e^{-xy}\sin x\;\dd y\,\dd x \quad\text{vs}\quad \int_0^\infty\!\!\int_0^\infty e^{-xy}\sin x\;\dd x\,\dd y

en [0,A]×[0,)[0, A] \times [0, \infty): muestra 0Asinxx ⁣dx=π20eAyysinA+cosA1+y2 ⁣dy\int_0^A \frac{\sin x}{x}\dd x = \frac\pi2 - \int_0^\infty e^{-Ay}\frac{y\sin A + \cos A}{1 + y^2}\dd y y recupera 0sinxx ⁣dx=π2\int_0^\infty \frac{\sin x}{x}\,\dd x = \frac\pi2, 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,)[0, A] \times [0, \infty) la función (x,y)exysinx(x, y) \mapsto e^{-xy}\sin xno es absolutamente integrable hasta y=y = \infty uniformemente en un sentido ingenuo, pero cada integral iterada converge y su igualdad se sigue de Fubini en [0,A]×[0,B][0, A] \times [0, B]más un límite BB \to \infty (la cola 0ABexysinx ⁣dy ⁣dx0AeBxsinxx ⁣dx0AeBx ⁣dx0\int_0^A\int_B^\infty e^{-xy}\abs{\sin x}\,\dd y\,\dd x \leq \int_0^A \frac{e^{-Bx}\abs{\sin x}}{x}\dd x \leq \int_0^A e^{-Bx} \dd x\to 0, utilizando sinxx\abs{\sin x} \leq x).

yy primero: 0exy ⁣dy=1x\int_0^\infty e^{-xy}\,\dd y = \frac1x para x>0x > 0, por lo que la primera integral es 0Asinxx ⁣dx\int_0^A \frac{\sin x}{x}\,\dd x.

xx primero: dos integraciones por partes (o tomando la parte imaginaria de 0Ae(iy)x ⁣dx\int_0^A e^{(i - y)x}\dd x) dar

0Aexysinx ⁣dx=1eAy(ysinA+cosA)1+y2.\int_0^A e^{-xy}\sin x\,\dd x = \frac{1 - e^{-Ay}(y\sin A + \cos A)}{1 + y^2} .

Integrando yy sobre [0,)[0, \infty), el término 0 ⁣dy1+y2=π2\int_0^\infty \frac{\dd y}{1 + y^2} = \frac\pi2 se escinde:

0Asinxx ⁣dx=π20eAyysinA+cosA1+y2 ⁣dy.\int_0^A \frac{\sin x}{x}\,\dd x = \frac{\pi}{2} - \int_0^\infty e^{-Ay}\,\frac{y\sin A + \cos A}{1 + y^2}\,\dd y .

El resto está limitado por 0eAyy+11+y2 ⁣dy0eAy1+y1+y2 ⁣dy0\int_0^\infty e^{-Ay}\frac{y + 1}{1 + y^2}\dd y \leq \int_0^\infty e^{-Ay}\cdot\frac{1+y}{1+y^2} \,\dd y \to 0 como AA \to \infty (convergencia dominada, o la crudo atado 1+y1+y232\frac{1 + y}{1 + y^2} \leq \frac32 dando 32A\frac{3}{2A}). Por lo tanto 0sinxx ⁣dx=π2\int_0^\infty\frac{\sin x}{x}\dd x = \frac\pi2 — 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 γ\gamma un simple curva C1\mathcal{C}^1 cerrada de longitud 2π2\pi, parametrizada por arco longitud en [0,2π][0, 2\pi], encerrando a área AA. Usando Corolario 20.22, Parseval y la desigualdad de Wirtinger (ejercicios de Capítulo 14), demostrar AπA \leq \pi, con igualdad para el círculo. (Normalize 02πx(s) ⁣ds=0\int_0^{2\pi} x(s)\dd s = 0; write 2A=x ⁣dyy ⁣dx2A = \oint x\,\dd y - y\,\dd xand bound 2A(x2+y2)2A \leq \int (x^2 + y'^2)carefully via 2A=02π(xyyx) ⁣ds2A = \int_0^{2\pi}(xy' - yx')\dd sand x2+y22xyx^2 + y'^2 \geq 2xy'.)

Solución

Solución de Ejercicio 20.10.

Parametrizar por longitud del arco s[0,2π]s \in [0, 2\pi], entonces x2+y2=1x'^2 + y'^2 = 1, y traducir para que 02πx(s) ⁣ds=0\int_0^{2\pi} x(s)\,\dd s = 0. Por Corolario 20.22,

2A=x ⁣dyy ⁣dx=02π(xyyx) ⁣ds.2A = \oint x\,\dd y - y\,\dd x = \int_0^{2\pi}\bigl(xy' - yx'\bigr)\dd s .

Integrando y ⁣dx\oint y\,\dd x por partes a lo largo del periodo (límite términos se cancelan por periodicidad), yx=yx-\int yx' = \int y'x, por lo que de hecho 2A=202πxy ⁣ds2A = 2\int_0^{2\pi}xy'\,\dd s. Entonces 2xyx2+y22xy' \leq x^2 + y'^2 da

2A02π(x2+y2) ⁣ds=02πx2+02π(1x2)=2π02π(x2x2) ⁣ds.2A \leq \int_0^{2\pi}\bigl(x^2 + y'^2\bigr)\dd s = \int_0^{2\pi} x^2 + \int_0^{2\pi}\bigl(1 - x'^2\bigr) = 2\pi - \int_0^{2\pi}\bigl(x'^2 - x^2\bigr)\dd s .

Desigualdad de Wirtinger (ejercicios Capítulo 14: para un 2π2\pi-función periódica C1\mathcal{C}^1 con media cero, x2x2\int x^2 \leq \int x'^2) hace que la última integral no sea negativa: AπA \leq \pi. La igualdad requiere igualdad en Wirtinger (x(s)=acoss+bsinsx(s) = a\cos s + b\sin s) y en 2xyx2+y22xy' \leq x^2 + y'^2(y=xy' = x puntualmente), lo que fuerza a y=asinsbcoss+cy = a\sin s - b\cos s + c: la curva es el círculo unitario (adecuadamente centrado). Desde una curva de longitud LL reescala a longitud 2π2\pi, la declaración general es AL24πA \leq \frac{L^2}{4\pi}: 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 BB de radio RR en R3\R^3, calcular Bz2 ⁣dx ⁣dy ⁣dz\iiint_B z^2\,\dd x\,\dd y\,\dd z en coordenadas esféricas, y deducir B(x2+y2+z2) ⁣dx ⁣dy ⁣dz\iiint_B (x^2 + y^2 + z^2)\,\dd x\,\dd y\,\dd z por simetría. Verifique el este último contra el cálculo del shell 0Rr24πr2 ⁣dr\int_0^R r^2\cdot4\pi r^2\,\dd r.

Solución

Solución de Ejercicio 20.11.

En coordenadas esféricas z=rsinφz = r\sin\varphi y  ⁣dx ⁣dy ⁣dz=r2cosφ ⁣dr ⁣dθ ⁣dφ\dd x\,\dd y\,\dd z = r^2\cos\varphi\,\dd r\,\dd\theta\,\dd\varphi:

Bz2=0Rr4 ⁣dr02π ⁣dθπ/2π/2sin2φcosφ ⁣dφ=R552π[sin3φ3]π/2π/2=4πR515.\iiint_B z^2 = \int_0^R r^4\,\dd r\int_0^{2\pi}\dd\theta \int_{-\pi/2}^{\pi/2}\sin^2\varphi\cos\varphi\,\dd\varphi = \frac{R^5}5\cdot2\pi\cdot \Bigl[\frac{\sin^3\varphi}3\Bigr]_{-\pi/2}^{\pi/2} = \frac{4\pi R^5}{15}.

Por la simetría de la pelota bajo coordenadas permutantes, Bx2=By2=Bz2\iiint_B x^2 = \iiint_B y^2 = \iiint_B z^2, entonces B(x2+y2+z2)=34πR515=4πR55\iiint_B(x^2 + y^2 + z^2) = 3\cdot\frac{4\pi R^5}{15} = \frac{4\pi R^5}5. Verificación de shell: 0Rr24πr2 ⁣dr=4πR55\int_0^R r^2\cdot 4\pi r^2\,\dd r = \frac{4\pi R^5}5— el integrando r2r^2 es constante en la esfera de radio rr, de área 4πr24\pi r^2.

Ejercicio 20.12 ★★

(Fórmula del cordón) Sea KK un polígono con vértices (x1,y1),,(xm,ym)(x_1, y_1), \dots, (x_m, y_m) en orden antihorario (índices modelo mm). Deduzca de Corolario 20.22 que

Area(K)=12i=1m(xiyi+1xi+1yi),\operatorname{Area}(K) = \frac12\sum_{i=1}^m \bigl(x_iy_{i+1} - x_{i+1}y_i\bigr) ,

y verifique la fórmula en el triángulo (0,0)(0,0), (1,0)(1,0), (0,1)(0,1).

Solución

Solución de Ejercicio 20.12.

Parametrizar el borde de (xi,yi)(x_i, y_i) a (xi+1,yi+1)(x_{i+1}, y_{i+1}) por γ(t)=((1t)xi+txi+1, (1t)yi+tyi+1)\gamma(t) = \bigl((1-t)x_i + tx_{i+1},\ (1-t)y_i + ty_{i+1}\bigr). Su contribución a 12(x ⁣dyy ⁣dx)\frac12\oint(x\,\dd y - y\,\dd x) es

1201(((1t)xi+txi+1)(yi+1yi)((1t)yi+tyi+1)(xi+1xi)) ⁣dt,\frac12\int_0^1\Bigl(\bigl((1-t)x_i + tx_{i+1}\bigr)(y_{i+1} - y_i) - \bigl((1-t)y_i + ty_{i+1}\bigr)(x_{i+1} - x_i)\Bigr)\dd t ,

y desde 01((1t)u+tv) ⁣dt=u+v2\int_0^1\bigl((1-t)u + tv\bigr)\dd t = \frac{u + v}2, esto equivale

14((xi+xi+1)(yi+1yi)(yi+yi+1)(xi+1xi))=12(xiyi+1xi+1yi),\frac14\Bigl((x_i + x_{i+1})(y_{i+1} - y_i) - (y_i + y_{i+1})(x_{i+1} - x_i)\Bigr) = \frac12\bigl(x_iy_{i+1} - x_{i+1}y_i\bigr),

los términos cruzados se cancelan. Sumar los bordes mm da la fórmula del cordón, de Corolario 20.22. Triangulo (0,0),(1,0),(0,1)(0,0), (1,0), (0,1): 12((0010)+(1100)+(0001))=12\frac12\bigl((0\cdot0 - 1\cdot0) + (1\cdot1 - 0\cdot0) + (0\cdot0 - 0\cdot1)\bigr) = \frac12, el área correcto.

20.6 Problema: el volumen de la pelota en la dimensión nn

Problema 20.1

Problema del fin de semana — Vn=πn/2/Γ(n2+1)V_n = \pi^{n/2}/\Gamma(\frac n2 + 1), y la extrañeza de la alta dimensiones

El disco tiene área π\pi, volumen de la bola 43π\frac43\pi — y entonces? Este problema calcula el volumen de la bola unitaria de Rn\R^n por cada nn, dos veces (mediante una recursividad de corte impulsada por el Integrales de Wallis, luego a través de la función Γ\Gamma 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\R^n, la integral se entiende como la integral iterada nn veces (cortando una coordenada a la vez, como en el capítulo n3n \leq 3); escribimos Bn(R)B_n(R)para la bola cerrada de radio RR centrado en 00, vn(R)v_n(R) para su volumen, y Vn=vn(1)V_n = v_n(1), con V0=1V_0 = 1 por convención.

Parte I — The slicing recursion.

  1. Sustituyendo xi=Ruix_i = Ru_i en cada uno de los nn iterados integrales, muestre vn(R)=VnRnv_n(R) = V_nR^n.
  2. Cortando Bn(1)B_n(1) a lo largo de su última coordenada, muestra

    Vn=Vn111(1t2)n12 ⁣dt.V_n = V_{n-1}\int_{-1}^{1}(1 - t^2)^{\frac{n-1}2}\,\dd t .
  3. Con t=sinθt = \sin\theta, identifique la integral como Integral de Wallis: 11(1t2)n12 ⁣dt=2Wn\int_{-1}^1(1 - t^2)^{\frac{n-1}2}\dd t = 2W_n, donde Wn=0π/2cosnθ ⁣dθ=0π/2sinnθ ⁣dθW_n = \int_0^{\pi/2}\cos^n\theta\,\dd\theta = \int_0^{\pi/2}\sin^n\theta\,\dd\theta.
  4. Demostrar las dos identidades de Wallis (integrar por partes; luego telescopio nWnWn1nW_nW_{n-1}):

    Wn=n1nWn2(n2),WnWn1=π2n(n1).W_n = \frac{n-1}nW_{n-2} \quad (n \geq 2), \qquad W_nW_{n-1} = \frac{\pi}{2n} \quad (n \geq 1).

Parte II — The recursion solved.

  1. Combine las preguntas 2 a 4 en la recursividad de dos pasos

    Vn=2πnVn2(n2).V_n = \frac{2\pi}{n}\,V_{n-2} \qquad (n \geq 2).
  2. Deduce el formas cerradas, para k0k \geq 0:

    V2k=πkk!,V2k+1=2k+1πk135(2k+1).V_{2k} = \frac{\pi^k}{k!}, \qquad V_{2k+1} = \frac{2^{k+1}\pi^k}{1\cdot3\cdot5\cdots (2k+1)} .
  3. Tabule V1,,V7V_1, \dots, V_7 numéricamente. Usando el relación Vn/Vn2=2π/nV_n/V_{n-2} = 2\pi/n y los valores de 2W52W_5 y 2W62W_6, demuestre que la secuencia (Vn)(V_n) aumenta hasta su máximo V5=8π2155.26V_5 = \frac{8\pi^2}{15} \approx 5.26 y disminuye a partir de ahí.
  4. Demuestre que Vn0V_n \to 0 es más rápido que cualquier patrón geométrico. secuencia, y que n1Vn\sum_{n\geq1} V_n converge: todos las bolas unitarias juntas tienen un volumen total finito.
  5. Probar la identidad generadora

    k0V2kx2k=eπx2(xR),\sum_{k\geq0} V_{2k}\,x^{2k} = \eu^{\pi x^2} \qquad (x \in \R),

    y deducir k0V2k=eπ23.14\sum_{k\geq0}V_{2k} = \eu^\pi \approx 23.14.

Parte III — Second route: Γ\Gamma and the Gauss integral.

  1. Muestra por Fubini (el integrando es un producto) que

    In=Rnex2 ⁣dx=(+et2 ⁣dt) ⁣n=πn/2,I_n = \int_{\R^n}\eu^{-\norm x^2}\dd x = \Bigl(\int_{-\infty}^{+\infty} \eu^{-t^2}\dd t\Bigr)^{\!n} = \pi^{n/2},

    la integral de Gauss nn-dimensional, entendida como una límite sobre cubos [R,R]n\intcc{-R}{R}^n.

  2. Recuperar Γ(s)=0ts1et ⁣dt\Gamma(s) = \int_0^\infty t^{s-1}\eu^{-t}\dd t (Definición 9.17). De Γ(s+1)=sΓ(s)\Gamma(s+1) = s\,\Gamma(s) (Teorema 9.18) y Γ(12)=π\Gamma(\tfrac12) = \sqrt\pi (sustituto t=u2t = u^2 e invocar la integral de Gauss), calcular

    Γ(k+1)=k!,Γ(k+32)=13(2k+1)2k+1π.\Gamma(k + 1) = k!, \qquad \Gamma\Bigl(k + \frac32\Bigr) = \frac{1\cdot3\cdots(2k+1)}{2^{k+1}}\,\sqrt\pi .
  3. Demostrar, por inducción mediante la recursividad de la pregunta 5, la fórmula única

    Vn=πn/2Γ(n2+1)(n1),V_n = \frac{\pi^{n/2}}{\Gamma\bigl(\frac n2 + 1\bigr)} \qquad (n \geq 1),

    y verifique que reproduzca tanto formas cerradas de la pregunta 6.

  4. Muestra 0er2rn1 ⁣dr=12Γ(n2)\int_0^\infty \eu^{-r^2}r^{n-1}\dd r = \tfrac12\Gamma\bigl(\tfrac n2\bigr) y deduce el identidad

    In=nVn0er2rn1 ⁣dr.I_n = n\,V_n\int_0^\infty \eu^{-r^2}\,r^{n-1}\,\dd r.

    Interpretarlo: la masa gaussiana de Rn\R^n es recogidos a lo largo de conchas esféricas cuyo “(n1)(n-1)-dimensional área” en el radio rr es nVnrn1nV_nr^{n-1} — ambas partes ahora están probadas de forma independiente, por lo que la interpretación no cuesta nada.

  5. Establezca sn1=nVns_{n-1} = nV_n (el área de la unidad de esfera Sn1S^{n-1}, consistente con vn(R)=0Rsn1rn1 ⁣drv_n(R) = \int_0^R s_{n-1}r^{n-1}\dd r). Tabular s0,,s3s_0, \dots, s_3 y marque s1=2πs_1 = 2\pi, s2=4πs_2 = 4\pi, s3=2π2s_3 = 2\pi^2.

Parte IV — High dimensions are strange.

  1. De la fórmula de Stirling (Teorema 6.13) aplicado a k!k!, mostrar incluso n=2kn = 2k:

    Vn1πn(2πen)n/2(n, n even),V_n \sim \frac{1}{\sqrt{\pi n}} \Bigl(\frac{2\pi\eu}{n}\Bigr)^{n/2} \qquad (n \to \infty, \ n \text{ even}),

    y explicar por qué la misma desintegración supergeométrica se limita se extiende al impar nn mediante la recursividad.

  2. La bola unitaria se encuentra en el cubo [1,1]n\intcc{-1}1^n de volumen 2n2^n. Calcule la relación de llenado Vn/2nV_n/2^n para n=2,3,10n = 2, 3, 10, y muestra que tiende a 00: en alto dimensión, esencialmente todo el cubo se encuentra en su esquinas.
  3. Demuestre que la fracción de vn(1)v_n(1) que se encuentra dentro la distancia ε\varepsilon de la esfera límite es 1(1ε)n11 - (1 - \varepsilon)^n \to 1; numéricamente, ¿qué fracción de una bola de dimensiones 100100 se encuentra en el capa exterior de espesor 1%1\%?
  4. Demuestre el Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}} (monotonicity of (Wn)(W_n), the ratio Wn/Wn21W_n/W_{n-2} \to 1, and WnWn1=π2nW_nW_{n-1} = \frac\pi{2n}) asintótico de Wallis y el límite inferior Wnπ2(n+1)W_n \geq \sqrt{\dfrac{\pi}{2(n+1)}} para todos los nn.
  5. (Concentración en una losa) La fracción de la unidad bola con la primera coordenada más allá de δ\delta es δ1(1x2)n12 ⁣dx/(2Wn)\int_\delta^1(1 - x^2)^{\frac{n-1}2}\dd x \,\big/\, (2W_n). Usando 1ueu1 - u \leq \eu^{-u} y la cola atado δeax2 ⁣dxeaδ22aδ\int_\delta^\infty \eu^{-a x^2}\dd x \leq \frac{\eu^{-a\delta^2}}{2a\delta}, muestra esto la fracción es como máximo

    e(n1)δ2/2(n1)δ/2πn+1\frac{\eu^{-(n-1)\delta^2/2}}{(n-1)\,\delta} \Big/ \sqrt{\frac{2\pi}{n+1}}

    y concluimos: para δ=s/n1\delta = s/\sqrt{n-1}, todos menos un la fracción O(es2/2/s)O(\eu^{-s^2/2}/s) de la pelota se encuentra en el losa x1s/n1\abs{x_1} \leq s/\sqrt{n-1}. una bola de alto dimensión es, estadísticamente, un panqueque delgado en cada dirección a la vez.

  6. Reúna las preguntas 16 a 19 en un párrafo: dónde el volumen de Bn(1)B_n(1) se encuentra (cerca del límite esfera, pero dentro de O(1/n)O(1/\sqrt 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.

  1. (Simplex) Deje Δn={xRn:xi0, xi1}\Delta_n = \{x \in \R^n : x_i \geq 0,\ \sum x_i \leq 1\}. Demuestre por corte e inducción que vol(Δn)=1n!\operatorname{vol}(\Delta_n) = \frac1{n!}.
  2. (Politopo cruzado) Deduzca que Cn={x:xi1}C_n = \{x : \sum\abs{x_i} \leq 1\} tiene volumen 2nn!\frac{2^n}{n!}, y verificar el sándwich CnBn(1)[1,1]nC_n \subseteq B_n(1) \subseteq \intcc{-1}1^n a nivel de volúmenes: 2nn!Vn2n\frac{2^n}{n!} \leq V_n \leq 2^n.
  3. Calcule V4V_4 de una tercera manera: corte R4=R2×R2\R^4 = \R^2 \times \R^2, integre el área del disco (z,w)(z, w) el disco (x,y)(x, y) en coordenadas polares y recuperar V4=π22V_4 = \frac{\pi^2}2.
  4. (Monte Carlo en problemas) Se sortea un punto uniformemente en el cubo [1,1]20\intcc{-1}1^{20}. Demuestre que el la probabilidad de que caiga en la bola inscrita es V20/2202.5108V_{20}/2^{20} \approx 2.5\cdot10^{-8}, más o menos Se necesitan cuarenta millones de sorteos antes del primer golpe. se espera: estimación de VnV_n mediante muestreo de rechazo colapsa en alta dimensión (la maldición de dimensionalidad).
  5. Síntesis. Dos derivaciones independientes se encontraron en Vn=πn/2/Γ(n2+1)V_n = \pi^{n/2}/\Gamma(\frac n2 + 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, Γ\Gamma, 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)B_n(R) se describe mediante límites iterados. RxnR-R \leq x_n \leq R, luego xn1R2xn2\abs{x_{n-1}} \leq \sqrt{R^2 - x_n^2}, y así sucesivamente; sustituyendo xi=Ruix_i = Ru_i en cada uno de los nn integrales de una variable multiplica cada una por RR y asigna los límites a los de Bn(1)B_n(1): vn(R)=Rnvn(1)=VnRnv_n(R) = R^n\,v_n(1) = V_nR^n.

2. Cortando a lo largo de xn=tx_n = t: la rebanada de Bn(1)B_n(1) es la pelota Bn1(1t2)B_{n-1}\bigl(\sqrt{1 - t^2}\bigr), entonces, por pregunta 1,

Vn=11vn1(1t2) ⁣dt=Vn111(1t2)n12 ⁣dt.V_n = \int_{-1}^1 v_{n-1}\bigl(\sqrt{1 - t^2}\bigr)\,\dd t = V_{n-1}\int_{-1}^{1}(1 - t^2)^{\frac{n-1}2}\,\dd t .

3. Con t=sinθt = \sin\theta,  ⁣dt=cosθ ⁣dθ\dd t = \cos\theta\,\dd\thetay(1t2)n12=cosn1θ(1 - t^2)^{\frac{n-1}2} = \cos^{n-1}\thetaen [π/2,π/2]\intcc{-\pi/2}{\pi/2}:

11(1t2)n12 ⁣dt=π/2π/2cosnθ ⁣dθ=20π/2cosnθ ⁣dθ=2Wn,\int_{-1}^1(1 - t^2)^{\frac{n-1}2}\dd t = \int_{-\pi/2}^{\pi/2}\cos^n\theta\,\dd\theta = 2\int_0^{\pi/2}\cos^n\theta\,\dd\theta = 2W_n,

y θπ2θ\theta \mapsto \frac\pi2 - \theta intercambia el seno y formas cosenos de WnW_n.

4. Escribe sinn=sinn2(1cos2)\sin^n = \sin^{n-2}(1 - \cos^2) y integrar sinn2coscos\int\sin^{n-2}\cos\cdot\cos por partes (v=sinn1n1v = \frac{\sin^{n-1}}{n-1}):

Wn=Wn2Wnn1Wn=n1nWn2.W_n = W_{n-2} - \frac{W_n}{n-1} \quad\Longrightarrow\quad W_n = \frac{n-1}{n}W_{n-2}.

De ahí nWnWn1=(n1)Wn1Wn2nW_nW_{n-1} = (n-1)W_{n-1}W_{n-2}: la secuencia (nWnWn1)(nW_nW_{n-1}) es constante, igual a 1W1W0=1π21\cdot W_1W_0 = 1\cdot\frac\pi2, por lo que WnWn1=π2nW_nW_{n-1} = \frac{\pi}{2n}.

5. Las preguntas 2–3 dan Vn=2WnVn1V_n = 2W_nV_{n-1}, dos veces:

Vn=2Wn2Wn1Vn2=4π2nVn2=2πnVn2.V_n = 2W_n\cdot 2W_{n-1}\,V_{n-2} = 4\,\frac{\pi}{2n}\,V_{n-2} = \frac{2\pi}n\,V_{n-2}.

6. De V0=1V_0 = 1: V2k=2π2kV2k2=πkV2k2V_{2k} = \frac{2\pi}{2k}V_{2k-2} = \frac\pi kV_{2k-2}, entonces V2k=πkk!V_{2k} = \frac{\pi^k}{k!} por inducción. De V1=2V_1 = 2: V2k+1=2π2k+1V2k1V_{2k+1} = \frac{2\pi}{2k+1}V_{2k-1}, entonces

V2k+1=2j=1k2π2j+1=2k+1πk135(2k+1).V_{2k+1} = 2\prod_{j=1}^k\frac{2\pi}{2j+1} = \frac{2^{k+1}\pi^k}{1\cdot3\cdot5\cdots(2k+1)} .

7. V1=2V_1 = 2, V2=π3.142V_2 = \pi \approx 3.142, V3=4π34.189V_3 = \frac{4\pi}3 \approx 4.189, V4=π224.935V_4 = \frac{\pi^2}2 \approx 4.935, V5=8π2155.264V_5 = \frac{8\pi^2}{15} \approx 5.264, V6=π365.168V_6 = \frac{\pi^3}6 \approx 5.168, V7=16π31054.725V_7 = \frac{16\pi^3}{105} \approx 4.725. La relación de un paso es Vn/Vn1=2WnV_n/V_{n-1} = 2W_n, y (Wn)(W_n) está disminuyendo (sinnsinn1\sin^n \leq \sin^{n-1} puntualmente). Ahora 2W5=24523=1615>12W_5 = 2\cdot\frac45\cdot\frac23 = \frac{16}{15} > 1 mientras que 2W6=2563412π2=5π16<12W_6 = 2\cdot\frac56\cdot\frac34\cdot\frac12\cdot\frac\pi2 = \frac{5\pi}{16} < 1: los ratios superan 11 hasta n=5n = 5 y están por debajo de 11 desde n=6n = 6 en — (Vn)(V_n) aumenta a su máximo V5V_5 y luego disminuye.

8. Para n13>4πn \geq 13 > 4\pi: Vn/Vn2=2π/n<12V_n/V_{n-2} = 2\pi/n < \tfrac12, entonces VnC2n/2V_{n} \leq C\cdot 2^{-n/2} con un fijo constante; mejor, para cualquier q>0q > 0, 2π/n<q22\pi/n < q^2 para grandes nn, entonces Vn/qn0V_n/q^n \to 0: la decadencia supera a toda geometría secuencia. La convergencia de Vn\sum V_n se deriva de la relación Vn/Vn20V_n/V_{n-2} \to 0 (compárese con una serie geométrica de algún rango en).

9. k0V2kx2k=k0(πx2)kk!=eπx2\sum_{k\geq0}V_{2k}x^{2k} = \sum_{k\geq0}\frac{(\pi x^2)^k}{k!} = \eu^{\pi x^2}, el serie exponencial (Capítulo 11), convergente para cada xx. En x=1x = 1: kV2k=eπ23.14\sum_kV_{2k} = \eu^\pi \approx 23.14.

10. En el cubo [R,R]n\intcc{-R}R^n el integrando es el producto iexi2\prod_i\eu^{-x_i^2}, por lo que la integral iterada factoriza: (RRet2 ⁣dt)n\bigl(\int_{-R}^R\eu^{-t^2}\dd t\bigr)^n. Dejar RR \to \infty y usar Ret2 ⁣dt=π\int_\R\eu^{-t^2}\dd t = \sqrt\pi (Ejemplo 20.17): In=πn/2I_n = \pi^{n/2}.

11. t=u2t = u^2 da Γ(12)=0t1/2et ⁣dt=20eu2 ⁣du=π\Gamma(\tfrac12) = \int_0^\infty t^{-1/2}\eu^{-t}\dd t = 2\int_0^\infty\eu^{-u^2}\dd u = \sqrt\pi. Iterando Γ(s+1)=sΓ(s)\Gamma(s+1) = s\Gamma(s): Γ(k+1)=k!Γ(1)=k!\Gamma(k+1) = k!\,\Gamma(1) = k!, y

Γ(k+32)=(k+12)(k12)12Γ(12)=(2k+1)(2k1)12k+1π.\Gamma\Bigl(k + \frac32\Bigr) = \Bigl(k + \frac12\Bigr)\Bigl(k - \frac12\Bigr)\cdots \frac12\cdot\Gamma\Bigl(\frac12\Bigr) = \frac{(2k+1)(2k-1)\cdots1}{2^{k+1}}\,\sqrt\pi .

12. Establecer Fn=πn/2/Γ(n2+1)F_n = \pi^{n/2}/\Gamma(\frac n2 + 1). Desde Γ(n2+1)=n2Γ(n2)=n2Γ(n22+1)\Gamma(\frac n2 + 1) = \frac n2\,\Gamma(\frac n2) = \frac n2\,\Gamma(\frac{n-2}2 + 1), obtenemos Fn=2πnFn2F_n = \frac{2\pi}nF_{n-2}: la misma recursividad que VnV_n (pregunta 5). Bases: F1=π/Γ(32)=π/(π2)=2=V1F_1 = \sqrt\pi/\Gamma(\frac32) = \sqrt\pi/(\frac{\sqrt\pi}2) = 2 = V_1y F2=π/Γ(2)=π=V2F_2 = \pi/\Gamma(2) = \pi = V_2. Por inducción Vn=FnV_n = F_n para todos nn; La pregunta 11 convierte esto nuevamente en los dos formas cerradas de pregunta 6.

13. Con r=tr = \sqrt t: 0er2rn1 ⁣dr=120tn21et ⁣dt=12Γ(n2)\int_0^\infty \eu^{-r^2}r^{n-1}\dd r = \frac12\int_0^\infty t^{\frac n2 - 1}\eu^{-t}\dd t = \frac12\Gamma(\frac n2). Por lo tanto

nVn0er2rn1 ⁣dr=Vnn2Γ(n2)=VnΓ(n2+1)=πn/2=In.n\,V_n\int_0^\infty\eu^{-r^2}r^{n-1}\dd r = V_n\cdot\frac n2\,\Gamma\Bigl(\frac n2\Bigr) = V_n\,\Gamma\Bigl(\frac n2 + 1\Bigr) = \pi^{n/2} = I_n .

Una vez probadas ambas partes, la identidad puede ser leer como la descomposición de capas de la integral gaussiana: la esfera de radio rr lleva área nVnrn1nV_nr^{n-1}, y el gaussiano El peso er2\eu^{-r^2} está integrado sobre las carcasas.

14. s0=V1=2s_0 = V_1 = 2 (la esfera 00 tiene dos puntos), s1=2V2=2πs_1 = 2V_2 = 2\pi, s2=3V3=4πs_2 = 3V_3 = 4\pi, s3=4V4=2π2s_3 = 4V_4 = 2\pi^2; y 0Rsn1rn1 ⁣dr=VnRn=vn(R)\int_0^R s_{n-1}r^{n-1}\dd r = V_nR^n = v_n(R): área es la derivada radial del volumen.

15. Para n=2kn = 2k, Stirling (Teorema 6.13) da k!2πk(k/e)kk! \sim \sqrt{2\pi k}\,(k/\eu)^k, entonces

V2k=πkk!(πe/k)k2πk=1πn(2πen)n/2(n=2k).V_{2k} = \frac{\pi^k}{k!} \sim \frac{(\pi\eu/k)^k}{\sqrt{2\pi k}} = \frac1{\sqrt{\pi n}}\Bigl(\frac{2\pi\eu}n\Bigr)^{n/2} \qquad (n = 2k).

Para nn impar: V2k+1=2W2k+1V2k2V2kV_{2k+1} = 2W_{2k+1}V_{2k} \leq 2V_{2k}, entonces se mantienen los mismos límites de desintegración supergeométrica (hasta un factor 22 y un desplazamiento de uno en el exponente) — por cada q>0q > 0, Vn=o(qn)V_n = o(q^n).

16.V2/4=π/40.785V_2/4 = \pi/4 \approx 0.785; V3/8=π/60.524V_3/8 = \pi/6 \approx 0.524; V10/210=π512010240.0025V_{10}/2^{10} = \frac{\pi^5}{120\cdot1024} \approx 0.0025. En general Vn/2nVn2/2n2=2π4n=π2n0\frac{V_n/2^n}{V_{n-2}/2^{n-2}} = \frac{2\pi}{4n} = \frac{\pi}{2n} \to 0: la relación tiende a 00 (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ε1 - \varepsilontiene volumen Vn(1ε)nV_n(1-\varepsilon)^n, por lo que la exterior shell lleva la fracción 1(1ε)n11 - (1 - \varepsilon)^n \to 1. Para n=100n = 100, ε=0.01\varepsilon = 0.01: (0.99)100=e100ln0.99e1.0050.366(0.99)^{100} = \eu^{100\ln0.99} \approx \eu^{-1.005} \approx 0.366: aproximadamente 63%63\% de la pelota se encuentra dentro de 1%1\% de su superficie.

18. (Wn)(W_n) disminuye, por lo que WnWn1Wn2=nn1WnW_n \leq W_{n-1} \leq W_{n-2} = \frac{n}{n-1}W_n: apretando, Wn1/Wn1W_{n-1}/W_n \to 1. Multiplicando por WnWn1=π2nW_nW_{n-1} = \frac\pi{2n}: Wn2π2nW_n^2 \sim \frac\pi{2n}, es decir Wnπ/(2n)W_n \sim \sqrt{\pi/(2n)}. Límite inferior: Wn2WnWn+1=π2(n+1)W_n^2 \geq W_nW_{n+1} = \frac{\pi}{2(n+1)}, entonces Wnπ/(2(n+1))W_n \geq \sqrt{\pi/(2(n+1))} por cada nn.

19. Numerador: 1x2ex21 - x^2 \leq \eu^{-x^2} da (1x2)n12e(n1)x2/2(1 - x^2)^{\frac{n-1}2} \leq \eu^{-(n-1)x^2/2}, y con a=n12a = \frac{n-1}2,

δ1(1x2)n12 ⁣dxδeax2 ⁣dxδxδeax2 ⁣dx=eaδ22aδ=e(n1)δ2/2(n1)δ.\int_\delta^1(1 - x^2)^{\frac{n-1}2}\dd x \leq \int_\delta^\infty\eu^{-ax^2}\dd x \leq \int_\delta^\infty\frac x\delta\,\eu^{-ax^2}\dd x = \frac{\eu^{-a\delta^2}}{2a\delta} = \frac{\eu^{-(n-1)\delta^2/2}}{(n-1)\delta}.

Denominador: 2Wn2π/(n+1)2W_n \geq \sqrt{2\pi/(n+1)} por la pregunta 18. Al dividir se obtiene el límite mostrado. Para δ=s/n1\delta = s/\sqrt{n-1} se convierte en n+12π(n1)  es2/2/s=O(es2/2/s)\sqrt{\tfrac{n+1}{2\pi(n-1)}}\;\eu^{-s^2/2}/s = O\bigl( \eu^{-s^2/2}/s\bigr), uniformemente en nn: fuera de la losa x1s/n1\abs{x_1} \leq s/\sqrt{n-1} casi no hay volumen, por ss 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)B_n(1) tiene norma cerca de 11 (pregunta 17: radial concentración cerca de la esfera), sin embargo, cada uno de sus nn las coordenadas son pequeñas, de orden 1/n1/\sqrt n (pregunta 19), lo cual es consistente desde las coordenadas nn del tamaño 1/n1/\sqrt n tiene norma del pedido 11. 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\Delta_n en xn=t[0,1]x_n = t \in \intcc01: el la porción es {xRn1:xi0, xi1t}=(1t)Δn1\{x' \in \R^{n-1} : x_i \geq 0,\ \sum x_i \leq 1 - t\} = (1-t)\Delta_{n-1}, de volumen (1t)n1vol(Δn1)(1-t)^{n-1}\operatorname{vol}(\Delta_{n-1}) por homogeneidad. entonces

vol(Δn)=vol(Δn1)01(1t)n1 ⁣dt=vol(Δn1)nvol(Δn)=1n!.\operatorname{vol}(\Delta_n) = \operatorname{vol}(\Delta_{n-1})\int_0^1(1 - t)^{n-1}\dd t = \frac{\operatorname{vol}(\Delta_{n-1})}{n} \quad\Longrightarrow\quad \operatorname{vol}(\Delta_n) = \frac1{n!}\,.

22. Los orantes de signo 2n2^n cortan CnC_n en 2n2^n copias de Δn\Delta_n (con superposiciones insignificantes en las coordenadas hiperplanos): vol(Cn)=2nn!\operatorname{vol}(C_n) = \frac{2^n}{n!}. si xi1\sum\abs{x_i} \leq 1 luego xi2(xi)21\sum x_i^2 \leq \bigl(\sum\abs{x_i}\bigr)^2 \leq 1: CnBn(1)C_n \subseteq B_n(1); y Bn(1)[1,1]nB_n(1) \subseteq \intcc{-1}1^n desde xix\abs{x_i} \leq \norm x. Por lo tanto 2nn!Vn2n\frac{2^n}{n!} \leq V_n \leq 2^n — consistente con la pregunta 15, que coloca a VnV_n entre los Escalas factoriales y geométricas.

23. Para (x,y)(x, y) en el disco unitario, la porción de B4(1)B_4(1) es el disco de radio 1x2y2\sqrt{1 - x^2 - y^2} en el (z,w)(z, w)-avión, de área π(1x2y2)\pi(1 - x^2 - y^2). en polares coordenadas:

V4=x2+y21π(1x2y2) ⁣dx ⁣dy=π02π ⁣ ⁣01(1ρ2)ρ ⁣dρ ⁣dα=π2π14=π22,V_4 = \iint_{x^2+y^2\leq1}\pi(1 - x^2 - y^2)\,\dd x\,\dd y = \pi\int_0^{2\pi}\!\!\int_0^1(1 - \rho^2)\rho\, \dd\rho\,\dd\alpha = \pi\cdot2\pi\cdot\frac14 = \frac{\pi^2}2 ,

de acuerdo con la pregunta 6.

24. La probabilidad es la relación de volumen. V20220=π1010!2200.025810485762.5108\dfrac{V_{20}}{2^{20}} = \dfrac{\pi^{10}}{10!\cdot2^{20}} \approx \dfrac{0.0258}{1\,048\,576} \approx 2.5\cdot10^{-8}. El número de sorteos hasta el primer acierto es de orden inverso, sobre 41074\cdot10^7: un muestreador de rechazo que funcionó maravillosamente para el disco (π/4\pi/4 de visitas) es inútil en la dimensión 2020 — 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 InI_n, el cambio polar de variables a través de la integral de Gauss de Ejemplo 20.17 y la función Γ\Gamma. ecuación funcional. Se encuentran en Vn=πn/2/Γ(n2+1)V_n = \pi^{n/2}/\Gamma(\frac n2 + 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.