Télécharger la plaquette Cast3M


fr en

Télécharger filc_test.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : filc_test.dgibi
  2. *
  3. *******************************************************
  4. ************************************************************************
  5. ************************************************************************
  6. *  Cas test de la procedure FILC
  7. *******************************************************
  8. *  Developpe par : 
  9. *     Alberto FRAU             (alberto.frau[at]cea.fr)     
  10. *     Benjamin RICHARD     (benjamin.richard[at]cea.fr)
  11. * 
  12. *******************************************************
  13. * Institution : 
  14. *     CEA/DEN/DANS/DM2S/SEMT/EMSI
  15. *******************************************************
  16. * Commentaires
  17. *
  18. *     On teste la procedure avec une signal composé de deux 
  19. *     harmoniques et on verifie que le signal filtré l'harmonique
  20. *     de frequence plus elevee est bien filtree
  21. *
  22. *******************************************************
  23. *
  24. * Parametres pour la definition du signal d'entree
  25. T0 = 0.;
  26. TMAX1 = 5.;
  27. DT1 = 1/((2.)*(40.));
  28. FO1 = 2.;
  29. FO2 = 10.;
  30. AMP1  = 2.;
  31. *
  32. GRAP1 = CHAINE 'N';
  33. *
  34. * Construction du signal par la somme des deux harmoniques (2Hz et 10Hz)
  35. LL1 = PROG T0 PAS DT1 TMAX1;
  36. VAL1 = ((((2)*(PI))*(FO1))/(PI))*(180.);
  37. VAL2 = ((((2)*(PI))*(FO2))/(PI))*(180.);
  38. LL2 = ((AMP1)*(SIN((VAL1)*(LL1)))) + ((SIN((VAL2)*(LL1))));
  39. EV1 = EVOL MANU LL1 LL2;
  40. EV1 = EV1 COUL BLEU;
  41. *
  42. * Filtrage à 5 Hz
  43. EV2 = FILC EV1 5. 0.1;
  44. EV2 = EV2 COUL ROUG;
  45. *
  46. SI ('NEG' GRAP1 'N');
  47.   DESS (EV1 ET EV2);
  48. FINSI;
  49. *
  50. * Test
  51. SI (((MAXI (EXTR EV2 ORDO)) - (AMP1)) > 0.01);
  52.   ERREUR 5;
  53. FINSI;
  54. *
  55. FIN;
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

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