FR 3011086 A1 20150327 FR 1302189 A 20130920 PROCEDE REALISANT CONJOINTEMENT LA SYNCHRONISATION, L'IDENTIFICATION, LA MESURE, L'ESTIMATION DU FILTRE DE PROPAGATION ET LA LOCALISATION D'EMETTEURS UTILES ET INTERFERANTS La présente invention concerne un procédé permettant de déterminer la position d'un émetteur et ses paramètres identifiants en utilisant la connaissance a priori de la forme d'onde. Elle est notamment utilisée pour réaliser une métrologie des signaux et canaux de propagation électromagnétiques ou acoustiques, lorsque celle ci fait intervenir des émetteurs sur la forme d'onde desquels une connaissance a priori est disponible, des capteurs et réflecteurs fixes ou mobiles, au moyen d'un critère conjoint portant sur des caractéristiques connues de la forme d'onde et sur des variables :espacedirection-d'arrivée/retard-distance/doppler-cinématique dénommées CCEA-RD-DC dans la suite, pour la détection et la détermination des caractéristiques de transmission et de réception des signaux, et, le cas échéant, pour la détermination de paramètres de position et de cinématique des émetteurs tels que les angles d'arrivée, le retard et le Doppler, le retard différentiel et le Doppler différentiel, etc. L'invention s'applique pour des systèmes actifs, radars, sonars, sondeur de canal, pseudolite en radiolocalisation, etc., réseau de communication et de réseau de radiolocalisation. Les systèmes conventionnels de métrologie des signaux radioélectriques et acoustiques, de mesure d'occupation spectrale, d'identification des émetteurs radioélectriques et acoustiques, de sondage de canaux de propagation, de détection et localisation d'interférences ou de brouillages, sont le plus souvent génériques et utilisent en général : des capteurs passifs ou algorithmes spécifiques à la détection et à l'identification des signaux radioélectriques ou acoustiques, que ces signaux soient considérés comme utiles ou interférents, des capteurs passifs ou algorithmes spécifiques à l'estimation de localisation des émetteurs des signaux radioélectriques ou acoustiques (que ces émetteurs soient considérés comme utiles ou interférents), par mesure d'angle d'arrivée (goniométrie), par mesure de retard ou de retard différentiel, par mesure de Doppler ou de Doppler différentiel, etc., des capteurs ou algorithmes spécifiques à l'estimation des caractéristiques de propagation (réponse impulsionnelle du canal de propagation (dénommé CIR pour « Channel Impulse Response »), soit passifs et exploitant alors les émetteurs déjà présents dans l'environnement, soit actifs (sondeurs de canal) avec utilisation d'un ou plusieurs émetteurs spécifiques). L'art antérieur décrit différentes structures de réception et traitement correspondantes. En général, la structure de réception d'un détecteur utilise des informations ou des modèles sur les canaux de propagation des signaux émanant des émetteurs et réflecteurs à détecter et sur le bruit global, lequel comprend le bruit thermique des récepteurs et les interférences potentielles. Les récepteurs conventionnels utilisés en détection électromagnétique ou acoustique, scrutant ou non l'espace par balayage électronique ou mécanique d'un faisceau, mettant en oeuvre ou non une analyse spatiale en amont et découplée de l'analyse Distance/Doppler, supposent, en général, implicitement ou explicitement: - des positions spatiales ou des directions d'arrivée des signaux en provenance des émetteurs à détecter et à mesurer, quitte à organiser le cas échéant une scrutation de l'espace angulaire ou un maillage spatial du terrain observé, - des cases Distance/Doppler, quitte à organiser le cas échéant une scrutation desdites cases Distance/Doppler, - un modèle de propagation, en espace libre ou autre, - un modèle de bruit global, le plus souvent gaussien, circulaire et de puissance inconnue, devant être calibré préalablement ou estimé in situ, - les signaux émanant de réflecteurs ou de multitrajets de propagation faibles devant les signaux principaux émanant des émetteurs les plus forts par trajet direct, - une phase inconnue des signaux émanant des émetteurs et des réflecteurs. Ces systèmes appliquent implicitement un découplage entre les analyses conduites à l'aide de différents senseurs et algorithmes. Ce découplage limite les performances du système global à la borne supérieure des performances individuelles de chaque capteur ou algorithme, ce qui rend lesdites performances globales tributaires des limitations intrinsèques à chaque senseur/algorithme ou à des conditions particulières de mise en oeuvre, sans renforcement possible par une exploitation des synergies entre les différents senseurs/algorithmes. Au contraire, la présente invention déroule des analyses conjointes qui se renforcent mutuellement en faisant jouer au maximum les synergies entre plusieurs algorithmes mis en oeuvre, typiquement les suivants : - détection et synchronisation en temps (ou retard) et en fréquence (ou Doppler) sur les signaux, - estimation des directions d'arrivées des émetteurs des signaux, estimation des retards de propagation et des retards différentiels de propagation entre émetteurs des signaux, estimation des Doppler différentiels entre émetteurs des signaux, algorithmes de reconnaissance des signaux et des émetteurs basés sur des analyses de structure des signaux : o caractéristiques d'enveloppe dans le domaine temps fréquence, o caractéristiques de stationnarité et de continuité, de périodicité et de récurrences des signaux impulsifs ; o dans le cas de signaux de communications, caractéristiques de trame, porteuses et canalisation en fréquence, rythmes et périodes symboles, etc. analyses des caractéristiques fines des signaux à des fins d'identification des émetteurs: o modulation et codages des signaux de communication, o jitter et structure interne des pulses radar, o etc., mesures de réponses impusionnelles de propagation (CIR), - etc. Définitions utilisées Pour un émetteur et une réception du signal de cet émetteur sur un réseau antennaire les notations suivantes seront utilisées : s(kTe) signal émis à l'instant k.Te, comprenant le signal de référence d(kTe), où Te est la période sur laquelle les signaux sont échantillonnés avec 01(e Le procédé selon l'invention repose notamment sur l'utilisation d'une information a priori contenue dans les signaux émis par des émetteurs et d'un traitement permettant d'obtenir conjointement au moins une information pour localiser les émetteurs et les identifier. Il consiste notamment à effectuer un couplage de la détection, de la synchronisation en temps et en fréquence, de l'analyse spatiale des champs reçus et de la séparation spatio-temporelle desdits signaux, de l'identification des signaux, de l'estimation des paramètres de direction d'arrivée et de localisation dans un processus conjoint, dénommé CC-EA-RD-DC, en se focalisant sur la présence et l'exploitation dans les signaux émis s(kTe) de séquences de références d(kTe) connues à priori. 5 L'invention concerne un procédé de réception, d'échantillonnage et de traitement des signaux échantillonnés à la période Te pour déterminer les paramètres de synchronisation, localisation et d'identification d'un émetteur s émettant une séquence de référence d connue de longueur K.Te ou les paramètres de synchronisation, localisation et d'identification de plusieurs émetteurs sj émettant des séquences de référence dj connues de longueur respectives Ki.Te,...,Kj.Te, au sein d'un système comprenant un ou plusieurs capteurs recevant le signal échantillonné x(kTe) provenant de l'émetteur s ou provenant des émetteurs multiples sj, caractérisé en ce qu'il comporte au moins les étapes suivantes : - détecter un premier signal si en utilisant la séquence de référence connue de Ki échantillons : di(kTe), - estimer le critère dl(/ ,m ,Ki) correspondant à la corrélation sur K1 échantillons du signal de référence di = di(kTe) inclus dans l'enveloppe complexe du signal émis s(kTe) de la A sortie d'un filtre wdi(I,m, K1), appliqué au signal reçu sur Ki échantillons et décalé en fréquence x((l+k)Te).exp{-2in.m.k/K1)} k=1,...,K1-1, où I est une hypothèse temporelle/retard-distance, m est une hypothèse fréquentielle/Doppler-cinématique, propres à la référence di, - rechercher les positions de synchronisation I (temporelle/retarddistance) et m1 (fréquentielle/Doppler-cinématique), propres à la séquence de référence di, qui maximisent le critère êdi (/,m,K1) , selon les hypothèses I et m, estimer la valeur du critère maximum, A êdl (11 ) et le filtre wdi(li,mi,Ki) correspondant, - estimer le canal hsi(11,m1) propre à la séquence de référence di aux positions de synchronisation (Il mi), - estimer l'incidence és d1(l1,mi) propre au signal de référence dl et aux positions de synchronisation (li, mi), et déterminer les paramètres de localisation de l'émetteur Ei émettant le signal s1 avec la séquence de référence di, - à partir de la valeur du canal hsi(Itmi) et de la sortie filtrée A Ysi(k-re)= wai(11,mi,K1)".x démoduler et décoder les messages de signalisation diffusés présents dans le signal si afin d'obtenir l'identification de l'émetteur E1, - à partir de la valeur du canal hsi(li,mi) et du signal démodulé si reconstituer le signal hsi(ltrni)*si, puis soustraire ce signal reconstitué au mélange de signaux x, pour constituer l'entrée x'=x- hsi(11,mi)*5i d'une deuxième étape du procédé, appliquée au signal X,, - réitérer ces étapes pour toutes les J séquences de référence d1, d2,..., dj ...dj détectées dans le mélange de signaux afin de déterminer successivement pour tous les émetteurs les paramètres de synchronisation (li,rni), j=1 .J, le filtre Wdi(li,mi,Ki) et le critère êdj(/,,m,,K,) associés, estimer le canal hsi(li,mi), puis démoduler le signal si et décoder les messages de signalisation diffusés par l'émetteur E, déterminer les paramètres d'incidences êàj,mj) j=1 .J et les paramètres de localisation associés à l'émetteur E, ainsi que les identifieurs associés à l'émetteur reconstituer le signal hsi(li,mi)*si propre à l'émetteur - reconstituer ensuite la somme des signaux hsi(11,ml)*sj +...+ hsj(Ij,mj)*sj associée à l'ensemble des émetteurs Ei...Ej afin de constituer l'entrée x(J)=bT=x-hsi(li,m1)*51-hs2(12,m2)*s2_ -hsJ(1J,n)*sJ d'une étape J+1 de traitement qui ne contient plus que la contribution des interférences résiduelles sans référence connue bT, - estimer la contribution des interférences résiduelles x(J)=bT sans référence connue, et déterminer la position des sources d'interférences résiduelles présentes dans bT, décoder les messages de signalisation diffusés présents dans b-r, et identifier les émetteurs des signaux interférents résiduels présents dans T. Selon une variante de réalisation, on détermine le critère és dl (11 Mi K1 ) aux positions de synchronisation (11, mi) en déterminant les valeurs des paramètres de l'hypothèse (li, mi) et les coefficients du filtre A Wd1(11111111K1), par comparaison d'une statistique êdi (11,mi,K1)=Ek-o...K1-1 (x((l1+k).Te)".Wd1(11,rni,K1).d(kTe)) à une valeur seuil f3, puis en recherchant les positions li m1 et la valeur précise du maximum de cette statistique au dessus du seuil p. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit d'un ou plusieurs exemples de réalisation annexés des figures qui représentent : - La figure 1, un schéma de dispositif pour la mise en oeuvre du procédé selon l'invention, - La figure 2, un schéma pour l'estimation d'un critère de synchronisation antibrouillé et du filtre spatial associé, particulièrement adapté aux cas de réception génériques à SINR faibles, - La figure 3, un synoptique représentant un exemple d'étapes mises en oeuvre par le procédé selon l'invention, - La figure 4, un schéma pour l'estimation d'un critère de synchronisation antibrouillé et du filtre spatial associé, dans le cas spécifique où les interférences présentent une propriété de non- circularité, et - La figure 5, un schéma pour l'estimation d'un critère de synchronisation antibrouillé et du filtre spatial associé, dans le cas spécifique d'un SINR fort et d'un vecteur spatial du signal utilise connu ou supposé.30 La description qui suit va être donnée pour illustrer le procédé selon l'invention pour un mélange de signaux en entrée d'une antenne de réception et d'un système de traitement. Le procédé selon l'invention consiste notamment à réaliser conjointement la détection de présence de 5 signaux émis s comprenant des séquences de références d connues a priori et présentes dans le mélange x de signaux en entrée d'antenne, l'estimation du filtre canal de propagation desdits signaux s et l'estimation des paramètres de localisation des sources desdits signaux. Il sera aussi possible de réaliser la détection et l'estimation des paramètres de 10 localisation des autres signaux présents dans le mélange en entrée, même s'ils ne présentent pas de signaux de référence. La figure 1 schématise un exemple de système permettant la mise en oeuvre du procédé selon l'invention. Plusieurs émetteurs E1, E3, En, sont reçus par un dispositif antennaire 1 comprenant un ou 15 plusieurs capteurs 6i de réception ou éléments antennaires. Le dispositif antennaire 1 est relié à un module 2 de traitement de signaux. Le module 2 comprend un processeur 3 adapté à exécuter les étapes du procédé selon l'invention qui vont être détaillées dans la suite de la description, ces étapes permettant notamment de localiser un émetteur, et de 20 l'identifier au moyen des séquences de références et du décodage de la signalisation qu'il diffuse par des procédés connus de l'Homme du métier. A un émetteur En est associé un signal d'émission s, et au signal une séquence de référence dn reconnaissable par le récepteur du système, ainsi que des messages diffusés au sein du signal sn qui sont décodables 25 et interprétables par le récepteur du système. La figure 3 schématise différentes étapes mises en oeuvre par le procédé selon l'invention. Avant d'expliciter cette figure, quelques éléments sur le traitement des signaux vont être exposés. La constitution et les récurrences du signal de référence d 30 présentes dans le signal émis s sont connues. Les algorithmes sont déroulés sur des signaux échantillonnés à la période Te, sur une certaine durée d'itération cohérente K.Te (i.e. pour des indices temporels k vérifiant 0 k K-1, et des instants vérifiant 0 tk=k.Te (K-1).Te sur une durée K.Te), pour des retards toTe (supposés, pour simplifier dans l'exemple non limitatif donné, multiples de la période d'échantillonnage), pour des décalages Doppler Afo = mo/(KTe). La résolution en fréquence/Doppler étant proportionnelle à 1/(KTe), on suppose pour simplifier les écritures, mais sans limiter la mise en oeuvre de l'invention, que les décalages Doppler sont des multiples de cette résolution. Sous forme échantillonnée, les signaux émis seront écrits s(kTe), les séquences de références d(kTe), le mélange de signaux en entrée d'antenne est noté x(k.Te), le signal en sortie de filtrage y(k.Te). Le vecteur filtre de canal de propagation us est lui aussi estimé à partir de la séquence de référence d(kTe) connu et aux récurrences connues, et ce pour tous les signaux reçus comportant un tel signal. Le vecteur hs peut le cas échéant être relatif à une réception multi-capteurs et correspondre alors : pour une propagation en espace libre, au vecteur directeur Ad(0d) de la source émettant d(kTe), vue sous la direction Od, pour une propagation à trajets multiples, au vecteur directeur d'un réflecteur, voire à un vecteur directeur moyen d'un cône de diffusion pour les cas plus complexes. La résolution en angle étant A® K/De, i.e. proportionnelle à l'inverse de l'ouverture De du réseau d'antennes utilisé, l'angle d'arrivée 0 est, par exemple, obtenu avec ou sans la connaissance du signal de référence d(kTe) présent dans le signal émis, à partir d'observation de versions décalées et translatées en fréquence des vecteurs x(kTe) des vecteurs signaux reçus sur le capteur. Ces observations comprennent pour certains des J signaux reçus dans le mélange les séquences de références d(kTe) (j=1,...,J), de constitution et de récurrence connues à priori, filtrées et bruitées durant la transmission. Pour d'autres signaux présents dans le mélange en entrée d'antenne ne comportant pas de références connues d(kTe), l'estimation angulaire est, par exemple, réalisée à l'aveugle sur le résidu x(kTe) du mélange en entrée x(kTe) après synchronisation, filtrage et soustraction des J signaux comportant les références connues di, j= 1,...,J. La figure 2 est un schéma décrivant un exemple de mise en oeuvre pour l'estimation ê d (11 ,Me ,K) du critère C(I',m') et pour l'estimation A A Wd(1' , rlY ,K) du filtre w(1',m'), en prenant en compte les interférences éventuelles, les paramètres de synchronisation étant réduits dans ce cas simplifié non limitatif aux paramètres temps et Doppler directement liés à la distance et à la vitesse de l'émetteur. Pour chaque émetteur E diffusant un signal faisant partie du mélange de signaux reçu sur le réseau d'antennes et qui comprend au moins un signal de référence d(kTe), propre à l'émetteur mais connu a priori, le procédé comporte par exemple les étapes suivantes (l'indice j est omis pour simplifier l'écriture, sans limitation du procédé), figure 2 : o on détermine au moins une statistique suffisante êd (i,m,K), 24, correspondant à la corrélation sur K échantillons du signal de référence d(kTe), 21, inclus dans l'enveloppe complexe du signal émis s(kTe) et de la sortie, y(I,m)= Wd(I,m)H.x(I.Te), 23, d'un filtre"Wd(I,m), 22, appliqué au signal total reçu x(kTe), i.e. êd(I,m,K) = y(I,m)H.d(kTe) = x(I.Te)H. (I, m).d(kTe) ; où I symbolise une hypothèse temporelle/retard-distance et où m symbolise une hypothèse fréquentielle/Doppler-cinématique, où H désigne la transposée conjuguée du vecteur signal, o pour cela, on détermine par exemple les valeurs des paramètres de l'hypothèse (I,m) et les coefficients du filtre"Wd(I,m) en comparant 24, 25, la valeur de la statistique êd(I,m,K) = x(I.Te)H.w(I,m).d(kTe) pour les paramètres de l'hypothèse (I, m) à une valeur seuil, notée 13 sur la figure 2 et sur les figures qui suivent, et en recherchant le maximum au dessus du seuil 13. Les valeurs lo et mo maximisant la statistique êd(lo,mo,K) seront dites positions de synchronisation temporelle et Doppler; o on détermine, aux différentes positions de synchronisation (10, mo) obtenues, au moins une estimation de canal Tic/ (10,mo) correspondant à la réponse du canal de propagation relativement au signal de référence d(kTe) intégré sur la durée K.Te; 0 on détermine, aux différentes positions de synchronisation (/0,mo) obtenues, au moins une estimation angulaire êd (IO5M0), 25, de la direction correspondant au signal de référence d(k.Te), 21, intégré sur la durée K.Te, d'après la composante vectorielle :4d d(Oc(io, mo)), 25, de l'estimation de canal hd(io, mo), 25, à la position de synchronisation (bo, mo) et un modèle de diagramme de rayonnement de l'antenne, ou une table de calibration de l'antenne La figure 3 représente un exemple de déroulement des étapes du procédé selon l'invention pour un mélange de J signaux 5i ...5j comportant des références connues di ...dJ, signaux auxquels s'additionne un bruit total b-r ne comportant pas de références connues, 20 en J+1 étapes successives : o étape 1: application du critère conjoint CC-EA-RD-DC à la référence di correspondant au signal si, dans le mélange initial x en entrée : 25 A - 311, estimation du critère êdi(ii,m1,1C1) du filtre w cm(11,mi,K1) associé au critère et des positions de synchronisation li m1 propres à la référence d1 connue sur Ki échantillons, - 312, estimation du canal hsi(11,mi) propre à di aux instants li . 313, estimation de l'incidence (5d1(l1,mi) propre à di aux instants li mi d'après la composante vectorielle du canal hs1(11,m1) aux instants li, ml, pour localiser 116 l'émetteur, A . 314, exploitation de la sortie filtrée ysi(kl-e)= w di(11,mi,Ki)".x et des estimations de canal hsi(li,mi) propres à di aux instants li mi pour égaliser démoduler et décoder le signal si afin d'obtenir 315 l'identification de l'émetteur correspondant, o en fin d'étape 1, reconstitution 317 du signal hsi(li,mi)*si, puis soustraction 318 au mélange en entrée x, pour constituer l'entrée x'=x12s1(11,mi)*5i de l'étape 2 de traitement, o étape 2: application du critère conjoint CC-EA-RD-DC à la référence d2 correspondant au signal s2, dans le mélange initial de signaux x' en entrée: . 321, estimation du critère êd2(/2,m2,K2) du filtre A associé wd2(I2,m2) associé au critère et des positions de synchronisation 12 m2 propres à la référence d2, . 322, estimation du canal hs2(I2,m2) propre à d2 aux instants 12 M2, . 323, estimation de l'incidence êd2(12,m2) propre à d2 aux instants 12 21 d'après la composante vectorielle As2(ê d202, M2 )) du canal hs2(I2,m2) aux instants 12 m2, cette valeur sera utilisée pour déterminer 326 la localisation de l'émetteur E2, A . 324, exploitation de la sortie filtrée ys2(kTe)= w d2(12,m2,K2)".x et des estimations de canal hs2(I2,m2) propres à d2 aux instants 12 m2 pour égaliser démoduler et décoder le signal s2 afin de parfaire l'identification, 325, de l'émetteur correspondant, o en fin d'étape 2, on va reconstituer, 326, le signal hs2(12,m2)*s2, puis on soustrait 327 au mélange de signaux en entrée x, pour constituer l'entrée x"=x-hsi(li,mi)*sihs2(12,m2)*s2 de l'étape suivante de traitement ... et ainsi de suite jusqu'à épuisement des références d1, d2,..., dj o étape J, application du critère conjoint CC-EA-RD-DC à la référence dj correspondant au signal sj puis à la fin de l'étape traitant la référence dj, reconstitution 330 de la somme des signaux hsi(11,m1)*5J +...+ hsJ(Ilmj)*sj associée à l'ensemble des émetteurs El...Ej afin de constituer l'entrée 331 x(J) = bT =x - hsi(11,mi)*51 - hs2(12,m2)*s2 - hsJ(1J,mj)*5j d'une étape finale J+1 de traitement, ne comprenant plus que le bruit total résiduels bT (= bruit + brouilleurs résiduels) sans références connues, o étape J+1 : on va appliquer des techniques alternatives d'estimation de la contribution des interférences résiduelles bT sans référence connue, 332: - Goniométrie des sources d'interférences résiduelles sans référence connue, présentes dans bT, 3001, localisation des émetteurs sans référence connue, 3003, - Identification, démodulation et décodage de la signalisation des signaux interférents résiduels sans référence connue, 3002, présents dans bT, 3004, identification des émetteurs, localisation 3005 des émetteurs. Une autre variante d'application de réalisation de l'invention, optimale lorsque les signaux utiles et les interférences présentent une propriété de non-circularité, est présentée figure 4. Le cas le plus favorable mais non limitatif correspond au cas où le signal utile émis et les interférences sont rectilignes, c'est-à-dire que leurs enveloppes complexes si(kTe) s2(kTe) sj(kTe) sont à valeurs réelles. C'est le cas en particulier pour des signaux émis modulés en amplitude (AM, ASK) ou 5 en phase à 2 états (BPSK). Dans ces conditions, il est intéressant de substituer à la statistique précédente 24 figure 2, êd (/',m',K), la statistique êd(/',m',K), 43 figure 4) , définie par: à ( Wrd ) d,m(/'Te)H fi 5e, In( /1Te)1. 3-éd,m(I'Te) H T OÙ .i'Akerre) [Lcm'((k+11)Te) , xm'((k+1')Te) ] , 10 et où Itcl 7 .icl,m'(1Te) et R im(lTe) sont définis par : K-1 d = 142 = E d(kt)12 k'=0 xd,m' (1') = L1 (r-Fle)d* (k').e K-1 2ur k'm' k'=0 K-1 Rxe,(1') (141e) k'=-0 où T désigne la transpositon et H désigne la transposition conjugaison. où m((k+I)Te) remplace x ((1+1')Te) dans les statistiques correspondant à la figure 2 précédente. Les étapes exécutées pour la 15 mise en oeuvre de ce filtre face à un signal utile ou interférent sont identiques à celles décrites à la figure 2. Les étapes exécutées pour la mise en oeuvre successive de ce filtre face à un mélange de signaux utiles et interférents sont identiques à celles décrites à la figure 3. Cette statistique s'applique aussi pour des signaux émis 20 devenant quasi-rectilignes après un prétraitement de dérotation. C'est en particulier le cas pour des signaux émis de typen/2-BPSK ou OQPSK ou encore certaines modulations de phase continue ou CPM (Continuous Phase Modulation) c'est-à-dire modulés en fréquence et à phase continue, comme par exemple les signaux MSK, GMSK d'indice 1/2 ' fréquemment utilisés dans les normes radiocellulaires et les transmissions militaires. Dans ces conditions, la statistique précédente doit être mise en oeuvre mais en substituant aux vecteurs spatiaux xm((k+I)Te, les vecteurs spatio-temporels dérotés à L prises temporelles :dm,st((k+i)Te) Udd(k+(L1)/2)Te)T,..., 3-c'dm((k -(L1)/2)Te)Tir si L est impair et3-é dm,st((k-Fi)Te) 4 dm((k.. dm((k- L/2)Te)rfilr si L est pair, où :i'dm((k+i)Te) [Lklm((k-Fi)Te)T xdm((k+1)Te)tT , xdin((k+I)Te) xcle-Fige) exri[j27cm(k+1)/K], xcl(0 4 t/Ts .x(t) et Ts est la durée symbole du signal émis. La figure 5 présente une mise en oeuvre particulière de l'invention lorsque le vecteur directeur As est connu ou supposé, les paramètres cinématiques étant réduits dans ce cas simplifié non limitatif au paramètres retards et Doppler directement liés à la distance et à vitesse relative de l'émetteur ou du réflecteur. Lorsque la propagation s'effectue en espace libre, le vecteur canal utile hs s'écrit sous forme condensée hs à hs eid's As où hs (1s et As correspondent respectivement à l'amplitude, à la phase et au vecteur directeur normalisé (première coordonnée réelle positive) du signal émanant de l'émetteur ou du réflecteur. Lorsqu'un balayage spatial de l'espace est envisagé, comme cela est le cas pour les structures conventionnelles de détection, on peut considérer que le vecteur As est connu lorsque la case spatiale de l'émetteur ou du réflecteur est scrutée. Dans les conditions, les structures conventionnelles de détection ne sont optimales que pour des signaux émanant d'émetteurs ou de réflecteurs de faible puissance, et peuvent devenir sous-optimales pour des signaux émanant d'émetteurs ou de réflecteurs de forte puissance. Lorsque le vecteur directeur As est connu ou supposé, une statistique suffisante pour la détection du signal connu s(kTe) sur la durée 0 k K-1 dans la case spatiale As, la case distance (ou encore hypothèse temporelle/retard-distance de l'objet) et la case Doppler m (ou hypothèse fréquentielle de l'objet mobile) est donnée par, référencée 51, sur la figure 5, AH]-1(Plaîxd,'',(PTe) F(P,m' \ ,As,K) = (1 I "cc/s) (AsH ^ V R l' Te ),1 A1 - m', K)) où Cd(P,m',K), 51, est défini en figure 2. L'expression précédente est une extension de Cd(/',m',K) lorsque le vecteur As est connu, m', l' étant respectivement la case Doppler et la case distance. Pour des SINRy faibles (SINRy Traitements supplémentaires Des traitements complémentaires peuvent être envisagés aux étapes décrites préalablement, par exemple : - des techniques de goniométrie à simple résolution, à commutation d'antenne ou parallèle haute résolution, visant à localiser angulairement 25 un émetteur correspondant à une position de synchronisation donnée peuvent être envisagées à partir des observations xm((k+I)Te), m((k+1)Te) ou dm,sd(k+0 Te) associées à cette synchronisation, - des techniques de goniométrie à haute résolution visant à localiser angulairement des émetteurs multiples correspondant à une position de 30 synchronisation donnée peuvent être envisagées soit à partir des observations xm((k-FOTe), m((k+1) Te) ou dm,st((k+4 Te) associées à cette synchronisation, soit à partir des résultats de corrélation des observations avec le signal de référence émis d(kTe) sur des intervalles sur multiples ou sous-multiples de K, c'est à dire à partir des vecteurs Pxdan(LK'Te), K' divisant K ou K divisant K', c'est le cas de la goniométrie à référence. Dans le cas où A, est connu, ces techniques peuvent être mises en oeuvre après un filtrage spatial des observations limitées à la case spatiale A. Exemples illustratifs et non limitatifs de mise en oeuvre : Un exemple illustratif et non limitatif de statistiques, d'estimées de canal et d'estimées angulaires relatives au signal de référence d(k.Te) est donné ci-après.: Posons les grandeurs suivantes : Le réseau antennaire comportant N éléments, chacune des grandeurs x, h, i, n, bT, est une vecteur colonne de taille N Dans un modèle spéculaire classique h, s'écrit sous la forme M-1 abrégée h, (1,m) = E K .6(1 - - l'=0 m'=0 I' est le retard du trajet, m' est le doppler du trajet, h,(1',m') est l'amplitude du trajet, (1),(1',m') est la phase du trajet, A,(0,(1',m')) est le vecteur directeur du trajet associé à la direction 0,(1',m'), 6 est la distribution de Dirac ; En contexte mono-source et dans un modèle spéculaire classique, x s'écrit sous la forme : M-1 L-1 2uc-m' x(1.Te)= E .e K .S((1 - 1')Te)+12-r(1.Te) r=0 'z,=0 Pour les statistiques à l'ordre 2 des signaux émis et reçus estimées sur la durée d'intégration cohérente K, on obtient : d =1 Ëld(k')12 k' =0 K-1 xd,m (1) =Ix(1 + k')d* (k? ).e k'=0 K-1 x (1) =EX(1k')''+ k') k'=0 Un critère de synchronisation peut s'écrire de la façon suivante : êd (1,m,K) =1 îxlid,m (1)â (1)il xd,rn (1) r d Ce critère est associé au filtre 11' d (1 ,n K) = 1 (1) xd (1) . d s'appliquant au vecteur signal en entrée x((l+k)Te) pour produire le signal en sortie i)d,', = d H K)-x d ,m(/ e) Pour des signaux contenant la séquence de référence d(k.Te), une estimation du canal de propagation peut s'obtenir, pour chaque trajet de retard l', de Doppler m', à partir de la séquence de référence d(k.Te), de la façon suivante : h,/ , = xHa , (1' = d (1' ,111')-ee .51,1 Ceci conduit à la réponse impulsionelle complète déterminée d'après la référence d(k.Te) : M-1 e s d (11 'ln') .8 - ) L-1 -m' f I d (1,M) = hd (I' , m' ).e K l'=0 m'=0 km 7T Le cas d'application précédemment détaillé de l'invention est particulièrement indiqué pour l'estimation de la réponse impulsionelle du canal de propagation des signaux en bande large et présentant des signaux de référence d(kTe) 0 k K-1 de durée K significative, procurant de ce fait une haute sensibilité, de bonnes résolutions en temps et en Doppler, et une forte capacité d'antibrouillage et de réjection des interférences. D'autres techniques de détection et de synchronisation connues de l'homme de l'art telles que le radiomètre, la détection cyclique des fréquences porteuse et des rythmes symboles etc. peuvent de même s'appliquer aux signaux présents dans le mélange en entrée d'antenne après filtrages et soustractions successives préalables des signaux présentant la référence d(k.Te), une fois estimées les réponses impulsionnelles /d associées aux dites références d(k.Te), c'est à dire s'appliquer au signal xs(I.Te) obtenu à partir du signal d'entrée x(k.Te) du filtre 1d et de l'estimée de canal hd relative à d par : 51 d (1.Te) = id (1, m,K).x(1.Te) .x' (1.Te) = x(1.Te)-(h d * 5d)(1.Te) 111-1 = x(1.Te) - Eh d 2ur- K .5 1 d((1 - 1').Te)) 1 .= 0 m=0 Ces techniques sont notamment intéressantes pour les applications suivantes : - lorsque tous les signaux du mélange ne comprenant pas de référence connue d(kTe) ont été filtrés et soustraits, - lorsque le mélange x de signaux en entrée présente un signal s dominant sans référence connue, - lorsqu'une des méthodes connues de l'homme de l'art s'avère en pratique plus performante dans la situation particulière de réception du mélange x(k.Te). Pour des signaux contenant la séquence de référence d(k.Te), une estimation du vecteur directeur de la source ou du trajet de la direction d'arrivée correspondant au retard I et au Doppler m, peut s'obtenir, à partir de la séquence de référence d(k.Te), à partir du vecteur directeur iîd(1,M) associé à l'estimation de h d(1,M), que l'on compare à des vecteurs calibrés tabulés pour en déduire l'incidence d(/,m) par un critère du type suivant (exemple illustratif et non limitatif) : arg max o' {Ad (1, .A Co' )} vecteurs normalisés tabulés en nombre J et de la forme /4.10=[5 AJ,2(") ri réel positif 4:(0' ) eVect{211(0;),...,A j(0 Cette façon de procéder pour l'estimation des directions d'arrivées des signaux en bande large et présentant des signaux de référence d(kTe) 0 k K-1 de durée K significative, procure une haute sensibilité, de bonnes résolutions en temps et en Doppler, et une forte capacité d'antibrouillage et de réjection des interférences. D'autres techniques d'estimation des signatures spatiales connues de l'homme de l'art telles que les méthodes de corrélation vectorielle en situation de réception mono source ou à source dominante, les méthodes à super résolution de type Capon connues de l'homme du métier, les méthodes à haute résolution de type Music, Jade, connues elles aussi de l'homme du métier, etc.... peuvent de même s'appliquer aux signaux présents dans le mélange en entrée d'antenne après filtrages et soustractions successives préalables des signaux présentant la référence d(k.Te), une fois estimée la réponse impulsionnelle d associée à ladite référence d(k.Te), c'est à dire s'appliquer au signal x'(I.Te) obtenu à partir du signal d'entrée x(I.te) du filtre 1;y'd (/,m,K)et de l'estimée de canal TI d relative à d par: d(1.Te)= ?d (1,m, K).x(1.Te) x'(1.Te) = x(1.Te) -(d * Yd)(1.Te) L-1 ir' = x(1.Te)- E 11 d(1' 2uK d((1 -1').Te)) r=o m=0 Ces techniques alternatives sont notamment intéressantes pour les applications suivantes : - lorsque tous les signaux du mélange ne comprenant pas de référence connue d(kTe) ont été filtrés et soustraits, - lorsque le mélange x de signaux en entrée présente un signal s dominant sans référence connue, - lorsqu'une des méthodes connues de l'homme de l'art s'avère en pratique plus performante dans la situation particulière de réception du mélange x(k.Te). La détermination de la statistique êd(1,m,K) fournit notamment : - la détection de l'émetteur diffusant le signal de référence d(kTe), ou l'identification de la nature de l'émetteur via la reconnaissance de son signal de référence d(kTe) (qui maximise le critère ê d(1,m,K)), ou - la synchronisation retard Doppler de l'émetteur (paramètres lo et mo qui maximisent le critère ê d(1,M,K)), ce qui permet d'en déduire, à partir d'un réseau de (plusieurs) senseurs en réception, - les retards absolus et/ou les retards différentiels correspondant à l'émetteur, ou - les estimées de localisation TOA (time of arrivai) et/ou TDOA (time difference of arrivai) à partir des estimées précédentes selon des procédés bien connus de l'homme de l'art, ou - les Doppler absolus et/ou les Doppler différentiels correspondant à l'émetteur, ou - les estimées de localisation FOA (frequency arrivai) et/ou FDOA (frequency differential of arrivai) à partir des estimées précédentes selon des procédés bien connus de l'homme de l'art. Le filtre i;svd (/, m,K) associé au critère êd (/,m,K) procure un signal en sortie y(k.Te)= Wd (/,m,K) .x(k.Te) purifié des interférences et des brouillages dus aux sources contenant la référence d(k.Te) présents dans le mélange en entrée x(I.Te). Outre la reconnaissance procurée par la détection du signal de référence d, le signal y peut alors être démodulé et décodé de manière à ce que le contenu décodé des messages identifie complètement et sans ambiguïté l'émetteur. La détermination du filtre de canal de propagation hd associé à d fournit notamment l'ensemble des caractéristiques de la propagation du signal incluant la référence d : - les niveaux des trajets principaux et réfléchis d'une part, qui permettent d'appliquer des méthodes de tracking en puissance (homing) pour favoriser la localisation des émetteurs depuis une plate-forme mobile, et/ou - les distributions statistiques d'amplitude, de retard et de Doppler correspondant à l'émetteur du signal de référence d, ce qui fournit des indications précieuses sur l'environnement de propagation et sur sa stationarité, aidant l'ingénierie de réseau d'une part et la post-intégration et la fusion des estimateurs de localisation d'autre part, et/ou des filtres à appliquer pour éliminer toutes les contributions de la source émettant ce signal dans le mélange x(I.Te) reçu sur antenne, afin d'itérer le processus aux signaux présentant d'autres séquences de références connues puis, le cas échéant, aux signaux ne présentant pas de référence connue. La détermination des incidences êd(/',m') de la source et des éventuels réflecteurs associés pour différents paramètres retard l' et Doppler m' fournit les paramètres d'entrée d'estimateurs de localisation basés par exemple : - sur des techniques de triangulation à partir d'estimations d'angles d'arrivées (AOA), ou sur des techniques de localisation exploitant le défilement angulaire et basées sur des estimations successives d'angles d'arrivées (VOAE), ou sur des techniques mixtes fusionnant à partir d'une plate forme mobile des estimations angulaires (AOA) et des estimations de puissance reçue entretenues au cours du temps (dites aussi AOA en mode « rapprochement par scrutation de la puissance reçue connu sous l'expression « mode homing » par les hommes du métier), ou sur des techniques mixtes fusionnant estimations angulaires, estimation de synchronisation temporelle (TOA ou TDOA), estimation de synchronisation en fréquence (FOA ou FDOA), estimations de puissances reçues, instantanées ou entretenues au cours du temps, etc. Exemples d'applications Le procédé selon l'invention est, par exemple, utilisé avec un système de réception de position connue, et face à des émetteurs ou des réflecteurs sur lesquels des informations de synchronisation, de position ou de cinématique sont disponibles sur l'émission ou la réflexion des signaux exploités par l'invention, ce qui permet de déduire de l'invention des informations temporelles (retard/distance) absolues et des information fréquentielles (Doppler vitesse/absolue) absolues à partir des informations relatives produites par le procédé objet de l'invention. Selon un exemple de mise en oeuvre, l'utilisation du procédé est dupliquée sur plusieurs systèmes de réception synchronisés et communicants entre eux ou avec un centre d'exploitation, pour produire des informations de synchronisation, de position ou de cinématique sur les émetteurs et réflecteurs afin d'obtenir des informations temporelles et fréquentielle absolues sur les émetteurs et les réflecteurs. Le procédé selon l'invention peut notamment être appliqué à la détection, l'identification et la goniométrie conjointes de signaux interférant des émissions de type GSM, UMTS, DVB, WiFi, WiMax, LTE, LTE-A, et bien d'autres. Le procédé selon l'invention est utilisé dans tout système dans lesquels il est possible d'utiliser une connaissance a priori des signaux. Des exemples de système sont donnés ci-après à titre illustratif et nullement limitatif. Les émetteurs de radiocommunications numériques ou acoustiques incluent dans leurs signaux des séquences particulières (dites séquences de référence) ainsi que des voies de signalisation, destinées à permettre la synchronisation des récepteurs, à favoriser la démodulation le décodage et l'interprétation des messages transmis. Dans les réseaux standardisés de radio transmission (GSM, 3GPP, 3GGP2, Wifi 802.11xx, Wimax 802.11yy, LTE, LTE-A, DVB-T/H, etc.), ces séquences et voies de signalisation sont précisément décrites. Elles constituent une information très complète dont tire profit l'invention en terme de synergie entres les traitements conjoints susmentionnés, s'appliquant tant aux signaux utiles qu'aux sources d'interférences éventuelles. Tout opérateur ou constructeur de réseau de communication ou de réseau de radiolocalisation connaît par principe les séquences de référence et voies de signalisation qui constituent la base du fonctionnement de son système, et se trouve de fait en mesure d'appliquer la présente invention. Tout opérateur ou constructeur de radar, sonar ou autre système actif (sondeur de canal, pseudolite en radio-localisation, etc.), maîtrise de même les signaux générés par son effecteur ou par son réseau d'effecteurs, et se trouve de même en mesure d'appliquer la présente invention. La présente invention s'applique à la métrologie des émetteurs 30 radioélectriques (voire acoustiques), au contrôle du spectre radioélectrique dans les bandes de fréquence utilisées par les radiocommnications et par les radars. Sans sortir du cadre de l'invention, le procédé selon l'invention peut aussi être utilisé dans de nombreux domaines de l'électromagnétisme et de l'acoustique, en détection, en transmission, en localisation et navigation, afin de mesurer précisément les signaux radioélectriques et acoustiques reçus, d'améliorer la connaissance de l'environnement et du milieu de propagation, d'améliorer en conséquence les traitements des signaux utiles, la mesure de l'occupation des fréquences dans les bandes des radiocommunications, dans les bandes radar, voire dans le domaine acoustique, de conduire un contrôle performant de l'utilisation du spectre radio électrique, d'améliorer le déploiement, la planification et l'ingénierie des réseaux de radiocommunication et de radiodiffusion. L'invention s'applique de même pour de nombreuses applications de surveillance et de contrôle : le contrôle aérien, le contrôle littoral, le contrôle des frontières, etc. Avantageusement, le procédé selon l'invention peut être exécuté pour un mélange de signaux avec dans ce mélange des signaux dont la puissance est inférieure au signal dominant, l'écart de puissance, traduisant à la fois un gain sur le bruit et un gain de rejection des signaux interférants, pouvant atteindre - 10 à 20 dB selon les cas dans le cas de signaux TDMA/GSM, - 20 à 25 dB dans le cas de signaux OFDM/DVB-T/H OFDM/VVifi OFDM/VVimax, OFDM/LTE, 20 à 30 dB dans le cas de signaux CDMA/3GPP et CDMA/2GPP2, voire plus dans certains cas. Entre autres avantages, le procédé selon l'invention permet notamment d'accéder directement aux positions des émetteurs et des réflecteurs sans avoir à balayer tout l'espace et donc sans avoir à scruter les directions dépourvues d'émetteur et de réflecteurs. Il permet aussi de relâcher les contraintes sur la maîtrise du diagramme du réseau de capteurs (par exemple, en cas de couplage entre éléments rayonnants, la calibration n'est plus nécessaire tant que l'estimation angulaire n'est pas réalisée), et il s'accommode parfaitement d'une propagation avec trajets multiples. Enfin, l'utilisation du critère couplé permet en pratique un réglage de seuil indépendant du niveau de bruit présent en entrée mais également une protection intrinsèque vis-à-vis des interférences sans risquer de rejeter les signaux utiles à la mesure et sans avoir à insérer de traitements pour augmenter la robustesse.10 REVENDICATIONS1 - Procédé de réception, d'échantillonnage et de traitement des signaux échantillonnés à la période Te pour déterminer les paramètres de synchronisation, localisation et d'identification d'un émetteur s émettant une séquence de référence d connue de longueur K.Te ou les paramètres de synchronisation, localisation et d'identification de plusieurs émetteurs sj émettant des séquences de référence dj connues de longueur respectives Ki.Te,...,Kj.Te, au sein d'un système comprenant un ou plusieurs capteurs (6i) recevant le signal échantillonné x(kTe) provenant de l'émetteur s ou provenant des émetteurs multiples caractérisé en ce qu'il comporte au moins les étapes suivantes : - détecter un premier signal si en utilisant la séquence de référence connue de Ki échantillons : di(kTe), - estimer, (311), le critère di(/ ,m ,K1)correspondant à la corrélation sur K1 échantillons du signal de référence di = di(kTe) k=0,...,K1-1 inclus dans l'enveloppe complexe du signal émis A s(kTe) de la sortie d'un filtre e(1,m, K1), appliqué au signal reçu sur K1 échantillons et décalé en fréquence x((l+k)Te).exp{- 21n.m.k/K1)} où lest une hypothèse temporelle/retard- distance, m est une hypothèse fréquentielle/Doppler-cinématique, propres à la référence di, - rechercher les positions de synchronisation li (temporelle/retarddistance) et m1 (fréquentielle/Doppler-cinématique), propres à la séquence de référence di, qui maximisent le critère Cd1(/,m,K1), selon les hypothèses I et m, estimer la valeur du critère maximum, A e'd1 (11 K1 ) et le filtre wci1(11,m1,K1) correspondant, - estimer, (312), le canal h51(11,mi) propre à la séquence de référence di aux positions de synchronisation (li mi),- estimer, (313), l'incidence (5d,(l1,mi) propre au signal de référence di et aux positions de synchronisation (11, mi), et déterminer les paramètres de localisation de l'émetteur Ei émettant le signal si avec la séquence de référence di, - à partir de la valeur du canal hsi(li,mi) et de la sortie filtrée Ysi (k-re)= A Wd1(117M17K1)H.X (314), démoduler et décoder les messages de signalisation diffusés présents dans le signal si afin d'obtenir l'identification de l'émetteur Ei, - à partir de la valeur du canal hsi(li,mi) et du signal démodulé si reconstituer le signal hsi(11,mi)*si, (317), puis soustraire (318) ce signal reconstitué au mélange de signaux x, pour constituer l'entrée x'=x-hs1(11,mi)*si d'une deuxième étape du procédé, appliquée au signal x', - réitérer ces étapes pour toutes les J séquences de référence di, d2,..., d .dj détectées dans le mélange de signaux afin de déterminer successivement pour tous les émetteurs E1 E, les paramètres de synchronisation (11,m1), j=1 .J, le filtre et le critère êdj(ii,mi,KJ) associés, estimer le canal hsi(li,mi), puis démoduler le signal si et décoder les messages de signalisation diffusés par l'émetteur E, déterminer les paramètres d'incidences êàj,mj) j=1...J et les paramètres de localisation associés à l'émetteur E, ainsi que les identifieurs associés à l'émetteur reconstituer le signal hsi(li,mi)*si propre à l'émetteur - reconstituer (330) ensuite la somme des signaux hsi(11,mi)*sj +...+ hsj(Ij,mj)*sj associée à l'ensemble des émetteurs Ei...Ej afin de constituer l'entrée x(J)=bT=x-hsi(11,m1)*si-hs2(12,m2)*s2_ -hsJ(1,,,ms)*ss, d'une étape J+1 de traitement qui ne contient plus que la contribution des interférences résiduelles sans référence connue bT,- estimer (331) la contribution des interférences résiduelles x(J)=bT sans référence connue, et déterminer la position des sources d'interférences résiduelles présentes dans IDT, décoder les messages de signalisation diffusés présents dans bT, et identifier les émetteurs des signaux interférents résiduels présents dans T. 2 - Procédé selon la revendication 1 caractérisé en ce qu'on détermine le critère e'd1 (11 K1 ) aux positions de synchronisation (li, mi) en déterminant les valeurs des paramètres de l'hypothèse (I1, mi) et les A coefficients du filtre wdi(l1,mi,K1), par comparaison (24) d'une statistique (x((li+k).Te)".1ki(li,mi,K1).d(kTe)) à une valeur seuil f3 (25), puis en recherchant les positions li m1 et la valeur précise du maximum de cette statistique au dessus du seuil 13. 3 - Procédé selon la revendication 1 caractérisé en ce que l'on réalise une étape d'estimation de canal hsi(11,m1) aux positions (Ii, mi), à partir des estimateurs statistiques ayant permis de construire les estimées du A filtre wdi(I,m,Ki) et du critère êdi(/,m,Ki) aux positions (11, mi). 4 - Procédé selon la revendication 1 caractérisé en ce que l'on réalise une étape d'estimation de goniométrie êdi(l1,m1) aux positions li et m1 à partir de la valeur estimée du canal hsi(11,m1) aux positions (11,m1), à partir A des estimées du filtre wd1(11,m ,K1) et du critère êdi(ii,mi,Ki) aux positions (11,mi). 5 - Procédé selon la revendication 1 caractérisé en ce que l'on réalise une étape d'estimation de canal hs1(11,mi) aux positions (lie mi), complétée par une estimation de canal hsi(I,m) pour les paramètres (I, m) variables au voisinage des positions (11, mi), à partir des estimateursA statistiques ayant permis de construire les estimées du filtre wci1(1,m,K1) et du critère êdi(1,m,KI) aux positions (li, mi) et aux positions (I, m) voisines de (li, mi). 6 - Procédé selon la revendication 3 caractérisé en ce que l'on utilise A directement les estimées du filtre Wd1 (l1 e ) ou critère êdl (119 Mi ,K1) pour estimer le canal hsi(li,mi) à chaque position (11,mi) réalisant le maximum dudit critère. 7 - Procédé selon la revendication 1 caractérisé en ce que les capteurs sont un réseau d'antennes, et pour un mélange des signaux reçus correspondant à un signal large bande pour le réseau d'antennes, on structure la bande de réception du signal en plusieurs sous-bandes étroites pour le réseau et on exécute les étapes de la revendication 1 sous-bande par sous-bande. 8 - Procédé selon la revendication 1 caractérisé en ce que le mélange des signaux reçus est à large bande devant la porteuse et en ce que l'on exprime le signal s correspondant à la transformation du mélange de signal s par les variables retard-position et Doppler-cinématique, compression retard/position et compression du Doppler-cinématique et on utilise ce signal transformé dépend de 4 paramètres au lieu de 2 pour exécuter les étapes de la revendication 1. 9 - Procédé selon la revendication 1 caractérisé en ce que lorsque la propagation s'effectue en espace libre, on exprime la valeur du canal utile ils sous forme condensée lis A hs eM's As où hs et As correspondent respectivement à l'amplitude, à la phase et au vecteur directeur normalisé du signal émanant de l'émetteur.30 10 - Utilisation du procédé selon l'une des revendications précédentes pour la détection, l'identification et la goniométrie conjointes de signaux interférant des émissions de type GSM, UMTS, IS 54 IS 136, IS 95/A/B LTE, LTE-A, 3GPP, 3GPP2, DVB-T/H, WiFi 802.11xx, SRC 802.15yy, Wimax 802.16zz 11 - Utilisation du procédé selon l'une des revendications précédentes pour déterminer des paramètres de localisation et d'identification d'un émetteur dans un système de radiocommunications numériques ou acoustique incluant dans un signal des séquences de référence ainsi que des voies de signalisation adaptées à permettre la synchronisation des récepteurs, à favoriser la démodulation, le décodage et l'interprétation des messages transmis. 12 - Utilisation du procédé selon l'une des revendications précédentes à la métrologie d'émetteurs radio-électriques dans des bandes de fréquences du domaine des radiocommunications.