FR 3087974 A1 20200501 FR 1859952 A 20181026 Procédé d'estimation d'un canal de communication sans fil, dispositif d'estimation d'un canal de communication sans fil et programme d'ordinateur associés 1. Domaine de l'invention Le domaine de l'invention est celui des systèmes de communication sans fil, en particulier les systèmes massivement multi-antennes dits « Massive MIMO » (pour « Massive Multiple Input Multiple Output » en anglais). 2. Présentation de l'art antérieur L'invention se place dans le cadre de l'étude des réseaux mobiles en particulier ceux de la Se génération (5G). Les débits de données envisagés pour ces réseaux dits 5G sont très élevés. Pour les atteindre, les systèmes mufti-antennes (MIMO) sont prometteurs, car leur débit théorique maximal évolue proportionnellement au nombre d'antennes utilisées, à puissance d'émission fixée. Cependant, atteindre un débit théorique maximal peut être problématique. Il faut pour cela mettre en oeuvre des traitements bien plus complexes que pour des systèmes à une seule antenne. Il est en outre préférable de connaître le canal de transmission (encore appelé canal de communication) qui sépare l'émetteur du récepteur. Pour cela on réalise une estimation du canal de communication. Pour un système avec Nt antennes à l'émission et N, antennes à la réception, l'estimation de canal MIMO consiste à déterminer un coefficient complexe par fréquence pour chacun des N'IV, couples d'antennes du système. Dans le cas où N, et/ou N, est très élevé (on parle alors de système massivement multi-antennes : « Massive MIMO »), il est très coûteux d'estimer chacun des coefficients de manière indépendante. Il est alors possible d'utiliser un modèle physique qui permet de prendre en compte des relations de liaison entre les coefficients et ainsi réduire le nombre de paramètres à estimer, et par conséquent le coût de l'estimation. Les modèles physiques utilisés dans l'état de l'art, postulent que le récepteur est suffisamment éloigné de l'émetteur par rapport à leurs tailles respectives pour que les fronts d'ondes sphériques soient bien approchés par des plans : c'est l'hypothèse ondes planes. Cette hypothèse est faite par toutes les méthodes d'estimation basées sur un modèle physique, voir par exemple l'article fondateur « Bajwa et al.: Compressed Channel Sensing: A New Approach to Estimating Sparse Multipath Channels » Vol. 98, No. 6, June 2010. Les systèmes massivement multi-antennes pour les réseaux mobiles de 5e génération sont amenés à avoir de plus en plus d'antennes. Les émetteurs/récepteurs seront de fait de plus en plus grands. Cela implique l'hypothèse ondes planes ne sera vérifiée que pour des distances de plus en plus grandes. 5 Ceci met en défaut les méthodes d'estimation actuelles, en particulier dans le cas de très grands réseaux d'antennes et/ou de récepteurs proches de l'émetteur. L'invention a notamment pour objectif de pallier à ces inconvénients de l'art antérieur. 3. Exposé de l'invention 10 Selon un premier aspect de l'invention, il est proposé un procédé d'estimation d'un canal de communication sans fil entre un émetteur et un récepteur, comprenant une pluralité de trajets permettant la propagation d'une onde, au moins l'émetteur ou le récepteur étant formé de plusieurs antennes. Selon une caractéristique générale de l'invention, le procédé comprend : 15 pour au moins un trajet, une détermination d'une matrice dite caractéristique, qui dépend d'un premier élément représentatif d'au moins une direction de propagation associée au trajet, et d'un deuxième élément représentatif d'une distance de propagation associée audit trajet, et une réalisation d'une estimation du canal de communication à partir de la au moins une 20 matrice caractéristique obtenue. En d'autres termes, les inventeurs proposent une nouvelle façon d'estimer un canal de communication sans fil, à partir d'une matrice dite caractéristique qui prend en compte non seulement une direction de propagation, mais également une distance de propagation. Cette matrice caractéristique permet une meilleure estimation du canal de communication, notamment lorsque la 25 configuration émetteur/récepteur (par exemple liée à la taille des antennes et/ou de la distance entre l'émetteur et le récepteur) ne permet pas de vérifier l'hypothèse ondes planes. De préférence, selon un mode de mise en oeuvre, la matrice caractéristique associée au trajet dépend en outre d'au moins l'un des deux éléments suivants : au moins un vecteur reliant le centre de gravité du réseau d'antennes de l'émetteur 30 et au moins une antenne de l'émetteur, au moins un vecteur reliant centre de gravité du réseau d'antennes du récepteur et au moins une antenne du récepteur. Ce mode de mise en oeuvre permet de prendre en compte toutes les configurations de réseaux d'antennes à l'émission et à la réception, quel que soit le nombre d'antennes et leurs 5 positions dans l'espace. Selon un mode de mise en oeuvre, la définition d'une matrice caractéristique dépend d'un modèle de canal parmi une pluralité de modèles. La matrice caractéristique est ainsi adaptable au modèle de canal choisi, selon l'hypothèse de propagation des ondes qui est faite. 10 Selon un mode de mise en oeuvre, la pluralité de modèles peut comprendre le modèle ondes sphériques et un modèle ondes dites paraboliques décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne du récepteur d'une part, et la distance existant entre les centres de gravité des réseaux d'antennes du récepteur et de l'émetteur d'autre part, 15 le modèle d'ondes dites paraboliques étant fonction uniquement : d'une grandeur dépendant d'au moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, un vecteur de direction d'arrivée de propagation associée au trajet, et d'un terme correctif fonction de la distance de propagation associée audit trajet, 20 ladite grandeur et ledit terme correctif correspondant respectivement audit premier et audit deuxième élément, mentionnés ci-avant. Le modèle d'ondes paraboliques a pour avantage de permettre une meilleure estimation du canal lorsque l'hypothèse ondes planes n'est plus valide, mais sans la complexité liée à l'hypothèse ondes sphériques, qui est particulièrement couteuse en temps de calcul. 25 De préférence, selon une première alternative d'un mode de mise en oeuvre de l'invention, le procédé inclut en outre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, qui comprend : une obtention de la distance entre les centres de gravité des réseaux d'antennes du récepteur d'une part et de l'émetteur d'autre part, 30 une comparaison de ladite distance à un premier et un second seuil, une validation d'un modèle ondes planes si la distance obtenue est supérieure au second seuil, du modèle ondes dites paraboliques si la distance obtenue est comprise entre le premier et le second seuil, et le modèle ondes sphériques si la distance obtenue est inférieure au premier seuil. Cette étape préalable permet de valider la pertinence du modèle de canal choisi pour une 5 estimation de canal qui sera réalisée ultérieurement, notamment pour le contexte émetteur/récepteur considéré. Avantageusement, le modèle validé lors de l'étape préalable est optimal pour l'estimation de canal à réaliser. De préférence : le premier seuil est fonction de la distance à la puissance 3/2 entre le centre de l'ensemble 10 des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur, et le second seuil est fonction de la distance au carré entre le centre de l'antenne de l'émetteur la plus éloignée du centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur. Le choix de ce premier et de ce deuxième seuil est particulièrement simple à implémenter. 15 De préférence, selon une deuxième alternative d'un mode de mise en oeuvre de l'invention, le procédé inclut en outre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, qui comprend une détermination d'une erreur relative d'approximation du modèle pour au moins un premier et un deuxième modèle, chaque erreur relative d'approximation du modèle étant fonction d'un écart 20 entre un modèle de référence et l'un desdits deux modèles une comparaison de chaque erreur relative d'approximation du modèle à une valeur prédéterminée de marge d'erreur, et une validation ou non du premier modèle de canal ou du deuxième modèle de canal en fonction du résultat de la comparaison. 25 Cette deuxième alternative permet de s'affranchir d'un écart de déphasage arbitraire dont dépend la première alternative présentée ci-avant. Elle permet également de mieux tenir compte des positions relatives et de l'orientation des antennes de l'émetteur et/ou du récepteur. Selon une première alternative d'un mode de mise en oeuvre de l'invention, la réalisation d'une estimation du canal de communication comprend de préférence, au moins pour un premier 30 trajet : 4 Une première estimation d'une première valeur de distance de propagation associée au premier trajet, et d'une première direction de propagation associée au premier trajet, et une détermination d'une première matrice caractéristique dépendant de ladite première valeur de distance de propagation et de ladite première direction de propagation estimées. 5 Ce mode de mise en oeuvre a pour avantage de proposer une estimation précise du canal. De préférence, la première alternative peut comprendre en outre pour un deuxième trajet : Une deuxième estimation d'une deuxième valeur de distance de propagation associée audit deuxième trajet, et d'une deuxième direction de propagation associée audit deuxième trajet, et une détermination d'une deuxième matrice caractéristique dépendant de ladite deuxième valeur 10 de distance de propagation et de ladite deuxième direction de propagation estimées, ainsi que d'un résidu résultant de la première estimation. Avantageusement, ce mode de mise en oeuvre permet de réaliser une estimation de canal dans le cas où plusieurs trajets existent, via une séquence d'estimations simples ne considérant qu'un unique trajet (algorithme dit « glouton »). 15 Selon une deuxième alternative de ce mode de mise en oeuvre de l'invention, au moins la première estimation ou la deuxième estimation peut être effectuée en une première sous-estimation, pour une première grandeur fixée et une deuxième grandeur, distincte de la première grandeur, pouvant prendre une pluralité de valeurs prédéterminées, puis 20 en une deuxième sous-estimation pour la deuxième grandeur fixée à une valeur résultant de la première sous-estimation et la première valeur pouvant prendre une pluralité de valeurs prédéterminées. Cette deuxième alternative propose de réaliser l'estimation de canal de manière séquentielle. Celle-ci est beaucoup moins couteuse en termes de temps de calcul que la première alternative 25 présentée ci-dessus. Par exemple : la première grandeur ou la deuxième grandeur est l'un des éléments suivants : un élément représentatif de la distance de propagation associée au trajet considéré, ou un élément représentatif de la direction de propagation associée au trajet considéré. 30 Choisir comme première grandeur un élément représentatif de la direction est particulièrement avantageux dans le cas d'une transmission directive. A l'inverse, choisir comme 6 première grandeur un élément représentatif de la distance de propagation est simple à mettre en oeuvre. Selon un deuxième aspect de l'invention, il est proposé un procédé de précodage d'un signal à transmettre via un canal de communication, comprenant : 5 une estimation du canal de communication selon l'une quelconques des revendications précédentes afin d'obtenir un canal de communication estimé, et un précodage du signal à transmettre en fonction du canal de communication estimé. Un précodage adapté permet de maximiser le débit de la communication. Le précodage 10 sera d'autant mieux adapté que l'estimation du canal est précise (or la prise en compte de la courbure des fronts d'onde améliore la précision de l'estimation de canal). Selon un troisième aspect de l'invention, il est proposé un procédé de démodulation d'un signal reçu après propagation via un canal de communication, comprenant : une estimation du canal de communication obtenue par l'un des procédés 15 d'estimation d'un canal de communication décrits ci-dessus afin d'obtenir un canal de communication estimé, et une démodulation du signal reçu en fonction du canal de communication estimé. Une démodulation adaptée permet de maximiser le débit de la communication. La démodulation sera d'autant mieux adaptée que l'estimation du canal est précise (or la prise en 20 compte de la courbure des fronts d'onde améliore la précision de l'estimation de canal). Selon un quatrième aspect de l'invention, il est proposé un procédé de détermination d'un modèle ondes dites paraboliques, en particulier implémenté par un ordinateur, ledit modèle d'ondes parabolique décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne du récepteur d'une part, et la 25 distance existant entre les centres de gravité des réseaux d'antennes du récepteur et de l'émetteur d'autre part, le modèle d'ondes dites paraboliques étant fonction uniquement : d'une grandeur dépendant tau moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, 30 un vecteur de direction d'arrivée de propagation associée au trajet, 7 et d'un terme correctif fonction de la distance de propagation associée audit trajet. Le modèle d'ondes paraboliques a pour avantage de permettre une meilleure estimation du canal lorsque l'hypothèse ondes planes n'est plus valide, mais sans la complexité liée à l'hypothèse ondes sphériques, qui est particulièrement couteuse en temps de calcul. 5 Selon un cinquième aspect de l'invention, il est proposé un procédé d'estimation d'un canal de communication sans fil entre un émetteur et un récepteur, comprenant une pluralité de trajets permettant la propagation d'une onde, au moins l'émetteur ou le récepteur étant formé de plusieurs antennes. Selon une caractéristique générale de ce cinquième aspect de l'invention, le procédé comprend : 10 pour au moins un trajet, une détermination d'une matrice dite caractéristique, qui dépend du modèle d'ondes paraboliques, et une réalisation d'une estimation du canal de communication à partir de la au moins une matrice caractéristique obtenue. La matrice caractéristique fonction du modèle d'ondes paraboliques permet une meilleure estimation du canal de communication, notamment lorsque la 15 configuration émetteur/récepteur (par exemple liée à la taille des antennes et/ou de la distance entre l'émetteur et le récepteur) ne permet pas de vérifier l'hypothèse ondes planes. Selon un sixième aspect de l'invention, il est proposé un dispositif d'estimation d'un canal de communication sans fil entre un émetteur et un récepteur, comprenant une pluralité de trajets permettant la propagation d'une onde, au moins l'émetteur ou le récepteur étant formé de 20 plusieurs antennes. Selon une caractéristique générale de ce sixième aspect de l'invention, le dispositif comprend : un moyen de détermination configuré pour mettre en oeuvre, pour au moins un trajet, une détermination d'une matrice dite caractéristique, qui dépend d'un premier élément représentatif tau moins une direction de propagation associée au 25 trajet, et d'un deuxième élément représentatif d'une distance de propagation associée audit trajet, et un moyen de réalisation configuré pour mettre en oeuvre une réalisation d'une estimation du canal de communication à partir de la au moins une matrice caractéristique obtenue. De préférence, la matrice caractéristique associée au trajet dépend en outre d'au moins 30 l'un des deux éléments suivants : 8 au moins un vecteur reliant le centre de gravité du réseau d'antennes de l'émetteur et au moins une antenne de l'émetteur, au moins un vecteur reliant centre de gravité du réseau d'antennes du récepteur et au moins une antenne du récepteur. 5 De préférence, la définition d'une matrice caractéristique dépend d'un modèle de canal parmi une pluralité de modèles. Selon un mode de réalisation, la pluralité de modèles peut comprendre le modèle ondes sphériques et un modèle ondes dites paraboliques décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne 10 du récepteur d'une part, et la distance existant entre les centres de gravité des réseaux d'antennes du récepteur et de l'émetteur d'autre part, le modèle d'ondes dites paraboliques étant fonction uniquement : d'une grandeur dépendant d'au moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, 15 un vecteur de direction d'arrivée de propagation associée au trajet, et d'un terme correctif fonction de la distance de propagation associée audit trajet, ladite grandeur et ledit terme correctif correspondant respectivement audit premier et audit deuxième élément. Selon une première variante d'un mode de réalisation, le dispositif peut inclure en outre en 20 outre un moyen de validation configuré pour mettre en oeuvre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, ledit moyen de validation comprenant un moyen d'obtention configuré pour mettre en oeuvre une obtention de la distance entre les centres de gravité des réseaux d'antennes du récepteur d'une part et de l'émetteur d'autre part, 25 un moyen de comparaison configuré pour mettre en oeuvre une comparaison de ladite distance à un premier et un second seuil, un moyen de décision configuré pour mettre en oeuvre une validation d'un modèle ondes planes si la distance obtenue est supérieure au second seuil, du modèle ondes dites paraboliques si la distance obtenue est comprise entre le premier et le second seuil, et le modèle ondes sphériques 30 si la distance obtenue est inférieure au premier seuil. 9 De préférence, le premier seuil est fonction de la distance à la puissance 3/2 entre le centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur, et le second seuil est fonction de la distance au carré entre le centre de l'antenne de 5 l'émetteur la plus éloignée du centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur. Selon une deuxième variante du mode de réalisation, le dispositif peut inclure en outre un moyen de validation configuré pour mettre en oeuvre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, qui comprend 10 un moyen de détermination configuré pour mettre en oeuvre une détermination d'une erreur relative d'approximation du modèle pour au moins un premier et un deuxième modèle, chaque erreur relative d'approximation du modèle étant fonction d'un écart entre un modèle de référence et l'un desdits deux modèles un moyen de comparaison configuré pour mettre en oeuvre une comparaison de chaque 15 erreur relative d'approximation du modèle à une valeur prédéterminée de marge d'erreur, et un moyen de décision configuré pour mettre en oeuvre une validation ou non du premier modèle de canal ou du deuxième modèle de canal en fonction du résultat de la comparaison. Selon un mode de réalisation, ledit moyen de réalisation est configuré pour mettre en oeuvre une réalisation d'une estimation du canal de communication et comprend: 20 un moyen d'estimation configuré pour mettre en oeuvre, pour au moins pour un premier trajet, une première estimation d'une première valeur de distance de propagation associée au premier trajet, et d'une première direction de propagation associée au premier trajet, et une détermination d'une première matrice caractéristique dépendant de ladite première valeur de distance de propagation et de ladite première direction de propagation estimées. 25 De préférence, le dispositif peut comprendre en outre : un deuxième moyen d'estimation configuré pour mettre en oeuvre, pour un deuxième trajet, une deuxième estimation d'une deuxième valeur de distance de propagation associée audit deuxième trajet, et d'une deuxième direction de propagation associée audit deuxième trajet, et une détermination d'une deuxième matrice caractéristique dépendant de ladite 30 deuxième valeur de distance de propagation et de ladite deuxième direction de propagation estimées, ainsi que d'un résidu résultant de la première estimation. 10 De préférence, au moins le premier moyen d'estimation ou le deuxième moyen d'estimation peut être configuré pour effectuer une première sous-estimation, pour une première grandeur fixée et une deuxième grandeur, distincte de la première grandeur, pouvant prendre une pluralité de valeurs 5 prédéterminées, puis une deuxième sous-estimation pour la deuxième grandeur fixée à une valeur résultant de la première sous-estimation et la première valeur pouvant prendre une pluralité de valeurs prédéterminées. De préférence, 10 la première grandeur ou la deuxième grandeur est l'un des éléments suivants : un élément représentatif de la distance de propagation associée au trajet considéré, OU un élément représentatif de la direction de propagation associée au trajet considéré. 15 Selon un septième aspect de l'invention, il est proposé un système de précodage d'un signal à transmettre via un canal de communication, comprenant : un dispositif d'estimation d'un canal de communication tel que l'un de ceux décrits ci-dessus, afin d'obtenir un canal de communication estimé, ledit système de précodage étant configuré pour mettre en oeuvre un précodage du signal à 20 transmettre en fonction du canal de communication estimé. Selon un sixième de l'invention, il est proposé un système de démodulation d'un signal reçu après propagation via un canal de communication, comprenant : un dispositif d'estimation d'un canal de communication tel que l'un de ceux décrits ci-dessus afin d'obtenir un canal de communication estimé, 25 ledit système de démodulation étant configuré pour mettre en oeuvre une démodulation du signal reçu en fonction du canal de communication estimé. Selon un huitième aspect de l'invention, il est proposé un dispositif configuré pour déterminer un modèle d'ondes dites paraboliques décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne 30 du récepteur d'une part, et la distance existant entre les centres de gravité des réseaux d'antennes du récepteur et de l'émetteur d'autre part, 1l le modèle d'ondes dites paraboliques étant fonction uniquement : d'une grandeur dépendant d'au moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, un vecteur de direction d'arrivée de propagation associée au trajet, 5 et d'un terme correctif fonction de la distance de propagation associée audit trajet. Selon un neuvième aspect de l'invention, il est proposé un système d'estimation configuré pour estimer un canal de communication sans fil entre un émetteur et un récepteur, comprenant une pluralité de trajets permettant la propagation d'une onde, au moins l'émetteur ou le récepteur étant formé de plusieurs antennes. Selon une caractéristique générale de ce neuvième aspect de 10 l'invention, le système comprend : Un moyen de détermination afin de déterminer pour au moins un trajet, une matrice dite caractéristique, qui dépend du modèle d'ondes paraboliques décrit ci-avant, et un moyen de réalisation pour réaliser une estimation du canal de communication à partir de la au moins une matrice caractéristique obtenue. 15 La matrice caractéristique fonction du modèle d'ondes paraboliques permet une meilleure estimation du canal de communication, notamment lorsque la configuration émetteur/récepteur (par exemple liée à la taille des antennes et/ou de la distance entre l'émetteur et le récepteur) ne permet pas de vérifier l'hypothèse ondes planes. Selon un dixième aspect de l'invention, il est proposé un programme d'ordinateur 20 comportant des instructions pour la mise en oeuvre d'un procédé d'estimation d'un canal de communication tel que décrit ci-dessus. Ce programme peut être exécuté par un processeur. Selon un onzième aspect de l'invention, il est proposé un programme d'ordinateur comportant des instructions pour la mise en oeuvre d'un procédé de précodage d'un signal à transmettre via un canal de communication tel que décrit ci-dessus ou d'un procédé de 25 démodulation d'un signal reçu après propagation via un canal de communication tel que décrit ci- dessus. Ce programme peut être exécuté par un processeur. Ces programmes peuvent utiliser n'importe quel langage de programmation. Il peut être téléchargé depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur. 30 Selon un douzième aspect de l'invention, il est proposé un support d'enregistrement, lisible par un processeur, caractérisé en ce qu'il est apte à mémoriser un programme d'ordinateur 12 comportant des instructions pour la mise en oeuvre d'un procédé de précodage d'un signal à transmettre via un canal de communication tel que décrit ci-dessus ou d'un procédé de démodulation d'un signal reçu après propagation via un canal de communication tel que décrit ci-dessus. Le support d'enregistrement est éventuellement amovible, mémorisant respectivement un 5 programme d'ordinateur mettant en oeuvre un procédé d'estimation de canal, ou de précodage ou de démodulation tels que décrits précédemment. 4. Liste des figures D'autres avantages et caractéristiques de l'invention apparaîtront plus clairement à la 10 lecture de la description suivante d'un mode de réalisation particulier de l'invention, donné à titre de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : la Figure 1 illustre de façon schématique un système comprenant un émetteur et des récepteurs pouvant intégrer un mode de réalisation de l'invention ; les Figures 2 et 3 illustrent respectivement des procédés de précodage et de démodulation 15 implémentant un mode de mise en oeuvre selon l'invention ; la Figure 4 illustre plus précisément et de façon schématique un système émetteur/récepteur pouvant implémenter un mode de réalisation de l'invention ; la Figure 5 illustre à titre d'exemple, un modèle général de représentation du système de communication schématisé sur la figure 4 ; 20 - la Figure 6 représente un mode de mise en oeuvre d'un procédé d'estimation du canal de communication selon l'invention ; les Figures 7a et 7b représentent des mises en oeuvre de deux alternatives d'une d'étape préalable de validation d'un modèle d'ondes, selon l'invention ; la Figure 8 illustre un exemple de mise en oeuvre de l'étape de réalisation d'une estimation 25 de canal, représenté sur de la figure 6 ; la Figure 9 illustre un autre exemple de mise en oeuvre de l'étape de réalisation d'une estimation de canal, représenté sur de la figure 6 ; et la Figure 10 illustre une manière particulière de réaliser un moyen de traitement pour un émetteur ou un récepteur. 30 13 5. Description d'un mode de réalisation particulier de l'invention La Figure 1 illustre un exemple de système de communication dans lequel s'inscrit l'invention. 5 Dans cet exemple, un émetteur Tx ou un émetteur/récepteur tel qu'une station de base, communique avec deux récepteurs Rx1 et Rx2 (ou émetteurs/récepteurs) respectivement formé d'une ou d'une pluralité d'antennes. Un canal de communication sans fil comprenant un ou plusieurs trajets (encore appelés chemins) relie l'émetteur Tx à chacun des récepteurs Rx1 et Rx2. Le canal de communication (ou canal de transmission) permet de transmettre des signaux depuis 10 l'émetteur vers chacun des récepteurs, le récepteur ayant accès au signal émis déformé par son passage par le canal, et auquel du bruit thermique s'est éventuellement ajouté. Il est connu que chaque trajet constituant le canal est associé à un gain complexe, défini par un déphasage et une atténuation. Le premier récepteur Rx1 de dimension Ll (par exemple une longueur caractéristique du premier 15 récepteur, ou alors le diamètre d'un cercle dans lequel s'inscrit l'ensemble des antennes du premier récepteur) est localisé à une distance D1 de l'émetteur Tx. Cette distance est suffisamment grande pour que l'hypothèse ondes planes soit applicable, malgré la taille du récepteur. A l'inverse, le second récepteur Rx2 de dimension L2 est situé à une distance D2 de l'émetteur Tx. La dimension L2 correspond par exemple au diamètre d'un cercle dans lequel s'inscrit un système 20 de 256 antennes associées entre-elles si l'on se place dans le cas d'un récepteur « massive MIMO ». La longueur D2 est insuffisante pour permettre l'application de l'hypothèse ondes planes. En effet la courbure des ondes symbolisées par les cercles concentriques, est trop importante du point de vue du récepteur Rx2. 25 Pour un signal à transmettre depuis un émetteur vers un récepteur, un précodage du signal à transmettre est effectué pour maximiser le débit de données sur le canal de communication sans fil. Comme illustré sur la Figure 2, une première étape 21 d'estimation du canal de communication sans fil est réalisée, afin de connaître les paramètres définissant ce canal de communication. Puis un précodage 22 du signal à transmettre via ce canal de communication est effectué à partir des 30 paramètres obtenus lors de l'étape d'estimation du canal de communication. 14 Pour un signal reçu par un récepteur depuis un émetteur, une démodulation du signal reçu via le canal de communication est effectuée de façon à maximiser le débit de données sur le canal. Comme représenté sur la Figure 3, une première étape 31 d'estimation du canal de communication sans fil est réalisée, afin de connaître les paramètres définissant ce canal de 5 communication, à partir de signaux émis connus. Puis une démodulation 32 du signal reçu du canal de communication est effectuée à partir du canal estimé, ou autrement dit des paramètres obtenus lors de l'étape d'estimation du canal de communication. La Figure 4 illustre très schématiquement un dispositif implémentant un mode de 10 réalisation de l'invention. Un émetteur Tx est relié à un récepteur RA via un canal de communication sans fil Cn symbolisé par un rectangle. L'émetteur Tx comprend un moyen de traitement MTx configuré entre autre pour mettre en oeuvre les étapes 21 et 22 décrites en référence à la Figure 2. Le récepteur Tx comprend un moyen de traitement MRA configuré entre autre pour mettre en 15 oeuvre les étapes 31 et 32 décrites en référence à la Figure 3. La Figure 5 illustre à titre d'exemple, un modèle général de représentation du système de communication schématisé sur la figure 4, dans un cas d'une transmission en ligne de vue directe (« line of sight » ou LoS en langue anglaise). Un point O, (respectivement 04 correspond au centre de gravité du réseau d'antennes d'émission 20 (respectivement de réception). Des vecteurs ;id a',,v).; et correspondent aux positions des antennes d'émission et de réception par rapport à leurs centres de gravité respectifs (des repères locaux différents sont utilisés pour l'émission et la réception). De plus, une distance de référence D correspond à la distance qui sépare les deux centres de gravité Of et O,. 25 La notation Dij désigne la distance qui sépare la j-ème antenne d'émission et la i-ème antenne de réception. Dans l'exemple de la figure 5, la distance 0,3 entre la deuxième antenne de l'émetteur et la troisième antenne du récepteur est représentée. ut correspond à la direction de propagation des ondes au départ du canal de communication dans le repère centré sur O, et e à la direction de propagation des ondes à l'arrivée du canal de 30 communication dans le repère centré sur 0,-. 15 En faisant l'hypothèse d'une atténuation et d'une phase proportionnelles à la distance Du parcourue par les ondes, le canal de transmission à la fréquence f entre la j-ème antenne d'émission et la i-ème antenne de réception s'exprime de la façon suivante : D rio 5 où :7 est la longueur d'onde, c correspondant la vitesse de la lumière, et h correspond au canal entre les points ()let Cl,. De plus, les distances D et Di' étant assez proches en pratique, l'hypothèse D- ï est faite. Enfin, DiJ la distance Di; s'exprime de la manière suivante e47e-Di ND2 2.11,7;1,7 - i72/.7t) Rj-à1 2 10 où R est une matrice de rotation permettant de mettre en correspondance les repères locaux utilisés à l'émission et à la réception. En considérant comme exacte cette expression pour la distance, les inventeurs ont pu établir une nouvelle définition du modèle ondes sphériques (SWM), telle qu'exprimé ci-dessous : Asw?,,Eu .11 Du - D V D2 +2D(ar,hur - at,ima - 15 On obtient donc pour expression du canal : Classiquement, du fait de l'hypothèse du modèle ondes planes, les modèles physiques de l'état de l'art utilisés pour l'estimation de canal approchent l'expression 'Aswmii par Apiikisttij - 20 Cette expression beaucoup moins complexe que celle utilisée dans le cas de l'hypothèse du modèle ondes sphériques, pose le problème du manque de précision dans le cas de récepteur et/ou émetteur comportant de larges antennes. Pour résoudre ce problème de complexité (et donc de coût en terme de temps de calcul) posé par le modèle d'ondes sphériques et d'imprécision posé par le modèle d'onde plane, les inventeurs ont 25 proposé un nouveau modèle d'ondes intermédiaire dit modèle ondes paraboliques (ParWM), tel que défini ci-dessous : 16 AParWM.ij = c -÷ -4 Fitkir> (au 0 at,pui 2D _1 r- dir at,j42 Ce nouveau modèle constitue un compromis entre le modèle ondes planes PWM et le modèle ondes sphériques SWM : il est plus précis que le PWM, et moins complexe que le SWM. En effet, l'expression de Aparwm.ij ne comporte que des produits et des additions alors que celle de 5 Assmij comporte une racine carrée. Dans la suite de la description, à des fins de simplification, l'invention sera placée dans le contexte particulier dit MISO (« Multiple In Single Out » en langue anglaise) multi-trajets, c'est-à-dire que l'on considère que le récepteur comprend une seule antenne. Bien entendu, l'homme du métier saura facilement étendre l'invention au contexte MIMO en considérant l'expression générale 10 de la différence de distance Du ou l'une de ses approximations proposées ci-avant. Dans le contexte MISO, le récepteur n'a qu'une seule antenne, )11,- = I. En conséquence, la matrice de Rotation R, facteur du vecteur représentatif de la direction d'arrivée du signal sur les antennes du récepteur, disparait de l'équation. Cela correspond par exemple à des réseaux 5G conçus pour des stations de bases ayant un grand nombre d'antennes et des terminaux mobiles 15 n'ayant qu'une seule antenne. Dans ce cas, les expressions données précédemment permettent de prendre en compte non seulement les canaux de communication constitués d'un unique trajet en vue directe mais aussi les trajets issus de réflexions sur des plans en vue indirecte. Chaque trajet est associé à un gain complexe correspondant au canal entre les centres de gravité des réseaux d'antennes d'émission et de réception. 20 Les inventeurs ont proposé une expression générale du canal de communication comprenant p trajets. Cette expression est valable pour les trois modèles considérés SWM, PWM et ParWM : h,,,, = 1JiVr , hkem(i7,71, DO. où At désigne le modèle physique choisi, 25 hk, -17,7"; et Dk sont respectivement le canal de transmission entre les centres de gravité des réseaux d'antennes d'émission et de réception, la direction de départ lors de la propagation des ondes pour le k-ième trajet et la distance parcourue par les ondes pour le k-ième trajet, et 17 em(ii7t E C' est un vecteur dit caractéristique du modèle physique M pour ce k-ième trajet. Dans le cas où l'émetteur et le récepteur comprennent chacun plusieurs antennes, le canal de communication est alors défini par une matrice caractéristique dont les dimensions sont données par le nombre d'antennes à l'émetteur et au récepteur. Dans l'exemple considéré où le récepteur 5 comporte une seule antenne, la matrice caractéristique comprend une seule colonne. On parle alors de vecteur caractéristique. Ce vecteur caractéristique s'exprime de manière générale comme suit : Avec pour les modèles ondes sphériques SWM et paraboliques ParWM : 10 ASWnik 1/1)i + 211)k ( dir,k) lai, je, - AParSVM,ik -a?,1,111,1 Dans le cas du modèle d'ondes planes PWM avec tspwm,ii D et du modèle d'ondes sphériques si D ,7( 10 Cette alternative est simple à mettre en oeuvre, mais elle a toutefois l'inconvénient de se baser sur un écart de déphase arbitraire de 1' De plus elle est indépendante de la position (une fois la distance obtenue) et de l'orientation relative du récepteur par rapport à l'émetteur. Une autre alternative illustrée sur la Figure 7b permet de résoudre ces problèmes. Cette 15 alternative repose sur une métrique appelée erreur relative d'approximation du modèle (rMAE), définie ci-dessous : ilh -140i/51(111W, où prois.i(u) A argminx," x112 , où Al est le modèle à valider et h est le canal de référence, ici le canal obtenu avec le modèles ondes sphériques. L'erreur relative d'approximation du modèle 20 quantifie une l'erreur d'approximation minimale qu'implique le modèle Al considéré. On peut considérer que le modèle Al est valide lorsque la rMAE est faible (par exemple si rMAE également 15 associée au premier trajet. Un résidu l'a' obtenu tel que décrit ci-dessus est délivré en entrée d'une deuxième étape 82. Cette deuxième étape comporte une deuxième estimation pour un deuxième trajet à partir des valeurs déterminées de la distance v, associée au deuxième trajet, de la direction également associée au deuxième trajet et du résidu 1-2'obtenu lors de l'étape 81. La matrice E(') est délivrée à l'issue de l'étape 82, permettant de déterminer le canal estimé 20 fl à partir de l'équation précédemment mentionnée : A noter que dans le cas où un seul trajet est considéré, seule l'étape 81 est effectuée. L'avantage de cette alternative est une grande précision dans le canal estimé obtenu. Cependant, la réalisation de l'estimation de canal est possible en construisant un dictionnaire de N--,N,r) vecteurs caractéristiques (correspondant à directions de départ et ND distances). La 25 complexité de la résolution approchée de ce problème est donc O(NI--)ND). En effet, cette alternative implique de tester toutes les distances pour chacune des directions, ce qui induit une forte complexité calculatoire. 22 Une autre mise en oeuvre moins coûteuse est possible en estimant le vecteur caractéristique du modèle choisi (ondes sphériques ou paraboliques) via la résolution du problème d'optimisation irglitaxi XeMNI.07 4 1,A Xem(uLk,D)1 #nruzlx D 11XCA1(1e:f.D){12 5 Cette autre mise en oeuvre correspond à une estimation séquentielle des directions de propagation respectivement associées à chaque trajet, puis des distances également associées à chaque trajet. Sa complexité est 0(N-, +ND). En effet, la stratégie revient à ne tester plusieurs distances que pour la meilleure direction de propagation. Cette meilleure direction de propagation peut être par exemple déterminée avec des fronts d'ondes plans. 10 La Figure 9 illustre un exemple de cette autre mise en oeuvre, par exemple pour la première estimation 81 de la figure précédente. Une première étape 91 de sous-estimation consiste à fixer la valeur de la distance DI et à déterminer la meilleure direction 71:r (celle qui maximise la fonction de coût). Lorsque cette meilleure direction de propagation est déterminée, 15 une deuxième sous-estimation 92 revient à tester plusieurs distances pour cette meilleure direction de propagation. A l'inverse, dans un mode de mise en oeuvre alternatif, la première sous-estimation peut consister à trouver la meilleure valeur de distance possible pour une direction de propagation donnée (par exemple dans le cas d'une propagation très directive), puis lors de l'étape de deuxième sous-estimation tester les directions de propagation pour cette meilleur valeur de 20 distance obtenue lors de l'étape de première sous-estimation. La Figure 10 illustre une manière particulière, parmi plusieurs possibles, de réaliser un moyen de traitement 100 (pour un émetteur ou un récepteur) configuré pour mettre en oeuvre un mode de réalisation d'un procédé selon l'invention. Le dispositif 100 comprend une mémoire vive 25 130 (par exemple une mémoire RAM), une unité de traitement 110 équipée par exemple d'un processeur, et pilotée par un programme d'ordinateur stocké dans une mémoire morte 120 (par Ir(c)11Xerwmtrit)1 23 exemple une mémoire ROM ou un disque dur). A l'initialisation, les instructions de code du programme d'ordinateur sont par exemple chargées dans la mémoire vive 130 avant d'être exécutées par le processeur de l'unité de traitement 110. La figure 10 illustrent seulement une manière particulière, parmi plusieurs possibles, de 5 réaliser le moyen de traitement 100 afin qu'il effectue certaines étapes du procédé selon l'invention. En effet, ces étapes peuvent être réalisées indifféremment sur une machine de calcul reprogrammable (un ordinateur PC, un processeur DSP ou un microcontrôleur) exécutant un programme comprenant une séquence d'instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre module 10 matériel). Dans le cas où le moyen de traitement est réalisé avec une machine de calcul reprogrammable, le programme correspondant (c'est-à-dire la séquence d'instructions) pourra être stocké dans un médium de stockage amovible ou non, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. 15 Il va de soi que les modes de mise en oeuvre et de réalisation qui ont été décrits ci-dessus ont été donnés à titre purement indicatif et nullement limitatif, et que de nombreuses modifications peuvent être facilement apportées par l'homme de l'art sans pour autant sortir du cadre de l'invention. REVENDICATIONS1- Procédé d'estimation d'un canal de communication sans fil entre un émetteur et un récepteur, comprenant une pluralité de trajets permettant la propagation d'une onde, au moins l'émetteur ou le récepteur étant formé de plusieurs antennes, le procédé comprenant : pour au moins un trajet, une détermination d'une matrice dite caractéristique, qui dépend d'un premier élément représentatif d'au moins une direction de propagation (e, re,-) associée au trajet, et d'un deuxième élément représentatif d'une distance de propagation associée audit trajet, et une réalisation d'une estimation du canal de communication (fi) à partir de la au moins une matrice caractéristique obtenue. 2 - Procédé selon la revendication 1, dans lequel la matrice caractéristique associée au trajet dépend en outre d'au moins l'un des deux éléments suivants : - au moins un vecteur (RI ii7X1) reliant le centre de gravité du réseau d'antennes de l'émetteur et au moins une antenne de l'émetteur, - au moins un vecteur (14, ... reliant centre de gravité du réseau d'antennes du récepteur et au moins une antenne du récepteur. 3- Procédé selon la revendication 1, dans lequel la définition d'une matrice caractéristique dépend d'un modèle de canal parmi une pluralité de modèles. 4- Procédé selon la revendication 3, dans lequel la pluralité de modèles comprend le modèle ondes sphériques et un modèle ondes dites paraboliques décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne du récepteur d'une part, et la distance existant entre les centres de gravité des réseaux d'antennes du récepteur et de l'émetteur d'autre part, le modèle d'ondes dites paraboliques étant fonction uniquement : 25 d'une grandeur dépendant d'au moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, un vecteur de direction d'arrivée de propagation associée au trajet, et d'un terme correctif fonction de la distance de propagation associée audit trajet, 5 ladite grandeur et ledit terme correctif correspondant respectivement audit premier et audit deuxième élément. 5- Procédé selon la revendication 4, incluant en outre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, qui comprend 10 une obtention (71a) de la distance entre les centres de gravité des réseaux d'antennes du récepteur d'une part et de l'émetteur d'autre part, une comparaison (72a) de ladite distance à un premier et un second seuil, une validation (73a) d'un modèle ondes planes si la distance obtenue est supérieure au second seuil, du modèle ondes dites paraboliques si la distance obtenue est comprise entre le premier et le 15 second seuil, et le modèle ondes sphériques si la distance obtenue est inférieure au premier seuil. 6- Procédé selon la revendication 5, dans lequel le premier seuil est fonction de la distance à la puissance 3/2 entre le centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur, et 20 le second seuil est fonction de la distance au carré entre le centre de l'antenne de l'émetteur la plus éloignée du centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur. 7- Procédé selon la revendication 4, incluant en outre une étape préalable de validation du modèle 25 de canal dont dépend la matrice caractéristique, qui comprend une détermination (71b) d'une erreur relative d'approximation du modèle pour au moins un premier et un deuxième modèle, chaque erreur relative d'approximation du modèle étant fonction d'un écart entre un modèle de référence et l'un desdits deux modèles 26 une comparaison (72b) de chaque erreur relative d'approximation du modèle à une valeur prédéterminée de marge d'erreur, et une validation ou non (73b, 74b) du premier modèle de canal ou du deuxième modèle de canal en fonction du résultat de la comparaison. 5 8- Procédé selon l'une quelconque des revendications précédentes, dans lequel la réalisation d'une estimation du canal de communication comprend au moins pour un premier trajet : Une première estimation (81) d'une première valeur de distance de propagation associée au premier trajet, et d'une première direction de propagation associée au 10 premier trajet, et une détermination d'une première matrice caractéristique dépendant de ladite première valeur de distance de propagation et de ladite première direction de propagation estimées. 9 - Procédé selon la revendication 8, comprenant en outre pour un deuxième trajet : 15 Une deuxième estimation (82) d'une deuxième valeur de distance de propagation associée audit deuxième trajet, et d'une deuxième direction de propagation associée audit deuxième trajet, et une détermination d'une deuxième matrice caractéristique dépendant de ladite deuxième valeur de distance de propagation et de ladite deuxième direction de propagation estimées, ainsi que d'un résidu résultant de la première 20 estimation. 10- Procédé selon l'une des revendications 8 ou 9, dans lequel au moins la première estimation ou la deuxième estimation est effectuée en une première sous-estimation (91), pour une première grandeur fixée et une deuxième 25 grandeur, distincte de la première grandeur, pouvant prendre une pluralité de valeurs prédéterminées, puis en une deuxième sous-estimation (92) pour la deuxième grandeur fixée à une valeur résultant de la première sous-estimation et la première valeur pouvant prendre une pluralité de valeurs prédéterminées. 11 - Procédé selon la revendication 10, dans lequel 27 la première grandeur ou la deuxième grandeur est l'un des éléments suivants : un élément représentatif de la distance de propagation associée au trajet considéré, ou un élément représentatif de la direction de propagation associée au trajet considéré. 5 12 - Procédé de précodage d'un signal à transmettre via un canal de communication, comprenant : - une estimation (21) du canal de communication selon l'une quelconques des revendications précédentes afin d'obtenir un canal de communication estimé (fi), et - un précodage (22) du signal à transmettre en fonction du canal de communication 10 estimé. 13 - Procédé de démodulation d'un signal reçu après propagation via un canal de communication, comprenant : une estimation (31) du canal de communication selon l'une quelconques des 15 revendications 1 à 11 afin d'obtenir un canal de communication estimé (fi), et - une démodulation (32) du signal reçu en fonction du canal de communication estimé. 14 - Dispositif d'estimation (100) d'un canal de communication sans fil entre un émetteur (Tx) et un récepteur (Rx:), comprenant une pluralité de trajets permettant la propagation d'une onde, au 20 moins l'émetteur ou le récepteur étant formé de plusieurs antennes, le dispositif comprenant : un moyen de détermination configuré pour mettre en oeuvre, pour au moins un trajet, une détermination d'une matrice dite caractéristique, qui dépend d'un premier élément représentatif d'au moins une direction de propagation 25 associée au trajet, et d'un deuxième élément représentatif d'une distance de propagation associée audit trajet, et un moyen de réalisation configuré pour mettre en oeuvre une réalisation d'une estimation du canal de communication (fi) à partir de la au moins une matrice caractéristique obtenue. 28 15 - Dispositif selon la revendication 14, dans lequel la matrice caractéristique associée au trajet dépend en outre d'au moins l'un des deux éléments suivants : au moins un vecteur (W,I retA,}) reliant le centre de gravité du réseau d'antennes de l'émetteur et au moins une antenne de l'émetteur, 5 au moins un vecteur (;a1,..'a«.1) reliant centre de gravité du réseau d'antennes du récepteur et au moins une antenne du récepteur. 16 - Dispositif selon la revendication 14, dans lequel la définition d'une matrice caractéristique dépend d'un modèle de canal parmi une pluralité de modèles. 10 17 - Dispositif selon la revendication 16, dans lequel la pluralité de modèles comprend le modèle ondes sphériques et un modèle ondes dites paraboliques décrivant un écart entre : la distance entre les positions d'au moins une paire d'antennes formée par une antenne de l'émetteur et une antenne du récepteur d'une part, et la distance existant entre les centres de gravité des réseaux 15 d'antennes du récepteur et de l'émetteur d'autre part, le modèle d'onde dite parabolique étant fonction uniquement : d'une grandeur dépendant tau moins l'un des éléments suivants : un vecteur de direction de départ de propagation associée au trajet, un vecteur de direction d'arrivée de propagation associée au trajet, 20 et d'un terme correctif fonction de la distance de propagation associée audit trajet, ladite grandeur et ledit terme correctif correspondant respectivement audit premier et audit deuxième élément. 18- Dispositif selon la revendication 17, incluant en outre un moyen de validation configuré pour 25 mettre en oeuvre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, ledit moyen de validation comprenant un moyen d'obtention configuré pour mettre en oeuvre une obtention de la distance entre les centres de gravité des réseaux d'antennes du récepteur d'une part et de l'émetteur d'autre part, 29 un moyen de comparaison configuré pour mettre en oeuvre une comparaison de ladite distance à un premier et un second seuil, un moyen de décision configuré pour mettre en oeuvre une validation d'un modèle ondes planes si la distance obtenue est supérieure au second seuil, du modèle ondes dites paraboliques si 5 la distance obtenue est comprise entre le premier et le second seuil, et le modèle ondes sphériques si la distance obtenue est inférieure au premier seuil. 19 - Dispositif selon la revendication 18, dans lequel le premier seuil est fonction de la distance à la puissance 3/2 entre le centre de l'ensemble 10 des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur, et le second seuil est fonction de la distance au carré entre le centre de l'antenne de l'émetteur la plus éloignée du centre de l'ensemble des antennes de l'émetteur et le centre de l'ensemble des antennes de l'émetteur. 15 20- Dispositif selon la revendication 17, incluant en outre un moyen de validation configuré pour mettre en oeuvre une étape préalable de validation du modèle de canal dont dépend la matrice caractéristique, qui comprend un moyen de détermination configuré pour mettre en oeuvre une détermination d'une erreur relative d'approximation du modèle pour au moins un premier et un deuxième modèle, 20 chaque erreur relative d'approximation du modèle étant fonction d'un écart entre un modèle de référence et l'un desdits deux modèles un moyen de comparaison configuré pour mettre en oeuvre une comparaison de chaque erreur relative d'approximation du modèle à une valeur prédéterminée de marge d'erreur, et un moyen de décision configuré pour mettre en oeuvre une validation ou non du premier 25 modèle de canal ou du deuxième modèle de canal en fonction du résultat de la comparaison. 21- Dispositif selon l'une quelconque des revendications 14 à 20, dans lequel ledit moyen de réalisation configuré pour mettre en oeuvre une réalisation d'une estimation du canal de communication comprend: 30 - un moyen d'estimation configuré pour mettre en oeuvre, pour au moins pour un premier trajet, une première estimation d'une première valeur de distance de propagation associée au 30 premier trajet, et d'une première direction de propagation associée au premier trajet, et une détermination d'une première matrice caractéristique dépendant de ladite première valeur de distance de propagation et de ladite première direction de propagation estimées. 5 22 - Dispositif selon la revendication 21, comprenant en outre: un deuxième moyen d'estimation configuré pour mettre en oeuvre, pour un deuxième trajet, une deuxième estimation d'une deuxième valeur de distance de propagation associée audit deuxième trajet, et d'une deuxième direction de propagation associée audit deuxième trajet, et une détermination d'une deuxième matrice caractéristique dépendant de ladite deuxième valeur de 10 distance de propagation et de ladite deuxième direction de propagation estimées, ainsi que d'un résidu résultant de la première estimation. 23 - Dispositif selon l'une des revendications 21 ou 22, dans lequel au moins le premier moyen d'estimation ou le deuxième moyen d'estimation est configuré pour effectuer 15 une première sous-estimation, pour une première grandeur fixée et une deuxième grandeur, distincte de la première grandeur, pouvant prendre une pluralité de valeurs prédéterminées, puis une deuxième sous-estimation pour la deuxième grandeur fixée à une valeur résultant de la première sous-estimation et la première valeur pouvant prendre une pluralité de valeurs 20 prédéterminées. 24 - Dispositif selon la revendication 23, dans lequel la première grandeur ou la deuxième grandeur est l'un des éléments suivants : un élément représentatif de la distance de propagation associée au trajet considéré, 25 ou un élément représentatif de la direction de propagation associée au trajet considéré. 25 - Système de précodage d'un signal à transmettre via un canal de communication, 30 comprenant : 31 un dispositif d'estimation du canal de communication selon l'une quelconques des revendications 14 à 24, afin d'obtenir un canal de communication estimé, ledit système de precodage étant configuré pour mettre en oeuvre un précodage du signal à transmettre en fonction du canal de communication estimé. 5 26 - Système de démodulation d'un signal reçu après propagation via un canal de communication, comprenant : un dispositif d'estimation du canal de communication selon l'une quelconques des revendications 14 à 24 afin d'obtenir un canal de communication estimé, 10 Ledit système de démodulation étant configuré pour mettre en oeuvre une démodulation du signal reçu en fonction du canal de communication estimé. 27 - Programme d'ordinateur comportant des instructions pour la mise en oeuvre des étapes d'un procédé selon l'une quelconques des revendications 1 à 11, lorsqu'il est exécuté par un processeur. 15 28 - Programme d'ordinateur comportant des instructions pour la mise en oeuvre des étapes d'un procédé selon l'une des revendications 12 ou 13 lorsqu'il est exécuté par un processeur. 29 - Support d'enregistrement, lisible par un processeur, caractérisé en ce qu'il est apte à 20 mémoriser le programme d'ordinateur selon l'une des revendications 27 ou 28. 25