Télécharger ftran.procedur

Retour à la liste

Numérotation des lignes :

  1. * FTRAN PROCEDUR AF221230 13/09/10 21:15:01 7820
  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 'ORDO') 2;
  40. DT_EV2 = 'EXTR' ('EXTR' EV2 'ORDO') 2;
  41. *
  42. 'SI' (NN_EV1 'NEG' NN_EV2);
  43. 'MESS' 'Erreur....les evolution n'ont pas les meme nobre des point';
  44. 'ERRE' 5;
  45. 'FINSI';
  46. 'SI' (DT_EV1 'NEG' DT_EV2);
  47. 'MESS' 'Erreur....les evolution n'ont pas les 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.  
  81.  
  82.  
  83.  

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