Ecuaciones Diferenciales Y Transformadas DE Laplace CON Aplicaciones PDF

Preview

Summary

te da toda la informacion ademas de ejercicios de las ecuaciones diferenciales y transformadas de Laplace con sus aplicaciones...

Description

Luis M. Sánchez Ruiz Matilde P. Legua Fernández

ECUACIONES DIFERENCIALES Y TRANSFORMADAS DE LAPLACE CON APLICACIONES

EDITORIAL UNIVERSITAT POLITÈCNICA DE VALÈNCIA

Para referenciar esta publicación utilice la siguiente cita: Sánchez Ruiz, Luis M; Legua Fernández, Matilde P. (2017). Ecuaciones diferenciales y transformadas de Laplace con aplicaciones. 5 HG. 9alencia: Editorial Universitat Politècnica de València

© Luis M. Sánchez Ruiz Matilde P. Legua Fernández

© 2017, Editorial Universitat Politècnica de València distribución: www.lalibreria.upv.es / Ref.: 0798_01_05_01 Imprime: Byprint Percom, sl

ISBN: 978-84-9048-649-8 Impreso bajo demanda

La Editorial UPV autoriza la reproducción, traducción y difusión parcial de la presente publicación con fines científicos, educativos y de investigación que no sean comerciales ni de lucro, siempre que se identifique y se reconozca debidamente a la Editorial UPV, la publicación y los autores. La autorización para reproducir, difundir o traducir el presente estudio, o compilar o crear obras derivadas del mismo en cualquier forma, con fines comerciales/lucrativos o sin ánimo de lucro, deberá solicitarse por escrito al correo [email protected]. Impreso en España

Índice General 1 Ecuaciones diferenciales

1

1.1 1.2

Introducción y definiciones básicas . . . . . . . . . . . . Ecuaciones de primer orden . . . . . . . . . . . . . . . .

1 2

1.3

Ecuaciones de variables separables . . . . . . . . . . . . . 1.3.1 Definición y resolución . . . . . . . . . . . . . . .

5 5

1.3.2

Ecuaciones diferenciales homogéneas . . . . . . .

7

1.3.3

Ecuaciones diferenciales reducibles . . . . . . . .

9

1.4

Ecuaciones diferenciales exactas.............................................11 1.4.1 Definición y resolución................................................11 1.4.2 Factores integrantes..................................................13

1.5

Ecuaciones lineales..................................................................14

1.6 1.7

1.5.1

Definición y resolución................................................14

1.5.2

Ecuación de Bernoulli..................................................16

1.5.3 Ecuación de Riccati.....................................................17 Aplicación: Problemas de mezclas..........................................18 Ecuaciones no lineales en y 0.................................................... 19 1.7.1 Ecuaciones resolubles en y0.............................................. 19 1.7.2

Aplicación: Cálculo de la envolvente..........................19

Ecuaciones resolubles en y..........................................21 1.7.4 Ecuaciones resolubles en x.......................................24 Aplicación: Trayectorias isogonales........................................26 Ejercicios resueltos..................................................................27 1.7.3

1.8 1.9

1.10 Ejercicios propuestos...............................................................33 i

ii 2 Ecuaciones de orden superior 39 2.1 Introducción.............................................................................39 2.2 Ecuaciones diferenciales incompletas......................................40 2.2.1 2.2.2 2.3

Ecuaciones donde falta la y..........................................40 Ecuaciones donde falta la x..........................................42

Ecuaciones lineales de orden n................................................43 2.3.1 Definición. El operador derivada.................................43 2.3.2 Ecuación lineal homogénea..........................................45 2.3.3 Ecuación lineal no homogénea.....................................49 2.3.4 Método de Lagrange o variación de parámetros . 54

2.4

2.3.5 Ecuación de Euler-Cauchy...........................................56 2.3.6 Reducción del orden.....................................................59 Aplicación: Circuitos eléctricos...............................................60

2.5 2.6

Ejercicios resueltos..................................................................61 Ejercicios propuestos...............................................................63

3 Sistemas de ecuaciones 67 3.1 Introducción.............................................................................67 3.2

3.3 3.4

Método matricial...................................................................69 3.2.1 Sistema lineal homogéneo............................................69 3.2.2 Sistema lineal no homogéneo.......................................77 Método de eliminación.............................................................79 Aplicación: Transformadores y redes......................................82

3.5 3.6

Ejercicios resueltos..................................................................84 Ejercicios propuestos...............................................................93

4 Métodos numéricos 95 4.1 Interpolación............................................................................95 4.2 4.3

4.4

Resolución numérica de PVI....................................................98 Métodos de un paso.................................................................99 4.3.1 Método de Euler........................................................99 4.3.2 Métodos de Runge-Kutta.........................................102 Métodos multipaso lineales....................................................103

iii

4.5 4.6

4.4.1 Métodos explícitos e implícitos..................................103 4.4.2 Generación de métodos lineales.................................104 4.4.3 Métodos predictor-corrector.....................................106 Ejercicios resueltos.................................................................107 Ejercicios propuestos..............................................................109

5 Transformadas de Laplace 111 5.1 Definición y conceptos básicos..............................................111 5.2 Propiedades y transformadas...............................................116 5.3 Aplicación: Integrales impropias...........................................130 5.4

Transformada inversa de Laplace........................................130 5.4.1 Definición y propiedades............................................130 5.4.2 Cálculo de algunas transformadas inversas..............135 5.4.3 Método de Heaviside..................................................138

5.5

5.6

Aplicaciones de las transformadas.........................................142 5.5.1 Resolución de PVI. Función de transferencia..........142 5.5.2 Ecuaciones con coeficientes variables........................147 5.5.3 Ecuaciones integro-diferenciales.................................148 5.5.4 Sistemas de ecuaciones diferenciales..........................148 5.5.5 Estabilidad de sistemas dinámicos.............................149 Ejercicios resueltos.................................................................154

5.7

Ejercicios propuestos..............................................................166

Prólogo Las ecuaciones diferenciales modelan casi todos los procesos que aparecen en la técnica en los cuales hay una relación de cambio entre las variables involucradas. En dichos procesos es habitual contar con unas condiciones iniciales de partida y, en la mayoría de ocasiones, se pueden emplear entonces dos herramientas específicas para la búsqueda de solución: por métodos numéricos cuya utilización aconseja conocer algunas nociones de interpolación y mediante la transformada de Laplace que permite obtener soluciones exactas en los casos más usuales como son por ejemplo los modelados por ecuaciones o sistemas de ecuaciones lineales. Con esta publicación se pretende satisfacer las necesidades básicas que puedan tener los alumnos de ingeniería en técnicas de resolución de ecuaciones diferenciales y problemas de valor inicial. Se han incluido los resultados teóricos que dan el soporte matemático necesario para abordar los problemas que aparecen, pero se ha evitado dar demostraciones que sean excesivamente tediosas o complicadas de desarrollar; solo se han in- cluido aquellas que pueden ayudar a conseguir una formació n adecuada para abordar tipos de ecuaciones o aplicaciones no tratadas en este texto. A lo largo del mismo hay una amplia exposición de ejemplos que facilitan la comprensión de los diferentes temas presentados, finalizando cada ca- pítulo con una selección de ejercicios cuya resolución se recomienda para verificar que se ha entendido la materia desarrollada. Los autores expresan su reconocimiento al profesor Manuel Legua (1924—99), catedrático desde 1964 a 1989 de la Escuela Universitaria de Ingeniería Técnica Industrial de Valencia —transformada en Escuela Técnica Superior de Ingeniería del Diseño en 2002—, que les transmitió la forma de enfocar la didáctica de las Matemáticas destinadas a cubrir las necesidades de los ingenieros. Asimismo desean expresar su agradecimiento a Isabel Morales por su ayuda en la elaboración de la primera versión de este texto, a los compañeros de profesión que han hecho sugerencias respecto de dicha y posteriores versiones, y a los alumnos que, con sus dudas y querer saber, les han hecho ver los temas en los cuales tenían una mayor dificultad. Esperamos que estas notas faciliten la labor de nuestros futuros alumnos. Los autores

Capítulo 1 Ecuaciones diferenciales 1.1

Introducción y definiciones básicas

Existen situaciones en que se desea determinar una función desconocida a partir de una ecuación, denominada diferencial, (ED) que contiene por lo menos una de sus derivadas respecto de una variable independiente. Ejemplo 1.1.1 Se lanza verticalmente una partícula P de masa m des- de la super cie de un planeta sin atmósfera cuya forma es una esfera de radio R. Conociendo el valor de la gravedad g sobre su super cie, establecer la velocidad v de P en función de su distancia x al centro del planeta. Hallar la velocidad que debemos imprimir a P para que escape del campo gravitacional si g = 9.81 m/s 2, R = 6500 km. Sol.: La fuerza gravitacional es F (x) = km. Como F (R) = km = —mg x R resulta que k = —gR2. Como es habitual denotamos a = dvd , entonces ma = — Integrando,

1

gR m 2 x

= md dv = md dv ddx = mvddv =, v dv = —x 2 dx.

v2 =

gR

2

gR

2 + C. Si dotamos a P con una velocidad inicial

x

v (R) = v0, entonces C = 1 v2 2 0— gR por lo que v2 =

2gR

x

2gR

2 + v2 — 2gR , v = q 0

2 + v2 — 2gR. x

0

Si v20 ≥ 2gR, P escapa ya que v no se anula. Para g = 9.81 m/s2 y R = 6500 km, la velocidad de escape es v 0 ≥ 11.293 km/s. 1

2

Capítulo 1

Una ecuación diferencial es ordinaria (EDO) si la incógnita es función de una sola variable, y en derivadas parciales (EDP) si es función de dos o más. En lo que sigue nos ocuparemos únicamente de las primeras. Una ecuación diferencial ordinaria es de orden n si involucra hasta la derivada n-sima de la función desconocida, diciéndose que está expresada en forma normal si viene dada por y(n) = f(x, y, y 0,. . ., y (n31)). Una función definida en un intervalo I ⊆R es solución o integral de la ecuación diferencial si al sustituirla en ella se obtiene una identi- dad. Si la solución contiene n constantes arbitrarias se llama integral general (IG). Cada solución obtenida dando valores a las constantes arbitrarias se llama integral particular (IP) . Y las soluciones que no pueden obtenerse a partir de una integral general se llaman integrales singulares (IS). Las gráficas de las soluciones se denominan curvas integrales. Se resuelve un problema de valor inicial (PVI) si, dados x 0, y0, y00 , . . . , y0(n31) 5 R, se busca una solución y(x) que satisfaga y(x 0) = y0 , y0 (x0) = y0 , . . . , y (n31) (x0) = y (n31) . 0

1.2

0

Ecuaciones de primer orden

Las EDO de primer orden se generan si entre la ecuación que representa una F familia de curvas F (x, y, C) = 0 y su derivada total respecto de x, F x(x, y, C )+ F y(x, y, C) y0 = 0, eliminamos el parámetro C. El resultado se denomina ecuación diferencial de F. Ejemplo 1.2.1 Hallar la ecuación diferencial de (x — c)2 + y2 = 1. Sol.: Eliminando c entre la ecuación de esta familia de circunferencias y su d erivada³ respecto´de x, 2 (x — c ) + 2 yy 0 = 0, obtenemos su diferencial ecuación y 2 y0 2 + 1 = 1. Aquí (x — c)2 + y2 = 1 es la integral general y cada una de las circunferencias es una integral particular. Es fácil ver que y = ±1 son soluciones singulares. Hay ecuaciones diferenciales, como 2 + y 02 = 0, que carecen de solución. Enunciaremos dos resultados que garantizan la existencia de soluciones, única en el primero de ellos con ayuda del siguiente concepto.

3

Ecuaciones diferenciales

2 Definición 1.2.2 Una función f : A⊆R $R veri ca una condición de Lipschitz con respecto a y en A si hay una constante L, denominada de Lipschitz, tal que |f ( x, y1) — f (x, y2)| ≤ L | y1 — y2 | en A.

Si f y es continua en A = [a, b] × [c, d], f satisface una condición de Lipschitz ya que eRntonces | fy ( x, y)| ≤ K ; (x, y) 5 A, por lo que y |f (x, y1) — f (x, y2)| ≤ y2 |fy (x, y)| dy ≤ K | y1 — y2| , ( x, y1) , ( x, y2) 5 A. 1

Teorema 1.2.3 (de Picard) Sea f una función continua en el rectángulo A = [a, b]×[c, d] donde satisface una condición de Lipschitz con respecto a y. Dado (x 0, y 0) del interior de A, existe un h > 0 tal que y0 = f (x, y), y(x0 ) = y 0, tiene una única solución en [x0 — h, x0 + h]. Teorema 1.2.4 (de Peano) Sea f una función continua en el rectán- gulo A×= [a, b] [c, d]. Si (x0, y0) está en el interior de A, entonces existe un h > 0 tal que y0 = f(x, y), y(x0) = y 0, tiene al menos una solución en [x0 — h, x 0 + h]. 1

Ejemplo 1.2.5 Analizar si f (x, y) = 2y 2 satisface una condición de 2 : |x| ≤ 2, |y| ≤ 1}. Estudiar Lipschitz en el rectángulo A = {(x, y) 5 R 1 0 la unicidad de solución del PVI y = 2y 2 , y(0) = 0. Sol.:

Como

f (0,y) 3f(0,0)

2y

=

32 1

no está acotado en las cercanías del

y30

origen, f no satisface una condición de Lipschitz en A. Por tanto no es de extrañar que, a pesar de ser f continua en A, el PVI dado presente más de una solución como lo son y 1 (x) = x2 , y2 (x) = 0. Resolución gráfica. Isoclinas Una ecuación diferencial y 0 = f (x, y), donde f es una función continua en D ⊆ R2, asocia a cada P (x, y) 5 D una dirección de coeficiente angular m = f (x, y), que es la dirección en P de cualquier solución que pase por P . Los puntos de D con su dirección se llama campo de direcciones, pudiéndose visualizar con ayuda de pequeños segmentos con punto medio en ciertos P 5 D y su correspondiente dirección. Si distribuimos los puntos P sobre curvas de ecuación m = f(x, y), m5 R, se obtienen las curvas isoclinas que contienen a los puntos a los que corresponde la misma dirección.

4

Capítulo 1

Ejemplo 1.2.6 Emplear el método de las isoclinas para hallar: a) La curva que pasa por— ( 1, 2) y en cada punto (x, y) la pendiente de la tangente es el cuadrado de la abscisa x. b) Las curvas tales que en cada punto (x, y) la pendiente sea el cociente entre ordenada y abscisa. c) La curva que pasa por el punto (0, 1) y veri ca que la pendiente de la tangente en cada punto (x, y) es igual al producto de ordenada y abscisa.

Sol.: a) La propiedad enunciada es y0 = x2. Las isoclinas son s s m = x 2 , x = m, x = — m, m ≥ 0. El campo de direcciones en un entorno de (—1, —2), por ejemplo en el cuadrado [—2, 0] × [—3, —1] , es el siguiente.

Este método no proporciona la expresión analítica de la curva integral que pasa por (—1, —2) pero sugiere su gráfica. b) La propiedad enunciada es y 0 = yx. Las isoclinas son m = yx con m 5 R, y representan los puntos de curvas integrales en los que la pendiente es m. En este caso coinciden con las curvas integrales. c) La propiedad enunciada es y 0 = xy por lo que las isoclinas son las hipérbolas m = xy. Dando a m valores negativos, positivos y 0, notamos que las curvas integrales son decrecientes en el segundo y cuarto cua- drantes, crecientes en el primero y tercero, tienen pendiente nula sobre OY (recta x = 0) y la recta y = 0 también es solución de la ecuación

5

Ecuaciones diferenciales

diferencial. La gráfica de la curva integral que pasa por (0, 1) es 5

4

3 y2

1 0 -2

1x

2

-1

.

1.3

Ecuaciones de variables separables

1.3.1

Definición y resolución

Una EDO y0 = f(x, y) se dice que es de variables separables si existen dos funciones continuas P y Q tales que, operando, puede escribirse como P (x) dx = Q(y) dy. Entonces se dice que tiene las variables separadas y una integral general es Z P (x) dx = Z Q(y) dy + C, C constante arbitraria.

Ejemplo 1.3.1 Resolver a)

(1 — y2) — 3xy y0 = 0.

b)

(1 — y 2) — 3xy y0 = 0,

c) d)

(1 — y 2) — 3xy y0 = 0, y( — I33 ) = 2. (1 — y2) — 3xy y0 = 0, y(4) = —1.

e)

(1 — y2) — 3xy y0 = 0,

y(4) = — I23 .

y(0) = 2.

Sol.: a) Separando las variables x e y, dx 3y dy = . x 1 — y2 Integrando en ambos miembros, 3

ln |x| = —23 ln |1 — y2| + ln C , |x| |1 — y2| 2 = C, C > 0. (IG)

6

Capítulo 1

Las soluciones de una EDO son funciones definidas en un intervalo donde la satisfacen. En este caso podemos expresar las soluciones en forma explícita p p y = ±q1 — 3 ( C/x)2, x > C, y = ±q 1 — 3 (C/x)2 , x < —C, correspondientes a soluciones en las que 1 — y2 > 0, e p p y = ±q 1 + 3 ( C/x )2, x > 0 , y = ± q1 + 3 (C/x) 2, x < 0, — y2 < 0. Por otra parte siemcorrespondientes a soluciones en las que 1 pre que separamos variables hemos de esudiar si hemos perdido alguna solución. En este caso hemos dividido por 1—y 2 por lo que hemos descartado que dicha expresión sea nula, Y efectivamente en este caso a) también admite las soluciones singulares y = ±1. (IS) ¯

I

³

I

´¯

3

=

—2

2

b) Para que y(4) = — —

2

3

= C. Por tanto la

2 2 , 4 ¯1 ¯¯ 3 ¯ 2 |x| ¯1 — y2 ¯ 2 = 1. 3

1

Analizando la condición inicial, que buscamos una solución en 4 > 0 donde y(4) < 0 y además 1 — y2 = 1 — 43 > 0, observamos que la función 1 que la satisface es (1/2x)2, x > . y = — q1 — 2 p I 3 , 3 |

3 c) Para que y (— 3 ) = 2 del PVI viene incluida en I

3

1 — 22| 32 = 3 = C. Por tanto la solución 3

|x| ¯1 — y2¯ 2 = 3. Analizando la condición inicial, que buscamos una solución en — I

I3

3

0 y además 1 — y = 1 — 4 < 0, observamos que la función que la satisface es p y = q 1 + 3 (3/x ) 2, x < 0.

7

Ecuaciones diferenciales 3

d) Si en la IG se plantea y(4) = —1 , 4 (1 — 1) 2 = 0 = C conduce a ¢ ¡ x 1 — y2 2 3 = 0 que no debe considerarse como solución ya que el proceso de obtención de la IG es válido para C > 0. La función y = —1 es la solución de es...