Calculo Diferencial

martes, 26 de junio de 2012

4.2.1 Metodo de Euler

Consiste en dividir los intervalos que va de $x_o\,$ a $x_f\,$ en $n\,$ subintervalos de ancho $h\,$ ; osea:

$h = {x_f - x_o \over n}\,$

de manera que se obtiene un conjunto discreto de $n+1 \,$ puntos: $x_o, x_1, x_2,.......,x_n\,$ del intervalo de interes $[x_o,x_f]\,$ . Para cualquiera de estos puntos se cumlple que:

$x_i = {x_0 + ih}, \,$ $0 \le i \le n \,$ .

La condición inicial $y(x_o) = y_o \,$ , representa el punto $P_o = (x_o, y_o)\,$ por donde pasa la curva solución de la ecuación de el planteamiento inicial, la cual se denotará como $F(x)= y \,$ .

Ya teniendo el punto $P_o\,$ se puede evaluar la primera derivada de $F(x)\,$ en ese punto; por lo tanto:

$F'(x) = {dy\over dx} \bigg\vert\begin{matrix}\\{P_o}\end{matrix} = f(x_o,y_o)\,$

Grafica A.

Con esta información se traza una recta, aquella que pasa por $P_o\,$ y de pendiente $f(x_o, y_o)\,$ . Esta recta aproxima $F(x)\,$ en una vecinidad de $x_o \,$ . Tómese la recta como reemplazo de $F(x) \,$ y localícese en ella (la recta) el valor de y correspondiente a $x_1\,$ . Entonces, podemos deducir segun la Gráfica A:

${y_1 - y_o\over x_1 - x_o} = f(x_o,y_o) \,$

Se resuelve para $y_1\,$ :

$y_1 = y_o+(x_1 - x_o) f (x_o,y_o) = y_o + h f(x_o, y_o) \,$

Es evidente que la ordenada $y_1 \,$ calculada de esta manera no es igual a $F (x_1)\,$ , pues existe un pequeño error. Sin embargo, el valor $y_1 \,$ sirve para que se aproxime $F' (x) \,$ en el punto $P = (x_1,y_1)\,$ y repetir el procedimiento anterior a fin de generar la sucesión de aproximaciones siguiente:

$\begin{array}{crl} y_1 = y_o + h f (x_o,y_o)\\ y_2 = y_1 + h f (x_1,y_1)\\ .\\ .\\ .\\ y_{i+1} = y_i + h f (x_i,y_i)\\ .\\ .\\ .\\ y_n = y_{n-1} + h f (x_{n-1},y_{n-1})\\ \end{array} \quad \,$

4.2 Metodos de solucion para sistemas de Ecuaciones diferenciales lineales

Métodos de solución de sistemas de ecuaciones diferencialeslineales

Un sistema de diferenciales lineales puede resolver las ecuaciones. Al igual que existen varias técnicas para resolver una ecuación diferencial lineal, también las hay para un sistema de ecuaciones diferenciales lineales. Como el método de eliminación de Gauss, método separable y reducible etc. Sea un sistema de ecuaciones diferenciales lineales representado como,

Entonces, la representación de la matriz equivalente de este sistema de ecuaciones diferenciales lineales será,

Ahora, la técnica más conveniente para resolver este sistema de ecuaciones diferenciales lineales es recurrir al método de multiplicación de matrices. La fórmula para resolverlo está dada de la forma,

Aquí A es la matriz que contiene los términos coeficientes de toda la ecuación del sistema, C, que es una matriz columna compuesta de elementos no homogéneos y, finalmente, la matriz X es la que contiene los elementos desconocidos. Cuando la matriz C es igual a cero, entonces el sistema dado es un sistema homogéneo de ecuaciones diferenciales lineales.

Todo el mundo está familiarizado con el hecho de que multiplicar unas cuantas matrices es mucho mejor que solucionar las ecuaciones algebraicas crípticas para cuatro variables desconocidas. Sin embargo, la técnica anterior nos da la solución en todo momento. Por lo tanto, tenemos que adoptar algunas otras técnicas para resolver las ecuaciones.

Dos sistemas de ecuaciones diferenciales lineales se llaman equivalentes en la situación de que ambos produzcan las mismas soluciones para variables desconocidas. Sin embargo, esto puede lograrse mediante aplicar algunas operaciones elementales como la multiplicación con una constante,la reordenación de las ecuaciones, etc.Sea un sistema de ecuaciones diferenciales lineales dado como,

x + y + z = 0

x – 2y + 2z = 4

x + 2y – z = 2

Las ecuaciones anteriores pueden resolverse mediante aplicar algunas operaciones elementales sobre las dos últimas ecuaciones y manteniendo la primera sin alterar. Esto nos ayudará a eliminar una de las variables y el resto de las ecuaciones pueden resolverse de forma simultánea para dos variables, lo cual es muy fácil de hacer.

x + y + z = 0

x – 2y + 2z = 4

      y – 2z = 2

Ahora resta la dos ecuacionesa partir de la primera como,

x + y + z = 0

   - 3y + z = 4

     y – 2z = 2

Continuando con el procedimiento, multiplicamos la tercera ecuación con tres y lo sumamos a la segunda como

x + y + z = 0

   - 3y + z = 4

           - 5z = 10

Esto nos da el valor de z = −2. Al sustituir este valor en el sistema de ecuaciones podemos eliminar los términos que contienen z y finalmente,las ecuaciones han sidoreducidas a dos variablesúnicas. Los valores de las otras dos variables pueden obtenerse mediante la solución de esas ecuaciones simultáneamente.

Por lo tanto, el valor de x es 4, y el de y es −2.

Otra técnica popular es la transformar la matriz aumentada AC de manera talque se convierta en triangular superior. Una vez más esto puede lograrse mediante realizar las operaciones elementales de transformación en las ecuaciones del sistema. Después que esto se ha hecho, es posible resolver fácilmente el sistemapara las variables desconocidas.

4.1.3 Solución general de sistemas de Ecuaciones diferenciales lineales

Solución general de los sistemas de ecuaciones diferenciales lineales y solución particular de los sistemas de ecuaciones diferenciales lineales

Para resolver un sistema de ecuaciones diferenciales lineales, es absolutamente esencial conocer los conceptos de valor propio y vector propio. Para una matriz M dada, es llamadolos valores propios, si la condición es verdadera,

Aquí x se llama vector propio de la matriz M.

Es decir, los vectores propios son aquellos vectores que luego de ser multiplicados por la matriz de entrada permanecen proporcionales a la matriz de entrada o resultan cero.

Sea A la matriz que contiene los valores propios, como [ 1, 2, 3 … n]. A continuación se indican los pasos para resolver un sistema de ecuaciones diferenciales.

1. Calcula las ecuaciones del sistema de ecuaciones diferenciales y construye la matriz que contiene los coeficientes de todas las ecuaciones en el orden en que aparecen en el sistema de entrada de la ecuación.

2. Los valores propios se obtienen de esta matriz que contiene los términos de los coeficientes.

3. Calcula todos los vectores propios de valores propios como los obtenidos en el paso anterior. Nómbralos en la secuencia a medida que son determinados como EV1, EV2, EV3 …EVn.

4. Calcula las ecuaciones correspondientes para cada uno de los conjuntos de valores propios y los vectores propios asociados. Repite este paso para cada par de valores y vectores propios.

5. Obtén la solución particular para un sistema de ecuaciones no homogéneo como,

Aquí X(t) se define como,X(t) = [x1 x2 x3 … xn]

Y en caso de que el sistema de entrada sea una ecuación diferencial homogénea, entonces la solución particular del sistema será dada de la forma,

La ecuación anterior nos da la relación,

En la relación anterior, y son valores propios y vectores propios, respectivamente.

6. Y la solución general para el sistema de ecuaciones diferenciales lineales es dada de la forma,

En general podemos decir que la solución de un sistema de ecuación diferencial es llamada solución general si los valores de las constantes no se obtienen en la solución final. La misma solución puede convertirse en una solución particular cuando tenemos el valor de las constantes determinadas. Esto se hace en el caso que el sistema de entrada de la ecuación diferencial sea un problema de valor inicial con las condiciones iniciales establecidas para la determinación de los términos constantes.

El ejemplo siguiente aclarará el procedimiento para resolver un sistema de ecuaciones diferenciales lineales.

Determina el conjunto de ecuaciones como xT(t) = [(x1(t), x2(t)] para el sistema de ecuaciones dx/ dt = A * x con las condiciones iniciales establecidas como x(0) = x0 = (x01, x02). El valor de la matriz A estádada como,

La ecuación característica de la matriz de coeficientes arriba es,

f( ) = 2 –trace(A) * + det(A) = 2 + 4 + 4

Y las raíces de esta ecuación nos dan repetidos valores propios de la matriz como, 1 = 2 = −2.

Entonces, el vector propio de la matriz es dado de la forma,

La matriz tiene un solo vector propio ya que ambos valores propios son los mismos. Por lo tanto, la solución general del problema se da como,

Por consiguiente, un vector propio generalizado puede ser calculado como, v = vp +s1* vh1 v1 v2 −1 0 +s1* 1 1

La solución particular de este problema sería vT(p) = (−1, 0) = v2, el cual es el valor propio generalizado de esta matriz, junto con los valores propios repetidos y vh es la solución homogénea dando el vector propiov1.

Y la solución general del problema es,

Estonos da,

x1(0) = c1 e0 + c2 (0 – e0) = c1 – c2 = x01

x1(0) = c1 e0 + c2 (0) = c1 = x02

4.1.2 Sistemas de Ecuaciones diferenciales lineales homogeneos

Sistemas de ecuaciones diferenciales lineales homogéneos

Sabemos que una ecuación diferencial lineal es de la forma,

Si esta misma ecuación se transforma en la forma,

Obtenemos una ecuación diferencial lineal homogénea. Esta se da cuando lafunción conocida no estápresente en la ecuación diferencial lineal, entonces se le llama una ecuación diferencial homogénea. Y si tenemos una gran cantidad de tales ecuaciones juntas, de manera tal que dependen unas de las otras, y definen colectivamente un problema común, entonces se les llama un sistema de ecuaciones diferenciales linealeshomogéneo.

Tales sistemas pueden ser resueltos de manera eficiente con la ayuda de las matrices, las cuales son denominadas matriz fundamental. Sean X1, X2 … X3 las soluciones de la matriz fundamental del sistema de entrada de ecuaciones diferenciales homogéneas, entonces puede representarse de manera condensada como,

En la ecuación anterior, las soluciones del sistema de ecuaciones diferenciales están definidas en algún intervalo, digamos I y la solución general del sistema de ecuaciones diferenciales es este

En la ecuación anterior, los términos que se mantienen dentro de los corchetes son los vectores fila, donde X1 = [xi1j], X2 = [xi2j] …Xn = [xinj]. Estas son las soluciones n fundamentales del sistema de entrada de ecuaciones diferenciales lineales homogéneas para el intervalo dado I.Entonces tenemos que la matriz fundamental para el sistema homogéneo de ecuaciones diferenciales lineales para el intervalo dado como I es,

Los pasos para resolver un sistema homogéneo de ecuaciones diferenciales lineales son los siguientes:

1. Construye la matriz de coeficientes para las ecuaciones del sistema dado.

2. Determina los valores propios de esta matriz de coeficientes del sistema dado de ecuaciones diferenciales lineales homogéneas.

3. Ahora, busca el vector propio inicial de este conjunto de valores propios y nómbralo como EV1.

4. Determina la primera ecuación de este vector en términos de constantes k1, k2.

5. Después de esto, determina el siguiente conjunto de valores propios, sus vectores propios correspondientes y su ecuación.

6. Anota la solución general para las ecuaciones en términos de constantes k1, k2.

7. Por último, deriva la solución general para el sistema de ecuaciones.

Aunque el procedimiento para solucionar el sistema de ecuaciones diferenciales lineales homogéneo es bastante fácil, se da un ejemplo ilustrativo que te ayudará a hacer los conceptos más claros.

dx/ dt = 2x + 3y

dy/ dt = 2x + y

Primero, escribamos la matriz constante para el sistema de ecuaciones diferenciales homogéneas dado. Esto es,