/* Simulation d'un BB et d'un AR1 */ data ar1; do i=1 to 200; retain y 0; u=rannor(100); ar1=0.6*y+u; y=ar1; output;end; run; /* représentation graphique dans le temps */ proc gplot; symbol i=join; plot u*i; plot y*i; run; /* Calcul et représentation de la densité spectrale */ PROC SPECTRA DATA=ar1 CENTER OUT=spectre P S WHITETEST ; VAR u ; WEIGHT PARZEN /*ou TUKEY ou BART */; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; RUN ; PROC SPECTRA DATA=ar1 CENTER OUT=spectre P PH S WHITETEST ; VAR y ; WEIGHT PARZEN ; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; RUN ; /* Simulation de processus MA */ data ma; u1=0; u2=0; do i = 1 to 200; u=rannor(1); ma1= u -0.8 * u1; ma2= u -0.8 * u1+0.2*u2; output; u1 = u; u2=u1; end; run; proc gplot; symbol i=join; plot u*i; plot ma1*i; plot ma2*i; run; PROC SPECTRA DATA=ma CENTER OUT=spectre P PH S WHITETEST ; VAR ma1 ma2 ; WEIGHT PARZEN ; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; PLOT S_02*FREQ; PLOT S_02*PERIOD ; RUN ; /* Simulation d'un processus non stationnaire (type déterministe) */ data ts; do i=1 to 200; u=rannor(17); ts=0.5+0.01*i+0.2*u; output;end; run; proc gplot; plot ts*i; run; PROC SPECTRA DATA=ts CENTER OUT=spectre P PH S WHITETEST ; VAR ts ; WEIGHT PARZEN ; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; RUN ; /* Simulation d'un processus non stationnaire (type stochastique) */ data ds; do i=1 to 200; retain y 0; u=rannor(40); ds=y+1.6*u; y=ds; output;end; run; proc gplot; plot ds*i; run; PROC SPECTRA DATA=ds CENTER OUT=spectre P PH S WHITETEST ; VAR ds ; WEIGHT PARZEN /*TUKEY BART 0.5 0 (ou 1 2 3 4 3 2 1) */; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; RUN ; /* Série d'electricité aux EU: nom de la base: elec dans la librairie a (que vous devez définir précédemment) */ proc gplot data=a.elec; plot elec*date; run; proc arima data=a.elec; identify var=elec; run; proc SPECTRA DATA=a.elec CENTER OUT=spectre P PH S WHITETEST ; VAR elec ; WEIGHT PARZEN ; RUN ; PROC GPLOT DATA=spectre ; PLOT S_01*FREQ; PLOT S_01*PERIOD ; RUN ;