INTRODUCCIN A LOS MODELOS DINMICOS

Tema 1
INTRODUCCIN A LOS MODELOS DINMICOS

Modelos matemticos y teoras

Un modelo constituye una representacin abstracta de un cierto aspecto de la realidad. En su estructura intervienen, por una parte, los elementos que caracterizan la realidad modelizada y, por otra parte, las relaciones existentes entre ellos. Un modelo matemtico es un tipo de modelo basado en la lgica matemtica, cuyos elementos son esencialmente variables y funciones, y las relaciones entre ellos vienen expresadas a travs de relaciones matemticas (ecuaciones, inecuaciones, operadores lgicos...) que se corresponden con las correspondientes relaciones del mundo real que modelizan (relaciones tecnolgicas, leyes fsicas, restricciones del mercado...).

 Importancia de los modelos matemticos en Economa:

La construccin de modelos revela, a veces, relaciones que no son evidentes a primera vista Una vez construido el modelo, es posible extraer de l propiedades y caractersticas de las relaciones que de otra forma permaneceran ocultas En aquellas situaciones econmicas del mundo real en las que no es posible experimentar con la realidad, ofrecen un marco terico para evaluar la toma de decisiones as como sus consecuencias

Modelos matemticos estticos vs modelos matemticos dinmicos

En un modelo esttico la variable tiempo no desempea un papel relevante. En un modelo dinmico, por el contrario, alguno/s de los elementos que intervienen en la modelizacin no permanecen invariables, sino que se consideran como funciones del tiempo, describiendo trayectorias temporales. El anlisis de un modelo dinmico tiene por objeto el estudio de la trayectoria temporal especfica de alguno/s de sus elementos
MODELOS DINMICOS DETERMINISTAS VS MODELOS DINMICOS ESTOCSTICOS
.

 Un modelo dinmico determinista es aquel en el que, tanto a los parmetros como a las variables temporales, se les asignan valores determinados con certeza absoluta En general existen pocos modelos deterministas en el campo de la Economa y las Finanzas, ya que en la mayor parte de los casos, las variables y parmetros involucrados en los modelos econmicos y financieros (tasas de inters, precios de activos, ....) son impredecibles. Habitualmente la modelizacin dinmica en modelos econmicofinancieros hace uso de modelos estocsticos. En un modelo estocstico, alguna variable (o parmetro) sigue un proceso estocstico, es decir, que los valores que toma a lo largo del tiempo no son determinados con certeza absoluta sino que siguen una distribucin de probabilidad. El estudio de los modelos dinmicos estocsticos (y sus aplicaciones econmico-financieras) constituye el contenido fundamental de la asignatura

Ejemplo: Modelo de capitalizacin compuesta

Consideremos un depsito financiero a 3 aos, con capital inicial C0 y tasa de inters anual del 6 %. Disea un modelo de capitalizacin compuesta considerando que la capitalizacin es:
(a) Anual Elementos del modelo
Variable tiempo t: variable discreta t {0, 1, 2, 3} Variable de estado C(t) que describe la evolucin del capital a lo largo del tiempo. Es funcin del tiempo y el estudio de su trayectoria temporal es el objetivo del modelo. t = 1 : incremento de tiempo transcurido entre dos valores de la variable t, es decir entre dos periodos Los modelos discretos suelen trabajar con valores de t equidistantes y, por tanto, con un incremento constante. n = 3; nmero de periodos. Se cumple n* t = intervalo temporal total.

t = 1
0 C(0) 1 C(1) 2 C(2) 3 aos C(3)

Relaciones
- Queremos establecer la relacin existente entre el capital en un instante t y el capital en el instante siguiente t + t. Aplicando la ley de capitalizacin compuesta se tiene que: C(t + t) = C(t) + C(t)*0.06 = C(t)*(1+0.06)

Resolucin del modelo

- Procediendo recursivamente obtenemos C(3): C(0) = C0 C(1) = C0*(1+0.06) C(2) = C(1)*(1+0.06) = C0*(1+0.06)*(1+0.06) = C0*(1+0.06)2 C(3) = C0*(1+0.06)3

(b) Mensual Elementos del modelo

Variable tiempo t: t {0, 1/12, 2/12, ..., 12/12, ..., 24/12, ..., 36/12} Variable de estado C(t) t = 1/12 (1 mes) n = 36

t = 1/12

0 C(0)

1/12 2/12 C(1/12) ....

3/12 ...

t C(t)

t + t ... C(t + t) ...

36/12=3 aos C(1)

Relaciones
- Como la tasa de inters es anual y el periodo de capitalizacin es mensual, debemos convertir la tasa anual en mensual. Para ello, sustituimos 0.06 por 0.06/12 = 0.06*(1/12) = 0.06*t C(t + t) = C(t) + C(t)*0.06* t

 Modelo general: Modelo dinmico discreto en diferencias finitas Capitalizacin compuesta C0 capital inicial r tasa de inters anual t expresado de modo que permita transformar la tasa de inters anual en la correspondiente a la duracin del periodo utilizado C(t t) tt C(t+ + t)= =C(t) C(t)+ +C(t) C(t)*r* *r* C(0) C(0)= =C C 00

 Modelo dinmico continuo

Supongamos que en el ejemplo anterior disminuimos la duracin del periodo y trabajamos con capitalizacin diaria. El modelo seguira siendo: C(t + t) = C(t) + C(t) *r* t C(0) = C0 pero t pasa a ser t=1/360 (transformando la tasa de inters anual en diaria). Cuanto menor sea t, menor ser el periodo de capitalizacin utilizado. Si hacemos que t 0, entonces r* t representa la tasa de inters instantnea. Para obtener el modelo de capitalizacin instantnea o modelo continuo de capitalizacin procedemos como sigue:

 Operamos en la ecuacin de relaciones del modelo reagrupando los trminos: C(t + t) - C(t) = C(t)*r* t
C (t + t ) C (t ) = C (t ) * r t

Tomamos lmites cuando t 0 en ambos miembros de la igualdad:

lim C (t + t ) C (t ) = lim C ( t ) * r t 0 t
No depende de t

t 0

Es la definicin de derivada de C(t): C(t)

C(t) C(t)= =C(t)*r C(t)*r C(0)=C 0 C(0)=C

0

Modelo continuo de capitalizacin

Obtenemos una ecuacin en la que se relaciona una funcin con su primera derivada. Este tipo de ecuaciones sern estudiadas en el tema siguiente y constituyen la esencia de las modelizaciones dinmicas en tiempo continuo. Son las denominadas ecuaciones diferenciales.

Tema 2
ECUACIONES DIFERENCIALES Y SISTEMAS

Definiciones bsicas
Una ecuacin diferencial es una ecuacin que relaciona una funcin desconocida, tanto en una como en varias variables, con sus derivadas (o derivadas parciales) hasta un cierto orden. Distinguimos dos tipos de ecuaciones diferenciales:
Ecuaciones diferenciales ordinarias (EDO) : son aquellas en las que la funcin incognita es una funcin de una sola variable y = y(t). Se representan como: F(t, y, y, y, ... , y(n)) = 0 o bien en su forma normal o explcita como: y(n) = H(t, y, y, y, ... y(n-1)) Ejemplos: C(t)=C(t)*r ; y = y+1

 Ecuaciones en derivadas parciales (EDP): son aquellas en las que la funcin incognita es una funcin de varias variables y = y(t1, t2, ..., tn), por lo que las derivadas que intervienen en la ecuacin son derivadas parciales. Se representan como:

y y n y , ..., , ..., )= 0 F ( t1 , ..., t m , r1 r2 rm t1 tm t1 t 2 ... t m

Ejemplo: z = z(t,x,y)

z 2z z +3 2 +4 = z +1 t x y

A continuacin nos centraremos en el estudio de las ecuaciones diferenciales ordinarias EDO. Consideremos las siguientes definiciones bsicas:

 Solucin de una EDO: es una funcin y(t) que al ser sustituida en la ecuacin, la convierte en una identidad. Ejemplo: La solucin de la ecuacin y(t)= -ky(t) es y(t)=Ce-kt, siendo C una constante arbitraria. Solucin general de una EDO: es el conjunto de todas las soluciones de la ecuacin diferencial. Se expresa como una funcin que depende de uno o ms parmetros. Dando valores concretos a dichos parmetros se obtienen las denominadas soluciones particulares. Ejemplos: 1. Demostrar que la solucin general de la ecuacin y(t) = y(t) + t es y(t)=Cet-t-1 2. Demostrar que dos soluciones particulares de la ecuacin y(t) = y(t) + t son y(t)=et-t-1 ; y(t)=-t-1 3. Demostrar que y(t)=et-1 no es solucin de y(t) = y(t) + t

 Orden de una EDO: se llama orden de la ecuacin diferencial al mayor orden de derivacin que aparece en la ecuacin. Ejemplos: y(t)=t2+y2 es una EDO de orden 1; 3y(t)+2y(t) = et es una EDO de orden 3; Condiciones iniciales de una EDO: se llaman condiciones iniciales a los valores asignados a la funcin y a sus derivadas hasta orden una unidad menos del orden de la ecuacin, en un instante dado t0: y(n) = H(t, y, y, y, ... y(n-1)) y(t0)=a0 y(t0)=a1 ... y(n-1)(t0)=an-1

 Problema de valores iniciales: dada la solucin general de una EDO, se llama problema de valores iniciales al problema de encontrar una solucin particular para la cual se verifiquen las condiciones iniciales dadas. Ejemplo: Resolver el problema de valores iniciales para
dy ( t ) = y (t ) + t dt y (0 ) = 1

sabiendo que y(t)=Cet-t-1 es la solucin general

Ejemplos con Mathematica

Representacin de la solucin general de y(t) = y(t) + t ; representacin de la solucin particular y(t) = 2et-t-1 Representacin de la solucin general de y(t) = -2x2t

Mtodos de resolucin de EDO

Mtodos analticos: consisten en la obtencin de la funcin solucin de la EDO expresada por medio de una formula, es decir de una expresin analtica, que involucra funciones elementales. Ahora bien, en general, las EDO (y los problemas de valores iniciales) no tienen solucin analtica, por lo que es necesario recurrir a los denominados mtodos numricos para la obtencin de soluciones. Mtodos numricos: consisten en la obtencin de soluciones numricas de problemas de valores iniciales. La solucin se expresa entonces por medio vectores de valores numricos; es decir que, aunque no conozcamos la expresin analtica de la solucin, podemos determinar en cada instante del tiempo cual sera su valor numrico.

Mtodos analticos
Para poder aplicar un determinado mtodo analtico de resolucin, es necesario que la ecuacin presente una estructura particular. Si cambia dicha estructura el mtodo deja de ser valido. A continuacin recogemos una clasificacin de EDO que admiten resolucin analtica.
No lineales
Variables separables Homogneas y reducibles a homogneas Exactas y reducibles a exactas Bernouilli Riccati

EDO
Lineales

Orden 1 Orden n con coeficientes constantes

Homogneas Completas

VARIABLES SEPARABLES
Una ecuacin diferencial se dice de variables separables si es de la forma: dy (t ) = f1 (t ) * f 2 ( y ) ( f2(y) 0 ) dt Mtodo de resolucin:
Se agrupa dt con la funcin en t, f1(t), y dy con la funcin en y, f2(y)

dy (t ) = f1 (t ) * dt f2 ( y)
Se integran ambos miembros de la igualdad respecto de su variable correspondiente

dy = f1(t)dt f2 ( y)

 Ejemplo: dy =
dt

2y t ( y 2) y2 dt dy = 2y t

dt t = ln t + C2

y2 dy = 2y

dt t
1 y ln y + C1 = ln t + C2 2

1 1 1 y2 2 y dy = ( 2 y )dy = 2 y ln y + C1

1 y ln y = ln t + K Solucin general 2

LINEALES HOMOGNEAS A COEFICIENTES CONSTANTES Una ecuacin diferencial se dice que es lineal, de orden n, con coeficientes constantes y homognea si es de la forma:
y ( n ) ( t ) + a1 y ( n 1) ( t ) + a 2 y ( n 2 ) ( t ) + ... + a n 1 y( t ) + a n y ( t ) = 0

donde ai son nmeros reales e y(i) es la derivada i-sima de y(t). Mtodo de resolucin:
Se construye el polinomio caracterstico asociado a la ecuacin diferencial:

P(k) = k n + a1k n1 + a2k n2 ++ an1k + an

Se obtienen las raices del polinomio caracterstico:

k n + a1k n1 + a2 k n2 + + an1k + an = 0 Distinguimos los siguientes casos: Las raices son todas reales y distintas dos a dos: {k1, k2, ..., kn}
Sol. general Las raices son reales mltiples {k1 m1, k2 m2 ,, k p mp }

y(t) = C1ek1t +C2ek2t ++Cneknt (C1,..., Cn constantes arbitrarias)

Sol. general
1 k1t k1t m k1t p p 1 y(t) =C e C te t C e C e C te + + + + + + + +t 10 11 1(m 1) p0 p1 1 kt kt mp 1

Cp(mp1)e p

kt

Cij constantes arbitrarias Si una raiz es compleja simple + i (y su conjugada - i), entonces a la solucin general construida como en los casos anterior se le aadir el sumando:

C1et cost + C2et sen t

 Ejemplos: 1. y ( 2) (t ) + 3 y(t ) + 2 y (t ) = 0
Polinomio caracterstico P (k ) = k 2 + 3k + 2 k1= -1, k2 = -2 Raices

y(t ) = C1e t + C2 e 2t

Sol. general

2. y (3) (t ) 3 y (2) (t ) + 3 y(t ) y (t ) = 0

Polinomio caracterstico P ( k ) = k 3 3k 2 + 3k 1 k1= -1 m1 = 3 (raiz triple) Raices

y(t ) = C1et + C2tet + C3t 2et

Sol. general

3. y (2) (t ) + 4 y(t ) + 5 y (t ) = 0
Polinomio caracterstico Raices k1= -2 i
P (k ) = k 2 + 4k + 5
Sol. general

y (t ) = C1e2t cos t + C2e2t sen t

Mtodos numricos
Los mtodos numricos se aplican a problemas de valores iniciales y se basan en la discretizacin de la ecuacin diferencial considerada, es decir, en la transformacin de la ecuacin diferencial (en tiempo continuo) en una ecuacin en diferencias finitas (en tiempo discreto). Estudiaremos los esquemas de discretizacin de Euler y Adams, aplicados a un problema de valores iniciales con una EDO de orden 1 (con f(y(t),t) no lineal en el caso general):
dy (t ) = f ( y (t ), t ) dt y (t 0 ) = y 0

Si la EDO es de orden n, se transforma sta en un sistema equivalente de n ecuaciones diferenciales de orden 1 y se aplican los mtodos estudiados.

ESQUEMA DE DISCRETIZACIN DE EULER

Consideramos el problema genrico de valores iniciales:
dy (t ) = f ( y (t ), t ) dt y (t 0 ) = y 0

Mtodo de resolucin:
Suponemos que nuestro intervalo de trabajo es [t0,T], y consideramos que dicho intervalo est dividido en subintervalos [ti,ti+1], de manera que ti+1=ti + h, siendo h la longitud del subintervalo
h

t0

ti

ti+1

Aplicando a la funcin y(t) el desarrollo de Taylor de orden 1 obtenemos que: Donde o(h2) representa el resto del polinomio de Taylor que depende de h2. Si h es muy pequeo (infinitesimal), a o(h2) se le llama infinitsimo de orden 2, y puede considerarse despreciable a efectos de aproximar.

y ( t i + h ) = y ( t i ) + hy( t i ) + ( h 2 )

Sustituyendo en la ecuacin diferencial sabemos que:

dy(t ) y '(ti ) = = f ( y(ti ), ti ) dt ti

Uniendo ambas expresiones obtenemos:

y (ti + h) = y (ti ) + hf ( y (ti ), ti ) + (h 2 )

Despreciando el infinitsimo de orden 2, o(h2), tenemos que:

y (ti + h ) y (ti ) + hf ( y (ti ), ti )

Diagrama de flujo del esquema de discretizacin de Euler:
y(t0) = y0 i=0
y0 es conocido, y(t1) se calcula a partir de y0 y se convierte en un valor conocido a partir del cual calcular y(t2), y as sucesivamente

y(ti + h) y(ti ) + hf ( y(ti ), ti )

i = i+1

Dado el problema de valores iniciales:

dy (t ) = 1 + ( y t )2 dt y (0) = 0.5

Compara el comportamiento de su solucin analtica y(t)=t + 1/(2-t), con el comportamiento de la solucin numrica proporcionada por el esquema de discretizacin de Euler, en el intervalo [0, 1.9999]
Construimos un programa en lenguaje FORTRAN con la siguiente estructura:

PROGRAMA EULER PROGRAM EULER

Nombre del programa

*PARAMETROS

Y VARIABLES

IMPLICIT REAL*8(A-H,O-Z) PARAMETER(T0=0.,Y0=0.5,A=1.9999,N=1000) DIMENSION T(N),Y(N),Z(N)

Dividimos el intervalo en 1000 subintervalos. La variable T almacena los cortes de los subintervalos La variable Y almacena los valores de la aproximacin numrica La variable Z almacena los valores de la solucin analtica
*FICHEROS OPEN(20,FILE='EULER.SOL')

Guarda en el fichero Euler.sol los valores de T, Y y Z en tres columnas

*INICIALIZACION

DE MATRICES. DO I=1,N T(I)=0. Y(I)=0. Z(N)=0. END DO

Asigna a las variables un valor inicial de 0

*RUTINA DE CALCULO. H=A/N T(1)=T0+H Y(1)=Y0+H*(1+(Y0-T0)**2)

Calcula la longitud H de los subintervalos y el valor de la aproximacin de Euler en el primer corte Y(1).

DO I=2,N T(I)=T(I-1)+H Y(I)=Y(I-1)+H*(1+(Y(I-1)-T(I-1))**2) END DO

Calcula los valores de la discretizacin de Euler en el resto de cortes Evalua la expresin analtica de la solucin en cada corte del intervalo Transcribe los valores obtenidos al fichero de soluciones Euler.sol creado antes

DO I=1,N Z(I)= T(I)+1 / (2.0 -T(I) ) END DO

*ESCRITURA

EN UN FICHERO DE LOS VALORES T(K), Y(K) , Z(K)

DO K=1,N WRITE(20,*) END DO

STOP END

Instrucciones de finalizacin de programa

Anlisis de la solucin
El fichero Euler.sol que genera el programa anterior al ser ejecutado contiene tres columnas de nmeros de 20 dgitos: La 1 columna contiene los 1000 cortes del intervalo temporal de trabajo. La 2 columna los valores de la aproximacin numrica en cada corte. La 3 columna los correspondientes valores de la solucin analtica. A continuacin presentamos las 11 primeras filas del fichero para poder comparar la precisin de la aproximacin numrica realizada:
Euler.sol Aproximacin Valor analtico
0.501250062539648988 0.502500375429675961 0.503750939046151869 0.505001753765901285 0.506252819966504175 0.507504138026297791 0.508755708324378886 0.510007531240605383 0.511259607155598039 0.512511936450743111 0.513764519508193462

0.999949991703033472E-03 0.501249937489628761 0.199989998340606694E-02 0.502500125016744681 0.299984997510910041E-02 0.503750562956604031 0.399979996681213389E-02 0.505001251685214259 0.499974995851516779E-02 0.506252191579335986 0.599969995021820170E-02 0.507503383016484677 0.699964994192123560E-02 0.508754826374932745 0.799959993362426951E-02 0.510006522033711440 0.899954992532730341E-02 0.511258470372612406 0.999949991703033732E-02 0.512510671772190118 0.109994499087333712E-01 0.513763126613763332

CONCLUSIONES
El mtodo de discretizacin de Euler nos permite obtener una aproximacin numrica de la solucin de la ecuacin diferencial del ejemplo anterior, con una precisin de 4 cifras decimales. Existen esquemas de discretizacin ms sofisticados que proporcionan mucho mejores aproximaciones. Por ejemplo, el mtodo de Runge-Kutta aplicado al ejemplo anterior tiene una precisin de 15 cifras decimales. El mtodo de Euler utiliza el polinomio de Taylor de orden 1 para disear el esquema de discretizacin. Si incrementamos el orden del polinomio de Taylor, mejoramos la pecisin de la aproximacin. Este hecho es el que incorpora Runge-Kutta. Asimismo, el esquema de discretizacin de Adams que describimos a continuacin utiliza el desarrollo de Taylor de orden 2.

ESQUEMA DE DISCRETIZACIN DE ADAMS

Consideramos el problema genrico de valores iniciales:
dy (t ) = f ( y (t ), t ) dt y (t 0 ) = y 0

Mtodo de resolucin:
Suponemos que nuestro intervalo de trabajo es [t0,T], y consideramos que dicho intervalo est dividido en subintervalos [ti,ti+1], de manera que ti+1=ti + h, siendo h la longitud del subintervalo
h

t0

ti

ti+1

Aplicando a la funcin y(t) el desarrollo de Taylor de orden 1 obtenemos que:

h2 y (t i + h ) = y (t i ) + hy (t i ) + y(t i ) ( h 3 ) 2
Donde o(h3) representa el resto de orden 3 del polinomio de Taylor.

Sustituyendo en la ecuacin diferencial sabemos que:

dy(t) = f ( y(ti ), ti ) dt ti

d2 y(t) = f ( y(ti ), ti ) 2 dt t
i

Uniendo ambas expresiones obtenemos:

h2 y (ti + h ) = y (ti ) + hf ( y (ti ), ti ) + f (ti ) + ( h 3 ) 2

Despreciando el infinitsimo de orden 3, o(h3), tenemos que:

Utilizando la definicin de derivada, aproximamos:

h2 y (ti + h ) y (ti ) + hf ( y (ti ), ti ) + f ( y (ti ), ti ) 2

f ( y(ti ), ti ) f ( y(ti 1 ), ti 1 ) f '( y(ti ), ti ) h

Sustituyendo y operando en la expresin anterior obtenemos:

h y(ti + h) y(ti ) + (3 f ( y(ti ), ti ) f ( y(ti1 ), ti1 )) 2

Notemos que para calcular el valor de la funcin en el instante y(ti+h), necesitamos conocer y(ti) e y(ti-1). Por tanto, el algoritmo del esquema de discretizacin de Adams arranca calculando y(2) a partir de y(0) e y(1). El valor de y(0) es conocido, pero y(1) deber obtenerse a partir de y(0). Para ello utilizamos el esquema de discretizacin de Euler:

y (1) y (0) + hf ( y (0), 0)

Diagrama de flujo del esquema de discretizacin de Adams
y(t0) = y0 y(t1) = y0 + h f(y0 ,0) i=2

h y (ti + h) y (ti ) + (3 f ( y (ti ), ti ) f ( y (ti 1 ), ti 1 )) 2

i = i+1