*       POT_POUR.PRG : POT-POURRI DE SERIES NON STATIONNAIRES
SEED 4
ALLOCATE 300
*
COMPUTE NBEG=1, NEND=300, SHOCK=150
SET DSHOCK NBEG NEND = T==SHOCK
SET DU     NBEG NEND = T>SHOCK
PRINT SHOCK-1 SHOCK+3 CONSTANT DU
*
ZERO Y1
EQUATION 1 Y1
# CONSTANT Y1{1}
ASSOCIATE(VAR=1.0) 1
# 0.0 0.5
SIMULATE 1 NEND 2
# 1 Y1
*
ZERO Y2
EQUATION 2 Y2
# CONSTANT DU Y2{1}
ASSOCIATE(VAR=1.0) 2
# -1.0 2.0 0.5
SIMULATE 1 NEND 2
# 2 Y2 2
*
ZERO Y3
EQUATION 3 Y3
# CONSTANT Y3{1}
ASSOCIATE(VAR=1.0) 3
# 0.0 1.0
SIMULATE 1 NEND 2
# 3 Y3
*PRINT NBEG NEND Y3
*
SET TREND NBEG NEND = T
SET DTS   NBEG NEND = TREND*DU
SET DTS   SHOCK+1 NEND = DTS-SHOCK
PRINT SHOCK-1 SHOCK+3 DU TREND DSHOCK DTS
*
SET Y4 NBEG NEND = 0.0 + 0.05*TREND(T) + Y1(T)
*
SET Y5 NBEG NEND = 0.0 + 0.05*TREND - 0.04*DTS + Y1
*
SET Y6 NBEG NEND = 0.0 + 0.05*TREND - 0.04*DTS - 4*DU + Y1(T)
*
OPEN PLOT C:\T837\POT_POU.RGF
SPGRAPH(HFIELDS=3,VFIELDS=2,HEADER='Pot-pourri de séries  non stationnaires')
 GRAPH(HEADER='Y1 série stationnaire') 1
 # Y1 NBEG+1 NEND
 GRAPH(HEADER='Y2 série non stationnaire',$
      SUBHEAD='Deux moyennes',SAMESIZE) 1
 # Y2 NBEG+1 NEND
 GRAPH(HEADER='Y3 série non stationnaire',$
      SUBHEAD='Phi=1',SAMESIZE) 1
 # Y3 NBEG+1 NEND
 GRAPH(HEADER='Y4 série non stationnaire',$
      SUBHEAD='Tendance',SAMESIZE) 1
 # Y4 NBEG+1 NEND
 GRAPH(HEADER='Y5 série non stationnaire',$
      SUBHEAD='Tendance coudée',SAMESIZE) 1
 # Y5 NBEG+1 NEND
 GRAPH(HEADER='Y6 série non stationnaire',$
      SUBHEAD='Deux tendances + saut',SAMESIZE) 1
 # Y6 NBEG+1 NEND
SPGRAPH(DONE)
*
*              Présentation en classe
*
SPGRAPH(HFIELDS=3,VFIELDS=2,HEADER='Pot-pourri de séries  non stationnaires')
 GRAPH(HEADER='Y1') 1
 # Y1 NBEG+1 NEND
 GRAPH(HEADER='Y2',SAMESIZE) 1
 # Y2 NBEG+1 NEND
 GRAPH(HEADER='Y3',SAMESIZE) 1
 # Y3 NBEG+1 NEND
 GRAPH(HEADER='Y4',SAMESIZE) 1
 # Y4 NBEG+1 NEND
 GRAPH(HEADER='Y5',SAMESIZE) 1
 # Y5 NBEG+1 NEND
 GRAPH(HEADER='Y6',SAMESIZE) 1
 # Y6 NBEG+1 NEND
SPGRAPH(DONE)