FR 2835071 A1 20030725 FR 0200789 A 20020122 de la dimension de variation de l'opacité :r 1 2835071 Procédé de ffltrage d'un signal numérique, synthétiseur de fréquence et programme d'ordinateur correspondants. Le domaine de l' invention est celui du filtrage de signaux numériques. Le procédé de filtrage selon l' invention a de nombreuses applications. I1 peut notamment, mais non exclusivement, être mis en _uvre dans les synthétiseurs de fréquence, utilisés dans des appareils électroniques comme les récepteurs (de radiocommunication ou autres), les générateurs de signaux, les modulateurs, les démodulateurs, les mélangeurs, etc. Plus généralement, il peut s'appliquer dans tous les dispositifs ou systèmes qui utilisent des signaux sinusoïdaux devant posséder une bonne pureté spectrale. On suppose que le signal d'origine à filtrer est un signal numérique sinusoidal, comprenant une composante sinusoïdale dominante à laquelle se superposent des signaux parasites. Ce signal numérique d'origine à filtrer est par exemple obtenu avec un système de synthèse de fréquence numérique directe (système SND). Ce système SND peut éventuellement être combiné à un convertisseur numérique / analogique (CNA), si on désire produire un signal analogique. I1 est connu que, lors de la synthèse numérique directe (SND), des signaux (raies) parasites sont générés en raison des différentes troncatures effectuées sur le nombre de bits utilisés (accumulateur de phase, adressage de la table du sinus et largeur du bus des données), ainsi que par les valeurs non entières du rapport fréquence d'horloge sur fréquence du signal synthétisé. Les solutions connues pour réduire le niveau des signaux parasites consistent soit à limiter l'utilisation du synthétiseur à des fréquences o les parasites ont des niveaux faibles, soit à augmenter le nombre des bits des différents éléments utilisés dans la SND. La première solution précitée est inacceptable dans les applications o la fréquence doit être variable. La seconde solution précitée est coûteuse et limite la vitesse de traitement. Selon encore une autre technique connue de filtrage, on ajoute une séquence pseudo-aléatoire à la phase du signal (méthode dite "dither method"). Cette autre technique connue est présentée par exemple dans le document ("data sheet") de la société GRAYCHIP INC intitulé: "Multistandard quad DDC chip" rev 0.3, April 20, 2 2835071 2000. Il y est montré que, dans certains cas, on peut réduire le niveau des raies parasites d'environ 14 dB. Cependant, le potentiel d'amélioration de cette autre technique connue est intrinsèquement limité car elle consiste à répartir l'énergie des parasites dans toute la bande de fréquence. L' invention a notamment pour objectif de pallier ces différents inconvénients de l'état de la technique. Plus précisément, l'un des objectifs de la présente invention est de fournir un procédé de filtrage d'un signal numérique qui soit simple à mettre en _uvre et peu coûteux en nombre d'opérations, de façon à ne pas nécessiter des temps de calcul élevés et pouvoir être mis en _uvre en temps réel. Un autre objectif de l' invention est de fournir un tel procédé permettant d'obtenir une amélioration supérieure à 14 dB, et pouvant aller jusqu'à 60 dB, voire au-delà. Ces différents objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention à l'aide d'un procédé de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons si, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoïdale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, le procédé comprenant une étape de filtrage numérique basée sur une simplification d'une matrice apparaissant dans la décomposition en valeurs singulières (SVD) de la matrice de Hankel du signal numérique d'origine. On notera qu'il n'était pas évident d'utiliser la décomposition SVD de la matrice de Hankel d'un signal, dans le présent contexte de filtrage d'un signal numérique sinusoidal. En effet, c'est une méthode de filtrage connue dans le domaine du traitement des signaux, pour extraire plu sieurs compos antes sinus odales d'un bruit. L'inconvénient majeur de cette méthode est qu'elle est gourmande du point de vue du nombre des opérations à effectuer. Il était donc difficile d'envisager son utilisation lorsqu'une mise en _uvre en temps rcel est souhaitée. 3 2835071 La présente invention repose sur une simplification de la décomposition en valeurs singulières de la matrice de Hankel du signal, et plus précisément sur une simplification d'une matrice apparaissant dans cette décomposition. Cette simplification s'appuie sur le fait que, dans le présent contexte, le signal d'origine ne comporte, outre le bruit, qu'une composante sinusoïdale dominante. Comme expliqué en détails par la suite, la matrice l; obtenue après simplification, Émodifée, permet d'obtenir une matrice de Hankel filtrée, Hfi,rée, dont les première et dernière lignes (qui mises bout à bout constituent la séquence de N échantillons "sfiltré;" du signal filtré) peuvent être calculées de façon simple et rapide. De façon avantageuse, la matrice I: est une matrice n x p s'écrivant: , O. O O 0 â2 0. O É- O 0 â3 À _.. O O.. O o les éléments 6j sont des valeurs singulières classées dans l'ordre décroissant et correspondent aux puissances des composantes présentes dans le signal numérique d'origine. La simplification de la matrice consiste à annuler toutes les valeurs de c,; sauf la plus grande a" correspondant à la composante sinusodale dominante du signal numérique d'origine. Ainsi, la matrice obtenue après simplification, Émodifiée, est très simple avec un seul élément non nul (a,). Avantageusement, la simplification de la matrice consiste en outre à poser a, = 1. I1 s'agit ici d'une simple question de normalisation. Dans un mode de réalisation préférentiel de l' invention, le procédé comprend les étapes suivantes: - obtention d'un vecteur v, dont les composantes Vk, k variant de 1 à N/2, sont les N/2 valeurs de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine; - calcul de scalaires u, et u: tels que: 4 2835071 N/2 * U1 = Évksk k=1 N12 * U2 = ÉVkSk+NI2 k =1 - calcul d'une séquence d'échantillons sfiltréj, i variant de 1 à N. représentant le signal numérique filtré, selon les équations suivantes: * sfiltrék = u vk k = 1, 2, 3 N/2 * sfiltrék+N,2 = u2 vk k = 1, 2, 3 N/2 De cette façon, le filtrage numérique appliqué au signal numérique d'origine est rapide et peu coûteux en puissance de calcul. Les équations de calcul des scalaires u, et u2, ainsi que les équations de calcul des N échantillons du signal filtré, résultent d'une simplification, telle que discutée ci dessus et présentée en détails par la suite, d'une matrice apparaissant dans la décomposition en valeurs singulières (SVD) de la matrice de Hankel du signal numérique d'origine. Dans l'ensemble de la présente description (et comme cela apparaît au vu de l'équation (e) (voir ci-après) des composantes vk du vecteur v), on entend par " fonction d'auto-corrélation translatée " la fonction d'autocorrélation translatée d'une unité et telle que le décalage nul correspond à k=1 (au lieu de k=0 dans le cas classique). De façon avantageuse, l'étape d'obtention du vecteur v consiste à calculer les composantes Vk, k variant de 1 à N/2, selon l'équation suivante: (N12)+1 Vk = SiS,;k-1 k = 1, 2, 3... N/2 i=l o si* est le complexe conjugué de s;. Ce calcul peut être appliqué dans tous les cas, que l'on connaisse ou non la fréquence du signal numérique d'origine, puisqu'il est effectué à partir des échantillons du signal numérique d'origine. Selon une variante avantageuse, dans laquelle la fréquence du signal numérique d'origine est connue, l'étape d'obtention du vecteur v consiste à calculer N/2 valeurs théoriques de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine. Ainsi, cette variante s'applique lorsqu'on connaît la fréquence du signal numérique d'origine. Ceci est notamment le cas lorsque le signal numérique d'origine résulte d'une synthèse de fréquence numérique directe (SND). Cette variante permet d'obtenir de meilleurs résultats en termes de qualité spectrale. De façon préférentielle, l'étape d'obtention du vecteur v est effectuée une fois pour toutes, pour une fréquence donnée du signal numérique d'origine. En d'autres termes, le vecteur v ne doit être recalculé qu'en cas de changement de fréquence du signal numérique d'origine. Dans un mode de réalisation particulier de l'invention, l'étape d'obtention du vecteur v est effectuée de façon préalable plusieurs fois, pour différentes fréquences possibles du signal numérique d'origine, de façon à disposer d'un jeu de vecteurs v associés chacun à une fréquence distincte. Le procédé comprend une étape de stockage du jeu de vecteurs v, et, lors de l'étape de filtrage numérique du signal numérique d'origine possédant une fréquence donnée, l'obtention du vecteur v consiste à lire dans les moyens de stockage le vecteur v associé à ladite fréquence donnée. Cette variante est particulièrement intéressante si les changements de fréquence sont fréquents (par exemple synthèse d'un signal "chirp"). En effet, elle permet de gagner du temps lors du filtrage du signal numérique d'origine, les vecteurs v ayant été calculés précédemment. L'invention concerne également un synthétiseur de fréquence comprenant des moyens permettant la mise en _uvre du procédé précité. L'invention concerne aussi un programme d'ordinateur, comprenant des séquences d' instructions adaptées à la mise en _uvre d'un procédé tel que précité lorsque ledit programme est exécuté sur un ordinateur. L'invention concerne encore un produit programme d'ordinateur de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoïdale dominante superposée à des signaux parasites, ledit procédé 6 2835071 permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, ledit produit programme d'ordinateur comprenant des instructions de code de programme enregistré sur un support utilisable dans un ordinateur comprenant des moyens de programmation lisibles par ordinateur pour effectuer une étape de filtrage numérique basée sur une simplification d'une matrice apparaissant dans la décomposition en valeurs singulières (SVD) de la matrice de Hankel du signal numérique d'origine. L'invention concerne aussi un produit programme d'ordinateur de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons si, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoidale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, ledit produit programme d'ordinateur comprenant des instructions de code de programme enregistré sur un support utilisable dans un ordinateur comprenant: - des moyens de programmation lisibles par ordinateur pour effectuer une étape d'obtention d'un vecteur v, dont les composantes Vk, k variant de 1 à N/2, sont les N/2 valeurs de la fonction d'autocorrélation translatée de la séquence d'échantillons si du signal numérique d'origine; - des moyens de programmation lisibles par ordinateur pour effectuer une étape de calcul de scalaires u, et U2 tels que: N/2 * U =ÉVkSk k=] N/2 * U2 = iVkSk+N/2 k=l - des moyens de programmation lisibles par ordinateur pour effectuer une étape de calcul d'une séquence d'échantillons sfiltré;, i variant de 1 à N. représentant le signal numérique filtré, selon les équations suivantes: * sfiltrék = u, vk k = 1, 2, 3 N/2 * sfiltrék+N,2 = u2 vk k = 1, 2, 3... N/2 (produit programme d'ordinateur). D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante d'un mode de réalisation préférentiel de l' invention, donné à titre d'exemple indicatif et non limitatif, et des dessins annexés, dans lesquels: - la figure 1 présente un organigramme simplifié d'un mode de réalisation particulier du procédé de filtrage de signal numérique selon l'invention; - la figure 2 est un synoptique illustrant la mise en _uvre d'un mode de réalisation particulier du procédé de filtrage selon l'invention; - la figure 3 présente le spectre d'un exemple de signal numérique d'origine, généré par synthèse numérique directe; - les figures 4 à 6 présentent chacune le spectre du signal numérique d'origine de la figure 3, après filtrage selon un mode de réalisation particulier de la présente invention, par blocs de 512 points (fig.4), 2048 points (fig.5) et 4096 points (fig.6) respectivement. L'invention concerne donc un procédé de filtrage d'un signal numérique d'origine. Comme illustré sur l'organigramme simplifié de la figure 1, on suppose dans la suite de la description que le signal numérique d'origine résulte d'une étape (1) de synthèse permettant d'obtenir un signal numérique d'origine (deux voies I et Q quadrature), représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter. Cette étape (1) de synthèse consiste par exemple en la mise en _uvre d'une synthèse de fréquence numérique directe (SND). La technique SND est bien connue de l'homme du métier et ne sera pas décrite plus en détails ci-après. On notera qu'avec cette technique, on connaît la fréquence du signal synthétisé. Il est clair cependant que la présente invention s'applique que l'on connaisse ou non la fréquence du signal numérique d'origine. Optionnellement, si on souhaite obtenir un signal analogique possédant une grande pureté spectrale, le filtrage selon le procédé de l'invention (2) est suivi d'une étape (3) de conversion numérique/analogique. 8 2835071 Dans le mode de réalisation particulier illustré sur la figure 1, le filtrage numérique (2) selon le procédé de l'invention comprend: - une étape (21) d'obtention d'un vecteur v, dont les composantes Vk, k variant de 1 à N/2, sont les N/2 valeurs de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine; - une étape (22) de calcul de scalaires u, et U2 tels que: N/2 * lll = ÉVkSk (a) k=l N/2 2 VkSk+N/2 (b) k=l - une étape (23) de calcul d'une séquence d'échantillons sfiltréi, i variant de 1 à N. représentant le signal numérique filtré, selon les équations suivantes: * sfiltrék = u, vk k = 1, 2, 3... N/2 (c) * sfiltrck+N,2 = u2 vk k = 1, 2, 3... N/2 (d) Les étapes (22) de calcul de scalaires u, et u2 et (23) de calcul d'une séquence d'échantillons sfiltré; nécessitent 2N multiplications complexes et N-2 additions complexes. Ainsi, pour traiter N = 512 points par exemple, il faut environ 1024 multiplications et 510 additions. Le temps de traitement d'une séquence de N points dépend linéairement de N et le procédé selon l'invention est donc particulièrement intéressant du point de vue du temps de calcul. L'étape (21) d'obtention du vecteur v peut être effectuce de deux façons différentes,. Dans le cas général, applicable notamment si on ne connaît pas la fréquence du signal synthétisé (et a fortiori si on la connaît), l'étape (21) d'obtention du vecteur v consiste à calculer les composantes Vk, k variant de 1 à N/2, selon l'équation suivante: (N/2)+1 Vk = SiS,+k-1 k = 1, 2, 3... N/2 (e) i=] o sj* est le complexe conjugué de s;. Dans ce cas général, le calcul du vecteur v nécessite (N2/4+N/2) multiplications complexes et N2/4 additions complexes. Ainsi, pour N = 512 par exemple, il faut environ 66 103 multiplications et 66 103 additions. Comme cela est indiqué ci-après, ce calcul peut être effectué une fois pour toutes. Une variante peut être envisagée dans le cas particulier o on connat la fréquence du signal synthétisé. C'est par exemple le cas si l'étape (1) de synthèse consiste en la mise en _uvre d'une synthèse de fréquence numérique directe (SND). Selon cette variante, l'étape (21) d'obtention du vecteur v consiste à calculer N/2 valeurs théoriques de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine. On obtient alors de meilleurs résultats en termes de qualité spectrale. On rappelle que la fonction d'auto-corrélation translatée du signal synthétisé, qui est un signal sinusoïdal, est donnée par: vk = exp [j2F(k-l)] o F est la fréquence du signal synthétisé et k variant de 1 à N/2, les N/2 instants d'échantillonnage. On présente maintenant, en relation avec le synoptique de la figure 2, la mise en _uvre d'un mode de réalisation particulier du procédé de filtrage selon l'invention. On observe que les équations (a) et (b) se calculent de la même façon, et qu'il en est de même pour les équations (c) et (d). Du point de vue réalisation, il sufft donc de concevoir un circuit calculant le signal filtré pour k=1 à N/2, et le même circuit servira pour k=(N/2)+1 à N. Dans le montage présenté en figure 2, on suppose qu'un compteur 8 est cadencé par une horloge H et possède une capacité N/2. Les sorties du compteur 8 permettent d'adresser une mémoire RAM 9 qui contient les N/2 composantes du vecteur v, qui sont supposées avoir été calculées au préalable (voir discussion détaillée ci-après). Le nombre u, est calculé par l'association d'un premier multiplicateur 4 et d'un accumulateur 5. Plus précisément, à l'instant i, la valeur v; est présentée en entrée du premier multiplicateur 4 et elle est multipliée par si. L'opération est effectuée pour i variant de 1 à N/2. Pour i=N/2, l'accumulateur 5 contient la valeur u (équation (a)). Cette valeur est stockée dans une mémoire 6 par le signal de sortie du décodeur 7, lorsque le compteur 8 a atteint sa valeur maximale. Le signal de sortie du décodeur 7, retardé par un élément référencé 11, provoque ensuite la remise à zéro de l'accumulateur 2835071 5. La valeur de u, est ensuite utilisée pour calculer le signal filtré "sfiltrck", pour k=1 à N/2 (équation (c)), dans le second multiplicateur 10. Durant cette phase, la valeur u2 est calculée par le premier multiplicateur 4 (équation (b). Lorsque la valeur de u2 est disponible, elle est utilisce dans le second multiplicateur 10 pour générer les N/2 éléments suivants du signal filtré, selon l'équation (d). Comme indiqué ci-dessus, les N/2 composantes du vecteur v peuvent être calculées au préalable, une fois pour toutes, pour une fréquence donnée du signal numérique d'origine. Dans ce cas, le vecteur v est stocké et il ne doit étre recalculé qu'en cas de changement de fréquence du signal synthétisé. Si les changements de fréquence sont fréquents (par exemple synthèse d'un signal "chirp"), on peut envisager une variante, pour gagner du temps, consistant à calculer au préalable plusieurs vecteurs v, associés chacun à une fréquence particulière possible du signal synthétisé, et à stocker ce jeu de vecteurs v dans une mémoire (RAM, ROM, EPROM,...). Ainsi, lors du filtrage numérique (2), l'obtention du vecteur v consiste à lire dans cette mémoire le vecteur v, parmi la pluralité de vecteurs stockés, associé à la fréquence du signal numérique d'origine que l'on souhaite filtrer. Les performances du procédé selon l'invention, dans le mode de réalisation présenté ci-dessus en relation avec les figures 1 et 2, ont été évaluces à l'aide d'un programme simulant un générateur de signal SND. La figure 3 illustre le spectre P1 l d'un exemple de signal généré à la fréquence F=11,33 MHz, pour une fréquence d'horloge FH=64 MHz. La séquence analysée comporte 16 348 points et elle est apodisée par une fenêtre de Kaiser. La raie parasite la plus forte se trouve à -60 dBc et le niveau moyen des parasites se situe entre -70 et -80 dBc. Les figures 4 à 6 présentent chacune le spectre du signal numérique d'origine de la figure 3, après filtrage selon l'invention, par blocs de 512 points (fg.4) ' 2048 points (fig. 5) et 4096 points (fig. 6) respectivement. On observe que le niveau des raies parasites diminue progressivement lorsque la fréquence s'éloigne de celle de la porteuse. On observe également que le signal devient de plus en plus pur lorsque le nombre de points traités augmente. Le degré d' améliorati on du spectre dépend légèrement de la fréquence choi sie, mais les simulations montrent que le traitement sur 2 048 points permet d'obtenir une pureté spectrale voisine de 100-110 dB quelle que soit la fréquence synthétisée. On présente ci-après le cheminement intellectuel de l'inventeur pour aboutir aux équations précitées. On suppose que le signal comprend N valeurs s={s, S2 S3 S4... SN}. La matrice de Hankel de ce signal est une matrice n x p, avec p=N/2 et n=(N/2)+1=p+1, formée de la façon suivante: s, s2 À À sp S2 S3. À H = S3 Sn À À À SN Cette écriture met en évidence le fait que l'on retrouve la séquence initiale en prenant la lèrC ligne de H. suivie de la dernière ligne. La décomposition SVD (décomposition en valeurs singulières) de cette matrice de Hankel du signal s'écrit: H = VI;VH o U et V sont des matrices n x n et p x p respectivement, formées de vecteurs singuliers, VH indique la transposée hermitienne de V et est une matrice n x p composée des valeurs singulières classées dans l'ordre décroissant: 12 2835071 cy, O. O O O (T2 0. O O 0 3 = À À P O O.. O Ul l U12 À À À Uln Vl I V12.. Vlp U21 U22 À À À À V21 V22.. . U=. V=. . U',, Un2 À À À Unn _ V plV p2. V pp Les valeurs 6j correspondent aux puissances des composantes présentes dans le Ssignal. Les valeurs les plus faibles correspondent au bruit. Le principe général de ce type de filtrage consiste à mettre les valeurs qui correspondent au bruit à zéro, afin de l'éliminer, puis d'effectuer le produit des matrices pour obtenir une matrice de Hankel filtrée: Hfi" = U dii V Cette méthode de filtrage, connue en elle-même, donne de bons résultats. Son inconvénient est qu'elle est lourde du point de vue temps de calcul. Dans le cadre de la présente invention, on simplifie la matrice en annulant toutes les valeurs de 6; sauf la plus grande 6" correspondant à la composante 1Ssinusoïdale dominante du signal numérique d'origine, et on pose 6, = 1. Dans ces conditions la matrice Émodifée devient simple avec un seul élément non nul égal à 1, les calculs se simplifient et on obtient: 13 2835071 * * u,lv" u"v2,. . u"vp, * * * U2IVIl u2,v2, À â21VPI Hfltrée =.. .. . .. UnlVII UnlV2. . UnlVp Pour retrouver le signal filtré, il suffit de calculer la première et la dernière ligne de Hfi,rée. On constate sur la matrice précédente que seules les deux composantes u" et un de la matrice U doivent étre connues et que seule la première ligne v"*, v2,*,...vp,* de la matrice VH doit étre connue. Les composantes v,, v2"...vp peuvent étre calculées à l'aide du produit matriciel suivant: v = HHH(:,1) o HH est la transposée hermitienne de H et H(:,1) représente la première colonne de H. Cette équation correspond à l'équation (e) précitée, utilisée lors de l'étape (21) d'obtention du vecteur v. Pour des raisons de vitesse de calcul, il a fallu trouver une astuce pour calculer u, et un rapidement, sans avoir à calculer la matrice U entière. Après étude, l'inventeur a trouvé que u, et Un pouvaient étre obtenus par: N12 Uri = ÉVkSk k=1 N12 k = 1,2,3 N /2 Unl = ÉVkSk+N12 k=l Ainsi, u" et un correspondent respectivement aux scalaires u, et u2 précités. Ces deux équations correspondent aux équations (a) et (b) utilisées lors de l'étape référencée 22 sur la figure 1. 14 2835071 On dispose maintenant de tous les elements pour passer à l'étape (23) de filtrage proprement dit, qui consiste à calculer la premiere et la dernière ligne de Hams pour obtenir le signal filtré par l'équation: Sf/tré {V UniV} Cette équation correspond aux équations (c) et (d) utilisées lors de l'étape S référencée 23 sur la figure 1. REVENDICATIONS 1. Procédé de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoïdale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, caractérisé en ce que le procédé comprend une étape de filtrage numérique basée sur une simplification d'une matrice l; apparaissant dans la décomposition en valeurs singulières (SVD) de la matrice de Hankel du signal numérique d'origine. 2. Procédé selon la revendication 1, caractérisé en ce que la matrice est une matrice n x p s'écrivant: a] O. O O O (T2 0. O É- O 0 â3 À _.. O O.. O o les éléments c'; sont des valeurs singulières classées dans l'ordre décroissant et correspondent aux puissances des composantes présentes dans le signal numérique d'origine, et en ce que la simplification de la matrice consiste à annuler toutes les valeurs de s; sauf la plus grande a" correspondant à la composante sinusoïdale dominante du signal numérique d'origine. 3. Procédé selon la revendication 2, caractérisé en ce que la simplification de la matrice l; consiste en outre à poser a, = 1. 4. Procédé de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoïdale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, caractérisé en ce que le procédé comprend les étapes suivantes: 16 2835071 obtention d'un vecteur v, dont les composantes Vk, k variant de 1 à N/2, sont les N/2 valeurs de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine; calcul de scalaires u, et U2 tels que: N/2 * U1 =ÉVkSk k=1 N/2 * U2 = 2.iVkSk+N,2 k=l - calcul d'une séquence d'échantillons sfiltré;, i variant de 1 à N. représentant le signal numérique filtré, selon les équations suivantes: * sfltrék = u1 vk k = 1, 2, 3 N/2 sfiltrék+N,2 = u2 vk k = 1, 2, 3 N/2 5. Procédé selon la revendication 4, caractérisé en ce que l'étape d'obtention du vecteur v consiste à calculer les composantes Vk, k variant de 1 à N/2, selon l'équation suivante: (N/2)+1 Vk SiS.;k-1 k = 1, 2, 3... N/2 i=1 1S o sj* est le complexe conjugué de s;. 6. Procédé selon la revendi cation 4, la fréquence du sig n al numéri que d' origine étant connue, caractérisé en ce que l'étape d'obtention du vecteur v consiste à calculer N/2 valeurs théoriques de la fonction d'auto-corrélation translatée de la séquence d'échantillons s; du signal numérique d'origine. 7. Procédé selon l'une quelconque des revendications 4 à 6, caractérisé en ce que l'étape d'obtention du vecteur v est effectuce une fois pour toutes, pour une fréquence donnée du signal numérique d'origine. 8. Procédé selon l'une quelconque des revendications 4 à 7, caractérisé en ce que l'étape d'obtention du vecteur v est effectuée de façon préalable plusieurs fois, pour différentes fréquences possibles du signal numérique d'origine, de façon à disposer d'un jeu de vecteurs v associés chacun à une fréquence distincte, en ce que le procédé comprend une étape de stockage du jeu de vecteurs v, 17 2835071 et en ce que, lors de l'étape de filtrage numérique du signal numérique d'origine possédant une fréquence donnée, l'obtention du vecteur v consiste à lire dans les moyens de stockage le vecteur v associé à ladite fréquence donnée. 9. Synthétiseur de fréquence, caractérisé en ce qu'il comprend des moyens permettant la mise en _uvre du procédé selon l'une quelconque des revendications 1 à 8. 10. Programme d'ordinateur, caractérisé en ce que ledit programme comprend des séquences d' instructions adaptées à la mise en _uvre d' un procédé selon l' une quelconque des revendications 1 à 8 lorsque ledit programme est exécuté sur un ordinateur. 11. Produit programme d'ordinateur de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusoïdale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, ledit produit programme d'ordinateur comprenant des instructions de code de programme enregistré sur un support utilisable dans un ordinateur comprenant des moyens de programmation lisibles par ordinateur pour effectuer une étape de filtrage numérique basée sur une simplification d'une matrice l; apparaissant dans la décomposition en valeurs singulières (SVD) de la matrice de Hankel du signal numérique d'origine. 12. Produit programme d'ordinateur de filtrage d'un signal numérique d'origine, représenté par une séquence d'échantillons s;, i variant de 1 à N. o N est le nombre d'échantillons à traiter, le signal numérique d'origine comprenant une composante sinusodale dominante superposée à des signaux parasites, ledit procédé permettant de diminuer le niveau des signaux parasites et obtenir un signal numérique filtré présentant une bonne qualité spectrale, ledit produit programme d'ordinateur comprenant des instructions de code de programme enregistré sur un support utilisable dans un ordinateur comprenant: 18 2835071 des moyens de programmation lisibles par ordinateur pour effectuer une étape d'obtention d'un vecteur v, dont les composantes Vk, k variant de 1 à N/2, sont les N/2 valeurs de la fonction d'auto-corrélation translatée de la séquence d'échantillons si du signal numérique d'origine; - des moyens de programmation lisibles par ordinateur pour effectuer une étape de calcul de scalaires u, et U2 tels que: N/2 * ui = ÉVkSk k=l N/2 2 VkSk+ N/2 k=l - des moyens de programmation lisibles par ordinateur pour effectuer une étape de calcul d'une séquence d'échantillons sfiltré;, i variant de 1 à N. représentant le signal numérique filtré, selon les équations suivantes: * sfiltrék = u, vk k = 1, 2, 3... N/2