Procedures | Cast3M

   1 : * MULTIREC  PROCEDUR  PICA      92/07/03    21:27:30     619            
   2 : *----------------------------------------------------------------------*
   3 : *                                                                      *
   4 : * PROCEDURE MULTIREC                                                   *
   5 : *                                                                      *
   6 : * CALCUL DE LA RECOMPOSITION D'UN SIGNAL A PARTIR D' UNE SEQUENCE DE   *
   7 : *        FILTRE ONDELETTE ORTHOGONALE                                  *
   8 : *                                                                      *
   9 : *----------------------------------------------------------------------*
  10 : *                                                                      *
  11 : * ENTREE                                                               *
  12 : *                                                                      *
  13 : * COARSE X/Y CONTENANT LE SIGNAL DE BASE (COARSE SIGNAL). DX DOIT ETRE *
  14 : *        CONSTANT DANS LE PROCESSUS DE DECOMPOSITION DIADIC.           *
  15 : *                                                                      *
  16 : * ONDELL X/Ys CONTENANT LA DECOMPOSITION EN ONDELETTE DANS L'ORDRE DES *
  17 : *        BANDES DE FREQUENCE CROISSANTE                                *
  18 : *                                                                      *
  19 : * H      Y CONTENANT LE FILTRE (NON FORCEMENT SYMETRIQUE).             *
  20 : * MMH    POINTS NEGATIFS (MMH<0--->SYMETRIQUE)                         *
  21 : *                                                                      *
  22 : * G      Y CONTENANT LE FILTRE ORTHOGONAL (NON FORCEMENT SYMETRIQUE).  *
  23 : * MMG    POINTS NEGATIFS (MMG<0--->SYMETRIQUE)                         *
  24 : *                                                                      *
  25 : *                                                                      *
  26 : * OPTIO  MOT CONTENANT LE MODE DE CORRELATION (DEFAUT='SYME')          *
  27 : *                                                                      *
  28 : * N      NOMBRE DE NIVEAU DE RECOMPOSITION DESIRE                      *
  29 : *                                                                      *
  30 : *----------------------------------------------------------------------*
  31 : *                                                                      *
  32 : * SORTIE                                                               *
  33 : *                                                                      *
  34 : * M      NOMBRE DE NIVEAU EFFECTIVEMENT ATTEINT                        *
  35 : *                                                                      *
  36 : * SIGNAL SIGNAL RECOMPOSE                                              *
  37 : *----------------------------------------------------------------------*
  38 : 'DEBPROC' MULTIREC COARSE*'EVOLUTION' ONDELL*'EVOLUTION'
  39 :                         H*'LISTREEL'   MMH*'ENTIER'
  40 :                         G*'LISTREEL'   MMG*'ENTIER'
  41 :                     OPTIO/'MOT'            N/'ENTIER';
  42 : *----------------------------------------------------------------------*
  43 : YPHI='COPIER' ('EXTRAIRE' COARSE 'ORDO' 1);
  44 : *
  45 : NPT='DIMENSION' YPHI;
  46 : NPTT='DIMENSION' ('EXTRAIRE' ONDELL 'ABSC' 1);
  47 : 'SI' (NPT 'NEG' NPTT);
  48 :   'MESSAGE' '*** nb de point COARSE neg nb de point ONDELL(1)';
  49 :   'MESSAGE' '--> pas de calcul';
  50 :   M=0;
  51 : 'SINON';
  52 :   M='DIMENSION' ONDELL;
  53 : 'FINSI';
  54 : *
  55 : 'SI' (MMH '<' 0); MMHT='TEXT' ' ';
  56 :                   MMH=('DIMENSION' H) - 1;
  57 :                   MPH= MMH;
  58 : 'SINON';          MMHT='TEXT' NPNE MMH;
  59 :                   MPH=('DIMENSION' H) - (MMH+1);
  60 : 'FINSI';
  61 : *
  62 : 'SI' (MMG '<' 0); MMGT='TEXT' ' ';
  63 :                   MMG=('DIMENSION' G) - 1;
  64 :                   MPG=MMG;
  65 : 'SINON';          MMGT='TEXT' NPNE MMG;
  66 :                   MPG=('DIMENSION' G) - (MMG+1);
  67 : 'FINSI';
  68 : *
  69 : NNPT=2*NPT;
  70 : 'SI' ((MMG '>' NNPT) 'OU' (MPG '>' NNPT));
  71 :   'MESSAGE' '*** nb de point insuffisant pour Phi au niveau 1';
  72 :   'MESSAGE' '--> pas de calcul';
  73 :   M=0;
  74 : 'FINSI';
  75 : 'SI' ((MMH '>' NNPT) 'OU' (MPH '>' NNPT));
  76 :   'MESSAGE' '*** nb de point insuffisant pour Psi au niveau 1';
  77 :   'MESSAGE' '--> pas de calcul';
  78 :   M=0;
  79 : 'FINSI';
  80 : *
  81 : 'SI' (M 'NEG' 0);
  82 :   'SI' (('NON' ('EXISTE' OPTIO)));
  83 :                OPTIO='TEXT' SYME; 'FINSI';
  84 :   'SI' (('NON' ('EXISTE' N))); N=1000; 'FINSI';
  85 :   'SI' (M > N); M=N; 'FINSI';
  86 : *
  87 :   I=0;
  88 :   'REPETER' LAB1 M;
  89 : *
  90 :     I=I+1;
  91 :     YPSI='EXTRAIRE' ONDELL 'ORDO' I;
  92 : *
  93 :     YYPSI='DIADIQUE' YPSI 'INVE';
  94 :     YYPHI='DIADIQUE' YPHI 'INVE';
  95 :     YPSI ='CVOLUTION' YYPSI G MMGT OPTIO;
  96 :     YPHI ='CVOLUTION' YYPHI H MMHT OPTIO;
  97 :     YYPHI=2 * (YPSI + YPHI);
  98 :     YPHI=YYPHI;
  99 : *
 100 :   'FIN' LAB1;
 101 : 'FINSI';
 102 : *
 103 : XXXX='EXTRAIRE' ONDELL 'ABSC' M;
 104 : XX  ='DIADIQUE' XXXX 'IVIN';
 105 : 'TITRE' 'Signal recompose sur ' M ' niveaux';
 106 : SIGNAL='EVOLUTION' 'MANUEL' 'temps' XX 'Y' YPHI;
 107 : 'FINPROC' M SIGNAL;