Brevet FR2998064A1 — publié le 20140516

ETAT DE L'ART Domaine technique de l'invention La présente invention concerne un dispositif et une méthode d'holographie 5 hétérodyne numérique hors-axe, notamment pour la mesure de vibrations mécaniques d'amplitudes nanométriques. Etat de l'art L'holographie est une méthode qui permet de reconstruire une image en trois dimensions (3D) d'un objet à partir de l'enregistrement d'une figure d'interférence entre 10 deux ondes cohérentes, une onde de référence et une onde diffusée par l'objet à observer. Dans le cas de l'holographie classique, l'enregistrement est effectué sur une plaque photographique. Une fois la plaque développée par des procédées connus de photographie, on obtient un hologramme. L'image 3D est reconstruite en éclairant l'hologramme avec l'onde de référence qui produit alors une onde transmise directement (ordre 0) et deux 15 ondes de lumière diffractées, une onde objet (ou onde signal) qui reconstitue une image virtuelle de l'objet et une onde conjuguée qui forme une image réelle parasite qui gène l'observation de l'image virtuelle et que l'on cherche à supprimer. En holographie classique, l'utilisation par exemple d'un matériau épais pour l'émulsion photographique utilisée pour l'enregistrement de l'hologramme permet de créer un hologramme de volume 20 qui, lors de la restitution, ne va diffracter que l'onde signal recherchée. On obtient ainsi des images de très grande qualité. Cependant, en holographie classique, les informations quantitatives enregistrées (amplitude et phase de l'onde signal) ne sont pas accessibles directement, ce qui en limite les applications. Afin de disposer d'informations quantitatives exploitables, une nouvelle technique 25 d'holographie est apparue, appelée holographie numérique, dans laquelle le film photographique est remplacé par un détecteur optoélectronique bidimensionnel, par exemple une caméra CCD. L'acquisition du signal d'interférence par le détecteur permet de former un interférogramme, c'est-à-dire une figure d'interférence échantillonnée sur le détecteur bidimensionnel. Le procédé de reconstruction, entièrement numérique, consiste à 30 calculer à partir de l'interférogramme le champ de l'onde diffractée dans un plan de l'objet, par exemple au moyen d'une transformée de Fresnel, qui peut, par exemple, comprendre un algorithme de transformées de Fourier rapides (FFT selon l'abréviation de l'expression anglo-saxonne « Fast Fourrier Transform »). Un exemple d'holographie numérique est exposé dans Schnars et Jüptner, « Direct recording of holograms by a CCD target and numerical reconstruction », Appl. Opt. 33, p. 179-181 (1994). L'holographie numérique présente de nombreux avantages par rapport à l'holographie classique, du fait notamment que le champ reconstruit de l'onde signal est une fonction complexe permettant d'accéder à l'amplitude et à la phase de l'onde signal. Il est ainsi possible d'établir la cartographie tridimensionnelle d'un objet, l'étude d'objets en mouvement et notamment l'étude d'un objet en vibration. Cependant, de la même façon qu'en holographie classique sur film mince, le champ de l'onde signal calculé à partir de l'interférogramme comprend simultanément l'ordre diffractif recherché, l'ordre 0 et l'ordre conjugué parasite. Il y a donc une incertitude sur l'amplitude complexe de l'onde objet que l'on cherche à déterminer. Un moyen pour supprimer l'image conjuguée parasite et l'onde non-diffractée est d'introduire un déphasage entre l'onde de référence et l'onde signal. Dans la demande EP 1 043 632 par exemple, une technique dite d'holographie hétérodyne numérique est décrite, notamment pour l'application à l'étude d'objets en vibration. Selon cette technique, la fréquence de l'onde de référence est décalée d'une quantité déterminée en fonction de la fréquence de vibration de l'objet et de la cadence d'acquisition du détecteur, par exemple au moyen de deux modulateurs acousto-optiques. Une démodulation numérique de l'interférogramme permet alors de supprimer les composantes dues à l'ordre 0 et à l'ordre conjuguée. Pour gagner encore en précision, il a été proposé (voir par exemple M. Gross et al., « Digital Holography with ultimate sensibility », Optics Letters, Vol. 32, (2007)), de combiner la technique d'holographie hétérodyne numérique avec une technique d'holographie hors-axe, dans laquelle un angle (typiquement quelques degrés) est introduit entre l'onde de référence et l'onde diffusée par l'objet. Couplée avec un filtrage spatial numérique, cette technique permet de filtrer l'onde transmise à l'ordre zéro ainsi que le bruit qui lui est attaché. Il est alors possible d'enregistrer et de reconstruire des images holographiques avec un meilleur rapport signal sur bruit, permettant la mesure de vibrations de très faible amplitude, à cadence d'acquisition vidéo. Dans l'article de V. Lédl et al. (« Frequency Shifted Digital Holography for the Measurement of Vibration with Very Small Amplitudes », 9th International Conference on Vibration Measurements by Laser and Noncontact Techniques, juin 2010), des mesures d'amplitudes de vibration de l'ordre du nanomètre ont ainsi été réalisées expérimentalement, en utilisant un montage d'holographie hétérodyne hors-axe. L'invention propose une méthode et un dispositif d'holographie hétérodyne hors-axe permettant, par rapport aux technologies connues, d'accéder à des rapports signal sur bruit supérieurs et de gagner de ce fait en précision de mesure, notamment sur les amplitudes de vibration inférieures au nanomètre. RE SUME DE L'INVENTION Selon un premier aspect, l'invention concerne une méthode d'holographie numérique pour détecter l'amplitude z de vibration d'un objet à au moins une fréquence de vibration co, comprenant : - la génération d'une onde laser de fréquence (pli, d'illumination de l'objet et d'une onde de référence cohérente avec l'onde d'illumination; - l'acquisition au moyen d'un détecteur optoélectronique bidimensionnel de bande passante temporelle cos d'une série d'interférogrammes résultant de l'interférence de l'onde de référence et d'une onde signal issue de l'objet, les ondes signal et de référence présentant un décalage angulaire, le champ de l'onde de référence comprenant une première composante de fréquence coi et au moins une deuxième composante de fréquence (02, les fréquences coi et c02 étant décalées de la fréquence coL de l'onde laser respectivement d'une quantité déterminée Ôi, Ô2, avec Ô1 = yicos et Ô2 = qw + y2cos, où q entier relatif et -0,5 La déposante a montré que grâce à la méthode ainsi décrite, un rapport signal sur bruit nettement amélioré des mesures d'amplitudes de vibration pouvait être obtenu, notamment grâce à la mesure simultanée des raies du spectre en fréquence des hologrammes. De très bons résultats sont obtenus en travaillant sur l'une des premières bandes de modulation de l'onde signal (q = ±1) et avec une onde de référence présentant deux composantes décalées en fréquence. Il est également possible d'augmenter le rapport signal sur bruit en augmentant le nombre de composantes de l'onde de référence et de travailler sur des bandes de modulation d'ordre plus élevé, par exemple pour accéder à des amplitudes de vibration plus grandes. De façon générale, yi et y2 peuvent s'écrire respectivement plIN - 1/2 et p2IN - 1/2 où N est le nombre d'interférogrammes dans une série d'interférogrammes acquis par le détecteur optoélectronique bidimensionnel et p1,2 = 1,... N. Selon une variante, yi et y2 sont de normes différentes pour éviter des effets de diaphonie (ou « cross-talk » selon l'expression anglo-saxonne) sur le spectre en fréquence des battements optiques entre les composantes fréquentielles de l'onde de référence et de l'onde signal. Par exemple, l'onde de référence pourra comprendre au moins une troisième composante EL03 de fréquence co3, la fréquence (o3 étant décalée de la fréquence coL de l'onde laser d'une quantité déterminée Ô3 = y3cos + qw, où q entier relatif et - 0,5 Selon une variante, la méthode comprend l'émission d'une onde laser de fréquence COL et la séparation de l'onde d'émission laser pour former l'onde d'illumination de l'objet et une deuxième onde envoyée vers un dispositif optoélectronique de décalage de fréquence optique pour former l'onde de référence, permettant ainsi de former des ondes d'illumination et de référence cohérentes. Selon une variante, la vibration de l'objet peut résulter de l'excitation forcée de l'objet à ladite fréquence w. Il est alors possible grâce à la méthode ainsi décrite, de déterminer, en plus de l'amplitude de vibration, le retard de phase entre la phase de vibration de l'objet et la phase du signal d'excitation, et ceci à partir du rapport des amplitudes de la première raie et de la deuxième raie du spectre en fréquence des battements optiques entre les composantes fréquentielles de l'onde de référence et de l'onde signal. Selon une variante, dans le cas d'une excitation forcée de l'objet, la fréquence d'excitation est variable temporellement de façon continue entre une fréquence d'excitation inférieure coi et une fréquence d'excitation supérieure OF pendant une durée (T) donnée, l'acquisition et le traitement des interférogrammes étant effectuées pendant ladite durée. La variation de la fréquence d'excitation est par exemple linéaire entre les valeurs inférieure et supérieure. Cette variation continue de la fréquence d'excitation, simultanément à une variation d'une des fréquences de décalage de l'onde de référence, permet un verrouillage du retard de phase de la vibration par rapport à l'excitation, permettant une mesure plus précise du retard de phase. Alternativement, la méthode peut s'appliquer aussi dans le cas d'un objet présentant un spectre de fréquences de vibrations. L'acquisition et le traitement des interférogrammes sont alors effectués pour obtenir l'amplitude z de vibration de l'objet à au moins une fréquence de vibration (o contenue dans le spectre. Cette méthode permet alors d'établir par exemple des cartographies d'objets en vibration et la détermination des points résonants, notamment dans le cadre de la vibrométrie laser pour le contrôle non destructif Selon un deuxième aspect, l'invention concerne un dispositif d'holographie numérique pour détecter l'amplitude z de vibration d'un objet à une fréquence de vibration (o, comprenant : - un détecteur optoélectronique bidimensionnel de bande passante temporelle (os; - des moyens de génération de deux ondes laser cohérentes entre elles, une onde laser d'illumination de l'objet de fréquence coL et une onde de référence, le champ de l'onde de 25 référence présentant au moins une première et une deuxième composantes de fréquences to2 décalées de la fréquence (oL de l'onde laser respectivement d'une quantité déterminée ôi, ô2, avec ô1 = yi(os et ô2 = cgo + y2cos, où q entier relatif et - des moyens de combinaison de l'onde de référence et d'une onde signal issue de 30 l'objet sur le détecteur optoélectronique bidimensionnel, les ondes de référence et signal présentant un décalage angulaire, l'acquisition du signal d'interférence des ondes signal et de référence par le détecteur optoélectronique bidimensionnel résultant en une série d' interférogrammes; et - des moyens de traitement des interférogrammes comprenant : - des moyens de calcul d'hologrammes de l'objet obtenus à partir de chacun des interférogrammes, un hologramme étant défini par l'amplitude complexe de l'onde signal dans un champ spatial donné d'un plan de l'objet; - des moyens pour calculer la transformée de Fourier temps - fréquence des amplitudes complexes desdits hologrammes, permettant d'obtenir en chaque point du champ spatial un spectre en fréquence des battements optiques entre les composantes fréquentielles de l'onde de référence et de l'onde signal, le spectre de battements optiques comprenant au moins deux raies pour des fréquences égales aux décalages de fréquence Ô1 et Ô2 ; et - des moyens pour calculer en chaque point du champ spatial à partir des amplitudes de la première raie et de la deuxième raie l'amplitude z de la vibration de l'objet. Selon une variante, le dispositif comprend une source d'émission laser à la fréquence coL, des moyens de séparation de l'onde d'émission laser pour former l'onde d'illumination de l'objet et une seconde onde et un dispositif optoélectronique de décalage de fréquence optique, destiné à recevoir ladite seconde onde pour former l'onde de référence. Selon une variante, les moyens de séparation de l'onde laser sont fibres. Par exemple, le dispositif optoélectronique de décalage de fréquence optique comprend deux modulateurs acousto-optiques travaillant à des fréquences prédéterminées sur des ordres de diffraction opposés pour former ladite onde de référence. Alternativement, le dispositif optoélectronique de décalage de fréquence optique comprend une cellule de Pockels. La cellule de Pockels présente l'avantage d'être moins coûteuse que des modulateurs acousto-optiques et de ne pas nécessiter d'électronique complexe. Selon une variante, le dispositif comprend en outre un moyen d'excitation de l'objet à la fréquence co de vibration donnée. Selon une variante, le détecteur optoélectronique bidimensionnel est une caméra de type CCD ou CMOS. BREVE DESCRIPTION DES DESSINS D'autres avantages et caractéristiques de l'invention apparaîtront à la lecture de la description, illustrée par les figures suivantes qui représentent : - la figure 1, un schéma d'un exemple de dispositif d'holographie numérique selon 1 ' invention ; - les figures 2A à 2C, des schémas illustrant des étapes de traitement des 5 interférogrammes dans un exemple de méthode d'holographie numérique selon 1 ' invention ; - les figures 3A à 3C, des diagrammes montrant dans un exemple de mise en oeuvre des mesures de spectres en fréquence du champ de l'onde signal (Figure 3A), du champ de l'onde de référence (Figure 3B) et des battements optiques entre les spectres en fréquence 10 des champs de l'onde signal et de référence (Figure 3C), obtenus par la méthode d'holographie numérique selon l'invention ; - la figure 4, un diagramme montrant des courbes comparatives de mesures d'amplitude z de vibration d'un objet, effectuées respectivement à l'aide de la méthode selon l'invention et selon une méthode de l'art antérieur; 15 - la figure 5, un schéma illustrant un exemple particulier de réalisation d'un dispositif selon l'invention ; - les figures 6A à 6D, des hologrammes d'un objet obtenus avec le dispositif de l'exemple de la figure 5 ; - la figure 7, un diagramme représentant, en fonction de la fréquence d'excitation de 20 l'objet, l'amplitude mesurée de vibration en un point donné de l'objet ; et - les figures 8A et 8B, des diagrammes représentant, en fonction de la fréquence d'excitation de l'objet, respectivement l'amplitude de vibration et la phase, les courbes mesurées étant comparées aux courbes théoriques. DESCRIPTION DETAILLEE 25 La figure 1 illustre un schéma d'un exemple de dispositif d'holographie numérique selon l'invention. Dans cet exemple, le dispositif est fibre. Le dispositif permet l'enregistrement et le traitement d'interférogrammes d'on objet présentant au moins une fréquence de vibration co, afin d'en déterminer l'amplitude et/ou la phase de vibration. Dans la suite de la demande, on parlera simplement de « fréquence de vibration » co pour 30 exprimer la fréquence angulaire (ou pulsation). La phase de vibration est plus précisément définie, dans le cas d'un objet soumis à une excitation forcée sinusoïdale, par le retard de phase par rapport à la phase d'excitation. La capacité à mesurer des amplitudes de vibration de très faibles amplitudes d'un objet ainsi que la phase de la vibration de l'objet soumis à une vibration forcée peuvent permettre de donner des informations très fines sur la nature de l'objet, telles que son intégrité mécanique, sa qualité et, le cas échéant, son état d'endommagement. Pour cela, le dispositif selon l'invention comprend des moyens 3, 7 de génération de deux ondes cohérentes temporellement, une onde d'illumination de l'objet W1 et une onde de référence WLO. Le dispositif comprend en outre un détecteur bidimensionnel 19 de bande passante temporelle os. Avantageusement, les moyens de génération du dispositif selon l'invention comprennent une source 3 cohérente temporellement, par exemple une source laser, émettant une onde d'émission WEM ayant une fréquence optique wu La source laser est, par exemple, une diode laser émettant à une longueur d'onde d'émission dans la bande spectrale visible ou proche infrarouge, par exemple une diode laser fibrée. Les moyens de génération du dispositif selon l'invention comprennent en outre un séparateur de faisceaux 7, permettant de former, à partir de l'onde d'émission WEm, l'onde d'illumination de l'objet W1 et une seconde onde WL, ayant chacune la fréquence optique wu Dans l'exemple de la figure 1, le séparateur de faisceaux 7 comprend un coupleur à fibres, permettant de coupler l'onde d'émission WEm issue d'une fibre optique 5 de sortie de la source d'émission 3 dans deux fibres optiques 9, 11 afin de former l'onde d'illumination W1 et la seconde onde WL. L'onde d'illumination W1 éclaire l'objet 15 que l'on souhaite examiner, soit en réflexion, soit en transmission. L'onde issue de l'objet, rétrodiffusée ou transmise, est appelée onde signal dans la suite de la demande et désignée Ws sur la figure 1. La vibration de l'objet 15 résultant par exemple d'une excitation forcée de l'objet à la fréquence co entraine une modulation de la phase (t) de la partie temporelle du champ E de l'onde signal Ws due à la variation du chemin optique entre l'objet 15 en vibration et le détecteur 19. Ceci conduit à l'apparition de bandes latérales optiques d'amplitudes complexes sr, aux harmoniques de la fréquence de vibration co de l'objet. La modulation de la phase (t) en fonction du temps t peut avoir la forme 4(t) = osin(cot + y), où 00 = 47cz/X est l'amplitude de modulation, exprimée en fonction de l'amplitude de vibration z de l'objet et de la longueur d'onde X de l'onde d'émission WEm, et y est le retard de phase mécanique du signal de vibration par rapport au signal d'excitation. Ainsi, le champ optique E de l'onde signal Ws peut s'écrire de la façon E =InEn , avec En et n entier relatif, n G [-00,00]. Les quantités s n = e n(0 0)em sont les amplitudes complexes des composantes fréquentielles d'ordre n de l'onde signal, avec s l'amplitude complexe du champ optique de l'onde signal et ./0 les fonctions de Bessel de première espèce. On parle de composante non décalée E0 du champ E de l'onde signal pour n = 0 et de bandes latérales pour les composantes En du champ E de l'onde signal, où n=1, -1, 2, - 2, etc. Il est alors possible de déduire des amplitudes complexes des composantes fréquentielles d'ordre 0 et 1 de l'onde signal, l'amplitude de vibration z de l'objet et le retard de phase y par l'équation : 27z- (1) 1 z e . so L'équation (1) est valable pour des amplitudes de vibration très inférieures à la longueur d'onde X de l'onde d'émission. Il en résulte notamment une expression de l'amplitude de vibration z de l'objet en fonction de l'amplitude relative de la première bande par rapport à la composante non décalée selon l'équation : 1E11 z (2) 27z- Dans l'exemple illustré sur la figure 1, les moyens de génération des ondes d'illumination de l'objet et de référence comprennent en outre un dispositif optoélectronique de décalage de fréquence 13, destiné à recevoir la seconde onde WL afin d'en décaler la fréquence optique coL pour former l'onde de référence WLO et un moyen de contrôle 23 du dispositif optoélectronique de décalage de fréquence 13. L'onde de référence WL0 décalée en fréquence est appelée oscillateur local dans la suite de la description. Le décalage de la fréquence coL de l'onde de référence permet une démodulation numérique des interférogrammes à une fréquence de démodulation contenue dans la bande passante temporelle cos du détecteur optoélectronique, qui correspond à sa cadence d'acquisition. Selon la présente invention, l'oscillateur local, ou onde de référence, présente une première fréquence coi et au moins une deuxième fréquence c02, les fréquences étant décalées de la fréquence d'émission coL d'une quantité déterminée ôj, j = 1...m, de telle sorte à ce que lors de la démodulation des interférogrammes, la composante non décalée du champ de l'onde signal et les bandes latérales utiles pour la détermination de l'amplitude de vibration z de l'objet soient contenues dans la bande passante cos de la caméra. Ainsi, le champ de l'oscillateur local peut s'écrire EL0=IjELoj ..., où ELoj est /st la composante j = 1.. .m ayant la forme ELoj = EL I'Ltei , avec ôj le décalage de fréquence. Comme cela sera expliqué plus en détails par la suite, la déposante a montré qu'un oscillateur local à deux fréquences (m = 2), appelée dans la suite de la description « oscillateur local dual », permet d'obtenir un rapport signal sur bruit nettement meilleur pour la mesure de l'amplitude de vibration ou du retard de phase que ceux obtenus dans l'art antérieur et donc d'atteindre notamment des mesures d'amplitudes de vibration très faibles. En choisissant pour l'oscillateur local un nombre supérieur de fréquences de décalage (m > 2), on peut augmenter encore le rapport signal sur bruit. Par exemple, dans le cas d'un oscillateur local dual, les fréquences décalées de l'onde de référence WLO pourront être données par (oi = coL + Ô1 et CO2 = COL ± Ô2, avec ô1 = yitos et Ô2 = cgo + y20)s, où q entier relatif et - 0,5 Selon une variante, on choisit yi, y2 de normes différentes pour éviter des effets de diaphonie (ou « cross-talk » selon l'expression anglo-saxonne) sur le spectre en fréquence des battements optiques entre l'onde de référence et l'onde signal. Selon une variante, le dispositif optoélectronique de décalage de fréquence 13 peut comprendre deux modulateurs acousto-optiques (AOM) fonctionnant à des fréquences déterminées sur des ordres de diffraction opposés pour obtenir les fréquences décalées de l'onde de référence. Le moyen de contrôle 23 peut comprendre un synthétiseur numérique permettant de programmer les fréquences synthétisées des signaux de commande des AOM. Les fréquences de décalage peuvent ainsi facilement être adaptées à la fréquence de vibration co de l'objet. Selon une variante, le dispositif optoélectronique de décalage de fréquence 13 peut comprendre une cellule de Pockels, permettant la génération de l'oscillateur local aux fréquences recherchées. Plus précisément, la cellule de Pockels est commandée par un signal destiné à créer au moins deux fréquences de décalage de l'onde de référence WLO. Ce signal peut par exemple être un balayage linéaire de tension de commande avec le temps, dont la pente commute séquentiellement entre deux valeurs au moins une fois pendant le temps d'acquisition d'une image, les deux pentes étant les deux décalages de fréquence ôi et ô2 de l'oscillateur local dual. L'interférence entre l'onde de référence ou oscillateur local WLO et l'onde signal Ws, recombinées à l'aide d'un moyen de combinaison 21, par exemple un cube séparateur, mène à des figures d'interférence qui sont acquises par le détecteur bidimensionnel 19. Le moyen de combinaison 21 est arrangé de manière à ce que les deux ondes interfèrent avec un angle 0 de l'ordre du degré. Ce montage hors-axe permet, lors du traitement des interférogrammes enregistrés, d'effectuer un filtrage spatial dans le domaine des fréquences spatiales k afin de supprimer l'image parasite conjuguée et l'ordre 0 qui apparaissent lors de la reconstruction des hologrammes. Grâce aux choix des fréquences de l'oscillateur local, l'interférence entre l'onde signal et l'onde de référence donne lieu à des battements optiques à des fréquences comprises dans la bande passante temporelle os du détecteur 19. Le dispositif d'holographie numérique selon l'invention comprend également des moyens de traitement 25 des interférogrammes enregistrés par le détecteur bidimensionnel 19. Les figures 2A à 2C illustrent schématiquement les étapes mises en oeuvre pour le calcul des hologrammes de l'objet, permettant ensuite d'accéder à l'amplitude z et/ou à la phase y de vibration de l'objet. Un nombre N d'interférogrammes /p est enregistré par intervalle de temps par le détecteur bidimensionnel, comme illustré par la figure 2A. Par exemple, 256 images sont acquises à une fréquence (os/(27c) = 20 Hz. Les fréquences décalées de l'onde de référence sont choisies telles que yi et y2 soient égaux respectivement à plIN - 1/2 et p2IN - 1/2 où p1,2 = 1,... N. Les moyens de traitement permettent alors de calculer lors d'une étape de reconstruction un hologramme de l'objet à partir de chaque interférogramme enregistré. Un interférogramme enregistré par le détecteur a la forme 2 1En -FIE LOj (3) La reconstruction de l'hologramme est faite par exemple de façon connue par transformée de Fresnel de l'interférogramme /p, permettant d'obtenir pour chaque interférogramme un hologramme à valeurs complexes ip = Fresnel j dans un champ spatial donné d'un plan de l'objet. Chaque hologramme est ainsi défini par l'amplitude complexe de l'onde signal dans le champ spatial donné de l'objet. De manière avantageuse, seulement une partie décentrée Hp (dite hors-axe) des images ip est conservée dans la suite, comme illustré sur la figure 2B, les hologrammes hors-axe Hp étant obtenus par filtrage numérique dans l'espace des fréquences spatiales, pour ne garder que l'ordre de diffraction utile lors de la reconstruction. Les hologrammes Hp sont stockés (figure 2C) puis une transformée de Fourier temps - fréquence des amplitudes complexes des hologrammes Hp est calculée, pour chaque point du champ, afin d'obtenir le spectre en fréquence fi FT{H}, montré schématiquement sur la figure 2C. La composante k de la transformée de Fourier peut être écrite comme suit : (C0k) pN p exp(-2ipkre /N) (4) On obtient ainsi un spectre fi(co,,) en fréquence des battements optiques entre les composantes fréquentielles de l'onde de reférence WL0 et de l'onde signal Ws pour chaque point du champ. On montre que chaque composante fi(co, = 7,cos)peut s'écrire : 5I-1(cok 7kcos) K k -1 e * LOk (5) où K est une constante. Par exemple dans le cas de l'oscillateur local dual, 1-1(71cos ) K eo E et 11(720,)= K * LO2 - Les figures 3A à 3C montrent des diagrammes représentant des spectres fréquentiels mesurés de l'amplitude E du champ de l'onde signal Ws (figure 3A), de 10 l'amplitude ELO du champ de l'onde de référence WLO (figure 3B) et des battements optiques entre les composantes fréquentielles de l'onde de référence WL0 et de l'onde signal Ws pour un point donné du champ (figure 3C) dans un exemple de mise en oeuvre. Dans cet exemple, l'oscillateur local est dual et comprend deux fréquences de décalage coi = coL + cos/4 et (02 = COL ± CO - cos/4 obtenus par exemple au moyen de deux 15 modulateurs acousto-optiques dont l'un fonctionne à la fréquence porteuse (oc, et le deuxième fonctionne aux fréquences (oc - cos/4 et (oc + co + cos/4. Les modulateurs acoustooptiques sont par exemple des cellules Bragg fonctionnant à la fréquence porteuse w/(2r) = 80 MHz. La source d'émission est une diode laser fonctionnant à 532 nm. Les ondes signal et de référence sont combinés sur le détecteur avec un angle de quelques dégrés. Le 20 détecteur est une caméra CCD de fréquence d'acquisition cos/(27c) = 20 Hz. L'objet est un actuateur piézoélectrique présentant une oscillation hors plan de fréquence co/(27c) = 10 kHz. Une tension variable comprise entre 10-2 V et 10 V est appliquée à l'actuateur piézoélectrique pour faire varier l'amplitude de vibration. Pour chaque tension, une série de 256 interférogrammes /p est acquise par la caméra CCD puis traitée selon la méthode 25 décrite ci-dessus. Les figures 3A et 3B illustrent les spectres en fréquences de l'onde signal et de l'onde de référence (oscillateur local). Sur la figure 3A, on observe la composante non décalée 1E01 de l'onde signal (ramenée à la fréquence 0) et les bandes latérales 1E11 et Is aux fréquences +co et - co respectivement. Sur la figure 3B, on observe les bandes latérales 30 de l'oscillateur local, respectivement IELoil et IEL02, correspondant respectivement à des décalages par rapport à la fréquence coL de l'onde d'émission laser de Ô1 = cos/4 et Ô2 =w - cos/4 Chaque interférogramme /p formé dans le plan du détecteur présente une expression donnée par l'équation (3) ci-dessus. Lors de la détection, seules les composantes fréquentielles comprises entre ± cos/2 sont visibles par le détecteur. Dans l'exemple ci- dessus, le détecteur enregistre la composante non décalée Iso et la première composante décalée 1E11 de l'onde signal modulée par la vibration de l'objet simultanément, grâce à la transposition des fréquences de l'onde signal Ws dans la bande passante temporelle du détecteur au moyen de l'onde de référence Ww. Ainsi, sur la figure 3C l'amplitude de la raie A est proportionnelle à le *L02E11 et l'amplitude de la raie B est proportionnelle àlE;01E01, la raie C correspondant à la contribution statique au spectre des battements. Les trois contributions A, B et C se trouvent à l'intérieur de la bande passante temporelle os du détecteur, les autres contributions de l'équation (3) ne sont pas détectées. Pour z La déposante a montré expérimentalement la faisabilité d'une cartographie spatiale de l'amplitude de vibration et de la phase de vibration, avec une excellente précision. Le montage expérimental est représenté sur la figure 5. Pour cette expérience, l'objet utilisé est un lamellophone à lamelles métalliques dont les lamelles présentent des fréquences de résonance différentes les unes des autres, fonction notamment de la masse et de la longueur de chaque lamelle. Un actuateur piézoélectrique permet d'exciter les lamelles du lamellophone à une fréquence de vibration variable o. Les éléments utilisés pour le montage expérimental de la figure 5 peuvent être sensiblement identiques à ceux décrits en référence à la figure 1. Dans cet exemple cependant, pour gagner en rapport signal sur bruit, l'onde de référence WLO (oscillateur local) présente trois fréquences de décalage (01, c02 et c03 décalées de la fréquence d'émission coL respectivement d'une quantité ôi, ô2 et ô3, les décalages en fréquence étant choisis pour permettre une démodulation des interférogrammes dans la bande spectrale temporelle os du détecteur, par exemple une caméra CCD. Par exemple, ô. = (n - 2)(o - cosy., avec os la fréquence d'acquisition du détecteur, n = 1, 2, 3 et - 0,5 En faisant varier la fréquence d'excitation co, on met successivement en vibration chaque lamelle du lamellophone à sa fréquence de résonance. Dans cette expérience, la fréquence d'excitation v = co/(27c) de l'actuateur piézoélectrique varie de 0 à 3 kHz, par pas de 1 Hz. 3000 séquences de N= 8 interférogrammes sont acquises à une fréquence cos/(27c) = 20 Hz, et traitées selon la méthode décrite ci-dessus. Les décalages en fréquence ygos, y2c0s et y3c0s sont respectivement de - 5 Hz, - 2,5 Hz et - 7,5 Hz. Les figures 6A à 6D montrent les composantes fréquentielles des hologrammes fi(co, ) telles que définies par l'équation (4) ci-dessus et calculées sur l'ensemble des points de l'objet. Plus précisément, la figure 6A montre la composante non décalée li/(71cos )1 de l'hologramme du lamellophone, proportionnelle à l'amplitude de la composante non décalée Iso de l'onde signal. Les figures 6B à 6D représentent la composante fréquentielle V--1(72cos) de l'hologramme du lamellophone proportionnelle à l'amplitude de la bande latérale 1E11 de l'onde signal, obtenues respectivement à trois fréquences d'excitation co différentes, correspondant respectivement aux fréquences de résonance des lamelles indiquées par une flèche sur chacune des figures. Les fréquences d'excitation y = co/(27c) sont respectivement 541 Hz, 1006 Hz et 2211 Hz pour les figures 6B, 6C et 6D. Une observation qualitative des fréquences de résonance des différentes lamelles est ainsi mise en évidence, en faisant varier la fréquence d'excitation du lamellophone et en adaptant une fréquence de décalage de l'oscillateur locale à la fréquence d'excitation. Comme expliqué précédemment dans le cas de l'oscillateur local dual, une mesure quantitative de l'amplitude de vibration z et du retard de phase y en un point de l'objet peut alors être faite à partir des amplitudes complexes des différentes raies du spectre en fréquence des hologrammes mesurées simultanément grâce aux fréquences de décalage multiples de l'oscillateur local. Dans l'exemple de la figure 5, l'oscillateur local présente 3 fréquences de décalage et l'amplitude de vibration z ainsi que le retard de phase y peuvent se déduire de deux des trois raies 1H(71cos H(72w) et V-I(73cos) du spectre en fréquence des hologrammes, ou des trois raies du spectre en fréquence des hologrammes, par exemple par les équations : ( Ii(72cos) fi(ricos) fi(73cos) fi(ricos) z 47z- (8) 1 arg(fi_072c0s 1107 co ) + arg (11-(r1cos) (9) "lo s ) en tenant compte des équations (1), (2) et (5). La figure 7 montre l'amplitude de vibration de la plus longue lamelle du lamellophone (indiquée par la flèche dans la figure 6B) en fonction de la fréquence d'excitation y. A la fréquence de résonance de 541 Hz de cette lamelle, l'amplitude de vibration 50 est maximale. Selon une variante, pour verrouiller la phase constante wo d'une mesure à l'autre et déterminer avec précision le retard de phase mécanique y, une variation temporelle continue de la fréquence d'excitation (« chirp » selon l'expression anglo-saxonne) est appliquée simultanément au signal d'excitation de l'objet et à l'un des décalages en fréquence Ô. de l'onde de référence. La variation temporelle d'excitation est par exemple linéaire et s'écrit : co(t) = coi + (coF - coi)t/T, (10) où coi est la fréquence d'excitation inférieure, (OF est la fréquence d'excitation supérieure, et T la durée de la variation. Les figures 8A et 8B montrent des résultats expérimentaux obtenus avec un montage expérimental selon l'exemple de la figure 5, tel que décrit ci-dessus. La fréquence d'excitation y varie continûment entre 536 Hz et 546 Hz pendant un temps de variation T= 99,2 s. Cette même variation de fréquence est appliquée à l'une des composantes de l'oscillateur local. A chaque fréquence co(t), N= 8 interférogrammes ont été acquis à une cadence de 20 Hz par la caméra CCD, et traités selon la méthode de la présente invention. L'amplitude de vibration de la lamelle la plus longue du lamellophone en fonction de la fréquence d'excitation y est montrée sur la figure 8A, où l'amplitude de vibration est maximale à la fréquence de résonance de la lamelle considérée. La figure 8B représente le déphasage y - wo en fonction de la fréquence d'excitation. La courbe met en évidence un saut de phase de 7C à la fréquence de résonance. Bien que décrite essentiellement à travers des exemples mettant en oeuvre une excitation forcée de l'objet, la présente invention s'applique également à l'étude de la vibration et plus particulièrement des fréquences de résonance d'objets présentant des spectres de vibration, par exemple dans le cadre de contrôle non destructif Les mesures précises de l'amplitude de vibration et du retard de phase telles que décrites précédemment pourront permettre d'analyser avec précision notamment les résonances de l'objet. Bien que décrite à travers un certain nombre d'exemples de réalisation détaillés, la méthode et le dispositif d'holographie hétérodyne hors-axe selon l'invention comprennent différentes variantes, modifications et perfectionnements qui apparaîtront de façon évidente à l'homme de l'art, étant entendu que ces différentes variantes, modifications et perfectionnements font partie de la portée de l'invention, telle que définie par les revendications qui suivent. REVENDICATIONS1. Méthode d'holographie numérique pour détecter l'amplitude z de vibration d'un objet (15) à au moins une fréquence de vibration w, comprenant : - la génération d'une onde laser de fréquence coL d'illumination de l'objet (W1) et 5 d'une onde de référence (WL0) cohérente avec l'onde d'illumination (Wi); - l'acquisition au moyen d'un détecteur optoélectronique bidimensionnel (19) de bande passante temporelle oh d'une série d'interférogrammes résultant de l'interférence de l'onde de référence (Wu)) et d'une onde signal (Ws) issue de l'objet (15), les ondes signal et de référence (Ws, Wu)) présentant un décalage angulaire (0), le champ de l'onde de 10 référence comprenant une première composante (ELoi) de fréquence w1 et au moins une deuxième composante (ELo2), de fréquence w2, les fréquences col et co2 étant décalées de la fréquence coL de l'onde laser respectivement d'une quantité déterminée 81, 82, avec 81 = Yiws et 82 = + y2ws, où q entier relatif et - 0,5 2. Méthode selon la revendication 1, dans laquelle les normes de yi et y2 sont différentes. 3. Méthode selon l'une des revendications 1 et 2, comprenant : - l'émission d'une onde laser (WEm) de fréquence coL ; etFR 12 60883 19 27658/028FR1 - la séparation de l'onde d'émission laser (WEm) pour former l'onde d'illumination de l'objet (WI) et une deuxième onde (WL) envoyée vers un dispositif optoélectronique de décalage de fréquence optique (13, 13') pour former l'onde de référence (Ww). 4. Méthode selon l'une quelconque des revendications précédentes, comprenant l'excitation forcée de l'objet (15) à ladite fréquence co. 5. Méthode selon la revendication 4, dans laquelle la fréquence d'excitation de l'objet est variable temporellement de façon continue entre une fréquence d'excitation inférieure coi et une fréquence d'excitation supérieure coF pendant une durée (T) donnée, et dans laquelle l'acquisition et le traitement des interférogrammes sont effectuées pendant ladite durée. 6. Méthode selon l'une des revendications 4 ou 5, comprenant en outre le calcul du retard de phase entre la phase de vibration de l'objet et la phase du signal d'excitation à partir du rapport des amplitudes de la première raie (A, B) et de la deuxième raie (A, B). 7. Méthode selon l'une quelconque des revendications 1 à 3, dans laquelle : - la vibration de l'objet (15) présente un spectre de fréquences ; et - l'acquisition et le traitement des interférogrammes sont effectuées pour obtenir l'amplitude z de vibration de l'objet à au moins une fréquence de vibration co contenu dans le spectre. 8. Méthode selon l'une quelconque des revendications précédentes, dans laquelle le champ de l'onde de référence (WL0) comprend au moins une troisième composante EL03 de fréquence w3, la fréquence w3 étant décalée de la fréquence coL de l'onde laser d'une quantité déterminée 83 = y3w, + qw, où q entier relatif et - 0,5 9. Méthode selon l'une quelconque des revendications précédentes, dans laquelle q = 1. 10. Dispositif d'holographie numérique pour détecter l'amplitude z de vibration d'un objet (15) à une fréquence de vibration co, comprenant : - un détecteur optoélectronique bidimensionnel (19) de bande passante temporelle (os;FR 12 60883 20 27658/028FR1 - des moyens de génératidn de deux ondes laser (3, 7) cohérentes entre elles, une onde laser d'illumination de l'objet (WI) de fréquence coL et une onde de référence (Ww); le champ de l'onde de référence (Ww) présentant au moins une première et une deuxième composantes (Ewi, ELo2) de fréquences (ni, co2 décalées de la fréquence coL de l'onde laser respectivement d'une quantité déterminée 81; 82, avec 81 = yi cos et 82 = qco + y20)' où q entier relatif et - 0,5 11. Dispositif d'holographie numérique selon la revendication 10, comprenant - une source d'émission laser (3) à la fréquence (DL; - des moyens de séparation (7) de l'onde d'émission (WEm) laser pour former l'onde d'illumination de l'objet (W1) et une seconde onde (WL) ; - un dispositif optoélectronique de décalage de fréquence optique (13, 13'), destiné à 30 recevoir ladite seconde onde (WL) pour former l'onde de référence (Ww). 12. Dispositif d'holographie numérique selon la revendication 11, dans lequel le dispositif optoélectronique de décalage de fréquence optique (13, 13') comprend deuxa . FR 12 60883 21 27658/028FR1 modulateurs acousto-optiques tràvaillant à. des fréquences prédéterminées sur des ordres de diffraction opposés pour former ladite onde de référence (Ww). 13. Dispositif holographique selon la revendication 11, dans lequel le dispositif optoélectronique de décalage de fréquence optique (13) comprend une cellule de Pockels. 14. Dispositif holographique selon l'une des revendications 11 à 13, dans lequel les moyens de séparation (7) de l'onde laser sont fibrés. 15. Dispositif holographique selon l'une quelconque des revendications 10 à 14, comprenant en outre un moyen d'excitation de l'objet (15) à la fréquence co de vibration donnée. 16. Dispositif holographique selon l'une quelconque des revendications 10 à 15, dans lequel le détecteur optoélectronique bidimensionnel (19) est une caméra de type CCD ou CMOS.15