Procedures | Cast3M

   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 :