Télécharger la plaquette Cast3M


fr en

Télécharger filc.procedur

Retour à la liste

Numérotation des lignes : 

  1. * FILC      PROCEDUR  BP208322  18/10/02    21:15:01     9943           
  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 (0.*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.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales