13  2024-12-20
Examen de Análisis III

Fecha de publicación

20 de diciembre de 2024

Ejercicio 13.1 La temperatura de una placa de inducción térmica en cada punto (x,y) de su superficie viene dada por la función T(x,y)=yex22y2.

  1. ¿Cuál es la tasa de variación instantánea de la temperatura en el punto (0,0) si empezamos a movernos en la dirección en la que x decrece un tercio de lo que aumenta y?

  2. ¿En qué puntos la temperatura será máxima y mínima? ¿Cuál será la temperatura en esos puntos?

  3. Usar el polinomio de Taylor de segundo grado en el punto (1,0) para aproximar la temperatura en el punto (0.9,0.1).

  1. La tasa de variación instantánea la da la derivada direccional de T en (0,0) siguiendo la dirección del vector v=(13,1).

    Calculamos primero el gradiente de la función.

    Tx=yex22y2(x22y2)x=yex22y2(2x)Ty=ex22y2+y(ex22y2)y=ex22y2+yex22y2(x22y2)y=ex22y2+yex22y2(4y)=ex22y2(14y2)

    Así pues, el gradiente es

    T(x,y)=(ex22y2(2xy),ex22y2(14y2)),

    y en el punto (0,0) vale

    T(0,0)=(e02202(20),e02202(1402))=(0,1).

    Por tanto, la derivada direccional es

    Tv(0,0)=T(0,0)(13,1)|(13,1)|=(0,1)(13,1)(13)2+12=(0,1)(110,310)=310.

  2. Para encontrar los puntos críticos de la función, igualamos el gradiente al vector nulo y resolvemos el sistema de ecuaciones.

    Tx=02xyex22y2=0x=0 o y=0,Ty=0(14y2)ex22y2=0y=±12.

    Por tanto, los puntos críticos son (0,1/2) y (0,1/2).

    Para determinar si son máximos o mínimos, calculamos el determinante de la matriz hessiana en cada punto crítico.

    2T(x,y)=(2Tx22Txy2Tyx2Ty2)=(x(2xyex22y2)x(14y2)ex22y2y(2xyex22y2)y(14y2)ex22y2)=(4x2yex22y22yex22y28xy2ex22y22xex22y28xy2ex22y22xex22y216y3ex22y212yex22y2).

    La matriz hessiana en el punto (0,1/2) es

    2T(0,1/2)=(40212e022(12)2212e022(12)280(12)2e022(12)220e022(12)280(12)2e022(12)220e022(12)216(12)3e022(12)21212e022(12)2)=(e1/2004e1/2).

    y su determinante vale H(0,1/2)=e1/2(4e1/2)=4e1, que es positivo, por lo que en el punto (0,1/2) hay un extremo relativo, y como 2Tx2(0,1/2)=e1/2<0, es un máximo relativo. La temperatura en este punto vale T(0,1/2)=1/2e022(12)2=1/2e1/20.3033.

    En el punto (0,1/2) la matriz hessiana es

    2T(0,1/2)=(40212e022(12)2212e022(12)280(12)2e022(12)220e022(12)280(12)2e022(12)220e022(12)216(12)3e022(12)21212e022(12)2)=(e1/2004e1/2).

    y su determinante vale H(0,1/2)=e1/24e1/2=4e1, que es positivo, por lo que en el punto (0,1/2) también hay un extremo relativo, y como 2Tx2(0,1/2)=e1/2>0, es un mínimo relativo. La temperatura en este punto vale T(0,1/2)=1/2e022(12)2=1/2e1/20.3033.

  3. La fórmula del polinomio de Taylor de segundo grado en el punto (1,0) es

    PT2(x,y)=T(1,0)+T(1,0)(x1,y)+12((x1,y)2T(1,0)(x1y)).

    Calculamos el gradiente en el punto (1,0)

    T(1,0)=(e12202(210),e12202(1402))=(0,e1),

    y la matriz hessiana en el punto (1,0)

    2T(1,0)=(4120e1220220e122028102e1220221e122028102e1220221e122021603e12202120e12202)=(02e12e10).

    Sustituyendo en la fórmula del polinomio de Taylor de segundo grado, obtenemos

    PT2(x,y)=0+(0,e1)(x1,y)+12((x1,y)(02e12e10)(x1y))=e1y2e1(x1)y.

    Y la aproximación de la temperatura en el punto (0.9,0.1) es

    T(0.9,0.1)PT2(0.9,0.1)=e10.12e1(0.91)0.1=0.1e1+0.02e10.0441455

Ejercicio 13.2 Un coche circula por un circuito elíptico cuya trayectoria viene dada por la función vectorial f(t)=(400cos(10t),100sen(10t)), donde t está dado en minutos y las coordenadas de f en metros.

  1. Calcular la rapidez del vehículo en el instante t=π.

  2. Calcular la curvatura de la trayectoria en ese instante.

  3. Calcular la componente tangencial del vector aceleración en ese instante.

  4. Calcular la componente normal del vector aceleración en ese instante.

  5. Suponiendo que los neumáticos no proporcionan ningún agarre (por ejemplo porque hay hielo en la carretera), ¿cuál es el mínimo ángulo que debería tener el peralte de la curva en este instante para que el coche no se salga del circuito? Tómese una aceleración debida a la gravedad de 9.8 m/s2.

  1. La rapidez del vehículo en el instante t=π es módulo del vector velocidad en ese instante, así que calculamos primero el vector velocidad.

    f(t)=(4000sin(10t),1000cos(10t))

    y en el instante t=π vale

    f(π)=(4000sin(10π),1000cos(10π))=(0,1000).

    por lo que la rapidez es |f(π)|=02+10002=1000 m/min.

  2. La curvatura de la trayectoria en el instante t viene dada por la fórmula

    κ(t)=|f(t)×f(t)||f(t)|3.

    Calculamos primero el vector aceleración.

    f(t)=(40000cos(10t),10000sin(10t))

    que en el instante t=π vale

    f(π)=(40000cos(10π),10000sin(10π))=(40000,0).

    Por tanto, la curvatura en el instante t=π es

    κ(π)=|(0,1000,0)×(40000,0,0)||(0,1000)|3=|(0,0,4107)|109=410710009=0.04.

  3. La componente tangencial del vector aceleración en el instante t=π es

    aT(π)=f(π)f(π)|f(π)|=(0,1000)(40000,0)1000=400000+100001000=0.

  4. La componente normal del vector aceleración en el instante t=π es

    aN(π)=|f(π)×f(π)||f(π)|=|(0,0,4107)|103=4107103=40000.

  5. La fuerza centrípeta que debe ejercer el peralte de la curva para que el coche no se salga del circuito es la misma pero opuesta a la fuerza centrífuga del movimiento del coche. Como la fuerza centrífuga es la componente normal del vector aceleración que hemos calculado antes, se tiene que la fuerza centrípeta es Fc=40000m kg m/min2=40000/3600 kgm/s2, donde m es la masa del coche. Por otro lado, la fuerza F que ejerce la carretera sobre el coche se descompone en la componente gravitatoria Fg=mg, que compensa la fuerza de la gravedad, y la componente normal a la componente gravitatoria, que es la fuerza centrípeta Fc. Si θ es el ángulo del peralte de la carretera, se cumple que Fg=Fcos(θ) y Fc=Fsin(θ). Así que la fuerza centrípeta es Fc=Fsin(θ)=40000m, y la fuerza gravitatoria es Fg=Fcos(θ)=mg. Resolviendo el sistema de ecuaciones, se tiene que el peralte de la curva en este instante para que el coche no se salga del circuito es

    θ=arctan(FcFg)=arctan(40000/3600mmg)=arctan(11.119.8)48.6.

Ejercicio 13.3 La ecuación (x2+y2+z2+R2r2)2=4R2(x2+y2) define la superficie de un toroide de radio mayor R y radio menor r.

  1. Para un toro de radio mayor 4 y radio menor 2, probar que la ecuación anterior define implícitamente a z como una función de x y y en un entorno del punto (4,0,2).

  2. Calcular el gradiente de z como función implícita de x y y en el punto (4,0) e interpretarlo.

  1. La ecuación del toro es

    (x2+y2+z2+4222)2=442(x2+y2)(x2+y2+z2+12)264(x2+y2)=0.

    Par demostrar que la ecuación anterior define a z como función implícita de x e y en un entorno del punto (4,0,2) debemos comprobar que se cumplen las condiciones del teorema de la función implícita.

    Sea F(x,y,z)=(x2+y2+z2+12)264(x2+y2). Entonces F(4,0,2)=(42+02+22+12)264(42+02)=0, de modo que el punto (4,0,2) pertenece a la superficie del toro.

    En segundo lugar, F es diferenciable en un entorno de (4,0,2), puesto que es una función polinómica y las derivadas parciales son

    Fx=2(x2+y2+z2+12)2x642x=4x(x2+y2+z220),Fy=2(x2+y2+z2+12)2y642y=4y(x2+y2+z220),Fz=2(x2+y2+z2+12)2z=4z(x2+y2+z2+12).

    Como todas las derivadas parciales son funciones polinómicas, son continuas en todo R3, y en particular en el punto (4,0,2).

    Finalmente, Fz(4,0,2)=42(42+02+22+12)=4232=2560, por lo que se cumplen las condiciones del teorema de la función implícita y la ecuación define implícitamente a z como función de x e y en un entorno del punto (4,0,2).

  2. Para calcular el gradiente de z como función implícita de x e y en el punto (4,0), tenemos que calcular las derivadas parciales de z con respecto a x e y. Para ello, primero calculamos las derivadas parciales de F en el punto (4,0,2).

    Fx(4,0,2)=44(42+02+2220)=0,Fy(4,0,2)=40(42+02+2220)=0,Fz(4,0,2)=42(42+02+22+12)=256.

    Por tanto, las derivadas parciales de z con respecto a x e y en el punto (4,0) valen

    zx=FxFz=0256=0,zy=FyFz=0256=0.

    Es decir, el punto (4,0), es un punto críttico de z y por tanto en este punto la función ni crece ni decrece en cualquier dirección.

Ejercicio 13.4 Una ventana como la de la figura de más abajo está formada por un rectángulo de base l y altura h y un semicírculo de diámetro l. Si el área de la ventana es fija, qué relación debe haber entre l y h para que el perímetro de la ventana sea mínimo. Calcular las dimensiones óptimas para un área de 1 m2. Usar el método de los multiplicadores de Lagrange.

Ventana

El área de la ventana es la suma de las áreas del rectángulo y le semicírculo, así que la función que da el área de la ventana es

A(l,h)=lh+π2(l2)2=lh+π8l2.

Por otro lado, el perímetro de la ventana es la suma de los tres lados del rectángulo y el perímetro del semicírculo, así que la función que da el perímetro de la ventana es

P(l,h)=l+2h+π2l.

Para minimizar el perímetro de la ventana podemos aplicar el método de los multiplicadores de Lagrange. Así que debemos resolver el sistema de ecuaciones

{A(l,h)=λP(l,h)A(l,h)=k

El gradiente de la función del área es

A(l,h)=(h+π4l,l)

y el gradiente de la función del perímetro es

P(l,h)=(1+π2,2).

Así que el sistema de ecuaciones a resolver es

{h+π4l=λ(1+π2)l=2λlh+π8l2=k

De la segunda ecuación obtenemos que l=2λ, sustituyendo en la primera ecuación obtenemos que

h+π42λ=λ(1+π2),

de donde se deduce que que h=λ, y por tanto la relación entre l y h es l=2h.

Finalmente, sustituyendo en la tercera ecuación, para un área de 1 m2, se tiene

2hh+π8(2h)2=1h2(2+π2)=1h=12+π20.5292 m,

y l=2h20.5292=1.0584 m.

Ejercicio 13.5 Calcular el centro de masas de una placa metálica semicircular de radio a en los siguientes casos:

  1. La densidad en cada punto es proporcional a la distancia al origen.
  2. La densidad en cada punto es proporcional a la distancia al eje x.
  1. La distancia al origen de un punto (x,y) es x2+y2, así que la densidad en cada punto es δ(x,y)=kx2+y2, donde k es una constante de proporcionalidad.

    La región de integración es el semicírculo de radio a, que en coordenadas rectangulares es

    R={(x,y)R2:x2+y2a2,y0}={(x,y)R2:axa,0ya2x2}.

    y en coordenadas polares es

    R={(r,θ)R+×R:0ra,0θπ}.

    Puesto que la región de integración es más simple en coordenadas polares, calculamos el centro de masas en estas coordenadas. Además la función de densidad se simplifica en coordenadas polares, puesto que x2+y2=r. Así que la densidad en cada punto es δ(r,θ)=kr. La masa de la placa es

    M=Rδ(r,θ)dA=0a0πkrrdθdr=0a0πkr2dθdr=k0ar2[θ]0πdr=kπ0ar2dr=kπ[r33]0a=kπa33.

    Por otro lado, los momentos respecto al eje x y al eje y son

    Mx=Ryδ(x,y)dA=0a0πrsen(θ)kr2dθdr=k0ar30πsen(θ)dθdr=k0ar3[cos(θ)]0πdr=k0ar3[(11)]dr=2k0ar3dr=2k[r44]0a=ka42.My=Rxδ(x,y)dA=0a0πrcos(θ)kr2dθdr=k0ar30πcos(θ)dθdr=k0ar3[sen(θ)]0πdr=k0ar3[sen(π)sen(0)]dr=k0ar3[00]dr=0.

    Así pues, el centroide tiene coordenadas

    x¯=MyM=0kπa3/3=0,y¯=MxM=ka4/2kπa3/3=3a2π.

  2. La distancia al eje x de un punto (x,y) es |y|, así que, en la región de integración, que cae en el semiplano y>0, la densidad en cada punto es δ(x,y)=ky, donde k es una constante de proporcionalidad.

    Al igual que antes trabajaremos en coordenadas polares ya que los cálculos son más simples. En coordenadas polares la función de densidad es δ(r,θ)=krsen(θ). La masa de la placa es

    M=Rδ(r,θ)dA=0a0πkrsen(θ)rdθdr=k0a0πr2sen(θ)dθdr=k0ar20πsen(θ)dθdr=k0ar2[cos(θ)]0πdr=k0ar2[(11)]dr=2k0ar2dr=2k[r33]0a=2ka33.

    Y los momentos respecto al eje x y al eje y son

    Mx=Ryδ(x,y)dA=0a0πrsen(θ)krsen(θ)rdθdr=k0ar30πsen2(θ)dθdr=k0ar30π1cos(2θ)2dθdr=k0ar3[θ2sen(2θ)4]0πdr=k0ar3(π20)dr=kπ20ar3dr=kπ2[r44]0a=kπa48.My=Rxδ(x,y)dA=0a0πrcos(θ)krsen(θ)rdθdr=k0ar30πsen(θ)cos(θ)dθdr=k0ar3[sen(θ)22]0πdr=k0ar3(sen(π)22sen(0)22)dr=k0ar3(00)dr=0.

    Así pues, el centroide tiene coordenadas

    x¯=MxM=kπa4/82ka3/3=3πa16y¯=MyM=02ka3/3=0.

Ejercicio 13.6 Calcular el volumen comprendido entre la gráfica de la función f(x,y)=xy y el plano z=0 en la región limitada por las curvas de la figura.

Región de integración

En primer lugar determinamos los puntos de corte de las dos gráficas para expresar la región de integración. Para ello, resolvemos el sistema de ecuaciones

{y=x3x=yy3

Despejando x en la primera ecuación obtenemos que x=y3, y sustituyendo en la segunda ecuación obtenemos que y3=yy3, de donde se deduce que y(12y2)=0, y por tanto tenemos tres soluciones y=0, y=1/2=21/2 y y=1/2=21/2. Sustituyendo en la primera ecuación obtenemos que x=0, x=23/2 y x=23/2 respectivamente. Por lo tanto, los puntos de corte son (0,0), (23/2,21/2) y (23/2,21/2).

Si integramos primero con respecto a y y luego con respecto a x, la región de integración resulta más compleja porque los límites de x no serían los puntos de corte de las dos gráficas, sino los extremos de la función x=yy3, pero además habría que descomponer la región de integración en cuatro partes. Por lo tanto, es más conveniente integrar primero con respecto a x y luego con respecto a y. En este caso, la región de integración puede descomponerse en las dos subregiones siguientes

R1={(x,y)R2:yy3xy3,21/2y0}R2={(x,y)R2:y3xyy3,0y21/2}

Como la función f(x,y)=xy es simétrica en el primer y tercer cuadrantes, para calcular el volumen encerrado entre la superficie y el plano z=0 en la región de integración, podemos calcular el volumen en la región R2 y multiplicarlo por 2. Así que el volumen es

Rf(x,y)dA=2021/2y3yy3xydxdy=2021/2[x22y]y3yy3dy=2021/2[(yy3)22yy62y]dy=2021/2[y22y4+y62yy52]dy=2021/2[y2y3+y52y52]dy=2021/2[y2y32]dy=2[y24y44]021/2=2(1421422)=2(18116)=2(116)=18.