0% encontró este documento útil (0 votos)
26 vistas30 páginas

Funciones Vectoriales y Campos en Matlab

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
26 vistas30 páginas

Funciones Vectoriales y Campos en Matlab

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Análisis II

Unidad Nº3: Funciones con


valores vectoriales
• Funciones vectoriales.

•Velocidad y aceleración. Longitud de arco.

• Campos vectoriales. Líneas de flujo.

•La divergencia y el rotacional.

•Aplicaciones con Matlab.

1
Funciones vectoriales

r
f :A⊆ R → R
n m

r r r r r r
x ∈ A ⇒ f ( x ) = ( f1 ( x ), f 2 ( x ),..., f m ( x ))

Es usual escribir los vectores en función de la


base canónica del espacio correspondiente. En
particular si m=3, los vectores unitarios son i, j,
k.
r
f :A⊆ R → R
n 3

r r r r r r
x ∈ A ⇒ f ( x ) = f1 ( x )i + f 2 ( x ) j + f 3 ( x )k
 Función vectorial: n = 1 y m ≥ 2
r
f : A ⊆ R → R2
r
f (t ) = (2 + t ,1 − 2t ); − 5 ≤ t ≤ 5
r
f (t ) = (2 + t )i + (1 − 2t )j
15

10

t 5

-5 5 0

-5

-10
-3 -2 -1 0 1 2 3 4 5 6 7
>> t=linspace(-5,5,1000);
plot((2+t),(1-2*t))
3
r
f : A ⊆ R → R2
r (
 t t 2 −1 2 t 2 −1
f (t ) =  2
) ( ) ;
, 2  −5 ≤ t ≤ 5
 t +1 t +1 
2
y 1.5

0.5

-0.5

-1

-1.5

-2
-5 -4 -3 -2 -1 0 1 2 3 4 5

x
>> t=linspace(-5,5,1000);
plot((t.*(t.^2-1))./(t.^2+1),(2*(t.^2-1))./(t.^2+1))

4
 Función vectorial: n = 1 y m = 3
r
f :A⊆ R→ R 3

r
f (t ) = (sen t , cos t , t ); 0 ≤ t ≤ 8π

z 30

25

20

15
t
10

0 8ᴨ 5

0
1
0.5 1
0 0.5
y -0.5
0
x
>> t=linspace(0,(8*pi),2000); -1
-0.5
-1
plot3(sin(t),cos(t),t),grid on
5
 Campo vectorial: n ≥ 2 y m ≥ 2
r
f : A ⊆ R2 → R2
r
f ( x, y ) = (− y, x )
3

-1

-2

-3
-3 -2 -1 0 1 2 3

6
Consideraciones sobre funciones vectoriales

• Una función vectorial es continua en una región D si


cada función componente es continua en D.
• Una función vectorial es derivable en una región D si
cada función componente es derivable en D.

Ejemplo: Para cada casa determinar el conjunto abierto


para el cual la función es continua
r  1 1 
a ) f ( x, y ) =  , 
 x+ y x− y
r
(
b ) f ( x, y ) = x + y , x − y
r
)
c) f (t ) = t ln(t + 1)
 Función vectorial de variable real
r
f :A⊆ R→ R , m≥2
m

En general se trabaja con m=2 o m=3, que son las


dimensiones de los espacios en que tienen lugar las
aplicaciones que interesa.
r
Definición: Si f = ( f1 , f 2 ,..., f m ) es una función
vectorial, definimos el límite, la derivada y la integral por

t →a
r
( t →a t →a t →a
)
lim f (t ) = lim f1 (t ), lim f 2 (t ),..., lim f m (t ) ,
r
f ′(t ) = ( f1′(t ), f 2′(t ),..., f m′ (t ) ),
br
 b b b

∫ a ∫
 a ∫ a ∫
f (t )dt =  f1 (t )dt , f 2 (t )dt ,..., f m (t )dt 
a 
siempre que los componentes de los segundos
miembros tengan sentido. 8
 Algunas propiedades de funciones escalares que son
válidos para funciones vectoriales se resumen en el
siguiente teorema:
r r
Si f , g y ur son r
r derivablesr ren un intervalo, lo mismo
ocurre con f + g , u. f y f .g , tenemos

( ) ( )
r r′ r r r′ r r
f + g = f ′ + g ′ , u. f = u′. f + u. f ′ y

( )
r r′ r r r r
f .g = f ′.g + f .g ′ .
r r
Si f y g tienen los valores en R3, también tenemos

( )
r r′ r r r r
f × g = f ′ × g + f × g′ .

( )

r r r
Si f [u (t )] entonces f [u (t )] = u ′(t ) f ′[u (t )] siempre
que las funciones sean derivables.
9
Velocidad y aceleración
Consideremos funciones vectoriales de variable real en el
espacio de 2 o 3 dimensiones. Este tipo de funciones son
modelos matemáticos utilizados para describir movimiento de
una partícula.
Suponiendo que las componentes son funciones con segundas
derivadas, se definen los siguientes objetos:

Vector posición r (t ) = ( x(t ), y (t ), z (t ) )



r

v (t ) = r ′(t ) = ( x′(t ), y′(t ), z ′(t ) )


r r
• Vector velocidad

r r
• Velocidad v = v (t ) = r ′(t )

• Vector aceleración a (t ) = v′(t ) = r ′′(t ) = ( x′′(t ), y ′′(t ), z ′′(t ) )


r r r
10
Aspectos geométricos

Si una función vectorial es derivable y de longitud


constante en un abierto I, entonces su derivada es un
vector perpendicular a la función en cada t en I.
r r
r (t ). v (t ) = 0 11
Ejemplos: Para cada trayectoria, determinar el vector
velocidad, la velocidad y el vector aceleración. Evaluar en to
indicado. Esbozar el gráfico correspondiente.
a) r (t ) = (cos t , sent ), t o = π / 4
r
r
b) r (t ) = 2t 2 , t , t o = 1 ( )
a)
1.5
b)
2
r
r v
v 1.5
1

1
r r
r 0.5 r a
0.5
r
0
r
0 a -0.5

-1
-0.5

-1.5

-2
-1 0 1 2 3 4 5 6 7 8
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1

12
 ( ) ′ ( )
  =   + ′   −  =  ( ) + ′ ( )  − 
 ( ) ′ ( )
Longitud de arco

• rLa longitud de una curva dada paramétricamente por


r (t ) = ( x(t ), y (t ), z (t ) ) para t en [a,b] está dada por
……………………
b
r
longitud =
∫ a
r ′(t ) dt
b
r

longitud = v (t ) dt
a
siempre que v (t ) = r ′(t ) = ( x′(t ), y′(t ), z ′(t ) ) sea
r r
continua en [a,b].

•En cualquier intervalo de tiempo [t1,t2] la partícula se


mueve una distancia a lo largo de su trayectoria dada
por t2
r
s=

t1
v (t ) dt
14
Componentes del vector aceleración. Vector
tangente unitario y vector normal
El vector aceleración se puede descomponer en dos vectores
perpendiculares, uno paralelo al vector velocidad, por lo tanto
tangente a la curva, y otro ortogonal a éste.
Se definen así los siguientes vectores:
r r
v (t ) r
• Vector tangente unitario T (t ) = r siempre que v (t ) ≠ 0
v (t )
r
r T ′(t ) r
• Vector normal principal N (t ) = r siempre que T ′(t ) ≠ 0
r T ′(t )
r T ′(t ) r
• Vector curvatura κ (t ) = N (t )
v (t ) r
T ′(t )
• Curvatura de la curva en t κ (t ) =
v (t )
15
•Con estos objetos se tiene que el vector aceleración se puede
escribir r
r
v (t ) = v(t ).T (t )
r r r
a (t ) = v′(t ).T (t ) + v(t ).T ′(t )
r r r r
a (t ) = v′(t ).T (t ) + v(t ). T ′(t ) .N (t )
r r r
a (t ) = v′(t ).T (t ) + v (t ).κ (t ).N (t )
2

Componente Componente
tangencial a T centrípeta o
normal a N
r r r
a (t ) = aT .T (t ) + a N .N (t )

16
17
Ejemplos: Para cada trayectoria, determinar los vectores
tangente unitario, normal principal y curvatura en to. Esbozar
el gráfico correspondiente.
a) r (t ) = (cos t , sent ), t o = π / 4
r
r
(
b) r (t ) = 2t 2 , t , t o = 1 )
1.5
a) b)
2
r r r
v =T r v r
1
1.5
T r aT T
1
a
r r
r r
0.5
0.5 r N
r aN N
0
y
r r
a=N
0
-0.5

-1

-0.5
-1.5

-2
0 1 2 3 4 5 6
x
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1

18
Campos vectoriales. Líneas de flujo

Si F es un campo vectorial, una línea de flujo para F es


una trayectoria c(t) tal que:
r´(t) = F(r(t)) (1)
Es decir, F da el campo de velocidades de la trayectoria
r(t).

Una línea de flujo es la trayectoria seguida por una


pequeña partícula suspendida en el fluido. Las líneas de
flujo también se llaman apropiadamente líneas de
corriente o curvas integrales

Geométricamente, una línea de flujo para un campo


vectorial dado F es una curva trazada sobre el dominio de
F de manera que el vector tangente a la curva en cada
punto coincide con el campo vectorial.
19
¿Cómo determinar líneas de flujo a partir del campo
vectorial?

Si r(t) = ( x(t) , y(t) , z(t) ) entonces r´(t) = (x´(t), y´(t), z´(t))


y F(x , y, z) = ( P(x, y, z) , Q (x, y, z) , R(x , y , z) ) ,
de acuerdo a la expresión (1) tenemos que:
( x´(t) , y´(t) , z´(t) ) = ( P(x, y, z) , Q (x, y, z) , R(x , y , z) ).
Igualando coordenadas tenemos un sistema de ecuaciones
diferenciales:  x´(t ) = P ( x(t ), y (t ), z (t ))

 y´(t ) = Q( x(t ), y (t ), z (t ))
 z´(t ) = R( x(t ), y (t ), z (t ))

En muchos casos es imposible encontrar fórmulas explicitas
para las líneas de flujo, lo que nos lleva a recurrir a métodos
numéricos.
21
Ejemplo de campos vectoriales r
y líneas de flujo f ( x, y ) = (− y, x )
6 r
r ′(t ) = f (x(t), y(t))
r
4 (x′(t), y′(t)) = (− y(t), x(t))
x′(t) = − y(t)
2 
 y′(t) = x(t )
0

x(t) = r cost

-2
 y(t ) = r sent

-4

-6
-6 -4 -2 0 2 4 6
La divergencia y el rotacional de un campo vectorial
r
Sea f ( x, y, z ) = f1 ( x, y, z ) i + f 2 ( x, y, z ) j + f 3 ( x, y, z )k

La divergencia de este campo vectorial es el campo escalar


dado por: r ∂f ∂f ∂f 3
div f = 1
+ 2
+ .
∂x ∂y ∂z
Si se piensa en el flujo de un fluido, la divergencia permite
medir la expansión, o divergencia, del fluido desde un punto.
El rotacional de este campo vectorial es el campo vectorial
definido por
r  ∂f ∂f   ∂f ∂f   ∂f ∂f 
rot f =  3 − 2  i +  1 − 3  j +  2 − 1  k
 ∂y ∂z   ∂z ∂x   ∂x ∂y  .
El rotacional del campo de velocidad mide qué tanto se
arremolina o gira el fluido en cierta dirección.
La divergencia y el rotacional de un campo vectorial pueden
escribirse en términos de las operaciones vectoriales producto
escalar y producto vectorial. ∂ ∂ ∂
Se define el operador nabla como: ∇ = i + j + k
∂x ∂y ∂z
.
Tratando a este operador como un vector es posible interpretar
a la divergencia de un campo vectorial como
r ∂ ∂ ∂ 
∇. f =  i + j + k ( f1 ( x, y, z )i + f 2 ( x, y, z ) j + f 3 ( x, y, z )k )
 ∂x ∂y ∂z 
r ∂f ∂f ∂f
∇. f = 1 + 2 + 3
∂x ∂y ∂z
El rotacional del campo vectorial se puede escribir como el
producto vectorial del operador con la función que define el
campo como sigue:
i j k
r r ∂ ∂ ∂
∇ × f = rot f =
∂x ∂y ∂z
f1 f 2 f 3
24
Ejemplo
Calcular la divergencia de  = −, −
 
 = ,−
 +  +

Ejemplo: calcular la
divergencia de la función
dada
Ejemplo
Calcular la divergencia de  = , 0
Ejemplo
Calcular la divergencia de  = , 
Ejemplo
Comprobar que un campo gradiente es irrotacional.

Comprobar que  , ,  = , −, 0 no es un campo


gradiente.

Ejemplo
Hallar el rotacional del campo vectorial  = (, − , 1)

Videos recomendados
[Link]
[Link]

También podría gustarte