Brevet FR3071619A1 — publié le 20190329

PROCÉDÉ DE DETERMINATION DE PARAMETRES CARACTERISTIQUES D'UN OSCILLATEUR DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne un procédé de détermination de paramètres caractéristiques d'un oscillateur, tels que la fréquence de résonance et le facteur de qualité. ARRIÈRE-PLAN TECHNOLOGIQUE DE L'INVENTION Un oscillateur à actionnement électrostatique de type microsystème ou 10 nanosystème électromécanique, respectivement MEMS (« microelectromechanical system » en anglais) ou NEMS (« nanoelectromechanical system » en anglais), est un dispositif comprenant un élément mobile pouvant être mis en oscillation. De tels oscillateurs consomment peu d'énergie et présentent un encombrement réduit. Ils trouvent donc une application particulièrement avantageuse 15 dans des objets nomades tels que des téléphones intelligents (smartphones) ou des tablettes numériques. Ces oscillateurs sont en particulier utilisés pour fabriquer des capteurs inertiels, comme des accéléromètres ou des gyromètres, destinés à équiper de tels objets. Un objet nomade peut par ailleurs être équipé d'un système de localisation par 20 satellite (GPS pour « Global Positionning System » en anglais) qui consomme beaucoup d'énergie. Les capteurs inertiels permettent de calculer en temps réel les déplacements de l'objet et donc de réduire le recours au système de localisation par satellite. Ce dernier reste cependant utile pour déterminer une position de référence de l'objet, par exemple à des intervalles de temps réguliers. 25 L'oscillateur à actionnement électrostatique comporte des moyens d'actionnement qui transforment un signal d'excitation appliqué en entrée de l'oscillateur en une force d'actionnement électrostatique, également appelée « force d'excitation », agissant sur l'élément mobile. Les déplacements de l'élément mobile sont mesurés par des moyens de détections qui génèrent un signal de 2 réponse en sortie de l'oscillateur. Le signal d'excitation et le signal de réponse se présentent par exemple sous la forme de tensions. L'oscillateur est caractérisé par deux paramètres essentiels, à savoir une fréquence de résonance et un facteur de qualité. En effet, les oscillations de 5 l'élément mobile sont régies par des équations de mouvement qui sont modifiées lorsqu'une force externe, par exemple due à une accélération ou à un effet gyroscopique, s'exerce sur l'élément mobile. La fréquence de résonance est alors également modifiée. En connaissant la fréquence de résonance initiale de l'oscillateur, il est donc possible de quantifier la force externe appliquée. De plus, il est préférable d'exciter l'oscillateur à une fréquence proche de la fréquence de résonance de manière à maximiser l'amplitude des oscillations. La force s'exerçant sur l'élément mobile peut ainsi être détectée plus facilement, et la sensibilité du capteur est améliorée. Le facteur de qualité de l'oscillateur a quant à lui une influence sur la précision du capteur inertiel. Plus le facteur de qualité est élevé, plus le capteur est précis et plus il est possible par exemple de limiter l'utilisation du système de localisation par satellite. Il est donc important de mesurer la fréquence de résonance et le facteur de qualité de l'oscillateur, notamment pour valider le développement et la fabrication du capteur inertiel. Une première méthode de mesure de la fréquence de résonance et du facteur de qualité consiste à exciter l'oscillateur au moyen d'une tension d'excitation sinusoïdale de fréquence Fo et à mesurer l'amplitude des oscillations à la fréquence 2F0. L'amplitude des oscillations est mesurée à la fréquence 2F0 car la force d'actionnement électrostatique est proportionnelle au carré de la tension d'excitation. En balayant une plage de fréquence, on obtient une densité spectrale d'amplitude des oscillations à partir de laquelle la fréquence de résonance peut être déterminée. En effet, la densité spectrale d'amplitude forme un pic de résonance présentant un maximum d'amplitude atteint à la fréquence de résonance. Le facteur de qualité peut quant à lui être déterminé à partir de la largeur du pic de résonance considérée à la moitié de sa hauteur. 3 Un inconvénient de cette méthode de mesure par balayage en fréquence est que pour mesurer correctement l'amplitude des oscillations, l'oscillateur doit retrouver une position de repos entre deux mesures successives. Le temps nécessaire pour que l'oscillateur retrouve sa position de repos après excitation est 5 approximativement égal à trois fois une constante d'amortissement qui est proportionnelle au facteur de qualité. En considérant par exemple un facteur de qualité de l'ordre de 106, la constante d'amortissement est de l'ordre de 10 secondes. Dans ce cas, il faut par exemple plus de 2 heures pour acquérir 250 points de mesure. 10 Par ailleurs, la précision avec laquelle le pic de résonance est défini dépend du nombre de points de mesure. Si le nombre de points n'est pas suffisant, le pic de résonance peut être élargi artificiellement et son maximum d'amplitude peut être décalé en fréquence, ce qui fausse la mesure de la fréquence de résonance et du facteur de qualité. 15 Le temps nécessaire pour mettre en oeuvre cette méthode de mesure est donc un facteur limitant la productivité, d'autant plus que la tendance actuelle est à l'augmentation du facteur de qualité des oscillateurs. Ce temps de mesure ne permet raisonnablement pas, en particulier d'un point de vue industriel, de caractériser tous les oscillateurs d'une plaque de silicium qui en comporte généralement plus d'une 20 centaine. Seuls quelques oscillateurs sont alors caractérisés, ce qui rend difficile d'estimer un rendement de fabrication des oscillateurs. Une deuxième méthode de mesure, dite « par impulsion à bande large », de la fréquence de résonance et du facteur de qualité consiste à exciter l'oscillateur au moyen d'une impulsion de tension et à acquérir le signal de réponse généré par 25 l'oscillateur. En effet, une impulsion dans le domaine temporel correspond à une amplitude constante dans le domaine fréquentielle. En pratique, l'impulsion présente plutôt une largeur de bande spectrale élevée. Par conséquent, en calculant la transformée de Fourier du signal acquis, il est possible de localiser le pic de résonance présent dans la bande spectrale de l'impulsion. Il s'agit d'une méthode de 30 mesure rapide puisque l'acquisition peut durer seulement quelques secondes. 4 Selon cette méthode, la durée d'acquisition détermine directement la résolution fréquentielle de la mesure. Par exemple, une durée d'acquisition de 1 seconde donne une résolution de 1 hertz, et une durée d'acquisition de 100 secondes donne une résolution de 0,01 hertz. 5 Un inconvénient de cette méthode de mesure par impulsion à bande large est que l'oscillateur à actionnement électrostatique comporte une capacité parasite entre les moyens d'actionnement et les moyens de détection. Il y a donc un couplage électrique direct de la tension d'excitation vers la sortie de l'oscillateur. La tension d'excitation étant à bande large, sa densité spectrale d'amplitude se superpose alors i.o à la densité spectrale d'amplitude du signal de réponse. Si le couplage capacitif est suffisamment important, une faible oscillation de l'oscillateur est alors noyée dans le spectre de la tension d'excitation, ce qui généralement rend cette méthode de mesure inutilisable. Pour pallier à ce problème, à la place d'une impulsion à bande large on peut 15 utiliser une impulsion à bande étroite pour exciter l'oscillateur. Cette impulsion à bande étroite a par exemple la forme d'une fonction sinus cardinal de largeur spectrale AF modulé par une fonction sinus de fréquence Fo. La tension d'excitation u(t) peut alors être définie par l'équation suivante : u(t) = sinc(AF. t). sin(2. Tr. Fo. t) (1) 20 La force d'excitation est proportionnelle au carré de la tension d'excitation u(t) et peut être exprimée par l'équation suivante : F(t) = k. [sinc(AF. t). sin(2. Tr. Fo. t)]2 (2) où k est un facteur de proportionnalité. L'équation (2) peut être développée pour obtenir l'équation suivante : 25 F(t) = k [î . sinc2 (AF. t) - î . cos(4. Tt. Fo. t) . sinc2 (AF. t)1 (2') La force d'excitation est donc la somme d'un premier terme proportionnel au carré du sinus cardinal et d'un deuxième terme proportionnel au carré du sinus cardinal modulé à la fréquence 2F0. 5 La figure 1 montre la densité spectrale d'amplitude (DSA) 110 de la tension d'excitation u(t). Cette densité spectrale 110 est une fonction porte de largeur AF centrée sur la fréquence Fo. La figure 1 montre également la densité spectrale d'amplitude de la force d'excitation F(t) qui comprend une première composante 121 5 et une deuxième composante 122 correspondant respectivement au premier terme et au deuxième terme de l'équation (2'). Ces composantes spectrales 121, 122 ont chacune une forme triangulaire de largeur 2AF, la première composante 121 étant centrée sur 0 Hz et la deuxième composante 122 étant centrée sur 2F0. La fréquence de résonance de l'oscillateur est recherchée dans une zone qui 10 correspond à la plage de fréquence sur laquelle s'étend la deuxième composante 122. En définissant de manière appropriée la fréquence Fo du sinus et la largeur spectrale AF du sinus cardinal, la densité spectrale d'amplitude 110 de l'impulsion excitatrice transmise par couplage capacitif ne se superpose pas avec la densité spectrale d'amplitude de la force d'actionnement électrostatique, comme illustré à la 15 figure 1. Ainsi, la densité spectrale d'amplitude résultant du couplage capacitif de l'impulsion excitatrice n'est pas présente dans la zone de recherche de la fréquence de résonance. Un inconvénient de cette méthode de mesure par impulsion à bande étroite est que la densité spectrale d'amplitude de la force d'excitation n'est pas constante. 20 Elle est en effet maximale à la fréquence 2F0 et décroît linéairement de part et d'autre jusqu'à s'annuler à la fréquence 2F0 - 2AF et à la fréquence 2F0 + 2AF. L'oscillateur est donc excité avec une force qui dépend de l'écart entre la fréquence de résonance et la fréquence 2F0. Ce phénomène est mis en évidence sur les figures 2A et 2B. 25 La figure 2A montre trois exemples de densités spectrales d'amplitude 211, 212, 213 de la force d'excitation lorsque la fréquence Fo est égale respectivement à 16,25 kHz, 16,5 kHz et 17,5 kHz. La fréquence de résonance de l'oscillateur est recherchée dans la fenêtre fréquentielle formée par chacune de ces densités spectrales. A chacune des densités spectrales 211, 212, 213 de la force d'excitation 30 correspond une densité spectrale d'amplitude 221, 222, 223 du signal de réponse acquis en sortie de l'oscillateur. 6 La figure 2B est un agrandissement de la figure 2A autour du pic de résonance. On constate que plus ce pic est proche du bord de la fenêtre fréquentielle, plus son amplitude est faible. La forme de la densité spectrale d'amplitude de la force électrostatique peut donc conduire à considérer que 5 l'oscillateur n'est pas fonctionnel alors que la force d'excitation est seulement trop faible pour mettre en oscillation l'élément mobile de l'oscillateur. RÉSUMÉ DE L'INVENTION Il ressort de ce qui précède qu'il existe un besoin de disposer d'une méthode à la fois rapide et fiable pour déterminer les paramètres caractéristiques d'un 10 oscillateur à actionnement électrostatique. Pour répondre à ce besoin, un premier aspect de l'invention concerne un procédé de détermination de paramètres caractéristiques d'un oscillateur à actionnement électrostatique, le procédé comportant les étapes suivantes : - générer une première tension d'excitation définie comme étant la somme 15 d'une première tension sinusoïdale et d'une impulsion de tension ; - appliquer la première tension d'excitation en entrée de l'oscillateur ; - acquérir dans le domaine temporel une première tension de réponse présente en sortie de l'oscillateur lorsque la première tension d'excitation est appliquée en entrée de l'oscillateur ; 20 - obtenir, par transformation dans le domaine fréquentiel, une première densité spectrale d'amplitude de la première tension de réponse ; - déterminer les paramètres caractéristiques de l'oscillateur à partir de la première densité spectrale d'amplitude. Grâce à l'invention, la caractérisation de l'oscillateur s'effectue de manière 25 rapide en obtenant une densité spectrale d'oscillation complète en effectuant une seule acquisition du signal de réponse, cette acquisition durant quelques secondes. Il n'est dans ce cas pas nécessaire d'exciter l'oscillateur pour chaque fréquence de la densité spectrale, ni d'attendre que l'oscillateur retrouve une position de repos entre deux points de la densité spectrale. De plus, l'utilisation de la somme d'une tension 30 sinusoïdale et d'une impulsion de tension pour exciter l'oscillateur permet d'obtenir 7 une force d'actionnement électrostatique présentant une densité spectrale d'amplitude constante. Les paramètres caractéristiques de l'oscillateur comprennent une fréquence de résonance. Selon un mode de mise en oeuvre, le procédé comportant les étapes 5 initiales suivantes : - estimer la fréquence de résonance de l'oscillateur ; - définir une plage de fréquence, dite « zone de recherche », dans laquelle la fréquence de résonance est recherchée, la zone de recherche étant centrée autour d'une fréquence médiane égale à la fréquence de résonance estimée, 10 la zone de recherche présentant une largeur prédéterminée. Un avantage est de savoir sur quelle plage de fréquence on souhaite chercher la fréquence de résonance afin de paramétrer la tension d'excitation convenablement. Selon un mode de mise en oeuvre, l'impulsion de tension est une réponse impulsionnelle finie d'un filtre passe-bas. 15 Selon un mode de mise en oeuvre, l'impulsion de tension présente une largeur spectrale égale la largeur de la zone de recherche, la première tension sinusoïdale présentant une première fréquence égale à la fréquence médiane. Un avantage est de faire coïncider la densité spectrale d'amplitude de la force électrostatique avec la zone de recherche de la fréquence de résonance. 20 Selon un mode de mise en oeuvre, l'impulsion de tension se présente sous la forme d'une fonction sinus cardinal. Selon un mode de mise en oeuvre, l'impulsion de tension est une réponse impulsionnelle finie d'un filtre passe-bande, le procédé comportant également les étapes suivantes : 25 - générer une deuxième tension d'excitation définie comme étant la somme d'une deuxième tension sinusoïdale et de l'impulsion de tension ; - appliquer la deuxième tension d'excitation en entrée de l'oscillateur ; - acquérir dans le domaine temporel une deuxième tension de réponse présente en sortie de l'oscillateur lorsque la deuxième tension d'excitation est 30 appliquée en entrée de l'oscillateur ; 8 - obtenir, par transformation dans le domaine fréquentiel, une deuxième densité spectrale d'amplitude de la deuxième tension de réponse ; - reconstituer une densité spectrale d'amplitude à partir de la première densité spectrale d'amplitude et de la deuxième densité spectrale d'amplitude ; 5 les paramètres caractéristiques de l'oscillateur étant déterminés à partir de la densité spectrale d'amplitude reconstituée. Un avantage est de supprimer la densité spectrale d'amplitude de la tension sinusoïdale transmis par couplage capacitif à la densité spectrale d'amplitude de la tension de réponse. Selon un mode de mise en oeuvre, l'impulsion de tension se présente sous la 10 forme d'une fonction sinus cardinal modulée par une fonction sinus. Selon un mode de mise en oeuvre, la fonction sinus cardinal présente une largeur spectrale, la fonction sinus présentant une fréquence de modulation égale à la largeur spectrale de la fonction sinus cardinal, la première tension sinusoïdale présentant une première fréquence et la deuxième tension sinusoïdale présentant 15 une deuxième fréquence égale à la somme de la première fréquence et de la largeur spectrale de la fonction sinus cardinal. Un avantage est d'obtenir une densité spectrale reconstituée continue. Selon un mode de mise en oeuvre, la largeur spectrale de la fonction sinus cardinal est égale à un quart de la largeur de la zone de recherche, la fréquence 20 médiane étant égale à la moitié de la somme de la première fréquence et de la deuxième fréquence. Un deuxième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions pour la mise en oeuvre du procédé selon le premier aspect de l'invention, lorsque ces instructions sont exécutées par un 25 processeur. BRÈVE DESCRIPTION DES FIGURES L'invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent, parmi lesquelles : 9 la figure 1 représente schématiquement la densité spectrale d'amplitude d'une tension d'excitation appliquée en entrée d'un oscillateur et la densité spectrale d'amplitude de la force électrostatique qui en résulte, dans le cas d'une mesure par impulsion à bande étroite selon l'art antérieur ; 5 la figure 2A représente schématiquement trois densités spectrales d'amplitude de la force électrostatique agissant sur un oscillateur excité au moyen d'un sinus cardinal modulé par un sinus, ainsi que les trois densités spectrales d'amplitude de la tension de réponse correspondante ; la figure 2B est une vue agrandie de la figure 2A autour des pics de 10 résonances ; la figure 3 représente schématiquement la densité spectrale d'amplitude d'une tension d'excitation appliquée en entrée d'un oscillateur, selon un premier mode de mise en oeuvre du procédé de l'invention ; la figure 4 représente schématiquement la densité spectrale d'amplitude de la 15 force d'excitation produite par la tension d'excitation de la figure 3 ; la figure 5 est un exemple de densité spectrale d'amplitude du signal de réponse d'un oscillateur de gyromètre, obtenu avec le premier mode mise en oeuvre du procédé de l'invention ; la figure 6 représente schématiquement les densités spectrales d'amplitude 20 d'une première tension d'excitation et d'une deuxième tension d'excitation, selon un deuxième mode de mise en oeuvre du procédé de l'invention ; la figure 7 représente schématiquement les densités spectrales d'amplitude de la force d'excitation produite par chacune des tensions d'excitation de la figure 6 ; 25 la figure 8 est un exemple de densité spectrale d'amplitude du signal de réponse de l'oscillateur de la figure 5, obtenu avec le deuxième mode mise en oeuvre du procédé de l'invention ; la figure 9A représente schématiquement trois densités spectrales d'amplitude de la force électrostatique agissant sur l'oscillateur de la figure 2A, ainsi que 30 les trois densités spectrales d'amplitude de la tension de réponse correspondante, selon le deuxième mode de mise en oeuvre du procédé de l'invention ; 10 la figure 9B est une vue agrandie de la figure 9A autour des pics de résonances ; la figure 10 illustre une comparaison entre une mesure par balayage de fréquence selon l'art antérieur et une mesure selon le procédé de l'invention. 5 Les figures ne sont présentées qu'à titre indicatif et nullement limitatif de l'invention. Pour plus de clarté, les éléments identiques ou similaires sont repérés par des signes de référence identiques sur toutes les figures. DESCRIPTION DÉTAILLÉE DE MODES DE MISE EN OEUVRE DE L'INVENTION 10 Un oscillateur de type MEMS ou NEMS comporte une partie fixe et un élément mobile par rapport à la partie fixe. L'oscillateur comporte également une entrée et une sortie entre lesquelles est disposé l'élément mobile. L'entrée de l'oscillateur est pourvue de moyens d'actionnement électrostatiques configurés pour mettre en vibration l'élément mobile. Les moyens 15 d'actionnement transforment un signal électrique d'entrée, également appelé « signal d'excitation », appliqué en entrée de l'oscillateur en une force d'actionnement électrostatique, également appelée « force d'excitation », agissant sur l'élément mobile. Le signal d'excitation se présente par exemple sous la forme d'une tension. La force d'excitation entraîne des déplacements oscillatoires de l'élément mobile. Les 20 moyens d'actionnement comprennent par exemple une première électrode fixe appartenant à la partie fixe et une première électrode mobile appartenant à l'élément mobile. Ces électrodes peuvent chacune avoir la forme d'un peigne et former ensemble une structure à peignes capacitifs interdigités. La sortie de l'oscillateur est pourvue de moyens de détection des 25 déplacements de l'élément mobile. Les moyens de détections sont configurés pour délivrer un signal électrique de sortie, également appelé « signal de réponse », qui dépend des déplacements détectés de l'élément mobile. Par exemple, les moyens de détection sont de type capacitif et comportent des moyens de mesure de capacité entre une deuxième électrode fixe appartenant à la partie fixe et une deuxième 30 électrode mobile appartenant à l'élément mobile. Ces électrodes peuvent chacune n avoir la forme d'un peigne et former ensemble une structure à peignes capacitifs interdigités. L'oscillateur est caractérisé par une fréquence de résonance et un facteur de qualité. Ces deux paramètres sont essentiels pour valider la conception de 5 l'oscillateur et pour son intégration dans un dispositif tel qu'un capteur inertiel, par exemple un accéléromètre ou un gyromètre. Lorsque l'oscillateur est excité à la fréquence de résonance, l'amplitude des oscillations de l'élément mobile est maximale. Lorsque le signal d'excitation n'est plus appliqué, les oscillations de l'élément mobile s'amortissent jusqu'à ce que 10 l'élément mobile retrouve une position de repos dans laquelle l'élément mobile cesse de se mouvoir. L'amortissement dans le temps t des oscillations dépend d'un terme exponentiel e-, où T est une constante d'amortissement, ou temps de relaxation, définissant la vitesse avec laquelle l'amplitude des oscillations décroît après excitation. 15 Le facteur de qualité Q de l'oscillateur est défini par la relation suivante : Q= Tr. FR. T OÙ FR est la fréquence de résonance de l'oscillateur. Le facteur de qualité traduit notamment la précision du capteur dans lequel l'oscillateur est intégré. Plus le facteur de qualité est élevé, plus la précision l'est aussi. La présente invention a pour objectif principal de proposer un procédé de 20 détermination des paramètres caractéristiques de l'oscillateur, ce procédé étant à la fois rapide et fiable. Le procédé selon l'invention vise en particulier à déterminer la fréquence de résonance et le facteur de qualité de l'oscillateur. Dans ce procédé, l'oscillateur est excité par une tension d'excitation qui résulte de la somme d'une tension sinusoïdale et d'une impulsion de tension. 25 Un premier mode de mise en oeuvre du procédé selon l'invention va maintenant être décrit, en référence aux figures 3 à 5. Le procédé comporte une étape de génération de la tension d'excitation. La tension d'excitation est par exemple fournie par un générateur de signal connecté électriquement à l'entrée de 12 l'oscillateur. Selon le premier mode de mise en oeuvre, la tension d'excitation est la somme d'une tension sinusoïdale de fréquence Fo constante et d'une réponse impulsionnelle finie d'un filtre passe-bas. La réponse impulsionnelle a de préférence la forme d'une fonction sinus cardinal de largeur spectrale AF, la tension d'excitation 5 u(t) pouvant alors être définie par l'équation suivante : u(t) = A. sinc(AF. t) + B. sin(2.Tr. Fo. t) (3) où A et B sont l'amplitude respectivement du sinus cardinal et de la tension sinusoïdale. Alternativement, la réponse impulsionnelle peut être celle de tout type de filtre 10 passe-bas tel qu'une fenêtre de Hamming, de Hanning, de Tchebychev, de Bartlett, ou de Blackman. La figure 3 représente schématiquement la densité spectrale d'amplitude 300 de la tension d'excitation u(t), par exemple calculée au moyen d'une transformée de Fourier. Une densité spectrale d'amplitude (DSA) est exprimée en unité par racine de 15 hertz, c'est-à-dire en volt par racine de hertz (V.Hz-°5) lorsqu'il s'agit de la densité spectrale d'amplitude d'une tension d'excitation, et en newton par racine de hertz (N.Hz-°'5) lorsqu'il s'agit de la densité spectrale d'amplitude d'une force d'excitation. La densité spectrale d'amplitude 300 de la tension d'excitation comporte : 20 - une première composante 301 correspondant à la densité spectrale du sinus cardinal, la première composante 301 étant une fonction porte centrée sur 0 Hz et de largeur spectrale AF ; - une deuxième composante 302 correspondant à la densité spectrale de la tension sinusoïdale, la deuxième composante 302 étant un pic de Dirac situé 25 à la fréquence Fo. La tension d'excitation est appliquée en entrée de l'oscillateur pour le mettre en oscillation. La force d'excitation qui en résulte est proportionnelle au carré de la tension d'excitation. La force d'excitation F(t) peut ainsi être exprimée par l'équation suivante : 13 F(t) = k [A. sinc2(AF. t) + B2 - B27. cos(4. Tt. Fo. t) + 2. A. B. sinc(AF. t). sin(2. n. Fo. t)1 (4) où k est un coefficient de proportionnalité. La force d'excitation est donc la somme 5 d'un premier terme proportionnel au carré du sinus cardinal, d'un deuxième terme égal à une constante, d'un troisième terme proportionnel à une fonction sinusoïdale de fréquence 2F0, et d'un quatrième terme proportionnel au sinus cardinal modulé à la fréquence Fo. La figure 4 représente schématiquement la densité spectrale d'amplitude 400 10 de la force d'excitation F(t) obtenue par transformation dans le domaine fréquentiel, par exemple au moyen d'une transformée de Fourier. La densité spectrale d'amplitude 400 comprend : - une première composante 401 correspondant à la densité spectrale d'amplitude du premier terme de l'équation (4), la première composante 401 15 ayant une forme triangulaire centrée sur 0 Hz et de largeur spectrale 2AF ; - une deuxième composante 402 correspondant à la densité spectrale d'amplitude du deuxième terme de l'équation (4), la deuxième composante 402 étant un pic de Dirac situé à 0 Hz ; - une troisième composante 403 correspondant à la densité spectrale 20 d'amplitude du troisième terme de l'équation (4), la troisième composante 403 étant un pic de Dirac situé à la fréquence 2F0 ; - une quatrième composante 404 correspondant à la densité spectrale d'amplitude du quatrième terme de l'équation (4), la quatrième composante 404 étant une fonction porte centrée sur la fréquence Fo et de largeur 25 spectrale AF. La fréquence de résonance de l'oscillateur est recherchée dans une plage de fréquence dite « zone de recherche » définie par la fréquence Fo de la tension sinusoïdale et par la largeur spectrale AF de la fonction sinus cardinal. La zone de recherche s'étend entre une fréquence minimale égale à Fo - '2 et une fréquence 30 maximale égale à Fo + '.2 . 14 Initialement, il est possible d'estimer la fréquence de résonance de l'oscillateur, en se basant notamment sur sa conception. La fréquence Fo de la tension sinusoïdale peut alors être définie de sorte que la fréquence de résonance estimée se trouve dans la zone de recherche, et en particulier que la zone de 5 recherche soit centrée sur la fréquence de résonance estimée. De même, la largeur spectrale AF de la fonction sinus cardinal peut être définie pour adapter la largeur de la zone de recherche. La tension d'excitation appliquée en entrée de l'oscillateur est transmise en sortie de l'oscillateur par couplage capacitif dû à une capacité parasite présente 10 entre l'entrée et la sortie de l'oscillateur. La largeur spectrale AF de la réponse impulsionnelle est de préférence inférieure à FO-AF/2 de manière à éviter que ce signal parasite, qui a une amplitude plus élevée, ne recouvre le signal issu de la résonance. Lorsque la tension d'excitation est appliquée à l'oscillateur, on procède à une 15 étape d'acquisition dans le domaine temporel du signal de réponse de l'oscillateur fourni par les moyens de détection. Pour ce faire, des moyens d'acquisitions sont connectés électriquement à la sortie de l'oscillateur. Les moyens d'acquisition peuvent être une carte électronique d'acquisition ou un oscilloscope numérique. Le signal de réponse se présente par exemple sous la forme d'une tension. 20 L'acquisition est effectuée avec une fréquence d'échantillonnage configurée pour respecter le critère de Shannon. La fréquence d'échantillonnage dépend ici de la fréquence Fo de la tension sinusoïdale et de la largeur spectrale AF de la réponse impulsionnelle, la fréquence d'échantillonnage Fe devant respecter l'inégalité suivante : AF) Fe 2. (F0 + 2) 25 La tension de réponse acquise est ensuite transformée dans le domaine fréquentiel pour obtenir sa densité spectrale d'amplitude. La densité spectrale d'amplitude de la tension de réponse est par exemple calculée au moyen d'une transformée de Fourier. A partir de cette densité spectrale d'amplitude, il est possible de déterminer la fréquence de résonance et le facteur de qualité de l'oscillateur. 15 La figure 5 montre un exemple de la densité spectrale d'amplitude 500 de la tension de réponse d'un oscillateur d'un gyromètre, obtenue grâce au premier mode de mise en oeuvre qui vient d'être décrit. Dans ce cas, la fréquence Fo de la tension sinusoïdale est égale à 32,5 kHz et la largeur spectrale AF du sinus cardinal est 5 égale à 5 kHz. La zone de recherche de la fréquence de résonance s'étend donc de 30 kHz à 35 kHz. La figure 5 met en évidence un pic de résonance 509 à la fréquence de résonance FR approximativement égale à 32.9 kHz. Une fois la fréquence de résonance connue, le facteur de qualité peut être déterminé par une mesure du décrément logarithmique qui consiste à mesurer le 10 temps d'amortissement après une excitation de l'oscillateur par une tension sinusoïdale à la fréquence de résonance. Le facteur de qualité peut également être déterminé en calculant le rapport de la fréquence de résonance sur la largeur du pic de résonance 509 mesurée à la moitié de sa hauteur. La tension d'excitation appliquée en entrée de l'oscillateur est transmise en 15 sortie de l'oscillateur par couplage capacitif. Un inconvénient de ce premier mode de mise en oeuvre est donc que le pic de Dirac 302 appartenant à la densité spectrale d'amplitude 300 de la tension d'excitation se retrouve dans la zone de recherche de la fréquence de résonance, ce pic 302 étant en effet présent à la fréquence Fo. Il s'agit d'un inconvénient mineur car ce pic de Dirac 302 présente une faible largeur 20 spectrale, inférieure à quelques dixièmes de hertz, alors que la zone de recherche s'étend sur une largeur de préférence supérieure à 1 kHz. Le phénomène de couplage capacitif est visible sur la figure 5 sous la forme d'un pic d'amplitude 508 présent à la fréquence Fo ici égale à 32,5 kHz. Toutefois, un deuxième mode de mise en oeuvre du procédé de l'invention 25 permet de remédier à cet inconvénient. Selon ce deuxième mode de mise en oeuvre, l'impulsion de tension utilisée pour former la tension d'excitation n'est plus une réponse impulsionnelle d'un filtre passe-bas, mais celle d'un filtre passe-bande. De préférence, cette réponse impulsionnelle a la forme d'une fonction sinus cardinal de largeur spectrale AF modulée par une fonction sinus de fréquence de 30 modulation AF'. 16 Selon ce deuxième mode de mise en oeuvre, on génère ainsi une première tension d'excitation et une deuxième tension d'excitation en additionnant l'impulsion de tension respectivement à une première tension sinusoïdale ayant première fréquence F1 et à une deuxième tension sinusoïdale ayant une deuxième 5 fréquence F2. De préférence la première tension sinusoïdale et la deuxième tension sinusoïdale ont la même amplitude de manière à obtenir des forces d'excitation ayant la même amplitude. La première tension d'excitation ui(t) et la deuxième tension d'excitation u2(t) peuvent être définies par les équations suivantes : Ill (t) = A. sinc(AF. t)sin(2. n. AF'. t) + B. sin(2. n. F1. t) (6.1) 10 u2(t) = A. sinc(AF. t)sin(2. n. AF'. t) + B. sin(2. n. F2. t) (6.2) où A est l'amplitude du sinus cardinal et B est l'amplitude des tensions sinusoïdales. La figure 6 représente schématiquement les densités spectrales d'amplitude de la première tension d'excitation ui(t) et de la première tension d'excitation u2(t), par exemple calculées au moyen d'une transformée de Fourier. Ces densités 15 spectrales comportent chacune : - deux premières composantes 605, 606 correspondant à la densité spectrale d'amplitude du sinus cardinal modulé à la fréquence AF', ces composantes 605, 606 étant des fonctions porte de largeur spectrale AF centrées respectivement sur -AF' et sur +AF' ; 20 - une deuxième composante 612, 622 correspondant respectivement à la densité spectrale d'amplitude de la première tension sinusoïdale et à la densité spectrale d'amplitude de la deuxième tension sinusoïdale, la deuxième composante 612, 622 étant un pic de Dirac situé à la fréquence F1 dans le cas de la première tension d'excitation et à la fréquence F2 dans le 25 cas de la deuxième tension d'excitation. La première tension d'excitation et la deuxième tension d'excitation sont successivement appliquées en entrée de l'oscillateur. Il en résulte respectivement une première force d'excitation Fi (t) et une deuxième force d'excitation F2(t) qui s'exercent successivement sur l'élément mobile de l'oscillateur et qui peuvent être 30 exprimées par les équations suivantes : 17 Fi (t) = k [A2. sinc2(AF. t). (1-cos(4.Tr.F1.0) + - -B2. cos(4.n. Fi. t) + 2 2 2 2. A. B. sinc(AF. t). sin(2. ii. AF'. t) . sin(2. n. Fi. t)1 (7.1) k [A2. sinc2(AF. t). (1-cos(4.Tr.F2.0 + B2 - B2 F2 (t) = . cos(4. T. F2. t) + 2 I 2 2 2. A. B. sinc(AF. t). sin(2. ii. AF'. t) . sin(2. n. F2. t)1 (7.2) 5 où k est un coefficient de proportionnalité. Ces équations qui expriment les forces d'excitation comportent chacune un dernier terme proportionnel au produit de l'impulsion de tension et de la tension sinusoïdale. La figure 7 représente schématiquement et partiellement la superposition des densités spectrales d'amplitude de la première force d'excitation Fi(t) et de la 10 deuxième force d'excitation F2(t). La figure 7 montre seulement les composantes spectrales à partir desquelles la zone de recherche de la fréquence de résonance est reconstituée. La densité spectrale d'amplitude de la première force d'excitation comporte deux composantes 711, 712 correspondant à la densité spectrale d'amplitude du dernier terme de l'équation (7.1), ces composantes 711, 712 étant 15 des fonctions porte de largeur spectrale AF centrées respectivement sur les fréquences F1 - AF' et sur F1 + AF'. Similairement, la densité spectrale d'amplitude de la deuxième force d'excitation comporte deux composantes 721, 722 correspondant à la densité spectrale d'amplitude du dernier terme de l'équation (7.2), ces composantes 721, 722 étant des fonctions porte de largeur spectrale AF centrées 20 respectivement sur les fréquences F2 - AF' et sur F2 + AF'. La zone de recherche ainsi obtenue est centrée sur une fréquence médiane Fo égale à la moitié de la somme de la première fréquence et de la deuxième fréquence (F0 - F1+2F2). Avantageusement, la fréquence de modulation AF' de la fonction sinus est égale à la largeur spectrale AF de la fonction sinus cardinal (AF' = 25 AF), et la deuxième fréquence F2 est égale à la somme de la première fréquence F1 et de la largeur spectrale AF (F2 = F1 + AF). La zone de recherche est ainsi continue et présente une largeur égale à quatre fois la largeur spectrale AF. 18 Si la fréquence de modulation AF' de la fonction sinus est inférieure à la largeur spectrale AF de la fonction sinus cardinal (AF' AF) et/ou si la deuxième fréquence F2 est supérieure à la somme de la première fréquence F1 et de 10 la largeur spectrale AF (F2 > F1 + AF), la densité spectrale reconstituée n'est plus continue et est donc nulle dans certaines bandes de fréquences, ce qui peut conduire à manquer la fréquence de résonance. Selon le deuxième mode de mise en oeuvre, le procédé comporte une première étape d'acquisition dans le domaine temporel d'une première tension de 15 réponse présente en sortie de l'oscillateur lorsque la première tension d'excitation est appliquée en entrée de l'oscillateur, et une deuxième étape d'acquisition dans le domaine temporel d'une deuxième tension de réponse présente en sortie de l'oscillateur lorsque la deuxième tension d'excitation est appliquée en entrée de l'oscillateur. Pour avoir une même résolution d'une acquisition à l'autre, les étapes zo d'acquisition sont de préférences effectuées avec une même fréquence d'échantillonnage Fe telle que : Fe > 2. (Fo + 2.AF ) La première tension de réponse et la deuxième tension de réponse sont transformées dans le domaine fréquentiel pour obtenir respectivement une première densité spectrale d'amplitude et une deuxième densité spectrale d'amplitude. Ces 25 densités spectrales d'amplitude sont par exemple calculées au moyen d'une transformée de Fourier. La première densité spectrale et la deuxième densité spectrale sont filtrées pour ne conserver que les portions de spectre correspondant aux composantes spectrales représentées aux figures 7A et 7B. Il s'agit par exemple d'un filtrage 19 numérique. Une densité spectrale d'amplitude est ensuite reconstituée en superposant les portions de spectre conservées. Les pics de Dirac présents aux fréquences F1 et F2, qui sont dus au couplage capacitif, sont ainsi éliminés de la densité spectrale d'amplitude reconstituée. La densité spectrale d'amplitude 5 reconstituée est équivalente à celle qu'on obtiendrait avec une force électrostatique non perturbée par le couplage capacitif. La figure 8 montre un exemple de la densité spectrale d'amplitude 800 de la tension de réponse du même oscillateur que pour la figure 5, obtenue grâce au deuxième mode de mise qui vient d'être décrit. Dans ce cas, la fréquence médiane 10 Fo de la zone de recherche est égale à 32,5 kHz. La largeur spectrale AF de la fonction sinus cardinal et la fréquence de modulation AF' de la fonction sinus sont égales à 1,25 kHz. Par conséquent la première fréquence F1 est égale à 31,875 kHz et la deuxième fréquence F2 est égale à 33,125 kHz. La zone de recherche de la fréquence de résonance s'étend donc de 30 kHz à 35 kHz. Un pic de résonance 809 15 est toujours présent à la fréquence de résonance FR approximativement égale à 32.9 kHz. Par contre, les pics de Dirac 612, 622 appartenant aux densités spectrales d'amplitude 610, 620 de la première tension d'excitation et de la deuxième tension d'excitation sont bien absents. La figure 9A montre trois exemples de densités spectrales d'amplitude 911, 20 912, 913 d'une force d'excitation formées selon le deuxième mode de mise en oeuvre. La fréquence médiane Fo de la zone de recherche définie par ces densités spectrales 911, 912, 913 est respectivement égale respectivement à 32,5 kHz, 33 kHz et 34 kHz. A chacun de ces spectres 911, 912, 913 correspond une densité spectrale d'amplitude reconstituée 921, 922, 923 de la tension de réponse de 25 l'oscillateur. Grâce à l'utilisation de la somme du sinus cardinal et de la tension sinusoïdale en tant que tension d'excitation, l'amplitude des densités spectrales 911, 912, 913 est constante sur toute la zone de recherche. Par conséquent, les densités spectrales reconstituées 921, 922, 923 comportent chacune un pic de résonance 30 ayant la même amplitude, comme mieux visible sur la figure 9B. L'amplitude du pic ne dépend donc pas de sa position par rapports aux bords de la zone de recherche, contrairement au cas de l'art antérieur illustré sur les figures 2A et 2B. 20 La figure 10 montre une densité spectrale d'amplitude 1001 de la tension de réponse d'un oscillateur obtenue par la méthode de balayage en fréquence de l'art antérieur. Dans ce cas, la mesure a duré environ 81 minutes pour une plage de fréquence de 75 Hz et un pas en fréquence de 0,02 Hz. La figure 10 montre 5 également, pour le même oscillateur, une densité spectrale d'amplitude 1002 de la tension de réponse obtenue par le procédé de l'invention. Dans ce cas, la mesure a duré environ 40 secondes pour une plage de fréquence de 500 Hz et un pas en fréquence de 0,25Hz (cette courbe 1002 est en fait le résultat de la moyenne de 10 courbes correspondant chacune à un temps d'acquisition de 4s à la fréquence 10 d'échantillonnage de 25 kHz ; en effectuant la moyenne de plusieurs acquisitions, le bruit est réduit). On constate que les pics de résonances des deux densités spectrales d'amplitude 1001, 1002 sont superposés par contre le procédé de l'invention permet de caractériser l'oscillateur beaucoup plus rapidement. Le procédé selon l'invention est par exemple mis en oeuvre par le biais d'une 15 application. Une application désigne un programme d'ordinateur comprenant des instructions exécutables par un processeur. Naturellement, l'invention n'est pas limitée aux modes de mise en oeuvre décrits en référence aux figures et des variantes pourraient être envisagées sans sortir du cadre de l'invention. REVENDICATIONS1. Procédé de détermination de paramètres caractéristiques d'un oscillateur à actionnement électrostatique, le procédé comportant les étapes suivantes : - générer une première tension d'excitation définie comme étant la somme d'une première tension sinusoïdale et d'une impulsion de tension ; - appliquer la première tension d'excitation en entrée de l'oscillateur ; - acquérir dans le domaine temporel une première tension de réponse présente en sortie de l'oscillateur lorsque la première tension d'excitation est appliquée en entrée de l'oscillateur ; - obtenir, par transformation dans le domaine fréquentiel, une première densité spectrale d'amplitude (500) de la première tension de réponse ; - déterminer les paramètres caractéristiques de l'oscillateur à partir de la première densité spectrale d'amplitude (500). 2. Procédé selon la revendication précédente, dans lequel les paramètres caractéristiques comprennent une fréquence de résonance, le procédé comportant les étapes initiales suivantes : - estimer la fréquence de résonance de l'oscillateur ; - définir une plage de fréquence, dite « zone de recherche », dans laquelle la fréquence de résonance est recherchée, la zone de recherche étant centrée autour d'une fréquence médiane égale à la fréquence de résonance estimée, la zone de recherche présentant une largeur prédéterminée. 3. Procédé selon l'une quelconque des revendications 1 et 2, dans lequel l'impulsion de tension est une réponse impulsionnelle d'un filtre passe-bas. 4. Procédé selon la revendication 3 lorsqu'elle dépend de la revendication 2, dans lequel l'impulsion de tension présente une largeur spectrale (OF) égale la largeur de la zone de recherche, la première tension sinusoïdale présentant une première fréquence (Fo) égale à la fréquence médiane. 5. Procédé selon l'une quelconque des revendications 3 et 4, dans lequel l'impulsion de tension se présente sous la forme d'une fonction sinus cardinal. 22 6. Procédé selon l'une quelconque des revendications 1 et 2, dans lequel l'impulsion de tension est une réponse impulsionnelle finie d'un filtre passe-bande, le procédé comportant également les étapes suivantes : - générer une deuxième tension d'excitation définie comme étant la somme 5 d'une deuxième tension sinusoïdale et de l'impulsion de tension ; - appliquer la deuxième tension d'excitation en entrée de l'oscillateur ; - acquérir dans le domaine temporel une deuxième tension de réponse présente en sortie de l'oscillateur lorsque la deuxième tension d'excitation est appliquée en entrée de l'oscillateur ; 10 - obtenir, par transformation dans le domaine fréquentiel, une deuxième densité spectrale d'amplitude de la deuxième tension de réponse ; - reconstituer une densité spectrale d'amplitude (800) à partir de la première densité spectrale d'amplitude et de la deuxième densité spectrale d'amplitude ; 15 les paramètres caractéristiques de l'oscillateur étant déterminés à partir de la densité spectrale d'amplitude reconstituée (800). 7. Procédé selon la revendication précédente, dans lequel l'impulsion de tension se présente sous la forme d'une fonction sinus cardinal modulée par une fonction sinus. zo 8. Procédé selon la revendication précédente, dans lequel la fonction sinus cardinal présente une largeur spectrale (AF), la fonction sinus présentant une fréquence de modulation (AF') égale à la largeur spectrale (AF) de la fonction sinus cardinal, la première tension sinusoïdale présentant une première fréquence (F1) et la deuxième tension sinusoïdale présentant une deuxième fréquence (F2) égale à la 25 somme de la première fréquence (F1) et de la largeur spectrale (AF) de la fonction sinus cardinal, la densité spectrale d'amplitude reconstituée (800) étant continue. 9. Procédé selon la revendication précédente lorsqu'elle dépend de la revendication 2, dans lequel la largeur spectrale (AF) de la fonction sinus cardinal est égale à un quart de la largeur de la zone de recherche, la fréquence médiane étant 30 égale à la moitié de la somme de la première fréquence (F1) et de la deuxième fréquence (F2). 23 10. Produit programme d'ordinateur comprenant des instructions pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 9, lorsque ces instructions sont exécutées par un processeur.