FR 3047578 A1 20170811 FR 1650945 A 20160205 1 DOMAINE TECHNIQUE GENERAL L'invention concerne le domaine des télécommunications et notamment le domaine des systèmes et procédés permettant d'estimer certains paramètres constitutifs d'un ou de plusieurs signaux ou porteuses modulées contenus dans une même bande de fréquences. ETAT DE LA TECHNIQUE La multiplication des systèmes de télécommunications engendre un besoin de réutilisation des canaux afin de satisfaire la demande en termes de bande passante et de qualité de service. Une telle réutilisation est possible grâce aux réseaux cellulaires ou spots. La densification de ces réseaux associée à une réutilisation de certains canaux entre cellules proches ou adjacentes engendre des situations où un terminal reçoit un signal d'intérêt entaché d'un ou plusieurs signaux interférents dont la puissance se rapproche de celle du signal utile. Cette concomitance de signaux utilisant la même ressource fréquentielle peut selon les cas être involontaire (interférence co-canal) ou délibérée (la connaissance a priori de l'utilisation de la même fréquence par un autre émetteur est utilisée pour réceptionner correctement les télécommunications et optimiser l'utilisation du spectre). Le taux d'erreur des terminaux d'architecture classique est alors très dégradé et il faut employer des algorithmes de séparation ou de réduction d'interférences en vue d'isoler le signal utile. La connaissance de certaines caractéristiques des signaux présents dans la bande considérée est nécessaire au bon fonctionnement de ces algorithmes. Ces caractéristiques sont, par exemple, la modulation utilisée, la puissance, la fréquence porteuse ou encore le rythme symbole de chacun des signaux présents dans la bande considérée. On connait des méthodes d'estimation des caractéristiques du signal principal lorsque les signaux interférents sont de très faibles niveaux relatifs : Pour l'estimation des modulations, on pourra se référer au papier de journal [Swami & Sadler, 2000] ou également au papier de conférence [Reichert, 1992]. Ce dernier s'appuie sur la transformée MPT (en anglais, « Mth Power nonlinear Transform ») précédemment proposée dans le document [Mazo, 1978] en vue d'estimer la fréquence porteuse et le rythme symbole d'un signal d'intérêt en l'abscence de signaux interférents. Pour l'estimation d'autres paramètres (puissances, fréquences porteuses, rythmes symboles...), on pourra se référer à l'ouvrage de référence [Proakis, 2007]. 3047578 2 D'autres techniques ont été proposées pour l'estimation de certains paramètres lorsque plusieurs signaux sont présents sur la même bande de fréquence. Celles-ci utilisent des transformées dans le domaine spectral ou des statistiques d'ordres supérieurs. A ce titre, on pourra se référer : 5 - au document [Zaerin & Seyfe, 2012] qui présente une méthode d'estimation des modulations basée sur des « cumulants » d'ordres supérieurs. Cette technique est efficace pour estimer les constellations des signaux utiles et des interférents. Toutefois, elle se révèle inefficace lorsque les SINR (en anglais, « Signal to Interference plus Noise Ratio », qui définissent le rapport entre la puissance du 10 signal d'intérêt et la puissance combinée du signal interférent et du bruit additif) sont faibles, soit quand la puissance de l'interférent est comparable à celle du signal d'intérêt ; - au document [Guenther, 2010] dans lequel il est utilisé une fonction spectrale de corrélation (en anglais « Spectral Correlation Function », (SCF)) pour l'estimation 15 des caractéristiques. Cependant, le caractère coopératif de la méthode présentée (on connaît certains paramètres a priori) en rend son utilisation impossible dès lors que les signaux interférents présentent un fort niveau et qu'aucune information a priori ne peut être utilisée. De plus, il n'exploite pas les valeurs de l'amplitude des pics de la SCF. 20 Ainsi, les techniques connues présentent de nombreux inconvénients. PRESENTATION DE L'INVENTION L'invention propose de pallier les inconvénients précités. A cet effet, l'invention propose un procédé d'estimation de paramètres constitutifs 25 d'un signal seul ou de signaux mélangés en fréquence constitué des étapes suivantes : - réception d'un signal seul ou d'un mélange de signaux y(t) dans une bande de fréquences donnée ; - détermination pour L, nombre réel strictement positif pris dans un ensemble de q, q strictement positif, valeurs prédéterminées {L1; ;Lq} de la fonction de corrélation 30 spectrale définie ci-dessous : LSCFy(f, RyL (T, (f) RilL(T, Y où est la fonction d'autocorrélation -,j2]ra cyclique de l'élévation à la puissance L de la version échantillonné y(n) du signal 3047578 3 y(t) et où FI.] est la transformée de Fourier et où, pour chaque q, a peut indépendamment prendre une ou plusieurs valeurs discrètes et/ou parcourir une ou plusieurs plages de valeurs ; traitement de chaque tranformation LSCFy(f,a) L étant pris dans un ensemble de q, 5 q strictement positif, valeurs prédéterminées {L1; ;Lq} pour en extraire un ensemble de pics p,, c'est-à-dire de maxima locaux détectés selon les variables f et OC ou selon la seule variable f pour chaque valeur discrète de a, chacun de ces pics étant caractérisé par un triplet {LAp,, LFp,, Lapa} où LFp, et Lap, sont respectivement la fréquence et la fréquence cyclique à laquelle il se situe et LAp, son amplitude ; 10 identification de groupes de triplets {LAp,, LFp,, Lap,} comme étant la signature laissée par un ensemble de paramètres constitutifs de signaux présents dans la bande de fréquence donnée. Dans le contexte de l'invention, pour M réel, il est utile d'utiliser le cas particulier de la fonction MPT définie comme suit : 15 11IPT (.;f) A où F4.] est la transformée de Fourier, ryM(r) est la fonction d'autocorrélation de l'élévation à la puissance M du signal échantillonné y(n) que l'on aura éventuellement normalisé et centré au préalable. Lorsque l'ordre M et l'ordre L sont identiques, la fonction MPT correspond à la fonction LSCF prise pour a nul : LSCFy(f, 0) = MPTy(f) \L = M. 20 La formulation générale utilisant la fonction LSCF inclut donc le cas particulier de la fonction MPT. Lorsque l'ordre M et l'ordre L sont identiques, la fonction MPT correspond à la fonction LSCF prise pour a nul : LSCFy(f, 0) = MPTy(f) L = M. La formulation générale utilisant la fonction LSCF inclut le cas particulier de la fonction MPT. 25 La phase d'identification peut s'appuyer sur plusieurs méthodes de classification des pics détectés. Toutes reposent sur les propriétés intrinsèques de chaque fonction LSCF (respectivement MPT) qui génère des raies spectrales caractéristiques du mélange, qu'on peut ainsi qualifier de signature. Les fréquences auxquelles apparaissent ces raies sont des combinaisons des rythmes et des fréquences porteuses résiduelles, alors que leurs profils 30 d'amplitudes selon l'ordre L (respectivement M) des transformations spectrales sont principalement fonction des modulations (constellation et filtre de mise en forme) et des amplitudes des signaux du mélange. 3047578 4 L'invention se fonde sur la richesse des combinaisons de transformations spectrales observables et l'exploitation conjointe des positions et des amplitudes des pics. En effet, le choix de certains ordres de transformation (L ou M) associé à la possibilité d'explorer une dimension supplémentaire de variable (la fréquence cyclique a) rend possible une 5 détermination séparée, sans connaissance a priori, des paramètres constitutifs du mélange de signaux. En particulier, les valeurs d'intérêt pour L (ou M) dépendent des constellations envisagées, et peuvent typiquement être sélectionnées dans l'ensemble {2, 4, 8} pour caractériser les constellations les plus courantes. L'invention est avantageusement complétée par les dispositions suivantes, 10 prisesseules ou en une quelconque de leur combinaison techniquement possible. L'étape de détermination est mise en oeuvre dans un premier temps, en considèrant la fréquence cyclique 0( comme étant la valeur discrète nulle (à savoir la fonction MPT) pour identifier des groupes de triplets {LAp,, LFp,, Lap, = 0}, et dans un deuxième temps, en sélectionnant un ensemble de valeurs de l'ordre L et de valeurs discrètes ou de plages de 15 variation de a#0 déterminés en fonction des résultats de la première application du procédé. La racine Lième du module de la transformée de Fourier est remplacée dans les formules de transformations spectrales par toute autre puissance ou un logarithme ou toute autre fonction réelle monotone du module de la transformée de Fourier. Le procédé comprend en outre un prétraitement du signal reçu comprenant les 20 sous-étapes suivantes : - détermination de la densité spectrale de puissance de y(n) ; - traitement de la densité spectrale de puissance pour estimer le niveau de bruit contenu dans le mélange ; - traitement de normalisation et de centrage du mélange de signaux à partir des 25 caractéristiques précédemment obtenues, la détermination des fonctions LSCF pour certaines valeurs de L prises dans {L1, ... Lq}, étant alors fondée sur le mélange normalisé et centré. L'identification comprend une étape de comparaison de groupes de triplets {LAp,, LFp,, Lap,} obtenus à des signatures prédéterminées à l'avance, c'est à dire des triplets ou 30 couples correspondant à un jeu d'hypothèses de paramètres constitutifs des signaux présents dans le mélange. L'identification permet l'estimation des modulations et des amplitudes des signaux en mettant en oeuvre un algorithme de classification des pics selon la fréquence à laquelle ils sont situés puis une recherche par calcul de paramètres constitutifs du mélange créant des 35 signatures ressemblantes aux groupes ainsi formés. 3047578 5 L'estimation des modulations et des amplitudes du mélange de signaux met en oeuvre un algorithme de recherche de maximum de vraisemblance des pics obtenus utilisant des modèles analytiques de signatures et de probabilité des écarts observés. L'estimation des modulations de chacun des signaux met en oeuvre un algorithme 5 qui recherche le minimum de l'Erreur Quadratique Moyenne (EQM) entre les différents triplets ou couples et une base de signatures préablement établie à partir de valeurs théoriques ou expérimentales. Les valeurs de L sont choisies en fonction des modulations présentes dans des normes de télécommunications considérées. 10 L'identification permet l'estimation des fréquences porteuses résiduelles et des rythmes symboles à partir des valeurs LFp,. L'identification permet l'estimation du nombre de signaux Nu à partir des valeurs LFp,. PRESENTATION DES FIGURES 15 D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : - la figure 1 illustre schématiquement des étapes d'un procédé selon l'invention ; - la figure 2 illustre schématiquement des détails des étapes du procédé de la figure 20 1. Sur l'ensemble des figures les éléments similaires portent des références identiques. DESCRIPTION DETAILLEE DE L'INVENTION Dans le cadre d'un procédé d'estimation de paramètres constitutifs d'un ou de plusieurs 25 signaux partiellement ou totalement superposés en fréquence dans une bande donnée, décrit en relation avec les figures 1 et 2, un signal y(t) est reçu EO par un récepteur 1 radio comprenant une ou plusieurs antennes 2. Le signal composite y(t) reçu est un signal seul ou un mélange de plusieurs signaux contenus dans une bande de fréquences et s'exprime temporellement de la manière 30 suivante : y(t 2rriu(t)-t+0,) u=1 Avec - : le bruit blanc gaussien additif; 3047578 6 - N. : le nombre de signaux contenus dans la bande de fréquences considérée ; - : l'amplitude du signal u; - fu : la fréquence porteuse du signal u; - Ou : la phase à l'origine du signal u ; 5 - su : les symboles émis par l'émetteuru; - Ru : le rythme symboles de l'émetteur u ; - Tu : le retard initial du signal u ; - hu : le filtre de mise en forme du signal u. L'objectif du procédé est d'estimer les paramètres constitutifs de chaque signal 10 présent dans la bande de fréquences considérée. On s'interesse en particulier, pour un signal d'indice u, à son rythme R' à son décalage de fréquence porteuse f' à son amplitude a' à sa constellation C, dont sont issus les symboles su et éventuellement à son filtre de mise en forme de l'impulsion h' Dans une première étape, le signal composite y(t) est numérisé El via une chaine 15 de réception classique. Une telle chaine de réception comprend classiquement un filtre radio-fréquence associé à un amplificateur faible bruit, suivi d'une mélangeur permettant le passage d'une bande comprenant la bande d'intérêt en bande de base, et enfin d'un échantillonneur. On obtient alors un signal échantillonné y(n) où n représente un indice d'échantillon. y(n) est notamment défini par : N y(n) au(n) - ei (2'1 h 20 u=i avec ri, correspondant à la fréquence porteuse résiduelle normalisée du signal u. De manière facultative pour faciliter les traitements ultérieurs, le signal numérisé y(n) est, dans une seconde étape, pré-traité E2. Ce prétraitement E2 comprend une sous-étape de centrage et de normalisation E21. 25 Ce pré-traitement E2 permet d'obtenir un signal ycn(n) tel que sa densité spectrale de puissance est centrée en 0 et tel que : Nü E (1,2 = 1 àu Gin (n) u=1 , avec n=1 Cette sous-étape est éventuellement effectuée à partir de l'estimation E23 du niveau de bruit ac..)2 contenu dans le mélange, et éventuellement, pour chaque signal du mélange, de 30 sa bande utile Bu, ces paramètres étant possiblement estimés E23 au moyen d'une détermination E22 de la densité spectrale de puissance du signal y(n). Ceci est autorisé par 3047578 7 la détection du rapport signal à bruit. La densité spectrale de puissance permet également de déterminer la bande de fréquences utile et la fréquence centrale de cette bande, comme explicité dans le document [Proakis, 2007]. De manière très générale, on détermine dans une troisième étape E3, pour 5 différentes valeurs de L prises dans un ensemble {L1, ... Lq} de q (q strictement positif) nombres réels strictement positifs qui calcule la fonction de corrélation spectrale définie ci-dessous : LSCFy(f, a) Fy RyL (7, a)_ (f) où R 1 yL(n)y*L(n -j2nan est la fonction d'autocorrélation 10 cyclique (en anglais, « cyclic autocorrelation function », (CAF)) du signal échantillonné y(n) élevé à la puissance L que l'on aura éventuellement normalisé et centré, où FI.] est la transformée de Fourier et où, pour chaque q, 0( peut indépendamment prendre une ou plusieurs valeurs discrètes et/ou parcourir une ou plusieurs plages de valeurs. En particulier, une troisième étape E3 simplifiée consiste à considérer tout d'abord 15 une fréquence cyclique a nulle. Ainsi, pour M, nombre réel strictement positif pris dans un ensemble de q (q strictement positif) valeurs prédéterminées {M1, Mq}, cela revient alors à calculer la fonction spectrale : MPTy(f) _ /(7) f où rym(r) est la fonction d'autocorrélation de l'élévation à la puissance M du signal 20 échantillonné y(n) que l'on aura éventuellement normalisé et centré au préalable. En fonction des résultats obtenus, l'analyse peut être complétée par le calcul de fonctions LSCFy(f,aA) limitée n'utilisant qu'un sous-ensemble pertinent de valeurs de L parmi celles de M utilisées au cours de la troisième étape. Les ordres des transformations spectrales M ou L considérés dépendent des 25 modulations envisagées dans le mélange, typiquement M = L = 2, 4, ou 8 sont suffisantes pour caractériser correctement les paramètres des différents signaux pour les modulations courantes telles que BPSK, QPSK, OQPSK, 8-QAM, 16-QAM, 8-PSK, 16-APSK, 32-APSK. Ensuite, dans une quatrième étape, chaque tranformation LSCFy(f,a) obtenue pour chaque valeur de L prise dans l'ensemble {L1, Lq} est traitée E4 pour en extraire des 30 triplets {LAp,, LFp,, Lapa} désignant ses pics p, selon les deux variables f et a où LAp, est l'amplitude du pic, LFp, la fréquence f où il se situe et Lap, la fréquence cyclique 0( à laquelle il apparaît. 3047578 8 Lorsqu'on a employé des valeurs discrètes de 0( pour certaines valeurs de L alors, à chaque valeur de 0( considérée, les pics ne sont détectés que selon la seule variable f et il s'en suit que Lap, vaut a. En particulier, quand on utilise la transformation MPTy(f) seule pour M pris dans {M1, Mq}, ces triplets ont tous une valeur a nulle, ils peuvent donc être 5 réduits à un simple couple {MAp,, MFp,}. On considère qu'un pic est présent dès lors que la fonction spectrale présente un maximum local, valable sur un intervalle d'une largeur qui peut être, selon les mises en oeuvre, déterminée à l'avance ou de manière adaptative, maximum dont l'amplitude est supérieure ou égale à un seuil qui peut, selon les mises en oeuvre, aussi être déterminé à 10 l'avance ou de manière adaptative. Par exemple, la valeur du seuil peut être choisie en fonction du nombre d'échantillons considérés et/ou du niveau de bruit estimé dans le mélange. On remarque que les fonctions LSCFy(f) et MPTy(f) sont extraites par une racine Lième ou Mième du module de la transformée de Fourier. L'utilisation d'une autre puissance 15 de la transformée de Fourier, voire même de son logarithme ou d'une quelconque fonction monotone sur les nombres réels positifs, fournit les mêmes positions de pics avec des amplitudes respectivement modifiées par la fonction employée. Il est donc possible d'utiliser ces formules alternatives pour extraire les triplets. Une fois les différents triplets obtenus, dans une cinquième étape on identifie E5 des 20 groupes de triplets {LAp,, LFp,, Lap,} observés comme des signatures caractéristiques laissées par la présence de certains types de signaux dans la bande de fréquences considérée, c'est-à-dire que ces groupes de pics situés à des fréquences particulières suivent un profil d'amplitude en fonction de l'ordre L qui est caractéristique des paramètres constitutifs du mélange: en particulier, cet ensemble de signatures peut indiquer la 25 puissance et/ou la fréquence porteuse et/ou le rythme symbole et/ou le type de modulation des signaux présents dans le mélange ainsi que leur nombre. Par convention, on considèrera que le signal d'intérêt est celui ayant l'amplitude la plus importante. On notera cependant que certaines applications s'intéressent plutôt à des signaux qui se trouveraient masqués par un ou plusieurs signaux utilisant la même ressource 30 fréquentielle mais avec plus forte puissance. En particulier, il est possible d'utiliser, en vue d'estimer les modulations et les amplitudes de chacun des signaux, différents algorithmes de reconnaissance des signatures les plus probables: - un algorithme de classification de pics apparaissant dans plusieurs transformations à 35 des fréquences reliées entre elles par une relation simple (identiques, harmoniques 3047578 9 ou régulièrement espacées) et ayant pour référence de comparaison une base de données de signatures préalablement établie ; ou un algorithme d'optimisation basé sur des formules analytiques donnant l'amplitude théorique et la relation de position des pics d'une signature. Dans ce cas, 5 les amplitudes et les modulations des signaux sont estimées de manière conjointe sous la forme du problème d'optimisation suivant : = , où est une fonction de coût strictement monotone selon la vraissemblance de l'ensemble [r,..... des hypothèses étant donnés des couples {MAp,, MFp,} ou des triplets {LAp,, LFp,, 10 Lapa} observés, et est l'ensemble des modulations et des amplitudes estimées ; ou bien encore, lorsque seules les modulations sont à estimer, un algorithme basé sur la minimisation de l'EQM (Erreur Quadratique Moyenne) entre les différents triplets {LAp,, LFp,, Lapa} ou couples {MAp,, MFp,} obtenus à partir du signal y(n), et 15 la base de signatures préablement établie à partir de valeurs théoriques ou expérimentales : ref Où { représente l'ensemble des modulations estimées pour les Nu Cl... signaux présents dans la bande de fréquences considérée, ref représente 20 la base de références des signatures, c'est à dire les valeurs des triplets {LAp,, LFp,, Lapa} établie préalablement pour différentes valeurs des paramètres constitutifs des modulations considérées en hypothèse { Nu } et fi représente l'ensemble des valeurs des triplets {LAp,, LFp,, Lapa} obtenus à partir du signal y(n). Similairement, les fréquences porteuses résiduelles et les rythmes symboles moyens 25 sont estimés à l'aide d'une méthode basée sur l'observation de l'ensemble des valeurs et LFp, (respectivement MFp,). Conjointement, le nombre de signaux Ni, est éventuellement déterminé. Références : 30 - [Guenther, 2010] : Guenther, B. (2010). Multi-User Signal Classification Via Cyclic (MsC Thesis) ; 3047578 10 - [Mazo, 1978] : Mazo, J. (1978). Jitter Comparison of Tones Generated by Squaring and by Fourth-Power Circuits. The Bell System Technical Journal ; - [Proakis, 2007] : Proakis, J. (2007). Digital Communications 5th Edition. McGraw-Hill Education ; 5 - [Reichert Reichert, 1992] : J. (1992). Automatic Classification of Communication Signal using Higher Order Statistics. IEEE International Conference on Speech, Acoustics and Signal Processing ; - [Swami & Sadler, 2000] : Swami, A., & Sadler, B. (2000). Hierarchical digital modulation classification using cumulants. IEEE Transactions on Communications ; 10 - [Zaerin & Seyfe] : Zaerin, M., & Seyfe, B. (2012). Multiuser modulation classification based on cumulants in additive white Gaussian noise channel. IET Signal Processing. REVENDICATIONS1. Procédé d'estimation de paramètres d'un signal seul ou de signaux mélangés en fréquence dans une bande donnée, le procédé comprenant les étapes suivantes : - réception (EO) d'un signal seul ou d'un mélange de signaux y(t) dans une bande de fréquences donnée ; - détermination (E3) pour L, nombre réel strictement positif pris dans un ensemble de q, q strictement positif, valeurs prédéterminées {L1; ;Lq} de la fonction de corrélation spectrale définie ci-dessous : LSCFy(f, a) RyL (T, a)_ (f) où r^=0 est la fonction d'autocorrélation cyclique de l'élévation à la puissance L de la version échantillonné y(n) du signal y(t) et où FI.] est la transformée de Fourier et où, pour chaque q, a peut indépendamment prendre une ou plusieurs valeurs discrètes et/ou parcourir une ou plusieurs plages de valeurs ; - traitement (E4) de chaque tranformation LSCFy(f,a) L étant pris dans un ensemble de q, q strictement positif, valeurs prédéterminées {L1; ... ;Lq} pour en extraire un ensemble de pics p,, c'est-à-dire de maxima locaux détectés selon les variables f et a ou selon la seule variable f pour chaque valeur discrète de a, chacun de ces pics étant caractérisé par un triplet {LAp,, LFp,, Lapa} où LFp, et Lap, sont respectivement la fréquence et la fréquence cyclique à laquelle il se situe et LAp, son amplitude ; - identification (E5) de groupes de triplets {LAp,, LFp,, Lap,} comme étant la signature laissée par un ensemble de paramètres constitutifs de signaux présents dans la bande de fréquence donnée. 2. Procédé d'estimation selon la revendication 1, dans lequel l'étape de détermination (E3) est mise en oeuvre dans un premier temps, en considèrant la fréquence cyclique OC comme étant la valeur discrète nulle pour identifier des groupes de triplets {LAp,, LFp,, 0}, et - dans un deuxième temps, en sélectionnant un ensemble de valeurs de l'ordre L et de valeurs discrètes ou de plages de variation de a#0 déterminés en fonction des résultats de la première application du procédé. 3047578 12 3. Procédé selon l'une des revendications 1 ou 2, dans lequel la racine Lième du module de la transformée de Fourier est remplacée dans les formules de transformations spectrales par toute autre puissance ou un logarithme ou toute autre fonction réelle monotone du module de la transformée de Fourier. 5 4. Procédé d'estimation selon l'une des revendications 1 à 3, comprenant en outre un prétraitement du signal reçu comprenant les sous-étapes suivantes : - détermination de la densité spectrale de puissance de y(n) ; - traitement de la densité spectrale de puissance pour estimer le niveau de bruit 10 contenu dans le mélange ; - traitement de normalisation et de centrage du mélange de signaux à partir des caractéristiques précédemment obtenues, la détermination des fonctions LSCF pour certaines valeurs de L prises dans {L1, ... Lq}, étant alors fondée sur le mélange normalisé et centré. 15 5. Procédé d'estimation selon l'une des revendications 1 à 4, dans lequel l'identification (E5) comprend une étape de comparaison de groupes de triplets {LAp,, LFp,, Lapa} obtenus à des signatures prédéterminées à l'avance, c'est à dire des triplets correspondant à un jeu d'hypothèses de paramètres constitutifs des signaux présents dans le mélange. 20 6. Procédé d'estimation selon l'une des revendications 1 à 5, dans lequel l'identification (E5) permet l'estimation des modulations et des amplitudes des signaux en mettant en oeuvre un algorithme de classification des pics selon la fréquence à laquelle ils sont situés puis une recherche par calcul de paramètres constitutifs du mélange créant des signatures 25 ressemblantes aux groupes ainsi formés. 7. Procédé d'estimation selon la revendication 5 ou 6, dans lequel l'estimation des modulations et des amplitudes du mélange de signaux met en oeuvre un algorithme de recherche de maximum de vraisemblance des pics obtenus utilisant des modèles analytiques 30 de signatures et de probabilité des écarts observés. 8. Procédé d'estimation selon la revendication 5 ou 6, dans lequel l'estimation des modulations de chacun des signaux met en oeuvre un algorithme qui recherche le minimum de l'Erreur Quadratique Moyenne (EQM) entre les différents triplets ou couples et une base de signatures préablement établie à partir de valeurs théoriques ou expérimentales. 3047578 13 9. Procédé d'estimation selon l'une des revendications précédentes, dans lequel les valeurs de L sont choisies en fonction des modulations présentes dans des normes de télécommunications considérées. 5 10. Procédé d'estimation selon l'une des revendications 1 à 9, dans lequel l'identification (E5) permet l'estimation des fréquences porteuses résiduelles et des rythmes symboles à partir des valeurs LFp,. 10 11. Procédé d'estimation selon l'une des revendications 1 à 10, dans lequel l'identification (E5) permet l'estimation du nombre de signaux Ni, à partir des valeurs LFp,.