Télécharger la plaquette Cast3M


fr en

Télécharger aftm.eso

Retour à la liste

Numérotation des lignes : 

aftm

  1. C AFTM      SOURCE    BP208322  20/09/18    21:15:03     10718          
  2.         SUBROUTINE AFTM(NFFT,NHBM,GAM,IGAM,DL) 
  3. c
  4. C=======================================================================
  5. * AFTM: Construction des operateurs fréquentiels
  6. * GAM  : Transformée de Fourier Directe
  7. *
  8. * GAM = [ 1 , COS(th1) , SIN(th1) , ... , , COS(H*th1) , SIN(H*th1) ]
  9. *       [                           ...                             ] 
  10. *       [ 1 , COS(thN) , SIN(thN) , ... , , COS(H*thN) , SIN(H*thN) ]
  11. *
  12. * IGAM : Transformée de Fourier Inverse
  13. *
  14. * IGAM = (1/N) [     2      , ... ,       2    ]
  15. *              [  COS(th1)  , ... ,   COS(thN) ] 
  16. *              [  SIN(th1)  , ... ,   SIN(thN) ] 
  17. *              [              ...              ] 
  18. *              [ COS(H*th1) , ... , COS(H*thN) ] 
  19. *              [ SIN(H*th1) , ... , SIN(H*thN) ]
  20. *
  21. * DL : Dérivée
  22. *
  23. * DL = diag(0,DL1,...,DLH)
  24. *      pour tout j=1,...,H  DLj = j[ 0  1 ]
  25. *                                  [-1  0 ] 
  26. C=======================================================================
  27. c        
  28.        IMPLICIT INTEGER(I-N)
  29.        IMPLICIT REAL*8(A-H,O-Z)
  30.  
  31.         INTEGER NFFT,NHBM
  32.         REAL*8 PI,XNFFT
  33.         REAL*8 GAM(NFFT,2*NHBM+1), IGAM(2*NHBM+1,NFFT), DL(2*NHBM+1) 
  34.         PARAMETER(PI=3.141592653589793D0)
  35.  
  36. *======= GAM ===========================================================
  37. *        Boucle sur les harmoniques
  38.          DO I = 1,NFFT   
  39.                 GAM(I,1) = 1.
  40.          ENDDO
  41. *        Boucle sur les harmoniques
  42.          DO J = 2,2*NHBM,2
  43. *           Boucle sur le temps 
  44.             DO I = 1,NFFT
  45.                 GAM(I,J)   = COS(J*PI*(I-1)/NFFT)
  46.                 GAM(I,J+1) = SIN(J*PI*(I-1)/NFFT)
  47.             ENDDO
  48.          ENDDO
  49.  
  50. *======= IGAM ==========================================================
  51.          XNFFT = 1.D0/NFFT
  52. *        Boucle sur le temps 
  53.          DO J = 1,NFFT
  54.             IGAM(1,J) = XNFFT
  55. *           Boucle sur les harmoniques
  56.             DO I = 2,2*NHBM,2
  57.                 IGAM(I,J)   = 2.*COS(I*PI*(J-1)/NFFT)/NFFT
  58.                 IGAM(I+1,J) = 2.*SIN(I*PI*(J-1)/NFFT)/NFFT
  59.             ENDDO
  60.          ENDDO
  61.  
  62. *====== DL =============================================================
  63.         DL(1)=0.
  64.         DO J = 2,2*NHBM,2
  65.                 DL(J) = (J/2)
  66.                 DL(J+1) = -(J/2)
  67.         ENDDO
  68.         END
  69.  
  70.  
  71.  

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