Télécharger aftm.eso

Retour à la liste

Numérotation des lignes :

  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