FR 3046312 A1 20170630 FR 1502711 A 20151229 3046312 EMULATEUR DE CANAL DE TRANSMISSION L'invention concerne notamment un émulateur pour canal de propagation comportant une entrée pour la réception d'un signal et une sortie de signal correspondant à l'émulation du signal reçu en entrée, filtré après émulation du canal de transmission. Dans de nombreuses applications, les signaux sont transmis via un canal de transmission ou médium de transmission. Lors de cette étape, des phénomènes de propagation se produisent et peuvent affecter le signal. Un système de communication transmet des données sur un medium appelé canal de transmission. Le canal de transmission apporte des distorsions au signal transmis dues, notamment, à l'atténuation en espace libre, la diffraction et la diffusion, ou encore à la présence de trajets multiples de propagation dont les recompositions s'effectuent de manière constructive ou destructive. Du bruit d'origines diverses peut également impacter cette transmission. Le bruit est, par exemple, généré par des équipements connectés ou situés à proximité du canal de transmission. Lors de la conception d'un système de communication ou de tout autre système dans lequel il y a propagation de signaux au sein d'un medium, on cherche à établir un modèle ou des modèles de canaux représentatifs des phénomènes physiques de transmission. Plusieurs types de modèles existent. Les modèles empiriques statistiques reproduisent les fonctions caractéristiques du canal de transmission de manière aléatoire en utilisant un certain nombre de statistiques observées sur un ensemble de mesures. Pour les modèles déterministes, on utilise une description fine du milieu de transmission. La publication intitulée «FPGA-based Wideband Channel Emulator for Evaluation of Wireless Sensor Networks in Industrial Environments », 2014 IEEE Emerging Technology and Factory Automation (ETFA) Iriaki Val, Félix Casado, Pedro Manuel Rodriguez, Aitor Arriola, décrit la mise en oeuvre d'un émulateur numérique de canaux de propagation 2 3046312 « DWCE », digital wireless channel emulator. Celui-ci s'articule autour d'une architecture temps réel (Transposition en fréquence, convertisseurs analogiques/numériques, FPGA) avec un algorithme permettant la génération de multitrajet, de bruit et de spectre Doppler. Les principales 5 limites de ces travaux sont un temps de calcul long (>900ns), une faible résolution temporelle (>8ns), ainsi que le nombre de macro-trajets émulables très limité ou encore l'effet Doppler appliqué en moyenne sur tout le profil retard puissance. La littérature offre aussi différents algorithmes de traitement basés sur des méthodes de calcul de la réponse du canal de 10 propagation avec tirage aléatoire des coefficients selon des distributions prédéfinies. A priori, aucun algorithme de traitement numérique, à implémenter dans un émulateur et connu du demandeur, ne propose à ce jour de reproduire la réponse simultanée d'un canal de propagation en temps réel 15 avec discrimination de l'effet Doppler. Le profil retard-puissance-Doppler émulé est fixe en dehors de l'évanouissement ou fading aléatoire qui lui est appliqué. Le Doppler est appliqué sous la forme d'un spectre aléatoire ou moyen. Dans les cas les plus sophistiqués, des modèles de canaux dynamiques peuvent être émulés. Cela se fait au détriment de certains 20 aspects des performances, comme le retard minimal ou le nombre maximal de macro-trajets émulables. Ces deux variables jouent pourtant un rôle crucial dans le réalisme de l'émulation et la variété d'équipements testables. Les modèles présentés dans l'art antérieur, malgré les résultats et les avantages qu'ils procurent, sont en général simplistes et un peu éloignés 25 de la réalité. L'objet de l'invention concerne notamment un émulateur de canal de propagation permettant une émulation de canaux variant dans le temps avec une grande variété de trajets et leurs Doppler associés. L'originalité de l'émulateur selon l'invention repose sur un préconditionnement des données 30 à fournir à l'émulateur afin de simplifier le modèle de canal de propagation de trois dimensions vers deux dimensions, ainsi que de procéder à sa 3 3046312 linéarisation. Ceci permet en outre de simplifier les opérations à effectuer côté hardware en vue d'émuler le canal et par ce fait, de réduire la puissance de calcul nécessaire. D'autre part, les calculs étant plus rapides, le retard minimal émulable sera beaucoup plus faible qu'avec les méthodes 5 conventionnelles. Pour finir, cette méthode dynamique a l'avantage de permettre de recréer des environnements comportant un grand nombre de macro-trajets, tout en contrôlant indépendamment et de manière précise, leur retard, leur décalage doppler ainsi que leur atténuation. Les modèles de canaux émulés s'en retrouvent d'autant plus fidèles à la réalité qu'avec les 10 émulateurs conventionnels aux retards importants, au nombre de macrotrajets émulés faibles et à l'effet Doppler approximé. L'invention concerne un procédé mis en oeuvre dans un émulateur pour canal de propagation connecté à un équipement, l'émulateur comprenant au moins une entrée recevant un signal E(t) émis par ledit 15 équipement dans une direction m qui induit un retard nAT avec une fréquence Doppler d. Of , et une sortie émettant un signal R(t) après traitement, un processeur adapté à exécuter au moins les étapes suivantes : - Choisir une modélisation discrète du canal de propagation sur les temps d'arrivée du signal et sur les fréquences Doppler, 20 - Définir une ou plusieurs zones correspondant à un intervalle de temps et à une valeur de fréquence Doppler et déterminer une matrice de canal de propagation DT dont les éléments sont obtenus en sommant toutes les valeurs complexes des trajets situés dans un intervalle de temps et une valeur de fréquence-Doppler choisis, 25 - Numériser le signal E(t) émis, - Combiner un échantillon numérisé à une ligne de la matrice DT matrice retard-Doppler-coefficient du canal afin d'obtenir le signal R(t) correspondant à un ensemble de trajets propagés, - Transmettre le signal R(t) à l'équipement. 30 Selon une variante de réalisation, le processeur est adapté à insérer dans un registre à décalage, à chaque temps symbole, un échantillon 4 3046312 du signal E(t), à faire subir à la matrice DT une rotation circulaire dans le sens des lignes et à réaliser une convolution entre la dernière ligne de la matrice et les échantillons mémorisés afin d'obtenir Ri échantillons qui sont ensuite convertis dans le domaine analogique avant transmission vers 5 l'équipement sous test. Le signal R(t) correspond, par exemple, à une somme sur toutes les valeurs de retard possibles n, du produit d'un échantillon du signal émis par la somme sur les valeurs possibles de Doppler de l'impact d'un trajet sur le signal qui suit ce trajet, Cn,d. ei.21.r.d.M(t-nàt) R(t) = (E(t - nAT). Dn(t - nAT)) 10 avec : (t)= LId u- rcn d - ej.2-rr.d.M.(01, I où d. Of. (t) est la valeur de fréquence Doppler pour un retard nAT, Cn,d un coefficient complexe représentant l'atténuation liée à la propagation en espace libre, (E(t nAT)) une valeur d'échantillon du signal reçu avec un 15 retard nAT. Selon une variante de réalisation, le procédé comporte une étape où l'on précalcule la matrice Dn et où l'on élimine les fréquences Doppler ayant un impact négligeable sur le canal de propagation. Selon une autre variante, on définit l'intervalle de temps pour 20 regrouper les trajets des signaux égaux à 1/Fs, où Fs est la fréquence d'échantillonnage de l'émulateur, et on détermine le signal R(t) de la manière suivante : R (i. At) = E ((i - n). At). Dn((i - n). At) n=0 n est un indice borné, i correspond à l'instant où l'on reçoit le signal. La taille d'une zone temps-décalage-Doppler peut être choisie en 25 fonction de la fréquence d'échantillonnage des convertisseurs et de la 5 3046312 fréquence d'activation des coefficients du canal dans le circuit logique programmable FPGA. L'invention concerne aussi un émulateur pour canal de propagation connecté à un équipement, l'émulateur comprenant un premier 5 module de transposition en fréquence, relié à un module de conversion analogique/numérique et numérique/analogique, lui-même relié à un module de traitement du signal, au moins une entrée recevant un signal E(t) émis par ledit équipement dans une direction m qui induit un retard nAT avec une fréquence Doppler d. Af et une sortie émettant un signal R(t) après 10 traitement, un processeur de signaux, une mémoire de stockage des signaux, l'émulateur est caractérisé en ce que le processeur est adapté à exécuter au moins les étapes du procédé selon l'invention afin de générer un signal émulé R(t) transmis à l'équipement. L'émulateur peut comporter un circuit logique programmable ou 15 FPGA (Field Programmable Gate Array). L'émulateur peut aussi comporter un processeur de signal numérique de type DSP. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit d'exemples de 20 réalisation donnés à titre illustratif et nullement limitatif, annexée des figures qui représentent : - Figure 1, un exemple de structure d'un émulateur numérique selon l'invention, - Figure 2A, un schéma d'un avion et de la propagation de signaux, 25 figure 2B, une représentation d'une réponse de canal calculée, figure 2C, la représentation des macro-trajets dans un diagramme retard-Doppler, - Figure 3A, une représentation des macro-trajets selon le procédé selon l'invention, figure 3B les données conditionnées dans un profil 30 retard-Doppler-puissance des macro-trajets, figure 3C une représentation de la matrice simplifiée de canal de propagation, figure 6 3046312 3D, une représentation de l'étape de convolution d'un échantillon du signal avec la matrice de canal, - Figure 4A, une représentation selon l'art antérieur de la combinaison des échantillons du signal avec la matrice de canal, figure 4B, la 5 transformation de la matrice de canal selon l'invention, figure 4C, une représentation des étapes de combinaison de l'échantillon de signal avec la matrice de canal pour obtenir le signal résultant final, et - Figure 5, un résumé des étapes mises en oeuvre par le procédé selon l'invention. 10 L'exemple qui suit donné à titre illustratif et nullement limitatif concerne un émulateur utilisé par un équipementier de radar. La figure 1 schématise un exemple de structure d'un émulateur numérique selon l'invention connecté à un système sous test. Le système sous test 100 comporte une sortie 101 émettant un signal E(t) et une entrée 15 102 qui reçoit le signal R(t) traité par l'émulateur 110. L'émulateur selon l'invention 110 comporte par exemple un premier module de transposition en fréquence 120, relié à un module de conversion analogique numérique et numérique analogique 130, en liaison avec un module de traitement du signal de type FPGA 140 dans lequel est 20 implémenté l'algorithme ou procédé selon l'invention. Le signal E(t) émis par le système sous test est transmis vers le module de transmission en fréquence qui comporte une première entrée 121 connectée à une première chaine de transposition en fréquence 122. Après transposition en fréquence, le signal transposé est transmis par une sortie 25 123 vers une entrée 131 du module convertisseur 130 et vers un convertisseur analogique numérique 135. Le signal transposé et numérisé ETN(t) est émis par une sortie 136 vers une entrée 141 du FPGA 140. Le signal transposé et numérisé est traité par le procédé selon l'invention, puis émis par une sortie 142 du FPGA vers une entrée 137 et un convertisseur 30 numérique analogique 138. Le signal analogique est ensuite émis par une sortie 139 du convertisseur vers le module de transposition en fréquence via 7 3046312 une entrée 124. Le signal traité par l'émulateur est alors retransposé en fréquence via une deuxième voie 125 et émis via une sortie 126 vers l'entrée du système sous test. Le FPGA 140 comprend notamment un processeur 143 adapté à 5 exécuter les étapes du procédé selon l'invention et une mémoire 144 ou buffer permettant de mémoriser les échantillons e(t) du signal qui sont combinés avec la matrice de canal. Les échantillons en provenance du convertisseur analogique numérique 135, par le biais de l'entrée 141 du FPGA, sont synchronisés avec le processeur de traitement numérique 143 10 par le biais d'un registre à décallage (FIFO) 145. Le processeur de traitement numérique 143 reçoit les coefficients ci de la dernière ligne de la matrice de canal DT qui sont chargés dans la mémoire 144 (RAM) par le biais d'un processeur MicroblazeTM 146, par exemple. Le processeur 146 réalise les permutations circulaires de DT et actualise la mémoire 144 avec les bons 15 coefficients. Le processeur de traitement numériqué 143 réalise l'opération de convolution entre les échantillons en provenance de l'entrée 141 et les coefficients issus de la mémoire 144. Le résultat de ce calcul est envoyé vers un deuxième registre à décalage (FIFO) 147, dont le rôle est de synchroniser la sortie du processeur 143 avec l'entrée du convertisseur numérique 20 analogique 138. En résumé, le signal émis par l'équipement en test est filtré selon le procédé de l'invention, pour émuler le canal de propagation, puis restitué en direction de l'équipement. L'équipement en test utilisera le signal reçu R(t) afin d'effectuer les mesures propres à sa fonction, par exemple une mesure 25 de hauteur pour un radioaltimètre, une détection de cibles pour un radar, etc. Sans sortir du cadre de l'invention, le FPGA peut être remplacé par tout type de dispositif comportant au moins une mémoire ou buffer et un processeur dans lequel les étapes du procédé selon l'invention sont implémentées, par exemple un ordinateur PC, un processeur de signal 30 numérique DSP, etc. 8 3046312 La figure 2A schématise les signaux émis, transmis et réémis entre un avion 200 équipé d'un radar 201 et un objet au sol 202. Un signal E(t) émis par le radar se propage via le medium de transmission, il est réfléchi par une des facettes 203 de l'objet au sol 202. 5 La figure 2B est une représentation des modèles de canaux non linéairement distribués dans le temps (l'axe des ordonnées correspond à l'amplitude d'un macro-trajet, les deux autres axes correspondant au Doppler (Hz) et aux retards (secondes), et la figure 2C une mise en évidence de l'irrégularité des données en fréquence et en temps (retard-Doppler) ce qui 10 rend leur utilisation dans une chaîne de traitement numérique difficile du fait de l'existence d'une période d'échantillonnage fixe. La figure 3A illustre une étape du procédé qui permet de régulariser la distribution des données en fréquence et en temps. Le procédé selon l'invention, lors d'une première étape, va préconditionner les signaux 15 en regroupant, par exemple, les trajets en zones de retard et Doppler linéairement espacées. Le choix de la taille d'une zone se fera en lien avec la fréquence d'échantillonnage des convertisseurs analogique/numérique et la fréquence d'actualisation des coefficients du canal dans le FPGA comme il est explicité ci-après. La figure 3A schématise dans un diagramme (retard, 20 Doppler) le résultat de cette étape où l'on a illustré plusieurs macro-trajets 310j au sein d'une même zone retard-décalage-Doppler 320. Le regroupement permet de construire une matrice de la forme retard-dopplercoefficient du canal modélisé qui est adaptée à la configuration matérielle de l'émulateur, notamment à sa fréquence d'échantillonnage Fs et à la 25 fréquence d'actualisation Fa des coefficients de filtrage dans le FPGA. D'autre part, le regroupement permet de répartir linéairement les échantillons en retard et en Doppler, ce qui diminue fortement la complexité de l'implémentation et de la charge de calcul dans le FPGA. La figure 3B schématise les données conditionnées dans un profil 30 retard-Doppler-puissance des macro-trajets. 9 3046312 La figure 3C illustre la matrice simplifiée du modèle. La matrice DT est une matrice diagonale, chacun des éléments de la matrice de canal correspond à la somme de toutes les valeurs complexes des trajets situés dans un intervalle de temps (fixé par 1/Fs, avec Fs par exemple la fréquence 5 d'échantillonnage des convertisseurs analogique/numérique et numérique analogique) et une valeur de Doppler choisis. La figure 3D schématise une étape de convolution, c'est-à-dire une représentation de la mise en oeuvre de l'algorithme de traitement numérique qui considère les échantillons arrivant dans le FPGA en les 10 combinant avec les lignes de la matrice avec une rotation circulaire selon les étapes décrites ci-après. Le processeur va insérer dans le registre à décalage, à chaque temps symbole un échantillon e(t) du signal E(t) transposé et numérisé, faire subir à la matrice DT une rotation circulaire dans le sens des lignes et réaliser une convolution entre la dernière ligne de la 15 matrice et le buffer d'échantillons, afin d'obtenir Ri échantillons qui sont ensuite convertis dans le domaine analogique avant transmission vers l'équipement sous test. Les étapes du procédé mis en oeuvre par l'émulateur selon l'invention vont maintenant être détaillées en relation avec les figures 4A à 20 4C et la figure 5. On choisit une modélisation du canal discrète aussi bien sur les temps d'arrivée que sur les fréquences Doppler. Soit E(t) le signal émis par l'équipement à l'instant t. Ce signal est émis dans une direction m qui va induire un retard n AT avec une fréquence Doppler égale à d * M (Modèle 25 discret). Le signal va percuter un élément de surface 203 situé dans la direction m (figure 2A). On note Rf(t) le signal reçu par cet élément de surface : Rf(t) = 5 (t n A (E(t). ei'2ind2:22et) (1) 2 T) 10 3046312 avec Ô() la fonction de Dirac, n * AT le retard du signal dans une direction donnée de l'espace, e(t) le signal émis à l'instant t, d*Af la fréquence de décalage doppler liée aux mouvements entre l'avion et l'élément sol. Ceci veut dire que la facette 203 de l'objet percuté par le signal émis par 5 l'avion reçoit un signal déphasé par l'effet Doppler (terme en ei2nd*Aft) et retardé de12-*-3T (terme en 8 2 - On considère que l'atténuation liée à la 2 propagation en espace libre intervient en même temps que la modification du signal par sa réflexion sur l'élément de surface considéré, le tout étant représenté par un coefficient complexe Cn. Soit Ef(t) le signal réfléchi par 10 l'élément de surface : *A Ef(t) = Cn. [ (t n T) * (E(t). ei2ltd4I (2) avec Cn un coefficient inhérent au modèle de canal choisi, qui reflète l'atténuation et le déphasage du signal ayant emprunté le macro-trajet considéré. Soit Rn(t) le signal réfléchi par la facette de l'objet dans la direction associée 15 à ce macro-trajet: Rn (t) = [5 n *21 * (Cg. [5 n * AT) * (E(t). d' t)1)1 n±e't) (3) 2 Le signal réfléchi par l'objet est de nouveau retardé par sa propagation et subi à nouveau l'effet Doppler. L'équation 3 se simplifie en développant : n * AT) Rn(t) = [6 (t 2 . d*if n*AT, )1 cl*Af n*Én, e 2 2 e. 2 2 puis : n * OT n * AT) . cl*Af n*Éa &hl n*AT Rn(t) = C. E ((t 2 ) 2 . Soit : n * AT) * (Cn. E 2 Rn(t) = Cn. E(t - n * AT). ej.2n.d*Af.(t-n*AT) (6) 11 3046312 Ce qui revient à écrire le signal reçu sur le radar et après réflexion sur une des facettes, selon la direction n par la formule suivante : Rn(t) = Cn. 5(t - n * AT) * E(t). eilltd*Af.t (7) En pratique, le signal reçu R(t) par le radar après propagation dans le medium de transmission et diffraction par l'objet est composé d'un 5 ensemble de trajets qui se sont propagés dans plusieurs directions de l'espace, soit une somme de trajets Rn. Pour obtenir le signal R(t), on va sommer les trajets à la fois en retard et en Doppler suivant l'équation : R(t) = / 1 [(Cn,d. 6(t - n. AT)) * E(t)ei.2'd.°1 (8) n d pour toutes les directions de l'espace considérées. 10 L'ensemble des trajets est ici distribué temporellement sur n. AT et fréquentiellement sur d. M, correspondant bien à une représentation discrète du canal de propagation. Le signal R(t) est représenté par : R(t) = E [C En 7 £-Jdt-n,d- E(t - nAT). ei.2-rt.d.M.(t-nAT)1 J (9). 15 Le signal R(t) correspond à une somme sur toutes les valeurs de retards possibles n, du produit d'un échantillon du signal émis par la somme sur les valeurs possibles de Doppler de Cn,d. ei.2n.d.M.(t-nAT) (impact d'un trajet sur le signal qui l'emprunte) : R(t) = En(E(t - nAT). Ed[Cn,d. ej.2n.d.M.(t-nAT)1 \ j i (10) . 20 On pose la matrice de canal correspondant à un retard n: = ., n,d. d[- D n (t) 7 C ej.2n.d.M.(t)1 i (11). Ce qui permet d'écrire : R(t) = 1(E(t - nAT). Dn (t - nAT)) (12). n La matrice Dn est analogue à une transformée de Fourier discrète à un coefficient prés, qui est réalisée sur l'espace fréquentiel occupé par le 25 Doppler, les variables n. AT et d. M sont linéairement espacées. 12 3046312 La matrice Dr, est une fonction périodique de période max(D, Dn est linéairement distribué dans le domaine fréquentiel, ce qui équivaut à une périodicité dans le domaine temporel et simplifie les calculs. Selon un mode de réalisation, la matrice Dn peut être précalculée 5 pour limiter les temps de calcul dans l'émulateur. Sa taille est bornée du fait de sa périodicité. Il est alors possible d'éliminer les fréquences Doppler trop basses, qui ont un effet négligeable par rapport aux autres fréquences plus élevées, pour réduire sa taille et limiter son occupation mémoire. Ceci revient à dire que l'on peut adapter la valeur de M pour réduire la taille de Dn. 10 Du point de vue discret (cas de l'émulateur), cette simplification permet avantageusement de réduire les temps de calculs à l'aide du pré-conditionnement des variables. Dans le domaine du temps discret, l'équation (6) devient : R (i. At) = E (i. At - ?LOT). Dn(i. At - nAr) (13) n=0 avec i correspondant à l'instant auquel le signal est reçu. 15 L'indice n est désormais borné pour préserver la causalité de l'expression. En effet, il est impossible de recevoir en R à l'instant i. At une composante provenant d'un signal E à un temps plus récent. En posant At = Ai = 1 - avec FS la fréquence d'échantillonnage de l'émulateur, Fs l'intervalle de temps choisi pour le regroupement des trajets est égal à At. 20 Aucune contrainte particulière n'est imposée lors du choix de l'intervalle de fréquence Lf pour les fréquences Doppler. Cette condition posée, l'équation (7) se simplifie pour devenir : R (i. At) = En=0 E((i - n). At). D7,((i - n). At) (14). En simplifiant avec des indices pour alléger la notation : Ri = (i - n). Dri(i - (0) n=0 13 3046312 somme de n multiplications entre un échantillon de signal et un élément de la matrice. La figure 5 schématise les étapes mises en oeuvre par l'invention. La matrice de coefficients du filtre est représentée en 500, cette matrice est 5 diagonale et ces éléments correspondent à la somme des valeurs complexes des trajets situés dans un intervalle de temps et une valeur de fréquence-Doppler choisis. Le signal E(t) est échantillonné à un instant t, 501. Puis le procédé combine les échantillons du signal 501 obtenus à un instant donné t avec une ligne de la matrice 502 afin d'obtenir le signal R(t) avec t=0, et ainsi 10 de suite pour tous les instants t, dans l'exemple t=0 à 4, 503, 504, 505 afin d'obtenir le signal Ri(t). L'émulateur selon l'invention est notamment utilisé par les équipementiers de radars et de produits de radiocommunication. Il peut aussi être mis en oeuvre pour des instruments de mesure radiofréquence RF, et 15 dans les domaines pour lesquels la connaissance de l'impact du canal de propagation sur un produit est importante. L'émulateur selon l'invention offre notamment comme avantages de diminuer les temps de calcul nécessaires au traitement numérique tout en offrant la possibilité d'associer à chaque trajet un Doppler qui lui est propre. 20 D'autre part, l'actualisation des valeurs du filtre permet de faire varier la réponse impulsionnelle à émuler et permet ainsi de reproduire fidèlement des scénarios de fonctionnement d'un équipement. La profondeur du filtre de convolution est variable et peut être réduite à quelques échantillons lorsque le délai minimum à émuler est court, et augmentée à plusieurs centaines 25 d'échantillons pour des scénarii mettant en jeu des distances et donc des retards plus importants. 14 REVENDICATIONS1 - Procédé mis en oeuvre dans un émulateur (110) pour canal de propagation connecté à un équipement (100), l'émulateur comprenant au moins une entrée recevant un signal E(t) émis par ledit équipement (100) dans une direction qui induit un retard nAT avec une fréquence Doppler d. M, et une sortie émettant un signal R(t) après traitement, un processeur (143) adapté à exécuter au moins les étapes suivantes : - Choisir une modélisation discrète du canal de propagation sur les temps d'arrivée et sur les fréquences Doppler, - Définir une ou plusieurs zones (320) correspondant à un intervalle de temps et une valeur de fréquence Doppler, et déterminer une matrice de canal de propagation DT dont les éléments sont obtenus en sommant toutes les valeurs complexes des trajets situés dans un intervalle de temps et une valeur de fréquence-doppler choisis, - Numériser le signal E(t), - Combiner un échantillon numérisé du signal E(t) à une ligne de la matrice DT matrice retard-Doppler-coefficient du canal afin d'obtenir le signal R(t) correspondant à un ensemble de trajets propagés, - Transmettre le signal R(t) audit équipement (100). 2 - Procédé selon la revendication 1 caractérisé en ce que le processeur est adapté à insérer dans un registre à décalage, à chaque temps symbole, un échantillon du signal E(t), à faire subir à la matrice DT une rotation circulaire dans le sens des lignes et à réaliser une convolution entre la dernière ligne de la matrice et le buffer d'échantillons afin d'obtenir Ri échantillons qui sont ensuite convertis dans le domaine analogique avant transmission vers l'équipement sous test. 3 - Procédé selon l'une des revendications 1 à 2 caractérisé en ce que le signal R(t) correspondant à une somme sur toutes les valeurs de retard 15 3046312 possibles n du produit d'un échantillon du signal émis par la somme sur les valeurs possibles de Doppler de Cmd. ei.2n.d.M.(t-nAT): R(t) = DE(t - nAT). Dr,(t - MT)) avec : D,i(t) = d LC n,c1. ej.217.d.M.(t)1, impact d'un trajet sur le signal qui suit ce trajet, 5 d. M. (t) , la valeur de fréquence Doppler pour un retard nAT , Cn,d , un coefficient complexe représentant l'atténuation liée à la propagation en espace libre, (E(t - nAT)) une valeur d'échantillon du signal reçu avec un retard nAt. 10 4 - Procédé selon l'une des revendications 1 à 3 caractérisé en que l'on précalcule la matrice Dn et on élimine les fréquences Doppler ayant un impact négligeable sur le canal de propagation. 5 - Procédé selon l'une des revendications précédentes caractérisé en ce 15 que l'on définit l'intervalle de temps pour regrouper les trajets des signaux égal à 1/Fs, où Fs est la fréquence d'échantillonnage de l'émulateur, et on détermine le signal R(t) de la manière suivante : R (i. At) = ((i - n). At). Dn((i - n). At) n=0 n est un indice borné, i correspond à l'instant où l'on reçoit le signal. 20 6 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce que l'on choisit la taille d'une zone retard-Doppler en fonction de la fréquence d'échantillonnage des convertisseurs et de la fréquence d'activation des coefficients du canal dans un circuit logique programmable. 25 7 - Emulateur (110) pour canal de propagation connecté à un équipement (100), l'émulateur comprenant un premier module de transposition en 16 3046312 fréquence (120), relié à un module de conversion analogique/numérique et numérique/analogique (130) et à un module (140) de traitement du signal, au moins une entrée (141) recevant un signal E(t) émis par ledit équipement dans une direction qui induit un retard nAt avec une fréquence Doppler d. M 5 et une sortie (142) émettant un signal R(t) après traitement, un processeur de signaux (143), une mémoire (144) de stockage des signaux, caractérisé en ce que le processeur (143) est adapté à exécuter au moins les étapes du procédé selon l'une des revendications 1 à 6 afin de générer un signal émulé R(t) transmis audit équipement (100). 10 8 - Emulateur selon la revendication 7 caractérisé en ce qu'il comporte un circuit logique programmable ou FPGA. 9 - Emulateur selon la revendication 7 caractérisé en ce qu'il comporte un 15 processeur de signal numérique de type DSP.