Télécharger la plaquette Cast3M


fr en

Télécharger ftran.procedur

Retour à la liste

Numérotation des lignes : 

  1. * FTRAN     PROCEDUR  GF238795  17/11/10    21:15:06     9606           
  2. *
  3. DEBPROC FTRAN EV1*'EVOLUTIO' EV2*'EVOLUTIO' 
  4.            F0*'FLOTTANT' VAL1/'ENTIER  ';
  5. *
  6. *********************************************************************
  7. * 
  8. *  Procedure FTRAN
  9. *    Cette procedure calcule la fonction de transfert à partir de 
  10. *    deux evolutons
  11. *                          Alberto FRAU (alberto.frauetcea.fr)
  12. *                                        DEN/DANS/DM2S/SEMT/EMSI
  13. *                          Benjamin RICHARD (benjamin.richardetcea.fr)
  14. *                                        DEN/DANS/DM2S/SEMT/EMSI
  15. *    INPUT:
  16. *        EV1   premiere evolution
  17. *        EV2   deuxieme evolution
  18. *        F0    frequence de coupure
  19. *        VAL1  nombre des points à effacer pour la transformée
  20. *
  21. *    OUTPUT:
  22. *        EV3   fonction de transfert (module et phase)
  23. *
  24. *               EV3 = fft(EV1)/fft(EV2)
  25. *
  26. *
  27. *********************************************************************
  28. * 
  29. LIST VAL1;
  30. 'SI' (('EXIST' VAL1) 'ET' ('EGA' ('TYPE' VAL1) 'ENTIER'));
  31.   KK1 = VAL1;
  32. 'SINON';
  33.   KK1 = 1;
  34. 'FINSI';
  35. LIST KK1;
  36. *
  37. NN_EV1 = 'DIME' ('EXTR' EV1 'ABSC');
  38. NN_EV2 = 'DIME' ('EXTR' EV2 'ABSC');
  39. DT_EV1 = 'EXTR' ('EXTR' EV1 'ABSC') 2; 
  40. DT_EV2 = 'EXTR' ('EXTR' EV2 'ABSC') 2; 
  41. *
  42. 'SI' (NN_EV1 'NEG' NN_EV2);
  43.   'MESS' 'Erreur....les evolutions n ont pas le meme nombre de points';
  44.   'ERRE' 5;
  45. 'FINSI';
  46. 'SI' (DT_EV1 'NEG' DT_EV2);
  47.   'MESS' 'Erreur....les evolutions n ont pas le meme pas de temps';
  48.   'ERRE' 5;
  49. 'FINSI';
  50. *
  51. NN_FRT = 'ENTI' (('LOG' NN_EV1)/('LOG' 2));
  52. 'SI' (((2)**(NN_FRT)) < NN_EV1);
  53.   NN_FRT = NN_FRT + 1;
  54. 'FINSI';
  55.  
  56. TFR1 = 'TFR' NN_FRT EV1 'MOPH' 'FMAX' F0;
  57. TFR2 = 'TFR' NN_FRT EV2 'MOPH' 'FMAX' F0;
  58. *
  59. E_TR = TFR1/TFR2;
  60. *
  61. LL1 = 'EXTR' E_TR 'ABSC';
  62. LL2 = 'EXTR' E_TR 'ORDO' 1;
  63. LL3 = 'EXTR' E_TR 'ORDO' 2;
  64. 'REPE' IND1 (KK1);
  65.   LL1 = 'ENLE' LL1 1;
  66.   LL1 = 'ENLE' LL1 ('DIME' LL1);
  67.   LL2 = 'ENLE' LL2 1;
  68.   LL2 = 'ENLE' LL2 ('DIME' LL2);
  69.   LL3 = 'ENLE' LL3 1;
  70.   LL3 = 'ENLE' LL3 ('DIME' LL3);
  71. 'FIN' IND1;
  72. *
  73. E_TR = 'EVOL' 'COMP' 'MOPH' LL1 LL2 LL3;
  74. *
  75. EV3 = 'EVOL' 'MANU' 'Freq [Hz]' LL1 'Modu' LL2;
  76. *
  77. 'FINPROC' EV3;
  78. *********************************************************************
  79. * 
  80.  

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