12. Teorema de Green (en el plano).

El siguiente teorema, llamado “Teorema de Green en el plano”, aplica para regiones planas limitadas por curvas cerradas y simples, regulares a trozos. Una idea intuitiva, en términos de "circulación", se puede ver en la sección 14..

Teorema 47 — Teorema de Green en el plano.

Sean P y Q campos escalares derivables con continuidad en un conjunto abierto S del plano XY. Sea C una curva simple cerrada regular a trozos y sea D la región encerrada por C (es decir, C = ∂D). Si D está contenida en S, se tiene la identidad

CPdx + Qdy = DQ ∂x P ∂y dA

donde C es recorrida en sentido contrario a las agujas del reloj.

Intuitivamente, C es recorrida en sentido contrario a las agujas del reloj si al caminar a lo largo de C la región D está siempre a la izquierda. Notar que C = ∂D indica que C es la frontera de D.

Área como una integral de línea. Si entonces Qx Py = 1, aplicando el teorema de Green (si se cumplen las condiciones) obtenemos otra manera para calcular el área de AD siendo la frontera de la región D una curva orientada contra-reloj. Por ejemplo si P(x,y) = 0 y Q(x,y) = x entonces

C0dx + xdy = D1dA = AD

Lo cual puede ser conveniente si la integral de línea no ofrece gran dificultad.

Teorema 48 — Área con el Teorema de Green..

Si D es una región plana limitada por una curva C, cerrada simple, regular a trozos y orientada contra-reloj. Si F = P,Q es un campo vectorial tal que ∂Q ∂x ∂P ∂y = 1 en D, entonces

AD = D1dA = CPdx + Qdy

En particular (para ciertas escogencias de F)

a
AD = Cxdy
b
AD = Cydx
c
AD = 1 2Cxdy ydx
d
Si P es un polígono simple de n vértices {(x0,y0),(x1,y1),,(xn1,yn1)}, con (xn,yn) = (x0,y0), entonces (fórmula de "shoelace" o "lazada")

A = 1 2 k=0n1(x kyk+1 xk+1yk).

Hay varias situaciones que hacen que la fórmula preferida sea

AD = 1 2Cxdy ydx

porque trata a las coordenadas x e y de manera equitativa. Las otras dos fórmulas dependen de proyectar la curva sobre un solo eje. Si la curva tiene comportamientos "extraños" o es muy alargada en una dirección, una de las fórmulas puede ser numéricamente menos estable o más propensa a errores si la parametrización no es uniforme.

La fórmula de ’lazada’ es el algoritmo base que utilizan los topógrafos en sus instrumentos para calcular el área de un terreno a partir de sus vértices

Propósito: Visualizar el teorema de Green para calcular área. El widget es matemáticamente honesto, en el recorrido "calcula" el área de la región encerrada entre el trozo de curva y los rayos al origen.

Ejemplo297Teorema de Green
Calcular Cy2dx + x2dy si C es la curva de la Figura 9.62.

Solución.

En este caso, P(x,y) = y2 y Q(x,y) = x2. Como se cumplen las condiciones del teorema de Green entonces,

Cy2dx + x2dy = D∂Q ∂x ∂P ∂y dA = 010x2 2x 2ydydx = 012x3 x4dx = 3 10

PIC

Fig. 9.62: C = C1 C2 C3

Ejemplo298Teorema de Green
Calcular I = C(x + y)dx + (3x + arctan y)dy si C es la curva de la Figura 9.63.

Solución. En este ejemplo,

P(x,y) = x + y

Q(x,y) = 3x + arctan (y).

Como se cumplen las condiciones del teorema de Green, entonces

I = D∂Q ∂x ∂P ∂y dA = 11x21(33x)/23 1dydx = 115 3x 2x2dx = 26 3

PIC

Fig. 9.63: C = C1 C2 C3

Ejemplo299Teorema de Green
Calcular
I = C(x +  arcsen x)dx + (2x + ln (y2 3))dy

si C es la curva de la Figura 9.64.


Solución. En este ejemplo,

P(x,y) = x +  arcsen x

Q(x,y) = 2x + ln (y2 3).

Como se cumplen las condiciones del teorema de Green podemos poner

PIC

Fig. 9.64: C = C1 C2

I = D∂Q ∂x ∂P ∂y dA = 221x244x2 2dydx = 226 3x2 2 dx = 16.

Ejemplo300Área Teorema de Green

Calcular el área de la región encerrada por la elipse x2 a2 + y2 b2 = 1.


Solución. Debemos calcular AD = Cxdy.

Parametrizamos la elipse con

r(t) = acos ti + bsen tj con t [0,2π[.

Esta parametrización orienta la elipse contra-reloj. En este caso

x(t) = acos t

y(t) = bsen t

dy = bcos tdt,

AD =Cxdy = 02πacos t bcos tdt = πab.

Ejemplo301Área de un polígono simple

Sea P un polígono simple de n vértices {(x0,y0),(x1,y1),,(xn1,yn1)}, con (xn,yn) = (x0,y0).

1
Verifique que su área es
A = k=0n1(xk+1 + xk)(yk+1 yk) 2 .
2
Verifique que la fórmula anterior es equivalente a la clásica de shoelace:
A = 1 2 k=0n1(x kyk+1 xk+1yk).


Solución.

1
Asumimos que C = C0 + C1 + + Cn1 está orientada en sentido antihorario.

Cada lado Ck (del vértice (xk,yk) al vértice (xk+1,yk+1)), se parametriza por

rk(t) =((xk+1 xk)t + xk,(yk+1 yk)t + yk),t [0,1].

Por el teorema de Green en el plano,

A = Cxdy = D1dA.

Entonces

A = k=0n1C kxdy = k=0n101[(xk+1 xk)t + xk](yk+1 yk)dt = k=0n1(yk+1 yk) [(xk+1 xk)t2 2 + xkt]01 = k=0n1(yk+1 yk) (xk+1xk 2 + xk)

Si el recorrido es horario, la fórmula entrega A.

2
La fórmula obtenida es equivalente a la clásica de shoelace:
A = 1 2 k=0n1(x kyk+1 xk+1yk).

En efecto, desarrollando el sumando de nuestra fórmula,

(xk+1 + xk)(yk+1 yk) 2 = xk+1yk+1 xk+1yk + xkyk+1 xkyk 2 ,

y al sumar sobre k los términos xk+1yk+1 y xkyk producen una suma telescópica que se anula porque (xn,yn) = (x0,y0):

k=0n1(x k+1yk+1 xkyk) = xnyn x0y0 = 0.

Por tanto,

k=0n1(xk+1 + xk)(yk+1 yk) 2 = 1 2 k=0n1(x kyk+1 xk+1yk).

13.
Ejercicios

13.1
Sea F un campo vectorial dado por
F(x,y) = (x + y)i (x2 + y2)j.

La curva C es la frontera del trapecio limitado por las curvas y = 0, x = 1, x = 3 y y = x como se muestra en la Figura 9.65.

PIC

Fig. 9.65:

1
Calcular la integral CF dr usando el teorema de Green.
64 3
2
Calcular la integral CF dr sin utilizar el teorema de Green.
C 1F dr +C 2F dr +C 3F dr +C 4F dr = 4 36 + 283 + 43 = 64 3

13.2
Considere el campo vectorial
F(x,y) = xi + (x + y2)j.

Calcular CF dr donde C = C1 + C2 tal y como se muestra en la Figura 9.66.

PIC

Fig. 9.66:

Como se cumplen las condiciones para aplicar el teorema de Green en el plano, excepto la orientación de la curva, entonces

CF dr = 223y2431 0dxdy = 8.

13.3
Considere el campo vectorial F(x,y,z) = yi + x2j

Calcule la integral de línea CF dr donde C es la curva que se muestra en la Figura 9.67.

PIC

Fig. 9.67:
Por el teorema de Green:

nCF dr = 010x2+1 (∂Q ∂x ∂P ∂y )dydx = 010x2+1 (2x 1)dydx = 01 (2x 1)(x2 + 1)dx = 1 6.

13.4
Calcule la integral CF dr si
F(x,y) = (xy2 + 2 + cos x)i + (yx2 ye sen y)j

C = C1 C2 C3,y (Figura 9.68).

PIC

Fig. 9.68:
Se omite.

13.5
Calcule C(4y + arctan (x/5)dx + (x2 + ln (y + 1))dy.
C es el camino representado en la Figura 9.69.

PIC

Fig. 9.69:
Se omite

13.6La finca “La Esperanza” tiene la forma poligonal que se muestra en la Figura 9.70. Sus vértices, medidos en kilómetros respecto a un origen en la esquina suroeste, son:

k 0 1 2 3 4 5 6 7 8
x k 2 8 15 18 17 14 9 5 1
yk 1 1 2 2 6 12 16 15 13 8

PIC

Fig. 9.70: Mapa

(a)
Estime el área de la finca contando cuadrados completos y fracciones de los cuadrados parciales. (Cada cuadrado de la malla mide 2 × 2 km.)
(b)
Calcule el área exacta aplicando el teorema de Green en la forma
A = 1 2 C(ydx + xdy) = 1 2 k=0n1(x kyk+1 xk+1yk),

con (xn,yn) = (x0,y0). Exprese el resultado en km2 y en hectáreas.

La curva C que recorre el contorno de la finca está formada por los segmentos rectos C0,C1,,C8 que unen los vértices Pk = (xk,yk) en sentido antihorario. Como cada lado es un segmento recto, la integral de línea se reduce a la suma de las contribuciones de cada lado, es decir, a la fórmula de shoelace:
A = 1 2| k=08(x kyk+1 xk+1yk)|,(x9,y9) = (x0,y0).

Calculamos cada término:

P0 P1 : 2 1 2 8 1 = 7, P1 P2 : 8 2 15 1 2 = 17 2 , P 2 P3 : 15 6 18 2 = 54, P 3 P4 : 18 12 17 6 = 114, P 4 P5 : 17 16 14 12 = 104, P 5 P6 : 14 15 9 16 = 66, P 6 P7 : 9 13 5 15 = 42, P 7 P8 : 5 8 1 13 = 27,

Sumando:

k=08(x kyk+1 xk+1yk) = 7 + 17 2 + 54 + 114 + 104 + 66 + 42 + 27 15 = 787 2 .

Por tanto,

A = 1 2 787 2 = 787 4 = 196.75km2.

Como 1km2 = 100 hectáreas,

A = 196.75km2 = 19675ha .

Una estimación por conteo de cuadrados debería dar un valor cercano a 196.75km2 (aproximadamente entre 180 y 210 km2 según la precisión en la aproximación de los cuadrados parciales).