Télécharger filc.procedur

Retour à la liste

Numérotation des lignes :

  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.  
  97.  

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