UTEC · FICA

Matemática IV

Uso de Coordenadas Polares

Índice

Teoría matemática

Cambio de variable

Para integrar sobre regiones con simetría circular se utiliza la transformación estándar de coordenadas polares:

x=rcos⁡θ,y=rsin⁡θ,x=r\cos\theta,\qquad y=r\sin\theta,

con r≥0r\ge0 y θ∈[0,2π)\theta\in[0,2\pi).

Jacobiano

El área elemental en polares no es dr dθdr\,d\theta; al derivar la transformación se obtiene

dA=∣det⁡∂(x,y)∂(r,θ)∣dr dθ=∣det⁡(cos⁡θ−rsin⁡θsin⁡θrcos⁡θ)∣dr dθ=r dr dθ.dA=\left|\det\frac{\partial(x,y)}{\partial(r,\theta)}\right|dr\,d\theta=\left|\det\begin{pmatrix}\cos\theta & -r\sin\theta\\ \sin\theta & r\cos\theta\end{pmatrix}\right|dr\,d\theta=r\,dr\,d\theta.

Por tanto

dA=r dr dθ.\boxed{dA=r\,dr\,d\theta.}

El factor rr refleja que un anillo de radio rr y grosor drdr tiene área aproximadamente 2πr dr2\pi r\,dr, y no 2π dr2\pi\,dr.

Cuándo usar polares

  • Frontera circular, anillo o sector. ¿Conviene polar? Sí.
  • Integrando que depende sólo de rr, por ejemplo x2+y2x^2+y^2 o x2+y2\sqrt{x^2+y^2}. ¿Conviene polar? Sí.
  • Expresiones que se simplifican al escribir x=rcos⁡θx=r\cos\theta, y=rsin⁡θy=r\sin\theta. ¿Conviene polar? Sí.
  • Regiones rectangulares y polinomios generales en xx, yy. ¿Conviene polar? Mejor Cartesiano.
  • Asimetrías que rompen la invariancia angular. ¿Conviene polar? Mejor Cartesiano.

Idea-fuerza

Antes de plantear la integral, dibuja la región: si la frontera es un arco de circunferencia y el integrando se reduce a una función de rr o de θ\theta, el cambio a polares convierte la integral en un producto de integrales separadas.

Visualización

  • Figura 1 — Malla polar sobre el plano xyxy.

    Malla polar con líneas radiales y circulares

    Figura 1. Cada celda de la malla polar tiene área r Δr Δθr\,\Delta r\,\Delta\theta, de modo que el Jacobiano del cambio aparece como factor natural de la integral.

  • Figura 2 — Región entre dos curvas polares.

    Región entre <span class="katex"><span class="katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>r</mi><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">r=f_1(\theta)</annotation></semantics></math></span><span class="katex-html" aria-hidden="true"><span class="katex-base"><span class="katex-strut" style="height:0.4306em;"></span><span class="mord mathnormal" style="margin-right:0.0278em;">r</span><span class="mspace" style="margin-right:0.2778em;"></span><span class="mrel">=</span><span class="mspace" style="margin-right:0.2778em;"></span></span><span class="katex-base"><span class="katex-strut" style="height:1em;vertical-align:-0.25em;"></span><span class="mord"><span class="mord mathnormal" style="margin-right:0.1076em;">f</span><span class="msupsub"><span class="vlist-t vlist-t2"><span class="vlist-r"><span class="vlist" style="height:0.3011em;"><span style="top:-2.55em;margin-left:-0.1076em;margin-right:0.05em;"><span class="pstrut" style="height:2.7em;"></span><span class="katex-sizing reset-size6 size3 mtight"><span class="mord mtight">1</span></span></span></span><span class="vlist-s">​</span></span><span class="vlist-r"><span class="vlist" style="height:0.15em;"><span></span></span></span></span></span></span><span class="mopen">(</span><span class="mord mathnormal" style="margin-right:0.0278em;">θ</span><span class="mclose">)</span></span></span></span> y <span class="katex"><span class="katex-mathml"><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>r</mi><mo>=</mo><msub><mi>f</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">r=f_2(\theta)</annotation></semantics></math></span><span class="katex-html" aria-hidden="true"><span class="katex-base"><span class="katex-strut" style="height:0.4306em;"></span><span class="mord mathnormal" style="margin-right:0.0278em;">r</span><span class="mspace" style="margin-right:0.2778em;"></span><span class="mrel">=</span><span class="mspace" style="margin-right:0.2778em;"></span></span><span class="katex-base"><span class="katex-strut" style="height:1em;vertical-align:-0.25em;"></span><span class="mord"><span class="mord mathnormal" style="margin-right:0.1076em;">f</span><span class="msupsub"><span class="vlist-t vlist-t2"><span class="vlist-r"><span class="vlist" style="height:0.3011em;"><span style="top:-2.55em;margin-left:-0.1076em;margin-right:0.05em;"><span class="pstrut" style="height:2.7em;"></span><span class="katex-sizing reset-size6 size3 mtight"><span class="mord mtight">2</span></span></span></span><span class="vlist-s">​</span></span><span class="vlist-r"><span class="vlist" style="height:0.15em;"><span></span></span></span></span></span></span><span class="mopen">(</span><span class="mord mathnormal" style="margin-right:0.0278em;">θ</span><span class="mclose">)</span></span></span></span>

    Figura 2. Cuando la frontera interior y la exterior se expresan en polares, la integral doble se reduce a una integral iterada en rr y θ\theta con límites separados.

  • Figura 3 — Región polar tipo sector.

    Sector circular, caso más simple de integración polar

    Figura 3. El sector circular entre dos rayos y dos circunferencias es la región canónica sobre la que se monta toda integral polar.

Ejemplo Forma 1: Integral sobre un disco

Problema

Calcular

I=∬D(x2+y2) dA,I=\iint_D(x^2+y^2)\,dA,

donde DD es el disco de radio R=3R=3.

Resolución

Analítica

Con el cambio a polares, x2+y2=r2x^2+y^2=r^2 y la región es 0≤r≤30\le r\le3, 0≤θ≤2π0\le\theta\le2\pi. Así,

I=∫02π ⁣∫03r2⋅r dr dθ=∫02π ⁣∫03r3 dr dθ=∫02π[r44]03dθ=∫02π814 dθ=81π2.\begin{aligned} I&=\int_0^{2\pi}\!\int_0^3r^2\cdot r\,dr\,d\theta \\[0.45em] &=\int_0^{2\pi}\!\int_0^3r^3\,dr\,d\theta \\[0.45em] &=\int_0^{2\pi}\left[\frac{r^4}{4}\right]_0^3d\theta \\[0.45em] &=\int_0^{2\pi}\frac{81}{4}\,d\theta=\frac{81\pi}{2}. \end{aligned}

Explicación. La integral interior suma contribuciones radiales, y la exterior aprovecha la independencia de θ\theta para multiplicar por la longitud angular total.

Resultado. I=81π2≈127,235 u4I=\dfrac{81\pi}{2}\approx127{,}235\ \text{u}^4.

Diferenciador (Ejemplo Forma 1 · Analítica)
  • Dato. Se entrega una región circular y un integrando que depende sólo de x2+y2x^2+y^2.
  • Herramienta. Se usa x2+y2=r2x^2+y^2=r^2 y el Jacobiano dA=r dr dθdA=r\,dr\,d\theta.
  • Bifurcación. En el primer paso se reescribe la región como producto de un intervalo radial y otro angular.
Verificación (Ejemplo Forma 1 · Analítica)

Por simetría rotacional, el promedio de x2+y2x^2+y^2 sobre el disco es R2/2=9/2R^2/2=9/2. Multiplicado por el área πR2=9π\pi R^2=9\pi se obtiene R22πR2=πR42=81π2\frac{R^2}{2}\pi R^2=\frac{\pi R^4}{2}=\frac{81\pi}{2}.

En coordenadas cartesianas, ∫−33∫−9−x29−x2(x2+y2) dy dx\int_{-3}^{3}\int_{-\sqrt{9-x^2}}^{\sqrt{9-x^2}}(x^2+y^2)\,dy\,dx se reduce al mismo resultado tras evaluar las dos cuadraturas. La coincidencia confirma la corrección.

Errores comunes (Ejemplo Forma 1 · Analítica)
  • Olvidar el factor rr en dAdA y subestimar el resultado por un factor medio R/2R/2.
  • Restringir θ\theta a [0,π][0,\pi] cuando la región cubre todo el plano.
  • Tratar al integrando como función sólo de rr sin comprobar que la dependencia de θ\theta se anule.

Ejemplo Forma 2: Integral sobre un anillo

Problema

Calcular el área del anillo con radio interior a=1a=1 y radio exterior b=3b=3, y luego evaluar

I=∬anillo(x2+y2) dA.I=\iint_{\text{anillo}}(x^2+y^2)\,dA.

Resolución

Analítica

Paso 1 — área. El integrando es 11, por lo que la altura local vale rr:

A=∫02π ⁣∫13r dr dθ=∫02π[r22]13dθ=∫02π9−12 dθ=4⋅2π=8π.\begin{aligned} A&=\int_0^{2\pi}\!\int_1^3r\,dr\,d\theta \\[0.45em] &=\int_0^{2\pi}\left[\frac{r^2}{2}\right]_1^3d\theta \\[0.45em] &=\int_0^{2\pi}\frac{9-1}{2}\,d\theta=4\cdot2\pi=8\pi. \end{aligned}

Paso 2 — momento polar. Reescribiendo x2+y2=r2x^2+y^2=r^2,

I=∫02π ⁣∫13r3 dr dθ=∫02π[r44]13dθ=∫02π81−14 dθ=20⋅2π=40π.\begin{aligned} I&=\int_0^{2\pi}\!\int_1^3r^3\,dr\,d\theta \\[0.45em] &=\int_0^{2\pi}\left[\frac{r^4}{4}\right]_1^3d\theta \\[0.45em] &=\int_0^{2\pi}\frac{81-1}{4}\,d\theta=20\cdot2\pi=40\pi. \end{aligned}

Explicación. La única diferencia con la Forma 1 es el límite inferior de rr, que ya no es 00. El integrando se conserva como r3r^3 porque f dA=r2⋅r dr dθf\,dA=r^2\cdot r\,dr\,d\theta.

Resultado. A=8π u2A=8\pi\ \text{u}^2 y I=40π u4I=40\pi\ \text{u}^4.

Diferenciador (Ejemplo Forma 2 · Analítica)
  • Dato. La región es un anillo y se mantiene la dependencia radial del integrando.
  • Herramienta. Se conserva la descripción en polares con r∈[a,b]r\in[a,b].
  • Bifurcación. La integral de r2r^2 se reduce a una diferencia de potencias cuartas evaluada en los radios.
Verificación (Ejemplo Forma 2 · Analítica)

El área puede obtenerse por diferencia de discos: πb2−πa2=9π−π=8π\pi b^2-\pi a^2=9\pi-\pi=8\pi. El momento polar cumple la misma relación: π2(b4−a4)=π2(81−1)=40π\frac{\pi}{2}(b^4-a^4)=\frac{\pi}{2}(81-1)=40\pi. Ambos resultados coinciden.

Para validar la integral en cartesianas, basta notar que el anillo es el disco grande menos el disco pequeño, por lo que

I=81π2−π2=40π,I=\frac{81\pi}{2}-\frac{\pi}{2}=40\pi,

que reproduce el valor polar.

Errores comunes (Ejemplo Forma 2 · Analítica)
  • Integrar sobre [−3,3]2[-3,3]^2 y olvidar la condición r≥1r\ge1, recuperando el disco lleno.
  • Intercambiar aa y bb y obtener un resultado negativo.
  • Asumir densidad uniforme sin comentar que, si se busca la masa, debe multiplicarse por ρ(r,θ)\rho(r,\theta).

Cierre opcional: sector con integrando angular

La Forma 3 del temario integra un integrando que depende sólo del ángulo, por ejemplo f(θ)=θf(\theta)=\theta, sobre un sector. Se modela con la misma estructura dA=r dr dθdA=r\,dr\,d\theta y se reduce a una integral sobre θ\theta multiplicada por una constante radial. Aparece en distribuciones angulares de campo en el entrehierro de un alternador trifásico. En este apunte sólo se enuncia: el desarrollo numérico se aborda en la versión con programación.

Errores comunes transversales

  • Olvidar el Jacobiano rr. Consecuencia: resultado menor por un factor medio R/2R/2. Prevención: escribir siempre r dr dθr\,dr\,d\theta en polares.
  • Usar θ∈[0,π]\theta\in[0,\pi] cuando la región cubre ambos medios planos. Consecuencia: área y momento reducidos a la mitad. Prevención: identificar si la región cruza el eje xx.
  • Confundir disco con anillo. Consecuencia: diferencia de área πa2\pi a^2. Prevención: verificar si el origen pertenece a la región.
  • Aplicar polares a un integrando asimétrico. Consecuencia: integral más difícil que en cartesianas. Prevención: cambiar a polares sólo si la región o ff lo justifican.
  • Asumir densidad uniforme en aplicaciones de ingeniería. Consecuencia: se pierde la dependencia angular del campo real. Prevención: verificar si BB o ρ\rho dependen de θ\theta.