DINÁMICA DISCRETA

DINÁMICA DISCRETA

Henri Claude Thonon (CV)

Volver al índice

C.– MODELOS DE ORDEN MAYOR QUE 2.

1.–     El Modelo del ciclo económico de Hicks.

          Este modelo es una generalización del modelo B.2.

El modelo:
1.–     Se supone que la renta nacional (Y) consta de tres corrientes: Consumo (C), Inversión (I) y Gasto Público (G):
                             Yt = Ct + It + Gt
2.–     El consumo no se define como una función de la renta actual, sino de la renta de los periodos anteriores; y para simplificar se supone de que este consumo es proporcional a dicha renta:
                             Ct = C(Yt–1, Yt–2, …..,Yt–n)
                             Ct = γ1•Y t–1 + γ2•Y t–2 +…..+ γn•Y t–n            (0 < åγi < 1)
3.–     La inversión se mantiene en una relación fija con respecto a los incremento de la renta de n periodos anteriores:
                    It = k1•∆Y t–2+ k2•∆Y t–3 + ........+ kn•∆Y t–n–1          (k1 ,k2,…, kn > 0)
                    It = k1•(Yt–1 − Y t–2) + k2•(Yt–2 − Y t–3)+ ........+ kn•(Yt–n − Y t–n–1)  (k1 ,k2,…, kn > 0)

4.–     El gasto Público se supone que aumenta en el tiempo a la tasa constante g:
                             Gt = Go•(1 + g)t

          De esta manera el modelo queda representado con el siguiente sistema de ecuaciones:
                    Yt = Ct + It + Gt
                    Ct = γ1•Y t–1 + γ2•Y t–2 +…..+ γn•Y t–n            (0 < åγi < 1)
                    It = k1•(Yt–1 − Y t–2) + k2•(Yt–2 − Y t–3)+ ........+ kn•(Yt–n − Y t–n–1)  (k1, k2,…, kn > 0)
                    Gt = Go•(1 + g)t

Solución del modelo (para el caso n = 2):

                    En este caso se tiene el sistema de ecuaciones:
                    Yt = Ct + It + Gt
                    Ct = γ1•Y t–1 + γ2•Y t–2                             (0 < γ1 + γ2 < 1)
                    It = k1•(Yt–1 − Y t–2) + k2•(Yt–2 − Y t–3)                             (k1, k2  > 0)
                    Gt = Go•(1 + g)t

          Sustituyendo las tres últimas ecuaciones en la primera:
Yt = γ1•Y t–1 + γ2•Y t–2 + k1•(Yt–1 − Y t–2) + k2•(Yt–2 − Y t–3) + Go•(1 + g)t
reagrupando términos y normalizando:
          Yt+3 – (γ1 + k1)•Yt+2  – (γ2 + k2 − k1)•Yt+1 + k2•Yt = Go•(1 + g)t+3.
Por lo tanto la solución vendrá dada por las raices de la ecuación cúbica:
x3 – (γ1 + k1)• x2  – (γ2 + k2 − k1)•x + k2 = 0
y la solución particular por:
Yp = Go•(1 + g)t+3 / [(1 + g)3 – (γ1 + k1)•(1 + g)2 – (γ2 + k2 − k1)•(1 + g) + k2]

          Análisis de la solución:
                    Dado que el término independiente (k2) es mayor que cero, y el coeficiente del término cuadrático (– (γ1 + k1)) es negativo tenemos los siguientes casos:
          a)       Dos raices reales positivas y una negativa, todas de módulo menor que 1. En este caso la solución converge a la solución particular:

γ1 =

0,7

 

m1=

0,62072308

 

γ2 =

0,01

 

m2=

0,27736095

0

I

Go =

4

 

m3=

-0,05808404

0

I

Co =

2

 

 

 

 

Io =

2

 

 

 

 

g =

0,2

 

 

 

 

k1 =

0,14

 

 

 

 

k2 =

0,01

 

 

 

 

 

 

 

 

t

Y

C

I

G

Yp

0

8

2

2

4

10,2795955

1

11,9

5,6

1,5

4,8

12,3355146

2

14,716

8,41

0,546

5,76

14,8026175

3

17,76544

10,4202

0,43324

6,912

17,763141

4

21,3324496

12,582968

0,4550816

8,2944

21,3157692

5

25,5935249

15,1103691

0,52987574

9,95328

25,578923

6

30,7049485

18,1287919

0,63222063

11,943936

30,6947076

7

36,8403325

21,7493992

0,75821007

14,3327232

36,8336492

8

44,2046181

26,0952822

0,91006799

17,1992678

44,200379

9

53,0431112

31,311636

1,09235382

20,6391214

53,0404548

10

63,6502016

37,572224

1,31103189

24,7669457

63,6485457

:

:

:

:

:

:

20

394,095033

232,625548

8,11908503

153,3504

394,095019

 

b)       Dos raices reales positivas y una negativa, alguna de módulo mayor que 1. En este caso la solución diverge a la solución particular:

γ1 =

0,8

 

m1=

1,58891618

 

γ2 =

0,1

 

m2=

1,07834917

0

i

Go =

4

 

m3=

-1,16726536

0

i

Co =

2

 

 

 

 

Io =

2

 

 

 

 

g =

-0,1

 

 

 

 

k1 =

0,7

 

 

 

 

k2 =

2

 

 

 

 

 

 

 

 

t

Y

C

I

G

Yp

0

8

2

2

4

11,480315

1

11,5

6,4

1,5

3,6

10,3322835

2

15,69

10

2,45

3,24

9,29905512

3

25,401

12,552

9,933

2,916

8,36914961

4

38,1229

20,3208

15,1777

2,6244

7,53223465

5

61,18761

30,49832

28,32733

2,36196

6,77901118

6

92,664949

48,950088

41,589097

2,125764

6,10111006

7

144,208704

74,1319592

68,1635573

1,9131876

5,49099906

8

216,124139

115,366963

99,0353066

1,72186884

4,94189915

9

327,877307

172,899311

153,428314

1,54968196

4,44770924

10

485,754647

262,301846

222,058087

1,39471376

4,00293831

:

:

:

:

:

:

20

21194,6427

11752,2974

9441,859

0,48630662

1,39573829

c)       Dos raices complejas conjugadas y una real negativa, todas de módulo menor que 1. En este caso la solución converge de manera oscilatoria a la trayectoria particular:

γ1 =

0,2

 

m1=

-0,83519329

 

γ2 =

0,5

 

m2=

0,71759665

0,28934212

i

Go =

4

 

m3=

0,71759665

-0,28934212

i

Co =

2

 

 

 

 

Io =

2

 

 

 

 

g =

0,2

 

 

 

 

k1 =

0,4

 

 

 

 

k2 =

0,5

 

 

 

 

 

 

 

 

t

Y

C

I

G

Yp

0

8

2

2

4

10,7329193

1

7,9

1,6

1,5

4,8

12,8795031

2

11,3

5,58

-0,04

5,76

15,4554037

3

14,432

6,21

1,31

6,912

18,5464845

4

19,7836

8,5364

2,9528

8,2944

22,2557814

5

24,83264

11,17272

3,70664

9,95328

26,7069376

6

31,49768

14,858328

4,695416

11,943936

32,0483252

7

38,2391152

18,715856

5,190536

14,3327232

38,4579902

8

46,625025

23,396663

6,02909408

17,1992678

46,1495882

9

55,8087655

28,4445626

6,7250815

20,6391214

55,3795059

10

67,1076624

34,4742656

7,8664511

24,7669457

66,4554071

:

:

:

:

:

:

20

411,441883

211,412154

46,6793298

153,3504

411,474364

d)                          Dos raices complejas conjugadas y una real negativa, algunas de módulo mayor que 1. En este caso la solución diverge de manera oscilatoria de la trayectoria particular:

γ1 =

0,4

 

m1=

-1,34909172

 

γ2 =

0,4

 

m2=

1,12454586

0,26387046

i

Go =

4

 

m3=

1,12454586

-0,26387046

i

Co =

2

 

 

 

 

Io =

2

 

 

 

 

G =

0,2

 

 

 

 

k1 =

0,5

 

 

 

 

k2 =

1,8

 

 

 

 

 

 

 

 

t

Y

C

I

G

Yp

0

8

2

2

4

36

 

1

9,5

3,2

1,5

4,8

43,2

 

2

13,51

7

0,75

5,76

51,84

 

3

20,821

9,204

4,705

6,912

62,208

 

4

32,9003

13,7324

10,8735

8,2944

74,6496

 

5

50,64125

21,48852

19,19945

9,95328

89,57952

 

6

75,973771

33,41662

30,613215

11,943936

107,495424

 

7

109,578702

50,6460084

44,5999705

14,3327232

128,994509

 

8

153,82126

74,2209892

62,4010034

17,1992678

154,793411

 

9

208,609262

105,359985

82,6101551

20,6391214

185,752093

 

10

276,76976

144,972209

107,030606

24,7669457

222,902511

 

:

:

:

:

:

:

 

20

1584,67735

1065,28365

366,043298

153,3504

1380,1536

 

30

6605,17549

3795,2227

1860,44754

949,505255

8545,5473

 

40

63863,649

36174,4243

21810,1385

5879,08627

52911,7764

 

50

372517,113

198431,433

137683,928

36401,7526

327615,773

 

60

3021269,59

1070376,95

1725502,58

225390,057

2028510,52

 

2.–     El Modelo de de inventario de Metzler.

            A diferencia del Modelo B3, el segundo enunciado en el cual se suponia que las espectativas de ventas son iguales al consumo del periodo anterior, se supone que las mismas son iguales al consumo del periodo anterior más una proporción del incremento del periodo anterior:
Ut = Ct−1 + r(Ct−1 ­− Ct−2) = βYt-1 + rβ(Yt−1 ­− Yt−2)
          Por lo cual la producción para existencias quedaría como:
St = (Ut − U t–1) = β •(Yt–1 − Y t–2) + rβ[(Yt−1 ­− Yt−2) ­− (Yt−2 ­− Yt−3)]
          Por lo tanto queda el siguiente sistema de ecuaciones:
                    Yt = Ut + St + V
                    Ut = βYt-1 + rβ(Yt−1 ­− Yt−2)                           (0 < β < 1)
                    St =  β(1+ r)Yt-1 − β•(1+2r)Y t–2 + rβ Yt−3   (0 < β < 1)

Solución del modelo:
          Sustituyendo las dos últimas ecuaciones en la primera:             
Yt = β•Y t–1 + rβ(Yt−1 ­− Yt−2) + β(1+ r)Yt-1 − β•(1+2r)Y t–2 + rβ Yt−3 + V
          reagrupando términos y normalizando:
                              Yt+3 − 2β•(1 + r)Y t+2 + β•(1 + 3r)Y t+2 − βr•Y t = V
Por lo tanto la solución vendrá dada por las raices de la ecuación cúbica:
x3 – 2β•(1 + r)• x2  + β•(1 + 3r)•x − βr = 0

y la solución particular por:
Yp = V/(1– β)

          Análisis de la solución:

                    Dado que el término independiente (βr) es negativo, y el coeficiente del término cuadrático (–2β•(1 + r)) es negativo también, tenemos los siguientes casos:
a)                          Dos raices complejas conjugadas y una real positiva, todas de módulo menor que 1. En este caso la solución converge de manera oscilatoria a la trayectoria particular:

b =

0,5

 

m1=

0,27330117

 

r =

0,5

 

m2=

0,61334941

0,73385575

i

U0 =

2

 

m3=

0,61334941

-0,73385575

i

U1 =

3

 

 

 

S0 =

2

 

 

 

V =

2

 

 

 

Y* =

4

 

 

 

 

 

 

 

 

t

Y

U

S

V

0

6

2

2

2

1

5,5

3

0,5

2

2

4,4375

2,625

-0,1875

2

3

3,6171875

1,953125

-0,3359375

2

4

3,42871094

1,60351563

-0,17480469

2

5

3,69909668

1,66723633

0,03186035

2

6

4,042099

1,91714478

0,12495422

2

7

4,20162773

2,10680008

0,09482765

2

8

4,15764403

2,14069605

0,01694798

2

9

4,03139111

2,06782609

-0,03643498

2

10

3,94228544

1,98413233

-0,04184688

2

:

:

:

:

:

20

4,00186251

1,99989594

0,00196657

2

b)                          Dos raices complejas conjugadas y una real positiva, una al menos de módulo mayor que 1. En este caso la solución diverge de manera oscilatoria a la trayectoria particular:


b =

0,7

 

m1=

0,34278286

 

r =

0,8

 

m2=

1,08860857

0,66979001

i

U0 =

2

 

m3=

1,08860857

-0,66979001

i

U1 =

3

 

 

 

S0 =

2

 

 

 

V =

2

 

 

 

Y* =

6,66666667

 

 

 

 

 

 

 

 

t

Y

U

S

V

0

6

2

2

2

1

5,7

3

0,7

2

2

6,3974

3,822

0,5754

2

3

7,6014308

4,868724

0,7327068

2

4

8,78383317

5,99525881

0,78857437

2

5

9,38172737

6,81082855

0,57089882

2

6

8,96587086

6,90202991

0,06384095

2

7

7,44206999

6,04322996

-0,60115997

2

8

5,17514389

4,35612051

-1,18097662

2

9

2,95102323

2,35312211

-1,40209888

2

10

1,74716929

0,82020869

-1,07303939

2

:

:

:

:

:

20

-9,42244632

-7,59102154

-3,83142478

2

          c)       Tres raices reales positivas, alguna de módulo mayor que 1. En este caso la solución diverge de solución particular:

b =

0,85

 

m1=

1,66194825

r =

1

 

m2=

1,36274446

U0 =

2

 

m3=

0,37530729

U1 =

3

 

 

 

S0 =

2

 

 

 

V =

2

 

 

 

Y* =

13,3333333

 

 

 

 

 

 

 

 

t

Y

U

S

V

0

6

2

2

2

1

5,85

3

0,85

2

2

8,41325

4,845

1,56825

2

3

15,1422963

9,330025

3,81227125

2

4

28,4621648

18,5906411

7,87152371

2

5

51,9002026

35,5147284

14,3854742

2

6

90,2818638

64,0375043

24,2443595

2

7

149,891514

109,363996

38,5275182

2

8

238,481185

178,07599

58,4051949

2

9

364,954329

278,010228

84,9441017

2

10

538,461009

417,713352

118,747656

2

:

:

:

:

:

17

2010,19803

2008,92326

-0,72522168

2

18

1144,95674

1540,83325

-397,876506

2

19

-867,855715

237,758133

-1107,61385

2

20

-4729,94511

-2448,56795

-2283,37717

2