Procedures | Cast3M

   1 : * FILC      PROCEDUR  AF221230  13/09/10    21:15:00     7820           
   2 : DEBPROC FILC EV1*'EVOLUTIO' F0*'FLOTTANT' DF1*'FLOTTANT';
   3 : *
   4 : **************************************************************
   5 : *  COMMENTAIRE:    Alberto FRAU (alberto.frauetcea.fr)
   6 : *                                    DEN/DANS/DM2S/SEMT/EMSI
   7 : *                  Benjamin RICHARD (alberto.frauetcea.fr)
   8 : *                                    DEN/DANS/DM2S/SEMT/EMSI
   9 : *
  10 : * La procedure filtre une signal pour une frequence F0 donnÃ©
  11 : * avec une intervalle DF1 de frequence sur lequel effectuer 
  12 : * couper le contenue frequentielle.
  13 : * 
  14 : * INPUT:
  15 : *    EV1 Evolution temporelle Ã  filtrer
  16 : *    F0  Frequence de coupure
  17 : *    DF1 intervalle autour de F01 pour aducir le filtrage
  18 : *
  19 : * OUTPUT:
  20 : *    EV3 Evolution temporelle filtrÃ©e
  21 : *    
  22 : *    | Modu               Fonction de filtrage                              
  23 : *    |                    
  24 : *    |                          DF1                          
  25 : *    |                    |------------|                         
  26 : *    |                                              
  27 : *    |____________________1
  28 : *    |                     * *                        
  29 : *    |                         *                     
  30 : *    |                           *              
  31 : *    |                           *                   
  32 : *    |                             *  * _____________0            
  33 : *    -------------------------------------------------
  34 : *    |---------------------------|                Freq
  35 : *                F0
  36 : **************************************************************
  37 : *
  38 : *
  39 : LX1 = 'PROG' 0. 'PAS' 0.05 1.;
  40 : LY1 = ((2.)*((LX1)**(3))) - ((3.)*((LX1)**(2))) 
  41 :          + ((1.)*('PROG' ('DIME' LX1) * 1.));
  42 : LX1 = LX1*DF1;
  43 : LX2 = 'PROG' 0. (F0 - (DF1/2));
  44 : LX2 = LX2 'ET' (LX1 + ('PROG' ('DIME' LX1) * F0));
  45 : LY2 = 'PROG' 1. 1.;
  46 : LY2 = LY2 'ET' (LY1);
  47 : *
  48 : NN1 = 'DIME' ('EXTR' EV1 'ABSC');
  49 : NN2 = 'ENTI'(('LOG' NN1)/('LOG' 2));
  50 : 'SI' (((2)**(NN2)) < NN1);
  51 :   NN2 = NN2 + 1;
  52 : 'FINSI';
  53 : *
  54 : TFR1 = 'TFR' NN2 EV1 'MOPH';
  55 : LF1 = 'EXTR' TFR1 'ABSC';
  56 : LM1 = 'EXTR' TFR1 'ORDO' 1;
  57 : LP1 = 'EXTR' TFR1 'ORDO' 2;
  58 : *
  59 : LX3 = LF1;
  60 : LY3 = 'IPOL' LF1 LX2 LY2;
  61 : FIL1 = 'EVOL' 'COMP' 'MOPH' LX3 LY3 LY3;
  62 : *
  63 : TFR2 = TFR1*FIL1;
  64 : *
  65 : EV2 = 'TFRI' TFR2;
  66 : *
  67 : TITX = 'EXTR' EV1 'LEGX' 1;
  68 : TITY = 'EXTR' EV1 'LEGY' 1;
  69 : *
  70 : NN3 = ('DIME' ('EXTR' EV2 'ABSC')) - NN1;
  71 : 'SI' (NN3 'NEG' 0);
  72 :   L_X1 = 'COPI' ('EXTR' EV2 'ABSC');
  73 :   L_Y1 = 'COPI' ('EXTR' EV2 'ORDO');
  74 :   'REPE' IND1 (NN3);
  75 :     L_X2 = ('ENLE' L_X1 ('DIME' L_X1));
  76 :     L_Y2 = ('ENLE' L_Y1 ('DIME' L_Y1));
  77 :     'DETR' L_X1;
  78 :     'DETR' L_Y1;
  79 :     L_X1 = 'COPI' L_X2;
  80 :     L_Y1 = 'COPI' L_Y2;
  81 :     'DETR' L_X2;
  82 :     'DETR' L_Y2;
  83 :   'FIN' IND1;
  84 : 'SINON';
  85 :   L_X1 = 'COPI' ('EXTR' EV2 'ABSC');
  86 :   L_Y1 = 'COPI' ('EXTR' EV2 'ORDO');
  87 : 'FINSI';
  88 : *
  89 : EV3 = 'EVOL' 'MANU' TITX L_X1 TITY L_Y1;
  90 : *
  91 : FINPROC EV3;
  92 : *
  93 : 
  94 : 
  95 : 
  96 :