Brevet FR3114464A1 — publié le 20220325

On propose un procédé de détection (Det) d’un signal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication. La propagation dans lesdits canaux de communication est caractérisée par au moins une variable. L’ensemble des valeurs pouvant être prises par ladite variable est découpé en une pluralité d’intervalles. Ce procédé de détection comprend les étapes de : - réception (E4) des signaux respectivement transmis dans la pluralité des canaux de communication, et - détection (E20) d’un signal correspondant à une valeur de ladite variable comprise dans un intervalle parmi la pluralité d’intervalles par comparaison, à un seuil prédéterminé, d’une valeur prise par un corrélateur lié à l’intervalle concerné et calculée en fonction des valeurs bruitées. Un procédé d’estimation de canal, un dispositif de détection, un dispositif d’estimation de canal et un programme d’ordinateur sont également décrits. Figure pour l’abrégé : 3 Procédés et dispositifs de détection d’un signal et d’estimation de canal, et programme d’ordinateur associé Domaine technique de l'invention La présente invention concerne le domaine technique des télécommunications. Elle concerne en particulier des procédés et dispositifs de détection d’un signal et d’estimation de canal, ainsi qu’un programme d’ordinateur associé. Etat de la technique Afin d’optimiser la transmission de données dans un système de communication, il est connu de mettre en œuvre un procédé d’estimation de canal visant à avoir une meilleure connaissance de l’état courant des canaux de communication utilisés par le système de communication. Dans ce contexte, l’article ” OMP with Grid-less Refinement Steps for Compressive mmWave MIMO Channel Estimation ” de L. Weiland, C. Stöckle, M. Würth, T. Weinberger and W. Utschick, 2018 IEEE 10th Sensor Array and Multichannel Signal Processing Workshop (SAM) , Sheffield, 2018, pp. 543-547, doi: 10.1109/SAM.2018.8448789, propose à ce sujet une méthode de type OMP (pour « Orthogonal Matching Pursuit ») permettant de caractériser les canaux de propagation dans le cas d’un système de communication de type MIMO (pour « Multiple-Input Multiple-Ou t put »). Présentation de l'invention Dans ce contexte, la présente invention propose un procédé de détection d’un signal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, la propagation dans lesdits canaux de communication étant caractérisée par au moins une variable, l’ensemble des valeurs pouvant être prises par ladite variable étant découpé en une pluralité d’intervalles, ledit procédé comprenant les étapes de : - réception des signaux respectivement transmis dans la pluralité des canaux de communication, et - détection d’un signal correspondant à une valeur de ladite variable comprise dans un intervalle parmi la pluralité d’intervalles par comparaison, à un seuil prédéterminé, d’une valeur prise par un corrélateur lié à l’intervalle concerné et calculée en fonction des valeurs bruitées. Ainsi, l’étape de détection du signal repose sur l’utilisation d’un corrélateur qui permet d’améliorer les performances de détection notamment en ce qui concerne la distinction par rapport au niveau de bruit présent. De plus, cette étape est peu coûteuse en calculs. D’autres caractéristiques non limitatives et avantageuses du procédé de détection conforme à l’invention, prises individuellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : - l’étape de détection comprend une étape de comparaison, audit seuil prédéterminé, d’une valeur maximale parmi une pluralité de valeurs prises par le corrélateur; - l’étape de détection comprend une étape de détermination, pour chaque intervalle, d’une somme de fonction de corrélateurs fenêtrés associés à l’intervalle concerné et fonction desdites valeurs bruitées ; - l’étape de détection comprend une étape de comparaison de la somme de fonction de corrélateurs fenêtrés audit seuil prédéterminé ; et - le seuil prédéterminé est fonction d’un niveau de bruit associé à l’intervalle concerné. L’invention concerne également un procédé d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit procédé comprenant les étapes de : - détection d’un signal par la mise en œuvre d’un procédé de détection tel que défini précédemment, et - estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. L’invention concerne également un procédé d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit procédé comprenant, tant qu’une condition d’arrêt n’est pas obtenue, la répétition d’étapes de : - détection d’un signal par la mise en œuvre d’un procédé de détection tel que défini précédemment, et - estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté, ladite condition d’arrêt correspondant à l’absence de détection d’un signal lors de la mise en œuvre du procédé de détection tel que défini précédemment. Ce procédé d’estimation de canal présente l’avantage de reposer sur des calculs qui ont déjà été effectués dans le cadre du procédé de détection précédemment introduit. Cela permet là aussi de réduire le coût d’exécution de ce procédé d’estimation. D’autres caractéristiques non limitatives et avantageuses du procédé d’estimation conforme à l’invention, prises individuellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : - l’étape d’estimation comprend une étape de détermination d’une valeur intermédiaire de ladite variable comme égale au centre de l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté ; - l’étape d’estimation comprend une étape de détermination d’une valeur intermédiaire de ladite variable par comparaison de corrélateurs fenêtrés correspondant à l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté ; - l’étape d’estimation comprend en outre une étape de détermination d’une fonction ayant une propriété de convexité associée à un produit scalaire d’un vecteur formé des valeurs bruitées et d’un vecteur de dépointage dépendant de ladite variable, le maximum de ladite fonction ayant une propriété de convexité et dudit produit scalaire étant obtenu pour la même valeur de ladite variable ; et - ladite valeur de la variable est déterminée par une étape d’optimisation sur la base de ladite valeur intermédiaire de la variable et utilisant ladite fonction ayant une propriété de convexité. L’invention concerne également un dispositif de détection d’un signal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, la propagation dans lesdits canaux de communication étant caractérisée par au moins une variable, l’ensemble des valeurs pouvant être prises par ladite variable étant découpé en une pluralité d’intervalles, ledit dispositif de détection comprenant : - un module de réception des signaux respectivement transmis dans la pluralité des canaux de communication, et - un module de détection d’un signal correspondant à une valeur de ladite variable comprise dans un intervalle parmi la pluralité d’intervalles par comparaison, à un seuil prédéterminé, d’une valeur prise par un corrélateur lié à l’intervalle concerné et calculée en fonction des valeurs bruitées. L’invention concerne également un dispositif d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit dispositif d’estimation de canal comprenant : - un dispositif de détection d’un signal tel que défini précédemment, et - un module d’estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. L’invention concerne également un dispositif d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit dispositif d’estimation de canal étant conçu pour activer, tant qu’une condition d’arrêt n’est pas obtenue : - un dispositif de détection d’un signal par la mise en œuvre d’un procédé de détection tel que défini précédemment, et - un module d’estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté, ladite condition d’arrêt correspondant à l’absence de détection d’un signal par le dispositif de détection tel que défini précédemment. L’invention propose enfin un programme d’ordinateur comprenant des instructions exécutables par un processeur et conçues pour mettre en œuvre un procédé tel que proposé ci-dessus lorsque ces instructions sont exécutées par le processeur. Bien entendu, les différentes caractéristiques, variantes et formes de réalisation de l'invention peuvent être associées les unes avec les autres selon diverses combinaisons dans la mesure où elles ne sont pas incompatibles ou exclusives les unes des autres. Description détaillée de l'invention De plus, diverses autres caractéristiques de l'invention ressortent de la description annexée effectuée en référence aux dessins qui illustrent des formes, non limitatives, de réalisation de l'invention et où : représente un système de télécommunication mis en jeu dans la présente invention, représente un dispositif d’estimation de canal adapté à mettre en œuvre un procédé d’estimation de canal conforme à l’invention, représente, sous forme d’un logigramme, un exemple de procédé d’estimation de canal conforme à l’invention, représente, sous forme d’un logigramme, un exemple de procédé de détection tel qu’il peut être mis en œuvre dans le cadre du procédé d’estimation de canal représenté sur la figure 3, représente, sous forme d’un logigramme, un exemple de procédé de d’estimation tel qu’il peut être mis en œuvre dans le cadre du procédé d’estimation de canal représenté sur la figure 3, représente l’évolution d’une fonction test utilisée dans le cadre du procédé d’estimation représenté sur la figure 5, représente la forme d’un corrélateur conforme à l’invention, représente la forme d’un noyau de convolution utilisé dans le cadre du procédé d’estimation représenté sur la figure 5, et représente la forme d’une fonction modifiée obtenue lors de la mise en œuvre du procédé d’estimation représenté sur la figure 5. La figure 1 représente un système de télécommunications comprenant un ensemble de Nt antennes d’émission T 1 , …, T Nt et un ensemble de Nr antennes de réception R 1 , …, R N r . Chacune des antennes d’émission T 1 , …, T Nt transmet des signaux électromagnétiques (représentant généralement des données à transmettre, codées par des symboles) dans un canal de communication C associé, où ces signaux sont reçus par les différentes antennes de réception R 1 , …, R N r . Dans l’exemple décrit ici, l’ensemble des antennes d’émission T 1 , …, T Nt est un réseau d’antennes linéaire uniforme (ou ULA pour " Uniform Linear Array "). La transmission des signaux électromagnétiques entre les antennes d’émission T 1 , …, T Nt et les antennes de réception R 1 , …, R N r est ici modélisée physiquement par une onde plane propagée (en ligne directe) selon une direction définie entre le centre de masse des antennes d’émission T 1 , …, T Nt et le centre de masse des antennes de réception R 1 , …, R N r . Cette situation correspond notamment au cas des systèmes dits " MIMO " (pour " Multiple-Input Multiple-Output ") qui comprennent une pluralité d’antennes d’émission T 1 , …, T Nt et une pluralité d’antennes de réception R 1 , …, R N r . Dans ce contexte des systèmes MIMO, la propagation des signaux dans les canaux de communication peut être caractérisée par un angle de départ θ tx , un angle d’arrivée θ r x , et un retard de propagation τ entre la pluralité d’antennes d’émission T 1 , …, T Nt et la pluralité d’antennes de réception R 1 , …, R N r . La présente invention se place dans un contexte de canal de communication multi-trajets (par exemple ici P trajets) entre les antennes d’émission T 1 , …, T Nt et les antennes de réception R 1 , …, R N r . Les signaux reçus y s’écrivent sous la forme suivante : avec X t le vecteur des signaux émis, θ tx l’angle de départ, n un vecteur caractérisant un bruit thermique et la notation T correspondant à l’opérateur de transposition de matrices. Dans ce cas, le canal de communication C est modélisé par le vecteur h s’exprimant sous la forme : où e t est le vecteur de dépointage caractéristique du réseau d’antennes, β p le gain des signaux électromagnétiques et t p la variable définie par t p =cos(θ tx ). Dans une version étendue de modèle MIMO, le vecteur h modélisant le canal de propagation s’exprime sous la forme : avec e t , e r et e f respectivement les vecteurs de dépointage associés aux directions t, r et τ, les variables t=cos(θ tx ), r=cos(θ r x ), τ le retard de propagation et ⊗ le produit de Kronecker. Les canaux de propagation sont donc caractérisés par l’intermédiaire de ces trois variables. L’objectif de la présente invention est de les déterminer. L’exposé suivant décrit la détermination de l’une de ces variables, par exemple, l’angle de départ θ tx , en considérant une modélisation selon une seule des trois variables. En variante, une modélisation à deux variables ou à trois variables peut être utilisée. La figure 2 représente sous forme fonctionnelle un dispositif d’estimation de canal 1 conforme à un exemple de mise en œuvre de l’invention. Ce dispositif d’estimation de canal 1 comprend une unité de commande 2. L’unité de commande 2 comprend également un dispositif de détection 5 et un dispositif d’estimation 7. L’unité de commande 2 comprend un processeur 20 et une mémoire 22. Le dispositif de détection 5 et le dispositif d’estimation 7 sont formés par un ensemble de modules fonctionnels. Par exemple, le dispositif de détection 5 comprend un module de réception de signaux et un module de détection d’un signal. Le dispositif d’estimation 7 comprend quant à lui le dispositif de détection 5 et un module d’estimation. Chacun des différents modules décrits est par exemple réalisé au moyen d’instructions de programme d’ordinateur conçues pour mettre en œuvre le module concerné lorsque ces instructions sont exécutées par le processeur de l’unité de commande 2. La mémoire de l’unité de commande 2 est par exemple conçue pour mémoriser par exemple des signaux pilotes, ici des signaux linéaires, utilisés pour tester les canaux de communication. La figure 3 est un logigramme représentant un exemple de procédé d’estimation de canal mis en œuvre dans le contexte décrit précédemment. Comme le montre la figure 3, le procédé d’estimation de canal débute à l’étape E2, lors de laquelle les antennes d’émission T 1 , …, T Nt émettent des signaux électromagnétiques X t à destination des antennes de réception R 1 , …, R N r . A l’étape E4, les signaux y sont reçus par les antennes de réception R 1 , …, R N r . Ces signaux sont transmis dans la pluralité des canaux de communication. Ces signaux reçus y sont donnés par l’expression précédemment introduite. Le procédé d’estimation de canal se poursuit ensuite par l’étape E6. Lors de cette étape, l’unité de commande 2 détermine des valeurs bruitées z . Ces valeurs bruitées z sont fonction des signaux émis X t et des signaux reçus y . Elles sont représentatives de la transmission des signaux à travers les canaux de communication. En revanche, ces valeurs bruitées z sont calculées en s’affranchissant des signaux pilotes. Elles s’expriment par exemple selon la formule : avec la notation T correspondant à l’opérateur de transposition de matrices. Ces valeurs bruitées z vont servir de base pour la détermination de la valeur t associée à l’angle de départ θ tx . Afin d’estimer le canal de communication, l’unité de commande 2 initialise à la valeur 0 un indice i (étape E8). Cet indice i désigne l’itération courante. Lors de cette étape, l’unité de commande 2 initialise également la valeur d’une variable dite « résidu » r es i (étape E8). Cette variable res i est ici initialisée à partir des valeurs bruitées déterminées à l’étape E6 : res 0 = z . Chaque itération permet de rechercher une valeur prise par ce résidu res i . Comme décrit précédemment, la propagation des signaux dans les canaux de communication est caractérisée par la variable t associée à l’angle de départ θ tx , la variable r associée à l’angle d’arrivée θ r x et le retard de propagation τ. Le procédé d’estimation de canal comprend alors des étapes d’estimation des valeurs de ces variables. Dans la suite, les étapes présentées s’attachent à décrire la détermination de la valeur prise par l’une de ces variables (ici la variable t associée à l’angle de départ θ tx ). Les différentes étapes s’appliquent de la même manière pour la détermination des valeurs des autres variables caractérisant la propagation des signaux dans les canaux de communication. Pour caractériser la propagation des signaux dans le canal de communication, il convient tout d’abord de détecter un signal dans le système de communication. Le procédé d’estimation de canal comprend alors un procédé de détection Det d’un signal dans le système de communication. La figure 4 est un logigramme représentant un exemple de procédé de détection tel qu’il peut être mis en œuvre dans le cadre du procédé d’estimation de canal selon l’invention. Comme le montre la figure 4, à l’étape E10, le procédé de détection comprend une étape de découpage de l’intervalle des valeurs pouvant être prises par la valeur t en une pluralité de sous-intervalles, par exemple K t sous-intervalles j. De par sa définition, la valeur t associée à l’angle de départ θx est comprise entre -1 et 1. A l’étape E10, l’intervalle [-1, 1] est donc découpé en une succession de sous-intervalles de manière à couvrir l’ensemble de l’intervalle [-1, 1]. Le nombre Kt est par exemple prédéterminé en amont de l’exécution du procédé de détection, par exemple en fonction d’un niveau de performance recherché. La suite de procédé de détection consiste alors à tester la présence d’un signal sur chacun des sous-intervalles constituant l’intervalle des valeurs pouvant être prises par la valeur t associée à l’angle de départ θ tx . Pour cela, à l’étape E12, le processeur du dispositif du détection détermine, pour chaque sous-intervalle j, une fonction dite « corrélateur » f( res i , j). Ce corrélateur dépend des valeurs bruitées z . Plus particulièrement, chaque corrélateur associé au sous-intervalle j concerné dépend de la valeur courante du résidu res i (qui est égal aux valeurs bruitées z lors de la première itération, où i=0). Il dépend également du vecteur de dépointage e t associé à la direction correspondant au centre t j du sous-intervalle j. En pratique, ce corrélateur peut être interprété comme un filtre spatial associé au sous-intervalle concerné, permettant de filtrer les signaux afin par exemple de les distinguer du bruit sur ce sous-intervalle j. Selon un mode de réalisation, le corrélateur est défini par l’expression suivante : Avec t j le centre du sous-intervalle concerné, e t le vecteur de dépointage associé à la direction correspondant au centre t j du sous-intervalle concerné, res i la valeur courante du résidu. Selon un autre mode de réalisation, sur chaque sous-intervalle, chaque corrélateur est défini comme une somme de fonction de corrélateurs fenêtrés associé au sous-intervalle concerné. Dans ce cas, il est défini par l’expression suivante : avec e t le vecteur de dépointage associé à la direction t, t j le centre du sous-intervalle concerné, res i la valeur courante du résidu, λ k des scalaires et r tx, k des vecteurs et la notation ⊙ un opérateur symbolisant le produit terme à terme entre les différents éléments des vecteurs concernés (également appelé produit de Hadamard). Les scalaires λ k et les vecteurs r tx, k correspondent respectivement les valeurs propres et vecteurs propres associés à la matrice R 0 définie selon l’expression : avec Δ t la largeur du sous-intervalle concerné, || a tx || la distance entre deux antennes d’émission, λ la longueur d’onde des signaux émis et la notation sinc correspondant à la fonction sinus cardinal définie par sinc(x)=sin(x)/x. Les valeurs propres λ k et les vecteurs propres r tx, k ne dépendent que de la largeur Δ t du sous-intervalle j concerné. Les valeurs propres λ k décroissent vers 0. L’expression est appelée dans cette description « fonction de corrélateurs fenêtrés ». Comme le montre la figure 4, le procédé de détection se poursuit ensuite avec l’étape E14. Préalablement, il est possible de noter ici que, dans les deux modes de réalisation décrits précédemment, les corrélateurs observent un maximum lorsque le résidu courant res i est colinéaire au vecteur de dépointage e t associé à la direction correspondant au centre t j du sous-intervalle j concerné. Lors de cette étape, le processeur du dispositif de détection identifie le sous-intervalle correspondant au maximum de vraisemblance du corrélateur déterminé à l’étape E12. Plus précisément, l’unité de commande 2 détermine la valeur maximale atteinte par les corrélateurs parmi les corrélateurs déterminés pour chaque sous-intervalle j. Une fois la valeur maximale des corrélateurs déterminé, le processeur du dispositif de détection compare cette valeur maximale déterminée à un seuil prédéterminé. Ce seuil prédéterminé est fonction d’un niveau de bruit associé à l’intervalle étudié. Il s’agit par exemple ici d’un bruit gaussien réparti dans toutes les directions. Selon le premier mode de réalisation décrit précédemment, le processeur du dispositif de détection compare donc ici au seuil prédéterminé la valeur maximale parmi chacune des valeurs du corrélateur calculée selon le premier mode de réalisation introduit précédemment. En variante, selon le second mode de réalisation, le processeur du dispositif de détection compare, au seuil prédéterminé, la somme maximale de fonction de corrélateurs fenêtrés déterminée parmi les différentes sommes de fonctions de corrélateurs déterminées pour chaque sous-intervalle j. En pratique ici, lors de cette étape, l’unité de commande 2 identifie le sous-intervalle sur lequel la puissance du signal émis est la plus importante et telle que ce signal ne peut pas être considéré comme du bruit. Si à l’étape E14, le maximum de vraisemblance du corrélateur (soit la valeur maximale parmi les corrélateurs selon le premier mode de réalisation, soit la somme maximale parmi les sommes de fonction de corrélateurs fenêtrés déterminées selon le deuxième mode de réalisation) est supérieur au seuil prédéterminé, le processeur du dispositif de détection identifie le sous-intervalle contenant la valeur t de l’angle de départ recherchée (étape E16). On considère donc ici qu’un signal a été détecté dans le système de communication. Comme le montre la figure 3, le procédé d’estimation de canal se poursuit ensuite à l’étape E20. Lors de cette étape, l’unité de commande 2 reçoit du dispositif de détection l’information quant à la détection d’un signal dans le système de communication. Dans le cas où un signal a été détecté, le procédé d’estimation de canal comprend alors un procédé d’estimation Est (décrit ci-dessous et représenté sur la figure 5) de la valeur associée à la variable caractérisant la propagation des signaux dans le canal de communication ; plus particulièrement ici, le procédé d’estimation concerne l’estimation de la valeur t associée à l’angle de départ θ tx concernant le signal détecté. Si à l’étape E14, le maximum de vraisemblance du corrélateur est inférieur au seuil prédéterminé, on considère qu’aucun signal n’a été détecté. A l’étape E20, l’unité de commande 2 reçoit alors l’information qu’aucun signal n’a été détecté par le dispositif de détection. Le procédé d’estimation de canal se poursuit à l’étape E22. Cette absence de détection forme ici une condition d’arrêt. Lorsque la condition d’arrêt a été obtenue, le canal de communication est caractérisé sur la base de l’ensemble des signaux détectés précédemment à l’indice courant i. Le vecteur h est alors défini par l’expression suivante : Avec β p ,i l’estimation du gain des signaux électromagnétiques associés à l’estimation τ p de la valeur t obtenue lors des itérations précédentes du procédé (précédant l’indice courant i). On a donc ici P = i-1. Chaque paramètre β p ,i a donc été obtenu lors des itérations précédentes. En variante, une autre condition d’arrêt peut être définie, par exemple en déterminant la norme du résidu et en identifiant lorsque celle-ci est inférieure à un seuil prédéfini. La figure 5 est un logigramme représentant un exemple de procédé de d’estimation tel qu’il peut être mis en œuvre dans le cadre du procédé d’estimation de canal. Dans le cas où un signal a été détecté dans le système de communication (et donc un sous-intervalle a été identifié), le procédé d’estimation Est débute à l’étape E30. Lors de cette étape, le processeur du dispositif d’estimation détermine, à partir du sous-intervalle identifié à l’étape E14 du procédé de détection, une première estimation t 0 de la valeur t de l’angle de départ θ tx . Selon un mode de réalisation, cette première estimation t 0 correspond au centre t j du sous-intervalle identifié. Un autre mode de réalisation se base sur les valeurs des corrélateurs fenêtrés déterminés sur le sous-intervalle identifié. Plus particulièrement, la première estimation t 0 de la valeur t de l’angle de départ θ tx se base sur une fonction de test f t déterminée sur chaque sous-intervalle j et dépendant des corrélateurs fenêtrés déterminés précédemment : La détermination de la première estimation t 0 se base alors sur une comparaison des corrélateurs fenêtrés associés au sous-intervalle j identifié à l’étape E14, correspondant donc au signal détecté lors du procédé de détection Det. La figure 6 représente l’évolution de cette fonction test f t en fonction des différentes valeurs possibles pour la valeur t de l’angle de départ θ tx . La courbe f correspond à cette fonction de test en tenant compte de l’ensemble des corrélateurs fenêtrés. La courbe a correspond à la contribution du premier corrélateur fenêtré à la fonction f t . La courbe b correspond à la contribution du deuxième corrélateur fenêtré à la fonction f t . Cette figure permet alors d’évaluer la contribution de chaque corrélateur fenêtré dans la fonction de test f t (res i , j). L’étude de chacune de ces contributions de chaque corrélateur fenêtré permet alors l’évaluation de la première estimation t 0 de la valeur t de l’angle de départ θ tx . En particulier, la différence de phase associée à chaque corrélateur fenêtré permet de localiser la partie du sous-intervalle contenant la valeur t recherchée (par exemple à gauche ou à droite du centre du sous-intervalle). L’amplitude de chaque corrélateur fenêtré permet d’identifier la position de la valeur t recherchée dans le sous-intervalle concerné. En pratique, cette première estimation t 0 de la valeur t associée à l’angle de départ θ tx correspond à une estimation approchée de cette valeur. A ce stade du procédé d’estimation, on dispose d’une estimation approchée de la valeur t associée à l’angle de départ θ tx (on sait également, grâce à la mise en œuvre du procédé de détection dans quel sous-intervalle se trouve cette valeur t). Les étapes suivantes du procédé d’estimation ont alors pour but d’affiner cette première estimation t 0 . Une solution classique consiste à utiliser des méthodes d’optimisation telles que la méthode de Newton-Raphson ou la descente du gradient. Cependant, ces méthodes reposent sur l’utilisation de fonctions présentant des propriétés de convexité. Ici, la fonction étudiée (nommée corrélateur précédemment) définie par présente un comportement ressemblant à une fonction Dirac au voisinage de son maximum (Figure 7). Les méthodes d’optimisation classiques ne sont pas adaptées pour traiter ce type de fonction. Il est donc nécessaire dans un premier temps de modifier le corrélateur afin de pouvoir ensuite y appliquer les méthodes d’optimisation classiques. Le procédé d’estimation se poursuit donc à l’étape E32 lors de laquelle le processeur du dispositif d’estimation détermine une fonction modifiée f mod _i . Cette fonction modifiée f mod _i est déterminée sur la base d’un produit scalaire entre le vecteur associé aux valeurs bruitées et le vecteur de dépointage. Plus précisément ici, le vecteur associé aux valeurs bruitées z correspond ici au vecteur associé à la valeur courante du résidu res i . La fonction modifiée f mod _i présente des propriétés de convexité, c’est-à-dire que cette fonction modifiée f mod _i est soit convexe, soit concave. Ces propriétés vont donc permettre de mettre en œuvre les méthodes d’optimisation classiques. De plus, cette fonction modifiée f mod _i atteint une valeur maximale pour la valeur t correspondant au maximum du corrélateur initial (c’est-à-dire la forme non modifiée). En pratique ici, la détermination de cette fonction modifiée f mod _i revient à déterminer le résultat d’une opération de convolution entre le corrélateur | e t T (t). res i |² et une fonction f n dite « noyau de convolution ». Ce noyau de convolution f n présente par exemple ici la forme d’une section de parabole (Figure 8). La fonction modifiée f mod _i , à laquelle va être appliquée la méthode d’optimisation, dépend ici d’une somme de corrélateurs fenêtrés (tels qu’introduits lors du procédé de détection) : avec λ k des scalaires et r tx, k des vecteurs dépendant du noyau de convolution f n choisi. La figure 9 représente la variation de cette fonction modifiée f mod _i en fonction des différentes valeurs possibles de la variable t. La fonction modifiée f mod _i a ici la forme d’une parabole centrée sur la valeur t associée à l’angle de départ recherché. Une méthode d’optimisation classique est donc appliquée à cette fonction modifiée f mod _i . On applique ici la méthode de Newton-Raphson, qui se décompose selon les étapes suivantes E34 à E44. Le but de cette méthode est de déterminer la position du maximum de la fonction modifiée f mod _i . Comme le montre la figure 6, le procédé d’estimation Est se poursuit à l’étape E34 d’initialisation à la valeur 0 d’un indice l. Cet indice l désigne le tour courant de cette méthode d’optimisation. Lors de cette étape, l’unité de commande 2 initialise également la valeur d’une variable t l . On utilise ici comme valeur d’initialisation la première estimation t 0 déterminée à l’étape E30. A l’étape E36, le processeur du dispositif d’estimation détermine, pour le tour courant, les valeurs des dérivées première et seconde de la fonction modifiée f mod _i en la valeur de la variable courante t l . En d’autres termes, en utilisant les notations classiques, le processeur du dispositif d’estimation détermine les valeurs f mod _i ’(z, t l ) et f mod _i ’’(z, t l ). Puis, à l’étape E38, le processeur du dispositif d’estimation détermine la valeur de la variable t l +1 définie par l’expression suivante : Le procédé se poursuit à l’étape E40 lors de laquelle le processeur du dispositif d’estimation évalue si la valeur t l +1 déterminée correspond au maximum de la fonction modifiée f mod _i . Pour cela, le processeur calcule la quantité |t l +1 -t l |. Si cette quantité est supérieure à une valeur prédéterminée ε (|t l +1 -t l |>ε), il ne peut pas être considéré que la convergence vers le maximum est atteinte. Le procédé se poursuit alors à l’étape E42 lorsque de laquelle la valeur t l est actualisée par la valeur t l +1 déterminée à l’étape E38. L’indice l est également incrémenté. Une nouvelle itération est alors mise en œuvre et le procédé reprend à l’étape E36. En revanche, si, à l’étape E40, la quantité |t l +1 -t l | est inférieure à la valeur prédéterminée ε (|t l +1 -t l | l +1 déterminée peut être considérée comme représentant la valeur du maximum de la fonction modifiée f mod _i (étape E44). En d’autres termes, cette valeur t l +1 correspond à une estimation plus fine de la valeur t associée à l’angle de départ θ t x recherché. Dans l’exemple d’un noyau de convolution parabolique, une seule itération est suffisante pour atteindre la valeur du maximum avec la méthode de Newton-Raphson. En variante, d’autres méthodes d’optimisation peuvent être utilisées telle que par exemple la descente de gradient. On note ici τ p la valeur associée à l’angle de départ obtenue à cette étape. Une fois l’estimation de la valeur τ p associée à l’angle de départ obtenue, le procédé d’estimation de canal se poursuit à l’étape E50 (Figure 3). A partir de l’estimation de la valeur t associée à la direction de départ obtenue, le paramètre β p ,i représentant l’estimation du gain des signaux électromagnétiques associés à l’estimation τ p de la valeur t obtenue par l’estimation issue de l’itération courante i du procédé est déterminé. Plus précisément, leur détermination se fonde sur le calcul d’un pseudo-inverse. En effet, en utilisant une notation matricielle, le canal de propagation h s’écrit : h = E t . b avec E t la matrice contenant les vecteurs de dépointage e t et b la matrice associée aux estimations de gain β p ,i . La détermination de la matrice b associée aux estimations de gain β p (et donc des estimations de gain β p ,i elles-mêmes) utilise alors une pseudo inversion selon la formule suivante : b =( E t H .E t ) -1 . E t H h . Ce paramètre β p ,i est ici recalculé à chaque itération du procédé pour toutes le valeurs p inférieures ou égales à l’indice courant i. Le signal détecté est alors entièrement caractérisé et le procédé se poursuit en actualisant la valeur du résidu r es i précédemment introduit (étape E52) pour tenir compte du dernier signal détecté et de l’estimation associée à l’angle de départ correspondant. En d’autres termes, le signal, parmi les signaux restant dans le résidu, dont la puissance était la plus élevée est retranché des valeurs bruitées z pour obtenir un nouveau réside res i+1 : Le nouveau résidu res i+1 est recalculé complètement à partir des valeurs bruitées z à chaque itération car les estimations des gains β p ,i sont actualisées à chaque itération. L’indice i est alors incrémenté à l’étape E54 et le procédé reprend ensuite en amont du procédé de détection Det, tant que la condition d’arrêt définie à l’étape E20 n’est pas obtenue. Ces valeurs estimées sont utilisées, dans un contexte de démodulation, afin de retirer l’influence du canal de propagation dans les données. Ces valeurs estimées peuvent aussi être utilisées par l’unité de commande 2 pour configurer des circuits de traitement des signaux électromagnétiques reçus par les antennes R 1 , …, R N r du réseau d’antennes (ces circuits de traitement étant inclus dans l’unité de commande 2 mais non représenté pour simplifier l’exposé). Ces valeurs estimées peuvent également être utilisées pour configurer des précodeurs conçus pour effectuer un précodage des signaux électromagnétiques à émettre au moyen des antennes R 1 , …, R N r du réseau d’antennes (lorsque ces antennes fonctionnent également en émission comme mentionné plus haut). Annexe : Démonstration de la formule « Math. 6 » La présente invention comporte une étape de détection de rayon dont les paramètres (ici la variable t associée à la direction de départ θ tx ) sont inconnus, basée sur les échantillons reçus : avec h=β.e t * (t). Alors la distribution de y sachant la direction de départ t est donnée par : La direction de départ appartient à l’intervalle [-1, 1]. Cet intervalle est découpé en une succession de sous-intervalles. On teste la probabilité de présence du rayon sur chaque sous-intervalle. La probabilité de présence du rayon dans l’intervalle est donnée par : De façon générale, on va calculer Cette probabilité peut être réécrite sous la forme : On effectue un développement limité de l’exponentielle au premier ordre au voisinage de 0 : e x ≈1+x. La probabilité se réécrit maintenant en introduisant des constantes K 0 et K 1 comme : On pose : Alors l’intégrale peut se réécrire sous la forme : Où |a| représente la distance entre deux antennes et λ représente la longueur d’onde associée à la fréquence de la porteuse utilisée. On remarque que l’intégrale fait intervenir la transformée de Fourier de la fonction noyau g(t), appelée La probabilité de présence de rayon peut se réécrire : Où et Dans le cas présenté ici où on a : Pour chaque sous-intervalle j, de taille Δ t et centré en t j , on pose : la partie dépendant de z de la probabilité de présence du rayon dans l’intervalle j. La matrice R 0 peut être décomposée en vecteurs propres et valeurs propres de la façon suivante : ce qui permet d’écrire la fonction f i ( z , j) comme Lors de la phase de détection, on calcule la probabilité de présence de rayon P j pour chaque intervalle j : On recherche ensuite le sous-intervalle le plus susceptible de contenir le rayon. Cette opération revient à retenir le sous-intervalle qui a relevé la probabilité maximale. Pour cela, le calcul complet de P j n’est pas nécessaire, il est suffisant de calculer f i ( z , j), la partie variable de P j . La valeur atteinte par f( z , j) sur le sous-intervalle retenu est testée par rapport à un seuil. Le résultat du test rend la décision du procédé de détection. Procédé de détection (Det) d’un signal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, la propagation dans lesdits canaux de communication étant caractérisée par au moins une variable, l’ensemble des valeurs pouvant être prises par ladite variable étant découpé en une pluralité d’intervalles, ledit procédé comprenant les étapes de : - réception (E4) des signaux respectivement transmis dans la pluralité des canaux de communication, et - détection (E20) d’un signal correspondant à une valeur de ladite variable comprise dans un intervalle parmi la pluralité d’intervalles par comparaison, à un seuil prédéterminé, d’une valeur prise par un corrélateur lié à l’intervalle concerné et calculée en fonction des valeurs bruitées. Procédé de détection (Det) selon la revendication 1, dans lequel l’étape de détection comprend une étape de comparaison (E14), audit seuil prédéterminé, d’une valeur maximale parmi une pluralité de valeurs prises par le corrélateur. Procédé de détection (Det) selon la revendication 1, dans lequel l’étape de détection (E20) comprend une étape de détermination (E12), pour chaque intervalle, d’une somme de fonction de corrélateurs fenêtrés associés à l’intervalle concerné et fonction desdites valeurs bruitées. Procédé de détection (Det) selon la revendication 3, dans lequel l’étape de détection (E20) comprend une étape de comparaison (E14) de la somme de fonction de corrélateurs fenêtrés audit seuil prédéterminé. Procédé de détection (Det) selon l’une quelconque des revendications 1 à 4, dans lequel le seuil prédéterminé est fonction d’un niveau de bruit associé à l’intervalle concerné. Procédé d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit procédé comprenant les étapes de : - détection (Det) d’un signal par la mise en œuvre d’un procédé de détection (Det) selon l’une quelconque des revendications 1 à 5, et - estimation (E44) de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. Procédé d’estimation de canal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit procédé comprenant, tant qu’une condition d’arrêt n’est pas obtenue, la répétition d’étapes de : - détection (E20) d’un signal par la mise en œuvre d’un procédé de détection selon l’une quelconque des revendications 1 à 5, et - estimation (E44) de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté, ladite condition d’arrêt correspondant à l’absence de détection d’un signal lors de la mise en œuvre du procédé de détection (Det) selon l’une quelconque des revendications 1 à 5. Procédé d’estimation selon la revendication 6 ou 7, dans lequel l’étape d’estimation (Est) comprend une étape de détermination (E32) d’une valeur intermédiaire de ladite variable comme égale au centre de l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. Procédé d’estimation selon la revendication 6 ou 7, dans lequel l’étape d’estimation (Est) comprend une étape de détermination (E32) d’une valeur intermédiaire de ladite variable par comparaison de corrélateurs fenêtrés correspondant à l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. Procédé d’estimation selon l’une quelconque des revendications 7 à 9, dans lequel l’étape d’estimation (Est) comprend en outre une étape de détermination (E34) d’une fonction (f mod ) ayant une propriété de convexité associée à un produit scalaire d’un vecteur formé des valeurs bruitées et d’un vecteur de dépointage dépendant de ladite variable, le maximum de ladite fonction (f mod ) ayant une propriété de convexité et dudit produit scalaire étant obtenu pour la même valeur de ladite variable. Procédé d’estimation selon la revendication 10 prise dans la dépendance de la revendication 8 ou 9, dans lequel ladite valeur de la variable est déterminée par une étape d’optimisation sur la base de ladite valeur intermédiaire de la variable et utilisant ladite fonction (f mod ) ayant une propriété de convexité. Dispositif de détection (5) d’un signal dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, la propagation dans lesdits canaux de communication étant caractérisée par au moins une variable, l’ensemble des valeurs pouvant être prises par ladite variable étant découpé en une pluralité d’intervalles, ledit dispositif de détection (5) comprenant : - un module de réception des signaux respectivement transmis dans la pluralité des canaux de communication, et - un module de détection d’un signal correspondant à une valeur de ladite variable comprise dans un intervalle parmi la pluralité d’intervalles par comparaison, à un seuil prédéterminé, d’une valeur prise par un corrélateur lié à l’intervalle concerné et calculée en fonction des valeurs bruitées. Dispositif d’estimation de canal (1) dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit dispositif d’estimation de canal comprenant : - un dispositif de détection (5) d’un signal selon la revendication 12, et - un module d’estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté. Dispositif d’estimation de canal (1) dans un système de communication comprenant une pluralité de canaux de communication, à partir d’une pluralité de valeurs bruitées respectivement représentatives de la transmission à travers lesdits canaux de communication, ledit dispositif d’estimation de canal étant conçu pour activer, tant qu’une condition d’arrêt n’est pas obtenue : - un dispositif de détection (5) d’un signal par la mise en œuvre d’un procédé de détection selon la revendication 12, et - un module d’estimation de la valeur de la variable dans l’intervalle comprenant la valeur de ladite variable correspondant au signal détecté, ladite condition d’arrêt correspondant à l’absence de détection d’un signal par le dispositif de détection selon la revendication 12. Programme d’ordinateur comprenant des instructions exécutables par un processeur et conçues pour mettre en œuvre un procédé selon l’une des revendications 1 à 11 lorsque ces instructions sont exécutées par le processeur.