Procedures | Cast3M

   1 : * COURSPEC  PROCEDUR  PICA      92/07/03    21:18:21     619            
   2 : ********************************************************************************
   3 : * PROCEDURE COURSPEC
   4 : * 
   5 : * spectre de puissance "stationaire" associe a une decomposition
   6 : * en ondelette
   7 : *
   8 : * STASPE=COURSPEC SPECT DELTT;
   9 : *
  10 : ********************************************************************************
  11 : * ENTREE
  12 : *
  13 : * SPECT  : valeur su spectre dans les bandes (calculer par VALSPE ou RESPOWNS)
  14 : * DELTT  : pas de temps du residu
  15 : *
  16 : ********************************************************************************
  17 : * SORTIE
  18 : *
  19 : * STASPE : spectre
  20 : *
  21 : ********************************************************************************
  22 : 'DEBPROC' COURSPEC SPECT*'LISTREEL' DELTT*'FLOTTANT';
  23 : *
  24 : M=('DIME' SPECT) - 1;
  25 : *
  26 : XIN=1./((2**M)*(2*DELTT));
  27 : XPROUT='EXTR' SPECT 1;
  28 : *
  29 : VALS='PROG' XPROUT XPROUT 0.;
  30 : XVAL='PROG' .01001 (XIN-0.0001) XIN;
  31 : *
  32 : I=1;
  33 : 'REPETER' LAB1 M;
  34 :  I=I+1;
  35 :  XPROUT='EXTR' SPECT I;
  36 :  VALS=VALS 'ET' ('PROG' XPROUT XPROUT 0.);
  37 :  XXIN=2*XIN;
  38 :  XVAL=XVAL 'ET' ('PROG' (XIN+0.0001) (XXIN-0.0001) XXIN);
  39 :  XIN=XXIN;
  40 : 'FIN' LAB1;
  41 : *
  42 : STASPE=EVOL MANU 'Freq' XVAL 'Spec' VALS;
  43 : PUISS=EXTR (SOMM STASPE) 1;
  44 : TITRE 'Generic Shape of the Stationary Power Spectrum: Power=' PUISS;
  45 : STASPE=EVOL MANU 'Freq' XVAL 'Spec' VALS;
  46 : *
  47 : 'FINPROC' STASPE;