Análisis Matemático
III
PROLOGO
Uno de los principales objetivos en la elaboración de este texto, es dotar como instrumento de consulta y guía en el estudio del Análisis III, a los estudiantes de Ingeniería de Mantenimiento de la Espoch, como también a otros que necesiten esta herramienta en su formación básica en las diferentes carreras dentro y fuera de nuestra institución. El ordenamiento de los contenidos del presente texto, corresponden al programa de análisis III de la escuela de Ingeniería de Mantenimiento, en los que van reflejados los muchos años de experiencia de dictar este asignatura, completando lógicamente con la consulta y recopilación de varios argumentos en los diferentes textos de matemática superior disponibles. Al ser esta una materia muy importante dentro de la formación básica de cualquier carrera de ingeniería, contiene la teoría y los distintos procedimientos analíticos para poder resolver los diferentes temas de cada uno de los capítulos, así también ejemplos resueltos y propuestos los mismos que servirán como una herramienta de consulta para los estudiantes de Ingeniería. Expreso mi deseo de que este esfuerzo, pueda contribuir al mejoramiento de la enseñanza de análisis III, así como como despertar en el lector lector el interés por profundizar los temas expuestos aquí, ya que son una herramienta herramienta importante en la solución de muchos problemas problemas
que aparecen aparecen con frecuencia en el estudio estudio de fenómenos
naturales y dentro de la vida cotidiana.
El autor
1
Análisis Matemático
III
CONTENIDO TEMAS: Página:
UNIDAD I 1. Conceptos básicos y terminología.
5
1.1 Definición.
5
1.2 Clasificación de las ecuaciones diferenciales.
5
1.3 Orden y grado de una ecuación diferencial.
6
1.4 Soluciones de las ecuaciones diferenciales ordinarias.
6
1.5 Origen de las ecuaciones diferenciales ordinarias.
8
1.6 Ejercicios propuestos. 13 UNIDAD II 2. Ecuaciones de primer orden.
16
2.1 Ecuaciones de Variables Separables.
16
2.2 Ecuaciones Homogéneas.
21
2.3 Ecuaciones Exactas o en diferenciales totales.
28
2.4 Ecuaciones Lineales.
36
2.5 Ecuaciones de 1er orden no resueltas con respecto a la Derivada. 2.6 Ecuaciones de Lagrange y de Clairaut.
42 49
UNIDAD III 3. Aplicaciones de las Ecuaciones Diferenciales de 1er orden.
52
3.1 Trayectorias Ortogonales.
53
3.2 Familias Isógona les.
55
3.3 Aplicaciones Geométricas.
56
3.4 Aplicaciones Físicas.
61
3.5 Ejercicios propuestos.
73
2
Análisis Matemático
III
CONTENIDO TEMAS: Página:
UNIDAD I 1. Conceptos básicos y terminología.
5
1.1 Definición.
5
1.2 Clasificación de las ecuaciones diferenciales.
5
1.3 Orden y grado de una ecuación diferencial.
6
1.4 Soluciones de las ecuaciones diferenciales ordinarias.
6
1.5 Origen de las ecuaciones diferenciales ordinarias.
8
1.6 Ejercicios propuestos. 13 UNIDAD II 2. Ecuaciones de primer orden.
16
2.1 Ecuaciones de Variables Separables.
16
2.2 Ecuaciones Homogéneas.
21
2.3 Ecuaciones Exactas o en diferenciales totales.
28
2.4 Ecuaciones Lineales.
36
2.5 Ecuaciones de 1er orden no resueltas con respecto a la Derivada. 2.6 Ecuaciones de Lagrange y de Clairaut.
42 49
UNIDAD III 3. Aplicaciones de las Ecuaciones Diferenciales de 1er orden.
52
3.1 Trayectorias Ortogonales.
53
3.2 Familias Isógona les.
55
3.3 Aplicaciones Geométricas.
56
3.4 Aplicaciones Físicas.
61
3.5 Ejercicios propuestos.
73
2
Análisis Matemático
III
UNIDAD IV 4. Ecuaciones Diferenciales de Orden Superior
77
4.1 Caso de Integración Inmediata
77
4.2 Ecuaciones que se reducen a un Orden Inferior
80
4.3 Ecuaciones Diferenciales Lineales de Orden Superior
84
4.4 Ecuaciones Diferenciales Lineales Homogéneas con Coeficientes constantes
87
4.5 Ecuaciones Diferenciales Lineales no Homogéneas con Coeficientes constantes
89
4.6 Ecuaciones de Euler
95
4.7 Ejercicios Propuestos
97
4.8 Aplicaciones de las Ecuaciones de 2do Orden
100
UNIDAD V 5. Sistema de Ecuaciones Diferenciales con Coeficientes Constantes
114
5.1 Reducción de un Sistema a una Ecuación de enésimo orden 116 5.2 Método Matricial (Euler)
119
5.3 Ejercicios Propuestos
123
UNIDAD VI 6. Series de Fourier
124
6.1 Introducción
124
6.2 Definiciones Preliminares
126
6.3 Integrales Requeridas para calcular los Coeficientes de La serie
129
6.4 Cálculo de los Coeficientes
130
6.5 Serie de Fourier de Funciones pares e impares
134
6.6 Series de Fourier de Medio Rango
140
6.7 Armónicas impares y pares
144
6.8 Suma de una serie de Fourier en un Punto de Discontinuidad 146
3
Análisis Matemático
III
6.9 Período General
147
6.10 Análisis Armónico Numérico
154
6.11 Ejercicios Propuestos
157
UNIDAD VII 7. La Transformada de Laplace
158
7.1 Trasformadas de Laplace de algunas funciones básicas
158
7.2 Demostración de algunas de ellas
159
7.3 Principales Propiedades de la Trasformada de Laplace
161
7.4 Aplicaciones de las Trasformadas de Laplace
165
7.5 Ejercicios Propuestos
175
8. Bibliografía
178
4
Análisis Matemático
III
UNIDAD I 1. Conceptos básicos y terminología. 1.1 Defi ni ción.- Una ecuación donde la incognita es la misma función y que depende de una o más variables independientes, de una o más variables dependientes y de sus derivadas, se denomina “ Ecuación Diferencial”.
Ejemplos:
1.
d 2 y dx
2
2 xy
dy dx
e y una variable independiente y una dependiente.
2 F 2 F 2. 0 x 2 y 2
3. 4.
u 2 u 0 x y
d 2 u dx 2
5. m 6. L
dos variables independientes y una dependiente.
2uv
d 2 x dt 2 d 2 q dt 2
du dv dx dy
una variable dependiente y 2 variables independientes
d 2 v dy 2
0
dos variables independientes y 2 dependientes
kx una variable independiente y una dependiente. R
dq dt
1 C
q 0 una variable independiente y una dependiente
1.2 Clasificación de las ecuacion es dif erenci ales. Las ecuaciones diferenciales se clasifican en dos tipos: La una que depende de una sola variable independiente a la que se lo denomina Ecuación
diferencial ordinaria o (EDO); y la otra que depende de dos o más variables independientes a la que se lo denomina Ecuación Diferencial en Derivadas parciales o (EDP). En los ejemplos anteriores tenemos: 1. EDO 2. EDP (Ecuación de Laplace) 3. EDP
4. EDP 5. EDO ( Ecuación diferencial del M.A.S) 6. EDO ( Ecuación diferencial de la corriente)
En este documento estudiaremos y analizaremos las ecuaciones diferenciales Ordinarias EDO de primer orden, es decir ecuaciones de la forma: 5
Análisis Matemático F (x, y, y´) = 0
III
(I)
1.3 Orden y grado de un a ecuación di ferenci al. Orden: El orden de una ecuación diferencial nos da la derivada de mayor orden. Grado: El grado de la ecuación es el exponente a que esta elevada la derivada de mayor orden
después de realizar cualquier análisis algebraico. Así: En los ejemplos anteriores tenemos: 1. 2do orden, primer grado.
4. 2do orden, primer grado.
2. 2do orden, primer grado
5. 2do orden , primer grado
3. 1er orden, segundo grado
6. 2do grado, primer orden.
1.4 Soluciones de una ecuación dif erencial ordinar ia. Sea F x, y, y ,...... y n 0 (I) una EDO sus posibles soluciones serán: a.- Una función y = f(x) siempre que está, transforme en identidad la ecuación (I). En este caso la Solución viene dada por: y f ( x) C donde C es una constante arbitraria que se llama familia de un parámetro.
b.- A menudo, cuando resolvemos una EDO las soluciones podemos expresarles en las formas: g(x, y) = 0 y x = x (t), y = y (t) a estas expresiones se denominan soluciones implícitas y soluciones paramétricas respectivamente.
Las curvas que describen estas ecuaciones se denominan curvas integrales y representan las gráficas de las soluciones, las mismas que pueden ser generales o particulares dependiendo si existen o no condiciones iniciales o de frontera. Ejemplo 5. Ej emplo 1.
Verificar que la función y = c1 cos x + c2 sen x, es solución de la ecuación y´´ + y = 0 (A) Comprobación: Derivamos dos veces: y c1 senx c2 cos x
y = c1 cos x + c2 sen x.
y
y c1 cos x c2 sen. x , reemplazando el valor de y e y´´
en (A) comprobamos que la expresión dada si es solución. Ej emplo 2.
Dada la relación: x2 + y2 = R2 dy dx
(1). , comprobar si es Solución de la ecuación
x y
6
Análisis Matemático
III
Solución: Derivando implícitamente (1) tendremos: 2 x 2 y
dy dx
0
dy dx
x
; Reemplazando en la y
ecuación dada comprobamos que se transforma en identidad, por lo que sí es solución. Las soluciones en este caso también se podrían presentar así: y R 2 x 2 , y R 2 X 2
Ej emplo 3. Dada las siguientes ecuaciones paramétricas:
x t 2 2t 2 2 3 2 y 3 t t c
Comprobar que son Solución de la ecuación diferencial x y 2 2 y 2 . Desarrollo: dy dy dx
2t 2 2t 2t t 1 dt t dt dx dx 2t 2 2t 1
dy dt
dt
Reemplazando tenemos: t 2 2t 2 t 2 2t 2 : por lo que al llegar a tener una identidad, se concluye que si es Solución de la ecuación. Ej emplo 4. Verificar que la función y ( x) e x
2
x
0
e t dt e x , es Solución de la ecuación 2
2
diferencial y 2xy 1. Desarrollo: Derivando y ( x) e x
2
x
0
e t dt e x y ( x) 2 xe x 2
2
2
e t dt 1 2 xe x . 2
2
Reemplazando en la ecuación dada tenemos: y 2 xy 2 xe x
2
y 2 xy 2 xe x
2
x
0 x
0
e t dt 1 2 xe x 2 x e x 2
2
x
2
2
2
e t dt 1 2 xe x 2 xe x
2
0
x
0
2 2 e t dt e x
e t dt 2 xe x 1, y 2 xy 1. /// 2
2
Lo cual queda demostrado que si es solución.
7
Análisis Matemático
III
Ej emplo 5.
Si y c 2 cx (A) es la primitiva de la ecuación diferencial 4 x y´2 2 xy´ y 0 (B); Hallar la solución para la condición inicial y (1) = 2
Solución: Derivamos (A): 2 y c y´ c y´
c 2 y c
(C)
Ponemos (C) en (B) 4 x
c2 4 y c
2 x
c 2 y c
y 0
c 2 x cx y c y y c
0 c 2 5c 4 0 c 1
y c y c = 4, por lo que las soluciones particulares son: y 12 x y 2
2
y 42 4 x
Nota: Las constantes arbitrarias dependen del orden de la ecuación diferencial: Si es de primer orden tendremos una sola constante, si es de segundo dos, y así sucesivamente.
1.5 Or igen de las ecuaci ones difer encial es. Las ecuaciones diferenciales aparecen no solo a partir de las familias de curvas geométricas, sino también del intento de describir una variedad de problemas sobre geometría, física, economía, ingeniería, ciencias sociales, etc. 1.5.1 Ecuación dif er encial de un a fami li a de cur vas.
Si se tiene la ecuación de una familia de curvas, se puede obtener su ecuación diferencial eliminando las constantes o parámetros aislándolas a un solo miembro y derivando. También se puede eliminar la constante derivando la ecuación dada, tantas veces como constantes arbitrarias tenga, y se resuelve el sistema formado con la ecuación original. Ej emplo 1.
Encontrar la ecuación diferencial cuya solución general es: y x 2 C 1e x C 2 e 2 x Solución: y x 2 C 1e x C 2 e 2 x (1.)
Derivamos (1) y 2 x C 1e x 2C 2 e 2 x (2.) (1.)
+ (2.)(-1), 8
Análisis Matemático
III
y y x 2 2 x 3C 2 e 2 x (3.)
Derivando (3.)
y y 2 x 2 6C 2 e 2 x (4.)
(3.)*2 + (4.) 2 y y y 2 x 2 2 x 2 y y 2 y 21 x x 2 . /// Esta es la ecuación diferencial
pedida. Ej emplo 2.
Encontrar la ecuación diferencial cuya solución general es: y C 1 x C 2 e x Desarrollo: y C 1 x C 2 e x (1.)
Derivamos (1.) y C 1 C 2 e x (2.) (1.)
+ (2.).
y y C 1 x 1 (3.)
Derivando (3.)
y y C 1 (4.)
(3.) en (4.) y 0. /// Esta es la ecuación diferencial pedida y y y y x 1 y x 1 y x
Ej emplo 3.
Encontrar la ecuación diferencial cuya solución general es la familia de circunferencias
x a 2 y b 2 r 2 en el plano xy, siendo a, b y r constantes arbitrarias. Solución:
x a 2 y b 2 r 2 (1). Derivamos (1).
2 x a 2 y b y 2r x a y y b y r (2)
Volviendo a derivar (2) 1 y y y 2 b y 0
y 2 y y 1 y
b (3)
Derivando (3) y simplificando tenemos: y y y y y 2 y y 1 y y 2 y y
y 2
0
y y 2 y 2 y y y y y 0 3 y y 2 y y 2 1. /// 2 y y 2 y y y
9
Análisis Matemático
III
Ej emplo 4.
Se quiere definir una curva dada por la condición de que en cada uno de sus puntos (x, y), su pendiente
dy dx
sea igual al doble de la suma de las coordenadas del punto. Y
RT. P(x, y) 0
X
Figura 1.1 La ecuación diferencial que representa esta condición es:
dy dx
2 x y
1.5.2 Ecuaci ones dif erenci ales de probl emas físicos.
Las Ecuaciones diferenciales de problemas físicos provienen de diferentes campos tales como: la mecánica, eléctrica, química, etc. Ej emplo 1.
Una partícula de masa m se mueve a lo largo de una línea recta (el eje x) estando sujeta: 1) una fuerza proporcional a su desplazamiento x desde un punto fijo O en su trayectoria y dirigida hacia O y 2) Una fuerza resistente proporcional a su velocidad. Expresar la fuerza total como una ecuación diferencial.
Datos: Fuerza uno = - kx Fuerza dos = -k 1v, pero como v = dx /dt tenemos que F2 = -k 1 dx / dt. Según la segunda ley de Newton se tiene: Que la suma de las fuerzas en la dirección del movimiento = (masa)*(aceleración). Como la aceleración es la segunda derivada del espacio respecto al tiempo tenemos:: kx k 1
dx dt
m
d 2 x 2
dt
m
d 2 x 2
dt
kx k 1
dx dt
donde k y k 1 son factores de proporcionalidad,
esta es la ecuación diferencial pedida. 10
Análisis Matemático
III
Ej emplo 2.
Una lancha que pesa W kg se desliza por un plano inclinado a °. Si la fuerza de rozamiento que se opone al movimiento es Fr kg y la resistencia del aire expresado en kilogramos equivale a 0,05 veces la velocidad en centímetros por segundo, hallar la ecuación del movimiento.
Figura 1.2 Solución: Del gráfico tomamos los siguientes datos: F = Componente del peso en la dirección del movimiento = Wsen Fr = Fuerza de rozamiento Fa = Resistencia del aire = 0.05v, De acuerdo a la segunda ley de Newton, y recordando que la aceleración es la primera derivada de la velocidad respecto al tiempo tenemos: F F r F a ma m
dv dt
Wsen F r 0.05v
dv dt
0.05v m
Wsen F r m
///.que
es
la
ecuación pedida. Ej emplo 3
Considérese un circuito simple conectado en serie como el que se muestra en el gráfico, y que consta de un inductor, un resistor y un capacitor. Queremos encontrar la ecuación diferencial que nos permita determinar la carga q (t).
11
Análisis Matemático
III
Figura 1.3 Para poder encontrar la solución al problema, tenemos que recurrir a la segunda ley de Kirchoff que dice: La suma de las caídas de voltaje a través de cada uno de los componentes del circuito es igual a la tensión E (t) aplicada. Como q (t) es la carga del capacitor en un instante cualquiera, entonces la corriente i (t) está dado por i (t) = dq / dt, además sabemos que las caídas de voltaje en cada uno de los componentes son: En un inductor = L
di dt
L
d 2 q dt 2
.
En un capacitor = (i / C)*q. En un resistor = i * R R
dq dt
.
Con estos datos y según la segunda ley de Kirchoff tenemos L
d 2 q dt 2
R
dq dt
1 C
q E (t ) , que es la ecuación que nos permite calcular la carga q (t), en esta
ecuación L, C y R son constantes, llamadas inductancia, capacitancia, y resistencia respectivamente.
1.6 Ej ercicios propuestos . Determinar el orden y grado de las siguientes ecuaciones ordinarias. 1. 2.
d 2 q dt 2
R
dq dt
q C
0 2
dy * y0 dx 2 dx dx
d 2 y
dy
12
Análisis Matemático 3. x
3
dy dx
x
d 2 y dx 2
III y
4
3
d 3 y dx 3
5
d 3 y d 2 y 4. 3 2 y 0 dx dx
5. cos x * y 2 senx * y 4 1 6. x * y 4 y 3 y 0 7. x 2 y xy y 6 y 2 0 8. y 2 y x y 2 y 3 Comprobar, si son soluciones de las ecuaciones diferenciales dadas, las expresiones que se indican a continuación. 1. y´´ x 2 y 2 ;
y
1 x
.
2. x y dx xdy 0 ; y
c 2 x 2
.
2 x
3. y y 0 ; y 3 senx 4 cos x 4. y l 1 l 2 y l 1l 2 y 0 ; y c1e l x c2 e l x . 1
2
5. x 2 y y 2 x y ; x 2 xy y 2 c 2 6. 2 x 3 y y y 2 3x 2 ; 1 x y 2 x 3 7. Verificar que la función y = arcsen (x,y), satisface a la ecuación diferencial x y y y 1 x 2 y 2
8. Comprobar que la función
x
x y sen.t 2 dt , satisface a la ecuación diferencial 0
y x y y 2 sen. x 2
9. Probar que la función
y
2
0
log sen 2 x 2 cos 2 d , satisface a la ecuación
x 1 diferencial 1 x 2 y 1 x y y log . 2
Hallar, para las familias de curvas que se dan, las líneas que satisfagan las condiciones iniciales indicadas.
13
Análisis Matemático
III
1. x 2 y 2 c ; y(0) = 5 2. y c1 sen x c2 ; y() = 1; y´() = 0 3. y c1e x c2 e x c3 e 2 x ; y(0) = 0; y´(0) = 1; y´´(0) = -2 Formar las ecuaciones diferenciales de las familias de curvas que se dan a continuación. 1. y 2 2cx 2. x 2 y 2 c 2 3. x 3 c x 2 y 2 4. y c1e 2 x c 2 e x 5. y 2
1
x
2 ce
y2 2
6. Hallar la ecuación diferencial de la familia de curvas que cumplen con la siguiente propiedad: “Si por un punto cualquiera de una curva de la familia se trazan las rectas
tangente y normal a ella, el área del triangulo formado por dichas rectas con el eje y es igual a
x 2 y 0 2
, donde y0 es la coordenada del punto en que la tangente corta al eje y.
7. Hallar la ecuación diferencial de la familia de circunferencias de radio variable r, cuyos centros están sobre el eje X. 8. Por un punto p(x,y) de una curva que pasa por el origen, se trazan paralelas a los ejes coordenados, las que determinan un rectángulo con dichos ejes. Hallar la ecuación diferencial de la curva, de modo que ésta divida al rectángulo formado en dos regiones, donde el área de la parte derecha sea el cuádruple del área de la parte izquierda. 9. Hallar la ecuación diferencial de todas las tangentes a la parábola x 2 2 y 1 Determinar las Ecuaciones diferenciales de los problemas físicos siguientes: 1. El radio se desintegra a una velocidad proporcional a la cantidad Q. 2. Para cierta substancia el cambio de presión de vapor (P) respecto a la temperatura (T) es proporcional a la presión de vapor e inversamente proporcional al cuadrado de la temperatura.
14
Análisis Matemático
III
3. En la práctica, un cuerpo B de masa m que va cayendo (tal como un hombre que desciende en un paracaídas) encuentra una resistencia del aire proporcional a su velocidad instantánea v (t). Encuentre la velocidad del cuerpo en cualquier instante. 4. Un circuito en serie contiene un resistor y un inductor, tal como se muestra en la figura. Determinar la ecuación diferencial de la corriente i (t) si la resistencia es R, la inductancia es L y la tensión aplicada es E (t).
Figura 1.4 5. Un circuito en serie contiene un resistor y un capacitor, tal como se muestra en la figura, encuentre la ecuación diferencial para la carga q (t) del capacitor si la resistencia es R, la capacitancia es C y el voltaje aplicado es E (t).
Figura 1.5
15
Análisis Matemático
III
UNIDAD II 2. Ecuaciones de primer orden. Las ecuaciones diferenciales de primer orden son de la forma: F(x, y, y´) = 0
(1) o
M (x, y) dx + N(x, y) dy = 0
Estas ecuaciones se clasifican de la siguiente manera: 1. Ecuaciones de variables separables 2. Ecuaciones homogéneas 3. Ecuaciones exactas 4. Ecuaciones lineales.
2.1 Ecuaci ones de variables separables. Se llaman ecuaciones diferenciales de variables separables, a las ecuaciones diferenciales ordinarias de primer orden y grado uno, que mediante cualquier procedimiento algebraico se les puede dar las siguientes formas: y f ( x). g ( y ) . (1) o bien: f ( x). g ( y )dx f 1 ( x). g 1 ( y )dy 0 : (1ª)
Técnicas de solución: En (1), dividiendo por g (y) y multiplicando por dx, tendremos: dy g ( y )
f ( x)dx. Luego de integrar esta expresión obtendremos la Solución general de la
ecuación dada, en la forma.
dy
g ( y) f ( x)dex C .
De la misma manera, dividiendo los dos miembros de la ecuación (1ª) por f 1 ( x). g ( y) e integrando tendremos la Solución general de la ecuación (1ª) en la forma: f ( x)
f ( x) 1
dx
g 1 ( y )
g ( y) dy C
Nota: Si para un valor determinado de y = y0, tenemos que g (y0) = 0, la función
y = y0
también es Solución de la ecuación (1). De la misma manera, las rectas x = a y y = b serán soluciones de la Ecuación (1ª), si son raíces de f 1(x) = 0 y g (y) = 0.
16
Análisis Matemático
III
Ej emplo 1. Resuelva la ecuación diferencial
4 y yx dy 2 x xy dx 0 2
2
Solución:
En este caso podemos separar variables de la siguiente manera: y 4 x 2 dy x 2 y 2 dx 0 , Dividiendo por: x,y obtenemos:
y x dy 2 y 2 4 x 2 dx 0
Integrando:
y 1 1 x 2 2 2 y 2 dy 4 x 2 dx 0 = 2 ln2 y 2 ln4 x ln(c)
2 y 2 ln(c) , de donde la Solución viene dada por. 2 y 2 c4 x 2 // ln 2 4 x Ej emplo 2. Resolver
dy dx
1 y 2
1 x 2 xy
Solución: 1 x 2 xydy 1 y 2 dx 1 y 2 dx 1 x 2 xydy 0
A fin de separar variables, dividiremos por x1 x 2 1 y 2 quedando la ecuación anterior así: dx x1 x
2
ydy
1 y 2
0
Integrando Por fracciones parciales, tenemos: dx
x
xdx 1 x
2
ydy 1 y
2
ln(c)
ln(1 x 2 )(1 y 2 ) ln Cx 2 ; C=2c.
Ej emplo 3. Resolver: Solución:
1 x 3
dy dx
ln x
1 2
ln1 x 2
1 2
1 x 1 y Cx 2
2
ln1 y 2 ln(c) . De donde
2
.////
x 2 y x 2 .
1 x 3 dy x 2 ( y 1)dx
A fin de separar variables, dividiremos todo por
1 x 3 .( y 1) quedando la ecuación anterior así: dy y 1
x 2 dx 1 x 3
. Integrando:
x 2 dx
dy
y 1
1 x 3
C ln( y 1)
2 3
1 x 3 C /// .
Ej emplo 4. Resolver: e x y senxdx (2 y 1)e y dy 0 2
17
Análisis Matemático
III
e x .e y senxdx (2 y 1)e y dy 0 2
Solución:
A fin de separar variables, dividiremos
todo por e y quedando la ecuación anterior así: e x senxdx (2 y 1)e ( y
2
y)
dy 0
Integrando el primero por partes y el segundo por
cambio de variable tenemos: e x senxdx (2 y 1)e ( y e x ( senx cos x) 2e ( y
2
y)
2
y)
dy C
C /// .
2.1.1 Ecu aciones Di fer enci ales que pueden r educi rse a variables separables.
Las ecuaciones diferenciales de la forma
y f ax by c
b 0 se reducen a
ecuaciones en variables separables por medio de la sustitución u ax by c donde u es la nueva función que se busca. Ej emplo 1. Resolver:
dy dx
1 x y 1
Solución: Se puede ver que esta ecuación no es de variables separables, si tomamos el cambio
de variables: u x y Se tiene: du dx
1
dy dx
Así la ecuación se transforma en: du dx
1
1 u 1
Integrando:
O bien
u 1du u2
u 1du u2
dx , que es una ecuación en variables separables.
1 dx ln c 1 du dx ln c 1 u
Resulta: u ln 1 u x ln c Reemplazando: x y ln x y 2 x ln c Lo que equivale a: x y 2 c1e y //. Ej emplo 2. Resolver: xy 2 x y y a 2 Solución: Se puede ver que esta ecuación no es de variables separables, si tomamos un
cambio de variable adecuado, en este caso: u xy y
u x
y
x
du dx x 2
u
18
Análisis Matemático
III
Reemplazando en la ecuación dada tenemos: du x u 2 2 u 2 dx u a 2 , Simplificando x du u u x a u 2 du a 2 xdx. , que es una x dx x x x 2 u2
ecuación en variables separables. 1
1
3
2
Integrando: u 3 a 2 x 2 C . /// 3
Reemplazando: x 3 y 3 a 2 x 2 C 1 . /// 2
Ej emplo 3. Resolver: (1 x 2 y )dx x 2 ( y x)dy 0 Solución: Se puede ver que esta ecuación no es de variables separables, si tomamos un
cambio de variable adecuado, en este caso: z x y y x z y dy dx dz Reemplazando en la ecuación dada tenemos: 1 x 2 ( x z ) dx x 2 ( z )(dx dz ) 0 Resolviendo los paréntesis y simplificando tenemos: (1 x 3 x 2 z x 2 z )dx x 2 zdz 0 (1 x 3 )dx x 2 zdz 0 que es una ecuación en variables
separables. Separando variables e integrando:
1 x 2
dx xdx zdz 0
(1 x 3 ) x 2 1 x
x 2 2
dx zdz 0
z 2 2
(1 x 3 ) x 2
c z 2 x 2
2 x
dx zdz 0
C /// .
Reemplazando z = x – y, y simplificando tenemos: x 2 2 xy y 2 x 2 y 2 2 xy
Ejemplo
2 x
2 x
C
C /// .
4.
Resolver:
y (8 x 2 y) 2 2(8 x 2 y ) 1 :
decodificando
tendremos:
y 4(4 x y ) 2 4(4 x y ) 1
Solución. Se puede ver que esta ecuación no es de variables separables, si tomamos un cambio de variable adecuado, en este caso: z 4 x y y z 4 x y
dy dx
dz dx
4
Reemplazando en la ecuación dada tenemos:
19
Análisis Matemático dz dx
III dz
4 4 z 2 4 z 1 . Simplificando y separando variables tenemos:
z z 2
Integrando:
dz z z 2
5
4 dx C completando cuadrados en la primera
5
4dx
4
integral
4
tendremos:
dz z 2 z
1
1
4
4 dx C
dz 2 z 1 4 dx C arctan g ( ) 4 x C 1 2 2 ( z ) 1 2
Reemplazando z = 4x + y, y simplificando tenemos: arctan g (
8 x 2 y 1 ) 4 x C /// . 2
2.1.2 Ej ercici os propuestos..
Resolver las ecuaciones diferenciales: 1. tgxsen 2 ydx cos 2 xctgydy 0 2. x y y y 3 3. xy y 1 x 2 4. y x y a1 x 2 y 5. y tagx y 6.
1 x 2 dy 1 y 2 dx 0
7. y (8 x 2 y 1) 2 8. (2 x 3 y 1)dx (4 x 6 y 5)dy 0 9. (2 x y)dx (4 x 2 y 3)dy 0 10. senx 11. x 2 y
dx dy
dy dx
cosh 2 y 0 1 x csc y
12. 1 e x y y e x ; y 1 parax 0 y; y 1 parax 13. y senx y ln
2
20
Análisis Matemático 14.
1 y e y
III y 2
y
0
x ln x
15. x 2 y 3 y x 2dx x 3 y 2 x dy 0 16.
dy dx
x y
x y 2
17. y sen( x y) 18. ( x 2 y 2 1)dx 2 x 2 dy 0 19. y 20.
dr d
e y 2 y xe y
sen e 2r sen
, r ( ) 0. 2 3e r e r cos 2
2.2 Ecu acion es homogé neas. 2.2.1 F unciones homogé neas. Definición.- Una función f(x, y) se dice que es homogénea de grado n si: f (tx, ty) = t n f(x, y) Ej emplo 1. Comprobar si la siguiente función es homogénea. f(x, y) = x + y
Solución:
f (tx, ty) = t(x+ y) = t f(x, y), es homogéneo de grado uno.
Ej emplo 2. Comprobar si la siguiente función es homogénea. f (x, y) = sen (x + y)
Solución: f (tx, ty) sen(tx ty) sen.t ( x y) tf ( x, y), no cumple con la condición, por la
tanto no es homogénea. Ej emplo 3. Comprobar si la siguiente función es homogénea. f ( x, y) x. sen
Solución: f (tx, ty ) tx. sen
y x
ty tx
y
y. sen . x
ty. sen
y y t x. sen y. sen tf ( x, y). es tx x x
ty
homogénea de
grado uno.
21
Análisis Matemático
III
Si M (x, y) dx + N (x. y) dy = 0 es una ecuación diferencial ordinaria en la que los coeficientes M (x. y) y N(x, y) son funciones homogéneas del mismo grado, se llaman ecuaciones
diferenciales homogéneas. Té cnica de soluci ón:
Para que estas ecuaciones puedan transformarse
en ecuaciones de variables Separables u
podemos realizar los siguientes cambios de variable.
y x
O u
x y
Ej emplo 1. Resolver la ecuación 2 x 3 ydx ( x 4 y 4 )dy 0 Solución:
Cada coeficiente, es una función homogénea de grado 4, por lo tanto la ecuación dada es homogénea del mismo grado, por lo que se puede utilizar el cambio de variable: u
x y
o x = uy diferenciando tendremos dx udy ydu
Después de efectuar las sustituciones, la ecuación se transforma en: 2u 3 y 4 udy ydu u 4 y 4 y 4 dy 0
Realizando algunas operaciones se reduce a: Integrando
se
obtiene:
3 x 4 y 2 y 6 c , con
2u 3 du
3u
c c1
6
4
1
dy y
2u 3 du 3u 4 1
ln c
dy y 1 6
0
ln(3u 4 1) ln y ln c1 :
o
bien
que es la solución.
Ej emplo 2. Resolver la ecuación ydx 2 xy x dy 0 Solución:
Cada coeficiente, es una función homogénea de grado 1, por lo que la ecuación dada es homogénea del mismo grado, para resolverle utilizamos el cambio de variable: u
x y
o x = uy derivando tendremos dx udy ydu
Después de efectuar las sustituciones, la ecuación queda: y udy ydu 2 uy 2 uy dy 0
22
Análisis Matemático
III
Realizando algunas operaciones se reduce a: Integrando se obtiene:
du
2
u
dy y
1 du 2 u
C
dy y
0
u ln y ln C u ln
C y
///.
Volviendo a variables originales y realizando operaciones tendremos: 2
C ln x y ln . Que es la solución. y y y
x
C
Ej emplo 3. Resolver la ecuación
( x y y) arctan
y x
x; con la condición de que y = 0 para
x = 1. Solución:
Cada coeficiente, es una función homogénea de grado 1, por lo que la ecuación dada es homogénea del mismo grado, para resolverle utilizamos el cambio de variable: u
y x
o y = u·x derivando tendremos y u x
du dx
.
Después de efectuar las sustituciones, la ecuación queda: du x ( u x ) ux arctan u x dx
Realizando algunas operaciones se reduce a: arctan u.du Integrando:
arctan u.du
dx
x
dx x
ln C u arctan u 1 / 2 ln 1 u 2 ln x ln C
u arctan u ln Cx 1 u 2
Volviendo a Variables originales tenemos: y arctan
y x
x ln C x 2 y 2 /// . Que es la Solución
general. Con la condición de que y = 0 para x = 1 encontramos la Solución particular. 0 ln C C 1.
Reemplazando este valor encontramos la Solución pedida.
23
Análisis Matemático y arctan
y x
III
x ln x 2 y 2 /// /// .
Nota: A veces es conveniente intentar los dos cambios de variables sugeridos, ya que uno de ellos puede conducir a una ecuación más fácil de integrar. 2.2.2 E cuaci cu acion ones es r educi bles bl es a hom ogé n eas. eas.
a x b1 y c2 : y si el determinante Las ecuaciones del tipo: y 1 a x b y c 2 2 2
a b 1 1 0, a1 b2
haciendo la sustitución x u , y v , donde las constantes y se determinan por el sistema de ecuaciones: a1 b1 c1 0 y a2 b2 c2 0 , obtenemos una ecuación diferencial homogénea respecto a las variables u y v . Si 0 poniendo en la ecuación tipo a1 x b1 y u , obtenemos una ecuación en variables separables.
Otra forma de transformar
una ecuación diferencial no homogénea, en homogénea, es
mediante la sustitución y z , ocurriendo esto cuando todos los términos de la ecuación son del mismo grado, dándole el grado 1 a la variable x, el grado a la variable y, y el grado -1 a la derivada dy/dx. Además podemos transformarlas a ecuaciones homogéneas utilizando sustituciones adecuadas. Ej emplo 1. Resolver la ecuación
2 x y 4dy x 2 y 5dx 0
Solución:
Primero hallamos el determinante de los coeficientes: 1 2 3 Como 3 0 realizamos las sustituciones respectivas. 2 1 x u Y
y v , encontramos los valores de y por medio de las siguientes
ecuaciones: 1. 2 4 0
2. 2 5 0
Resolviendo el sistema: = -1
y =2
Entonces la sustitución sería: x u 1 y
y v 2 ,
dx = du
y dy = dv.
24
Análisis Matemático
III
Reemplazando en la ecuación original queda:
2u 1 v 2 4du u 1 2v 2 5 dv 0 . 2u v du u 2v dv 0 . Que es una ecuación homogénea de
Simplificando.
grado uno, por lo que realizamos la sustitución: u vz y du vdz zdv Reemplazando y simplificando tendríamos: dv v
z 1 z 2 1
dz 0
v
1
2 z
2 z 2 1
v 2 ( z 2 1)( z 1) 2 ( z 1) 2
dz 2
1
z 1
2
z 1
ln z 2 1 ln
Integrando: ln v ln
dv
ln c
v 2 ( z 1) 3 z 1
z
dz
2
1
0
ln C
c , Sustituyendo: z
u v
y
Volviendo a variables originales Tendremos la Solución: x y 13 C 1 x y 3. /// Ej emplo 2. Resolver:
x y dx 6 xy 3
2
dy 0
Solución: y z dy z 1 dz ,
Si hacemos
x z dx 6 xz 3
reemplazando en la ecuación dada tenemos:
. z 1dz 0 x z 3 dx 6 xz 3 1. dz 0. (1)
2
Luego: 1
es el grado de x
3 es el grado de z 3 3 es el grado de 6 xz 3 1 , Entonces para que la ecuación (1) sea homogénea debe cumplirse: 3 1
1 3
, remplazando en la ecuación (1) y realizando algunas simplificaciones,
obtendremos la ecuación homogénea de grado uno en en x y z. x z dx 2 xdz 0 // A continuación realizamos la sustitución
z ux dz udx xdu , remplazando en la
ecuación homogénea tenemos: x uxdx 2 xudx xdu xdu 0, De donde realizando algunas simplificaciones y separando variables se tiene:
2du 1 3u
dx x
0,
integrando 2
du
1 3u
dx x
2
0 ln 1 3u ln x ln C 1 3u 3
2
3
C x
.///.
25
Análisis Matemático
Como u
z x
III 1
, y y z z y , obtenemos la respuesta: x 3 y 3
3
dy
Ej emplo 3. Resolver:
dx
2 3 3
1
Cx 3 /// /// .
y 2 x 2 xy
Solución:
Si hacemos z y 2 y z , derivando
dy dx
1
.
dz
2 z dx
, reemplazando todo en la ecuación
dada, esta se transforma en ecuación una homogénea de grado uno. 1
.
dz
2 z dx
z x 2 x z
Sea z = u.x, se tiene: u x
du dx
Simplificando tenemos:
dz dx
u x ux x x
du dx
dz z x dx
x
.
, reemplazando (2) en (1), simplificando y separando variables
x
du dx
u 1 u du
dx x
0, integrando du
dx x
0 , de
/// . donde u ln x C ///
Volviendo a variables originales la respuesta quedaría así: y 2 x(C ln x) /// /// . 2.2.3 Ej er cici os propues pr opuestos. tos.
Determinar cuál de las siguientes funciones son homogéneas x
1. f ( x, y) y 2 tg
y
2. f ( x, y ) 4 x 4 y 4 x y 3. f ( x, y) x.arctag y.arctg x y 2 x
4. f ( x, y) x y e y 4 xy. yx 3 x 2 y 2 . 5. f ( x, y) 3 5 x 2 y 4 2
2
Resolver Las siguientes ecuaciones diferenciales: 1. x 2 y 2 dx xydy 0
26
Análisis Matemático
III
2. x 2 dy y 2 xy dx 0 y
3. xe y dx xdy 0 x
4. y
y x
cos
y x x
y y y 5. xsen y cos dx x cos dy 0
6.
x
x
x y
x
x y dx
x y
x y dy 0
2 2 y y 7. x y y.arcsen x dx x.arcsen x dy 0
8. xdy ydx x 2 y 2 dx 9. Hallar la solución particular de la ecuación x 2 3 y 2 dx 2xydy 0 , con la condición de que y = 1 para x = 2 10. y
1 3 x 3 y 1 x y
.
x 2 y 1
11. y
2 x 4 y 3
12. x y 1dx y x 3dy 0 13. x 2 y x
dy dx
xy 2 y 0
14. x 2 2 y 3 dx 3 xy 2 dy 0 12 12 15. 4 xy 6 y dx 4 y 3 x dy 0, sugerencia
16.
dy dx
17. 2
dy dx
18. y
3 x 2 y x 5 y x
3
1
u y 2 .
, sugerencia z = x3.
y 4 x x 2 y x
y 2 2 xy x 2 y 2 2 xy x 2
.
; y = -1, para x = 1.
27
Análisis Matemático 19. y(
dy dx
III
) 2 2 x
dy dx
y 0; y 5 , para x = 0.
20. (3 x 5 3 x 2 y 2 )dx (2 y 3 2 x 3 y)dy 0
2.3 Ecuaci ones exactas o en di fer enci ales total es. Una ecuación M(x, y) dx + N (x, y) dy = 0 se dice que es una ecuación exacta o un diferencial exacto, si cumple con la condición necesaria y suficiente de que:
M N . y x
Técnica de Solución: La técnica para resolver esta clase de ecuaciones es por integración respecto a x e y respectivamente, determinando así una Solución: f ( x, y) c. Ej emplo 1. Resolver: e 2 y y cos xy dx 2 xe 2 y x cos xy 2 y dy 0 Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: Así tenemos que: M ( x, y) e 2 y y cos xy
(1)
N ( x, y) 2 xe 2 y x cos xy 2 y
De (1)
M 2e 2 y xysenxy cos xy y
De (2)
N 2e 2 y xysenxy cos xy x
entonces
y
(2)
M N y x
lo que nos indica que es
una ecuación exacta, entonces si existe una función “f” tal que:
f dx
M f x ( x, y) ,
f dy
N f y ( x, y ).
Integrando respecto a x e y (1) y (2) respectivamente tendríamos: f ( x, y ) xe 2 y senxy . f ( x, y ) xe 2 y senxy y 2
Comparando, se tiene: f ( x, y ) xe 2 y senxy y 2 .
Así la Solución es: xe 2 y senxy y 2 c
28
Análisis Matemático Ej emplo 2. Resolver:
III
2 xe
y
e x dx x 2 1e y dy 0
Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: Así tenemos que: M ( x, y) 2 xe y e x
(1)
N ( x, y ) x 2 e y e y
(2)
De (1).
M 2 xe y y
De (2).
N 2 xe y x
Entonces
y
M N . y x
Lo que nos indica que es una ecuación
exacta, entonces si existe una función “f” tal que:
f dx
M f x ( x, y) ,
f dy
N f y ( x, y).
Integrando respecto a x e y, (1) y (2) respectivamente: f ( x, y ) x 2 e y e x f ( x, y ) x 2 e y e y
Comparando, se tiene: f ( x, y ) x 2 e y e x e y
Así la Solución es: x 2 e y e x e y c Ej emplo 3. Resolver: (3 x tan y 2
2 y 3 x 3
)dx ( x sec y 4 y 3
2
3
3 y 2 x 2
)dy 0
Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: Así tenemos que: M ( x, y) 3 x tan y 2
2 y 3 x 3
N ( x, y) x sec y 4 y 3
De (1).
2
3
(1) 3 y 2 x 2
y
(2)
M 6 y 2 2 2 3 x sec y 3 y x
29
Análisis Matemático
III
6 y 2 N 2 2 3 x sec y 3 x x
De (2).
Entonces
M N . y x
Lo que nos indica que es una
ecuación exacta, entonces si existe una función “f” tal que:
f dx
f
M f x ( x, y) ,
dy
N f y ( x, y).
Integrando respecto a x e y, (1) y (2) respectivamente: f ( x, y) x tan y 3
f ( x, y) x tan y 3
y 3 x 2 y 3 2
x
y4
Comparando, se tiene: f ( x, y) x tan y 3
y 3 2
x
y4
Así la Solución es: x tan y 3
y 3 x
2
y 4 c
Nota: Puede probarse que cada ecuación diferencial de primer orden: M(x, y) dx + N (x, y) dy = 0, puede transformarse en exacta, multiplicando por un factor de
integración ( x, y) adecuado, tal factor tiene la propiedad de que u x, y M x, y u x, y N x, y. y x
Desafortunadamente, no es fácil determinar dicho factor de integración u, si no en algunos casos especiales, porque en otros podría resultar una experiencia frustrante. 2.3.1 Casos especiales en que se determi nan factor es integran tes.
Dada la ecuación M(x, y) dx + N (x, y) dy = 0 que no es exacta, entonces:
1.-
2.-
si:
si:
M N y x N
M N y x M
f ( x), ( x, y) e
f ( x ) dx
g ( y), ( x, y ) e
g ( y ) dy
30
Análisis Matemático
3.-
si:
III
f ( x ) dx g ( y ) dy M N Nf ( x) Mg ( y), ( x. y) e y x
Ej emplo 1. Resolver:
( x y)dx x ln xdy 0
Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: M 1 y N 1 ln x x
Entonces
M N . Lo que nos indica que la ecuación no es exacta, y x
entonces buscaremos el factor de integración adecuado.
1.-
1 1 ln x x ln x
1 x
que es solo función de x
De donde e
F ( x ) dx dx
e
x
1
e ln x e ln x x 1
De donde el factor integrante es
1 x
1
x
Multiplicando la ecuación por el factor: y 1 dx ln xdy 0 x
1 x
x y dx x ln xdy 0 Es exacta.
Entonces si existe una función “f” tal que: dx
M ,
dy
N
Integrando respecto a x e y respectivamente, y comparando, se tiene: ( x, y ) x y ln x
31
Análisis Matemático
III
Así la Solución es: x y ln x c Ej emplo 2. Resolver: ydx – (2x+y) dy= 0 Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: M 1 y N 2 x
Entonces
M N . y x
Lo que nos indica que la ecuación no es exacta,
entonces buscaremos el factor de integración adecuado. 1 2
1.-
2 x y
1
que no es solo función de x, por lo que probamos con el segundo
2 x y
caso. 2 1
2.-
y
ue
3
como vemos es solo función de y, por lo que su factor integrante es
y
3
y dy
e ln y 3
1 y 3
De donde la ecuación 1 y3
ydx
1 y3
2 x ydy 0
es exacta.
Entonces si existe una función tal que: dx
M ,
dy
N
Integrando respecto a x e y respectivamente, y comparando, se tiene: ( x, y )
x y 2
1 y
Así la Solución es:
x y 2
1 y
c
32
Análisis Matemático
III
Ej emplo 3. Resolver: y2 dx+ (x2+xy) dy =0 Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: M 2 y y N 2 x y x
Entonces
M N . y x
Lo que nos indica que la ecuación no es exacta,
entonces buscaremos el factor de integración adecuado. M N 2 y 2 x y f ( x) N g ( y) M y x
3.-
3 y 2 x f ( x)( x 2 + x.y)- g(y) (-y2) = f(x) (x2+ x y) + (g (y) (y2) = xy f(x) +y2 g (y)
Haciendo: - 3y = xy f (x) + y2 g (y),
-2x = x2 f (x)
Tenemos: f ( x)
2
2 3 y xy y 2 g y x
y
x
1
de donde g ( y) , por lo que el factor y
de integración será: e
f x dz g y dy
2
1
f x dx g y dy x dx y dy 1 e e 2
x y
Es un factor de integración, de donde la ecuación 1 2
y dx
1
2
x y
2
x y
x
2
xy dy 0
Es exacta
Entonces si existe una función tal que: dx
M ,
dy
N
Integrando respecto a x e y respectivamente, y comparando, se tiene: 33
Análisis Matemático ( x, y)
y x
III
ln y
Así la Solución es:
y
ln y c
x
Nota: Al calcular el factor integrante, no es necesario poner la constante de integración. Ej emplo 4. Resolver: x 2 y dx xdy 0 Solución:
Comprobamos si la ecuación dada es exacta aplicando la condición: M 2 y N 1 x
Entonces
M N . y x
Lo que nos indica que la ecuación no es exacta,
entonces buscaremos el factor de integración adecuado. Podrá ser?
1
x 2
Compruébalo.
2.3.2 EJERCI CI OS PROPUESTOS.
Resolver las siguientes ecuaciones: xy dy 0 1. x y 1 dx y 2 y 1
2
2. xdy ydx x 3 dx 0 3. e x ln y
y x
x
y
dx ln x seny dy 0
y 2 dx 2 yarctagx 2 x senhydy 0 {ñ 2 y 4. 2 1 x
5. x 3 3 xy 2 2dx 3 x 2 y y 2 dy 0 6.
2 xdx y 3
y 2 3 x 2 y 4
dy 0
7. Hallar la integral particular de la ecuación satisfaga a la condición inicial
x x x y x e dx e y 1 dy 0 que y
y (0) = 2 34
Análisis Matemático 8.
y x
III
dx y 3 ln x dy 0
9. y1 xy dx xdy 0 10. x cos y yseny dy xseny y cos y dx 0 11. ( xy 2 y x)dx x( xy 1)dy 0 1 2 x 2 2 y y 12. ye 2 x 3 xe 2 y dx e 3 x e e dy 0; y(1) = 0. 2
13. Determinar el valor de a y b para que la ecuación diferencial sea exacta y resolverla. ( y x 3 )dx (ax by 3 )dy 0
14. ( xy 3 1)dx x 2 y 2 dy 0 x 2 x. senx )dy 0 15. ( senx x. cos x)dx 2( 2 y y
16. Demostrar que
1 x. M y. N
, donde xM + yN ≠ 0, es un factor integrante para la
ecuación diferencial Homogénea M(x,y) +N(x,y) = 0 17.
1 x
dx (1 xy 2 )dy 0.
18. senx(2 3 ysen 2 x)dx sec xdy 0 19. ( x 2 y)dx xdy 0 20.
y senx cos 2 ( xy) cos 2 ( xy)
dx
x cos 2 ( xy)
dy seny.dy 0
2.4 Ecuaciones lineales. Una ecuación de primer orden de la forma: M x, y dx N x, y dy 0 se dice que es lineal si mediante cualquier procedimiento algebraico podemos escribirles en las formas: dy dx
p x y q x ( I ) o
dx dy
p y x q y .
Si q(x) = 0, o q(y) = 0, la ecuación lineal toma el nombre de ecuación lineal homogénea, p ( x ) dx. . que es una ecuación en variables separables y su Solución es: y Ce
35
Análisis Matemático
III
Técnicas de Solución: a.) Utilizando la sustitución De donde
dy dx
u
dv dx
v
y = u.v (a), donde u y v son funciones de x;
du dx
;
reemplazando en (I)
du dv du dv u v p( x)u.v q( x) p( x)u v u q( x) , (D) dx dx dx dx
Si se exige que du dx
P ( x ) dx p( x)u 0 , resolviendo por variable separables tenemos: u e
v q( x)e
Reemplazando “u” en (D) hallamos v:
Reemplazando
u
y
v
P ( x ) dx
dx + c
en (a) tendremos la Solución general
de
( I ) dada por:
p ( x ) dx P ( x ) dx p ( x ) dx p ( x ) dx y e ( q( x).e c) ye q( x)e dx c
y.u q( x).udx c siendo u e
p ( x ) dx
o
. // Ejemplo 1.
b.) Por medio del método de la Variación de la constante arbitraria. Dada la ecuación Lineal en x:
dy dx
p x y q x ( I )
1. Encontramos la Solución de su respectiva ecuación lineal Homogénea es decir: p ( x ) dx. y Ce . (2).
2. Tratando en la expresión (2) a C como una función de “x” , derivamos esta y encontramos y
d p ( x ) dx Ce .(3) dx
3. Reemplazando (2) y (3) en (I), encontramos C(x). Este valor reemplazamos en (2) p ( x ) dx. . Ejemplo 2 y encontramos la Solución general de la forma: y C ( x)e
Ej emplo 1.
4
Resolver: y y x
o
x
4 y y x . (I) x
Solución:
Realizamos la sustitución Derivando tenemos:
dy dx
u
y = u.v (a), dv dx
v
du dx
; (b), reemplazando (a) y (b) en la ecuación (I)
36
Análisis Matemático
III
du 4 dv du 4 dv u v u.v x u v u x , (D) dx x dx dx x dx
Si se exige que du dx
4
u 0 , resolviendo por variable separables tenemos: u e
4
dx
x
u x 4 //
x
v x. x 4 dx c v
Reemplazando “u” en (D) hallamos v:
Reemplazando y x 4 (
u
1
y
v
1
c
2 x 2
en (a) tendremos la Solución general
de
(I) dada por:
1
c) y x 2 cx 4 //
2 x 2
2
Comprobación: 4 p ( x) , x
Tenemos que integración es
p( x)dx 4
dx x
4
4 ln x ln x , de donde el factor de
u x 4 , por lo que la Solución general de la ecuación dada es:
yx 4 x. x 4 dx c y x 4 x 3 dx c
y x 4 .
1 2 x 2
1
cx 4 y x 2 c.x 4 // 2
Resolver: xdy 2 ydx x 2e x dx.
Ej emplo 2. Solución:
Ponemos la ecuación de la forma:
dy dx
2
x 2
x
x
. y
e x . (I)
1. Encontramos la Solución de la respectiva ecuación lineal homogénea.
dy dx
2 x
y0
2
y C ( x)e
x dx
2
y C ( x)e ln x y C ( x). x 2 . (2)
2. Derivando (2), y tratando a C como función de x tenemos: y
dy dx
C ( x). x 2 2 xC ( x) (3).
3. Remplazando (2) Simplificando: C ( x)
y (3)
x 2.e x x
3
en (I).
C ( x). x 2 2 xC ( x)
2 x
C ( x). x x 2 .e 2
x
x
.
. Integrando: C ( x) x 2 e x dx 2 x 3 e x dx K
37
Análisis Matemático
III
C ( x) x 2 e x 2 x 3 e x dx 2 x 3 e x dx K C ( x) x 2 e x K . // , este valor reemplazamos
en (2) y encontramos la Solución general. y x 2 e x K . x 2 y e x Kx 2 . /// 2
dx
x
x
Comprobando tenemos: p( x) p( x)dx 2 de integración es y. x
2
x 2e x x
2
2 ln x ln x , de donde el factor
u x 2 , por lo que la Solución general de la ecuación dada es: 2
. x dx c y x
2
x 2 e x x
3
dx c
y x 2 x 2 e x dx 2 x 3 e x dx c
y x 2 x 2 e x 2 x 3 e x dx 2 x 3 e x dx c y x 2 . x 2 e x cx 2 y e x cx 2 // Ej emplo 3. Resolver:
dr 2rctg sen2 d 0
Solución:
Ponemos la ecuación de la forma:
dr d
2ctg .r sen2 , Ecuación lineal en r.
Tenemos que p( ) 2ctg p( )d 2 ctg .d 2 ln sen ln sen , de donde el 2
factor de integración es u sen 2 , por lo que la Solución de la ecuación dada es: sen 4 c rsen 2 sen2 . sen 2 .d c rsen 2 2 sen3 . cos .d c 2. 4
De donde rsen 2 sen 4 c1 // . Siendo c1 =2C. 2.4.1 Ecuaciones de Bern oulli .
Una ecuación diferencial se dice que es de Bernoulli de grado n si puede escribirse en la forma: dy dx
p x y q x y n
Los casos n 0 y n 1 no los analizamos, ya que están contenidos en el caso de la ecuación lineal de primer orden. Esta ecuación se reduce a lineal valiéndose de la sustitución z y1n . o también se puede utilizar directamente la sustitución y = u.v.
38
Análisis Matemático
III dy
Ej emplo 1 . Resolver:
dx
y x
xy 2
Solución: 2
Dividimos por y la ecuación original.
y
2
dy dx
dz z dx
x
x
dz z dx
x
1
x
x . (A)
y y 2
Realizamos la sustitución y 12 z y 1 z
y
dy dx
dz dx
en la ecuación (A)
x (D), es una ecuación lineal en z, por lo que procedemos a
resolverla como tal. 1
dx
x
x
Entonces p( x) p( x)dx
1
ln x ln x , de donde el factor de integración
es u x 1 , entonces la Solución general de la ecuación (D) es:
zx 1 x. x 1 dx c dx c x c z x 2 cx // ;
Volviendo a variables originales
tendríamos: y
1 x 2 cx
Ej emplo 2.
//
Resolver: 2
dx dy
x y
x 3 cos y 0
Solución: 3
Dividimos por x a la ecuación dada: 2 x Realizando la sustitución x 2 z y
3
dx dy
2 x 3
dx dy
x
2
y
cos y (B). dz dy
, en
(B)
Tendremos una ecuación lineal en z: dz z dy
y
cos y (E), aplicando el método respectivo tendremos:
Como p( y)
i y
p( y)dy
dy y
ln y . , de donde el factor de integración será: u y
Entonces la Solución general de la ecuación (E) es: 39
Análisis Matemático
III
zy y. cos ydy c yseny senydy c yseny cos y c zy yseny cos y c //.
Volviendo a variables originales: x 2 y yseny cos y c // dy
Ej emplo 3 . Resolver:
dx
y
1
x 1
( x 1) 3 y 2 2
Solución: 2
Dividimos por y la ecuación original.
y
2
dy dx
dz dx
z
1
dz
2
dx
( x 1) 3
x 1
z x 1
x 1
1
( x 1) 3 . . (A)
y y 2
Realizamos la sustitución y 12 z y 1 z
y 1
2
dy dx
dz dx
en la ecuación (A)
1
( x 1) 3 . (D), es una ecuación lineal en z, por lo 2
que procedemos a resolverla como tal. Entonces p( x)
1 x 1
p( x)dx
dx x 1
1
ln x 1 ln x 1 , de donde el factor de
integración es u ( x 1) 1 , entonces la Solución general de la ecuación (D) es: z ( x 1) 1
1 2
( x 1) 3 .( x 1) 1 dx c
1 2
( x 1) 2 .dx c
1
1 ( x 1) 3 c z ( x 1)( ( x 1) 3 c). 6 6
Volviendo a variables originales tendríamos: y
6 ( x 1)(( x 1) 3 6c)
/// .
2.4.2 Ej ercici os propuestos. Hallar las soluciones de las siguientes ecuaciones: 1.2.3.-
dy dx di dt dy dx
y 2 2 x 6i 10 sen2t
2 y x
x 3
4.- y 2 dx 2 xy 3dy 0
5.- 1 y 2 dx 1 y 2 seny xy dy.
40
Análisis Matemático
III
6.- x y y e x 0.; y b para x = a. 7.- y y.tgx 8.- 2 xy
dy dx
1 cos x
.; y 0 para x = 0
y 2 x 0
1
9.- ydx x x 3 y dy 0 2
10.- 3 xdy y1 xsenx 3 y 3 senx dx 11.- 2 xy 3 y dx 2 xdy 0 12.- xdy ydx x 3 y 6 dx 13.- 8 x y y
1 y 3 x 1
14.- y cos x y 3 x cos x senx dx 2senxdy 0 15.-
dy dx
y 3 1 2 xy 2
: y (0) = 1.
Ecuaciones diferenciales de 1er orden, no resueltas respecto a la
2.5
derivada. er
2.5.1 Ecu aciones dif erencial es de 1 orden, de grado superior.
Su forma general es: f ( x, y, y , y 2 ,..., y n ) 0 (A)
Técnica de Solución:
Resolvemos Algebraicamente con respecto a y´, (Podemos reemplazar y´ por z) obteniendo n soluciones algebraicas de la forma:
y f 1 ( x, y ),
y f 2 ( x, y), ......
( B ),
De cada Solución Algebraica que es una ecuación diferencial de primer orden y primer grado, encontramos las soluciones.
La Solución general de la ecuación diferencial dada será la multiplicación de estas soluciones. F ( x, y, C ) F 1 ( x, y, C ). F 2 ( x, y, C )..... 0, ( C ).
41
Análisis Matemático
III
Además, para la ecuación (A) puede existir una Solución Singular, que geométricamente representa la envolvente de la familia de curvas (C), y se puede obtenerse eliminando C del sistema de ecuaciones: F C ( x, y, C ) 0
F ( x, y, C ) 0,
O eliminando p = y´ del sistema de ecuaciones: f p ( x, y, p) 0.
f ( x, y, p) 0,
Nota: Debe comprobarse dichas soluciones singulares ya que no todas las veces son Solución de la ecuación (A). Ej emplo 1. Resolver:
x y 2 2 x y y 0
Solución:
Dividiendo todo por x y luego resolviendo con respecto a y´ (como si se tratase de una ecuación algebraica de 2do grado), tenemos: 2 4 y1 2
2
y1 1 1
y x
y x
2 2 1 y1 2
, (a)
y x
y1 2 1 1
2
y 2 1 1
y
y x
y x
, de done tenemos que:
, (b) que son ecuaciones homogéneas, por lo
dy
que tendremos que utilizar la sustitución y = u x y
dx
u x
du dx
,
Sustituyendo en la ecuación (a) tendremos: u x
du dx
1 1 u x
du dx
1 u 1 1 u
du
1 u 1 1 u
dx x
Integrando: 2 ln 1 u ln x ln C ln 1 u
2
ln
C x
1 u 2
C x
, Volviendo a variables originales
2
y C 1 1 ///. x x
Resolviendo la ecuación (b) de la misma manera, su Solución será: 42
Análisis Matemático
III
2
y C 1 1 ///. x x
Multiplicando entre sí, obtenemos la Solución general de la ecuación dada.
2 x y C 2 4 x 2 xy 0 , simplificando
y C 2 4Cx /// . (c).
(Familia de parábolas). Para encontrar la Solución singular derivamos respecto a C la Solución general. 2 y C 4 x C 2 x y.(d ). Poniendo (d) en (c).
y 2 x y 2 42 x y x x y 0 // . (e) Solución singular. Comprobación: Derivando (e) y reemplazando en la ecuación dada Tendremos: 1 y 0 y 1,
x 1 2 x 1 x 0 0 0 // . Lo que demuestra que si es una 2
Solución de la ecuación.
Ej emplo 2. Resolver: y 2
2 y x
y 0
Solución:
Resolviendo con respecto a y´ (como si se tratase de una ecuación algebraica de 2do grado), tenemos: 2 y1 2
y1
y x
y x
4
y 2 x 2
4
2
y 2 x 2
1, (a)
y1 2
y
y x
y1
y x
y 2 x 2
1 , de done tenemos que:
y 2 x 2
1, (b) que son ecuaciones homogéneas, por lo
que tendremos que utilizar la sustitución y = u x y
dy dx
u x
du dx
,
Sustituyendo en la ecuación (a) tendremos:
43
Análisis Matemático du
u x
dx
III
u u 2 1 x
du dx
du
u2 1
u2 1
dx x
Integrando:
ln u u 2 1 ln x ln C ln u u 2 1 ln Cx u u 2 1 y
Volviendo a variables originales
x
y 2 x
2
Cx,
1 Cx y 2 x 2 Cx 2 y // .
Resolviendo la ecuación (b) de la misma manera, su Solución será: y 2 x 2 C y // .
Multiplicando entre sí, obtenemos la Solución general de la ecuación dada.
x C 2
2
1 2Cy x 2 C 2 2Cy 0 , (c).
Para encontrar la Solución singular, reemplazamos en la ecuación dada p = y´, luego esta nueva ecuación derivamos respecto a
p, resolviendo este sistema hallamos la Solución
buscada. p 2 2 y 2 x 2
y x
p 1 0 : (¡)
2 p 2
y x
0 p
y x
.(2)
y y
2 . 1 0 y 2 2 y 2 x 2 0 x 2 y 2 0 // . Solución singular. x x
Comprobación: Derivando la Solución singular y reemplazando en la ecuación dada Tendremos: 2 x 2 y y 0 y
x y
2
,
x y x 2 . 1 0 x 2 y 2 0 Lo que demuestra que es x y y
una Solución de la ecuación, por lo tanto ésta es la Solución singular. 2.5.2 Ecu aciones con ausenci as de la variabl e in dependi ente (x) o dependi ente (y). 2.5.2.1 Ecuacion es de la forma, f ( x, y ) 0 ausente “y” .
Técnica de Solución: Introducimos la sustitución y´ = p, y realizamos el siguiente procedimiento:
44
Análisis Matemático Ponemos: y´ = p o
III dy dx
p , (¡); de donde tenemos que x = H(p) y dx = H´(p)dp.
De (1) dy pdx dy p H ( p)dp . Integrando tendremos: y p H ( p )dp. Por lo que la Solución general de la ecuación dada es: y p H ( p)dp que es de forma paramétrica. Ejemplo 1. x H ( p). 2.5.2.2 Ecuacion es de la forma, f ( y, y) 0, ausente (x).
Técnica de Solución: Introducimos la sustitución y´ = p, y realizamos el siguiente procedimiento: Ponemos: y´ = p o
dy dx
De (1) dy pdx dx
p , (¡); de donde tenemos que y = g(p) y dy = g´(p)dp. dy p
Integrando tendremos: x
dx g ( p) dp p
g ( p)dp p
.
Por lo que la Solución general de la ecuación dada es:
y g ( p ) x g ( p )dp . que es de forma paramétrica. Ejemplo 2. p
2.5.2.3 Ecuacion es de la forma:
Técnica de Solución:
.
Introducimos la sustitución y´ = p. Para determinar las Variable “x” y
1 p
h h dp . , p p dy
x h( y, p)
2.5.2.4 Ecuaciones de la form a:
Técnica de Solución:
“y” se determina por el sistema de ecuaciones:
Ejemplo 3.
Introducimos la sustitución y´ = p Para determinar las Variable “x” y
p
h h dp . x p dx,
“y” se determina por el sistema de ecuaciones:
y h( x, p) .Ejemplo 4.
45
Análisis Matemático . Resolver: Ej emplo 1
III x cos y ln y .
Solución: 1
Reemplazamos y´ = p (1), de donde tendríamos x = cos p + ln p, y dx ( senp )dp . p
1 y p senp dp p. senp.dp dp y p cos p senp p C . p y p cos p senp p C Entonces su Solución general es x cos p ln p. Ej emplo 2. Resolver:
y y 2 2 ln y .
Solución:
Reemplazamos y´ = p (1), de donde tendríamos y p 2 2 ln p , x
2
y dy (2 p )dp . p
1 2 dp 2 2 p dp 2 dp 2 2 x 2 p . p p p p
y p 2 2 ln p Entonces su Solución general es 2 x 2 p . p Ej emplo 3 . Resolver: 3 x
y y
6 y y 2 .
Solución:
Remplazando y p , tendremos 3 x
3
dx dy
3 p
y p
6 py 2 (1), derivando respecto a “y”, inclusive p p y
1 p
2 p y
dp dy
p 2
y dp 2
p dy
dp dy
dp dx 1 6 2 yp y 2 , reemplazando y simplificando tenemos: dy dy p p 12 yp 6 y
2
dp dy
2 p 2 p
y
dp dy
6 y 2 p 2 12 yp 3
6 y 2 p 2 12 yp 3 0 1 6 p 2 y (2 p y
dp dy
)0
46
Análisis Matemático
III
Igualando a cero el segundo factor, 2 p y 2
dy y
dp p
0 2
dy y
dp p
dp
0 separando variables e integrando tenemos:
dy
ln(C ) ln( y) 2 ln( p) ln(C ) ln( y 2 p) ln(C ) y 2 p C . //
Despejando p y sustituyendo en la ecuación (1), se halla la Solución general y C 2 y 3 6 2 . y 3 x 6C y 3 3Cx 6C 2 . /// 3 x C C y 2 y
.Ejemplo 4. Resolver: 4 y x 2 y 2 Solución:
Reemplazando y´ = p, tendremos: 4 y x 2 p 2 , Derivando respecto a “x”, inclusive p. 4
dy dx
2 x 2 p
dp dx
2 p x p
dp dx
(2 p x)dx pdp 0 (a): que es una ecuación
homogénea de 1er orden, por lo que realizaremos la sustitución
x = p.z: de donde
dx = pdz + zdp. Sustituyendo en (a),
(2 p p. z )( pdz zdp) pdp 0 (2 z ) pdz ( z 2 2 z 1)dp 0 z 2 ( z 1)
2
dp
dz
p
0
z 1 ( z 1)
2
dz
dz ( z 1)
2
dp p
0
dz z 1
dz ( z 1)
2
dp p
0,
Integrando, ln z 1
1 z 1
ln x p
ln p C ln p( z 1)
p x p
1 z 1
C , reemplazando z
x p
y simplificando
C .///
y 1 ( x 2 p 2 ) 4 Entonces la Solución sería: p C . ln x p x p
47
Análisis Matemático
III
2.6 Ecuaciones de Lagrange y de Clairaut. 2.6.1 Ecuación de Lagr ange.
La ecuación de la forma y xf ( p) g ( p). (a) donde p = y´, recibe el nombre de ecuación de Lagrange.
Técnica de Solución: Por medio de la derivación, y teniendo en cuenta que dy = p.dx, la ecuación de Lagrange se transformas en lineal con respecto a x. pdx f ( p)dx x f ( p) g ( p)dp . (b)
Si p f ( p), De las ecuaciones (a) y (b) se obtiene la Solución general en forma paramétrica x Cf 1 ( p) g 1 ( p),
y Cf 1 ( p) g 1 ( p) f ( p) g ( p).
Donde p es un parámetro y f 1 (p) y g1 (p) son funciones conocidas determinadas. Además, puede existir Solución singular, que se calcula por el procedimiento general. Ej emplo 1. Resolver: y 1 y x y 2 (A) Solución:
Sustituyendo
y´ = p, en (A) tendremos: y 1 p x p 2 (I) derivando, y teniendo en
cuenta que dy = p.dx: pdx 1 p dx xdp 2 pdp ( p 1 p)dx ( x 2 p)dp
dx dp
x 2 p. (B) que es una
ecuación lineal en x, Resolviendo (B) tenemos: u e
dp
u e p . Por lo que su Solución será:
xe p 2 pe p dp c xe p 2e p p 1 c x ce p 2 p 2. //
Sustituyendo este valor en (I) tenemos: y (1 p)ce p 2 p 2 p 2 simplificando y c(1 p)e p p 2 2. // Por lo que la Solución general de la ecuación de Lagrange dada
es. x ce p 2 p 2 y c1 p e p p 2 2 Ej emplo 2. Resolver: y 2 x y ln y (A)
48
Análisis Matemático
III
Solución:
Sustituyendo y´ = p, en (A) tendremos: y 2 xp ln p (I) derivando, y teniendo en cuenta que dy = p.dx: dy 2 pdx 2 xdp
1 p
dp pdx 2 pdx 2 x
1
dx
1
dx
2
1
dp p 2 x 0 x 2 p dp p dp p p
(B) que es una ecuación lineal en x, Resolviendo (B) tenemos: ue
2
dp
p
xp 2
2
u e ln p u p 2 . Por lo que su Solución será: 1
p
2
. p 2 dp c xp 2 p c x
c p
2
1 p
. //
c 1 Sustituyendo este valor en (I) tenemos: y 2 2 p ln p simplificando p p y ln p
2c p
2 , por lo que la Solución general de la ecuación de Lagrange dada es.
c 1 x p 2 p y ln p 2c 2. p
2.6.2 Ecuación de Clai r aut .
Una ecuación que se puede poner en la forma y xp ( p) donde p y , después de realizar algunas operaciones algebraicas es una ecuación de Clairaut, su método de resolución es el mismo que el de las ecuaciones de Lagrange y su Solución general viene dada por: y cx (c) , que no es más que la sustitución en la ecuación dada de c por p.
Esta ecuación puede tener también una Solución: singular, que se obtiene de eliminar p entre las ecuaciones: 1. y xp ( p) y
2. x P 0
Ej emplo 1. Resolver y x y
a y 2
. Donde a es una constante.
Solución:
Ponemos en la ecuación dada y p de donde tenemos: y xp
a p 2
(I).
49
Análisis Matemático
III
Derivando y reemplazando dy pdx tendremos: dy xdp pdx
2a
dp pdx pdx x
p
2a
2a dp dp x 0. p p
Examinamos los dos factores:
dp 0 dp o p c. Por lo que la Solución general de la ecuación dada es:
y cx
a c2
. // (Familia de rectas).
Igualando el segundo factor a cero tendremos: x
2a p
y x *
0 x 2a x
2a p
a
2a
2
(II). Eliminando p entre esta y la ecuación (I)
y 2a
x
Ejemplo 2.
ax 2 4a 2
y
8a 2 x 2 4a
Resolver: y 3 x y 6 y 2 y 2
Solución:
Ponemos en la ecuación dada y p de donde tenemos: y 3 xp 6 y 2 p 2 (I). Multiplicando la ecuación por y2 se obtiene y 3 3 xy 2 p 6 y 4 p 2 . Mediante la transformación y3 = z y 3 y 2 y z 3 y 2 p z , se obtiene: 2 z x z z 2 , que es una ecuación de Clairaut, de donde su Solución es: 3 z Cx
2 3
C 2 y 3 Cx
2 3
C 2 . ///
2.6.3 Ej ercici os propuestos.
1. x 2 y 2 xy y 6 y 2 0 2. x y 2 2 y y 4 x 0 3. y 2 2 y y y 2 e x 1 4. y 2 2 x y 8x 2 0 5. x y 2 2 y 2 6. x1 y 2 1
50
Análisis Matemático
III
7. y y ln y 8. y arcsen. y ln1 y 2 9. y y 1 y cos . y 10. x y sen. y 11. 2 y x y y ln y 12. y 2 x y sen y 3
13. y x y e y 2
14. y x y y 2 15. x y 2 y y y 1 0
51
Análisis Matemático
III
UNIDAD III 3. Aplicaciones de las ecuaciones diferenciales de 1er orden. Los diversos problemas que pueden resolverse con la ayuda de ecuaciones diferenciales, vienen dados por las interpretaciones de la función derivada, así por ejemplo: rapidez de variación, velocidad, tasa de cambio, pendiente, marginal, etc.
3.1 Tr ayectori as ortogonales. Dos familias de curvas, se dice que son ortogonales, si cada miembro de la una familia es ortogonal a cada miembro de la otra familia. Dos curvas son ortogonales entre sí, cuando se intersecan ortogonalmente, es decir sus rectas tangentes son ortogonales entre sí, en el punto común. Dos rectas son ortogonales entre sí, cuando el producto de sus pendientes es – 1. Así entonces, si dos familias de curvas están dadas por las ecuaciones: G1 x, y, c1 0 , G2 x, y, c2 0
Y si la pendiente de la primera familia está dada por: y f x, y entonces la pendiente de la segunda familia (la ortogonal) debe ser: y
1 f x, y
Nota. En la práctica, se da una familia de curvas y lo que se busca es la familia ortogonal. Ej emplo 1. Encontrar la familia ortogonal a la familia de hipérbolas y
c x
Solución:
En primer lugar derivando encontramos la ecuación diferencial de la familia dada (pendiente): Derivando y
c x 2
Como c xy Reemplazando tenemos: y
y x
Por lo tanto la ecuación diferencial (o pendiente) de la familia ortogonal es: y
x y
Resolviendo por variables Separables tenemos: ydy xdx ydy xdx C y 2 x 2 C///.
52
Análisis Matemático
III
Esta es la ecuación de la familia ortogonal a la familia de hipérbolas dada. Ej emplo 2. Encontrar el miembro de la familia ortogonal a la familia de curas dada por la
ecuación: 3 xy 2 2 3cx 0 que pasa por el punto 0,10. Solución:
Determinemos la pendiente a la familia dada: Derivando: 3 y 2 6 xy y 3c 0 Como: c
3 xy 2 2
y
se tiene que
3 x
1 3 x 2 y
es la ecuación diferencial de la familia
dada, por tanto la ecuación diferencial de la familia ortogonal es: dy dx
x y , integrando llegamos: y ce 2
3 2 x 2
utilizando la condición dada se tiene que c 10 ,
por lo tanto la curva buscada está dada por la ecuación: y 10e
3 2 x 2
Nota: Si la familia de curvas planas se da por un ecuación en coordenadas polares r , , a 0 , (A),
d
(A) y
d
donde a es un parámetro. Eliminando el parámetro “a” entre la ecuación
0,
encontramos
la
ecuación diferencial de la familia
(A) F r , , r 0 . Sustituyendo en esta r´ por
r 2 r
obtenemos
la ecuación diferencial de la familia de las
r 2 trayectorias ortogonales. F r , , 0. r Ejemplo 1.
Encontrar la familia de curvas ortogonales de la familia de cardioides:
r a 1 cos (A)
Solución:
Derivando (A) tenemos:
dr d
asen a
1
*
dr
sen d
(B).
Poniendo (B) en (A) r
1
*
dr
sen d
1 cos
dr d
rsen 1 cos
.
53
Análisis Matemático
III
Reemplazando en esta expresión
r 2 r
rsen 1 cos
1 r
dr d
r por
sen r 1 cos
d dr
r 2 r
tenemos:
sen r 1 cos
Separando variables: 1 cos sen
* d
dr r
Integrando: dr 1 cos . d sen r lncsc ctg ln sen ln r lnc
de donde: r c 1 cos Es la ecuación de la familia ortogonal.
3.2 F amil ia isogonales. Dos familias se dice que son isogonales a , si cada miembro de la una familia interseca a
todo miembro de la segunda familia a .
Si F x, y, y 0 es la ecuación diferencial de la familia dada, si tan k
F x, y,
y k
0 , es la ecuación diferencial de la familia isogonal.
1 k y
Ej emplo 1. Encuentre la familia isogonal a 60 de la familia de rectas que pasan por el
origen: Solución:
Tenemos: 60 , tan60 3
y c , y
c
y x
y x
De donde, la ecuación de la familia isogonal es: y
y 3 x x 3 y
54
Análisis Matemático
III
Resolviendo (ecuación homogénea), se tiene que la Solución es:
2 3
y 2 2 ln c x y 0 x
tan 1
3.3 Aplicaci ones geomé tr icas. Básicamente las aplicaciones geométricas están relacionadas con la derivada, a continuación pondremos algunas propiedades de ésta relacionadas con las curvas y que nos servirán para realizar los diversos ejemplos. En coordenadas rectangulares.
Ecuación de la tangente a la curva en cualquier punto: si dy/dx es la pendiente de la tangente a la curva en (x,y) tenemos: Y y
dy dx
( X x).
Ecuación de la normal a la curva en cualquier punto: si – dx/dy es la pendiente de la recta normal a la curva en (x,y) tenemos: Y y
dy dx
dx dy
dy dx
dx dy
dy
;
dy dx
;
. respectivamente. dx dy
) 2 : y
) 2 , respectivamente.
Longitudes de la normal entre el punto (x,y) y los ejes X y Y. Son: y 1 ( x 1 (
dx
. respectivamente.
Longitudes de la tangente entre el punto (x, y) y los ejes X y Y. Son: y 1 ( x 1 (
( X x).
Segmentos de la normal al interceptarse con el eje X y Y, son: X x y Y y x
dy
Segmentos de la tangente al interceptarse con el eje X y Y, son: X x y Y y x
dx
dy dx
)2 : y
) 2 , respectivamente.
55
Análisis Matemático
III
Las longitudes de la subtangente y la subnormal son: y
Elemento de longitud de Arco: ds (dx) 2 (dy) 2 .
Elemento de Área: ydx o xdy.
dx dy
y y
dy dx
, respectivamente.
Ej emplo 1. Hallar una curva que tenga un segmento de tangente cuya longitud sea igual a la
distancia desde el punto de contacto hasta el origen de coordenadas
Figura 3.1 Solución: d 1 d 2
Condición del ejemplo: y 1 (
dx dy
) 2 x 2 y 2 , elevando al cuadrado ambos lados de la igualdad y simplificando
tenemos: ( y 1 (
y 2 (
dx dy
dx dy
dx
dy
) 2 ) 2 ( x 2 y 2 ) 2 y 2 1 (
) 2 x 2
dx dy
x y
..(1)..;
dx dy
dx
dy
) 2 x 2 y 2 y 2 y 2 (
) 2 x 2 y 2
x
..(2). y
Como observamos hemos obtenido dos ecuaciones en variables separables. Resolviendo (1) tenemos: dx dy
x y
dy y
dx x
dy y
dx x
ln C ln y ln x ln C ln y ln Cx Y Cx. /// , la
Solución de uno es una familia de rectas.
56
Análisis Matemático
III
Resolviendo (2) tenemos: dx dy
x
y
dy y
dx x
dy y
dx x
ln C ln y ln x ln C ln y ln
C x
Y
C x
. ///
La Solución de (2) es una familia de hipérbolas. Por la tanto las dos familias son las soluciones del ejemplo. Ej emplo 2. Determine las curvas para las cuales la subtangente es igual a la suma de las
abscisas y ordenadas en el punto de tangencia.
Figura 3.2 Solución:
La subtangente se define como la proyección sobre el eje 0 x de la parte de la tangente comprendida entre el punto de tangencia y la intersección de esta con el eje 0 x. Según el gráfico, la subtangente es TK , si es el ángulo de la recta tangente, tenemos: tan y
y OM
De donde: OM
y y y
Por la condición del problema, tenemos:
x y,
y
x
x y y
(Ecuación Homogénea de
primer grado). Con el cambio de variable: z z y
dz dy
zy y y
y
dz dy
x y
y
dx dy
z 1 z y
z y
dz dy
dz dy
tenemos:
1 dz
dy y
57
Análisis Matemático
III
Integrando: dy
dz y
ln C z ln y ln C ln
y C
z
y C
e z y Ce z ///. x
Volviendo a Variables originales obtenemos la Solución: y Ce y ///. Ej emplo 3. Un muchacho está en la esquina de una piscina rectangular, en el otro extremo
que está a 20 m tiene atado a un bote. El muchacho camina a lo largo del lado opuesto al bote sujetando la cuerda que ata al bote, encuentre la trayectoria que sigue el bote y determine la posición del muchacho cuando el bote está a 12 m del borde.
Figura 3.3 Solución:
Sea el ángulo de inclinación de la cuerda, entonces: tan
De donde: dx
400 y 2 y
dy dx
y 400 y 2
dy
Integrando:
dx
400 y 2 y
20 400 y 2 c dy C x 400 y 20 ln y 2
Utilizando las condiciones iniciales: x 0, y 20 obtenemos el valor de C = 0. Así la ecuación que sigue el bote es: 20 400 y 2 . (Ecuación de la trayectoria del bote) x 400 y 20 ln y 2
58
Análisis Matemático
III 20 400 y y
Ahora bien: AE x 400 y 2 AE 20 ln
(I)
Para responder a la pregunta, reemplazamos y 12 en ecuación (I) y encontramos que el muchacho está a 22 m de su posición inicial. Ej emplo 4 . Encuentre la forma que debe tener un espejo de rotación para que un rayo de luz,
que parta de una fuente, se refleje siempre en forma paralela al eje de rotación.
Figura 3.4 Solución:
Ubiquemos la fuente de luz en el origen del plano xy , y el eje de rotación hagámoslo coincidir con el eje Ox . Según el gráfico, si el ángulo SPT es , entonces el ángulo OPA también es , puesto que el ángulo de incidencia es igual al ángulo de reflexión, de donde el triángulo AOP es isósceles, es decir, A O = O P. Además tan y
PM
AM AO OM ,
AM OP x,
AM
PM y,
OP x 2 y 2
De donde la ecuación para esta familia de curvas es: y =
y x x y 2
2
. (Ecuación
Homogénea de grado uno) 59
Análisis Matemático
III
Con el cambio de variable: z z y
dz dy
zy z 2 . y 2 y 2 y
x
y
y
z y
dz dy
dx dy
z y
dz dy
tenemos:
z z 2 1 y
dz dy
z 2 1
dz z 2 1
dy y
Integrando:
dz z 1 2
dy y
ln C ln z z 2 1 ln y ln C ln z z 2 1 ln yC
z z 2 1 yC ///,
Volviendo a variables originales y realizando simplificaciones algebraicas llegamos a:
y 2 c 2 2cx , y 2 2c x
c
que es la ecuación de una familia de parábolas.
2
3.4 Apl icaci ones físicas. 3.4.1 M ecánicas.
Las aplicaciones a la mecánica están relacionadas básicamente con la segunda ley de Newton que dice: F m *a, donde “m” es la masa y “a” la aceleración, cabe recordar también que: a)
a
b) v
dv dt dx dt
d 2 x dt
2
, donde “v” es la velocidad, “t” el tiempo y “x” el espacio y
, donde “t” es el tiempo y “x” el espacio.
, Ej emplo 1. Un paracaidista, está cayendo con una velocidad de 176 pies/seg = 53,65 m seg . 2 cuando se abre su paracaídas. Si la resistencia del aire es Wv 256 lb, donde W es el peso
total del hombre y del paracaídas, hallar su velocidad como una función del tiempo t después de abierto el paracaídas. Solución:
En este caso, tendríamos la siguiente condición: Fuerza neta sobre el sistema = peso del sistema - resistencia del aire. W dv Wv 2 W . 256 g dt
60
Análisis Matemático
III
Simplificando, separando variables e Integrando entre los límites t = 0, v = 176 y t = t y v = v. Tenemos dv
dt 8 v 2 256
ln .
v 16 v 16
1 17 6 v 2 256 8
dv
v
5
v 16
6
v 16
ln . 4t
5
e
1 v 16 dt ln 0 32 v 16
t
4t
6
v 16
6 5e 4t 6 5e 4t
t 8 17 6 v
t
0
/// .
Como conclusión podemos observar que el paracaidista alcanza rápidamente una velocidad aproximadamente constante. Ej emplo 2. Una bala que se mueve con una velocidad de 400 m/s, atraviesa una pared de 20 cm de ancho y sale con una velocidad de 100 m/s. Si la fuerza de resistencia de la pared es
proporcional al cuadrado de la velocidad, determine en que tiempo atraviesa la bala la pared. Solución:
Sean: V 0 = la velocidad inicial V f = la velocidad final h = ancho de la pared
Por la segunda ley de Newton tenemos: m
dv dt
kv 2
el signo negativo lo ponemos para indicar que es fuerza de resistencia. Separando variables e integrando tenemos: dv v2
k m
dt
dv
v
2
Con la condición inicial 1 v
k 1t
1 v0
k m
dt C
1 v
k m
t C // .
de v = v0 para t = 0, encontramos C
... donde.....k 1
1 v0
. de donde:
k m
Si T es el tiempo necesario para que atraviese la pared, se tiene:
1 v1
k 1T
Si x = x (t) es el desplazamiento recorrido hasta el instante t, se tiene: v
1 v0
dx dt
entonces:
61
Análisis Matemático dx dt
III
v0 k 1v 0 t 1
Separando variables e integrando, tenemos: dx
v0 k 1v0 t 1
dt dx
v0
k v t 1
dt C 1 x
1 0
1 k 1
ln( k 1v0 t 1) C 1
Con la condición de que para t = 0, x = 0, encontramos C1 = 0. y con x = h para t = T , Reemplazando se obtiene T
1 1 v0 v1 v0 ln v1 h
Reemplazando los datos del problema: v0 = 400, v1 = 100 y h = 0.2, obtenemos: que el tiempo que necesita la bala para cruzar la pared es de: T = 0.00108 (seg). 3.4.2 El é ctr icas.
Así como la mecánica tiene como base fundamental las leyes de Newton, el tema de la electricidad también tiene una ley que describe el comportamiento de los circuitos eléctricos conocido como Ley de Kirchoff, que dice lo siguiente: La suma algebraica de todas las caídas de voltaje alrededor de un circuito eléctrico es cero (otra manera de enunciar esto es decir que el voltaje suministrado ( f .e .m.) es igual a la suma de las caídas de voltaje). Ej emplo 1.- Un generador con una f e m de 120 voltios se conecta en serie a una resistencia
de 12 Ohmios y un inductor de 3 Henrios. Si el interruptor K se cierra en el tiempo t = 0, establezca una ecuación diferencial para la corriente, y determine la corriente en el tiempo t.
62
Análisis Matemático
III
Figura 3.5 Solución:
Llamando I a la corriente (en amperios) que fluye como se indica en la figura, tenemos: Datos: a.-
Voltaje suministrado = 120 voltios.
b.-
Caída de voltaje a través de la resistencia: R I = 12 I
c.-
Caída de voltaje a través del inductor: L
Por la ley de Kirchoff, se tiene: 3
dI dt
150 sen 20t
dI dt
3
DI dt
3
12I 120
Puesto que el interruptor se cierra en t = 0, debemos tener I (0) = 0. Separando variables, integrando y utilizando la condición inicial se tiene: dI 40 4 I
1
dt .
dI
4 10 I
dt ln C ln 10 I ln C 4t 10 I Ce 4t I 10 Ce 4t
Utilizando la condición inicial encontramos C = 10. Reemplazando tenemos:
I 101 e 4t //.
Ej emplo 2.- Una f e m de cayente E 250e 4t se conecta en serie con una resistencia de 25
ohmios y un condensador de 0.02 Faradios. Asumiendo Q = 0 en t = 0 , encuentre la carga y la corriente en cualquier tiempo, Demuestre que la carga alcanza un máximo, calcúlelo y halle cuando se obtiene.
63
Análisis Matemático
III
Figura 3.6 Solución:
Refiriéndonos a la figura, tenemos:
1.-
Voltaje suministrado: E 250e 4t
2.-
Caída de voltaje en la resistencia: R I = 25 I .
3.-
Caída de voltaje a través del condensador:
Q C
Q 0.02
Por la ley de Kirchoff: 25 I 50Q 250e 25
dQ
dQ dt
dt
4t
y puesto que I
dQ dt
50Q 250e 4t , dividiendo todo para 25,
2Q 10e 4t , esta es una ecuación Lineal en Q, entonces su Solución es:
2 dt 2 dt Qe 10e 4t .e dt C Qe 2t 10 e 4t .e 2t dt C Qe 2t 10 e 2t dt C
Q e 2t . 5e 2t C Q 5e 4t Ce 2t .
Con la condición: Q = 0 en t = 0 encontramos C = 5, de donde: Q 5e 2t 5e 4t ///
Ahora como
I
dQ dt
, se tiene: I 25e 2t 5e 4t 10e 4t I 20e 4t 10e 2t // . que es
la corriente. Para determinar cuándo Q es máxima, debemos hacer
dQ dt
I 0 en la ecuación resultante.
64
Análisis Matemático
III
0 20e 4t 10e 2t 10e 2t 2e 2t 1 0 e 2t
1 2
2t ln1 ln 2 t 0.35 seg .
3.4.3 Quimicas. Ej emplo 1. Dos químicos A y B reaccionan para formar otro químico C. Se encuentra que la
tasa a la cual C se forma varía con las cantidades instantáneas de los químicos A y B presentes. La formación requiere 2 lb de A por cada 1 lb de B. Si 10lb de A y 20lb de B están presentes inicialmente y 10 lb de C se forman en 20 minutos, encuentren la cantidad de químico C en cualquier tiempo. Solución:
Sea X la cantidad de libras de C formadas en el tiempo t (en horas). Luego su formación. Para formar X libras de C se requieren
2 3
X libras de A y
dX dt
1 3
es la tasa de
X libras de B
(por las condiciones del problema). Luego en el instante t cuando se forman X libras de C , la 2
1
3
3
cantidad de A es 10 X y la cantidad de B es 20 X , así: 2 1 k 10 X 20 X , con las condiciones dt 3 3
dX
t = 0, X = 0;
t = 1/3,
X = 10,
resolvemos la ecuación encontrada: Separando variables, integrando y utilizando las condiciones, se tiene: dX 2 1 10 X 20 X 3 3
kdt .
dx
15 x 60 x
2 9
kdt
dX
15 X 60 X
2 9
k dt ln C
B 2 1 1 2 1 A dX kdt ln C ln15 X ln 60 X kdt ln C 9 45 45 9 45 15 x 60 X
60 X 15 X
Ce10kt . ///
Utilizando las condiciones iniciales y despejando X tenemos: 15e 0.17t 1 X . /// 1 0.17 t e 4
65
Análisis Matemático
III
Ej emplo 2. El aire en un cuarto de dimensiones 150*40*10 pie 3 contiene 0,3 % de CO2 ,
aire fresco, conteniendo un 0,04 % de CO2, es admitido a través de un ventilador a una tasa de 6000 (pie3 /min.). ¿Cuál es el porcentaje de CO2 después de 30 minutos? Solución:
Sea X = cantidad de CO2 en pie3 en el cuarto en el instante t. X
Entonces la concentración es:
60.000
.
Durante el intervalo dt la cantidad de CO2 que entra en el cuarto es 6.000 (0.0004) dt ( pie3) y la cantidad que sale es: 6.000
X 60.000
dt ( pie3 )
Entonces, el cambio dX en el intervalo dt es:
dX 6.000 0.0004
X X 24 dt . dt dX 60.000 10
Separando Variables e integrando tenemos: dX X 24
1 10
dt
X 24 Ce
Condiciones:
1 10
t
dX
X 24
10
t ln C X 24 Ce
1 10
t
con t = 0, X = 0.003(60.000) = 180 encontramos C = 156
1 10
t
//.
Cuando t = 20: X 24 156e
60.000
10
dt ln C ln X 24
1
. ///
Por lo tanto: X 24 156e
31.77
1
3
X 31.77. , entonces el porcentaje de CO2 es:
0.0005 0.05% // .
Ej emplo 3. Un tanque de 100 Dl está lleno con salmuera que contiene 60 Kg de sal disuelta,
entra agua en el tanque a una velocidad de 2 Dl por minuto y la mezcla; conservada uniforme mediante agitación, sale a la misma velocidad. ¿Cuánta sal queda en el tanque después de una hora?
66
Análisis Matemático
III
Solución:
Sea S el número de Kg. De sal en le tanque después de t minutos, entonces la concentración es: S 100kg .Dl . Durante el intervalo dt entran en el tanque 2.dt. Dl de agua y salen 2.dt. Dl de mezcla que contienen
2S 100
dt
S 50
dtKg . de sal.
Entonces el cambio dS en la cantidad de sal en el tanque es: dS
S 50
dt .
Separando variables e integrando tenemos: dS S
1 50
dt
dS S
1 50
dt ln C ln S
1 50
t ln C S Ce
1 50
t
.
Con la condición S = 60 para t = 0 encontramos C = 60. Por lo tanto S 60e
1 50
t
// .
Para t = 1 hora = 60 minutos
S 60e
1 50
*60
S 18 Kg . //
3.4.4 Cambio de temperatur a.
La ley de enfriamiento (calentamiento) de Newton establece, que la rapidez de cambio de un cuerpo en cualquier tiempo t, es proporcional a la diferencia de las temperaturas del cuerpo y del medio circundante en el tiempo t. Consideremos: T = temperatura del cuerpo en el tiempo t Tm = temperatura del medio circundante To = temperatura inicial del cuerpo (t = 0). Luego la Ecuación que gobierna este fenómeno es: dT dt
k (T Tm) , ya sea que aumente o disminuya la temperatura, con k como factor de
proporcionalidad.
67
Análisis Matemático
III
Ej emplo 1. Un cuerpo que está a 110 grados centígrados, tarda 12 minutos en enfriarse a 70
grados centígrados en un medio que está a 18 grados centígrados. Después de qué tiempo el cuerpo estará a 30 grados centígrados? Solución: Sea T = T (t) la temperatura del cuerpo en el instante t y T 0 la temperatura del
medio, tenemos las condiciones: Tm =18º, T o = 110, T (12) = 70.
La ecuación diferencial para este ejemplo.- es:
dT dt
k (T Tm)
Resolviendo la Ecuación por Variables separables tenemos: dT T Tm
kdt
dT T Tm
k dt ln C ln T Tm ln C kt
T Tm C
e kt Tm Ce kt // .
Con la condiciones de que T = 110º, t = 0 y Tm = 18º, encontramos el valor de C = 92, reemplazando en la Solución encontrada tenemos: T Tm 92e kt . y con T = 70º, Tm = 18º y t = 12, hallamos el valor de k = -0.05. Con estos valores y tomando T = 30 tenemos: 12 0 18 92e 0.05t 92e 0, 05t 12 0.05t ln t 40,7 min . // . 92 Ej emplo 2. Un termómetro que marca 75º F se lleva fuera donde la temperatura es de 20º F.
Cuatro minutos después el termómetro marca 30º F. Encontrar: a) La lectura del termómetro siete minutos después de que este ha sido llevado al exterior y , b) El tiempo que le toma el termómetro caer desde 75º F hasta más o menos medio grado con respecto a la temperatura del aire. Solución: Sea T = T (t) la temperatura del cuerpo en el instante t y Tm la temperatura del
medio, tenemos las condiciones: Tm =20º F, T o = 75º F, T (4) =30º F.
La ecuación diferencial para este Ejemplo.- es:
dT dt
k (T Tm)
Resolviendo la Ecuación por Variables separables tenemos:
68
Análisis Matemático dT T Tm
kdt
III dT
T Tm
k dt ln C ln T Tm ln C kt
T Tm C
e kt Tm Ce kt // .
Con la condiciones de que T = To = 75º, t = 0 y Tm = 20º, encontramos el valor de C = 55º, reemplazando en la Solución encontrada tenemos: T Tm 55e kt . y con T = 30º, Tm = 20º y t = 4, hallamos el valor de k = -0.426. Reemplazando este valor en la ecuación encontrada tenemos: T Tm 55e 0, 426t /// . a) Con t = 7 minutos, T 20º 55e .426( 7) T 23º F . 0.5 b) 20.5 20 55e 0.426t 55e 0, 426t 0.5 0.426t ln t 11min . // . 55
3.4.5 Descomposición y crecimiento. Ej emplo 1. Un cultivo de bacterias tiene una población inicial Q0 , después de una hora el
número de bacterias medido es
4 3
Q0 . Si la rapidez con la cual crece la población es
proporcional al número de bacterias presente, determine en que tiempo se triplicará el número de bacterias: Solución: Sea Q Qt el número de bacterias en el instante t dQ dt
dQ dt
Es la rapidez de crecimiento, entonces: Q
Como debe existir una constante de proporcionalidad, sea esta k , así la ecuación diferencial para este problema es: dQ dt
kQ
Y la condición inicial es: Q0 Q0
Separando variables e integrando: dQ Q
kdt
dQ Q
k dt ln c ln Q ln c kt Qt ce kt 4
Utilizando las condiciones: Q0 Q0 y Q1,5 Q0 , se determina que: c Q0 3
69
Análisis Matemático
III
4 3 3 2
Y: k ln
Para determinar el tiempo necesario para que se triplique, hacemos: Q 3Q0 2
Reemplazando: 3Q0 Q0 e
3
4 t 3
ln
despejando el tiempo y calculando tenemos: t 5.71h ///.
Que es la Solución buscada. Ej emplo 2. Un líquido transporta droga dentro de un órgano de volumen V (cm3 ) a una tasa 3
3
de “m” (cm /s) y sale a una tasa “n” (cm /s). La concentración de la droga en el líquido que
entra es p (g/cm3 ). Encuentre la ecuación diferencial para la concentración de la droga en el órgano en cualquier tiempo, junto con las condiciones apropiadas y resuelva dicha ecuación. Solución: La situación está descrita en forma breve por la figura, la cual nos muestra un
compartimiento simple de volumen V junto con la entrada y salida.
Figura 3.7 Si X representa la concentración de la droga en el órgano, (esto es, el número de gramos de droga por cm3), entonces la cantidad de droga en el órgano en cualquier tiempo t está dada por: V (cm3 ) X ( gr / cm3 ) XV ( gr )
El número de gr/s que entran el órgano están dado por: m(cm3 / s) * p( gr / cm3 ) mp( gr / s )
El número de gr/s que salen del órgano está dado por: n(cm3 / s) * X ( gr / cm3 ) nX ( gr / s) .
Con la condición que la tasa de cambio de la cantidad de droga en el órgano es igual a la tasa que entra menos la tasa que sale, formamos la ecuación diferencial pedida.. 70
Análisis Matemático d ( XV ) dt
III
mp nX
Si asumimos que la concentración de la droga en el órgano en el tiempo t = 0 es X0, entonces la condición inicial es: X (0) = X 0. Resolviendo por el método de separación de variables:: dX
mp nX
X
mp n
dt V
C n
dX mp nX
1
n
dt ln C ln mp nX ln C t V V
mp nX C
nt
e V
nt
e
V
. ///
Utilizando la condición inicial encontramos C.: C mp nX 0 ,
Reemplazando C en la Solución: nt mp V X 0 X e /// . n n
mp
Ej emplo 3. Se sabe que un cierto material radioactivo decae a una velocidad proporcional a su
cantidad de material presente, Un bloque de ese material tiene originalmente una masa de 120 gramos y cuando se le observa después de 25 años, su masa ha disminuido a 100 gramos. Encuentre una expresión para la masa de ese material como función del tiempo. Solución:
Sea x (t) = cantidad de sustancia radioactiva en cualquier tiempo (t) la descripción matemática es:
dx(t ) dt
kx(t ).
Resolviendo la ecuación por variables separables tenemos: dx(t ) x(t )
kdt
dx(t )
x(t )
k dt ln C ln x(t ) kt ln C ln
Integrando: x(t ) C
kt x(t ) Ce kt /// .
Para determinar el valor de C tenemos la siguiente condición: para t = 0, x (t) = 120 gr. De donde
C = 120. Luego reemplazando en la Solución encontrada tenemos:
x(t ) 120e kt // .
71
Análisis Matemático
III
Ahora para determinar el valor de k tenemos que t = 25 años, y x (t) = 100gr. De donde: 100 120e 25k k
1
12 ln( ) // . 25 10
Luego la ecuación buscada quedaría así: x(t ) 120e
1 20
12 ln( ) t 10
/// .
3.5 Ej ercicios propuestos . Hallar las trayectorias ortogonales para las siguientes familias de curvas: 1. y 2 2ax a 2, a es un parámetro 2. y ae ax , a const . 1
3. x 2 y 2 a 2 2
4. x k y k a k 1
5. x 2 y 2 a 2 3
6. y 2 4 x a
7. Encontrar las trayectorias ortogonales de la familia de elipses con centro en (0,0) y de vértices en (1,0) y (-1,0). 1
8. Encontrar las trayectorias ortogonales de la familia de curvas y x tan g ( y k ). 2
9. La temperatura de una placa delgada está dad por T ( x, y) e y cos x. Encontrar la ecuación de las líneas de flujo de calor.
10. Hallar las trayectorias ortogonales de la familia de curvas que satisfacen la siguiente propiedad: “La recta tangente a una de las curvas en un punto cualquiera P, es la bisectriz del ángulo determinado por la recta vertical que pasa por P y la recta que une P con el origen de coordenadas.
11. Hallar las trayectorias isógona les de la familia de parábolas y 2 4ax, el ángulo formado es α = 45º.
12. Hallar la familia de trayectorias que cortan las líneas x 2 2a( y x 3), formándose un ángulo α = 60º.
72
Análisis Matemático
III
En los ejercicios dados a continuación, primero encontrar la ecuación diferencial de acuerdo a las condiciones de los ejemplos y luego hallar su solución. 13. Hallar una curva que pase por el punto 0, 2 de modo que el coeficiente angular de la tangente en cualquiera de sus puntos sea igual a la ordenada del mismo punto, aumentada tres veces.
14. Hallar una curva que posea la propiedad de que la magnitud de la perpendicular bajada del origen de coordenadas a la tangente sea igual a la abscisa del punto de contacto.
15 Hallar la curva que tiene la propiedad de que el segmento de la tangente a la curva comprendida entre los ejes de coordenadas se divide por la mitad en el punto de contacto.
16. Un barco retrasa su movimiento por la acción de la resistencia del agua, que es proporcional a la velocidad del barco. La velocidad inicial del barco es 12 m/sg, después de 7sg su velocidad será de 9 m/sg. ¿Después de cuánto tiempo su velocidad será de 2 m/sg?
17. Determinar la trayectoria S recorrida por un cuerpo durante el tiempo t, si su velocidad es proporcional al trayecto, sabiendo que en 15 sg el cuerpo recorre 120 m y en 18 sg 200m.
18. Un punto material de masa igual a 1 g se mueve en línea recta debido a la acción de una fuerza que es directamente proporcional al tiempo, calculado desde el instante t 0 , e inversamente proporcional a la velocidad del punto. En el instante t 10 s la velocidad era igual a 50 cm/s, y la fuerza, igual a 4 dinas. ¿Qué velocidad tendrá el punto al cabo de un minuto del comienzo del movimiento?
19. Una bala se introduce en una tabla de h 10 cm de espesor con la velocidad v0 200 m/s traspasándola con la velocidad v1 80 m/s. Suponiendo que la resistencia de la tabla al movimiento de la bala es proporcional al cuadrado de la velocidad, hallar el tiempo del movimiento de la bala por la tabla.
20. Demostrar que la curva para la cual la pendiente de la tangente en cualquier punto es proporcional a la abscisa del punto de contacto, es una parábola.
21. Hallar la curva para la cual la pendiente de la tangente en cualquier punto es n veces mayor que la pendiente de la recta que une este punto con el origen de coordenadas.
22. El fondo de un depósito, de 300 litros de capacidad, está cubierto de sal. Suponiendo que la velocidad con que se disuelve la sal es proporcional a la diferencia entre la 73
Análisis Matemático
III
concentración en el instante dado y la concentración de la disolución saturada (1 kg de sal para tres litros de agua) y que la cantidad de agua pura dada disuelve 1/3 de kg de sal por min., hallar la cantidad de sal que contendrá la disolución al cabo de una hora.
23. Cierta cantidad de una substancia indisoluble que contiene en sus poros 2 kg de sal se somete a la acción de 30 litros de agua. Después de 5 min. Se disuelve 1 kg de sal. ¿Dentro de cuánto tiempo se disolverá el 99% de la cantidad inicial de sal?
24. Si la temperatura del aire es de 25º C y el cuerpo se enfría en 17 minutos desde 110º C hasta 75º C. ¿Dentro de que tiempo su temperatura descenderá hasta 45º C?
25. El radiactivo tiene una vida promedio de 5600 años aproximadamente. ¿En cuántos años desciende el 20% de su cantidad original? ¿Al 10%? Rpta: t
5600 Ln(0.20) Ln(2)
26. Cierta cantidad de sustancia, que contenía 3 kgrs. De humedad, se colocó en una habitación de 100m3 de volumen donde el aire tenía al principio el 25% de humedad. El aire saturado, a esta temperatura, contiene 0.12 kgr. De humedad por 1m3, si durante el primer día la sustancia perdió la mitad de su humedad. ¿Qué cantidad de humedad quedará al finalizar el segundo día? Rpta: 0.82 Kg. Sug.
ds dt
kss 6
27. El Presidente y el primer Ministro piden café y reciben tazas a igual temperatura y al mismo tiempo. El Presidente agrega inmediatamente una pequeña cantidad de crema fría; pero no se toma café hasta 10 min. Después. El primer Ministro espera 10 min. Y, luego añade la misma cantidad de crema fría y comienza a tomarse su café. ¿Quién tomará el café más caliente?
Rpta. El Presidente
28. Si un circuito eléctrico contiene una resistencia R (ohmios) y un condensador de capacidad C (faradios) en serie y una f.e.m. E (voltios), a la carga del condensador q está dada por: R
dq dt
q C
E , encontrar q(t).
29. La ecuación diferencial de un circuito que contiene una resistencia R, capacidad C y f.e.m. e E . sen.(wt ) es R
di dt
i C
de dt
Suponiendo constantes R, C, E, W, hallar la
corriente i en el instante t.
74
Análisis Matemático
III
30. Deducir para el sistema indicado en la figura, tres ecuaciones aplicando las leyes de Kirchoff, supóngase que la carga del condensador es nula cuando t = 0 y deduzca que siempre i2 = 2 amperios y que i1
1 10
e 10t , tendiendo rápidamente a cero.
Figura 3.8 31. Hallar i en función del tiempo t en el sistema representado en la figura indicada, si son cero todas las corrientes iniciales y la carga del condensador.
Figura 3.9 32. Empleando coordenadas rectangulares, hallar la forma del espejo si los rayos que parten de un punto dado, al reflejarse, son paralelos a una dirección dada
75
Análisis Matemático
III
UNIDAD IV 4. Ecuaciones diferenciales de orden superior Las ecuaciones diferenciales de orden superior o de n-ésimo orden tienen la forma: F ( x, y, y , y ,.........., y n ) 0 Vale recordar también que las soluciones de estas ecuaciones pueden ser generales las mismas que deben tener un número igual de constantes arbitrarias al orden de la ecuación y las particulares dependiendo si hay o no condiciones iniciales o de frontera, éstas también deben ser en número igual al orden de la ecuación dada. Así por ejemplo: Si tenemos una ecuación de tercer orden deben existir tres constantes arbitrarias en caso de las generales y tres condiciones iniciales para las particulares. y x ln( x), y(1) y (1) y (1) 0 o lo que quiere decir lo mismo que y, y`, y y``, son
iguales a cero (0), cuando x = 1.
4.1 Caso de integración inmediata. Si tenemos la ecuación diferencial y ( n) f ( x). Donde f(x) sólo es función de “x”. Para encontrar la Solución general de este tipo de ecuación tenemos que integrar n veces. Así:
y ..... f ( x)dx......dx C 1
x
n 1
(n 1)!
C 2
x
n 2
(n 2)!
.........C n1 x C n .
Ej emplo 1. Hallar la Solución: general de la siguiente ecuación. y xe x .
Solución: Integrando sucesivamente la ecuación dada tenemos: y xe x e x C 1 y ( x 1)e x C 1 y ( x 1)e x e x C 1 x C 2 y ( x 2)e x C 1 x C 2 . y ( x 2)e e C 1 x
x
x 2 2
C 2 x C 3 y ( x 3)e C 1 x
x 2 2
C 2 x C 3 /// .
Ejemplo 2. Hallar la Solución particular de la siguiente ecuación. y x ln x. Con las condiciones: y(1) y (1) y (1) 0
Antes de encontrar la Solución indicada vale recordar que para encontrar una Solución particular, es necesario tener condiciones iniciales, las mismas que tienen que ser igual en número al orden de dicha ecuación. Así en nuestro ejemplo tenemos una ecuación de 3er orden, por lo que el número de condiciones iniciales tienen que ser 3 lo cual observamos. Solución: Integrando sucesivamente la ecuación dada tenemos: 1 1 y x 2 ln x x 2 C 1 . 2 4 1
Con la condición de que y 0 , cuando x = 1, encontramos que C 1 . Por lo que 4
76
Análisis Matemático
III
1 1 1 y x 2 ln x x 2 . 2 4 4
Integrando por segunda vez tenemos: 1 1 1 1 1 5 1 y x 3 ln x x 3 x 3 x C 2 . y x 3 (ln x ) x C 2 . 6 18 12 4 6 6 4 1
Ahora con la condición de que y 0 , cuando x = 1, encontramos que C 2 . Por lo que: 9
1 5 1 1 y x 3 (ln x ) x . 6 6 4 9
Integrando otra vez más:
5 1 1 1 1 4 13 1 1 x 4 (ln x ) x 4 x 2 x C 3 y x (ln x ) x 2 x C 3 . 24 6 96 8 9 24 12 8 9 1 Por último con la condición de que y 0 , cuando x = 1, encontramos que C 3 . Por lo 32 1 13 1 1 1 que: y x 4 (ln x ) x 2 x . ///. Que es la Solución pedida. 24 12 8 9 32 y
1
4.1.1 Ecuación de la f orm a: d 2 y g ( y ). Donde g(y) es función de “y” únicamente. dx 2
Para encontrar su Solución procedemos de la siguiente manera: Ponemos la ecuación de la siguiente manera: d y g ( y )dx. Multiplicamos ambos miembros por y , y d y g ( y) y dx. y g ( y )dy. que es una ecuación en Como y dx dy , la ecuación se convierte en: y d variables separables. 1
Integrando tenemos: y 2 g ( y)dy C 1 . Donde el 2do miembro es un función de y
únicamente. Extrayendo la raíz cuadrada, las variables x y quedan separadas y podremos integrar otra vez. Encontrando así su Solución general.
2
Ej emplo 1. Resolver:
Solución:
d 2 y 2
d 2 y 2
d x
a2 y 0
a 2 y Multiplicando todo por y dx .
d x y d y a y. y dx ; recordando que y dx dy , tenemos: y d y a 2 y.dy. 2
Integrando:
y d y
1
a y.dy. y 2
2
2
a 2 y 2 2
C y 2 a 2 y 2 C 1 . con C1= 2C.
77
Análisis Matemático
III
Despejando y´ tenemos una ecuación en variables separables, por lo que luego de integrar y haciendo algunas simplificaciones tenemos: dy
y C 1 a 2 y 2
1 a
ay
arc. sen
C 1 a y 2
2
x C 2 y
C 1
dy
dx C 1 a
C 1 a y 2
2
dx C 2 .
sen(ax aC 2 ).
Utilizando la identidad trigonométrica del seno de la suma de ángulos y realizando algunas sustituciones tenemos: y c1 senax c2 cos ax /// . que es la solución. Ej emplo 2. Resolver: d 2 y
Solución: y d y
b y
2
d x
d 2 y d 2 x
b y
3
b y 3
.
Multiplicando todo por y dx .
. y dx ; recordando que y dx dy , tenemos: y d y
3
Integrando:
b
yd y y
3
dy C
1 2
y
a
2
2 y 2
b y 3
C y 2
dy.
C 1 y 2 a y 2
. con
C1= 2C.
Despejando y´ tenemos una ecuación en variables separables, por lo que luego de integrar y haciendo algunas simplificaciones tenemos: y
1 y
1
C 1 y 2 a
ydy C 1 y a 2
dx
ydy C 1 y a 2
dx C 2 .
C 1 y 2 a x C 2 C 1 y 2 a C 1 x C 1C 2 C 1 y 2 a (C 1 x C 3 ) /// .
C 1
C1.C2 = C3. 4.1.2 E jercici os propuestos:
1. 2. 3. 4.
y xe x , y(0) y (0) y (0) 0 .
y IV x 1 y x ln( x), y(1) y (1) y (1) 0 y x senx x 5. y ( x 1) 4
6. 7.
d 2 x 2
dt d 2 x dt 2
t 2 3 t
78
Análisis Matemático 8. y 9. y d 2 s
10.
2
dt
d 2 s
11.
2
dt
d 2 s
12.
2
dt
III
b y
3
d
0
y 2
1
( s 2) 3
1
as
k s 2
, Hallar “t” , si para t = 0 es
s = a y ds/dt = 0.
4.2 Ecuaciones que se r educen a un orden in ferior. 4.2.1 Ecuación diferencial que no contiene “y” de forma explicita.
Esta ecuación es de la forma: F ( x, y , y , y ,......) 0 en la que realizando la sustitución de p y
dy dx
, p y , etc. obtenemos una ecuación de orden inferior en una unidad la misma
que quedaría así: F ( x, p, p, p,......) 0 . Ej emplo 1. Resolver: x y y 0
Solución: Sustituyendo p y ; p
dp dx
y , en la ecuación dada tenemos: x
Despegando variables e integrando:
dp p
dx x
0
dp p
dx x
ln( p) ln( x) ln(C 1 ) ln( px) ln(C 1 ) px C 1 p
Como p y
dy dx
dp dx
p0
0
C 1 x
. ///.
, sustituyendo en la Solución anterior, despegando e integrando obtenemos
la Solución general de la ecuación dada. dy dx
C 1 x
dy C 1
dx x
dy C 1
dx x
C 2 y C 1 ln( x) C 2 /// .
Ej emplo 2. Encontrar la Solución particular de: x y 1 y 2 con las condiciones de que: y 0 para x = 1; y = 1 para x = e2.
79
Análisis Matemático
III
Solución: Sustituyendo p y ; p
dp dx
y , en la ecuación dada tenemos: x
Despegando variables e integrando:
dp 1 p
2
dx x
dp 1 p
2
dx x
dp dx
1 p2
ln(C 1 )
ln p 1 p 2 ln( x) ln(C 1 ) ln p 1 p 2 ln C 1 x p 1 p 2 C 1 x . 1 p 2 C 1 x p 1 p 2 C 1 x 2 2C 1 xp p 2 C 1 x 2 2C 1 xp 1. ///. 2
2
Como y p 0 para x = 1, tenemos que x 2 2 xp 1 /// . Como p y obtenemos: x 2 x 2
dy
dy dx
C1 = 1, por lo que la ecuación queda:
, sustituyendo en la Solución anterior, despejando variables e integrando
1 2dy
x 2 1
dx 2 dy
dx x 1 y x 2 ln x C 2 /// . (a) 4
x 2 1 x
1 dx C 2 2 y x 2 ln( x) C 2 . 2
2
Con la condición de que y = 1 para x = e , tenemos que C 2 1
encontramos lo solicitado. y x ln x 4
2
8 e4 4
8 e4 4
. Reemplazando en (a)
/// .
4.2.2 Ecuación diferencial que no contiene “x” de forma expl ícita.
Esta ecuación tiene la forma: F ( y, y , y , y ,.....) 0 realizando la siguiente sustitución obtenemos una ecuación de orden inferior en una unidad. Sustitución: y p, y p p y p 2 p p p 2 etc. Con p
dp dy
, la ecuación dada queda
de la siguiente manera: F ( y, p, p , p 2 p pp 2 ,......) 0 . Ej emplo 1. Encontrar la Solución: de y y y 2 1 para las condiciones de que: y = 1; y´ = 1
para x = 0. Solución:
Como la ecuación dada es de orden dos la sustitución es: p y ; p
dp dy
y .
80
Análisis Matemático
III
Realizando la sustitución, separando variables, integrando y simplificando tenemos: dp
y. p
dy
p 2 1. y. p
1
dp dy
p 2 1 0
pdp p 1 2
ln p 1 ln( y) ln(c) ln p 1 ln 2
2
2
c2 y 2
dy y
pdp p 1
p 1 2
2
c2 y 2
dy y
p 2
0
c 2 y 2 y 2
/// .
Utilizando la condición p = y´= 1; para y = 1, encontramos el valor de c = 0, reemplazando en la ecuación anterior tenemos: p2 = 1, p = 1 //. Reemplazando p dy dx
dy dx
, separando variables e integrando tenemos:
1 dy dx dy dx C 1 y x C 1 /// .
Por último con la condición y = 1, para x = 0, encontramos que C1 = 1. Reemplazando, tenemos la Solución pedida y = x + 1///. Ej emplo 2. Encontrar la Solución: de y 2 y y y 2 y . Solución:
Como la ecuación dada es de orden dos la sustitución es: p y ; p
dp dy
y .
Realizando la sustitución, y simplificando tenemos: p 2 yp
dp dy
y 2 p y
dp dy
y 2 p y
dp dy
p y 2
(a) es una ecuación lineal en “p”, donde p( y )
1 y
dp dy
1
p y. (a). y
y Q( y) y por lo que su Solución
general en fórmula es: (b) p.u u.Q( y)dy C 1 . Con u e En nuestro caso u e
dy
y
p ( y ) dy
1
.
u e ln y u
1 y
. Reemplazando todo en (b), e integrando
tenemos: p.
1 y
1
. y.dy C 1 p y( dy C 1 ), p y( y C 1 ). ///. y
81
Análisis Matemático
III
Para terminar el Ejemplo.- reemplazamos a p y ; luego de separar variables e integrar llegamos a la Solución solicitada. dy dx
y( y C 1 )
dy y( y C 1 )
dx
A B dy x C 2 . dx C 2 y( y C 1 ) y y C 1 dy
1 1 1 y dy x C 2 ln y ln y C 1 x C 2 ln x C 2 . C y C ( y C ) C y C 1 1 1 1 1
1
y y C 1 .e x C /// . C 1
2
4.2.3 Ej ercici os propuestos.
1. y 2 5 y 6 0 y 2. x y y ln( ) x 3. y y 2e y 2
4. x y 2 1 y 2 2 5. y y y x 6. y 1 y 2
2
7. x 2 y xy 1 8. ( x 1) y ( x 2) y x 2 0 9. y y y 2 y 2 y y y 10. y y y 2 y y 2 y 2 11. y y y 2 y 3 ln( y) 0 12. x y y 0 Hallar las soluciones que satisfagan a las condiciones que se indican: 3 1,......; y 1,... y 1 , para x = ½. 13. y y
14. x y 1 y 2 ;........ y 0,... para... x 1;.... y 1,.... para.... x e 2 y y 15. y (1 ln );.... y 1 / 2,... y 1,.. para...x 1 x
x
82
Análisis Matemático
III
4.3 Ecuaciones dif erenciales li neales de orden superi or. 4.3.1 Definiciones.
Una ecuación de la forma:
y ( n ) p1 ( x) y n1 .......... .. pn ( x) y Q( x). (a),
donde
p1 ( x),....., pn ( x) , son funciones de “x” solamente o constantes se llama ecuación diferencial Lineal de orden “n”. Si en (a) q( x) 0 , la ecuación lineal toma el nombre de Homogénea.
Antes de ver la manera de resolver estas ecuaciones es necesario recordar las definiciones de independencia y dependencia lineal de funciones.
Dependencia lineal.- Las funciones y1 f 1 x , y 2 f 2 x ,....... y n f n x se llaman 2 ,...., C n , tales, que sin ser linealmente dependientes, cuando existen unas constantes C 1 , C todas iguales a cero, se tiene C 1 y1 C 2 y 2 ..... C n y n 0; En el caso contrario estas funciones reciben el nombre de linealmente independientes. Una de las formas de comprobar la dependencia e independencia lineal de las funciones es mediante un determinante especial llamado WRONSKIANO el mismo que tiene la siguiente definición.
W ( y1 , y 2 ,......, y n )
y1
y 2
y3 ...
y n
y1
y 2
y3........
y n
........
.........
y 3n 1 ....
y nn 1
....... ....... y1n 1
y 2n 1
Como podemos observar en su definición será necesario que cada función tenga al menos su derivada de orden (n-1) en un intervalo Al resolver el wronskiano vamos a obtener dos soluciones la una nula y la otra diferente de cero, las mismas que nos determinaran si las funciones analizadas son linealmente dependientes o independientes respectivamente. Nota: Que el Wronskiano W 0 , para que las funciones sean linealmente independientes es una condición necesaria pero no suficiente. Ejemplos : Mediante el Wronskiano demostrar que el siguiente conjunto de funciones son
linealmente independientes.
1.- 1, e x ,2e 2 x . Solución: 1 W 0 0
e x
e x e x
2e 2 x 4e 2 x 8e x 4e x 12e x /// . 8e 2 x
83
Análisis Matemático
III
Como el w 0 estas funciones son linealmente independientes.
2.- x 2 , x 4 , x 8 . Solución: x 2
x 4
x 8
W 2 x
4 x 3
8 x 7 (224 x11 24 x11 16 x11 ) (8 x11 96 x11 112 x11 ) 48 x11 /// .
2
12 x 2
56 x 6
Como el w 0 estas funciones son linealmente independientes.
3.- x,2 x,4 x. Solución: x
2 x
4 x
W 1
2
4 0 /// .
0
0
0
Como el w 0 estas funciones son linealmente dependientes. 4.3.2 Ecuaci ones homogé neas.
Al haber ya, recordado los conceptos de dependencia e independencia lineal de funciones podemos decir que la Solución general de la ecuación diferencial lineal homogénea de orden n (1) y n P 1 x y n1 .... P n x y 0 Con coeficientes continuos P i x i 1,2,...., n tiene la forma y C 1 y1 C 2 y 2 .... C n y n , Donde y1 , y 2 ,... y n , son soluciones linealmente independientes de la ecuación (1) (también se llama sistema fundamental de soluciones).
4.3.3
Ecuaci ones no homogé neas .
La Solución general de la ecuación diferencial lineal no homogénea y n P 1 x y n1 .... P n x y Q x ,
(2)
Siendo los coeficientes P i x y el segundo miembro Q x funciones continuas, tiene la forma
Donde es la solución general de la correspondiente ecuación homogénea (1) e una solución particular de la ecuación no homogénea dada. (2). Si se conoce un sistema fundamental de soluciones y1 , y 2 ,...., y n de su correspondiente ecuación homogénea (1), la Solución general de la ecuación no homogénea (2) se puede hallar por la fórmula 84
Análisis Matemático
III y C 1 x y1 C 2 x y 2 ....C n x y n ,
Donde las funciones C i x i 1,2,...., n se obtienen del sistema de ecuaciones C 1 x y1 C 2 x y 2 ........ C n x y n 0,
C 1 x y1 C 2 x y 2 ........ C n x y n 0 .................................................................................... (3) .................................................................................... n 2 n 2 n 2 C 1 x y1 C 2 x y 2 ........ C n x y n 0 n 1 n 1 n 1 C 2 x y 2 ........ C n x y n Q x C 1 x y1 (Este procedimiento se llama método de variación de las constantes arbitrarias).
4.4 Ecuaciones diferenci al es lineales homogé neas con coeficientes constantes. Las ecuaciones diferenciales lineales homogéneas con coeficientes constantes son de la forma: bn y n bn1 y n 1 .......... . b1 y b0 y 0 , donde b0, b1,…..bn son constantes. Para resolver este tipo de ecuaciones diferenciales primero consideremos el polinomio característico de la siguiente forma: P (m) bn m n bn1m n 1 .......... . b1m b0 0
Como el polinomio característico P (m) = 0 es de grado n, entonces al encontrar su Solución puede obtenerse n raíces las mismas que pueden ser: reales distintas, reales iguales con multiplicidad s o imaginarias, las que nos permitirán dar diferentes formas a la Solución de esta ecuación diferencial.
1º Caso.- Cuando las raíces del polinomio característico P (m) 0 son reales distintas, el sistema fundamental de soluciones tendrá la forma siguiente: m x m x m x e , e ,......., e , y la Solución general de la ecuación diferencial homogénea es: m x m x m x y g c1e c2 e ....... cn e , ///. 1
2
n
1
2
n
2º Caso.- Cuando las raíces del polinomio característico P (m) 0 son reales y alguna de multiplicidad s, consideremos m1 m2 ...... m s m , donde m es la raíz de multiplicidad s, y n – s son las demás raíces y distintas, el sistema fundamental de soluciones tendrá la forma siguiente: e m1 x , xe mx , x 2 e mx ......., x s 1e mx , e m s 1 x ,.....e mn x .
En este caso la Solución general de la ecuación diferencial homogénea será: y g c1e mx c2 xe mx c3 x 2 e mx ....... c s x s 1e mx c s 1e m x ..... cn e m x /// . s 1
n
3º Caso.- Cuando las raíces del polinomio característico P (m) 0 alguno de estas son imaginarias: m1 a1 ib1 ; m2 a1 ib1 ; m3 a 2 ib2 ; m4 a2 ib2 , y las demás raíces supongamos que sean reales y distintas, luego el sistema fundamental de soluciones tendrá la forma siguiente: 85
Análisis Matemático
III
e a1 x cos .b1 x; e a1 x sen.b1 x; e a2 x cos .b2 x; e a2 x sen.b2 x; e m5 x ,.......e mn x .
Por lo que en este caso la Solución general de la ecuación diferencial homogénea vendrá dada por: y g c1e a x cos .b1 x c2 e a x sen.b1 x c3 e a x cos .b2 x c4 e a x sen.b2 x c5 e m x ...... cn e m x /// . 1
1
2
2
5
n
Ejemplos: Resolver las siguientes ecuaciones diferenciales:
1.- y 3 y 2 y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 2 3m 2 0 resolviendo encontramos sus raíces m1 1; y m2 = 2; por lo que su sistema fundamental de soluciones es: ex, e2x, y la Solución general y g c1e x c2 e 2 x /// .
2.- y 3 y 6 y 8 y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: m2 = -4 y m3 P (m) m 3 3m 2 6m 8 0 ; resolviendo encontramos sus raíces m1 2; = -1; por lo que su sistema fundamental de soluciones es: e2x, e-4x, y e-x, y la Solución general y g c1e 2 x c2 e 4 x c3 e x /// .
3.- y y 16 y 20 y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 3 m 2 16m 20 0 ; resolviendo encontramos sus raíces m = 2 de multiplicidad 2 y m1 = - 5; por lo que su sistema fundamental de soluciones es: e2x, xe2x, y e-5x, y la Solución general y g c1e 2 x c2 xe 2 x c3 e 5 x /// .
4.- y 10 y 25 y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 2 10m 25 0 ; resolviendo encontramos sus raíces m = 5 de multiplicidad 2 ; por lo que su sistema fundamental de soluciones es: e5x, xe5x, y la Solución:general y g c1e 5 x c 2 xe5 x /// .
86
Análisis Matemático
III
5.- y y 5 y 7 y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 3 m 2 5m 7 0 ; resolviendo encontramos sus raíces m1 = -1; m2 1 6i y m3 1 6i ; por lo que su sistema fundamental de soluciones es: e-x, e x cos . 6 x; e x sen. 6 x , Y la Solución general y g c1e x c2 e x cos . 6 x c3 e x sen. 6 x /// .
6.- y1V 2 y y 0 Solución: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 4 2m 2 1 0 ; resolviendo encontramos sus raíces m1 = i; de multiplicidad 2 y m2 i , de multiplicidad 2; por lo que su sistema fundamental de soluciones es: cos x; x.cos x; sen x; x.sen x, Y la Solución general y g c1 cos x c2 x cos x c3 sen. x c4 xsenx /// .
4.5 Ecuaciones diferenciales lineales no h omogé neas con coefi cientes constantes. Las Ecuaciones diferenciales lineales no homogéneas son de la forma: bn y n bn1 y n 1 .......... . b1 y b0 y Q( x) (A),
donde b0, b1,…..bn son constantes reales Para hallar la Solución general de la ecuación diferencial lineal no homogénea (A), primero se determina la Solución general de su respectiva ecuación diferencial lineal homogénea yg, y luego se busca una Solución particular cualquiera de la ecuación diferencial no homogénea y p, y la Solución general de la ecuación (A) es la suma de las dos. y y g y p .
Para determinar la función y p existen dos procedimientos: el uno por medio del método de los coeficientes indeterminados y el otro por la variación de la constante arbitraria o del parámetro.
4.5.1 Casos donde se puede encontrar indeterminados.
, uti lizando el mé todo de los coeficientes
1er Caso. Cuando el segundo miembro de la ecuación (A) es la función Q x P n ( x) , donde Pn(x) es un polinomio de grado n entonces:
87
Análisis Matemático
III
a) Si m 0 , no es raíz del polinomio característico P (m) 0 , la Solución particular es de la forma Y p Rn ( x) , donde R n(x) es un polinomio grado n con coeficientes indeterminados. b) Si m 0 , es raíz del polinomio característico P (m) 0 , entonces la Solución particular es de la forma Y p x s Rn ( x) , donde s es la multiplicidad de m 0 .
2do Caso. Si la función del segundo miembro de (A) es: Q x e ax Pn x donde Pn x es un polinomio de grado n, y a es real, tenemos: a) Si a no es raíz de la ecuación característica P (m) 0 , entonces su Solución particular es: y p e ax Rn x donde Rn x es un polinomio de grado n con coeficientes indeterminados. b) Si a es raíz de la ecuación característica P (m) 0 , se considera como Solución particular y p x s e ax Rn x , donde s es el grado de multiplicidad de la raíz a.
3er
Caso.
Si
tenemos que el segundo miembro de la ecuación (A) es Q x e P n x cos bx S m x senbx, donde Pn(x) y Sm(x) son polinomios de grado n y m respectivamente, entonces: ax
a) Si a bi , no es raíz de la ecuación característica P (m) 0 , su solución particular es: Y p e ax T N x cos cos bx U N x senbx senbx, Donde T N x y U N x son polinomios de grado N =máx =máx n, m. b) Si por el contrario, a bi , es raíz de la ecuación característica P (m) 0 , entonces la Solución particular sería: Y p x s e ax T N x cos bx U N x senbx senbx, Donde s es el grado de multiplicidad de la raíz a bi . (Para la ecuación de 2° orden s = 1 En el caso general, para resolver la ecuación las ecuaciones diferenciales lineales de cualquier orden se emplea el método de variación de las constantes arbitrarias. Hallar las soluciones de las siguientes ecuaciones: Ejemplos: Hallar
1.- y 4 y 4 y x 2 Solución: El polinomio característico, de la ecuación diferencial lineal homogénea correspondiente a la ecuación dada es:
88
Análisis Matemático
III
P (m) m 2 4m 4 0 ; resolviendo encontramos sus raíces m1 = 2; de multiplicidad 2; por
lo que su sistema fundamental de soluciones es: e 2 x ; xe 2 x y su Solución es: y g c1e 2 x c2 xe 2 x /// . Y de acuerdo a la parte a del 1er caso la Solución particular es: y p Ax 2 Bx C . Para determinar los valores de A, B y C, formamos un sistema de ecuaciones derivando esta Solución cuantas veces sean necesarias y reemplazando en la ecuación dada así: y p Ax 2 Bx C y p 2 Ax B; y y p 2 A .
Como y 4 y 4 y x 2 , es la ecuación dada tenemos: 2 A 4(2 Ax B) 4( Ax 2 Bx C ) x 2
4 A 1 8 A 4 B 0 2 A 4 B 4C 0
Resolviendo el sistema y reemplazando los valores en la Solución particular tenemos: A
1 4
; B
1 2
; y..C
3 8
1
1
4
2
y p x 2
x
3 8
/// /// .
Luego la Solución general de la ecuación no homogénea es: 1 1 3 y c1e 2 x c 2 xe 2 x x 2 x /// /// . 4 2 8
2.- y iv 3 y 4 y 2 x 4 8x 2 Solución: Sg. Y = Yg + Y p Para Hallar Yg
El polinomio característico, de la ecuación diferencial lineal homogénea ASOCIADA a la ecuación dada es: P (m) m 4 3m 2 4 0 ; resolviendo encontramos sus raíces m1 = 2; m2 = -2; m3 = i y m4 = -i; por lo que la Solución es: y g c1e 2 x c2 e 2 x c3 cos cos x c4 senx /// /// . Y de acuerdo a la parte y p Ax 4 Bx 3 Cx 2 Dx E .
“a”
del
1er
caso
la
solución
particular
es:
Para determinar los valores de A, B, C, D, y E, formamos un sistema de ecuaciones derivando esta Solución cuantas veces sean necesarias y reemplazando en la ecuación dada así:
89
Análisis Matemático
III
y p Ax 4 Bx 3 Cx 2 Dx E y p 4 Ax 3 3 Bx 2 2Cx D; y p 12 Ax 2 6 Bx 2CX
.
y p 24 Ax 6 B; y iv 24 A
Como y iv 3 y 4 y 2 x 4 8x 2 , es la ecuación dada tenemos: 24 A 3(12 Ax 2 6 Bx 2C ) 4( Ax 4 Bx 3 Cx 2 Dx E ) 2 x 4 8 x 2
4 A 2 4 B 0 36 A 4C 0 18 B 4 D 8 24 A 6C 4 E 2
Resolviendo el sistema y reemplazando los valores en la Solución particular tenemos: 1 9 1 9 A ; B 0; C ; D 2; y.. E 8 y p x 4 x 2 2 x 8 /// /// . 2 2 2 2
Luego la Solución general de la ecuación no homogénea es: 1 9 y c1e 2 x c2 e 2 x c3 cos x c4 senx x 4 x 2 2 x 8 //// . 2 2
4.5.2 Ejemplos de ecuaciones diferenciales lineales de orden “n” con coeficientes constan con stantes tes no h omogé omo gé n eas, eas, uti u ti l i zando el m é todo tod o de l a var i ació aci ón del par ámetr o.
Hallar las soluciones de las siguientes ecuaciones: Ejemplos: Hallar
1.- y y csc x Solución: a.- Encontramos la Solución de su respectiva ecuación homogénea. y y 0
Polinomio característico: p(m) 0 m 2 1 0 m i y g C 1 cos x C 2 senx. 90
Análisis Matemático
III
b.- La Solución particular de la ecuación diferencial viene dada por: C 1 ( x) cos x C 2 ( x) senx 0 ; y p C 1 ( x) cos x C 2 ( x) senx, Tal que C ( x ) senx C ( x ) cos x 0 2 1
De donde:
C 1 ( x)
0
senx
csc x
cos x
cos x
senx
1 1
1. C 1 ( x) 1 C 1 ( x) x.
senx cos x
cos x C 2 ( x)
0
senx csc x cos x
ctgx ctgx. C 2 ( x) ctgx C 2 ( x) ln senx . senx 1
senx cos x
Entonces la Solución general de la ecuación es: Y y g y p Y C 1 cos x C 2 senx x cos x senx ln senx .
2.- y 4 y 4 y x 2 e 2 x Solución: a.- Encontramos la Solución de su respectiva ecuación homogénea. y 4 y 4 y 0
Polinomio característico: p(m) 0 m 2 4m 4 0 m12 2 y g C 1e 2 x C 2 xe 2 x . b.- La Solución particular de la ecuación diferencial viene dada por: y p C 1 ( x)e 2 x C 2 ( x) xe 2 x , Tal que
C 1 ( x)e 2 x C 2 ( x) xe 2 x 0 ; Dividido todo para e 2 x 2 x 2 x 2 2 x 2C 1 ( x)e C 2 ( x)e (1 2 x) x e C 1 ( x) C 2 ( x) x 0 2 2C 1 ( x) C 2 ( x)(1 2 x) x
91
Análisis Matemático
III
De donde:
C 1 ( x)
0
x
x 2
1 2 x
1
x
x 1 1
x 1 . C 1 ( x) x 1 C 1 ( x) ln(x).
2 1 2 x
1
0
2 x 2
C 2 ( x)
1
x
x 2 1
x 2 . C 2 ( x) x 2 C 2 ( x) x 1 .
2 1 2 x
Entonces la Solución general de la ecuación es: Y y g y p Y e 2 x (C 1 C 2 x ln( x) 1) /// .
3.- y y y y 4 x e x Solución: a.- Encontramos la Solución de su respectiva ecuación homogénea. y y y y 0
Polinomio característico: p(m) 0 m 3 m 2 m 1 0 m12 1....m3 1 y g C 1e x C 2 xe x C 3 e x . b.- La Solución particular de la ecuación diferencial viene dada por: y p C 1 ( x)e x C 2 ( x) xe x C 3 ( x)e x , Tal que
C 1 ( x)e x C 2 ( x) xe x C 3 ( x)e x 0 x x x ; C 1 ( x)e C 2 ( x)e (1 x) C 3 ( x)e 0 C ( x)e x C ( x)e x (2 x) C ( x)e x 4 xe x 2 3 1
De donde:
92
Análisis Matemático 0
xe x
e x
0
e x (1 x)
e x
e (2 x)
e
x
C 1( x)
C 1( x)
C 1( x)
III
x
4 xe
x
e x
xe x
e x
e x
e x (1 x)
e x
e x
e x (2 x)
e x
e x
0
e x
e x
0
e x
e x
4 xe x
e x
x
x
e
xe
e
x
x
e
e (1 x )
e x
e x
e x ( 2 x)
e x
e x
xe x
0
e x
e x (1 x)
0
e x
e x ( 2 x)
4 xe x
e x
xe x
e x
e x
e x (1 x )
e x
e x
e x ( 2 x)
e x
x
1
1
3
2
C 1( x) 2 x 2 x C 1 ( x) x 3 x 2 .
C 2 ( x) 2 x C 1 ( x) x 2 .
1
1
3
2
C 1( x) 2 x 2 x C 1 ( x) x3 x 2 .
C 2 ( x) x 2 C 2 ( x) x 1 .
La solución general de la ecuación es: y y g y p y e2 x (C 1 C 2 x ln( x) 1) /// .
4.6 Ecu aciones de Eu ler. La ecuación lineal de la forma:
ax bn y n A1 ax bn 1 y n1 .... An1 ax b y An y f x , (1) Donde a, b, A1 ,...., An1 , An , son constantes, se llama ecuación de Euler. Para el recinto ax + b > 0, introducimos una nueva variable independiente t , poniendo: ax b e t . De donde encontrando las derivadas:
93
Análisis Matemático
III y ae t
dy dt
,
y a e n
2
2t
d 2 y dy 2 , dt dt
d 3 y d 2 y dy y a e 3 3 2 2 y así sucesivamente, dt dt dt Y reemplazando en la ecuación de Euler, ésta se transforma en una ecuación lineal con m
3t
3
coeficientes constantes en “t”. Si el valor de “b” es igual a cero (0) en la ecuación de Euler, se puede buscar directamente soluciones particulares de la forma y x k , obteniendo para k una ecuación que coincide con la ecuación característica de su respectiva ecuación lineal homogénea en “t”.
Ejemplos: Resolver las siguientes ecuaciones diferenciales de Euler:
1.- (2 x 1) 2 y 2(2 x 1) y 4 y 0 (1) Solución: Haciendo en la ecuación (1) la sustitución 2 x 1 e t , se tiene: y 2e
t
dy dt
...;.... y 4e
2t
(
toma la forma: 2 t
e .4e
2t
(
d 2 y dt 2
dy dt
) 2e .2e t
d 2 y
2
dt
t
dy dt
dy dt
). Entonces sustituyendo y simplificando la ecuación
4 y 0
d 2 y dt 2
2
dy dt
y 0 (2).
De donde su ecuación característica y sus raíces son: m 2 2m 1 0 m12 1. y la Solución general de (2) es: y C 1e t C 2 te t pero como 2 x 1 e t t ln(2 x 1) resulta y C 1 (2 x 1) C 2 (2 x 1) ln( 2 x 1) /// .
2.- ( x 2) 2 y 3( x 2) y 3 y 0 (1) Solución: Haciendo en la ecuación (1) la sustitución x 2 e t , se tiene: y e t
dy dt
...;.... y e 2t (
la forma: 2 t
e .e
2t
(
d 2 y dt 2
dy dt
) 3e .e t
d 2 y dt 2
t
dy dt
dy dt
). Entonces sustituyendo y simplificando la ecuación toma
3 y 0
d 2 y dt 2
2
dy dt
3 y 0 (2).
94
Análisis Matemático
III
De donde su ecuación característica y sus raíces son: m 2 2m 3 0 m1 1.... y.....m2 3 y la Solución general de (2) es: y C 1e t C 2 e 3t pero como x 2 e t , resulta y C 1 ( x 2) C 2 ( x 2) 3 /// .
3.- x 2 y 2 x y 6 y 0 (1) Solución: Haciendo en la ecuación (1) la sustitución x e t , se tiene: y e
t
dy dt
...;.... y e
2t
(
la forma: 2 t
e .e
2t
(
d 2 y dt 2
dy dt
) 2e .e t
d 2 y dt 2
t
dy dt
dy dt
). Entonces sustituyendo y simplificando la ecuación toma
6 y 0
d 2 y dt 2
dy dt
6 y 0 (2).
De donde su ecuación característica y sus raíces son: m 2 m 6 0 m1 2.... y.....m2 3 y la Solución general de (2) es: y C 1e 2t C 2 e 3t pero como x e t , resulta y C 1 x 2 C 2 x 3 /// .
Segundo método: Se busca una solución de la ecuación dad de la forma y x k , donde k es un número desconocido. Se tiene: y kx k 1 ...;.... y k (k 1) x k 2 Entonces sustituyendo y simplificando la ecuación (1) toma
la forma: x 2 k (k 1) x k 2 2 xkx k 1 6 x k 0 k (k 1) x k 2kx k 6 x k 0 / x k 0
Tenemos: k (k 1) 2k 6 0 k 2 k 6 0 k 1 2.... y.....k 2 3 . El sistema fundamental de soluciones son: y1 x 2 .... ..... y 2 x 3 y su Solución general será: y C 1 x 2 C 2 x 3 /// .
4.- x 2 y 2 y x 2 y 2 y 0 , todo multiplicando por “x” x 3 y 2 xy 0 (1) Por el segundo método: Se busca una Solución de la ecuación dada de la forma y x k , donde k es un número desconocido. Se tiene: y kx k 1 ...;.... y k (k 1) x k 2 ....;.... y k (k 1)(k 2) x k 3 .
Entonces
sustituyendo
y
simplificando la ecuación (1) toma la forma: x 3 k (k 1)(k 2) x k 3 2 xkx k 1 0 k (k 1)(k 2) x k 2kx k 0 / x k 0
Tenemos: 95
Análisis Matemático
III
k (k 1)(k 2) 2 0
k 3 3k 2 0 k 12 0.... y.....k 3 3 . El sistema fundamental de
soluciones son: y1 1....;.... y 2 ln( x).... ... y3 x 3 y y C 1 C 2 ln( x) C 3 x 3 /// .
su
solución
general
será:
4.7 Ej er cicios propuestos. En los siguientes ejemplos determine mediante el Wronskiano la independencia o dependencia lineal del siguiente grupo de funciones. 1. 4, x 2. 1, 2, x, x2 3. x, 2x, x2 4. sen x, cos x, cos 2x 5. 5, arctg x, arctg x 6. ex, 2ex, e-x 7. π, arcsen x, arccos x 8. 4, sen2x, cos 2x 9. e-3xsen 2x, e-3xcos 2x 10. 1/x, e1/x. Hallar la Solución de las siguientes ecuaciones lineales homogéneas con coeficientes constantes. 11. y IV y 0 12. y IV 2 y 4 y 2 y 5 y 0 13. y V 4 y IV 5 y 5 y 6 y 4 y 0 14. y 3 y 2 y 0 15. y 6 y 11 y 6 y 0 16. y 3 y 3 y y 0;... y(0) 1,.... y (0) 2,... y (0) 3. 17. y 4 y 3 y 0;... y(0) 6,... y (0) 10 18. y 2 y 2 y 0 19. y 2 y y 0 20. 3 y 2 y 8 y 0 21. y 2 y 3 y 0;... y(0) 1,... y (0) 3 22. y 4 y 3 y 0;... y(0) 6,... y (0) 10 23. y 8 y 0 24. y 2 y 0 25. y 5 y 4 y 0;... y(0) 5,... y (0) 8 Indicar la forma de las soluciones particulares de las siguientes ecuaciones lineales no homogéneas. 26. y 4 y x 2 e 2 x
96
Análisis Matemático
III
27. y 4 y 4 y sen2 x e 2 x 28. y 4 y 3 y xe x cos 2 x x 2 e x sen2 x. 29. y 4 y senxsen2 x 30. y y x 31. y 4 y 3 y 1 32. y y 3 Hallar la Solución de las siguientes ecuaciones lineales no homogéneas con coeficientes constantes. 33. y 2 y y 2 34. y 4 y 4 y 2 x 2 35. y 2 y 2 y e x cos x x e x 36. y y senx 37. y y x 2 senx. 38. y 2 y y 3 sen2 x cos x 39. y y y 1 senx x x 2 40. y 6 y 9 y 9 xe 3 x 41. y y y y xe x 42. y IV 2 y y x 3 Hallar la Solución de las siguientes ecuaciones lineales no homogéneas por el método de variación de las constantes arbitrarias. 43. y y tgx 44. y 4 y 45. y y
1 cos 2 x 1
e x 1 e x 46. y 2 y y x e x 47. y 2 y y 2 x 1 2 x 48. y y e cos e x
49. Hallar la Solución de la ecuación y 4 y senx , que satisface las condiciones y = 1, y`= 1 para x = 0. 50. Hallar la Solución de la ecuación y 2 y e 2 x x 2 1 , que satisface las condiciones y = 1/8, y`= 1 para x = 0. Hallar la Solución de las siguientes ecuaciones de Euler. 51. x 2 y x y 3 y 0 97
Análisis Matemático
III
52. x 3 y 3 x 2 y 6 x y 6 y 0 53. x 2 y x y 4 y 0 54. (1 x) 2 y 3(1 x) y 4 y 0 55. (1 x) 2 y 3(1 x) y 4 y (1 x) 3 56. Hallar la Solución particular de la ecuación x 2 y x y y x , cuyas condiciones de frontera son: y = 0, y`= 1 para x = 1. 57. x 2 y x y y x(6 ln x).
4.8 Apl icaciones de las ecuaci ones 4.8.1 Aplicaciones físicas.
orden.
4.8.1.1 Movimiento armónico simple. Supongamos que un cuerpo de masa “m” está sujeto a un resorte flexible suspendido de un
soporte rígido, según la Ley de Hooke que dice que el resorte ejerce una fuerza de restitución opuesta a la dirección del alargamiento y proporcional a su magnitud. Y de acuerdo con la Ley de Newton. La ecuación diferencial que rige este movimiento es: m d 2 x 2
dt
k m
x 0:
d 2 x dt 2
kx. o
donde k es la constante de proporcionalidad.
Remplazando 2
k m
, la ecuación toma la forma:
d 2 x 2
dt
2 x 0 (1).
Introduciendo operadores tendríamos ( D 2 2 ) x 0 D
d
.: dt
Para resolver este problema tenemos que “x” es el desplazamiento de la masa desde la posición de equilibrio. Las condiciones iniciales obvias son dos: x= α y dx/dt = β, para t = 0,
también adoptándose la convención de que los desplazamientos medidos hacia abajo desde la posición de equilibrio se consideran positivos y los que van hacia arriba negativos, tendremos que su Solución general tendrá la forma: x(t ) C 1 cos( wt ) C 2 sen(wt ). (2).
Donde las constantes se determinan con las condiciones iniciales. Luego podremos determinar el periodo de las vibraciones y la frecuencia de las mismas utilizando las siguientes fórmulas: T
2
, y
f
1
T
2
.
Una forma práctica de poner la Solución de la ecuación anterior, es de la forma: x(t ) Asen( .t ) (3). 98
Análisis Matemático
III
Donde A es la amplitud de las oscilaciones y se relaciona con las constantes arbitrarias por la formula: A C 1 2 C 2 2 y es un ángulo de fase y se define por: sen
C 1 A
; cos
C 2 A
.
Esta forma de expresar la Solución es muy conveniente, ya que nos permite calcular los valores del tiempo en los cuales la gráfica de x (t) atraviesa el eje t positivo (la recta x = 0). Observando también que sen( .t ) 0 cuando .t n donde n es un entero positivo. En la ecuación (1) se supone que en el sistema no hay rozamiento. Considérese ahora el caso en que el movimiento del peso sea retardado por una fuerza de rozamiento β (dx/dt) proporcional a la velocidad, donde β es una constante positiva. Entonces, la ecuación diferencial es: m
d 2 x dt 2
kx
O bien
d 2 x 2
dt
2b
dx dt
dx dt
,
2 x 0, (4)
2b
m
y
2
k m
.
Si introducimos el operador D = d/dt, la ecuación (2) se convierte en: ( D 2 2bD 2 ) x 0.
De aquí tenemos que la ecuación característica es: m 2 2bm 2 0, con raíces m b b 2 2 . Surgen así tres casos, dependiendo del las magnitudes relativas de b y ω.
Caso 1.- Si b = ω, las raíces de la ecuación son iguales, y la Solución es: x(t ) (C 1 C 2 t )e
t
. (5)
Observamos que al transcurrir el tiempo, x tiende a cero. El movimiento no es oscilatorio. Caso 2.- Si b > ω, las raíces de la ecuación son reales pero desiguales, y la Solución es: x(t ) C 1e m1t C 2 e m2t . (6)
También podemos observar que al transcurrir el tiempo, x va ha tender cero. El movimiento tampoco es oscilatorio.
99
Análisis Matemático
III
Caso 3.- Si b < ω, sea ω2 - b2 = a2. Entonces las raíces serán imaginarias del tipo: m1 b ai,
y
m2 b ai, y la Solución es:
x(t ) e bt (C 1 cos at C 2 senat ). (7).
Esta ecuación realizando algunas acotaciones podemos también poner como: x(t ) Ae bt sen(at ) . (8)
Esta ecuación representa un movimiento vibratorio amortiguado. Es análogo al movimiento armónico simple, de periodo T 2 a , excepto que en su amplitud no es constante, si no que esta dado por Ae-bt. Como esta última expresión tiende a cero
Mientras t crece, las vibraciones tienden a desaparecer con el transcurso del tiempo. Observamos sin embargo que la ecuación (8) se reduce a la (3) en ausencia de rozamiento. El efecto del rozamiento es doble: 1. b = c / (2m) aparece como coeficiente en el factor de amortiguamiento exponencial e-bt. Cuanto mayor sea b, con más rapidez tienden las vibraciones a hacerse imperceptibles. 2. El periodo T 2 / a 2 / 2 b 2 es mayor que el T0 = 2π / ω en el sistema sin rozamiento. En las figuras que se encuentran a continuación se representan las soluciones de la ecuación. La magnitud de b con respecto a ω, determina:
Ej emplo 1. Un cuerpo que pesa 3lb estira un muelle 6plg (1/2 pie). Dicho cuerpo se suelta en
t = 0 desde un punto que esta a 8plg (2/3 pie) bajo la posición de equilibrio, con una velocidad dirigida hacia arriba de 4/3 pie/sg. Determine la función x (t) que describe el movimiento libre resultante.
Solución: Por la ley de Hooke tenemos: 1 lb . F kx , reemplazando valores 3 k ( ) k 6 2 pie
Además pasamos las unidades de peso a unidades de masa: m
W g
3 32
slug .
Estos valores reemplazando en la ecuación que rige el movimiento tenemos: d 2 x
6 d 2 x 0 2 64 x 0 . Esta es una ecuación lineal homogénea 3 dt 2 dt 32
100
Análisis Matemático
III
Resolviendo: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 2 64 0 ; resolviendo encontramos sus raíces m 8.i ; por lo que su sistema fundamental de soluciones es: cos .8t ; sen.8t y la Solución general es: x(t ) c1 cos .8t c2 sen.8t /// . Derivando esta Solución: x(t ) 8c1 sen.8t 8c2 cos .8t /// . Y utilizando las condiciones iniciales: x (0) = 2/3, dx/dt (0) = -4/3 (El signo negativo por ser opuesta al movimiento (hacia arriba)), en estas ecuaciones encontramos los valores de c1 y c2. c1
2 3
,
1
y c2 . 6
2
1
3
6
Reemplazando en la Solución general tenemos: x(t ) cos .8t sen.8t /// . Ej emplo 2. Un cuerpo que pesa 64lb sujeto al extremo de un muelle lo estira 0.32 pie. El
cuerpo ocupa una posición que esta a 8plg (2/3 pie) sobre la posición de equilibrio, y desde ahí se le comunica una velocidad dirigida hacia abajo de 5 pie/sg. a) Encuentre la ecuación del movimiento. b) Calcule la amplitud y el periodo del movimiento. c) Cuál es la posición y velocidad instantánea en 0,3 segundos.
Solución: a) Por la ley de Hooke tenemos: F kx , reemplazando valores 64 k (0.32) k 200
lb pie
.
Además pasamos las unidades de peso a unidades de masa: m
W g
64 32
2. slug .
Estos valores reemplazando en la ecuación que rige el movimiento tenemos: d 2 x dt 2
100 x 0 . Esta es una ecuación lineal homogénea
Resolviendo: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 2 100 0 ; resolviendo encontramos sus raíces m 10i ; por lo que su Solución general es: x(t ) c1 cos .10t c2 sen.10t /// . Derivando esta Solución: x(t ) 10c1 sen.10t 10c2 cos.10t /// . Y utilizando las condiciones iniciales: x (0) = 2/3, encontramos los valores de c1 y c2. c1
2 3
,
dx/dt (0) = 5; en estas ecuaciones
1
y c2 . 2
101
Análisis Matemático
III 2
1
3
2
Reemplazando en la Solución general tenemos: x(t ) cos .10t sen.10t /// . Que es la ecuación del movimiento.
2 2 b) La amplitud es igual: A c12 c2 2 A 2 3 1 2 0.8333.
El periodo es: T
2
, T
2 10
0.628
c) La posición después de 0,3 segundos seria: 1 cos .10.(0.3) sen.10(0.3) x(0.3) 0,4645 pie /// . 3 2 20 Y la velocidad: x(t ) sen.10(0.3) 5 cos .10(0.3) 4.02 pie seg /// . 3 Ej emplo 3. Un cuerpo que pesa 16lb se sujeta al extremo de un muelle de 5 pie de largo. En estado de equilibrio, el resorte mide 8.2 pie. Si el peso se empuja hacia arriba y se suelta, a partir del reposo, desde un punto que esta 2 pie sobre la posición de equilibrio. Determinar los desplazamientos x(t) sabiendo además que el medio ofrece una resistencia numéricamente x(0,3)
2
igual a la velocidad instantánea.
Solución: El alargamiento experimentado por el resorte después de que se le sujeta el peso es: 8.2 - 5 = 3.2 pie, luego por la ley de Hooke se tiene: F kx , reemplazando valores 16 k (3.2) k 5
lb pie
.
Además pasamos las unidades de peso a unidades de masa: m
W g
16 32
1
. slug . 2
Estos valores reemplazando en la ecuación que rige el movimiento tenemos: 1 d 2 x 2 dt 2
5 x
dx dt
0
d 2 x dt 2
2
dx dt
10 x 0 . Esta es una ecuación lineal homogénea
Resolviendo: El polinomio característico, correspondiente a la ecuación diferencial dada es: P (m) m 2 2m 10 0 ; resolviendo encontramos sus raíces m 1 3i ; por lo que su Solución general es: x(t ) e t (c1 cos .3t c2 sen.3t /// . Derivando esta Solución: x(t ) e t (3c1 sen.3t 3c2 cos .3t ) e t (c1 cos .3t c2 sen.3t ) /// . Y utilizando las condiciones iniciales: x (0) = -2, encontramos los valores de c1 y c2. c1 2 , y c 2
2 3
dx/dt (0) = 0; en estas ecuaciones
.
102
Análisis Matemático
III 2
Reemplazando en la Solución general tenemos: x(t ) e t (2 cos .3t sen.3t ) Que es la 3
ecuación del movimiento.
En conclusión podemos decir que como las raíces del polinomio característico son imaginarias el sistema es sub amortiguado.
4.8.1.2 El cable colgante. Considere un cable o una cuerda que cuelga de dos apoyos A y B (figura 4.1a), no necesariamente al mismo nivel. Asuma que el cable es flexible de modo que debido a su carga (la cual puede ser debida al propio peso o fuerzas externas actuantes, o a una combinación de estas) toma la forma como la de la figura. Sea C la posición más baja del cable, y escoja los ejes x e y, donde el eje y pasa por C. Considere aquella parte del cable entre el punto más bajo y cualquier punto P en el cable con coordenadas (x,y). Esta parte estará en equilibrio debido a la tensión T en P (figura 4.1b), la fuerza horizontal H en C, y la carga vertical total en el segmento CP del cable denotada W (x) la cual se asume que actúa en algún punto Q, no necesariamente en el centro del arco CP.
Figura 4.1a
Figura 4.1b
Descomponiendo la tensión T en sus dos componentes T x = T .cosӨ y T y = T.senӨ y aplicando las ecuaciones de equilibrio en la figura b,
F 0 Tenemos:
F
y
x
y
F 0 x
F
y
0
0 . H - T.cosӨ = 0 T cos H (1) T.senӨ - W = 0 Tsen W (2)
Dividiendo (2) para (1), y usando la definición de que tg dy dx
W H
dy dx
, se tiene:
. (3)
103
Análisis Matemático
III
En esta ecuación H es constante, puesto que es la tensión en el punto más bajo, pero W puede depender de x. Derivando la ecuación (3) respecto a x tendremos: d 2 y dx
2
Ahora
1 dW H dx dW dx
.
representa el incremento en W por unidad de incremento en x, esto es, la carga
por unidad de distancia en la dirección horizontal. Ejemplo 1: Determinar la forma que adopta un cable uniforme, suspendido de sus dos
extremos, y sometido únicamente a la acción de su propio peso w kg/m de longitud.
Figura 4.2 Como ya dijimos H es constante, debiéndose aquí a la parte OQ del cable, mientras W ws, donde s es la longitud de OP. Luego d 2 y dx
2
1 dW H dx
w ds H dx
Las condiciones son:
w H
1 (dy / dx) 2 . (a)
y = 0 cuando x = 0 dy/dx = 0 cuando x = 0.
Como la ecuación (a) no es lineal, tenemos que hacer el cambio de variable que la ecuación quedaría: dp dx
w H
1 p 2
dp 1 p 2
w H
dy dx
dx. Integrando entre los límites x = 0; p = 0 y
p , por lo
x = x;
p = p tenemos:
104
Análisis Matemático p
0
dp 1 p 2
Como p
w
x
H 0
dy dx
III
dx arcsen.h. p
dy dx
w
w x p sen.h.( x). H H
w x). dy sen.h( x)dx , integrando entre los límites x = 0, H H
sen.h.(
w
y = 0 y x = x, y = y, tenemos: y
x
w H w Cos.h( x) 1 ///. Esta curva es una catenaria. Si se dy sen.h( x)dx y H w H o 0
hubiera tomado el origen de coordenadas a una distancia H/w por debajo del punto más del cable, la ecuación de la curva tomaría la forma: y
H
w Cos.h( x). w H
4.8.1.3 Defl exi ón de vigas.
El problema consiste en determinar la curva que forma la elástica de la viga al ser sometida a flexión por acción de su peso propio y de cargas externas o la combinación de ambas, estas cargas pueden ser distribuidas, puntuales etc. Aquí se estudia únicamente las vigas que son uniformes, tanto en forma como en material, y es conveniente considerarlas como si estuviesen constituidas por fibras dispuestas longitudinalmente.
Figura 4.3 En la figura 3 podemos observar, que las fibras de la parte superior están comprimidas y las de la mitad inferior alargadas; a ambas mitades los separa una superficie cuyas fibras ni están comprimidas ni alargadas. La fibra, que en un principio coincidía con el eje horizontal de la viga, está ahora en esta superficie de separación dispuesta según una curva (la curva elástica o de flexión). A continuación vamos a tratar de obtener la ecuación de esta curva. Consideremos una sección transversal de la viga a una distancia x de un extremo. Sea CD su intersección con la superficie de separación y P su intersección con la curva elástica. Se demuestra en mecánica que el momento M con respecto CD de todas las fuerzas exteriores que actúan sobre cualquiera de los dos segmentos en que la sección transversal divide a la viga. 1.-
105
Análisis Matemático
III
Es independiente del segmento considerado y 2.- La expresión que relaciona esta condición es: (A)
EI R
M .
Donde E es el modulo de elasticidad de la viga, I el momento de inercia de la sección transversal con respecto a CD con constantes asociadas a la viga, R es el radio de curvatura de la curva elástica en P y M es el momento flector en la sección transversal (punto P de la curva elástica) y es igual al suma algebraica de los momentos de las fuerzas exteriores que actúan sobre el segmento de viga. Supondremos que fuerzas hacia arriba producen momentos negativos y las fuerzas hacia abajo dan momentos positivos Al estar definido R por: dy 2 1 ( dx ) R 2 d y
3/ 2
, y sabiendo que en la práctica la pendiente de
dx 2
la curva elástica en todos los puntos es numéricamente pequeña, podemos tener una buena aproximación tomando a R
1 d 2 y
, por lo que reemplazando en la formula (A) tendremos:
dx 2
(B) EI
d 2 y dx 2
M /// .
Ej emplo 1. Una viga horizontal de 2l metros de longitud esta simplemente apoyada en sus
extremos (figura 4.4). Hallar la ecuación de su curva elástica y su máxima deformación vertical (flecha máxima) cuando tiene una carga uniformemente distribuida de w Kg/m.
Figura 4.4 Solución: Puesto que la viga es simplemente apoyada en O y R, cada uno de estos soportes lleva la mitad del peso de la viga, es decir wl. Hallando el momento flector M(x) en el punto
P, escogiendo el lado izquierdo de P, podemos observar que actúan dos fuerzas: 1. La fuerza hacia arriba wl , a una distancia x de P , produciendo un momento negativo. 2. La fuerza hacia abajo wx, a una distancia x/2 (centro de gravedad de OP ) de P , produciendo un momento positivo. Entonces el momento resultante total en P seria: 106
Análisis Matemático
III
wl x wx 2 wlx M ( x) x wx M ( x) /// . 2 2 2 2
Con este valor la ecuación (B) quedaría: EI
d 2 y dx
2
wx 2 2
wlx 2
. (1) dy
wlx 2
C 1 . (2) 6 2 Con la condición de que dy/dx = 0 para x = l en el punto medio de la viga encontramos 1 C 1 wl 3 . Por lo que la ecuación (2) quedaría: 3 dy wx 3 wlx 2 1 3 wl . (3). EI dx 6 2 3 wx 4 wlx 3 1 3 wl x C 2 . (4) Integrando (3) tenemos: EIy 24 6 3 Con la condición de que y = 0 para x = 0, encontramos C 2 0 . Por lo que la ecuación
Integrando (1) una vez tenemos: EI
dx
wx 3
(4) quedaría: EIy
wx 4 24
wlx 3 6
1
w
3
24 EI
wl 3 x y
( x 4 4lx 3 8l 3 x) /// . Que es la ecuación pedida.
Para encontrar la flecha máxima debemos tomar a x = l , ya que esta se produce en el punto medio de la viga, por lo que tendríamos y máx
w 24 EI
(l 4l 8l ) 4
4
4
5wl 4 24 EI
/// . :
Ej emplo 2. Una viga horizontal de l metros de longitud esta empotrada en un extremo y libre
en otro (figura 4.5). Hallar la ecuación de su curva elástica y su máxima deformación vertical (flecha máxima) cuando tiene una carga uniformemente distribuida de w Kg/m.
Figura 4.5 107
Análisis Matemático
III
Solución: Puesto que la viga es en voladizo tomamos el origen en el punto fijo O y sean (x,y) las coordenadas del punto P, Considerando el segmento PR, vemos que la única fuerza que tenemos es el peso w(l – x)kg en el punto medio de PR, es decir a ½(l-x) metros de P, tomando el momento flector M(x) en el punto P tenemos: 1 1 M ( x) w(l x). (l x) w(l x) 2 /// . 2 2
Con este valor la ecuación (B) quedaría: EI
d 2 y dx
2
1
w(l x) 2 . (1) 2
Integrando (1) una vez tenemos: EI
dy dx
1
w(l x) 3 C 1 . (2) 6
1
Con la condición en O: dy/dx = 0 para x =0 encontramos C 1 wl 3 . Por lo que la 6
ecuación (2) quedaría: EI
dy dx
1
1
6
6
w(l x) 3 wl 3 . (3).
Integrando (3) tenemos: EIy
1 24
w(l x) 4
1 6
wl 3 x C 2 . (4)
Con la condición en O: y = 0 para x = 0, encontramos C 2
1 24
wl 4 . . Por lo que la
ecuación (4) quedaría: EIy
1 24
w(l x) 4
ecuación pedida.
1 6
wl 3 x
1 24
wl 4 EIy
w 24 EI
(4lx 3 6l 2 x 2 x 4 ) /// . Que
es la
Para encontrar la flecha máxima debemos tomar a x = l , ya que esta se produce en el extremo de la viga, por lo que tendríamos y máx
w 24 EI
(4l 6l l ) 4
4
4
wl 4 8 EI
/// .
4.8.1.4 Ci r cuitos Elé ctr icos.
Es la aplicación de la segunda ley de Kirchov que dice: La suma de las caídas de tensión a través de los elementos de un circuito cerrado es igual a la fuerza electromotriz total E en el circuito. Así tenemos que la caída de tensión a través de una resistencia de R ohmios es Ri, a través de una bobina de L henrios de inductancia es L di/dt y a través de un condensador de C faradios de capacidad es q/C. También tenemos que la corriente i amperios y la carga q culombios están relacionadas por i = dq/dt. R, L y C son constante
108
Análisis Matemático
III
Figura 4.6 La ecuación diferencial de un circuito eléctrico que contiene una inductancia L, una resistencia R, un condensador C y una fuerza electromotriz E (t) (figura 4.6) es: L
di dt
L
Ri
d 2 q dt 2
R
q C
E (t ). (1) o bien, reemplazando i = dq/dt, di/dt = d 2q/dt 2, tenemos:
dq dt
q C
E (t ). (2). De donde podemos encontrar q = q(t).
Derivando (2) y empleando dq/dt = i tendremos: L
d 2 i dt 2
R
di dt
i C
E (t ). De donde determinamos i = i (t).
Ej emplo 1. Un circuito eléctrico consta de una inductancia de 0.2 henrios, una resistencia de -6 25 ohmios y un condensados cuya capacidad es de 30 microfaradios (un microfaradio = 10 faradios), (figura 4.7). Hallar la carga q y la corriente i en el tiempo t, siendo las condiciones iniciales q = 0,05 culombios, i = dq/dt = 0 para t = 0.
Figura 4.7 Solución: Como L = 0.2, R = 25, C = 30 x 10-6 y E (t) = 0, tenemos: 0.2
d 2 q dt 2
25
dq dt
q 25 x10 6
0
d 2 q dt 2
125
dq dt
200000q 0 . Esta es una ecuación lineal
homogénea Resolviendo: El polinomio característico, correspondiente a la ecuación diferencial dada es:
109
Análisis Matemático
III
P (m) m 2 125m 200000 0 ; resolviendo encontramos sus raíces m 62,5 442.8i ; por
lo que la Solución general es: q(t ) e 62.5t (c1 cos 442.8t c2 sen442.8t ) //
Derivando esta Solución: q(t ) i e 62.5t (442.8c2 62.5c1 ) cos 442.8t (62.5c2 442.8c1 ) sen442.8t ) //
Utilizando las condiciones iniciales: q (0) =0.05, dq/dt (0) = i (0) = 0 en estas ecuaciones encontramos los valores de c1 y c2. c1 0.05 , y c 2 0.141 . Por lo que la carga q y la intensidad i son:
q(t ) e 62.5t (0.05 cos 442.8t 0.141 sen442.8t ) // q(t ) i e 62.5t 30.9 sen442.8t ) q(t ) i 30.9e 62.5t sen442.8t ///
4.8.2 Ej er cicios propuestos.
1. Una masa de ½ slug está suspendida de un resorte cuya constante es 4 lb/pie. El sistema se pone en movimiento en un medio que opone una fuerza de amortiguación numéricamente igual a dos veces la velocidad instantánea. Encuentre la Solución estacionaria si se aplica al sistema, a partir de t = 0, una fuerza exterior f (t) = 20sen 2t.
2. Un péndulo de ½ pie (15.24 cm) de longitud se suelta con una velocidad de ½ radián/seg. Hacia la vertical, desde una posición determinada por la distancia de 1/5 rad respecto de la vertical. Hallar la ecuación del movimiento.
3. Una partícula de masa m es repelida de O con una fuerza igual a k>0 veces la distancia desde O. Si la partícula parte del reposo a una distancia a de O, hallar su posición t segundos más tarde.
4. Si en el problema 3 k=m y a=12 pies=365.76 cm, determinar a) la distancia desde O y la velocidad cuando t=2 seg, b) cuándo estará a 18 pies=548.64 cm de O y la velocidad que tendrá entonces.
5. Una cadena colocada sobre una clavija pulida pende 8 dm de un lado y 10 dm del otro. Si la fuerza de rozamiento es igual al peso de 1 dm de cadena, hallar el tiempo que tarda la cadena, al resbalar, en caerse.
110
Análisis Matemático
III
6. Si el interior de dos superficies esféricas concéntricas de radios r 1 y r 2 , r 1 < r 2 , tiene una carga eléctrica, la ecuación diferencial para el potencial V en un punto cualquiera entre las dos superficies esféricas y a una distancia r de su centro común es d 2V dr 2
2 dV r dr
0
Resolver para V , siendo V=V 1 cuando r=r 1 y V=V 2 cuando r=r 2.
7. Se tiene un resorte que se alarga 3 pulgadas (7,62 cm) si se le aplica un peso de 9 libras (4,082 kg.). Se suspende de él un peso de 24 libras (10,886 kg.) y se interrumpe el estado de reposo en que se encontraba. Hallar la ecuación del movimiento si una vez suspendido dicho peso a) Se tira hacia abajo 4 pulgadas (10,16 cm) y se suelta después. b) Se tira hacia abajo 2 pulgadas (5,08 cm) y se da una velocidad hacia arriba de 24 pulgadas/seg (60,96 cm/seg). c) Se tira hacia abajo 3 pulgadas (7,62 cm) y se da una velocidad hacia debajo de 48 pulgadas/seg (121,92 cm/seg). d) Se empuja hacia arriba 3 pulgadas (7,62 cm) y se suelta después. e) Se empuja hacia arriba 4 pulgadas (10,16 cm) y se da una velocidad hacia arriba de 60 pulgadas/seg (152,40 cm/seg).
8. Las características de un resorte son tales que se alarga 3 pulgadas (7,62 cm) bajo la acción de un peso de 30 libras (13,608 kg). Se suspende de él un peso de 64 libras (29,030 kg) y se interrumpe su estado de reposo. La resistencia del medio es, numéricamente, igual a 8 dx/dt libras. Hallar la ecuación del movimiento del peso si a) Comienza a moverse hacia abajo con una velocidad de 10 pies/seg (3,048 m/seg). b) Se tira hacia abajo 6 pulgadas (15,24 cm) y se da una velocidad hacia arriba de 10 pies/seg.
9. Se tiene un resorte que se alarga 6 pulgadas (15,24 cm) bajo la acción de un peso de 3 libras (1,361 kg). Se suspende de él un peso de 3 libras y se interrumpe su estado de reposo. Se tira, entonces, hacia abajo del peso hasta desplazarle 3 pulgadas y se suelta. Determinar la ecuación del movimiento si a) Actúa en el resorte una fuerza fija 3/2 sen 6t .
111
Análisis Matemático
III
b) Actúa en le resorte una fuerza fija 3/2 sen 8t.
10. Sea una ménsula de longitud l m (viga empotrada en un extremo y libre en el otro). Hallar la ecuación de la curva elástica y la flecha máxima si hay una carga uniforme de w kg/m y una carga de W kg en el extremo libre.
11. Una viga de 2l m de longitud está apoyada en ambos extremos y tiene una carga uniforme de w kg/m. Tómese el origen en el punto medio (punto más bajo) de la viga y hállese la ecuación de la curva elástica y de la flecha máxima. I y wl l x 1 / 2wl x 2 1 / 2wl 2 x 2
Sugerencia: y y 0
y
x0
12. Una viga de 3l m de longitud está apoyada en ambos extremos. Hay una carga uniforme de w kg/m y cargas concentradas de W kg a una distancia de l m de cada extremo. Tomando el origen como en el problema 11, hállese la flecha máxima. 3l 3l 1 x 2 W x , x ; y 2 4 2 2 2 Sugerencia w 9l 2 1 x 2 Wl , 0 x . M 2 4 2 M
w 9l 2
13. Un circuito consta de una inductancia de 0.05 henrios, una resistencia de 5 ohmios y un condensador de 4104 faradios de capacidad. Si q=i=0 para t=0, hallar q e i en función de t cuando a) hay una f.e.m. constante = 110 voltios,
b) hay un f.e.m.
alterna = 200 cos 100t . Hallar las soluciones de régimen permanente en b).
14. Resolver el problema 13 después de sustituir la resistencia de 6 ohmios por una resistencia de 60 ohmios.
112
Análisis Matemático
III
UNIDAD V 5. Sistema de ecuaciones diferenciales con coeficientes constantes Se llama sistema lineal de n ecuaciones diferenciales con coeficientes constantes a la expresión de la forma: dx dt f 1 (t , x1 x2 ......... xn ) dx2 f (t , x x ......... x ) 2 1 2 n dt ................................... dxn dt f n (t , x1 x2 ......... xn ) 1
Esto también se puede escribir en la forma: dxi dt
aij xi f i (t )
(i = 1, 2, 3, 4, 5,……n),
Donde aij son constantes y f i(t), funciones dadas. Si f i (t ) 0, el sistema se llama homogéneo y si f i (t ) 0, el sistema se llama no Homogéneo. Notación vectori al: X x1 , x 2 ....... xn V n,
f f 1 , f 2 ........ f n ,
dx dt
f t , x
x1 1 (t ),....... x n n (t ) determinadas y diferenciables, con El conjunto de funciones derivadas continuas en el intervalo (a, b), se llaman Solución del sistema en este intervalo, si las funciones convierten en identidades a todas las ecuaciones del sistema.
El problema de la búsqueda de la Solución: x1 1 (t ),....... xn n (t ) que satisfaga a las condiciones iniciales x1 (t 0 ) x10 ,....... xn (t 0 ) x10 , se denominada problema de valor inicial (Cauchy)
113
Análisis Matemático
III
Existen algunos métodos para encontrar su solución, a continuación explicaremos dos procedimientos muy usados. 5.1 Reducci ón de un sistema, a una ecuación de n-é simo orden. Este es el método más sencillo para hallar la Solución del sistema (1), ya que simplificando se obtiene una ecuación diferencial lineal con coeficientes constantes. Ilustraremos el método para un sistema de dos ecuaciones. Consideremos el siguiente sistema: dx dt ax by f (t )..............(1.1) dx cx dy g (t )..............(1.2) dt
Donde a, b, c, d, son constantes, f (t), g (t) son funciones conocidas: x (t), y (t) son funciones 1 dx
ax f (t ) (2) b dt
incógnitas de la ecuación (1) despejamos y Reemplazamos “y” en (1.2) se obtiene:
d 1 dx
d dx ax f (t ) cx ax f (t ) g (t ) dt b dt b dt
2
1 d x b dt 2
adx dt
d dt
d dx f (t ) cx ax f (t ) g (t ) 0
Simplificando se tiene: A
b dt
d 2 x dt 2
B
dx dt
Cx Rt
Donde A, B, C, son constantes. De esta ecuación diferencial de coeficientes constantes, hallamos x = x (t, C 1 , C 2 ) y poniendo en (2) el valor hallado de x, así como dx/dt , hallamos y. Ej emplo 1. Integrar el siguiente sistema de ecuaciones:
dx dt x 2 y......(1.1) dy x 3 y.......(1.2) dt
114
Análisis Matemático
III
Solución: 1 dx De (1.1) se tiene y x , (2) reemplazando en (1.2) 2
dt
d 1 dx 3 dx x x x , hallando la derivada dt 2 dt 2 dt 1 dx 2 dt d 2 x dt 2
1 d 2 x 2 2dt 2
4
dx dt
3
3 dx
2
2 dt
x x
, simplificando se tiene:
5 x 0 , que es una ecuación Diferencial lineal de segunda orden homogénea de
coeficientes constantes. El polinomio característico es: P m m 2 4m 5 0 de donde m1 2 i m2 2 i
La Solución general es: x g (t ) e at (c1 cos bt c2 senbt ) x(t ) e 2t (c1 cos t c2 sen.t ) (3).
Es decir: Derivando (3) tenemos:
x(t ) 2e 2t (c1 cos t c 2 sen.t ) e 2t (c1 sen.t c2 cos t )
x(t ) e 2t (2c1 c2 ) cos t (2c2 c1 ). sen.t )
(3.1)
Reemplazando (3) y (3.1) en (2) y simplificando tenemos: y(t )
1 2
e 2t (c 2 c1 ) cos t (c1 c2 ) sen.t (4).
Por lo que la Solución del sistema es: x(t ) e 2t (c1 cos t c2 sen.t ) 1 2t y ( t ) e (c 2 c1 ) cos t (c1 c 2 ) sen.t 2 Ej emplo 2. Integrar el siguiente sistema de ecuaciones:
dx dt y 1......(1.1) dy x 2t .......(1.2) dt
115
Análisis Matemático
III
Solución: dx
De (1.1) se tiene y
dt
1 , (2) reemplazando en (1.2)
d dx ( 1) x 2t , hallando la derivada dt dt d 2 x dt 2
x 2t
d 2 x dt 2
x 2t (3), que es una ecuación Diferencial lineal de segunda orden
no homogénea de coeficientes constantes, cuya Solución general es: x(t ) x g x p .
Para hallar x g. El polinomio característico es: P m m 2 1 0 de donde m1 1 m2 1
La Solución: x g (t ) c1e t c2 e t (4) Y la Solución particulares: x p = At + B (5) Derivando tenemos: x p A
y
x p 0 (5.1)
Reemplazando (5) y (5.1) en (3) y simplificando tenemos: 0 ( At B) 2t A 2 Y B = 0, estos valores reemplazamos en (5), esa es la Solución
particular buscada. x p = 2t (6).
Por lo que la Solución general de la ecuación es: x (t ) c1e t c2 e t 2t (7). Derivando (7) tenemos: x(t ) c1e t c2 e t 2 (7.1) Reemplazando (7) y (7.1) en (2) y simplificando tenemos: y (t ) c1e t c2 e t 1
(8).
Por lo que la Solución del sistema es: 116
Análisis Matemático
III x(t ) c1e t c 2 e t 2t y (t ) c1e t c 2 e t 1
5.2 M é todo matr icial o de Euler para sistemas lineales con coefi cientes constantes.. Este procedimiento sirve para encontrar la Solución de sistemas de primer orden con coeficientes constantes homogéneos. Consideremos el siguiente sistema: dx1 dt a11 x1 a12 x 2 ... a1n x n dx 2 a x a x ... a x 21 1 22 2 2n n dt ...................... ......................... dx n a n1 x1 a n 2 x 2 ... a nn x n dt
(1).
Las soluciones de este sistema de ecuaciones es de la forma: 1. x1 = 1 en, x2 = 2 en,...., xn = n en , Donde x 1 1e rt rt x 2 2 e
Como:
dx 1 dt dx 2 dt
1 re rt 2 re rt
…........ .............. rt x n n e
dx n dt
n re rt
i , i = 1, ...., n son constantes.
(2).
Reemplazando (2) en (1), simplificando por ert y ordenando se obtiene un sistema de ecuaciones para determinar β1, β2,………βn (a11 r ) 11 a12 2 ... a1n n 0 a 21 1 (a 22 r ) 2 .... a 2 n n 0 ........(3) ...................... ......................... a n1 1 a n 2 2 ... (a nn r ) n 0 117
Análisis Matemático
III
El sistema (3) tiene una Solución no nula si y sólo si su determinante Δ es igual a cero.
a11 r ∆=
a 21
a12
...........
a1n
a 22 r ...........
a2n
......
.......
a n1
an2
.......
.....
=0
(4).
........... a nn r
De este determinante encontramos la ecuación característica P(r) = 0 valores propios del sistema, y cada valor de r determinará β diferente. En caso de tener tres ecuaciones.
dx dt ax by cz dy a1 x b1 y c1 z .......(1) dt dz dt a 2 x b2 y c2 z
Las soluciones son: rt
rt
rt
x = e , y = ue , z = ve , Donde , u, v y r son constantes.
x e rt y e rt z e rt
dx dt dx dt dx dt
.rert .rert .rert
…. (2)
Reemplazando (2) en (1), simplificando por ert y (a r ) b c 0 ordenando se obtiene un sistema de ecuaciones para a1 (b1 r ) c1 0 ........(3) determinar λ, μ, y ν. a b (c r ) 0 2 2 1 El sistema (3) tiene una Solución no nula si y sólo si su determinante Δ es igual a cero.
118
Análisis Matemático
∆=
III
a r
b
c
a1
b1 r
c1
a2
b2
c2 r
=0
(4).
Si r 1, r 2 y r 3 son raíces del polinomio característico P(r) = 0 y son reales y distintas tenemos: Si r = r 1, se resuelve el sistema (3) y se obtiene:1, μ1, ν1 Si r = r 2, se resuelve el sistema (3) y se obtiene:2, μ2, ν2 Si r = r 3, se resuelve el sistema (3) y se obtiene:3, μ3, ν3 Obteniéndose tres sistemas de soluciones particulares: x1 1e r 1t ,
y1 1e r 1t ,
z 1 1e r 1t .
x 2 2 e r 2t ,
y 2 2 e r 2t ,
z 2 2 e r 2t .
x3 3 e 3 ,
y 3 3 e r 3t ,
r t z 3 3 e 3 .
r t
De donde la Solución general del sistema tiene la forma: x = C 1 x1 + C 2 y1 + C 3 z 1 y = C 1 x2 + C 2 y2 + C 3 z 2 z = C 1 x3 + C 2 y3 + C 3 z 3.
Ej emplo 1. Integrar el siguiente sistema de ecuaciones:
dx dt 3 x y z dy x 5 y z dt dz dt x y 3 z
(1).
Solución: Escribimos el sistema para la determinación de λ, μ y ν.
(3 r ) 0 (5 r ) 0........(2) (3 r ) 0
De donde la ecuación característica es:
119
Análisis Matemático 3 r
1 1 3
III
1 1 5 r 1 = 0, desarrollando el determinante obtenemos: 1 3 r 2
r - 11r + 36r – 36 =0, de donde sus raíces son: r 1 = 2, r 2 = 3, r 3 = 6 .
Estos valores reemplazamos en el sistema (2). r = r 1 = 2. Resolviendo tenemos: λ 1 = 1, μ1 = 0, ν1= -1. De esta manera (3 2) 0 encontramos los otros valores. (5 2) 0 (3 2) 0 Si r = r 2 =3. λ 2 = 1, μ2 = 1, ν2= 1. Y cuando r = r 3 = 6. λ 3 = 1, μ3 = -2, ν3= 1. Con estos valores escribimos las soluciones particulares.
x1 e 2t ,
y1 0,
z 1 e 2t .
x 2 e 3t ,
y 2 e 3t ,
z 2 e 3t .
x3 e 6t ,
y3 2e 6t ,
z 3 e 6t .
Y la Solución del sistema será: 2t
3t
6t
x = C 1e + C 2e + C 3 e 3t
6t
y = C 2e - 2C 3 e 2t
3t
6t
z = - C 1e + C 2e + C 3 e .
El caso en que las raíces del polinomio característico sean imaginarias o múltiples queda de deber para el lector.
5.3 Ej ercicios propuestos. dx dt y z dy 1.- x z dt dz dt x y
dx dt y z dy 2.- 3 x z dt dz dt 3 x y
120
Análisis Matemático
III
dx dt 8 y dy 3. - 2 z dt dz dt 2 x 8 y 2 z
dx dt 2 x y 4.- dy 9 x 2 y dt
dx t dt x 3 y e 5.- dy 2 y e t dt
dx dt 4 x y t 1 6.- dy 2 x y t 1 dt
dx 2 2 dt 6 x y 6t t 3 ,... x(0) 2;... y (0) 3 7.- dy 2 y 2t 1 dt
dx t dt x 4 y 3te 8.- dy x y e t dt
dx dy 2 dt dt sen.t 2 y x(0) 2 9.- , y ( 0 ) 1 . dx dy t y dt dt e
dx x (0) 1 dt y 10,- , dy x y (0) 1 dt
121
Análisis Matemático
III
UNIDAD VI 6. Series de Fourier 6.1 I ntr oducción. oducción. Para hablar de este tema recordemos que ya estamos familiarizado con el uso de la series de Taylor y de Maclaurin, donde se expresa una función f x como una serie de potencias en “x”. Hacemos esto porque, con frecuencia, es más fácil tratar con potencias de “x” que con las funciones originales. Desarrollos de las funciones en series de Taylor y Maclaurin: x
x 2
1!
2!
Serie de Maclaurin: f ( x) f (0) f (0) f (0) En torno a x = 0 Serie de Taylor: f ( x) f (a) f (a)
( x a) 1!
f (a)
( x a) 2 2!
......... f (0) ( n)
......... f (a)( ( n)
En torno a x = a, a cualquier valor entero.
( x a) n n!
x n n!
.....
.....
Ejemplos de funciones desarrolladas en serie de potencias: ln(1 x). (1) 1. f ( x) ln( Derivando (1) algunas veces tenemos: f ( x)
1 1 x
f ( x)
f ( x)
1 (1 x) 2 2
(1 x) 3
f v ( x)
6
(1 x) 4 …………………..
Reemplazando en estas derivadas y en la función x = 0 tenemos en serie de Maclaurin: ln1 x x
x 2 2
x 3 3
x4 4
............
Y con x a ; a cualquier valor entero en serie de Taylor: ( x a) ( x a) 2 ( x a) 3 ( x a) 4 ln1 x ln(1 a) .......... .. (1 a) 2(1 a) 2 3(1 a) 3 4(1 a) 4 2. f ( x) cos x. (2) 122
Análisis Matemático
III
Derivando (2) algunas veces: f ( x) sen. x f ( x) cos cos . x f ( x) sen. x f v ( x) cos cos x ……………….
Reemplazando en estas derivadas y en la función x = 0 tenemos en serie de Maclaurin: Maclaurin: cos x 1
x 2
Y con x
2!
2
4!
x6 6!
............
; en serie de Taylor:
( x cos x
x 4
1
)
2
( x
)3
2 ......... 1* 2 * 3
Por ejemplo, en algunos problemas que tratan con oscilaciones, resulta más conveniente usar una serie de senos. Como los senos son periódicos, una serie así sólo puede representar una función periódica. En el párrafo siguiente se recordará al estudiante algunas de las características de las funciones periódicas.
Defi ni ciones ciones prelimin prelim in ares. ares. 6.2 Defini Defi ni ción ción de fun ciones per per iódicas iódicas: :
Una función f (t) es periódica o tiene un período T si para todo t , f (t +T ) = f (t ), ), donde T es una constante positiva. El valor más pequeño de T >0 se llama llama el período principal o período fundamental o simplemente el período de f ( (t ). ). Ej emplo 1.
La función sen(x) tiene períodos, 2π, 4π, 6π ,…, ,…, y las funciones sen(x+2π), sen(x+4π), sen(x+6π),.., van a tener el mismo periodo de sen(x), porque sen(x+nπ) donde n es un número entero par va a ser siempre igual a sen(x). Y = se n.x Y
1,5 1 0,5 0 -0,5 0 -1 -1,5
X
2
4
6
8
F igura 6.1 Ej emplo 2.
123
Análisis Matemático
III
Dada la gráfica, determinar si la función es periódica, que periodo tiene y cual su ecuación. Y 3 2 1 0 -6
-4
-2
-1 0
X 2
4
6
8
-2 -3
F igur a 6.2
La función si es periódica, de periodo 4 y corresponde a la función f(x) = x Ej emplo 3. Si una función tiene periodo 2π, entonces F (t ) f ( .t ) f (
2 T
.t ), lo que significa que para
calcular el periodo de cualquier función periódica debemos utilizar la formula
2 T
Ej emplo 4.
Si a R y f tiene período T entonces entonces se cumple:
∫ ∫ ∀ ∈ Ej emplo 5.
Encontrar el periodo de la función
2 2
Si la función f (t) es periódica con un periodo T, entonces, de tiene que:
por lo que se
Puesto que cos ( + 2m) = cos para cualquier entero m se tiene que
Donde m y n son enteros, Por consiguiente T = 6 m; cuando m = 4 y n = 3, se obtiene el mínimo valor de T. (esto se puede ver mediante el procedimiento de ensayo error y error). De donde, T = 24 en general, si la función f(t) = cos co s ω 1t + cos ω2t es periódica con período T, entonces es posible encontrar dos enteros m y n tales que:
124
Análisis Matemático
III
1T
= 2nm 2 T = 2mn el cociente es
Es decir, la relación 1 / 2 debe ser un numero racional. Ej emplo 6.
Decir si la función f (t) = Cos10t + Cos (10 + π)t es una función periódica. Aquí tenemos que ω1 = 10 y ω2 = 10 + π. Aplicando el enunciado anterior tenemos:
Como observamos no es un número racional, por lo que es imposible encontrar un valor T que satisfaga f (t) = f (t+T) por consiguiente f (t) no es una función periódica. Ej emplo 7.
Encontrar el periodo de la función f (t) = (10cost)2 Si aplicamos la identidad trigonométrica cos2t = 1/2 (1+cos2t) se tiene:
10 100100 1cos25050cos2
Puesto que una constante y una función periódica de periodo T para cualquier valor de T, el período de cos 2t es, se concluye que el periodo de f(t) es . Demostrar que si f(t + T) = f(t), entonces.
∫ ∫ ∫ ∫ Definición: Una función periódica f x cuyo período es 2 , puede expresarse en la forma f x c0 c1 sen x 1 c2 sen2 x 2 c3 sen3x 3 ....... (2.1) siempre que f x satisfaga otras condiciones determinadas. Para las funciones que probablemente encuentre el estudiante, de las únicas condiciones que deberán estar seguros que se cumplen son: 1.- La función debe ser de valor único, es decir, para cada valor de x debe haber un solo valor de f x , es decir la función debe ser Biyectiva. 2.- La función tiene que ser finita y nunca infinita Una serie como la del segundo miembro de (2.1) se conoce como SERIE TRIGONOMÉTRICA DE FOURIER ya que existen muchos tipos de estas series.
125
Análisis Matemático
III
El término c1 sen x 1 se llama primera armónica o fundamental. Los términos c2 sen2 x 2 , c3 sen3x 3 , etc. se conocen como segunda armónica, tercera armónica, etc. Si el período de f x no es 2 , mediante un cambio apropiado de variable puede hacerse una conversión hacia una función de período 2 , como se verá posteriormente en el programa.
Forma alternativa de la serie (2.1): para ponerle en esta forma utilizamos la trigonometría como sigue: sen x 1 senx cos 1 cos xsen 1 Y, por tanto, puede escribirse c1 sen x 1 a1 cos x b1 senx
Donde a1 c1 sen 1 y b1 c1 cos 1 Convirtiendo los otros términos en una forma semejante, la serie queda: f x
1 2
a0 a1 cos x a 2 cos 2 x a3 cos 3x ........
b1 senx b2 sen2 x b3 sen3x ........
Posteriormente quedará clara la razón para escribir c 0 como
(3.1) 1 2
a0 .
Ahora, nuestro problema, en cualquier caso particular, es encontrar los coeficientes lo que es equivalente, los valores de a0 , a1 , a2 ,......, b1 , b2 ,....... o, c0 , c1 , c2 ,......, 1 , 2 ,...... En la práctica, es más fácil calcular los a y b, entonces, pueden hallarse los c y los , a partir de ellos, si se requiere la serie en la forma (2.1). Pero en general es posible que la forma (3.1) sea apropiada por sí misma. Una manera más concisa de representar la serie trigonométrica de Fourier es:
2 ∑ 12
Aplicaciones : Algunos procedimientos que utilizan los ingenieros en la solución de problemas, sólo son ajustables a funciones senoidales y, consecuentemente, a menudo resulta útil usar las series de Fourier para expresar otras funciones en términos de senos (y o cosenos).
a.- En el campo de la ingeniería eléctrica. Dos de estos casos son: El uso de números complejos al tratar con circuitos de c.a; depende de funciones senoidales.
126
Análisis Matemático
III
Si se requiere el valor promedio de la potencia en un circuito sobre el tiempo T, siendo T el periodo del voltaje aplicado, tiene que evaluarse la integral , que se simplifica, si el voltaje y la corriente pueden escribirse en funciones de cosenos de arcos múltiples.
∫
b.- En ingeniería civil. Son útiles al abordar algunos problemas de vigas, en la determinación y construcción de la curva de deflexión, cuando éstas soportan cargas que son de forma irregular, y que no resulta fácil hallar la solución de las ecuaciones estándar de las vigas:
c.- En la ciencia médica.
Una nota musical consiste de una armónica fundamental y armónicas superiores, y, en el estudio de la fisiología del oído, la teoría de Helmholtz supone que las diferentes fibras en el oído interno sólo son excitadas por aquellas armónicas que están en armonía. La actividad del corazón produce una corriente eléctrica. La variación del voltaje que resulta es periódica y puede medirse mediante un electrocardiógrafo. La onda resultante puede ser analizada en una serie de Fourier y los coeficientes mostrarán si existe alguna enfermedad o daño. La temperatura humana fluctúa periódicamente y, por tanto puede representarse como una serie de Fourier.
d.- Otras de las aplicaciones es su uso en ajustar condiciones en la frontera a las soluciones de ecuaciones parciales.
6.3 I ntegrales requeri dos para calcular los coeficientes de la serie (3.1) El cálculo de los coeficientes a0 , a1 , a2 ,......, b1 , b2 ,...... requerirá los valores de ciertas integrales definidas, los cuales citaremos ahora. En estas integrales tanto m como n son enteros. Dos de las integrales que se requerirán son cos nxdx y sennxdx . (I). Como el
estudiante sabe, la integral de un seno o un coseno sobre un período completo es cero. Esto es obvio observando sus gráficas (si el lector no está satisfecho sobre este punto, un esquema sencillo debe convencerlo). El período de cos nx es 2 n , por lo tanto, integrando desde hasta
, integramos sobre n períodos completos. De donde tendríamos:
1 cos nxdx sen.nx n
1
2
sen.n sen.n sen.n 0 , para cualquier n. y de la n
n
misma manera podríamos demostrar que sennxdx = 0, cuyo periodo es 2 n . 127
Análisis Matemático
III
Las otras integrales que se necesitarán para determinar los coeficientes de estas series son las combinaciones: , y cos mx cos nxdx senmxsennx dx senmx cos nxdx . (II)
Para la primera integral, se tiene
cos mx cos nxdx
1 2
cosm n x cosm n xdx
=
0
si m n , ya que ambos
cosenos se integran sobre un cierto número de períodos completos. 1 1 1 Pero si m n , la integral queda cos 2nx 1dx sen2nx x 2 2 2n Para el estudiante queda la resolución de las otras dos integrales, cabe recordar que se debe considerar que m n como un caso especial.
6.4 Cálcul o de l os coeficientes. Ahora ya estamos listos para encontrar las fórmulas para calcular los coeficientes a y b. Primero, mostraremos cómo hallar a 0 . Ahora bien, escribimos la serie en la forma: f x
1 2
a0 a1 cos x a2 cos 2 x ... an cos nx ... b1 senx b2 sen2 x ... bn sennx ... (3.1).
Y vamos a integrar ambos miembros desde hasta por lo que tenemos:
f x dx
1 a0 a1 cos x a2 cos 2 x .. an cos nx b1 senx b2 sen2 x ..bn sennx dx 2
1
2
f x dx
a0 dx
Las integrales de los demás términos serán cero, ya que cada una es un caso de (I) o bien, de (II). Resolviendo tenemos: Nótese que hasta
1
2
a0
f x dx
1 2
1 2
a0 x
1 2
a 0 2 a0 . a 0
1
f x dx
f x dx , es igual al valor medio de f x sobre el intervalo
.
Ahora vamos a encontrar los coeficientes de los términos cosenoidales.
128
Análisis Matemático
III
Para hallar a1 , multiplicaremos ambos miembros de (3.1) por cos x e integramos desde hasta . Obteniendo:
f x cos xdx
a1 cos 2 xdx
Una vez más, las integrales de todos los demás términos serán cero, ya que cada una es un caso de (I), o (II). Resolviendo:
f x cos xdx a1
1
a1
1
2
1
1
(1 cos 2 x)dx a1 x sen2 x a1 . 2 4
f x cos xdx
Ahora bien, Para hallar a2, a3, etc., multiplíquense ambos miembros de (3.1) por cos 2x, cos 3x, etc., respectivamente, e integramos entre y , en general, para encontrar an, multiplicamos ambos miembros de (3.1) por cos nx e integramos entre y por lo que tendremos:
f x cos nxdx
a n cos nxdx a n 2
1
2
(1 cos 2nx)dx a n
1 2
x
1 2n
an .
sen2nx
De donde: an
1
f x cos nxdx .
La misma que nos servirá para calcular todos los términos de a. A continuación vamos a encontrar los coeficientes de los términos senoidales. Para hallar b1 , multiplicaremos ambos miembros de (3.1) por sen x e integramos desde hasta . Obteniendo:
f x senxdx
b sen 1
2
xdx
Una vez más, las integrales de todos los demás términos serán cero, ya que cada una es un caso de (I), o (II). Resolviendo:
b1
f x senxdx b1
1
1
2
1
1
(1 cos 2 x)dx b1 x sen2 x b1 . 2 4
f x sen. xdx
Ahora bien, Para hallar b2, b3, etc., multiplíquense ambos miembros de (3.1) por sen 2x, sen 3x, etc., respectivamente, e intégrese entre y , en general, para encontrar bn,
129
Análisis Matemático
III
multiplicaremos ambos miembros de (3.1) por sen nx e integraremos entre y obteniendo:
f x sen.nxdx
b sen n
2
nxdx bn
1
2
(1 cos 2nx)dx bn
1 2
x
1 2n
bn .
sen2nx
De donde: bn
1
f x sen.nxdx .
Fórmula que nos servirá para calcular todos los términos de b. Como dijimos anteriormente que era más fácil calcular los términos a y b y a partir de estos calcular los c vamos a poner las siguientes equivalencias:
12 12 − 12 12 2 − 12 Ej emplo 1.
Calculando la serie de Fourier de la función f dada por:
0≤ −1 1 ≤2 F igur a 6.3
Solución: 1
2
a0
an
1
1 2 (
0
(
0
dx
2
cos nxdx
dx)
2
1 2
( x
o
cos nxdx )
x
2
)
1 2
1 1 ( sen.nx n
0
( 2 ) 0. 1
sen.nx n
2
) 0.
130
Análisis Matemático bn
1
( sen.nxdx 0
III
2
sen.nxdx )
1
1 ( cos .nx n
1 1 1 (cos n 1) (cos 2n cos n ) n n
1
cos .nx
0
n
2
)
4 ........n.impar n 0.............n. par
Por lo que reemplazando valores la serie será: f ( x)
4 sen.3 x sen.5 x ...... , sen. x 3 5
Ej emplo 2.
Calcular la serie de Fourier de la función f dada por: f ( x) e x , en el intervalo ,
F igur a 6.4
Solución: 1
2
a0
1 2
(
e x dx)
1 2
(e x )
2
(e e )
shx
.
1 a n ( e cos nxdx) (cos nx nsen . x ) e cos n e cos n 2 2 1 n (1 n ) 2 sh (1 n 2 ) ............ si..n..es.. par e e cos n . 2 sh (1 n 2 ) ........... si..n..es...impar . (1 n 2 ) 1
x
1 e x
1
x
1 e x
1 e n cos n e n cos n bn ( e sen.nxdx ) ( sen . x n cos nx ) 2 2 1 n (1 n ) 1
131
Análisis Matemático
III
2nsh (1 n 2 ) ............ si..n..es.. par n(e e ) cos n . 2 nsh (1 n 2 ) ........... si..n..es...impar . (1 n 2 )
Por lo que reemplazando valores la serie será: f ( x)
2 sh 1
1 1 1 1 2 3 cos x cos 2 x cos 3 x ........ senx sen2 x sen3x .... 2 2 5 16 2 5 10
6.5 Ser ie de Four ier de funciones pares e impares. Si nosotros de antemano pudiéramos analizar si la función de la que vamos a encontrar una serie de Fourier, es par o impar, podríamos reducir el trabajo a realizar en el cálculo de los coeficientes a la mitad, como sucede en el ejemplo uno. Entonces será útil, poder reconocer si una función es impar o par. Por lo que a continuación vamos a recordar algunas definiciones y propiedades importantes de las funciones pares e impares. 6.5.1 D efi ni ciones.
La función f(x) se define como par si f (-x) = f(x). De esta definición concluimos que la grafica de estas funciones es simétrica respecto al eje de las ordenadas “y”, como ejemplo podríamos citar a cualquier potencia par de “x”, del mismo modo cos x y, más generalmente, cos nx. La Función f(x) se define como impar si f (-x) = -f(x). Gráficamente podemos decir que estas funciones son simétricas respecto del origen, como ejemplo podríamos poner cualquier potencia impar de “x”, sen x y, más generalmente, sen nx. También vale recordar que el producto de dos funciones pares o dos impares es par y que el producto de una función par por una impar es impar. Ejemplos: Analizar y decir cuál de las siguientes funciones son pares, impares o de ninguno de estos tipos. 1. 2. 3. 4. 5.
f (x) = x cos x 2 f(x) = x cos 5x f(x) = (x+π) cos x f(x) = sen x sen 2x.
132
Análisis Matemático
III
F igura 6.5 6.5
6.
F igura 6.6 6.6 7.
F igura 6.7 6.7 8.
F igura 6.8 6.8
Solución: 1. f (x) = x cos x f (-x) = (-x) cos(-x) = -x cos x = - f(x), por lo que x cos x es una función impar. impar. 2 2. f(x) = x cos 5x
133
Análisis Matemático 2
III 2
f (-x) = (-x) cos(-5x) = x cos 5x = f(x), f (x), este producto es una función par. par. 3. De ninguno de los dos tipos. 4. Par 5. Par 6. Impar 7. Impar 8. Par
Al trabajar trabajar con las series de FOURIER FOURIER con frecuencia vamos a encontrarnos con integrales h de la forma f ( x)dx , por lo que vamos a definir algo importante referente a las funciones h
pares e impares. 1. Si f(x) es una función impar, 2. Si f(x) es una función par,
h
h
h
h
f ( x) dx 0.
h
f ( x)dx 2 f ( x)dx 2 0
0
h
f ( x)dx pero esta última no
es muy conveniente por tener límites negativos. Las series de Fourier para una función par solo puede contener términos que sean también funciones pares, es decir, solo tendrá términos cosenoidales y será de la forma: 1 2
a 0
a n cos nx
n 1
Por otro lado, la serie de Fourier para una función impar sólo puede contener términos que también sean impares, es decir, sólo tendrá términos senoidales y será de la forma:
b n sen nx
n 1
Ejemplo 1 Consideraremos la onda dientes de sierra representada en la siguiente gráfica.
F igur a 6.9
Esta onda esta descrita por la ecuación: f(x) = x en
- < x < , además esta función es
periódica ya que f (x+2π) = f(x)
Como f (-x) = -f(x) decimos que es una función impar. Por lo que al calcular esta función en serie de Fourier los coeficientes a 0 = a n = 0. Y solo calcularemos los coeficientes bn.
134
Análisis Matemático
III
− cos ∫ coscos 2 − cos 1 [1 ] − cos −
bn =
1
2
x. sen.nxdx , x. sen.nxdx = 0
función par
, haciendo un análisis para los diferentes valores que puede tomar n
Por lo que la serie de Fourier para esta onda es:
1 1 1 1 f(x) = 2 senx sen2 x sen3 x sen4 x sen5 x ...... 2 3 4 5
Ejemplo 2 Obtenga la serie de Fourier para la función periódica que se muestra en la figura.
F igur a 6.10 6.10
(− )
Esta Función está definida por: f(x) = 1 en el intervalo de los valores de b n = 0, y sólo calcularemos los términos de Solución: 1 2
=
,
a0 = valor medio de f (x) sobre un período 1 2
135
Análisis Matemático
a n =
1
f ( x) cos nxdx
=
1
2
1. cos nxdx
=
III
2
2
2
cos nxdx 0
=
n sen 2 n 2
0 − 1511 La serie de Fourier es
1 2
2 1 1 1 cos x cos 3 x cos 5 x cos 7 x....... 3 5 7
A veces, puede suceder que una función que no es impar ni par, es simétrica respecto a un punto sobre el eje, y que no sea el origen. Si este punto se toma como un origen temporal, entonces la función puede tratarse como impar. La función de onda cuadrada que se ve en la figura, puede tratarse de esta manera. El estudiante verá que es simétrica respecto al punto (0,
1 2
), si se toma como nuevo origen este
punto, como se muestra en diagrama que sigue, entonces tendremos una función impar g (x).
F igur a 6.11
Los ejes originales se indican mediante ---------------.
136
Análisis Matemático
III
g (x) está dada por las ecuaciones:
− − 0 2
Ahora vamos a encontrar la serie de Fourier para g(x), al observar la gráfica con el nuevo origen de esta función g (x) vemos que es impar por lo que los coeficientes a 0 = a n = 0 por lo que sólo procedemos a calcular los bn. b n = =
1
g ( x ) sennxdx
2
1 sennxdx ya que g (x) sen nx es una función par 0 2
1 1 = cos nx n 0 =
1 n
1 cos n
0 2
De donde: g (x) =
2 1 1 1 senx sen3 x sen 5 x sen7 x ..... . 3 5 7
Volviendo al origen inicial, f ( x )
1 2
g ( x )
= 1 2 senx 1 sen3 x 1 sen5 x 1 sen7 x .....
Para concluir diremos que en esta ocasión solo tuvimos que calcular b n y que el valor de 2
1 2
3
5
7
a 0 es la distancia del nuevo origen al original.
6.6 Series de F our ier de medio rango. En algunos problemas donde la función no es periódica dentro del intervalo de definición, resulta útil representarla mediante una serie de Fourier. Para sortear esta dificultad, usamos una función que es periódica y que coincida con la función dada sobre el intervalo de 137
Análisis Matemático
III
definición. El periodo que se elija, obviamente debe ser mayor o igual al intervalo dado. Tomando el intervalo como la mitad de un período, hace posible obtener una serie sencilla, definiendo la función para otro semiperiodo en tal forma que sea impar, dado por tanto una serie solo de senos, o bien par, dando una serie de cosenos (incluyendo el término que puede considerarse como
0
1 2
a 0 , ya
Por ejemplo, si tenemos la función definida entre 0 y como se muestra en la figura.
F igura 6.12
Puede elegirse cualquiera de las ondas siguientes para representarla.
F igur a 6.13
F igura 6.14
138
Análisis Matemático
III
Figura 6.13.- Una onda impar y su serie sólo contendrán senos (conocida como serie senoidal de Fourier de medio rango). Figura 6.14.- Es una onda par y su serie sólo contendrán cosenos (llamada serie cosenoidal de Fourier de medio rango). En la serie de Fourier de medio rango, el cálculo de los coeficientes se simplifica mediante la suposición de que sea una función impar, o bien par, la que representa la función requerida sobre el intervalo dado. Si suponemos una función impar a f (x), entonces a 0 = a n = 0, y calculamos solamente b n=
bn
1
f ( x ) sennxdx
=
2
f ( x) sennxdx
ya que f (x) sen nx es par.
0
Si por lo contrario suponemos una función par a f (x), entonces b n =0, y calculamos a 0 =
1
2
f ( x)dx f ( x)dx 0
y a n =
1
2
f ( x) cos nxdx f ( x) cos nxdx 0
ya que f (x) cos nx es par.
Volviendo ahora al ejemplo que hemos usado anteriormente, calculemos la serie para el caso Figura 6.13 En el intervalo 0 a , la función está definida por las ecuaciones 1 2 x f ( x) x
0 x
2 3
2 3
x
2 3 1 Por lo que: bn = f ( x ) sennxdx = xsennxdx ( x) sennxdx 0 2 2 0 3
2
139
Análisis Matemático
bn
2 1 x
III 2 / 3
= cos nx 0 2 n
1 2
2 / 3
0
1 x cos nx cos nxdx cos nxdx n n 2 / 3 2 / 3 n
1
2 2n 1 2n 2n 1 2n = cos 2 sen cos 2 sen 3 3 3n 3 3 3n 2n n
bn =
3 .n
sen 2
2n 3
1 ⁄ √ 2 −√ 0 ⁄2 25 De donde: f ( x)
3 3 1 1 1 senx 2 sen2 x 2 sen4 x 2 sen5 x ... 2 2 4 5
Ejemplo: Para una función f (x) entre 0 y , definida `por la ecuación f ( x) x . Realizar un esquema de las ondas que resultaran, si esta función se representa mediante una serie de Fourier de medio rango que contenga. (1)
senos,
(2)
cosenos
Y obtener la serie en cada caso.
F igura 6.15
140
Análisis Matemático
b
=
n
2
III
2 1
0
0
2 n ; dando valores a n
( x) sennxdx cos nx( x ) cos nxdx n 0
1 1 1 tenemos: f ( x) 2 senx sen2 x sen3 x sen4 x ... . 2 3 4
F igura 6.16 a 0 =
a n =
2 2
2
2
2 1
x cos nxdx sen.nx( x) n
0
0
0
De donde dando valores a n tenemos: f ( x)
2
1 n
2
n 2
sen.nxdx
(1 cos n ).
4 1 1 cos x cos 3 x cos 5 x ... 9 25
6.7 Ar mónicas im pares y pares.
Hemos visto que, si una función es impar o par, resulta útil observar esto a partir de la grafica antes de empezar a calcular los coeficientes. Otros tipos especiales de onda que pueden ocurrir, son aquellas cuyas series de Fourier sólo contienen armónicas pares o bien, armónicas impares. Nuevamente, resulta útil poder observar, a partir de la gráfica, que esto va a suceder. Si
f(x)
=
1 2
a 0 a1 cos x a 2 cos 2 x a3 cos 3 x ... b1 senx b2 sen2 x b3 sen3x ... ,
entonces, reemplazando x por x , tendremos
141
Análisis Matemático f ( x )
1 2
III
a0
a1 cos( x ) a2 cos(2 x 2 ) a3 cos(3 x 3 ) ... b1 sen( x ) b2 sen(2 x 2 ) b3 sen(3x 3 ) ...
1
a0 a1 cos x a2 cos 2 x a3 cos 3 x ... b1 senx b2 sen2 x b3 sen3x ... 2
1
Si la serie para f (x) sólo contiene armónicas pares (se incluye a 0 , ya que puede considerarse 2
como
1
a 0 cos 0x), puede verse que 2 f ( x ) f ( x)
Por otra parte, si sólo contiene armónicas impares, entonces f ( x ) f ( x)
A continuación mostramos un ejemplo de cada caso.
F igura 6.17
F igura 6.18
Como hemos podido darnos cuenta que en el primer caso, la función, de la misma manera en que se repite en intervalos de 2, también se repite en intervalos de, como puede concluirse a
142
Análisis Matemático
III
partir de la relación f ( x ) f ( x) . Esto pudiera tomarse también como un ejemplo del caso más general, donde el período es diferente de 2, lo cual consideramos brevemente. A continuación ponemos ejemplos para ver si las series que representan las funciones siguientes contendrán a) sólo armónicas pares, b) sólo armónicas impares, c) tanto armónicas pares como impares.
F igura 6.19
F igura 6.20
F igura 6.21
Solución: 1.Sólo armónicas impares. 143
Análisis Matemático 2.3.-
III
Sólo armónicas pares. (La onda senoidal rectificada) Ambas.
A continuación un ejemplo de aplicación: Un voltaje generado por una máquina rotatoria es un ejemplo de una onda en ingeniería cuya serie de Fourier sólo contiene armónicas impares. Como resultado de la semejanza entre los polos magnéticos n y s en cualquiera de las máquinas, el voltaje generado siempre tiene este tipo de simetría. 6.8 Suma de un a serie de F our ier en un pun to de disconti nuidad.
Si se sustituye un valor particular de x en una serie de Fourier, la suma que resulta es el valor correspondiente de f (x). Sin embargo, si se hace esto en un punto de discontinuidad x = k, digamos, la suma es el promedio de los dos valores límite de f (x) conforme x k desde la izquierda y derecha, respectivamente así si k es el punto de discontinuidad la suma de la serie de Fourier es: 1 2
f ( x 0) f ( x 0)
Ejemplo:
Sea la onda
F igura 6.22 2 1 1 1 senx sen3 x sen5 x sen7 x ... 2 3 5 7 Sustituyendo, por ejemplo, x en la serie tendremos:
Su serie de Fourier es:
1
1 1 2 1 1 1 1 sen sen sen sen ... , Lo cual es (1 0) . 2 3 5 7 2 2
6.9 Periodo general .
A continuación vamos a ver funciones que tienen el periodo diferente 2, ya que hasta el momento sólo consideramos funciones de dicho periodo, pero los problemas mas prácticos
144
Análisis Matemático
III
comprenden el uso de dichas funciones con un período diferente. Veremos ahora un ejemplo de una función cuyo período no es 2. Sea la onda:
F igura 6.23
Cuya ecuación es y f ( x) , definida por
1 2 −20 2 0 00
De hecho, esto puede convertirse al caso directo, donde el período es 2, mediante un cambio en la escala horizontal como se ilustra a continuación.
F igura 6.24
Tratemos ahora a t como una nueva variable independiente y, por tanto, la ecuación de la onda puede escribirse en la forma y F (t ) . A partir de la gráfica, escriba: (i) Las ecuaciones que definen F (t). (Ignore la escala x mientras hace esto.) (ii) La ecuación que da t en términos de x.
145
Análisis Matemático
III
1 f ( x ) (t )....... t 0 F (t ) 0...................0 t F (t 2 ) F (t )
t
x
2
Ahora tenemos una función F (t) la cual tiene periodo 2. Mostrar que la serie de Fourier para está es F (t )
a0
an
1
4
1
2
2 2
1 1 1 1 1 cos t 2 cos 3t 2 cos 5t ... sent sen2t sen3t ... 3 5 2 3
1 2
1
F (t ) cos ntdt 0
1 1 = (t ) cos ntdt
=
1 2 n 2
1 cos n
0
1 bn F (t ) sen.ntdt
=
1
0
1
.ntdt t sen
146
Análisis Matemático =
III
1 n
F (t )
1 2 2 4
1 1 1 1 1 cos t 2 cos 3t 2 cos 5t ... sent sen2t sen3t ... 3 5 2 3
Ahora bien, F (t ) y f ( x) y f ( x)
1 4
t
x
2 1
, de modo que la serie de Fourier para f (x) es
2 x 1 3 x 5 x 1 3 x 1 x 1 2 cos 2 cos ... sen sen x sen ... cos 2 2 3 2 2 2 2 3 2 5
Viendo cómo tratar un ejemplo de este tipo, en el cuadro siguiente consideraremos el caso general con período 2 , donde puede tomar cualquier valor. Cuando el período es 2 , tomemos nuevamente una nueva escala horizontal, pero en esta ocasión elegida de tal manera que t 2 corresponda a x 2 . 3 -3
-2 -2
- -
0 0
2 2
3 3
t x
Escriba la ecuación que da t en términos de x. t
x
y f ( x) Se convierte ahora en y F (t ) con período 2 y la serie de Fourier par F (t) será 1 2
a0 a1 cos t a2 cos 2t a3 cos 3t ... b1 sent b2 sen2t b3 sen3t ...
Donde a0 an
1
F (t )dt 1
F (t ) cos ntdt
1 bn F (t ) sen.ntdt
La serie y las formulas para los coeficientes pueden escribirse en términos de x, usando t
x
, para dar la serie de Fourier para f (x).
147
Análisis Matemático
III
De
donde:
f ( x)
a0
1
a0 a1 cos
2 1
x
a2 cos
2 x
a3 cos
3 x
... b1 sen
x
b2 sen
f ( x) dx ya que dt dx y x corresponde a t
1
2 x
b3 sen
3 x
...
f ( x)dx
an
1
bn
1
f ( x) cos
f ( x) sen
n x
n x
dx
dx
Estas formulas pueden considerarse como las estándar para los coeficientes cuando el período es general, de modo que no es necesario hacer la sustitución con t en cada ocasión. Anteriormente se señaló que para una función de periodo 2, las integrales usadas al evaluar los coeficientes pueden ser tomadas sobre cualquier intervalo de longitud 2. De modo semejante, aquí las integrales pueden tomarse sobre cualquier intervalo de longitud 2 . Por ejemplo, pueden usarse los límites 0 y 2 , en lugar de y . Como ejemplo, consideremos la onda que se muestra:
F igura 6.25
Primero, escribimos las ecuaciones que definen la función y luego hallamos el valor de este caso tenemos:
12 −0 0
en
0 12
148
Análisis Matemático
III
El período es de 6, de modo que el valor de es 3. Escribiendo, las formulas para a 0, a definiciones anteriores tenemos: a0
1
0
3
3
1
0
n x
1
3
1
3
12dx 4 xdx ( 3
0
n
y b
n,
para esta onda particular, usando las
x Área bajo la gráfica entre -3 y +3)
n x dx 4 x cos dx 12 cos 3 3 3 3 0 0 3 1 n x n x bn 12 sen dx 4 xsen dx 3 3 3 3 0 an
a0
3
1 4 xdx 12 dx ( x Área bajo la gráfica entre 0 y 6) 3 3 0 3 6
n x dx 12 cos dx 4 x cos 3 0 3 3 3 3 6 1 n x n x bn 4 xsen dx 12 sen dx 3 0 3 3 3 Procediendo a evaluar los coeficientes con el primer conjunto de ecuaciones tenemos: an
1
3
6
n x
1 a0 .54 18 3 n x dx x cos dx 3 cos 3 3 3 3 0 0 3 3 4 n x n x 3 n x 9 3 sen x sen sen dx n n 3 n 3 3 3 3 0 0
an
4
0
n x
3
3
3 n x cos n n 3 0 4
12 n 2 2
(cos n 1)
0 −
=
Vamos a demostrar ahora que bn
12 n
.
149
Análisis Matemático
bn
bn
1
0
4
0
III
n x dx 4 xsen dx 12 sen 3 3 3 3 0 3
n x
3 sen
n x
3 3
3
3
dx xsen
n x
0
3
dx
4 n x n x 3 n x 9 3 cos x cos cos dx 3 3 3 n 3 0 0 n 3 n 0
3
3
4
9
3 n 12
cos n 1 cos n
n
Por lo que dando los valores n en las expresiones encontradas, la serie de Fourier de la función del ejemplo es: 24 x 1 1 5 x 2 x 1 12 x 1 cos cos x cos ... sen sen sen x ... 2 3 32 52 3 3 2 3 3
f ( x) 9
Ahora veremos que fácilmente la idea de una serie de Fourier de medio rango puede extenderse a otro medio que no sea. De la siguiente manera. Para una función impar: bn
y
a 0 = a n = 0
2
0
f ( x) sen
n x
dx
Para una función par. bn 0
an
a0
,
2
2
f ( x) cos
f ( x)dx 0
n x
0
dx
Ej emplo 1.
Encuentre la serie senoidal de Fourier de medio rango para x 2 en el intervalo 0 < x < 3 y haga un esquema para la onda representada por esta serie sobre varios períodos. bn
2
3
x sen 3 2
0
n x 3
dx
150
Análisis Matemático
2 3 x 2
III
n x 3 cos 2 x cos dx 3 n 3 0 0 n 3 n x
3
3
2 27
6 3 x n x 6 3 n x cos n sen sen dx n n n 3 n n 3 3 0 0 2
3
3
18 3 n x cos n 2 2 cos 3 n n 3 n 0 2 27 54 cos n 2 2 (cos n 1) 3 n n 3
27
− − 1 La serie es de 18
x 4 1 2 x 1 4 1 4 x sen sen 3 3 sen x sen ... 3 3 2 3 4 3 3 3
F igura 6.26 6.10 Análisis armóni co..
Las formulas obtenidas para calcular los coeficientes en los apuntes anteriores, comprendieron integrales cuya evaluación depende de conocer la ecuación de y en función de x. Sin embargo, ésta puede no conocerse y, en su lugar, sólo se den o puedan medirse valores numéricos de y para ciertos valores de x. En tales casos, es necesario modificar las formulas para los coeficientes, de manera que pueda usarse un método numérico para su evaluación. Ej emplo 1.
Si se aplica un voltaje alterno a un circuito que consiste de un rectificador, una resistencia y una inductancia en serie, la corriente es de la forma 151
Análisis Matemático
III
F igura 6.27
En un caso particular, los valores de i, medidos en 12 valores equidistantes de que cubren un período de 2, fueron: 0 0
i
/6
/3
/2
2.3
5.5
8.9
2/3 5/6 10.8 11.4 9.9
Los coeficientes de la serie de Fourier, a 0, a siguiente manera: a0
an
bn
1
1
1
n
y b
7/6 4.8 n
4/3 0
3/2 0
5/3 0
11/6 0
en estos casos se encuentran de la
2
f ( )d
=2
valor medio de f ( ) sobre un período.
0
2
f ( ) cos n d
=2
valor medio de f ( ) cos n sobre un período.
=2
valor medio de f ( ) sen n sobre un período.
0
2
f ( ) senn d 0
No se conoce f ( ), de modo que se encuentra una aproximación para el valor medio de f ( ), tomando la media de los 12 valores de i en las tabuladas, es decir
i 12
De manera semejante, el valor medio de f ( ) cos n se toma como
i cos n 12
y el de f ( ) sen n como
152
Análisis Matemático
III
isenn 12
De donde a0 2 ; a1 i
12
i
6
i cos 6
; b1
isen ; etc. 6
El cálculo se plantea mejor en forma de tabla. Por conveniencia, cambiaremos a grados
i
i sen Cos i 2 Cos 2 0 1 0 0 0 1 0 2.3 0.866 1.99 0.5 1.15 0.5 1.15 5.5 0.5 2.75 0.866 4.76 -0.5 -2.75 8.9 0 0 1 8.9 -1 -8.9 10.8 -0.5 -5.4 0.866 9.35 -0.5 -5.4 11.4 -9.87 0.5 5.7 0.5 5.7 9.9 0.866 -9.9 0 0 1 9.9 4.8 -1 -4.16 -0.5 -2.4 0.5 2.4 0 0 0 -0.5 0 0 0.866 0 0.866 0 -1 0 0 -0.5 0 -1 0 -0.5 0 0 0 0 0 0.5 0 0.5 0.866 0.866 -0.5
0 30 60 90 120 150 180 210 240 270 300 330
Cos
53.6
24.59 -4.1
/6 8.9
a 0 = 8.9
a 1 = - 4.1
b 1 = 4.6
i Cos sen
a 2 = 0.3
Sen 2
i sen 2 0 1.99 4.76 0 -9.35 -9.87 0 4.16 0 0 0 0
0 0.866 0.866 0 0.866 0.866 0 0.866 0.866 0 0.866 0.866
27.46
2.10
-8.31
4.6
0.3
-1.4
b 2 = - 1.4
Continué ahora la tabla para encontrar a 3 , b 3 , a 4 y b 4. a 0 = -0.8
b 3 = 0.0
a 4 = 0.2
b 4 = 0.1
Otro calculo más da a 5 = - 0.1, b 5 = - 0.1 De donde la serie es: f(x) 4.5 - 4.1cos + 0.3 cos 2 - 0.8 cos 3 + 0.2 cos 4 - 0.1 cos 5 ... 4.6 sen - 1.4 sen 2 ... 0.1 sen 4 - 0.1sen 5 .....
153
Análisis Matemático
III
Esto también puede expresarse en la forma: F(x) = 4.5 + 6.2 sen ( - 40) + 1.4 sen (2 + 168) + 0.8 sen (3 - 90) + 0.2 sen (4 + 27) + 0.1 sen (5 - 135) +… Como se nota las amplitudes de las armónicas sucesivas decrecen, de modo que las armónicas superiores contribuyen muy poco al total. El método aproximado que hemos usado para calcular los a i y b i se vuelve menos exacto conforme las armónicas se hacen mayores. La única manera de mejorar la exactitud del método es decrecer el intervalo entre los valores tabulares, lo cual, por supuesto, significa incrementar el número de lecturas. Como una guía general, si el período se divide en k k intervalos, el método solo es confiable en la estimación de las primeras 1 armónicas.
2
Por ejemplo, en el presente problema, con 12 intervalos, sólo se puede esperar ir hasta la quinta armónica e incluso entonces, por lo regular, sólo consideraríamos las amplitudes de la cuarta y quinta armónicas relativas a la de la fundamental, en lugar de tratar cada una como exacta en sí misma. 6 .11 Ejercici os propuestos.
1.- Una función f(x) de periodo 2l está definida por:
− − 0 0 | | − −− −−ℎℎ cos (− 2 )
Hallar la serie de Fourier. 2.- Hallar la serie de Fourier que represente la onda senoidal rectificada 3.- Desarrollar las siguientes funciones en serie de Fourier en los intervalos dados: .
.
4.- Hallar la serie de Fourier senoidal que represente a la onda de la figura siguiente:
F igura 6.28
154
Análisis Matemático
III
arcsi arcsinnsin − | |
5.- Desarrollar en serie de Fourier las siguientes funciones periódicas:
02
6.- Desarrollar la función en serie de Fourier: a) de cosenos de arcos múltiples; b) de senos de arcos múltiples; c) en el intervalo Dibujar la gráfica de la función y las gráficas de las sumas de las series de Fourier para los casos a) y b). 7.- La temperatura de una persona es una función del instante del día y la curva de su temperatura normalmente se repite diariamente. Por ejemplo, su temperatura es la misma a las 2 p.m. de cada día. Los valores que siguen se extrajeron de una gráfica de temperatura tomada cada hora. Hallar una serie de Fourier hasta la tercera armónica para representar la curva de temperatura Tiempo
Media noche
2:00
4:00
6:00
8:00
Temp;°C
36·61
36·61
36·45
36·54
36·69
10:00
36·76
Media día
36·94
14:00
16:00
18:00
20:00
22:00
36·13
36·35
36·25
36·98
36·86
8.- La tabla muestra el desplazamiento x mm de una pieza deslizante con relación a un punto de referencia fijo para cada 30°de rotación de la manivela. Hallar una serie de Fourier pata x hasta la tercera armónica ɵ
x
0 298
30 356
60 373
90 337
120 254
150 155
180 80
210 51
240 60
270 93
300 147
330 221
UNIDAD VII 7. La tr ansformada de L aplace. Definición: Sea una función f definida en el intervalo 0, , se denota y define la
transformada de Laplace de f por:
L f (t )( s) e st f (t )dt 0
7.1 Transformadas de Laplace de algunas fun cion es básicas.
155
Análisis Matemático Lcos at ( s ) L sen at ( s )
s s 2 a 2 s s 2 a 2
III , s 0 , s 0
1
Le at ( s )
, s a s a n! Lt n ( s ) n1 , s 0, n s k Lk ( s ) ,.. s 0 s a ,.. s a L senh.at 2 s a 2 s ,.. s a Lcosh at 2 s a 2
Para el caso general de la función Lt ( s)
( 1)
s 1
, 1
, tenemos:
Donde la función Gamma ( ) se define por:
( x) e t t x 1dt y
tiene las siguientes propiedades:
0
( x 1) x ( x) ( x 1) n!, n N
1 2
7.2 Demostración de algunas de ellas.
1.- Hallar las transformadas de Laplace de las siguientes funciones: f (t ) t n
L t ( s) f ( s) e st t n dt lim n
0
b
e
b 0
st n
t dt , integrando por partes se tiene:
156
Análisis Matemático
III
dv e st u t 1 du nt n1dt v e st s n
b
st n
f ( s) lim e t dt b
t n e st b s
0
n
b
e s
st n 1
t dt
0 0 0 ∫ } 0 0
0
Para
Siguiendo el mismo procedimiento por algunas veces se llega a la siguiente respuesta:
Se tiene:
f (t ) senbt
L senbt ( s) f ( s) e st senbt .dt lim
b
e
st
senbt .dt , integrando por partes se tiene:
b 0
0
dv e st u senbt 1 st du b cos bt .dt v e s
b
st
f ( s) lim e senbt .dt b
senbt .e st b 0
s
0
b
b
e s
st
c cos bt .dt , volviendo a integrar por
0
partes y simplificando tenemos:
i 0 0 0 0 f (t ) e at
Le
at
( s) f ( s) e
b
b 0
e st e at dt lim
b
e
b 0
i − 0 0
e dt lim
st at
t ( s a )
dt , integrando se tiene:
f (t ) cosh at
157
Análisis Matemático
III
Lcosh at ( s) f ( s) e st cosh at .dt lim 0
b
e
st
b 0
cosh at .dt , integrando por partes se
tiene: dv e st u cosh at 1 st du a. senh.at .dt v e s
b
f ( s) lim e b
0
st
b cosh at .e st b a st cosh at .dt e senh.at .dt , volviendo a integrar 0 s 0 s
por partes y simplificando tenemos:
7.3
i − 0 0 0 ℎ || Principales propi edades de la transform ada de L aplace.
1. L y L-1 (operador inverso) son lineales. 2. L f ( n) (t ) s n L f (t ) s n1 f (0 ) ... f (n1) (0 ) 1 Ejemplo 1. Sea f (t)= senat , entonces f ' (t ) cos at , por la propiedad anterior: a L f ' (t ) sL f (t ) f (0)
Reemplazando:
( ) () 0
t 1 3. L f (t )dt L f (t ) 0 s
En términos de la transformada inversa:
158
Análisis Matemático
III
t
1 L ( s) L1 ( s)dt s 0 1
t
Ej emplo 1. Como
1 cosh at dt senh at tomando f (t) = cosh at , se tiene: a 0
t 1 L f (t )dt L f (t ) 0 s t 1 L senh at L cosh at dt a 0 1
1 L senh at Lcosh at a s as a L senh at 2 2 L senh at 2 2 s ( s a ) s a
4. Si L f (t ) ( s) , entonces Leat f (t ) ( s a) En términos de la transformada inversa: L1 ( s a) eat L1 ( s)
Ej emplo 1. Como Lcosh 4t
s s 2 16
, se tiene Le 2t cosh 4t
s 2 ( s 2) 2 16
2 s 3 Ej emplo 2. Calcular L1 2 s 4 s 20
Se tiene: 2 s 3 2 s 3 s 2 4 s 12 ( s 2) 2 16
s 2 7 4 2 2 2 ( s 2) 16 4 ( s 2) 16
Pero:
2t e cosh 4t ( s 2 ) 16 2t 4 L1 e senh.4t 2 ( s 2) 16 L1
s 2 2
159
Análisis Matemático
III
7 2t 2 s 3 2 t De donde: L1 2 2e cosh 4t e senh.4t 4 s 4 s 12
5. Si L f (t ) ( s) , entonces Lt f (t ) 1 n
n
d n ( s) ds n
En términos de la transformada inversa: d n ( s) n n 1 L (1) t L ( s) n ds 1
Ej emplo 1. Lt 2 senh.t Lt senh t 2
d 2 d d 1 L senh . t ( ) ds 2 ds ds s 2 1
16 s 2 2
s
2
1
3
1 2 s 2 1
Ej emplo 2. Calcular L1
Al comparar
1 ( s 2 1) 2
con
1 1 d 1 d 1 2 , se tiene que: 2 2 s ds s 2 1 ds s 1 ( s 1) 2
Así:
1 1 1 d 1 2 L ( s 1 ) 2 s ds s 1
L1
1
2
1
2
t
20
d 1 2 dt ds s 1
L1
t
1 1 tL 2 dt 2 0 s 1
1
− ∫ − 1
1
2
2
t cos t sen t
6. La función unitaria escalonada se denota y define por:
160
Análisis Matemático
III
1 si t a u a (t ) 0 si t a
Con la ayuda de esta función se tiene la siguiente propiedad: Lua (t ) g (t a ) e as L g (t )
En términos de la transformada inversa: L1e as ( s ) ua (t ) L-1 ( s)t t a
Ej emplo 1. Si f (t) = u a (t )sen t , entonces Lua (t )sen t Lua (t ) sen(t a) a
e as L sen(t a) eas L sen t cos a sen a cos t eas (cos a eas
1 s 2 1
sen.a
s
) s 2 1
cos a s. sen a s 2 1
e 5 s Ej emplo 2. Calcular L 2 s 10 s 26 1
Por la fórmula L1e as ( s) ua (t ) L1 ( s) t t a e 5 s 1 1 Tenemos: L1 2 u3 (t ) L 2 t t 5 s 10 s 26 s 10 s 26 1 u3 (t ) L1 t t 5 2 ( s 5) 1 u3 (t )(e 5t senh.t ) t t 5 u3 (t )e 5(t 5) senh(t 5)
7. La convolución de dos funciones f(t) y g(t), se denota y define por: t
f * g (t ) f ( ) g (t )d 0
Esta función tiene la siguiente propiedad: f * g = g * f
La transformada de Laplace de la convolución de dos funciones es: L f (t ) * g (t ) L f (t ) L g (t )
En términos de la transformada inversa: 161
Análisis Matemático
III
L1 ( s) ( s ) L1 ( s)* L1 ( s)
1 Ej emplo 1. Hallar L1 2 2 s ( s 4)
Tenemos: 1 1 1 1 L1 2 2 L 2 2 s s 4 s ( s 4) 1 1 L1 2 * L1 2 s s 4 t * senh.2t t
t
1
senh.2tdt cosh 2t 0
1
1
2
2
cosh 2t .
2
0
7.4 Apl icaciones de las tr ansformadas de Laplace. 7.4.1 En la solu ción de ecuacion es y sistemas de ecuacion es dif er enciales.
La transformada de Laplace, es una alternativa para resolver problemas de ecuaciones diferenciales con valores iniciales o de frontera, mediante esta técnica éstas se transforman en ecuaciones algebraicas. Ej emplo 1. Resuelva el problema de valores iniciales: y // 4 y / 13 y 2t 3e2t cos 3t ; y (0) 0, y / (0) 1
a.- Se Toma transformadas de Laplace en cada lado de la ecuación diferencial: L y 4 L y 13 L y 2 Lt 3 Le 2t cos 3t
s L y sy0 y' (0) 4 sL y y0 13 L y 2
( s 2 4 s 13) L y 1
2 s 2
3( s 2) ( s 2) 2 9
L y
2 s
2
3 s 2
s 22 9
1
2
s 2 4 x 13 s 2 s 2 4 s 13
3 s 2
s
2
4 s 13
2
b.- Tomando la transformada de Laplace inversa en cada término tenemos:
I.- Como:
1 s 2 4 s 13
1
3
3 s 2 2 9
Entonces su transformada inversa es: 162
Análisis Matemático
III
1 1 2t L1 2 e sen3t 3 s 4 s 13
II.-
2
A este término antes de encontrar su transformada inversa, le
s ( s 4s 13) 2
2
descomponemos utilizando fracciones parciales de la siguiente manera: 2 s 2 s 2 4 s 13
A s
B
s 2
Cs D s 2 4s 13
− −
Resolviendo tenemos: De donde: 2 s 2 s 2 4 s 13
8 1 169 s
2 1
13 s 2
8
s 2 10 3 169 s 22 9 507 s 22 9
Entonces su transformada inversa es: 2 8 2 8 2t 10 2t t L1 2 2 e sen3t e sen3t s s 4 s 13 169 13 169 507
III.- Como:
3 s 2
s
2
4 s 13
2
1 d
3 2 2 ds s 2 9
Su transformada inversa es: 3 s 2 1 1 3 1 1 3 d d L L 2 2 2 2 2 ds 2 ds s 2 9 s 2 9 s 4 s 13
L1
t 2t 1 3 1tL1 e sen3t . 2 2 s 2 9 2
Combinando los resultados, se tiene: y
179 507
e 2t sen3t
8 169
e 2t cos 3t
1 2
te 2t sen3t
2 13
t
Ej emplo 2. Resolver la ecuación Diferencial de valores iniciales: y 3 y 2 y 4t 12et ; y (0) 6, y(0) 1
a.- Se Toma transformadas de Laplace en cada lado de la ecuación diferencial: 163
Análisis Matemático
III
L y 3 y 2 y L4t 12et L y 3 L y 12 y 4 Lt 12 Let
s L y sy0 y' (0) 3 sL y y0 2 L y 2
( s 2 3 s 2) L y 6 s 1 18 ( s 3 s 2) L y 2
4 s
2
12 s 1
4 s
2
12 s 1
( s 2 3 s 2) L y 6 s 19
6 s 4 13 s 3 7 s 2 4 s 4 s 2 ( s 1)
L y
4 s
2
12 s 1
6 s 4 13 s 3 7 s 2 4 s 4 s 2 ( s 1)( s 2 3 s 2)
6 s 4 13 s 3 7 s 2 4 s 4 6 s 4 13 s 3 7 s 2 4 s 4 L y 2 s 2 ( s 1)( s 2 3 s 2) s ( s 1)( s 2)( s 1)
Descomponiéndole en fracciones parciales tenemos: L y
A B C D E 2 s s ( s 1) ( s 2) ( s 1)
Resolviendo tenemos: De donde: L y
3 s
2 s 2
2 ( s 1)
2 2 −2
2 ( s 2)
3 ( s 1)
b.- Tomando la transformada de Laplace inversa en cada término tenemos: 3 2 2 2 3 y L1 2 s s ( s 1) ( s 2) ( s 1)
De donde la Solución es: y 3 2t 2et 2e 2t 3et .
Ej emplo 3. Resuelva el problema:
2−1 − 20 00
Solución: recordemos que Ltf (t )
d
L f (t ),
Si L[y]= φ(s) entonces
2 −2 −2 0
d
ds d
ds
L y' ' 2
d
ds
L y' 2 L y' 2 Ly 0
d s s sy0 y' 0 2 s s y0 2 s s 2 s ds ds
s
2
2
2 s ' s 4 s 4 s 0
164
Análisis Matemático
III
Resolviendo esta ecuación con el método del factor integrante: 1
L y s c
1
2 s s
1 1 1 1 1 y L c s s 2 s 2 s 22 s 2 s 22 1
1
La Solución es: y = c (-1 + t + e-2t + te-2t). Ej emplo 4. Resuelva el problema de valores iniciales:
0001 Donde : f(t) =
0 si 0 ≤ t ≤ π 1 si π ≤ t ≤ 2 π 0 si t ≥ 2 π
Solución: Tenemos que f(t) = uπ(t) - u2π(t) Tomando transformadas de Laplace a cada lado: L[y’’] + L[y] = L[u π(t)]- L[u2π(t)]
1 1 s 2 L y sy0 y' 0 L y e s e 2 s s s
2 s 1 1 e 2 s s 2 1 s 2 1 s s 1
L y e s
1
1 1 s 2 s 1 2 s s L y e s e s 2 e e 2 2 s s 1 s s 1 s 1
Tomando transformadas inversas : y = uπ(t) - uπ(t)cos(t-π) – u2π(t) + u2π(t)cos(t-2π) + sen t y = uπ(t) - uπ(t)cos t - u2π(t) + u2π(t) cos t + sen t Sen t, si t ≤ π y=
1 + cos t + sen t, si π ≤ t ≤ 2π 2 cos t + sen t, si t > π
Ej emplo 5. Resolver el siguiente sistema de ecuaciones diferenciales:
−−25−2 −
Con las condiciones x (0) = 0, x´(0) = 0, y(0) = 0 y´(0) = 0.
Solución :
1.- Tomando la transformada de la place a cada una de las ecuaciones del sistema tenemos :
165
Análisis Matemático
III
− −−25−2 − −5 −−00− 0 5 0 −0 − −2 20 − −2 −1− 5−
2.- Resolvemos el sistema algebraico en s, por el método de CRAMER.
=
− ()− − ()− − − 2 −() −() 1−2 2
Tomando transformada Inversa tenemos :
=
Tomando transformada Inversa tenemos :
Ej emplo 6. Resolver el siguiente sistema de ecuaciones diferenciales:
166
Análisis Matemático
III
15 15 2 −
Con las condiciones x (0) = 35, x´(0) = -48, y(0) = 27, y´(0) = -55.
Solución :
1.- Tomando la transformada de la place a cada una de las ecuaciones del sistema tenemos :
1 5 − 15 2 15 − −0 −0−0 −0 −00 15 2 − 5−21 2−1 2.- Resolvemos el sistema algebraico en s, por el método de CRAMER.
− − 2 −15 − =
Tomando transformada
Inversa tenemos :
=
Tomando
transformada
Inversa tenemos :
− 0 −0 2 167
Análisis Matemático
III
7.4.2 Ej ercici os de apli cación de las ecuacion es dif er enci ales. Ej emplo 1. Un cuerpo que pesa 3lb estira un muelle 6plg (1/2 pie). Dicho cuerpo se suelta en
t = 0 desde un punto que esta a 8plg (2/3 pie) bajo la posición de equilibrio, con una velocidad dirigida hacia arriba de 4/3 pie/sg. Determine la función x (t) que describe el movimiento libre resultante.
Solución: Por la ley de Hooke tenemos: 1 lb F kx , reemplazando valores 3 k ( ) k 6 . 2 pie
Además pasamos las unidades de peso a unidades de masa: m
W g
3 32
slug .
Estos valores reemplazando en la ecuación que rige el movimiento tenemos: d 2 x
6 d 2 x 0 2 64 x 0 . Esta es una ecuación lineal homogénea 3 dt 2 dt 32
Condiciones de valores iniciales: x (0) = 2 / 3, x´ (0) = -4 / 3. 1.- Se Toma transformadas de Laplace en cada lado de la ecuación diferencial: L x(t ) 64 x(t ) L0 s 2 L x(t ) sx0 x' (0) 64 L x(t ) 0 2 4 2 4 2 s 2 ( s 2 64) L x(t ) s 0 ( s 2 64) L x(t ) s L x . 2 3 3 3 3 3 s 64
L x(t )
2
s
3 s 2 64 s 2 64 2
s
2
. 3 s 2 64 4( s 2 64) Tomando la transformada inversa: 2 s 8 2 1 x(t ) L1 2 . x(t ) cos 8t sen8t 2 3 s 64 4( s 64) 3 6 L x(t )
8
Que es la Solución del ejemplo propuesto.
168
Análisis Matemático
III
Ej emplo 2. Un inductor de 3 Henrios esta en serie con una resistencia de 30 ohmios y una
0
f.e.m. de 150 sen 20t voltios, suponiendo que en t = 0, la corriente i = 0, hallar la corriente en cualquier tiempo .
F igura 7.1
Estos valores reemplazando en la ecuación que rige el movimiento tenemos: L
di(t ) dt
Ri(t ) E
di(t ) dt
R L
i (t )
E L
. Condiciones iniciales: t = 0, i (0) = 0
1.- Tomando transformadas de Laplace a cada lado de la ecuación diferencial: R di(t ) R E 150 sen20t i(t ) L sLi(t ) i(0) Li(t ) L L L dt L L
L
( s
R L
) Li (t )
3000 L( s 2 400)
Li (t )
la inversa:
3000 , Reemplazando valores y tomando R 2 L( s )( s 400) L
10t i(t ) 2e 2 sen20t 2 cos 20t . 2 ( s 10)( s 400)
i(t ) 1000 L1
1
Ej emplo 3. Una cadena de 6m de longitud se desliza desde una mesa sin rozamiento. Si el
movimiento comienza desde el momento en que cuelga 1 m de la cadena. ¿Cuánto tiempo tardará en deslizarse toda la cadena? La ecuación que rige este movimiento es: d 2 x 2
dt
k m
x 0. Condiciones iniciales: t = 0, x (0) = 1 y x´ (0) = 0
Haciendo k 1 = k / m
d 2 x dt 2
k 1 x 0.
1.- Tomando transformadas de Laplace en cada lado de la ecuación diferencial: 169
Análisis Matemático
III
L x(t ) k 1 x(t ) L0
s L x(t ) sx0 x' (0) k L x(t ) 0 2
1
( s 2 k 1 ) L x(t ) s 0 ( s 2 k 1 ) L x(t ) s L x(t )
s s 2 k 1
Tomando la transformada inversa: s x(t ) L1 2 . x(t ) cosh k 1 t s k 1 Que es la Solución de la ecuación, pero nos piden el tiempo, por lo que tendremos que despejar t en función de x.
ℎ n√ −1 √ 5 Reemplazando los valores: Ln (6+
.
y x = 6 tenemos:
Que es la respuesta.
Ej emplo 4. (Problema Tautócrono)
Un alambre tiene la forma de una curva en el plano xy con su extremo más bajo en el origen. Asumiendo que no hay fricción, encuentre la forma que debe tener esta curva para que una bolilla bajo la influencia de la gravedad se deslice desde el reposo hasta 0 en un tiempo especificado T, independientemente de su posición inicial.
Solución: Sea A(x,y) el punto inicial del movimiento, y sea B (η,ξ) un punto de la curva buscada. Como se asume que no hay rozamiento, entonces no hay pérdida de energía, por lo tanto: EA = EB Si σ denota la trayectoria, se tiene:
d Mgy = mg ξ + m 2 dt 1
2
De donde: 170
Análisis Matemático
III
2
d 2 g y dt
Por lo tanto: d dt
2 g y
Escogemos el signo negativo, puesto que si t aumenta, la trayectoria disminuye, es decir es una función decreciente, así: d dt
2 g y
Por otro lado σ es una función que depende de (η, ξ), si suponemos que η = η (ξ), entonces σ depende solo de ξ, es decir: σ = F (ξ), dσ = F’(ξ) dξ
Reemplazando: F ' d
dt
2 g y
Integrando: t
0
1
dt
F ( )( y ) 2 g
0
'
1 2
d
y
De donde: T
1
'
2 g
F ( y ) * ( y )
1 2
(el símbolo * representa la convolución)
Tomando la transformada de Laplace, y teniendo en cuenta que:
L y n
(n 1)
s
n 1
1 y 2
Se tiene:
2 g T
L F ' ( y) L y 2 ; 2
1
De donde: L F ' ( y )
1
2gT
1 2
1 2
L F ' ( y )
1
s 2
2 g T 1
s 2
171
Análisis Matemático
III
Tomando la transformada de Laplace inversa: F ' ( y )
2 g
Ty
1 2
Además como: d 2 dx 2 dy 2 2
2
d d dx Tenemos: F ' ( y ) 1 y dy dt dt 2
dx 2 g Por lo tanto: 1 2 T 2 y 1 dy
De donde: dx
k y
1 donde k
2 gT 2 2
Haciendo y ksen2 Se obtiene: x 2k
2
sen2 4
Si 2 , se tiene que la Solución es: x
k
y
k
2 2
sen 1 cos
7.5 Ej ercicios propuestos. Resolver las siguientes ecuaciones y sistemas de ecuaciones diferenciales
1. 2. 3. 4.
− −−21 1 00 05 2 2−1 01 0 2−2 − − −− 0 00 −01 0 1 0 2 ,
172
Análisis Matemático
5.
6. 7. 8. 9. 10. 11. 12. 13.
14. 15.
III
− 2−0 0 0 0 −1 0 −20 2 2 −− 00 12 2 −2−02 002 1 0−1 0 2 0 01 − 52 0 0 0 1 − −− 0 00 5 0 − −21 0 0 0 0 1 2 −0 1 0 0 0 0 −1 −2012
16. Un cuerpo que pesa 5lb estira un muelle 8plg. Dicho cuerpo se suelta en t = 0 desde un punto que esta a 10plg arriba de la posición de equilibrio, con una velocidad dirigida hacia abajo de 4/3 pie/sg. Determine la función x (t) que describe el movimiento libre resultante. 17. Un inductor de 5 Henrios esta en serie con una resistencia de 40 ohmios y una f.e.m. de 200 sen 15t voltios, suponiendo que en t = 0, la corriente i = 0, hallar la corriente en cualquier tiempo .
0
F igura 7.2
173
Análisis Matemático
III
18. Un péndulo de ½ pie (15.24 cm) de longitud se suelta con una velocidad de ½ radián/seg. Hacia la vertical, desde una posición determinada por la distancia de 1/5 rad respecto de la vertical. Hallar la ecuación del movimiento 19. Una masa de ½ slug está suspendida de un resorte cuya constante es 4 lb/pie. El sistema se pone en movimiento en un medio que opone una fuerza de amortiguación numéricamente igual a dos veces la velocidad instantánea. Encuentre la Solución estacionaria si se aplica al sistema, a partir de t = 0, una fuerza exterior f (t) = 20sen 2t. 20. Un circuito consta de una inductancia de 0.05 henrios, una resistencia de 5 ohmios y un condensador de 4104 faradios de capacidad. Si q = i= 0 para t = 0, hallar q e i en función de t cuando a) hay una f.e.m. constante = 110 voltios, b) hay un f.e.m. alterna = 200 cos 100t . Hallar las soluciones de régimen permanente en b.
8. BIBL I OGRAFI A. 1. Ecuaciones Diferenciales por: KREIDER - KULLER - OSTBERG. 2. Problemas de Ecuaciones Diferenciales Ordinarias por: A. KISELIOV – M. KROSNOV - G. MAKARENKO.
3. Ecuaciones Diferenciales por: DONALD - L. KREIDER. 4. Ecuaciones Diferenciales por: RALPH PALMER AGNEV. 5. Ecuaciones Diferenciales Aplicadas por: M.R. - SPIEGEL. 6. Ecuaciones Diferencia2es Elementales por: L.M. KELLS. 7. Ecuaciones Diferenciales y Problemas con valores en la frontera por: WILLIAM E.
174