Brevet FR3131669A1 — publié le 20230707

Procédé de génération d’un signal de test pour capteur RF représentatif du signal provenant d’au moins un émetteur et reçu par un récepteur RF comportant les étapes suivantes : - échantillonnage du signal émis puis pour chaque échantillon : - calcul de son instant d’arrivée au récepteur - calcul des coordonnées de positions et vitesses du récepteur à l’instant d’arrivée du signal - calcul de l’amplitude du signal reçu en fonction des positions de l’émetteur et du récepteur puis à générer le signal de test à partir des valeurs calculées par : - échantillonnage des instants d’arrivée à la fréquence d’échantillonnage de l’émission ; - calcul par interpolation du signal en réception et l’atténuation correspondante ; - calcul de la variation de fréquence entre l’émission et la réception; - calcul par interpolation du décalage en fréquence applicable ; - construction du signal de test à partir des valeurs précédemment calculées. Procédé de génération d’un signal de test pour capteur RF L’invention concerne le domaine des capteurs de signaux de télécommunication et, plus particulièrement, les capteurs spatiaux défilants. Dans le cadre du développement, de la qualification, voire du réétalonnage d’un capteur, notamment, dans le domaine des télécommunications, il est utile d’émettre un signal RF de télécommunication connu de référence et de le traiter avec ce capteur puis de comparer les résultats de ce traitement avec ceux obtenus au préalable, par exemple lors de test sur un capteur semblable, cette comparaison étant d’autant plus importante que la mobilité des émetteurs et des capteurs a une influence importante sur les performances des traitements de ces derniers. Il est connu de l’homme du métier de mettre en œuvre trois étapes successives pour qualifier les performances d’un capteur, à savoir : - une première étape consistant à effectuer des vérifications sur chaque signal émetteur connu, supposé fixe et pris isolément, le capteur étant lui aussi supposé fixe ; - une seconde étape consistant à effectuer des vérifications des signaux résultant d’une somme pondérée de signaux émetteurs connus supposés fixes, le capteur étant supposé fixe ; - une troisième étape consistant à effectuer des vérifications sur des signaux acquis en environnement réel, le capteur étant mobile. Dans les deux premières étapes, la mobilité n’est pas prise en compte. Dans la troisième étape, une analyse au préalable des signaux acquis est indispensable, pour s’assurer que les conditions de qualification sont réunies, c’est-à-dire que les émissions sont conformes aux émetteurs prévus au cahier des charges. Cette analyse est lourde et parfois en échec quand la densité des signaux est importante. L’objectif de l’invention est de résoudre ces inconvénients en proposant un procédé de génération d’un signal de test pour capteurs RF conforme aux revendications. Ce procédé permet de générer le signal acquis par un capteur fixe ou mobile en présence de différents émetteurs fixes ou mobiles, en décrivant les données d’entrée ainsi que les équations à utiliser. La maîtrise du signal généré est assurée au travers de la représentativité du signal d’émission de chaque émetteur impliqué, et dans la prise en compte fine, au niveau de l’échantillon, de la mobilité émetteur-capteur dans la construction du signal reçu par le capteur. Un procédé selon l’invention comporte en outre les avantages suivants : Elle permet de calculer l’influence de la mobilité entre l’émetteur et le récepteur sur un signal reçu, échantillon par échantillon, ce qui permet d’aborder l’utilisation de ces signaux particuliers dans les algorithmes de traitement du signal ; Elle s’appuie sur une représentation des signaux de départ en échantillons complexes en bande de base, et donc de pouvoir utiliser des signaux réels très représentatifs des formes d’ondes utilisées par les émetteurs ; Elle permet de calculer uniquement l’influence de la mobilité entre l’émetteur et le récepteur sur un signal reçu, ce qui délimite le calcul et permet de le maîtriser, moyennant le respect de l’hypothèse de propagation en espace libre. montre une représentation schématique en géométrie plane d’un exemple d’évolution de la position d’un récepteur mobile entre l’instant d’émission du premier échantillon d’un signal et sa réception ; présente une représentation schématique en géométrie plane d’un exemple d’évolution de la position d’un récepteur mobile entre l’instant d’émission d’un n ième échantillon d’un signal et sa réception ; illustre, en géométrie plane, un exemple d’évolution de la position d’un récepteur mobile entre l’instant d’émission d’un n ième échantillon d’un signal et sa réception en cas de variation de vitesse du récepteur entre ces instants ; et montrent respectivement un exemple du spectre amplitude /fréquence et sa correspondance en temps/fréquence d’un signal reçu sans effet Doppler ; et montrent respectivement un exemple du spectre amplitude /fréquence et sa correspondance en temps/fréquence d’un signal reçu avec effet Doppler. Nous allons maintenant décrire plus particulièrement un mode de réalisation d'un procédé de génération d’un signal de test pour capteur RF comportant une partie réelle et une partie imaginaire et étant représentatif du signal provenant d’au moins un émetteur et reçu par un récepteur RF, ce récepteur et cet au moins un émetteur étant mobiles l’un par rapport à l’autre, cette mobilité pouvant engendrer une variation de la fréquence centrale, due à l’effet Doppler, entre le signal émis et le signal reçu. Dans ce mode de réalisation et afin de simplifier la démonstration, le capteur est mobile et les émetteurs sont fixes et ont la même date de référence. De plus, les hypothèses suivantes ont été posées : La vitesse vrad de déplacement entre l’émetteur et le récepteur sur l’axe émission-réception est comptée positive dans le sens de la propagation de l’onde ; Il n’y a pas de multi-trajet (utilisation des distances directes) ; L’échantillon est modélisé par une valeur complexe du type Aej2πt ; La largeur de bande des signaux émis est supposée plus faible que la largeur de bande de réception ; Les caractéristiques de mobilité des émetteurs et du capteur sont connues et décrites dans le temps par des données de positions/vitesses ; La vitesse radiale entre l’émetteur et le récepteur est très largement inférieure à la vitesse de la lumière ; Le type d’émetteur est quelconque mais son spectre d’émission respecte l’approximation bande étroite et plus particulièrement sa largeur de bande est inférieure à un dixième de sa fréquence centrale. Dans cet exemple de réalisation, les données d’entrée, pour chacun des émetteurs, sont : Une date de référence, (jj/mm/aaaa hh:mn:ss,0) fixée arbitrairement ici au 1 er janvier 2000 ; Sa ou ses positions (E j (m)) et vitesses correspondantes (V E j (m.s -1 )) à des temps successifs j, j variant de 1 à p en fonction du temps dont l’origine est la date de référence et dans un repère donné tel par exemple le centre de la terre (Earth Centered Earth Fixed) ; Un signal d’émission de référence échantillonné, sa fréquence centrale d’émission (MHz), sa puissance isotrope rayonnée équivalente (W), sa fréquence d’échantillonnage (Hz), ainsi que le nombre d’échantillons ; La valeur du temps de décalage (s) du début de l’émission par rapport à la date de référence ; La durée de l’émission du signal de référence (s). Les données d’entrée du récepteur, appelé capteur dans la suite, sont : sa fréquence centrale de réception (MHz) ; sa date de référence (jj/mm/aaaa hh:mn:ss,0) fixé arbitrairement ici au 1 er janvier 2000 ; ses positions (R k (m)) à des temps successifs k, k variant de 1 à q au cours du temps dont l’origine est la date de référence et dans un repère donné tel par exemple le centre de la terre (Earth Centered Earth Fixed) ; ses vitesses (V Rk (m.s-1)) aux temps successifs k, k variant de 1 à q dont l’origine est la date de référence et dans un repère donné tel par exemple le centre de la terre (Earth Centered Earth Fixed) ; Par ailleurs, la différence de vitesse, notamment radiale, entre l’émetteur et le récepteur, est modélisée, sur chaque échantillon du signal bande de base émis, par : un retard, conséquence de la prise en compte de la vitesse de propagation des ondes, constante et égale à la célérité de la lumière dans le vide, notée c ( c = 299 792 458 m.s-1 ) et de la vitesse radiale de déplacement v rad du capteur par rapport aux émetteurs ; un décalage en fréquence (effet Doppler) entre la fréquence d’émission f eme et la fréquence de réception f rec . Elle est calculée selon la formule classique sans prendre en compte la relativité : (1) un affaiblissement de l’amplitude de l’échantillon calculé en considérant une propagation en espace libre et en champ lointain. De plus, les hypothèses suivantes ont été posées : Entre deux échantillons de positions/vitesses, l’émetteur et le récepteur sont considérés comme subissant un déplacement rectiligne uniforme ce qui permet d’appliquer de la trigonométrie plane) ; La période entre deux échantillons de positions/vitesses est identique pour l’émetteur et le récepteur afin de permettre la synchronisation des cinématiques émetteur et récepteur ; La période entre deux échantillons de positions/vitesses est supérieure au temps qui sépare la réception de deux échantillons de signal successifs pour ne pas avoir de variation de vitesse systématique entre deux échantillons. Le signal de référence choisi est un signal représentatif d’un signal radiofréquence de communication vocale et du type d’émetteur et de la forme de l’onde émise. Afin d’avoir des positions et vitesses des émetteurs à des dates identiques à celles du récepteur, des tableaux des positions/vitesses des émetteurs sont construits en interpolant les valeurs initiales (E j , V Ej ) sur la datation des positions/vitesses (R k , V Rk ) des récepteurs. Cette procédure permet de synchroniser les positions et vitesses des émetteurs et du récepteur dans le temps. Au niveau de chaque émission arrivant sur le récepteur, l’effet Doppler aura pour effet de modifier progressivement la fréquence centrale du signal émis. La largeur de bande de réception qu’il faut prendre en compte pendant toute la durée de l’émission peut être significativement plus grande. C’est pourquoi cet accroissement de largeur de bande est calculé à partir la formule (1) où la vitesse radiale est égale à : = . La nouvelle largeur de bande est égale à (3) Le signal de chaque émetteur est rééchantillonné à une fréquence d’échantillonnage f2 telle que cette dernière soit plus grande que la nouvelle largeur de bande calculée et soit un sous multiple de 2 de la fréquence d’échantillonnage en réception afin de minimiser les effets de bords du suréchantillonnage. La présente un émetteur fixe occupant une position E1 et un récepteur mobile occupant des positions successives Rk, k ∈ [1,…,K] dans le temps. Le signal émis est quant à lui échantillonné à la fréquence d’échantillonnage f2 soit avec une période tech et possède N échantillons numérotés de 1 à N. A l’instant tech1E, l’émetteur E est à la position E1 et le premier échantillon ech1 du signal d’amplitude A ech1E , est émis par E mais toujours présent en E 1 tandis que le récepteur R est supposé être à la position R 1 connue et n’a pas encore reçu le premier échantillon ech1 du signal. Dans le cas où le récepteur ne se trouverait pas à une position connue à l’instant où le premier échantillon est émis, par exemple s’il se trouvait entre deux positions connues R k et R k+1, il serait nécessaire de calculer sa position par interpolation à partir de ces valeurs. A l’instant tech1R où le récepteur reçoit le premier échantillon ech1 du signal émis avec une amplitude atténuée A ech1R , l’émetteur E est toujours à la position E1 tandis que le récepteur, mobile, est à la position R tech1R . La première étape du procédé consiste à calculer, au regard de la , l’instant t ech1R où le récepteur R reçois le premier échantillon ech1 du signal émis par l’émetteur E sachant que (4) où d ech1R divisé par c est le décalage temporel entre l’instant d’émission du premier échantillon du signal par l’émetteur E et l’instant où le récepteur, mobile, le reçoit. Pour calculer la distance d ech1R , on utilise la règle du cosinus dans un triangle quelconque : (5) Comme la position de l’émetteur et du récepteur, lorsque ce dernier se trouve en R1, sont connues, la distance d1 peut être calculée par la formule suivante : (6) La distance d R1, Rtech1 est fonction de t ech1R et s’exprime par : (7) où les vitesses VX R1 , VY R1 , VZ R1 sont connues. Le cosinus de l’angle α est connu et s’exprime par le produit scalaire suivant : (8) En reportant (5) dans (4), on obtient : (9) En reportant (7) dans (9), on obtient : (10) En élevant au carré de chaque côté : Ce qui donne : (12) L’équation (12) est une équation du second degré à une inconnue de type : (13) Avec : (14) Les deux solutions sont : (15) Pour éliminer une des solutions, on utilise le fait que le retard tech1R – tech1E est positif dans tous les cas donc, notamment, dans le cas simple où tech1E est nul et le cosinus de l’angle α est nul (angle droit). Dans ce cas simple, on obtient : (16) Avec : (17) (18) (19) On a au numérateur : (20) Dans ce cas le numérateur est négatif et comme le dénominateur est positif, tech1R serait négatif, ce qui n’est pas possible. La solution est donc : (21) Avec : (14) Et : (22) (23) Une seconde étape consiste à calculer l’amplitude A tech1R du premier échantillon quand il arrive au niveau du récepteur R au temps t ech1R . Il a parcouru d’après l’équation (3) la distance : (24) L’affaiblissement de propagation sur cette distance est obtenu par la formule : (25) Où f est la fréquence centrale du signal qui porte l’échantillon. L’amplitude reçue est donc : (26) (27) Une troisième étape consiste à calculer les coordonnées de la position R tech1R du récepteur quand il reçoit le premier échantillon. Elles sont les suivantes : (28) Si la distance entre R 1 et R tech1R est plus grande que celle entre R 1 et R 2 , alors la position calculée R tech1R est au-delà de R 2 et le calcul de t ech1R doit être repris en remplaçant R 1 et R 2 par R 2 et R 3 et ainsi de suite jusqu’à ce que la distance entre R k et R tech1R soit plus petite que celle entre R k et R k+1 . Par conséquent, k est le plus petit entier positif ou nul tel que : (29) Quand l’inégalité est vérifiée, la vitesse du récepteur, lorsqu’il reçoit le premier échantillon, a l’expression suivante : Lorsque la vitesse du récepteur ne change pas durant le trajet d’un échantillon echn , la trajectoire du récepteur est rectiligne entre la position RtechnE du récepteur quand l’échantillon echn est en E1 et la position RtechnR quand l’échantillon echn arrive sur le récepteur. La configuration est alors celle présentée sur la . Pour généraliser l’exposé précédant au cas du nième échantillon, il suffit de reprendre les équations (21) et (27) en les généralisant soit : (30) (31) Avec : (32) Pour l’échantillon n , la position du récepteur R est celle à la position R technE au moment de l’échantillon n en E décalée de la distance parcourue selon la vitesse en Rk sur une durée égale à la différence entre t echnR et t echnE . (33) La distance d n à l’instant technE est : (34) Le cosinus devient : (35) Dans notre cas de figure, l’inégalité suivante est vraie : (36) Le principe de calcul consiste donc à calculer systématiquement R technR pour chaque échantillon echn successif et pour un premier couple de valeur k et k+1 . Pour chaque valeur de R technR , on conserve cette valeur si l’inégalité est vérifiée, et on passe au calcul du R technR suivant. Si l’inégalité (36) n’est pas vérifiée, cela signifie que la vitesse radiale du récepteur a varié durant le trajet de l’échantillon echn du signal émis. Dans ce cas, l’évolution du système est présentée sur la . L’équation 4 se réécrit de la façon suivante : (37) La distance l n est connu et s’exprime par : (38) En constatant que la distance dRk+1,RtechnR est fonction de la vitesse de R (connue) et de la différence du temps global technR - technE et du temps que R a mis pour parcourir dRtechnE, Rk+1 (connue), la distance dRk+1,RtechnR s’exprime par : (39) Où les vitesses VX Rk , VY Rk , VZ Rk , VX Rk+1 , VY Rk+1 , VZ Rk+1 , sont connues. Le cosinus de l’angle α est connu et s’exprime par le produit scalaire suivant : (40) Ainsi, l’équation générale (36) s’écrit alors : (41) L’équation (41) est une équation du second degré à une inconnue de type : (42) Avec : (43) La solution est donc : (44) (45) Pour l’échantillon n , la position du récepteur R est déduite de l’équation (39) : (46) Lorsque ce calcul est terminé, le procédé selon l’invention consiste à calculer la position R technR du récepteur pour l’échantillon suivant en incrémentant k de 1 et à partir des équations (21) et (31). Dans la suite, est présenté une autre variante de réalisation de l’invention dans laquelle l’émetteur et le récepteur sont mobiles. Pour cela, dans les équations (21) et (31), il suffit de remplacer la position E 1 par la position En de l’émetteur à l’instant t echnE correspondant au nième échantillon émis. Les équations (31) deviennent : (47) La distance d n devient : (48) Avec : (49) L’équation (35) est modifiée comme suit : (50) Si l’inégalité (36) n’est pas vérifiée, cela signifie que la vitesse radiale du récepteur a varié durant le trajet de l’échantillon echn du signal émis et tandis que A et C restent inchangés, la valeur de B dans les équations (43) devient : (51) Avec : (52) et (53) Que le ou les émetteurs soient fixes ou mobiles et que la vitesse radiale entre émetteur et récepteur soit constante ou varie, comme explicité précédemment, l’étape suivante d’un procédé selon l’invention consiste à construire le signal théoriquement reçu par le récepteur à partir du signal émis par l’émetteur et des calculs de décalage, de variation de fréquence et d’atténuation calculés précédemment. Pour cela, pour chaque couple émetteur-récepteur, le traitement est identique. Dans cet exemple de réalisation, le signal émis est de forme complexe en bande de base, fonction du temps et s’écrit : (54) Et l’atténuation attdeb a été considéré comme le rapport entre l’amplitude reçue et émise de l’échantillon : (55) Un observateur qui se situe au niveau du récepteur observe le signal dans une fenêtre temporelle qui commence à l’instant du premier échantillon présent au niveau du récepteur et qui dure la durée des instants t echnR qui peut être légèrement plus grande que la durée du signal émission si le récepteur s’éloigne de l’émetteur. En supposons que l’observateur échantillonne le signal à la même fréquence que la fréquence d’émission. Les instants des échantillons du signal retardé correspondent à l’interpolation des temps retardés ( t echnR ) par une séquence échantillonnée régulièrement à la fréquence f ech (inverse de t ech ) sur la durée totale de t echnR . (56) Le signal émis retardé est l’interpolation du signal émis au temps tempsecheemerec . (57) L’atténuation finale à prendre en compte est aussi l’interpolation de l’atténuation initiale au temps tempsecheemerec (58) Le décalage en fréquence est calculé selon les équations (1) et (2). L’équation (1) est obtenue à partir des équations (19), (46) et (47). Pour calculer v rad selon (2), les positions/vitesses du récepteur sont d’abord interpolées aux instants technR . Comme le signal est en bande de base, le décalage en fréquence est un décalage relatif et a pour expression : (59) Une fois le décalage en fréquence Δ f rec calculé, le décalage en fréquence définitif à la réception est l’interpolation de la fréquence Δ f rec au temps tempsecheemerec . (60) Finalement le signal reçu en bande de base se construit par le calcul suivant : (61) La mise en œuvre du procédé selon l’invention est réalisée au moyen d’un logiciel et d’un ordinateur de type « personal computer ». Un premier sous-programme lit le fichier signal de la forme d’onde d’un émetteur, l’échantillonne à une fréquence d’échantillonnage respectant l’équation 3. Les échantillons d’émission sont sauvegardés. Un second sous-programme calcule : les instants d’arrivée au récepteur de ces échantillons du signal émis en fonction des positions et vitesses respectives de l’émetteur et du récepteur ; l’atténuation du signal reçu à ces instants d’arrivée. Un troisième sous-programme : Rééchantillonne les instants d’arrivée à la fréquence d’échantillonnage de l’émission ; Calcule, par interpolation, à partir des instants d’arrivée rééchantillonnés le signal en réception, et l’atténuation correspondante ; Calcule, pour les instants d’arrivée, la variation de fréquence entre l’émission et la réception due à l’effet Doppler en fonction de la vitesse relative radiale entre l’émetteur aux instants de début de l’émission et le récepteur aux instants d’arrivée ; Calcule par interpolation, à partir des instants d’arrivée rééchantillonnés, le décalage en fréquence qui sera appliqué au signal de réception ; Génère le signal de test à partir de cette variation de fréquence et des valeurs interpolées du retard de réception et de l’atténuation. Des tests ont été réalisés à partir d’un récepteur et trois émetteurs émettant dans des lieux différents. Ces émetteurs émettent des signaux différents à savoir une forme d’onde NXDN, une forme d’onde TETRA et une forme d’onde APCO 25. Ces trois émissions, provenant de trois lieux différents, ont été regroupés en un seul signal bande de base en les plaçant à des fréquences centrales proches. Trois niveaux de vérification ont été réalisés : Avant construction des signaux reçus, avec un outil d’analyse technique pour examiner l’occupation temps-fréquence ; Après construction des signaux reçus, avec un outil d’analyse technique pour examiner l’occupation temps-fréquence ; Après exploitation d’un des signaux reçus par un outil de décodage pour vérifier les propriétés de ce signal notamment l’intelligibilité de la voix. Avant que le signal reçu ne soit calculé, il a été enregistré au format ‘.wav’ et visualisé avec un outil d’analyse de signaux. Un exemple de signal reçu en supposant l’émetteur et le récepteur fixe, donc sans effet Doppler, est fourni sur la . Sur la figure temps/ fréquence de la correspondant au spectre de la . En plaçant une ligne verticale 10 partant du centre de la fréquence de réception, on constate que la fréquence centrale ne varie pas dans le temps, donc qu’il n’y a pas de variation de fréquence due à l’effet Doppler. La montre un exemple d’un spectre de signal reçu lorsque le récepteur est mobile par rapport à l’émetteur. Sur la figure temps/fréquence de la correspondant au spectre de la , en plaçant une ligne verticale 11 partant du centre de la fréquence de réception, on constate une fréquence centrale plus élevée que la fréquence de réception au départ de l’acquisition et qui décroit progressivement. Ce décalage de fréquence est généré par l’effet Doppler et due au déplacement radial relatif entre l’émetteur et le récepteur. On remarque aussi sur les courbes que les conséquences de l’effet Doppler ne sont pas mises en évidence sur les représentations du spectre que le signal soit court ou long et/ou en forme de chute d’eau. Pour la vérification des bonnes propriétés des signaux reçus, les signaux sont enregistrés au niveau du récepteur en bande étroite de façon unitaire puis ensemble en bande large. Les signaux sont enregistrés sans effet Doppler et avec effet Doppler. La vérification a lieu avec un logiciel de décodage sur le signal APCO 25 qui dure le plus longtemps (13 secondes). Ce signal contient un message vocal de lecture des nombres de 1 à 10. Cet outil de décodage selon l’état de l’art, part du signal complexe bande de base, reconstruit et effectue plusieurs étapes, notamment de synchronisation, démodulation et décision pour reconstruire le signal. Le fichier ‘.wav’ au format 16 bits de l’émission est lu par l’outil et le message vocal est clairement reproduit. Le fichier ‘.wav’ en réception fixe en bande étroite est lu par l’outil WAVECOM et le message vocal est reproduit avec une qualité moindre mais intelligible. Cette diminution de la qualité est liée à la diminution du rapport signal à bruit causée par l’affaiblissement de propagation et l’accroissement de la bande d’échantillonnage. Le fichier ‘.wav’ en réception mobile en bande étroite est lu par l’outil et le message vocal est reproduit correctement uniquement sur la fin du message. Cette reproduction partielle est liée au fait que le logiciel de l’outil ne peut pas suivre la variation de fréquence centrale totale sur toute la durée du signal. Pour s’assurer que c’est bien un problème dû à l’outil de production et non pas à un défaut de propriétés du signal, plusieurs signaux sont construits à partir du signal reçu pour plusieurs décalages de fréquence Doppler (500, 700 et 1000 Hz). On vérifie alors que la production correcte est plus longue et décalée dans le signal. Les propriétés intrinsèques liées au protocole sont donc bien conservées dans le signal. La méthode est donc validée. Procédé de génération d’un signal de test pour capteur RF comportant une partie réelle et une partie imaginaire et étant représentatif du signal provenant d’au moins un émetteur et reçu par un récepteur RF, ce récepteur et cet au moins un émetteur étant mobiles l’un par rapport à l’autre et leurs positions et leurs vitesses étant connues en fonction du temps et les signaux émis par ledit au moins un émetteur étant eux aussi connus et comportant une partie réelle et une partie imaginaire, procédé comportant les étapes suivantes : - Une étape d’échantillonnage du signal émis par ledit au moins un émetteur à une fréquence fech ; puis, pour chaque échantillon de ce signal : - Une étape de calcul de l’instant d’arrivée au récepteur de cet échantillon du signal émis en fonction des positions et vitesses correspondantes de l’émetteur et du récepteur ; - Une étape de calcul des coordonnées de positions et vitesses du récepteur à l’instant d’arrivée au récepteur de l’échantillon du signal émis; - Une étape de calcul de l’amplitude du signal reçu en fonction des positions de l’émetteur à la date de l’émission de cet échantillon et du récepteur à sa date de réception ; puis à générer le signal de test à partir des valeurs précédemment calculée selon : - Une étape d’échantillonnage des instants d’arrivée à la fréquence d’échantillonnage de l’émission ; - Une étape consistant à calculer par interpolation, à partir des instants d’arrivée rééchantillonnés, le signal en réception et l’atténuation correspondante ; - Une étape de calcul, pour les instants d’arrivée, de la variation de fréquence entre l’émission et la réception due à l’effet Doppler en fonction de la vitesse relative radiale entre l’émetteur aux instants de début de l’émission et le récepteur aux instants d’arrivée; - Une étape consistant à calculer par interpolation, à partir des instants d’arrivée rééchantillonnés, le décalage en fréquence qui sera appliqué au signal de réception ; - Une étape de construction du signal de test à partir de cette variation de fréquence et des valeurs interpolées du retard de réception et de l’atténuation. Procédé selon la revendication 1, caractérisé en ce que, pour chaque échantillon de ce signal, il comporte une étape de calcul des positions et vitesses de l’émetteur et du récepteur à partir des valeurs connues de positions et vitesses. Procédé selon l’une quelconque des revendications 2 et 3, caractérisé en ce qu’il comporte une étape de calcul de la largeur de bande de réception en fonction des vitesses radiales minimale et maximale entre l’émetteur et le récepteur, ainsi que de la fréquence centrale et de la largeur de bande du signal émis. Procédé selon la revendication 3, caractérisé en ce que la largeur de bande de réception est calculée à partir de la formule suivante : Où : vradmin est la vitesse radiale minimale entre l’émetteur et le récepteur vradmax est la vitesse radiale maximale entre l’émetteur et le récepteur f eme est la fréquence centrale du signal émis lb eme est la largeur de bande du signal émis. Procédé selon l’une quelconque des revendications 3 et 4, caractérisé en ce qu’il comporte une étape supplémentaire de rééchantillonnage du signal d’émission à une fréquence f2 telle que cette dernière soit plus grande que la nouvelle largeur de bande calculée et soit un sous multiple de 2 de la fréquence d’échantillonnage en réception. Procédé selon l’une quelconque des revendications 1 à 5, caractérisé en ce que le signal émis échantillonné est de forme complexe en bande de base et en ce que le signal généré de réception échantillonné S RRF est du type : S RRF (ti) = A F (ti) ∗ S E (ti)∗exp (2iπ∗ti∗Frecfin) Où : ti est la date d’émission par l’émetteur de l’échelon i à laquelle est ajoutée le temps de retard t echnR d’arrivée au récepteur du premier échantillon du signal émis, i variant de 1 à m, m étant le nombre total d’échantillons du signal émis, A F est l’atténuation du signal calculée au temps ti S E :correspond au signal complexe reçu émis au temps ti Frecfin est la variation de fréquence due à l’effet Doppler au temps ti. Procédé selon l’une quelconque des revendications 1 à 6 caractérisé en ce que pour le n ième échantillon du signal émis, le calcul de l’instant t echnR de début de l’arrivée au récepteur de cet échantillon et l’amplitude correspondante du signal reçu sont effectuées à partir des équations respectives suivantes : (21) (31) Avec : (32) Procédé selon la revendication 7, caractérisé en ce qu’il comporte une étape de calcul de la position du récepteur lorsqu’il commence à recevoir le n ième échantillon du signal émis à partir de la formule suivante : Procédé selon la revendication 8, caractérisé en ce qu’il comporte une étape supplémentaire de vérification de l’exactitude de la position calculée du récepteur lorsqu’il commence à recevoir le n ième échantillon du signal émis, cette étape consistant à vérifier l’inégalité suivante : puis, lorsque l’inégalité n’est pas vérifiée, à recommencer les étapes des revendications 7 et 8 en remplaçant Rk+1 par Rk+2, voire Rk+3, et ce, jusqu’à ce que l’inégalité soit vérifiée.