FR 2738090 A1 19970228 FR 9509976 A 19950822 La présente invention concerne un procédé et un dispositif de multiplexageldémultiplexage spatial de signaux radioélectriques organisés en trames, de caractère non gaussien, et transmis en duplex sur un même canal fréquentiel et temporel ne nécessitant pas de calibration, et s'applique plus particulièrement aux radiocommunications cellulaires entre au moins une station de base et plusieurs mobiles. Avec l'accroissement constant de la demande de communications mobiles et la limitation du nombre de canaux alloués aux opérateurs des réseaux de radiocommunications cellulaires, les problèmes de saturation vont devenir cruciaux dans les années à venir. Des techniques classiques permettent actuellement de multiplexer les communications en fréquence (AMRF), en temps (AMRT), ou en code (AMRC). L'AMRF (Accès Multiple à Répartition de Fréquence) consiste à allouer à chaque communication une paire de fréquences (une pour la voie descendante et une pour la voie montante). Toute seule, cette technique très simple de multiplexage n'offre quasiment pas de perspective d'amélioration de l'efficacité spectrale. L'AMRT (Accès Multiple à Répartition dans le Temps) consiste à allouer à chaque communication un intervalle temporel précis. Seule, cette technique, déjà optimisée pour amélioration de l'efficacité spectrale dans le cadre de la transmission de phonie par exploitation des silences présents sur le signal de parole, laisse entrevoir peu d'amélioration supplémentaire sur ce plan. L'AMRC (Accès Multiple à Répartition dans les Codes) consiste à allouer à chaque communication un code, définissant une loi de saut en fréquence sur des intervalles de temps courts. Même si elle permet une certaine augmentation de l'efficacité spectrale (saturation progressive par dégradation du signal et non brutale par interruption du service comme dans le cas des deux techniques précédentes), cette méthode reste lourde et chère à mettre en oeuvre. Plusieurs voies sont actuellement explorées pour améliorer les capacités des systèmes de communication cellulaires à savoir la diversification des cellules et des formes d'ondes avec la mise en oeuvre de cellules de tailles de plus en plus réduites ("microcell" et "nanocell" du projet de normalisation ATDMA, abréviation anglosaxonne pour "Advanced Time Division Multiple Access", pour les zones urbaines), et surtout une nouvelle méthode de multiplexage pouvant se combiner aux techniques existantes l'Accès Multiple à Répartition Spatial (AMRS ou SDMA en terminologie anglo-saxonne). Cette technique est basée sur l'idée que les signaux radios correspondant à plusieurs liaisons entre une station de base et des mobiles peuvent être dans de nombreux cas décorrélés spatialement. L'exploitation de cette décorrélation permet de séparer les canaux spatiaux (même s'ils occupent tous un même canal fréquentiel et temporel) au moyen d'un système émission/réception multivoie associé à une base d'antennes adaptée à la gamme de fréquences utilisée. Plusieurs méthodes mettant en oeuvre ce concept ont déjà été proposées. Une première méthode est basée sur des techniques d'estimation de direction d'arrivée comme le système décrit dans le brevet de Richard ROY (Brevet US PCTIUS92II 0074, 12/1991, intitulé "Spatial Division Multiple Access Wireless Communication System"). Une deuxième méthode utilise des caractéristiques de cyclostationnarité des signaux de transmissions numériques pour les séparer par des traitements en aveugle (Brevet US N" 5260968, du 9/11/93, intitulé "Method and Apparatus for Multiplexing Communications Signals through Blind Adaptative Spatial Filtering" de William A. GARDNER et Stephan V. SCHELL). Pour la première méthode, l'estimation des directions d'arrivées de différentes sources radioélectriques à partir des signaux reçus sur une base d'antennes, nécessite la parfaite maîtrise du front d'onde ce qui entraîne une contrainte de calibration de la base d'antennes utilisée. Outre le surcoût induit par la calibration, l'écart type de l'estimation indirecte du vecteur directeur à partir de la direction d'arrivée estimée et de la table de calibration de la base d'antennes est nettement plus important que celui obtenu par des méthodes d'estimation directe. Le vecteur directeur correspond à la valeur prise par la fonction de transfert de la base d'antennes pour l'azimut, le site et la fréquence porteuse de l'onde électromagnétique incidente émise par une source. Cette information est ensuite exploitée par le traitement de démultiplexage spatial. Par ailleurs, les techniques de séparation à l'ordre 2 utilisant la cyclostationnarité imposent sur les signaux à séparer soit des différences de rythmes symbole soit des décalages de fréquences porteuses pour fonctionner. Ces limitations rendent ces techniques inadaptées aux configurations de canaux formés dans des radiocommunications entre une station de base et des mobiles présentant des multitrajets décorrélés qui sont non séparables par les méthodes précédentes et donc non opérationnelles dans la plupart des zones urbaines. Dans la suite de la description, les termes "capteurs " et "antennes" seront utilisés indifféremment ainsi que les termes "multicapteurs", "multivoie", et "réseau multicapteur" et base d'antennes" sachant d'autre part qu'une base d'antennes, ou réseau multicapteur, forme un nombre de voies de transmission égal au nombre d'antennes, ou de capteurs. L'invention a pour but de pallier les imperfections et limitations des méthodes précitées. A cet effet, l'invention a pour objet un procédé de multiplexage/démultiplexage spatial de signaux radioélectriques organisés en trames et transmis en duplex entre au moins une station de base et plusieurs mobiles communiquant sur un même canal fréquentiel et temporel au moyen d'un récepteur et d'un émetteur multivoie intégrés dans la station de base et couplés à une base d'antennes, caractérisé en ce qu'il consiste pour des signaux radioélectriques non gaussiens:: - a estimer l'information spatiale relative à chaque mobile à partir du signal reçu par le récepteur multivoie, pour les fréquences de réception et d'émission, au moyen de méthodes de séparation de sources en aveugle, et à partir de cette information, - a isoler en présence de multitrajets dans le canal, par filtrage spatial, les trajets respectifs à chaque mobile possédant un vecteur directeur dont la puissance est supérieure à un seuil déterminé pour permettre le démultiplexage spatial, et - à émettre simultanément en direction des trajets les plus puissants parmi les trajets isolés, les signaux à destination des mobiles correspondants, en protégeant chaque mobiles vis à vis des signaux émis vers les autres par filtrage spatial, pour permettre le multiplexage spatial. Un avantage de la présente invention est de permettre, pour l'estimation des vecteurs directeurs associés aux différents signaux mobiles, d'appliquer une méthode de séparation de sources en aveugle ne nécessitant ni la maîtrise du front d'onde ni connaissance à priori sur la forme d'onde (méthode transparente). Le procédé selon l'invention exploite pour cela le caractère non Gaussien des signaux sources. Cette hypothèse est vérifiée pour la plupart des modulations utilisées dans les communications radio analogiques ou numériques notamment pour les modulations à enveloppe constante. D'autres avantages et caractéristiques de la présente invention apparaîtront plus clairement à la lecture de la description qui suit et des figures annexées qui représentent: - la figure 1, les principales étapes du procédé selon l'invention, - la figure 2, un schéma fonctionnel d'un dispositif pour la mise en oeuvre du procédé selon l'invention, - la figure 3, un exemple de structure d'une "supertrame" prenant en compte les contraintes du mode AMRS, - la figure 4, un schéma fonctionnel précisant l'interface entre le contrôleur AMRS et les moyens d'estimation et de multiplexageldémultiplexage spatial, - la figure 5, un schéma de la décomposition d'un canal fréquentiel et temporel en sous-canaux spatiaux, - la figure 6, une illustration de la notion de canal fréquentiel et temporel, - la figure 7, une table de corrélation spatiale des mobiles présents dans la cellule (pris deux à deux), et - la figure 8, un schéma fonctionnel d'un module de commutation et de tarage du dispositif pour la mise en oeuvre du procédé selon l'invention. Sur les figures, les éléments homologues sont désignés par le même repère. Le procédé selon l'invention dont les principales étapes sont illustrées à la figure 1 consiste, à partir d'une émission/réception multivoie 1, à estimer 2 les informations spatiales concernant chaque source mobile active présente dans la cellule aux deux fréquences, supports de l'émission et de la réception altemativement, puis à partir de ces informations à réaliser un multiplexage/démultiplexage global 3 (temporel, fréquentiel et spatial) optimal des communications. L'information spatiale caractérisant un mobile est définie par les vecteurs directeurs des multitrajets du signal reçu sur la base d'antennes. Ces vecteurs directeurs représentent la réponse (amplitude et phase) de la base d'antennes à un signal incident suivant un azimut caractéristique de la position de la source. L'estimation 2 des vecteurs directeurs de chaque source utilise des signaux reçus de la source d'intérêt en configuration mono-émission. Autrement dit dans la forme d'onde organisée en trames, une plage temporelle périodique, appelée par la suite "burst", est allouée à chaque mobile pour la réactualisation de son information spatiale. La méthode appliquée pour cette estimation 2 exploite donc le caractère non gaussien des signaux sources à traiter et fonctionne parfaitement en l'absence de calibration de la base d'antennes. Par contre, compte-tenu de l'écart généralement important entre les fréquences des liaisons montantes (mobiles vers station de base) et descendantes (station de base vers mobiles), l'estimation réalisée sur les signaux reçus et utilisée pour un filtrage spatial en réception n'est pas transposable pour le filtrage spatial en émission.C'est pourquoi, dans le procédé selon l'invention une permutation périodique des groupes de fréquences des liaisons montantes et descendantes peut être effectuée de façon à permettre l'apprentissage et la réactualisation de l'information spatiale également aux fréquences d'émission. La plupart des stations de base gèrent plusieurs canaux fréquentiel en mode duplex et dans ce contexte le procédé selon l'invention peut permettre d'assurer une gestion globale de ces canaux à l'intérieur desquels le multiplexage spatial génère des sous-canaux AMRS de façon à regrouper préférentiellement sur un même canal fréquentiel les mobiles les moins corrélés spatialement. De plus, le procédé selon l'invention réactualise les coefficients des filtres spatiaux en réception et en émission, réalisant respectivement un démultiplexage spatial et un multiplexage spatial en prenant en compte les défauts d'appairage, conséquence d'erreurs résiduelles de tarage des voies en émission et en réception. Un dispositif pour la mise en oeuvre du procédé selon l'invention est schématisé par la figure 2 et comporte les sous-ensembles suivants: - une base d'antennes 4 couplée à des moyens 5 de commutation de fréquence et de tarage, -un récepteur multivoie 6 comportant des moyens de transpositionlamplification radio, et des moyens de conversion analogique/numérique multivoies non représentés, - un démultiplexeur spatial 7, -un nombre déterminé de démodulateurs 8i égal au nombre maximal de canaux spatiaux fixé par canal fréquentiel, par exemple trois démodulateurs pour une base comportant dix antennes, - un émetteur multivoie 9 comportant des moyens de transpositionlamplification radio et de conversion numériquelanalogique multivoies non représentés, - un multiplexeur spatial 10, -un nombre déterminé de modulateurs Ili égal au nombre de démodulateurs 8i du récepteur multivoie 6, - un module de traitement numérique 12, délimité par une ligne fermée en trait discontinu, comportant un contrôleur 13 et un estimateur 14 spécifiques au traitement AMRS, qui est au coeur du dispositif pour la mise en oeuvre du procédé selon l'invention. La mise en oeuvre du concept de I'AMRS au moyen des méthodes de séparation de sources en aveugle, par le matériel qu'il nécessite, les contraintes des traitements appliqués et les nouvelles possibilités qu'il offre, nécessite une prise en compte système. Ceci se traduit par un certain nombre de tâches de contrôle et de gestion indispensables à une intégration optimale du multiplexage spatial dans les réseaux de communications radiocellulaires. Ces tâches de contrôle et de gestion sont assurées dans le dispositif pour la mise en oeuvre du procédé selon l'invention par le module de traitement numérique 12. La séparation, au niveau de la station de base, des signaux en provenance ou à destination des différents mobiles "actifs" c'est-àdire en communication avec la station de base, dans une cellule déterminée, implique l'acquisition et la réactualisation de l'information spatiale spécifique de chaque mobile et nécessite pour cela une synchronisation des traitements. Pour réaliser cette fonction, le procédé selon l'invention s'appuie sur une forme d'onde de type AMRT avec une segmentation temporelle du canal: le signal reçu par la station de base, ainsi que le montre la figure 3, est organisé en trames , trame Ul à trame &num;K, comportant chacune un nombre déterminé M de "bursts". Chaque "burst", classiquement alloué à un utilisateur grâce au multiplexage spatial est partagé par plusieurs communications à l'exception d'un "burst" par trame qui reste alloué spécifiquement à un mobile.Ces "bursts" périodiques dénommés "bursts" de localisation AMRS sont donc alloués successivement à tous les canaux spatio-temporels potentiellement utilisables sur le canal fréquentiel (par canal fréquentiel est sous-entendu une paire de fréquences émission'réception) La structure d'une "supertrame", apparaissant sur la figure 3, correspond à la période de réactualisation de l'information spatiale de tous les canaux spatio-temporels soit un nombre K de trames égal au double du nombre de ces canaux (pour chaque canal spatio-temporel, l'information spatiale est acquise sur deux trames: une pour chacune des fréquences du canal fréquentiel). Suivant la durée de la "supertrame" en regard de la stationnarité de l'information spatiale, cette méthode doit être adaptée en allouant par exemple plus d'un "burst" par trame à la localisation AMRS pour réduire la période de réactualisation. Il est évident que le choix du nombre de "bursts" de localisation par trame se répercute directement sur les capacités du système. La synchronisation des traitements dans ce contexte consiste essentiellement à orienter les "bursts" reçus soit vers la mesure d'information spatiale ("burst" de localisation d'un mobile actif), soit vers le démultiplexage spatial ("burst" de trafic), ou, soit nulle part ("burst" de localisation inutilisé), et à distribuer ensuite l'information spatiale aux filtres de multiplexage/démultiplexage en émission/réception en fonction des affectations et réaffectations de canaux gérées par le contrôleur AMRS 13. La figure 4 illustre un schéma fonctionnel de l'interface entre le contrôleur AMRS 13, l'estimateur AMRS 14 et les multiplexeur/démultiplexeur spatiaux émission et réception 10 et 7. L'orientation des "bursts" est effectuée par l'intermédiaire d'un commutateur 15 recevant sur une première entrée le signal reçu du récepteur 6, délimité par une ligne fermée discontinue, après transposition BF. Le commutateur 15 reçoit également sur une deuxième entrée un signal de contrôle CTRL délivré par le contrôleur AMRS 13. Il comporte une première, une deuxième et une troisième sorties correspondant respectivement aux "burst" de localisation d'un mobile actif, au "burst" de trafic et au "burst" de localisation inutilisé (mobile inactif). La première sortie est couplée à l'entrée de l'estimateur AMRS 14 estimant le vecteur directeur Sj de la source i correspondant à un mobile actif i déterminé. Le vecteur directeur estimé S est injecté au contrôleur AMRS 13. La deuxième sortie est inactive et la troisième sortie est couplée à l'entrée du démultiplexeur spatial 7, délimité par une ligne fermée discontinue, constitué des filtres spatiaux W1, w; et VV+3 recevant respectivement le même nombre N de sous-canaux spatiaux contenu dans le "burst" trafic reçu. "+" en exposant correspond à la fonction de transposition conjuguée.Les sorties des filtres spatiaux Wï+, W2+ et Vs sont couplées respectivement à l'entrée des démodulateurs 8i DEMOD &num;1 à &num;3, délimités par une ligne fermée discontinue, et délivrent les informations relatives aux mobiles sur le réseau de stations de base représenté sur la figure par deux traits parallèles. Une réactualisation des poids en réception de chaque coefficient des filtres spatiaux W1 . W2 et w; est réalisée à partir du vecteur directeur S estimé par l'estimateur AMRS 14 et après contrôle par le contrôleur AMRS 13. En émission, les modulateurs Ili, MOD &num;1 à MOD &num;3, délimités par une ligne fermée discontinue, reçoivent respectivement des données issues du réseau de stations de base, à émettre en direction des mobiles. Chaque signal modulé est injecté respectivement à l'entrée de chaque filtre spatial W1. W2 et W3 en émission constituant le multiplexeur spatial 10 délimité par une ligne fermée discontinue. De la même façon qu'en réception, une réactualisation du poids en émission des coefficients des filtres spatiaux Wi , W2 et W3 est effectuée à partir du même vecteur directeur S estimé par l'estimateur AMRS 14 et contrôlé par le contrôleur AMRS 13. Les sous-canaux sortant respectivement de chaque filtre spatial W1. W2 et W3 sont sommés puis émis après transposition en fréquence intermédiaire Fl dans l'émetteur multivoie 9 délimité par une ligne fermée discontinue. Le filtre spatial adapté à l'émission vers un mobile déterminé est calculé à partir de son vecteur directeur S estimé sur le signal reçu. La nécessité d'éloigner suffisamment les fréquences d'émission et de réception notamment pour autoriser l'utilisation d'une base d'antennes commune pour les deux sens de transmission émission/réception peut rendre inutilisable en émission les vecteurs directeurs des sources estimés sur le signal reçu. Pour résoudre ce problème, une permutation des fréquences d'un même canal fréquentiel (émission/réception) est opérée par le procédé selon l'invention, aboutissant au niveau de la station de base à une permutation globale des sous-bandes allouées aux liaisons montantes et descendantes. Cette permutation réalisée au rythme trame permet l'estimation de l'information spatiale sur la fréquence utilisée en émission avec une trame de décalage par rapport à la réception, ce qui signifie que dans la "supertrame" deux "bursts" de localisation consécutifs sont utilisés pour effectuer une estimation complète de l'information spatiale sur chaque canal spatio-temporel. Un contrôleur de permutation 16 des fréquences est décrit dans la suite de la description en regard de la figure 8. Avec le multiplexage spatial, chaque canal, classiquement défini comme le support fréquentiel et temporel d'une transmission entre une station de base SB et un mobile, se décompose en sous-canaux spatiaux. La figure 5 illustre une schématisation de cette décomposition où trois canaux X, Y et Z supportent les transmissions avec respectivement trois mobiles MBxl, MBx2, MBx3, deux mobiles MByl, MBy2 et un mobile MBzl. La figure 6 donne une illustration de la notion de canal fréquentiel et temporel pour les canaux X, Y, Z dans le cas de liaisons montantes et dans le cas de liaisons descendantes. Pour effectuer une affectation de canal, le procédé selon l'invention se base sur une table de corrélation spatiale inter-mobile. Cette table dont un exemple est présenté sur la figure 7, contient les coefficients de corrélation spatiale CS de chacun des mobiles avec tous les autres à l'intérieur de la flotte gérée par la station de base. Le coefficient de corrélation spatiale CSij entre deux mobiles i et j est établi à partir de leurs vecteurs directeurs aib, âjb (estimés sur la fréquence basse Fb du canal) et âih, àjh (estimés sur la fréquence haute Fh du canal) suivant la formule: En présence de multitrajets, c'est le vecteur directeur du trajet le plus puissant qui représente l'information spatiale représentative du mobile au regard de l'affectation de canal. Sur la figure 7, les mobiles sont regroupés de façon à faire apparaître l'ensemble des corrélations spatiales à l'intérieur de chaque canal respectivement pour le canal X, le canal Y et le canal Z (zones entourées de pointillés) dont les valeurs sont plus particulièrement prises en compte dans la méthode d'affectation de canal. Les lignes et colonnes hachurées correspondent à des sous-canaux spatiaux non utilisés. A chaque nouvelle estimation d'un vecteur directeur relatif à un mobile actif (après l'analyse de chaque "burst" de localisation AMRS), la ligne et la colonne associées de la table sont réactualisées et des tests d'affectation de canal sont appliqués suivant un algorithme d'affectation déterminé. Son principe en est le suivant : une réactualisation de la répartition des mobiles sur les différents canaux de la cellule est effectuée toutes les deux trames après analyse des deux "bursts" de localisation alloués à chaque mobile successivement sur les deux fréquences du canal duplex. Si les dernières informations spatiales mesurées sont associées à un mobile entrant dans la cellule, le contrôleur AMRS applique au nouveau mobile p le test d'admission suivant: pour que le mobile p soit admis sur le canal non saturé Cn, il faut et il suffit que quelque soit le mobile x déjà présent sur l'un des sous-canaux spatiaux de Cn: CSxp seuil (2) si cette condition n'est remplie sur aucun canal non saturé (par rapport à ces capacités potentielles, c'est à dire au nombre maximal de sous-canaux spatiaux) alors le mobile ne peut être pris en charge par la station de base. Dans le cas contraire, le mobile est admis sur le canal Cn tel que: Max#cn(CSpx) = Mini=1,2-..,NbCan [Maxy e c; (CSpy )] (3) Si les dernières informations spatiales mesurées sont associées à un mobile p déjà présent sur la cellule, alors le contrôleur AMRS commence par vérifier si le mobile reste suffisamment décorrélés spatialement des autres mobiles présents sur le même canal par application du test suivant: CSxp seuil pour tout mobile x du canal Cn (4) si cette condition est remplie, l'affectation des canaux n'est pas modifiée. Dans le cas contraire, le moyen de contrôle AMRS teste la possibilité de permuter les affectations de canaux entre le mobile p et un mobile q sur un canal Cm différent.Pour que cette permutation soit possible, il faut et il suffit que les deux conditions suivantes soient simultanément réalisées: CSyp ( Seuil pour tout mobile y du canal Cm (5) CSqx seuil pour tout mobile x du canal Cn (6) Si aucune permutation n'est possible, le mobile p ne peut plus être pris en charge par la station de base. Si, et cela doit être le cas le plus fréquent, plusieurs permutations sont possibles le critère de choix optimal proposé est le suivant: le mobile q retenu pour réaliser la permutation avec le mobile p est celui qui minimise l'expression: où Cm désigne le canal affecté au mobile q. Ainsi que la plupart des techniques de filtrage spatial de signaux multivoies n'utilisant pas de séquence d'apprentissage, la méthode proposée n'est applicable que si les chaînes d'émission et de réception n'introduisent ni déphasage, ni gain différentiels et ni distorsions (non linéarités d'amplification). Il est très difficile et très coûteux de réaliser un dispositif satisfaisant ces contraintes, c'est pourquoi une méthode est décrite ci-après pour corriger par traitements numériques les dérives et distorsions sur les différentes voies émission et réception. La correction des déphasages et gain différentiels pour chaque voie formée par chaque antenne de la base d'antennes 4, est réalisée par une procédure de tarage se décomposant en deux étapes. Cette procédure est mise en oeuvre par un dispositif dont un schéma fonctionnel est illustré à la figure 8. Dans une première étape, le tarage de la chaîne de réception 17i de chaque voie i comportant des moyens de transposition de fréquence, de filtrage, d'amplification et de conversion analogique/numérique, est réalisée par injection en amont de la chaîne de réception 1 7i d'un signal de tarage, puis adaptation des coefficients d'un filtre d'égalisation 18i couplé à la sortie de la chaîne de réception 1 7i (sauf la première prise comme référence) permettant de compenser la dispersion des réponses impulsionnelles des différentes voies i, par un analyseur de tarage 19 commun à toutes les voies i formées par la base d'antennes 4. Dans une deuxième étape, le tarage de la chaîne d'émission 20i de chaque voie i comportant des moyens de conversion numérique/analogique, d'amplification, de filtrage et de transposition de fréquence, est réalisée de façon analogue, en injectant également un signal de tarage en amont de la chaîne d'émission 20i. Une partie du signal émis est réinjectée sur la chaîne de réception 17i, préalablement tarée, via un coupleur résistif 21 i. Enfin, un filtre d'égalisation 22i couplé en entrée de la chaîne d'émission 20i est adapté suivant le même principe que celui de la chaîne de réception 17i. Sur cette figure est également représenté le contrôleur 16 de la permutation des fréquences. Un commutateur deux voies assure l'émission et la réception respectivement à la fréquence basse Fb, et la fréquence haute Fh du signal reçu. Ces deux fréquences Fb et Fh étant séparées au moyen d'un duplexeur 24i disposé entre le coupleur résistif 21i et le commutateur deux voies 23i. Le contrôleur 16 de permutation des fréquences Fb et Fh est couplé aux chaînes émission et réception 20i et 17i de chaque voie i, à l'analyseur de tarage19 et au commutateur deux voies 23i. Un premier commutateur une voie 25i disposé à l'entrée de la chaîne de réception 17i reçoit sur une première entrée le signal de tarage réception, sur une deuxième entrée le signal prélevé par le coupleur résistif 21 i et sur une troisième entrée le signal réception à la fréquence basse Fb ou fréquence haute Fh en fonction de la position du commutateur deux voies 23i. Un deuxième commutateur une voie 26i est couplé en sortie du filtre d'égalisation de réception 18i et commande le tarage en réception. Un troisième commutateur une voie 27i, disposé à l'entrée de la chaîne d'émission 20i avant le filtre d'égalisation 22i reçoit sur une première entrée le signal à émettre et sur une deuxième entrée le signal de tarage en émission. Le signal de tarage utilisé est composé de N raies sinusoïdales équiespacées dans la bande utile. Le principe de l'adaptation des p coefficients des filtres d'égalisation 22i et 18i des chaînes d'émission et de réception 20i et 17i est basé sur la minimisation du critère suivant: m' - F. wiII2 (8) H1 où m'i = Hi(Hi réponse fréquentielle de la voie i estimée aux fréquences du signal de tarage), F.wj est la réponse fréquentielle du filtre d'égalisation de la voie i, avec F matrice Nxp telle que F(m,n) = e-2.j.##m.n , et une des fréquences composant le signal de tarage. La solution au sens des moindres carrés pour le filtre d'égalisation de la voie i est donc: wi=[F+.F]-1 F+.m' (9) La formation de voie en émission ainsi que l'indépendance des sources font que les signaux à émettre en mode AMRS n'ont pas une enveloppe constante et que leur dynamique peut dépasser la plage de linéarité des amplificateurs des chaînes d'émission 20i. Pour éviter des distorsions désappairées entre voies, un filtre de pré-correction 28i est disposé à l'entrée de la chaîne d'émission 20i et dont les caractéristiques sont les suivantes: en modélisant les distorsions d'amplitude et de phase d'un amplificateur comme suit: les fonctions de pré-correction f et g' établies après apprentissage des fonctions f et g sont définies par: où G désigne un éventuel gain fixe de l'ensemble filtre de pré-correction + amplificateur. La technique mise en oeuvre pour estimer l'information spatiale des sources mobiles est l'une des clés essentielles de cette invention. Cette technique très prometteuse dite de "séparation de sources en aveugle' exploite donc le caractère non gaussien et indépendant des signaux à séparer. La séparation de sources en aveugle permet de prendre en compte des multitrajets corrélés ou non temporellement dans l'estimation des vecteurs directeurs, ce qui n'est pas le cas des techniques classiquement appliquées en goniométrie haute résolution ni celui des techniques de séparation exploitant la cyclostationnarité. Les multitrajets corrélés (temporellement) sont intégrés dans un seul vecteur directeur composite tandis que les multitrajets décorrélés sont séparés par la méthode aboutissant à l'estimation d'un vecteur directeur par trajet indépendant. Il est rappelé que l'algorithme d'estimation de l'information spatiale est appliqué à des signaux en configuration mono-émission (sur les "bursts" de localisation). Les traitements se décomposent comme suit: - estimation de l'espace source (associés aux multitrajets décorrélés temporellement), - blanchiment des signaux sources, et - estimation des vecteurs directeurs blanchis. L'estimation de l'espace source est effectué via une décomposition en vecteurs propres, valeurs propres (méthode QR, de Jacobi ou autre) de la matrice Hermitienne de corrélation des observations x(t): Rx = E x(t).x(t)+ = U.D.U+ (12) où E[x(t).x(t)+] correspond à l'espérance mathématique, avec U matrice des vecteurs propres, et D matrice diagonale des valeurs propres. La matrice de corrélation des signaux sources est définie alors par: Rs = Us.Ds.Us+ (13) avec Us matrice des vecteurs propres de l'espace source, et D5 matrice diagonale des valeurs propres de l'espace source. L'isolation de l'espace source est opérée par application d'un seuil minimal sur les valeurs propres. La matrice de blanchiment des sources B s'écrit: Le blanchiment des signaux consiste ensuite simplement à appliquer la matrice de blanchiment au vecteur d'observation, ce qui donne pour les données blanchies la formule suivante: La dernière étape consiste à estimer H correspondant à la matrice unitaire des vecteurs directeurs orthonormalisés des sources. Plusieurs algorithmes existent pour effectuer cette estimation parmi eux est proposé celui de COMON et de CARDOSO-SOULOUMIAC les plus robustes testés, consistant à résoudre le problème d'optimisation suivant: Si on pose z = H,,+. y (16) COMON propose comme estimateur de H la matrice unitaire Hmc minimisant la somme de tous les carrés des modules des cumulants croisés d'ordre 4 de z, ce qui revient à trouver la matrice H mc maximisant le facteur de contraste définit par COMON à savoir: (méthode détaillée dans l'article de P. COMON : "Independent Component Analysis - Signal Processing", Vol. 36, N" 3, Special Issue On Higher Order Statistics, pp 287-314, Avril 1994). CARDOSO et SOULOUMIAC proposent comme estimateur de H la matrice Hmc minimisant la somme des carrés des modules des cumulants croisés d'ordre 4 de z ayant des premiers et seconds indices différents, ce qui revient à trouver la matrice Hmc maximisant la fonction Ccar.sou. (méthode détaillée dans l'article de J.F. CARDOSO et A. SOULOUMIAC "Blind Beamforming for Non Gaussian Signals", IEE Proc-F, Vol. 140, N" 6, pp 362-370, Décembre 1993). L'estimation finale des vecteurs directeurs À non blanchis utilisés au niveau du filtrage spatial s'obtient simplement par application de la matrice inverse de la matrice de blanchiment: Â = B-1.H = US.Ds1/2.H (19) Pour le démultiplexage spatial en réception, le filtrage spatial consiste à optimiser pour chaque mobile le diagramme d'antenne en fonction de son information spatiale, c'est dire en formant des lobes dans la direction des multitrajets du signal mobile visé et des trous dans les directions des autres signaux mobiles. La station de base en mode AMRS applique donc autant de filtres spatiaux qu'il y a de mobiles actifs sur le même canal. L'application du filtre spatial illustrée à la figure 4 se résume à la formule suivante: yj(t) = wi+.x(t) (20) où y(t) représente le signal filtré spatialement sur lequel la contribution du mobile a été amplifiée et celle des autres mobiles fortement atténuée. La méthode proposée d'adaptation des coefficients du filtre spatial de réception porte le nom de filtrage adapté spatial estimé, ou FAS estimé, et consiste à réactualiser le filtre spatial adapté au mobile i en fonction de la matrice de corrélation estimée des observations et de l'information spatiale suivant la formule: wi = (Rx + i21d)-1.ai (21) avec: : Wi le vecteur poids adapté à la séparation de la source i, I'estimation de la matrice de corrélation des observations, le vecteur directeur composite du mobile i construit par sommation des vecteurs directeurs des multitrajets pondérés par leur puissances relatives, cl;2, le terme de bruit fictif rajouté pour rendre l'estimation plus robuste vis à vis des risques de réjection du signal utile, liés aux erreurs commises sur âi (le rajout de ce terme de bruit n'est pas crucial pour la méthode de séparation de source en aveugle, puisque les vecteurs sources sont directement et donc précisément estimés à partir des signaux reçus alors que dans le cas de la goniométrie chaque vecteur directeur est déduit d'une estimation de direction d'arrivée via des tables de calibration), et Id, correspond à la matrice identité. Pour le multiplexage spatial en émission, le principe du filtrage spatial est de réaliser l'opération réciproque du filtrage spatial en réception, c'est à dire au lieu de séparer des signaux arrivant sur la base d'antennes suivant différents azimuts, il s'agit avec cette meme base de générer un signal composite dont chaque composante adressée à un mobile particulier se propage suivant la direction définie par son vecteur directeur principal. A la différence du filtrage en réception où, pour optimiser le rapport signal sur bruit, il est important de tenir compte de tous les multitrajets d'une même source, en émission il est préférable de ne focaliser l'énergie que dans la direction du multitrajet reçu le plus puissant. L'application des pondérations aux signaux à émettre issus des différents modulateurs MOD#1 à &num;3 Ili comme illustrée sur la figure 4, s'écrit comme suit: x= W.m (24) avec: m = (m1(t), m2(t),..., mM(t)) signaux issus des M modulateurs, x = (x1(t), x2(t) z , xN(t)) signaux émis sur les N antennes, W matrice des vecteurs wj (i=1 à M). La méthode d'adaptation des coefficients du filtre spatial en réception est théoriquement applicable en émission, mais en raison des erreurs résiduelles de tarage des voies en émission et en réception, des erreurs dites "de recopie" sont commises qui occasionnent des erreurs de pointage avec un mauvais positionnement des "trous" dans le diagramme de rayonnement et le risque au niveau de la réception par le mobile d'un mauvais rapport Signal/Brouilleur sans possibilité de réjection. C'est pourquoi, I'adaptation en émission est effectuée suivant une technique légèrement différente dénommée FAS synthétique. Le calcul de wi, filtre spatial en émission adapté au mobile i, est obtenu en appliquant à son vecteur directeur principal (associé au multitrajet le plus puissant) I'inverse, non plus de la matrice de corrélation estimée des observations sur les N antennes, mais d'une matrice de corrélation synthétique des signaux brouilleurs RJj: Cette matrice est construite à partir des vecteurs directeurs de deux brouilleurs fictifs positionnés de part et d'autre (en azimut) de chaque brouilleur réel (chaque mobile constituant un brouilleur pour tous les autres). L'intéret de cette méthode est de permettre une légère erreur sur les vecteurs directeurs des brouilleurs en élargissant les "trous" du diagramme de rayonnement autour de leurs directions estimées. La matrice de corrélation synthétique utilisée pour l'estimation du vecteur des poids wi (focalisation en direction du mobile i) s'écrit sous la forme: avec: ajk représentant les vecteurs directeurs des deux brouilleurs fictifs: ajk = âj(#+k.##) (27) En l'absence d'information sur l'azimut de la source e et sur la variation d'un vecteur directeur avec l'azimut de sa source (sans table de calibration), ces vecteurs sont obtenus par déphasage du vecteur directeur estimé âj(#) : âjk =D.aj() (28) avec Dk matrice diagonale de déphasage dont les coefficients diagonaux sont de la forme: dk (i) = ej.k(i-N).## (29) où:: A(p est un déphasage de l'ordre de quelques degrés à adapter en fonction de la base d'antennes 4, llj désigne la puissance de la source j considérée comme brouilleur dans ce cas, ce paramètre est fixé à une valeur proportionnelle à la profondeur du "trou" qu'il est nécessaire de former dans la direction du mobile j sur le diagramme de rayonnement (typiquement 30 dB au dessus du signal utile). ai2 est un terme de bruit ajouté pour régler la profondeur du "trou" dans le diagramme de rayonnement, son ajout a également pour conséquence de rendre la matrice RJj inversible et d'éviter le recours à un calcul coûteux de matrice pseudo-inverse. La présente invention s'applique plus particulièrement aux stations de base bien que certains aspects du procédé selon l'invention qui vient d'être décrit puissent trouver des applications au niveau des mobiles notamment en ce qui concerne l'augmentation de leur sensibilité en réception. Le principe de la présente invention n'est pas limité à des mobiles mais peut s'appliquer également à des stations secondaires fixes déportées par rapport à la station de base. REVENDICATIONS 1. Procédé de multiplexageldémultiplexage spatial de signaux radioélectriques organisés en trames et transmis en duplex entre au moins une station de base et plusieurs mobiles communiquant sur un meme canal fréquentiel et temporel au moyen d'un récepteur et d'un émetteur multivoie (1) intégrés dans la station de base et couplés à une base d'antennes, caractérisé en ce qu'il consiste, pour des signaux radioélectriques non gaussiens - a estimer (2) I'information spatiale relative à chaque mobile à partir du signal reçu par le récepteur multivoie, pour les fréquences de réception et d'émission, au moyen de méthodes de séparation de sources en aveugle, et à partir de cette information, - à isoler (3) en présence de multitrajets dans le canal, par filtrage spatial, les trajets respectifs à chaque mobile possédant un vecteur directeur dont la puissance est supérieure à un seuil déterminé pour permettre le démultiplexage spatial, et - a émettre simultanément en direction des trajets les plus puissants parmi les trajets isolés, les signaux à destination des mobiles correspondants, en protégeant chaque mobiles vis à vis des signaux émis vers les autres par filtrage spatial, pour permettre le multiplexage spatial. 2. Procédé selon la revendication 1, caractérisé en ce que l'estimation (2) de l'information spatiale est effectuée en configuration mono-émission, sur des plages de temps périodiques, aménagées dans une structure de trame déterminée, et allouées spécifiquement à l'estimation de l'information spatiale de chaque mobile sur deux trames consécutives afin de prendre en compte les deux fréquences de la transmission duplex. 3. Procédé selon l'une quelconque des revendications 1 et 2, caractérisé en ce qu'il consiste à permuter à chaque trame les fréquences émission et réception pour permettre l'acquisition de l'information spatiale sur les deux fréquences lorsqu'elles sont trop éloignées pour que cette information soit transposable de l'une à l'autre, 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'il consiste pour optimiser l'efficacité du multiplexage et du démultiplexage spatial, à gérer l'affectation des canaux au niveau de l'ensemble des ressources fréquentielles et temporelles allouées à la station de base en fonction des mouvements des mobiles, en utilisant un critère de corrélation spatiale entre mobiles intégrant les informations spatiales acquises sur les deux fréquences utilisées alternativement pour la transmission duplex. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que le filtrage spatial en émission du signal, à destination d'un mobile déterminé, est adapté à partir d'une matrice de corrélation synthétique de signal bruit+brouilleurs, établie à partir de paires de brouilleurs fictifs disposés spatialement de part et d'autre des mobiles en direction desquelles la station de base ne doit pas émettre ce signal, de façon à élargir les "trous" du diagramme de rayonnement en direction des brouilleurs et par conséquent à augmenter la tolérance du système aux erreurs de pointage consécutives aux erreurs résiduelles de tarage des voies d'émission et de réception. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il consiste à appliquer périodiquement une procédure de tarage des chaînes d'émission et de réception de chaque voie formée par la base d'antennes et de pré-correction d'amplification à l'émission afin de limiter les erreurs de pointage en émission, 7.Dispositif pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce qu'il comporte: - une base d'antennes émission/réception (4) couplée à des moyens de commutation et de tarage (5), -un récepteur multivoie (6) comportant des moyens de transposition de fréquence, d'amplification et de conversion analogique/numérique, couplé en sortie des moyens de commutation et de tarage (5), -un démultiplexeur spatial (7), couplé en sortie du récepteur multivoie (6), comportant un nombre déterminé de filtres spatiaux, chaque filtre spatial formant un canal spatial, - un nombre déterminé de démodulateurs (8i) égal au nombre de canaux spatiaux déterminés par canal fréquentiel et temporel, couplés en sortie du démultiplexeur (7), et délivrant sur le réseau des stations de base les données émises par les mobiles, - un nombre déterminé de modulateurs (iii) égal au nombre de démodulateurs (8i), recevant respectivement des données reçues du réseau des stations de base et à émettre par la base d'antennes (4) à destination des mobiles, - un multiplexeur spatial (10), couplé en sortie des modulateurs (ils), comportant un nombre déterminé de filtres spatiaux égal au nombre de modulateurs (11 i), - un émetteur multivoie (9) comportant des moyens de conversion numérique/analogique, de transposition de fréquence et d'amplification, couplé en sortie du multiplexeur spatial (10), et en ce qu'il comporte:: - un module de traitement numérique AMRS (12), comportant un contrôleur AMRS (13) et un estimateur (14) des vecteurs directeurs des trajets relatifs aux mobiles communiquant avec la station de base, permettant la réactualisation des poids des filtres spatiaux à l'émission et à la réception, la gestion des fréquences d'émission et de réception, et les commandes des moyens de communication et de tarage (3). 8. Dispositif selon la revendication 7, caractérisé en ce que les moyens de tarage comportent pour chaque voie (i), comportant respectivement une chaîne d'émission (22i) et une chaîne de réception (17i), et formée respectivement par chaque antenne de la base d'antennes (4): - un filtre d'égalisation du signal d'émission (22i) couplé en entrée de la chaîne d'émission (20i), -un filtre d'égalisation (18i) du signal de réception couplé en sortie de la chaîne de réception (17i), - un commutateur en entrée (25i) et en sortie (26i) de la chaîne de réception (17i), un commutateur (27i) en entrée de la chaîne d'émission (20i) et un commutateur deux voies (23i) en sortie de la chaîne émission (20i), permettant au dispositif de passer dans les différents modes de fonctionnement: tarage réception, tarage émission, émission/réception, - un coupleur (21 i) disposé entre l'antenne (Ai) et le commutateur deux voies (23i), permettant de prélever une partie du signal d'émission pour l'injecter à l'entrée de la chaîne de réception (17i) via le commutateur d'entrée (25i) de la chaîne de réception (17i), et en ce qu'il comporte en commun pour toutes les voies (i): - un analyseur de tarage (19) permettant l'analyse des différents signaux de tarage arrivant sur les commutateurs d'entrée (27i) et de sortie (26i) respectivement des chaines d'émission (20i) et de réception (17i), et permettant d'adapter les filtres d'égalisation (22i et 18i) des chaînes d'émission et de réception (20i et 17i). 9. Dispositif selon la revendication 8, caractérisé en ce qu'un filtre de pré-correction (28i) de la chaîne d'émission (20i) est disposé entre le filtre d'égalisation (22i) de la chaîne d'émission (20i) et la chaîne d'émission (20i). 10. Dispositif selon l'une quelconque des revendications 7 à 9, caractérisé en ce qu'il comporte en outre, pour la permutation des fréquences émission et réception d'un même canal fréquentiel, un contrôleur de permutation de fréquences (16i), couplé respectivement au récepteur (6) et à l'émetteur (9) multivoie, et aux moyens de commutation et de tarage (4).