CALV Guia2
CALV Guia2
CALV Guia2
Departamento de Matemáticas
ESCUELA
COLOMBIANA
DE INGENIERÍA Cálculo vectorial 2023 - 1
JULIO GARAVITO Guia de Trabajo II: Funciones Vectoriales
UNIVERSIDAD
1. Objetivos
Gráficar de funciones vectoriales y representar los vectores tangente unitario, normal unita-
rio y binormal en algunos puntos de una curva usando un asistente computacional. Utilizar
esta herramienta para obtener expresiones simplificadas para la torsión y la curvatura.
Revisar algunos elementos relativos a parametrización de regiones.
Para definir la función vectorial r(t) = hcos(t), sin (t), cos (2t)i, se puede ejecuta la siguiente
instrucción
r[t_]:={Cos[t],Sin[t],Cos[2*t]}
ParametricPlot3D[r[t],{t,0,2Pi},PlotStyle->{Thickness->0.02,Red}]
Graphics3D[{Black,Thick,Arrowheads[0.1],Arrow[{{0,0,0},r[2*Pi/3]}]}]
1
La siguiente instrucción permite ver el recorrido de la curva (animar), para esto se usa la
π π
función Manipulate; en ella la variable k satisface que 20 ≤ k ≤ 2π con incrementos de 20 .
El objeto que se manipula es la gráfica de la curva, en la cual t se define de 0 a k.
Manipulate[ParametricPlot3D[r[t],{t,0,k},PlotStyle->{Thickness->0.01,Red},PlotRange->{-1,1}],
{k,Pi/20,2*Pi,Pi/20}]
Manipulate[Show[{ParametricPlot3D[r[t],{t,0,k},PlotRange->{-1,1}],
Graphics3D[{Arrow[{{0,0,0},r[k]}]}]}],{k,Pi/20,2*Pi,Pi/20}]
T[t_]:=r’[t]/Norm[r’[t]]
T[t] T[Pi/3]
Es posible que el programa incluye unas funciones que dependen de la naturaleza del paráme-
tro. Para obviar estos valores absolutos, se puede usar Simplify y la función Assuming.
T[t_]:=Assuming[t>=0,Simplify[r’[t]/Norm[r’[t]]]]
T[t]
El gráfico del vector tangente unitario en t = 2π/3 se obtiene con la siguiente instrucción.
Graphics3D[{Blue,Thick,Arrow[{r[2*Pi/3],r[2*Pi/3]+T[2*Pi/3]}]}]
Con los siguientes comandos se obtiene la gráfica de r(t) junto con el vector posición y el
vector tangente en t = π/2.
Show[{ParametricPlot3D[r[t],{t,0,2*Pi},PlotStyle->{Thickness->0.01,Red},PlotRange->{-1,1}],
Graphics3D[{Arrow[{r[0],r[Pi/2]}]}],Graphics3D[{Blue,Thick,Arrow[{r[Pi/2],r[Pi/2]+T[Pi/2]}]}]}]
2
ContourPlot3D[{x^2-2*y^2==4*z,2*z==6-x^2-3*y^2},{x,-3,3},{y,-3,3},{z,-3,3},Mesh->False,
ContourStyle->Opacity[0.8]]
ParametricPlot3D[{2*Cos[t],Sqrt[3]*Sin[t],1/4(-1+5Cos[2t])},{t,0,2Pi},PlotStyle->{Red,
Thickness[0.015]}]
3. Parametrizaciones básicas
Presentamos algunas ejemplos sencillos sobre parametrización de curvas, superficies o sóli-
dos, pues, en los capı́tulos posteriores se hará uso de ellas.
Intercambiando a y b en esta fórmula se logra cambiar el punto inicial con el punto final. La
identidad cos2 t + sin2 t = 1 motiva el uso de la siguiente elección
También se puede asignar x = sin(t) y y = cos(t), de esta manera se describe la misma curva,
la diferencia está en los puntos terminales o en la orientación.
y y
1 1
x x
−1 1 −1 1
−1 −1
Figura 1: Gráficos de u(t) = hcos(t), sin(t)i y u(t) = hsin(t), cos(t)i
La ecuación x2 − y2 = 1 representa una hipérbola cuyas ramas abren hacia el eje x, la siguiente
idea es motivada por la relación cosh2 t − sinh2 t = 1, para todo t.
3
y y
1 1
x x
−2 −1 1 −2 −1 1
−1 −1
4
y
El cuadrilátero con vértices A(−2, −1),
B•4
B(−1, 4), C(3, 1), D (4, 3), se parametriza for- D
•
mando segmentos adecuadamente elegidos
(evitar crear un corbatı́n). Con las expresiones 2
P = (1 − t) A + tB y Q = (1 − t)C + tD, se de- •C
fine la función −
→r (t, u) = (1 − u) P + uQ, con
x
t ∈ [0, 1], u ∈ [0, 1]. En este caso, −→
r (t, u) = −2
A•
2 4
ht + 5u − 2, 5t + 2u − 3tu − 1i .
Con las siguientes instrucciones se consigue el gráfico
ParametricPlot[(1-u)*((1-t)*{-2,-1}+t*{-1,4})+u*((1-t)*{3,1}+t*{4,3}),{t,0,1},{u,0,1}]
ParametricPlot3D[{4Cos[\[Theta]]*Sin[\[Phi]],3*Sin[\[Theta]]*Sin[\[Phi]],
2*Cos[\[Phi]]},{\[Theta],0,2*Pi},{\[Phi],0,Pi},MeshStyle->Gray,PlotStyle->Orange]
5
Los gráficos se consiguen con las siguientes instrucciones.
a=ContourPlot3D[{2*z==x^2-y^2,2*x^2+y^2+2*y==8},{x,-5,5},{y,-5,5},{z,-9,5},Mesh->False,
ContourStyle->{Directive[Orange,Opacity[0.6]],Directive[Cyan,Opacity[0.2]]},PlotPoints->80,
Ticks->{{-3,0,3},{-3,0,3},{-8,-4,-2,2,4}}];
b=ParametricPlot3D[{3*u*Cos[t]/Sqrt[2],-1+3*u*Sin[t],27*u^2*Cos[t]^2/4+
3*u*Sin[t]-9*u^2/2-1/2},{t,0,2*Pi},{u,0,1},PlotStyle->{Cyan},PlotPoints->80];
Show[a,b]
ContourPlot3D[{z==9-(x^2+y^2)/4,x^2+(y-3)^2==4},{x,-6,6},{y,-6,6},{z,0,9},Mesh->False,
ContourStyle->{Orange,Cyan}]
ParametricPlot3D[{2*(2*u+1)*Sin[t],2*(2*u+1)Cos[t]-3*u+3,(u-1)*(12*(1+2u)*Cos[t]
-23-25u)/4},{u,0,1},{t,0,2*Pi}]
z
9 y
6
7
3 P•
5 •Q
x
3
-6 1 -6 −6 −3 3 6
-3 -3 −3
3 -1 3
6 6
x y −6
La base del cilindro se parametriza como x = 2 sin(t), y = 3 + 2 cos(t), t ∈ [0, 2π ]. La inter-
sección del paraboloide y el plano z = 0 se parametriza como x = 6 sin(t), y = 6 cos(t), con
t ∈ [0, 2π ].
Con el fin de parametrizar la proyección de σ sobre el plano xy, se pueden formar segmen-
tos entre dos puntos de la forma P (2 sin(t), 3 + 2 cos(t), 0) y Q (6 sin(t), 6 cos(t), 0) , es decir
(1 − u) P + uQ. Por tanto, la función que parametriza la región sobre el plano xy acotada por las
dos circunferencias presentadas antes es − →v (u, t) = (2 (2u + 1) sin t, 2 (2u + 1) cos t − 3u + 3, 0) .
La coordenada z sobre la superficie del paraboloide se consigue al remplazar
z = 14 (1 − u) (25u − 12(2u + 1) cos t + 23) . Por lo anterior, la superficie σ se parametriza, con
t ∈ [0, 2π ], u ∈ [0, 1], mediante la función
1
r(t, u) = h2(2u + 1) sin t, 2(2u + 1) cos t − 3u + 3, (1 − u)(25u − 12(2u + 1) cos t + 23)i
4
6
z
x
−2 2
la curva C,
La superficie sobre el paraboloide hiperbólico acotado por se puede parametri-
2 2
1 4u 4u
zar reemplazando x y z, para obtener el valor de y, y = 4 √ sin(t) − √ cos(t) =
5 3
4 2 32 2 2
5u − 15 u cos (t).Por lo tanto,
−
→ 4u 4 2 32 2 2 4u
r (t, u) = √ cos(t), u − u cos (t), √ sin(t)
3 5 15 5
con t ∈ [0, 2π ], u ∈ [0, 1], es una manera de parametrizar la mencionada superficie.
Para parametrizar la parte del paraboloide acotado por la curva C, el valor de y se consigue
2 2
4u 4u 32 2
sustituyendo en esta expresión y = 4 − √ cos(t) − √ sin(t) = 4 − 15 u cos2 (t) −
3 5
16 2
5u . Por tanto,
−
→ 4u 32 16 4u
r (t, u) = √ cos(t), 4 − u2 cos2 (t) − u2 , √ sin(t)
3 15 5 5
con t ∈ [0, 2π ], u ∈ [0, 1], representa la superficie deseada. Los gráficos se obtiene ejecutan-
do las siguientes instrucciones
a=Show[ContourPlot3D[{4*y==z^2-x^2,y==4-x^2-z^2},{x,-3,3},{y,-3,4},{z,-3,3},
Mesh->None,PlotPoints->90,ContourStyle->{{Gray,Orange},Opacity[0.6]}],
ParametricPlot3D[{4/Sqrt[3]*Cos[t],4-((4/Sqrt[5]*Sin[t])^2+(4/Sqrt[3]*Cos[t])^2),
4/Sqrt[5]*Sin[t]},{t,0,2*Pi},PlotStyle->{Purple,Thickness->0.015}]];
b=ParametricPlot3D[{4*u/Sqrt[3]*Cos[t],4-32/15*u^2*Cos[t]^2-16*u^2/5,
4*u/Sqrt[5]*Sin[t]},{t,0,2*Pi},{u,0,1},PlotPoints->90,Ticks->{{-2,0,2},{0,2,4},
{-1,0,1}},PlotStyle->RGBColor[0.9,0.1,0.1],MeshStyle->RGBColor[0.3,0.1,0.1]];
Show[a,b]
7
2
q con la esfera, el valor de z satisface que (1 + 3 cos(t)) +
Sobre la curva de intersección del cilindro
(2 sin(t)) 2 + z2 = 16, de donde z = ± 5 sin2 (t) − 6 cos(t) + 6, por tanto, una parametrización
q
−→ 2
para las curvas en mención es r (t) = 1 + 3 cos(t), 2 sin (t), ± 8 sin (t) − 6 cos(t) + 6 , t ∈
2
[0, 2π ] Para la superficie sobre la esfera dentro
q el cilindro, el valor de z está dado por (1 + 3u cos(t)) +
(2u sin(t)) 2 + z2 = 16, de donde, z = ± 15 − 6u cos(t) − 9u2 + 5u2 sin2 (t), por tanto,
q
2 2 2
r(u, t) = 1 + 3u cos(t), u sin (t), 15 − 6u cos(t) − 9u + 5u sin (t) ,
con t ∈ [0, 2π ], u ∈ [0, 1], es una parametrización de la parte de la esfera dentro del cilindro.
z
5
3P
•
1 -3
-3 -1
-1
3 1-11
x 5 3
5 y
-3
•
Q
Paraqrepresentar la parte del cilindro que está dentro de la esfera, denotando por
z1 = 6 − 6 cos(t) + 5 sin2 (t), podemos tomar dos puntos P = (1 + 3 cos(t), 2 sin(t), z1 ) y
Q = (1 + 3 cos(t), 2 sin (t), −z1 ), con ellos se parametriza el segmento rectilı́neo que los une.
Simplificando la expresión (1 − u) P + uQ, con u ∈ [0, 1], se obtiene la expresión la función
q
−→ 2
r (u, t) = 1 + 3 cos(t), 2 sin (t), (1 − 2u) 6 − 6 cos(t) + 5 sin (t) ,
con t ∈ [0, 2π ], u ∈ [0, 1], los gráficos se obtienen con ayuda de las siguientes instrucciones
Show[ContourPlot3D[{x^2+y^2+z^2==16,(x-1)^2/9+y^2/4==1},{x,-4,4},{y,-4,4},{z,-4,4},
ContourStyle->{Directive[Orange,Opacity[0.15]],Directive[Gray,Opacity[0.15]]},
Mesh->False],ParametricPlot3D[{{1+3*u*Cos[t],2*u*Sin[t],Sqrt[15-6*u*Cos[t]
-9*u^2+5*u^2*(Sin[t])^2]}},{u,0,1},{t,0,2Pi},PlotStyle->Red]
Show[ContourPlot3D[{x^2+y^2+z^2==16,(x-1)^2/9+y^2/4==1},{x,-4,4},{y,-4,4},{z,-4,4},
ContourStyle->{Directive[Orange,Opacity[0.15]],Directive[Gray,Opacity[0.15]]},
Mesh->False],ParametricPlot3D[{{1+3*Cos[t],2*Sin[t],(1-2*u)Sqrt[6-6*Cos[t]
+5*(Sin[t])^2]}},{u,0,1},{t,0,2Pi},PlotStyle->Red]]
8
Sólido acotado por dos paraboloides
Los paraboloides elı́pticos con ecuaciones z = 3 − 2x2 − y2 y z = x2 + 2y2 , determinan una
región que se parametriza hallando la curva de intersección que corresponde a x2 + y2 = 1. El
interior de esta región se representa con x = u cos v, y = u sin v con u ∈ [0, 1] y v ∈ [0, 2π ].
z
La coordenada z de dos puntos P y
2
z = 3 − 2x − y 2 3
Q
Q sobre los paraboloides, corresponde a
2
z1 = 21 (6 − 3u2 − u2 cos(2v)) y z2 = u2 (3 −
b
2
cos(2v)). El sólido se parametriza uniendo
2
1 los puntos P(u cos v, u sin v, u2 (3 − cos(2t))) y
z = x2 + 2y2
-1 Q(u cos v, u sin v, 21 (6 − 3u2 − u2 cos(2t))), me-
-1 b
1
P 1 diante segmentos rectilı́neos.
2 y
x
Es decir, la expresión (1 − w) P + wQ se simplifica y permite definir la función
u2
r(u, v, w) = u cos v, u sin v, 3(1 − w) + (6w − 3 − cos(2t)) ,
2
u ∈ [0, 1], v ∈ [0, 2π ], w ∈ [0, 1]. Una representación de esta regiónse obtiene con las siguientes
instrucciones.
RegionPlot3D[x^2+2*y^2<=z<=3-2*x^2-y^2&&x^2+y^2<=1,{x,-1,1},{y,-1,1},{z,0,3},PlotPoints->80,
Mesh->False,PlotStyle->Orange]
b) r(t) = hcos(2t), sin (2t), cos(4t)i d) r(t) = hcos2 (t), sin2 (t), sin (2t) + cos(t)i
c) r(t) = he−t cos(t), e−t sin(t), e−t sin(t)i e) r(t) = hcos(t), 2 sin (t), ti
2. Con las funciones vectoriales del numeral anterior, obtener una representación (dinámica)
gráfica de los vectores T, N y B, en puntos diferentes;
Sugerencia: Las funciones Manipulate o Animate permiten hacerlo.
3. Utilizar las funciones Dot y Cross que permiten evaluar el producto punto y el producto
cruz para hallar la curvatura y la torsión de las siguientes funciones vectoriales.
9
Parametrice y grafique en Mathematica las siguientes regiones.
x 2 y2 y2
1. La región del plano acotada por las elipses 16 + 9 = 1 y ( x − 1 )2 + 4 = 1, se proyecta
sobre el paraboloide hiperbólico 3z = y2 − x2 − 2y.
z z
3 3
-3 -3
6
3 y x
3
y
x -3 -3
2. La región del plano comprendida entre las 4. La porción del cilindro elı́ptico con ecua-
elipses x2
+ y2
= 1 y (x − 1 )2 + y2
= 1. ción 4x2 + 9y2 = 36 acotado por los planos
16 9 4
3z + x = 9 y 3z − y = −9.
y
x2 y2
3
16 + 9 =1
x
−4 −2 2 4
y2
−3 ( x − 1 )2 + 4 =1
10