FR 2806499 A1 20010921 FR 0003597 A 20000320 La présente invention concerne un procédé d'estimation d'une matrice de corrélation de signaux à caractéristiques inconnues, reçus par un réseau comportant un nombre N déterminé de capteurs, les capteurs recevant des signaux composés notamment de signaux ayant des caractéristiques connues ainsi que de signaux de caractéristiques inconnues telles que du bruit et/ou des brouilleurs. L'invention trouve son application dans différents domaines, notamment les suivants: * la goniométrie de différentes sources arrivant sur le réseau, c'est-à-dire la détermination de leurs directions d'arrivée, * le filtrage spatial afin de protéger une ou plusieurs liaisons de sources interférentes, * tout autre dispositif dans lequel il est souhaitable d'estimer la matrice de corrélation de signaux utiles et celle des signaux interférents. Dans de nombreuses applications notamment dans le domaine des télécommunications, I'utilisation des méthodes de traitement d'antenne se révèle particulièrement intéressante. Parmi les signaux arrivant sur le réseau de capteurs qui constituent l'antenne, certains sont de caractéristiques connues, les signaux utiles en télécommunications par exemple, d'autres présentent des caractéristiques inconnues, tels les interférences. Un grand nombre de ces méthodes sont basées sur l'estimation de la matrice de corrélation des signaux reçus par les capteurs Dans la suite de la description différents paramètres et notations sont utilisés, notamment les suivantes: * N. le nombre de capteurs d'un réseau, * x(n), le vecteur observation reçu par le réseau de N capteurs après numérisation, ou signaux capteurs, * xu(n), la composante du vecteur observation correspondant aux signaux dits signaux utiles, dont les caractéristiques sont connues du traitement, * M, le nombre de signaux utiles, * xb(n), la composante du vecteur observation correspondant aux signaux dits signaux interférents (bruit et/ou brouilleurs), dont les caractéristiques sont inconnues du traitement, * P. le nombre de signaux interférents, * Rx = E(x(n) x(n)t), la matrice de corrélation théorique des signaux capteurs, * Rs, la matrice de corrélation théorique du vecteur signal xu(n), * Rb, la matrice de corrélation théorique du vecteur signal xb(n). *,t, représente l'opération de transposition -conjugaison. En supposant que les signaux utiles sont décorrélés des signaux interférents, la relation (1) liant les différentes matrices de corrélation théoriques définies ci-dessus, peut s'exprimer de la manière suivante: Rx = E(x(n) x(n)t) = E((x,(n)+xb(n)) (xU(n)+xb(n))t) Rx = E(x, (n) x,(n)t) + E(xb(n) xb(n)t) + 2 Re{E(x,(n) xb(n)')} Rx= Rs+ Rb (1) En disposant de la connaissance de la matrice de corrélation des signaux utiles Rs et de la matrice Rx, la matrice Rb peut s'obtenir simplement à partir de la relation (1): Rb=Rx-Rs. (1') Appliquée aux estimées de ces matrices de corrélation, la relation (1') ne donne pas de bons résultats, comme il est démontré par exemple dans le livre " Introduction to adoptive arrays " ayant pour auteurs RA.Monzingo, W.Miller, publié aux éditions John Wiley and Sons, New York, 1980. En effet, lI'estimée de la matrice de corrélation Rx comporte également des termes relatifs à la corrélation entre les signaux utiles et les signaux interférents E(xu(n) xb(n)t), qui ne sont nuls qu'asymptotiquement. L'objet de l'invention concerne un procédé permettant d'estimer notamment la matrice de corrélation des signaux interférents à partir de la matrice de corrélation des signaux reçus sur un réseau comportant un nombre N de capteurs et d'une estimée de la matrice de corrélation des signaux utiles. L'invention concerne un procédé pour estimer la matrice de corrélation de signaux de caractéristiques inconnues dans un réseau comportant un nombre N de capteurs, lesdits capteurs recevant au moins les signaux Sb de caractéristiques inconnues et des signaux Su de caractéristiques connues, I'ensemble des Su et Sb formant les signaux capteurs Sx. Il est caractérisé en ce qu'il comporte au moins les étapes suivantes: a) déterminer la matrice /s ou Rs,(x de corrélation des signaux Su de caractéristiques connues, b) estimer la matrice de corrélation tx des signaux Sx capteurs, A AA c) former une matrice A égale à A = Rsi'2 Rx Rs"/2 d) décomposer la matrice A en éléments propres pour obtenir ses valeurs propres X; et ses vecteurs propres ui, avec I e) construire une matrice B ayant des vecteurs propres v; sensiblement égaux aux vecteurs propres ui de la matrice A, et des valeurs propres, égales aux valeurs propres X, à une constante 3 près, la valeur de la constante 3 étant déterminée à partir des N-P plus petites valeurs propres de la matrice A, P étant le nombre de brouilleurs ou signaux interférents, avec A, 5 f) déterminer la matrice estimée lR*b = /s2 B S12. Pour une matrice B de rang N-P, la valeur de la constante C est par exemple la plus petite valeur propre des Xi ou bien elle est déterminée à partir des N-P plus petites valeurs propres des Ri en prenant la valeur moyenne. Le procédé selon l'invention s'applique aussi au traitement d'antennes. L'invention concerne aussi un dispositif pour estimer la matrice de corrélation de signaux de caractéristiques inconnues dans un réseau comportant un nombre N de capteurs, lesdits capteurs recevant au moins les signaux Sb de caractéristiques inconnues et des signaux Su de caractéristiques connues, I'ensemble des Su et Sb formant les signaux capteurs Sx. Il est caractérisé en ce qu'il comporte des moyens adaptés pour mettre en oeuvre le procédé selon l'une des caractéristiques décrites ci-dessus. L'invention présente notamment les avantages suivants: * dans les applications par goniométrie, elle permet de diminuer le nombre de capteurs utilisés, la goniométrie étant effectuée uniquement sur certains signaux, et d'améliorer la précision de goniométrie des sources sélectionnées, * dans les applications OLS abréviation de l'expression "opposition des lobes secondaires" ou SLC en anglais pour (side lobe canceller), la minimisation de puissance recherchée peut être réalisée uniquement sur les signaux interférents tout en protégeant les signaux utiles, * dans les applications FAS abréviation de l'expression (filtre adapté spatial) ou SMF (Spatial Matched Filter), le procédé selon l'invention permet d'augmenter la vitesse de convergence du traitement et la robustesse par rapport aux erreurs. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite, à titre illustratif et nullement limitatif, en regard des dessins annexés qui représentent: ò la figure 1 schématise les différents signaux arrivant sur un réseau comportant un nombre N de capteurs, la figure 2 représente un algorithme comportant les étapes du procédé selon l'invention, * la figure 3 schématise un dispositif OLS, et * laI figure 4 représente une application du procédé selon l'invention au filtrage adapté spatial. Le procédé selon l'invention utilise notamment les estimées de la matrice de corrélation des signaux utiles et de la matrice de corrélation des signaux capteurs pour en déduire l'estimée de la matrice de corrélation des signaux interférents. La figure 1 schématise une antenne comportant un réseau de N capteurs Ci. Les capteurs Ci recevant tous les signaux utiles, bruit, brouilleurs, sont en liaison chacun avec un récepteur, l'ensemble des N récepteurs Di constituant un récepteur multivoies connus de l'Homme du métier. Le récepteur multivoie est lui-même relié à un dispositif 1 de numérisation et de traitement des signaux en réalisant les étapes du procédé selon l'invention. Pour cela le dispositif comporte par exemple un microprocesseur 2 ou un calculateur adapté pour mettre en oeuvre les étapes. Le procédé selon l'invention repose sur des concepts rappelés en utilisant les matrices théoriques ci-après, à partir de la relation Rx = Rs + Rb (1) Par Rs on désigne la matrice Rs qui peut être connue ou de manière plus générale Rs,ac la matrice Rs connue au coefficient a près, Rs,a = aRs L'idée consiste à a) " blanchir " la matrice d'observation Rs de la façon suivante: 1 - si la matrice Rs est inversible, le procédé consiste à multiplier à droite et à gauche les deux termes de l'équation (1) par (Rs,c)4"2, I'expression (2) suivante est obtenue: (Rs,aX)l/2 Rx(Rs,a()"2 = (Rs,a(X)112 RS (Rs,af)"/2 + (Rs,a)'112 Rb (Rsa)'112 (Rs, o)- 1/2 Rx (Rs,)-1/2 = ' Id + (Rs, ay/2 Rb (Rs, c-112 (2) ou encore A = c' Id + B avec Id la matrice identité 2 - si la matrice Rs n'est pas inversible, ce qui est le cas en particulier si le nombre M de signaux utiles est inférieur au nombre N de capteurs du réseau N. La matrice est rendue inversible en ajoutant un terme 6 Id, avec 6 suite des étapes du procédé correspond alors à Rs,c(x = c Rs + 8 Id. La matrice de corrélation Rb calculée par la suite comprendra la contribution intégrale des signaux interférents, et une contribution réduite du bruit de fond (de puissance 2 - 6), ce qui revient à augmenter artificiellement le INR (Interference to Noise ratio) des signaux interférents dans la matrice de corrélation estimée. Ainsi qu'il a été mentionné précédemment, I'expression (2) peut être de la forme A = cc Id + B b) décomposer en éléments propres matrice A blanchie obtenue au cours de l'étape a) N A = Rs'112 Rx Rs'1/2 = X Ui ui (3) i=1 O: Xi et ui représentent respectivement les valeurs propres et les vecteurs propres de la matrice A. A la fin de l'étape b) les valeurs propres xi et les vecteurs propres ui de la matrice A sont connus, avec l propres obtenus à l'étape b). Dans le cas o le nombre de signaux interférents est strictement inférieur au nombre N de capteurs, et en l'absence de bruit de fond, alors la matrice B = (Rs, O'-1/2 Rb (Rs,c)'1/2 est de rang N-P, (N étant le nombre de capteurs et P le nombre de signaux interférents) et d'après l'équation (2), la matrice A possède P valeurs propres supérieures à c-' et N-P valeurs propres identiques telles que: X1 = X2 =... = kN-P = o(x. Dans ces conditions, la matrice B élaborée a pour éléments propres: * des vecteurs propres vi qui sont sensiblement identiques aux vecteurs propres ui de la matrice A, et a des valeurs propres yi égales aux valeurs propres ci - X1 La matrice B s'exprime par exemple par la relation N B= (xi- X 1) uji U(4) i = 1 Avec yi = (Xi- X1) d) connaissant B et la matrice de corrélation Rs, il est possible d'obtenir la matrice de corrélation des signaux interférents Rb = Rs;12 B Rs12 (5) Les relations (2) à (5) précédentes, écrites à partir des matrices de corrélation théoriques, servent de base pour les étapes du procédé selon l'invention qui utilisent les estimées des matrices. Les étapes exposées ci-dessus utilisant les matrices de corrélation théorique, sont appliquées aux estimées de matrice. Cela implique la connaissance des estimées des matrices de corrélation des signaux utiles et des signaux capteurs pour déterminer l'estimée de la matrice des signaux interférents. Pour cela l'information disponible sur les signaux utiles doit permettre de construire une estimation suffisamment bonne de la matrice de corrélation Rs, à une constante multiplicative près. Cette estimation peut être obtenue de plusieurs façons décrites ci-après. Le nombre P de signaux interférents est de préférence inférieur au nombre N de capteurs constituant l'antenne de réception. Si cette dernière condition n'est pas vérifiée, un biais sera obtenu sur l'estimation de la matrice Rb. Le procédé met en oeuvre plusieurs étapes mentionnées sur l'algorithme de la figure 2. A) Déterminer la matrice blanchie A par les étapes suivantes A 1 ère étape - 10 - estimer la matrice de corrélation Rs à partir des informations sur les sigaux utiles. Par exemple, si le traitement connaît les vecteurs directeurs (qui peuvent être déterminés à partir des directions d'arrivée) et les puissances A des signaux utiles, la matrice de corrélation estimée Rs s'écrit: M Rs = Snm Smt (6) m= 1 avec 7m = la puissance et Sm = vecteur directeur Si la matrice de corrélation ainsi formée n'est pas inversible, la matrice est rendue inversible comme il a été exposé précédemment par ajout d'un terme 8 Id: M Rs = n7m Sm smt + Id (7) m= selon le cas, I'expression (6) ou l'expression (7) est utilisée dans la suite du procédé 2 ème étape - 20 estimation de la matrice de corrélation des signaux capteurs reçus sur K échantillons, à partir par exemple de I'estimateur classique non biaisé: K Rx = Y x(k) x(k)t (8) Kk=1 L'intégration de la matrice de corrélation Rx doit être effectuée sur un nombre d'échantillons K suffisant, afin de pouvoir négliger les termes de corrélation entre les signaux utiles et les signaux interférents. 3 ème étape - 30 - calcul de la matrice A = Rs"2 Rx Rs"/2 et décomposition de la matrice A en éléments propres N A = i ui uit (9) 1=1 avec R,'"2 = inverse racine carrée de R, Le calcul et la décomposition sont décrits par exemple dans le livre ayant pour titre "Numerical Recipes in C" et pour auteurs W.H. PRESS, B.R FLANNERY, S.A. TEUKOLSY, W.T. VET1ERLING, 1992 (Editeur: Cambridge University). A l'issue de cette 3 éme étape les vecteurs propres ui et les valeurs propres I, de la matrice A sont connus B) Déterminer une matrice B et l'estimée de la matrice des signaux interférents 4 ème étape - 40 formation de la matrice B à partir des éléments de A en déterminant ses éléments propres: vecteurs propres et valeurs propres * les vecteurs propres vi sont sensiblement égaux ou égaux aux vecteurs propres ui obtenus par décomposition de la matrice A, a les valeurs propres yi sont sensiblement égales ou égales aux valeurs propres Xi obtenues par décomposition de la matrice A corrigées d'une valeur ou constante C qui est choisie à partir des plus petites N-P valeurs propres de la matrice A ou égale à la plus petite valeur propre de la matrice A. la matrice B prend alors l'expression suivante N Y_ (Xi -) Ui UJt (1 0) i=1 Avantageusement, la valeur de la constante D est la plus petite valeur propre des Xi Selon un autre mode de réalisation, la valeur de la constante 3 est déterminée en moyennant les N-P valeurs propres de la matrice A. éme étape - 50 - déterminer la matrice de corrélation des vecteurs signaux interférents. A A partir de B et de la valeur estimée pour la matrice Rs, le procédé ^ permet d'obtenir la valeur estimée de la matrice de corrélation Rb correspondant aux signaux interférents. A A 2 A 2 Rb = Rs'2 B Rs'2 (11) La connaissance de la valeur estimée de la matrice de corrélation des signaux interférents peut être utilisée dans différents traitements de signaux dont certains exemples sont donnés ci- après à titre illustratif et nullement limitatif. Le procédé selon s'applique dans différents types d'applications dont certaines sont données à titre illustratif et nullement limitatifs. Application à la goniométrie Plusieurs méthodes de goniométrie sont basées sur l'exploitation de la matrice de corrélation des signaux reçus sur un réseau de capteurs. Certaines de ces méthodes, dites méthodes à haute résolution telle la méthode MUSIC décrite par exemple dans la référence (1) ou la méthode ESPRIT décrite dans la référence (2) par exemple. Ces méthodes permettent de faire la goniométrie de plusieurs sources arrivant simultanément sur un réseau de capteurs. (1) article intitulé " Multiple Emitter Location and Signal Parameter Estimation " de la revue IEEE Trans. Ant. Prop., Vol AP-34, n 3, pp 276280, March 1986 ayant pour auteur M.R.D.Schmidt. (2) article intitulé "Estimation of Signal Parameters via Rotational Invariance Techniques " de la revue IEEE Transaction ASSP, Vol ASSP-37, pp 984-995, July 1989 ayant pour auteurs MM.Roy, T. Kailath. Toutefois, ces méthodes possèdent certaines limitations: I'utilisation d'un réseau de N capteurs permet ainsi typiquement d'effectuer la goniométrie de N/2 sources. Lorsque le réseau reçoit des signaux utiles et des signaux interférents, il est possible de les séparer, en mettant en oeuvre les étapes du procédé décrites ci-dessus. Après avoir séparé les signaux interférents des signaux utiles, la goniométrie est par exemple effectuée selon des méthodes connues de l'Homme du métier, en mettant en oeuvre uniquement la A matrice de corrélation estimée Rb. La goniométrie détermine les directions d'arrivée des P signaux interférents, et non pas les directions d'arrivée de M+P sources, ce qui améliore ses performances. Elle est réalisée comme si les signaux utiles étaient absents. L'application du procédé selon l'invention nécessite un nombre de capteurs moindre que celui habituellement utilisé dans l'art antérieur. Application au filtrage spatial OLS Le schéma de principe de l'OLS est rappelé sur la figure 3. L'OLS est utilisé pour protéger les signaux utiles arrivant d'une direction donnée vis-à- vis de signaux interférents, Il est décrit par exemple dans les références (3) et (4) suivantes: (3) article intitulé "Adaptative Arrays " de la revue IEEE Trans. Ant.Prop., Vol AP 24, N 5, pp 585-598, September 1976 ayant pour auteur S.P.Applebaum,. (4) article intitulé " Explorations in fixed and adaptive resolution at GE and SURC ", de la revue IEEE Trans. Ant.Prop., Vol AP 24, N 5, pp 585-598, September 1976 ayant pour auteur P.W.Howells. Il est basé sur l'utilisation d'un capteur principal directif 60, pointant dans la direction des signaux utiles, et de L capteurs auxiliaires 61 permettant d'éliminer l'influence des signaux interférents présents dans la voie principale. Les L capteurs auxiliaire et le capteur principal constitue le réseau de N capteurs. Les pondérations wo de l'OLS sont calculées de manière à minimiser la puissance totale du signal de sortie yp(f - wot z(O. Elles s'expriment à partir de Rz, la matrice de corrélation des signaux des voies auxiliaires et de rz, le vecteur d'inter-corrélation entre z(f et yp(f: w = Rz' rzy (12) avec Rz =E(z(t, z(t) t) et rzy: E(z(), yp(f t) o z(t) est le signal de sortie des voies auxiliaires et yj(1 le signal de sortie de la voie principale. L'OLS doit être employé avec précaution, afin d'éviter le risque d'une élimination des signaux utiles. Ce risque est d'autant plus fort que les brouilleurs sont de puissance faible par rapport aux utiles, ou que le nombre de voies auxiliaires est supérieur au nombre de brouilleurs. Un bon conditionnement de l'OLS doit être mis en oeuvre. Ainsi les voies auxiliaires doivent être aussi peu directives que possible dans la direction des signaux utiles et aussi directives que possible dans la direction des brouilleurs. Afin d'éviter le risque d'une élimination des signaux utiles, plusieurs méthodes de robustification de 'OLS existent dans la littérature: * Ajout d'un bruit fictif de manière à "( masquer " la présence des signaux utiles, * Contrainte de norme sur les pondérations Wa: I I w, I I Le procédé selon l'invention consiste à mettre en oeuvre l'OLS comme si les signaux utiles étaient absents, en estimant les corrélations de l'expression (12) correspondant aux signaux interférents. Le risque de leur élimination est ainsi écarté et l'algorithme optimise l'élimination des signaux interférents et non des signaux utiles, A Pour cela la matrice de corrélation Rb est estimée en mettant en oeuvre les étapes du procédé selon l'invention et en posant x(X = (z()T yp(O)), La formule (12) permet alors d'obtenir les pondérations w, en considérant que la matrice Rz (respectivement le vecteur rzy) correspond aux L premières lignes et aux L premières colonnes (respectivement aux L premières lignes de A la dernière colonne) de la matrice estimée Rb. Filtre Adapté Spatial - figure 4 Afin de protéger un signal utile dont la direction d'arrivée est connue, il est possible de mettre en ceuvre le filtre adapté spatial sur direction d'arrivée. L'algorithme utilisé consiste à minimiser la puissance totale de sortie de l'antenne tout en maintenant un gain unité dans la direction du signal utile. Les notations suivantes sont utilisées * w, le vecteur de pondération correspondant au filtre spatial, * s, le vecteur directeur du signal utile, qui s'obtient par exemple à partir de la connaissance de la direction du signal utile, connu ou estimé par goniométrie par exemple, * y(t) = wt x(t9 la sortie du filtre spatial. Le vecteur w s'écrit: w = a Rx' s, o a est une constante multiplicative près. L'obtention du vecteur w est par exemple décrit dans la référence (5) ayant pour titre "Système mixte de goniométrie à haute résolution et de réfjection adaptative de brouilleurs dans la gamme HF", ayant pour auteurs F.Pipon, P.Chevalier, G.Multedo publié dans Proc.GRETSI, pp 685- 688, Septembre 1991. Le filtre spatial conduit à une maximisation du SNIR: rapport entre la puissance du signal arrivant de la direction correspondant au vecteur directeur s (signal utile) et la somme de la puissance du bruit de fond et de la puissance des signaux venant des autres directions (signaux interférents). Dans la pratique, w est estimé en utilisant la matrice de corrélation A A estimée Rx, et un vecteur s proche du véritable vecteur directeur s. Cette technique exige que le vecteur directeur s soit connu A parfaitement. En effet, lorsque s ne correspond pas au vecteur directeur s, l'algorithme risque d'éliminer le signal utile. Le procédé selon l'invention permet de s'affranchir de ce risque et d'améliorer les performances du filtre adapté spatial. A Le filtre spatial w est calculé à partir de la matrice estimée Rb, déterminée en mettant en oeuvre le procédé selon l'invention calculé selon A A I'invention, et de s. Le filtre spatial w est alors égal à a Rb- s. L'utilisation de cette méthode permet également avantageusement d'améliorer la vitesse de convergence de l'algorithme de calcul de w. Sans sortir du cadre de l'invention, le procédé s'applique notamment dans les domaines suivants: * quelque soit le type de sources: bande étroite, large bande, ou trajets multiples, * télécommunications spatiales, par exemple pour la protection d'un théâtre d'opérations vis-à-vis d'interférences, * écoute passive, par exemple après une étape de goniométrie et une intensimétrie des émetteurs reçus, * GE: localisation angulaire d'interférences en présence de signaux connus, 25. Radar, acoustique, radiocommunications. REVENDICATIONS 1 - Procédé pour estimer la matrice de corrélation de signaux de caractéristiques inconnues reçus par un réseau comportant un nombre N de capteurs, afin de traiter les signaux reçus par ledit réseau, lesdits capteurs recevant au moins les signaux Sb de caractéristiques inconnues et des signaux Su de caractéristiques connues, I'ensemble des Su et Sb formant les signaux capteurs Sx caractérisé en ce qu'il comporte au moins les étapes suivantes: A a) déterminer la matrice Rs ou Rs,a de corrélation des signaux Su de caractéristiques connues, A b) estimer la matrice de corrélation Rx des signaux Sx capteurs, c) former une matrice A égale à A = Rs'2 Rx Rs"2 d) décomposer la matrice A en éléments propres pour obtenir ses valeurs propres X. et ses vecteurs propres ui, avec l la valeur moyenne. A 4 - Procédé selon la revendication 1 caractérisé en ce que la matrice Rs déterminée à l'étape a) est déterminée à partir d'informations de puissance et de directions d'arrivée des signaux utiles selon l'expression: M Rs = nm Sm Sm (6) m = 1 avec -m = la puissance et sm = vecteur directeur - Procédé selon l'une des revendications 1 à 4 caractérisé en ce que la matrice Rs ou Rs,ox, est rendue inversible en ajoutant un terme 6 Id, avec 6 G2, o a2 est la puissance de bruit de fond. 6 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce qu'il comporte une étape de goniométrie effectuée sur les signaux interférents Sb A en utilisant la matrice estimée Rb obtenue au cours de l'étape f). 7 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce qu'il comporte une étape de détermination des pondérations wo dans un dispositif de type OLS o les N capteurs comportent L capteurs auxiliaires (61) et un capteur principal (60), en déterminant une matrice Rz correspondant A aux L premières lignes et aux L premières colonnes de la matrice estimée Rb et un vecteur Rzy correspondant aux L premières lignes de la dernière A colonne de la matrice estimée Rb et en utilisant l'expression wo = Rz' rZy avec Rz =E(z(t). z(t) t) et rzy = E(z(t), yp(t) t) 8 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce qu'il comporte une étape de détermination d'un filtre spatial w à partir de la A A matrice estimée Rb avec w = a Rb-'s. 9 - Application du procédé selon l'une des revendications 1 à 8 au traitement d'antennes. - Dispositif pour estimer la matrice de corrélation de signaux de caractéristiques inconnues reçus par un réseau comportant un nombre N de capteurs, afin de traiter les signaux reçus par ledit réseau, lesdits capteurs recevant au moins les signaux Sb de caractéristiques inconnues et des signaux Su de caractéristiques connues, I'ensemble des Su et Sb formant les signaux capteurs Sx caractérisé en ce qu'il comporte des moyens adaptés pour mettre en oeuvre le procédé selon l'une des revendications 1 à 5.