FR 3052615 A1 20171215 FR 1655322 A 20160609 1 Procédé de démodulation d'un signal reçu, produit programme d'ordinateur et dispositif correspondants. 1 DOMAINE TECHNIQUE Le domaine de l'invention est celui de la transmission de données via un lien radiofréquence basé sur la modulation d'une forme d'onde dite « chirp » telle qu'utilisée dans la technologie LoRa®. Plus précisément, l'invention concerne une méthode de démodulation d'une telle forme d'onde, qui présente des performances améliorées par rapport aux techniques existantes avec une complexité d'implémentation comparable. La technologie LoRa® étant dédiée à la transmission basse consommation par les objets connectés, l'invention a des applications dans tous les domaines de la vie personnelle et professionnelle dans lesquels les objets connectés sont présents, notamment, mais non exclusivement, dans les domaines de la santé, du sport, des applications domestiques (sécurité, électroménager, etc.), suivi d'objets, etc. 2 ARRIERE-PLAN TECHNOLOGIQUE Présentés comme la « troisième révolution de l'Internet », les objets connectés sont en train de s'imposer dans tous les domaines de la vie quotidienne et de l'entreprise. La plupart de ces objets sont destinés à produire des données grâce à leurs capteurs intégrés afin de fournir des services à valeur ajoutée pour leur propriétaire. De par les applications visées, ces objets connectés sont pour la plupart nomades. En particulier, ils doivent pouvoir transmettre les données produites, régulièrement ou à la demande, à un utilisateur déporté. Pour ce faire, la transmission radio longue portée du type radio mobile cellulaire (2G/3G/4G...) a été une technologie de choix. Cette technologie permettait en effet de bénéficier d'une bonne couverture réseau dans la plupart des pays. Cependant, l'aspect nomade de ces objets s'accompagne souvent d'un besoin en autonomie d'énergie. Or, même basés sur une des technologies radio mobile cellulaire les plus économes en énergie, les objets connectés actuels continuent de présenter une consommation rédhibitoire pour permettre un déploiement à grande échelle à un coût raisonnable. 3052615 2 Face à la problématique de la consommation du lien radio pour de telles applications nomades, de nouvelles technologies radio basse consommation et bas débit dédiées spécifiquement aux réseaux « Internet des Objets », c'est-à-dire des technologies radio pour des réseaux dits LPWAN (pour « Low-Power Wide-Area 5 Networks » en anglais), sont en train de voir le jour. Dans ce contexte, deux sortes de technologies peuvent être distinguées : d'un côté, il existe des technologies propriétaires comme par exemple la technologie de la société Sigfox®, ou bien la technologie LoRa®, ou encore la technologie de la société Qowisio®. En pratique, ces technologies non standardisées reposent toutes 10 sur l'utilisation de la bande de fréquences « Industriel, Scientifique et Médical », dite ISM, et sur la réglementation associée à son utilisation. L'intérêt de ces technologies est qu'elles sont déjà disponibles et permettent le déploiement rapide de réseaux sur la base d'investissements limités. En outre, ils permettent le développement d'objets connectés très économes en énergie et à faible coût ; 15 d'un autre côté, il existe plusieurs technologies promues par des organismes de normalisation. A titre d'exemple, on peut citer trois technologies en cours de standardisation auprès du 3GPP (pour « 3rd Generation Partnership Project » en anglais) : NB-loT (pour « Narrow Band - Internet of Things » en anglais), LTE MTC (pour « Long Term Evolution - Machine Type Communication » en anglais) et EC- 20 GPRS (pour « Extended Coverage - General Packet Radio Service » en anglais). Cependant, de telles solutions ne sont pas encore entièrement spécifiées et vont par ailleurs reposer sur des bandes de fréquences licenciées. Dans ce contexte, il apparaît que les technologies propriétaires basées sur l'utilisation de la bande ISM se présentent comme les solutions de choix sur le court 25 terme et l'une ou certaines d'entre elles peuvent alors s'imposer de facto comme la solution à utiliser. Par exemple, le brevet EP 2 449 690 B1 décrit une technique de transmission de l'information basée sur la modulation d'un signal chirp de base sur laquelle se base la technologie LoRa®. 3052615 3 Or, certains opérateurs, comme Bouygues® ou Orange® en France, se sont déjà portés sur la technologie LoRa® pour déployer leur réseau dédié aux objets connectés. Cependant, les premiers retours relèvent des expériences utilisateur peu satisfaisantes liées à des performances limitées du lien radio en conditions réelles. 5 Il existe donc un besoin pour améliorer les performances d'un récepteur implémentant la technologie LoRa® en conditions réelles, et notamment face à un canal de propagation radio mobile présentant des évanouissements. Il existe également un besoin qu'une telle amélioration n'induise pas une surconsommation en énergie du récepteur et ne pénalise donc pas l'autonomie de 10 l'objet connecté embarquant un tel récepteur. 3 RESUME Dans un mode de réalisation de l'invention, il est proposé un procédé de démodulation d'un signal reçu. Ce signal reçu résulte de la modulation d'un signal chirp de base, dont une fréquence instantanée varie linéairement entre une première 15 fréquence instantanée f0 et une deuxième fréquence instantanée f1 pendant un temps symbole Ts, et de la transmission du signal chirp modulé dans un canal de transmission. La modulation correspond, pour un symbole de rang s d'une constellation de N symboles, s étant un entier de 0 à N-1, à une permutation circulaire du motif de variation de ladite fréquence instantanée sur le temps symbole Ts, obtenue par un 20 décalage temporel de s fois une durée temporelle élémentaire Tc, telle que N*Tc=Ts. Un tel procédé comprend une étape d'estimation d'un symbole porté par le signal reçu, mettant en oeuvre les sous-étapes suivantes : détermination de N composantes de décision, à partir du signal reçu et d'un signal chirp de référence obtenu en modulant le signal chirp de base par un symbole de 25 référence correspondant à un symbole de rang r dans la constellation, une composante de décision d'indice I, notée composante D1, étant fonction d'un terme dont la phase dépend quadratiquement de I, avec I un entier de 0 à N-1 ; décision du rang k du symbole porté par le signal reçu, à partir de la composante de décision, d'indice k, notée composante Dk, présentant un extremum de valeur 30 parmi les N composantes de décision. 3052615 4 Ainsi, l'invention propose une solution nouvelle et inventive pour permettre l'estimation d'un symbole porté par un signal reçu résultant de la modulation d'un signal chirp de base présentant une variation linéaire de sa fréquence instantanée, ou de manière équivalente, présentant une variation quadratique de sa phase instantanée. 5 Pour ce faire, le procédé revendiqué propose de tenir compte de cette variation quadratique de la phase instantanée du signal reçu afin de mettre en oeuvre un récepteur optimal pour décider le rang du symbole reçu. Les performances de réception sont ainsi améliorées tout en conservant une complexité comparable à celle des récepteurs de l'art antérieur. 10 Selon un mode de réalisation, l'étape d'estimation d'un symbole comprend en outre les étapes suivantes, pour N échantillons du signal reçu et pour N échantillons du signal chirp de référence, pris aux mêmes instants multiples de Tc : conjugaison des N échantillons du signal chirp de référence, respectivement des N échantillons du signal reçu, délivrant N échantillons d'un signal chirp conjugué ; 15 multiplication terme à terme des N échantillons du signal chirp conjugué par les N échantillons du signal reçu, respectivement du signal chirp de référence, délivrant N échantillons d'un signal multiplié ; transformation de Fourier directe ou inverse du signal multiplié, délivrant N échantillons Y1 d'un signal transformé, avec I un entier allant de 0 à N-1 ; 20 et la composante Dk est en outre fonction d'un terme proportionnel à une amplitude de l'échantillon d'indice k, Yk, parmi les N échantillons Y1 du signal transformé, ainsi que de la phase de l'échantillon Yk. Ainsi, le procédé revendiqué propose de tenir compte de l'information complète (i.e. amplitude et phase) contenue dans les échantillons du signal en sortie de 25 transformée de Fourier, directe ou inverse, et non pas sur la seule base du module de ces échantillons comme fait dans l'état de la technique. Les performances de réception sont ainsi améliorées tout en conservant une complexité comparable. Selon un mode de réalisation, la composante Dk est en outre fonction d'un sous-ensemble de N' échantillons Yr, parmi les N échantillons Y1 du signal transformé, avec n 30 différent de ak, avec N' N, et avec a un paramètre appartenant à {-1,1}. 3052615 5 Ainsi, la méthode revendiquée permet de prendre en compte la dispersion du canal et l'interférence entre symboles qui en résulte pour décider du rang du symbole reçu, améliorant par là-même les performances de réception en présence d'un canal de transmission présentant des trajets multiples. 5 Selon un mode de réalisation, le procédé comprend une étape d'obtention (E45) de N coefficients de canal et un échantillon d'indice n du sous-ensemble d'échantillons Yn est pondéré par un coefficient de couplage proportionnel au coefficient de canal Hak-n[N] dépendant de la différence entre les indices ak et n, avec a un paramètre appartenant à {-1,1}, et à un terme dont l'argument dépend quadratiquement de 10 l'indice k. Le terme proportionnel à une amplitude de l'échantillon Yk est un coefficient de canal Ho indépendant de k. Ainsi, les termes pondérant les échantillons Yn ont une composante dépendant uniquement de l'écart entre les indices de ces échantillons considérés en sortie de transformée de Fourier. En effet, l'invariance dans le temps de la réponse impulsionnelle 15 du canal conduit à des termes représentatifs de l'interférence entre symboles dépendant uniquement de l'écart entre les indices des échantillons de signal considérés. Cependant, la variation quadratique de la phase du signal reçu impose que le couplage entre échantillons ne soit pas invariant dans le temps pour un écart donné entre indices d'échantillons considérés. 20 Ainsi, la prise en compte de ces deux effets dans la structure même d'une composante Dk considérée pour estimer le symbole reçu permet d'implémenter une réception aux performances améliorées en présence d'un canal de transmission présentant des trajets multiples, tout en permettant de travailler dans le domaine fréquentiel, i.e. en travaillant sur les échantillons en sortie de transformée de Fourier. 25 Selon différents modes de réalisation, la composante Dk est fonction d'un terme proportionnel à : rn k(o-k-n) la partie réelle de la somme vN v* ,2trriTi e2in N Hak-n[N]Sk, ou du complexe conjugué de la somme, lorsque la transformation de Fourier est une transformée de Fourier directe et lorsque le signal chirp conjugué correspond à la 30 conjugaison du signal chirp de référence ; ou 3052615 6 la partie réelle de la somme EiV, yn*e2in e-2in rn k(o-k-n) N Hak_n[N]Sk, ou du complexe conjugué de la somme, lorsque la transformation de Fourier est une transformée de Fourier inverse et lorsque le signal chirp conjugué correspond à la conjugaison du signal chirp de référence ; ou rn k(ak-n) 5 la partie réelle de la somme EnN=1, Yii*e 2in e2in N rLio.* k-n[NiSicK , ou du complexe conjugué de la somme, lorsque la transformation de Fourier est une transformée de Fourier directe et lorsque le signal chirp conjugué correspond à la conjugaison du signal reçu ; ou rn k(ak-n) la partie réelle de la somme Lio- YN -irri * =1v*e2k-n[N]Sk* ou du 10 complexe conjugué de la somme, lorsque la transformation de Fourier est une transformée de Fourier inverse et lorsque le signal chirp conjugué correspond à la conjugaison du signal reçu ; k2 avec Sk = (-1)1(eJTEOE N et avec a un paramètre appartenant à {-1,1}. Ainsi, la prise en compte sous forme analytique, i.e. dans la structure même de 15 l'organe d'estimation des symboles reçus, de la forme d'onde du signal considéré, par exemple la variation quadratique de sa phase instantanée, permet une implémentation simple et efficace du récepteur optimal au sens du maximum de vraisemblance en canal de transmission multitrajet dans le domaine fréquentiel, i.e. en travaillant sur les échantillons du signal en sortie de transformée de Fourier, directe ou inverse. 20 Par ailleurs, dans une variante, seuls N' coefficients de canal sont pris en compte parmi les N possibles, simplifiant par là-même les traitements embarqués dans le récepteur. Selon un mode de réalisation, les coefficients de canal 1-16k_n[N] sont nuls pour n différent de ak. 25 Ainsi, le procédé revendiqué permet d'implémenter le récepteur optimal au sens du maximum de vraisemblance dans le domaine fréquentiel, i.e. en travaillant sur les échantillons en sortie de transformée de Fourier en présence d'un canal qui se réduit à un bruit blanc Gaussien additif, ou canal AWGN (pour « Additive white Gaussian noise » en anglais), qui n'introduit donc pas d'interférences entre symboles. Les performances 3052615 7 du récepteur sont ainsi améliorées et présentent un critère d'optimalité en canal AWGN pour un surcoût calculatoire minime. Selon un mode de réalisation, l'étape d'obtention comprend en outre une estimation des coefficients de canal à partir des N échantillons Yr, du signal transformé 5 et d'au moins un symbole ki prédéterminé. Ainsi, le procédé revendiqué permet d'estimer les paramètres nécessaires à la prise en compte du canal de transmission en vue de l'implémentation d'un récepteur optimal pour l'estimation des symboles reçus en travaillant dans le domaine fréquentiel, i.e. en travaillant sur les échantillons en sortie de transformée de Fourier. Par ailleurs, la 10 prise en compte sous forme analytique de la forme d'onde du signal considéré, par exemple la variation quadratique de sa phase instantanée, permet de n'avoir à estimer que la partie invariante dans le temps des termes générant de l'interférence entre symboles, i.e. de la partie ne dépendant que de la différence entre les indices des échantillons considérés, conduisant par là-même à une implémentation efficace de 15 l'étape d'estimation des paramètres représentatifs de l'impact du canal de transmission, et donc du récepteur. Selon un mode de réalisation, les coefficients de canal estimés formant un Ho [ vecteur fi = 111- , l'estimation des coefficients étant effectuée sur la base de Ns FIN-1 symboles reçus, ki désignant le rang du i-ème desdits Ns symboles dans la constellation 20 de N symboles, ri désignant le rang d'un symbole de référence utilisé lors de la réception dudit i-ème symbole, Yi(i) désignant N échantillons dudit signal transformé obtenu lors de ladite réception dudit i-ème symbole, le vecteur R estimé de H = [Ho s H1 , exprime comme H = -Ns 2,i=0 y i avec - _ 1 -Ns-1_,ri, HN-1 3052615 8 ri(o-o-i(ki-ri)) (Ici-r i)-0 e+2jn Ne N 17` o-o-i(ki-ri)[N] e+2jrc ri(1-o-i(ki-ri)) 2 in(ki-r (0 ri(N-1-cri(ki-ri)) (ki-ri)-(N-1 2jrc )17(0 e+2jn Ne N N N y 1-o-i(lci-ri)[1^1] lorsque la transformation de Fourier correspond à une transformée de Fourier directe et lorsque le signal chirp conjugué correspond à la conjugaison du signal chirp de référence ; ou e N ej N Yo-i(ki-ri)-0[N] ri(o-dici-ri)-o) 2 .n(lci-rilo (i) ri(o-i(ki-ri)-1) 2j (0 2j7r N e N Ycci(ki-ri)-1[N] e 2 n ri(o-i(ki-ri)-N+1) (Ici-ri)-(N-1) e2jrc N y(i) i o-i(lci-ri)-N+1[N]- lorsque la transformation de Fourier correspond à une transformée de Fourier inverse et lorsque le signal chirp conjugué correspond à la conjugaison du signal chirp de référence ; ou -1) 2 in(k i-r (0* -17'(1) = N l'cKi-ri[N] 2jrc e 2jrc -0) 2i7,(ki-ri).0 Y(0* r (o-i(ki-ri) e N lorsque la 2 ;7tri(o-i(ki-ri)-N+1) 2 jn (Ici-ri)-(N-1) i\* e e N Y` o-i(ici-ri)-N+1[N] 10 transformation de Fourier correspond à une transformée de Fourier directe et lorsque le signal chirp conjugué correspond à la conjugaison dudit signal reçu ; ou 2i7rOci-ri).0 (0* e+2jn N e N 0-6t(ci-ri)[N] i-r i)) 2 ; (ki-ri).1 (0* e N e jiT N Yi_ 0 -(ki_ri)[N] ri(N-1-cri(ki-ri)) rc 2j(ki-r i)-(N-1)* e+2jn Ne N I %,N-1-o-i(lci-ri)[N] transformation de Fourier correspond à une transformée de Fourier inverse et o. e N toc lorsque le signal chirp conjugué correspond à la conjugaison dudit signal reçu; k2 15 avec Sk = (-1)keJTEOEN et avec a un paramètre appartenant à {-1,1}. - 0i) S* ici-ri[N] lorsque la 3052615 9 Ainsi, l'estimation des paramètres nécessaires à la prise en compte du canal de transmission correspond à l'erreur quadratique minimale entre le symbole émis et le symbole reçu, réduisant par là-même les erreurs d'estimation sur le symbole reçu. Par ailleurs, dans une variante, seuls N' coefficients de canal sont pris en compte 5 parmi les N possibles, simplifiant par là-même les traitements embarqués dans le récepteur. Selon un mode de réalisation, l'étape d'estimation des coefficients de canal comprend les sous-étapes suivantes : calcul de paramètres représentatifs du coefficient de canal Ho et d'un autre des 10 coefficients de canal ; obtention de paramètres représentatifs des coefficients de canal restants à partir des paramètres calculés. Ainsi la forme d'onde chirp ainsi que le choix de la valeur de Tc dans des systèmes effectifs comme le système LoRa® (8us) qui reste grand par rapport à la 15 dispersion temporelle maximale du canal conduisent au fait que seulement deux paramètres restent à estimer (e.g. Ho et un autre terme H1 avec I non nul) pour la détermination de l'ensemble des termes H1, conduisant ainsi à une grande simplicité d'implémentation de l'étape d'estimation des paramètres représentatifs de l'impact du canal de transmission, et donc du récepteur au final. 20 Selon un mode de réalisation, le coefficient de canal d'indice I non nul est i1 inversement proportionnel à sin N. Ainsi la forme d'onde chirp et le choix de la valeur de Tc dans des systèmes effectifs comme le système LoRa® qui reste grand par rapport à la dispersion temporelle maximale du canal, conduisent également à une décroissance exponentielle de 25 l'amplitude des termes H1 en fonction de I. Ceci montre qu'il peut être envisagé de n'utiliser qu'une quantité restreinte des termes H1 pour modéliser l'effet du canal, e.g. les termes correspondant à un indice I inférieur ou égal à 10, réduisant par là-même la complexité calculatoire du récepteur optimal au sens du maximum de vraisemblance. 3052615 10 Selon un mode de réalisation, le symbole prédéterminé est un symbole d'une séquence d'apprentissage ou un symbole reçu dont le rang k a été décidé lors d'une exécution précédente de ladite étape d'estimation de symbole. Ainsi, l'estimation des paramètres nécessaires à la prise en compte du canal de 5 transmission peut se faire sur la base de symboles connus, e.g. séquences d'apprentissage ou de synchronisation, permettant par là-même une estimation robuste de ces paramètres, ou sur la base de symboles de données préalablement reçus, permettant par là-même d'affiner cette estimation au cours de la réception. L'invention concerne également un programme d'ordinateur, comprenant des 10 instructions de code de programme pour la mise en oeuvre d'un procédé de démodulation d'un signal reçu résultant de la modulation d'un signal chirp de base tel que décrit précédemment, selon l'un quelconque de ses différents modes de réalisation, lorsque ledit programme est exécuté par un processeur. Dans un autre mode de réalisation de l'invention, il est proposé un dispositif de 15 démodulation d'un signal reçu, le signal reçu résultant de la modulation d'un signal chirp de base telle que décrit précédemment. Un tel dispositif de démodulation comprend une machine de calcul reprogrammable ou une machine de calcul dédiée, apte à et configurée pour : déterminer de N composantes de décision, à partir du signal reçu et d'un signal chirp 20 de référence obtenu en modulant le signal chirp de base par un symbole de référence correspondant à un symbole de rang r dans la constellation, une composante de décision d'indice I, notée composante D1, étant fonction d'un terme dont la phase dépend quadratiquement de I, avec I un entier de 0 à N-1 ; décider du rang k du symbole porté par le signal reçu, à partir de la composante de 25 décision, d'indice k, notée composante Dk, présentant un extremum de valeur parmi les N composantes de décision. Un tel dispositif de démodulation est notamment apte à mettre en oeuvre le procédé de démodulation d'un signal reçu résultant de la modulation d'un signal chirp de base selon l'invention (selon l'un quelconque des différents modes de réalisation 30 précités). 3052615 11 Ainsi, les caractéristiques et avantages de ce dispositif sont les mêmes que ceux du procédé de démodulation décrit précédemment. Par conséquent, ils ne sont pas détaillés plus amplement. 4 LISTE DES FIGURES 5 D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante, donnée à titre d'exemple indicatif et non limitatif, et des dessins annexés, dans lesquels : la figure 1 illustre les caractéristiques d'un chirp non modulé utilisé dans la technologie LoRa® ; 10 la figure 2 illustre les fréquences instantanées et les phases instantanées de différent chirps modulés selon la technologie LoRa® ; les figures 3a et 3b illustrent des structures de réception selon différents modes de réalisation de l'invention ; la figure 4 illustre les étapes d'un procédé de démodulation selon différents 15 modes de réalisation de l'invention ; la figure 5 illustre la décroissance des termes de couplage entre échantillons selon différents modes de réalisation de l'invention ; la figure 6 illustre les performances obtenues comparativement à celles obtenues par la technique de l'art antérieur dans un mode de réalisation 20 particulier de l'invention ; les figures 7a et 7b présentent des exemples de structures de dispositif de démodulation selon différents modes de réalisation de l'invention. 5 DESCRIPTION DETAILLEE DE L'INVENTION Sur toutes les figures du présent document, les éléments et étapes identiques 25 sont désignés par une même référence. Le principe général de l'invention repose sur l'estimation d'un symbole d'un signal reçu, correspondant à un signal chirp modulé transmis dans un canal de transmission, à partir de N composantes de décision représentatives du symbole, dans une constellation de N symboles. 3052615 12 Pour ce faire, la 1-ème composante parmi les N composantes de décision est une fonction de I via un terme complexe dont l'argument varie quadratiquement en fonction de I. L'indice Ti représentatif du symbole reçu dans la constellation de N symboles est alors déterminé comme fonction de l'indice k de la composante de décision qui présente 5 un extremum de valeur parmi les N composantes de décision. La solution proposée permet notamment de démoduler un signal généré en utilisant la technique décrite dans le brevet EP 2 449 690 B1 précité. Comme déjà indiqué, ce brevet EP 2 449 690 B1 décrit une technique de transmission de l'information basée sur la modulation d'un signal chirp de base. Comme 10 représenté sur la figure 1, la fréquence instantanée 102 du signal chirp de base varie linéairement entre une première fréquence instantanée f0 et une deuxième fréquence instantanée f1 pendant la durée Ts d'un symbole. Une telle fréquence instantanée est ici représentative de la vitesse de rotation dans le plan complexe du vecteur dont les coordonnées sont données par les signaux en phase 100 et en quadrature 101 15 représentant le signal modulant destiné à moduler la porteuse radio de manière à transposer le signal chirp de base sur des fréquences porteuses et ainsi générer un signal radiofréquence. Le signal chirp étant à enveloppe constante, le signal en phase 100, respectivement en quadrature 101, oscille entre deux valeurs extrêmales, 10 et 11, 20 respectivement QO et 01, sa fréquence variant linéairement dans le temps tout comme la fréquence instantanée 102 du signal chirp de base résultant. De par la variation linéaire de la fréquence instantanée 102, le signal chirp de base ainsi défini présente une phase instantanée 103 qui varie de manière quadratique entre deux valeurs 00 et 01 pendant la durée Ts, la fréquence instantanée étant la dérivée de la phase instantanée. 25 Les signaux chirp modulés sont alors obtenus par permutation circulaire du motif de variation de la fréquence instantanée du chirp de base, sur une durée Ts, obtenue suite à un décalage temporel de k fois une durée temporelle élémentaire, dite durée « chip », Tc. L'indice k est alors représentatif du rang d'un symbole dans une constellation de N symboles et l'on a alors N*Tc=Ts. À titre illustratif, la figure 2 30 représente la fréquence instantanée 102, 102', 102", 102'" et la phase instantanée 103, 3052615 13 103', 103", 103"' de différents chirps modulés correspondant respectivement à k=0, k=1, k=2 et k=3, i.e. permettant la transmission de l'information sur la base d'une constellation de quatre symboles. Le chirp de base, correspondant à k=0, est alors interprété dans ce cas comme portant le symbole de rang zéro dans la constellation. 5 Les inventeurs ont constaté que, selon cette technologie, la détermination de la valeur d'un symbole reçu via un tel signal, i.e. la détermination de son rang k dans la constellation de N symboles, est équivalente à la détermination de l'indice k ayant servi de base au calcul du décalage temporel utilisé pour engendrer le motif de phase et de fréquence instantanées du signal chirp modulé en question. 10 Il apparaît en effet que le signal chirp de base peut s'exprimer dans le domaine temporel, et sur la durée d'une période symbole, i.e. pour t allant de 0 à Ts comme s(t) = 0(t) où - 0(0 = 27 (fo + t) t + 00 2Tfos avec 00 la valeur initiale de la phase. En pratique, le signal LoRa® est tel que la bande passante du signal chirp, i.e. 15 If1 - foi, est ajustée à l'inverse de la durée chip Tc, et f1 est choisi tel que f1 = -fo. Sachant que Ts=N*Tc, l'expression de la phase instantanée du signal chirp peut alors se réécrire comme t °(t) T oe UNTc 2)f + avec a un paramètre appartenant à (-1,1) permettant de modéliser aussi bien des chirps montants (i.e. avec une fréquence instantanée croissante), que descendants (i.e. 20 avec une fréquence instantanée décroissante). L'expression analytique, sk(t), d'un chirp modulé par un symbole de rang k dans la constellation de N symboles (k étant donc compris entre 0 et N-1), et correspondant donc à une permutation circulaire du motif du chirp de base comme décrit ci-dessus, s'exprime alors comme 25 sk(t) = s(t - kTc[Ts]) = eicp(t-kTc[Ts]) (Eq-1) où [.] désigne la fonction modulo. 3052615 14 Cette équation peut alors se reformuler comme suit, pour t compris entre 0 et Ts=N*Tc : Cr (Eq-2a) s k (t) = e E/3' (/7(''c-k)) avec : 5 .i(u) = (u - 1)u pour u E [0,1] (Eq-2b) ^Ti(u) = (u + 1)u pour u E [-1,0] (Eq-2c) On décrit maintenant, en relation avec les figures 3a et 3b deux structures de réception permettant d'estimer un symbole porté par un signal reçu, correspondant à un signal chirp de base modulé selon la technique décrite ci-dessus, i.e. permettant de 10 décider de l'indice k utilisé pour engendrer le motif de variation des fréquence et phase instantanées de ce signal, selon différents modes de réalisation de l'invention. Plus particulièrement, ces figures illustrent les structures utilisées pour effectuer les traitements sur les signaux en phase, I, et en quadrature, Q, représentant le signal modulant obtenu après démodulation radiofréquence, ou démodulation RF, du signal 15 radiofréquence reçu (dans la suite de cette demande de brevet, le terme démodulation RF désigne la transposition en bande de base du signal reçu, cette transposition délivrant des signaux I et Q analogiques représentant le signal modulant la porteuse RF reçue, et le terme démodulation désigne les traitements effectués sur les signaux I et Q, souvent après échantillonnage et quantification, conduisant à la détermination de 20 l'information contenue dans le signal modulant). Lors de cette démodulation RF, il est toujours possible de choisir une fréquence porteuse de façon à ce que f1= -fo. En pratique, de tels signaux I et Q sont obtenus via l'utilisation d'un récepteur RF connu de l'homme du métier (par exemple un récepteur à conversion directe, superhétérodyne ou toute architecture équivalente), implémentant un démodulateur RF 25 en quadrature et délivrant deux voies analogiques I et Q. Les signaux I et Q sont alors échantillonnés par un convertisseur analogique numérique ADC 301 (par exemple un convertisseur flash, ou à base de modulateur sigma-delta, ou du type SAR pour « Successive Approximation Register » en anglais, ou tout autre équivalent) présent sur la voie de réception correspondante. Dans une chaine 30 de réception classique, un tel convertisseur travaillant à une fréquence 3052615 15 d'échantillonnage souvent élevée par rapport à la bande passante du signal utile, le signal délivré par l'ADC est décimé par un étage de décimation 302 (par exemple un filtre à phase linéaire du type CIC pour « Cascaded Integrator-Comb en anglais » ou tout autre équivalent) présent sur chacune des voies I et Q de manière à délivrer chacun N 5 échantillons, pouvant être interprétés comme les parties réelles et imaginaires de N échantillons complexes. Les N échantillons complexes sont alors délivrés à un dispositif de démodulation 300, 300' comprenant différents modules. Selon le mode de réalisation illustré sur la figure 3a, les N échantillons 10 complexes sont directement délivrés à un multiplieur complexe 303. Le multiplieur complexe 303 multiplie alors terme à terme les N échantillons complexes avec N échantillons complexes représentatifs d'un signal chirp de référence conjugué délivré par un module de génération 307, ici une table de correspondance, ou LUT (pour « LookUp Table » en anglais), stockant les échantillons précalculés correspondants. 15 Un tel signal chirp conjugué est ici défini comme un signal chirp dont la fréquence instantanée varie inversement de celle du signal chirp en question. Par exemple, reconsidérant le cas d'un signal chirp de base comme discuté ci-dessus en relation avec la figure 1, i.e. dont la fréquence instantanée varie linéairement de f0 à f1 sur une durée Ts, le signal chirp de base conjugué présentera alors une fréquence 20 instantanée qui varie linéairement de f1 à f0 sur la même durée Ts. De la sorte, la multiplication d'un signal chirp par son conjugué annule la variation linéaire de sa fréquence instantanée. Le résultat présente alors une fréquence instantanée constante. Dans un autre mode de réalisation illustré sur la figure 3b, le signe de la partie imaginaire des N échantillons complexes correspondant au signal reçu est inversé par un 25 module d'inversion 310. Ainsi, le module d'inversion 310 délivre des signaux correspondant aux signaux bande de base I et Q représentatif du signal chirp conjugué du signal chirp effectivement reçu. Les N échantillons complexes ainsi obtenus sont alors délivrés au multiplieur complexe 303 qui les multiplie terme à terme avec N échantillons complexes 30 représentatifs du signal chirp de référence délivré par le module de génération 307. 3052615 16 Les N échantillons complexes délivrés par le multiplieur complexe 303 sont donc, dans ce deuxième mode de réalisation, les complexes conjugués de ceux obtenus dans le mode de réalisation décrit ci-dessus en relation avec la figure 3a. Les N échantillons complexes délivrés par le multiplieur complexe 303 sont alors 5 délivrés à un module de transformée de Fourier discrète 304. Dans un mode de réalisation, la transformée de Fourier discrète implémentée est une transformée de Fourier discrète directe. Dans un autre mode de réalisation, la transformée de Fourier discrète implémentée est une transformée de Fourier discrète inverse. 10 Ainsi, quatre modes de réalisation apparaissent ici : dans un premier mode de réalisation, la conjugaison est appliquée au signal chirp de référence (cas de la figure 3a), et la transformée de Fourier discrète implémentée est une transformée de Fourier discrète directe ; dans un deuxième mode de réalisation, la conjugaison est appliquée au signal 15 chirp de référence (cas de la figure 3a), et la transformée de Fourier discrète implémentée est une transformée de Fourier discrète inverse ; dans un troisième mode de réalisation, la conjugaison est appliquée au signal chirp reçu (cas de la figure 3b), et la transformée de Fourier discrète implémentée est une transformée de Fourier discrète directe ; 20 dans un quatrième mode de réalisation, la conjugaison est appliquée au signal chirp reçu (cas de la figure 3b), et la transformée de Fourier discrète implémentée est une transformée de Fourier discrète inverse. Dans des variantes, N s'exprime comme une puissance de deux et la transformée de Fourier discrète en question est implémentée comme une transformée 25 de Fourier rapide. Les N échantillons complexes transformés délivrés par le module de transformée de Fourier discrète 304 sont alors fournis à un module de génération 305 de N composantes de décision représentatives du rang k, dans la constellation de N symboles, du symbole porté par le signal reçu. 3052615 17 Les N composantes sont alors délivrées à un module de décision 306 qui décide le rang k du symbole reçu en fonction de l'indice de la composante qui présente un extremum de valeur parmi les N composantes. Dans une variante, les N composantes représentatives du rang k du symbole 5 modulant le signal chirp de base prennent en compte l'effet du canal de propagation. Un estimateur de canal 308 estime alors des coefficients de canal sur la base d'échantillons fournis par le module de transformée de Fourier discrète 304 et du rang du symbole reçu correspondant décidé par le module de décision 306. On décrit maintenant, en relation avec la figure 4 un procédé de démodulation 10 d'un signal reçu, permettant notamment d'estimer un symbole porté par le signal reçu selon différents modes de réalisation de l'invention. Lors d'une étape E40, un signal chirp conjugué est obtenu. Comme décrit ci-dessus en relation avec les figures 3a et 3b, ce signal chirp conjugué peut correspondre soit au signal résultant de la conjugaison du signal bande de base sr(t) représentant le 15 chirp de référence d'une durée Ts délivré par le module de génération 307 (premier et deuxième modes de réalisation précités), soit au signal résultant de la conjugaison du signal bande de base y(t) représentant le signal chirp reçu (troisième et quatrième modes de réalisation précités), également d'une durée Ts. De manière générale, le signal chirp de référence correspond à un signal chirp de 20 base modulé par un symbole de référence de rang r dans la constellation de symbole. Dans une variante, r est pris égal à 0 est le signal chirp de référence est le signal chirp de base. Lors d'une étape E41, le multiplieur complexe 303 délivre le signal multiplié au module de transformée de Fourier discrète 304. 25 Dans les premier et deuxième modes de réalisation précités, ce signal multiplié s'exprime ainsi comme y(t)s,*(t), et dans les troisième et quatrième modes de réalisation précités, ce signal multiplié s'exprime ainsi comme y*(t)sr(t), i.e. comme le complexe conjugué du signal délivré par le multiplieur complexe 303 dans les premier et deuxième modes de réalisation. 3052615 18 Une expression analytique du produit y(t)sr(t) est tout d'abord dérivée ci- après. De manière générale, le signal chirp reçu a été propagé via un canal de propagation radioélectrique dont la réponse impulsionnelle h(t) peut s'exprimer 5 classiquement comme une somme de P trajets décalés dans le temps, chaque trajet pouvant alors être modélisé par une amplitude complexe Ap et un délai réel Tp de sorte que h(t) = EpP:,1)Ap 8(t - Tp) (Eq-3) avec 8(t) la distribution de Dirac. 10 Par ailleurs, le signal reçu est également entaché d'un bruit additif w(t) supposé Gaussien et centré de sorte que l'on peut écrire de manière générale que : y(t) = (h * s k)(t) + w(t) avec t E [O, TS T + - -max] et Tmax = TP-1, le support de la réponse impulsionnelle h(t) étant [0, Tmaxl. Une fois le récepteur synchronisé en temps, on peut alors écrire, en supposant 15 que le signal reçu correspond à un signal chirp de base modulé par un symbole de rang k dans la constellation de symboles, que +00 P-1 y(t) = h(r)sk (t - r)dr + w(t) = Ap sk(t - Tp) w(t) 0 P=° Ainsi, en sortie du multiplieur complexe 303 et dans les premier et deuxième modes de réalisation précités, il apparaît que Lors d'une étape E42, une transformée de Fourier est appliquée par le module 20 de transformée de Fourier discrète 304 afin de délivrer un signal transformé. Afin de simplifier les écritures, la suite du calcul est présentée pour le cas particulier où le symbole de référence correspond au chirp de base, i.e. pour r=0, quand bien même les résultats seront donnés pour le cas général. N T En notant u = L et u0 = N - + - et en définissant E comme T Ts 0 si (u - u0) E [0,1] (1 si (u - u0) E [-1,0[ (y (04 (t) = J h(r)sk(t - r)dr sr* (t) + w (t)sr* (t) 0 E = 3052615 19 on peut alors utiliser l'expression de sk(t) donnée par (Eq-2a) pour exprimer sk(t - r)s*(t) comme = e 2 - u 0) - c7) (u)] sk(t - 'Os* (t) = e2jrc-26[(u-uo+e)(u-u0+e-1)-(u-1)il = e 2 jrci 26[-2uu0 +2eu+4-2euo +uo +s2 ejrc1^lo-(4+uo)e2j7rN o-((s-uo)u-suo) Par application d'une transformée de Fourier discrète directe (DFT) sur le signal échantillonné ukn(r) = sk(nTc - 'Os* (nTc), il apparaît que DFT({ukn(r)} Uk/(T) 2 (_1)k ei (k+*) ejrc o-* e 2jrc o-N((s-uo)k-suo)e _2'.'N n=0 En notant q le terme e 2jN .7T ( que = e2j 0-reL N-1 6V'+'71.c)±'), il apparaît 2jrcNo-((s-u0)7711-suo)e ln Tc n=0 e -2_11'17 (0-(k+*)+1)7,4 N-1 + e-zirc(0-(k+*)+1» n=0 n=k+1 k N-1 = cr-Tc lqn)+ n=k+1 qn n=0 e e-Tc (1 qk+1) qk+1(1 te-k-1) qk+1 1-a (1 - e-2i"*) (1 - q) Ainsi sin (Tra 7T) = N (-1)k ei eek+r'c)2 (e(k±*)+1)(k Cette équation peut être reformulée de manière à faire apparaître les termes qui dépendent du canal de propagation et ceux liés à la forme d'onde utilisée. Ainsi N sin [714(o- (k + + 1)1 3052615 20 ic(o-k-F0 . re (o-k+1) N (-1yrk+1 e I N Ukler) Ne 'in sin (7 (o-k + / + o-TT )) ejirt e(GL)2-*)(-1)kei"ç N sin [kr (o-k + 1 + o- TT )] On peut alors finalement exprimer les échantillons du signal transformé comme P-1 DFT ((y (nTc.)s*(nTc.)))(1) = Yl = Ap Uki(rp) + DFT ((lm (nTc.)s* (nTc.))) P=0 ou sous une autre forme , i. k(o-k-F0 Yl = Ne 'n N I I o-k+I[N]Sk avec I et k de 0 à N-1 et / fit (up \ 2 -cp F11= (-1)1e FT EP-1 A e N TiON -P 1) (Eq-4a) P0 P Tc 5 ON (X) = sin(irx) (Eq-4b) N sin (le) (Eq-4c) Sk (-1)k ejne; (Eq-4d) W1 = DFT fin, (nTc)s*(nTc)} Dans le cas général où le signal chirp de référence correspond à un signal chirp de base modulé par un symbole de référence de rang r dans la constellation de 10 symboles, le calcul donne pour les N échantillons du signal transformé Y1 obtenus en sortie du module de transformée de Fourier 304 : dans le premier mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y(nTc.).sr* (nTc.) et sur w(nTc.).sr* (nTc.)) : . rl . (k-r)(o-(k-r)+1) Yl = Nee-2-InNe 2. (kN 116(k-r)+1[1\1]S k-r[N] + Wl(Ecl-Sa) 15 dans le deuxième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y(nTc.)sr*(nTc.) et sur w(nTc.)sr*(nTc.)) : . rl (k-r)(o-(k-r)-1) Yl = e+2Jn Ne 2 j reN H 0-(k-r)-1[1\1]S k-r[N] + Wl(Ec1-513) dans le troisième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y* (nTc.)sr(nTc.) et sur w* (nTc.)sr(nTc.)) : rl . (k-r)(o-(k-r)-1) 20 = N e±211T Ho-* (k-r)-1[1\1]Sk* -KM+ (Ecl-Sc) dans le quatrième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y* (nTc.)sr(nTc.) et sur w* (nTc.)sr(nTc.)) : rl . (k-r)(o-(k-r)+1) = e+2PriTre+21' Ho-* (k-r)+1[N]Sk* -r[N] + (51-5d) 3052615 21 Par ailleurs, afin de simplifier la lecture, les mêmes notations Y1, I/1 et W1 sont utilisées pour désigner les échantillons correspondants obtenus en sortie du module de transformée de Fourier 304 quel que soit le mode de réalisation précité considéré. Lors d'une étape E43, N composantes de décision ,01, I entier de 0 à N-1, pouvant 5 être interprétées comme représentant les N composantes d'un vecteur de décision et représentatives du rang du symbole porté par le signal reçu sont déterminées par un module de génération 305. Pour ce faire, il est proposé dans un mode de réalisation d'appliquer un critère de maximum de vraisemblance sur les N échantillons 1/1 délivrés par le module de 10 transformée de Fourier discrète 304. En effet, l'hypothèse Gaussienne pour le bruit additif w(nT,) reste vraie pour les échantillons W1 obtenus en sortie du module de transformée de Fourier discrète 304, la transformation de Fourier d'une distribution Gaussienne donnant une autre distribution Gaussienne. En reconsidérant par exemple le premier mode de réalisation précité 15 (correspondant à l'application d'une transformée de Fourier directe sur y(nTc)sr(nT,)), et en reconsidérant le cas particulier où le symbole de référence correspond au signal chirp de base, i.e. pour r=0, pour plus de clarté dans les écritures, les échantillons W1 peuvent s'exprimer comme suit sur la base de l'équation (Eq-5a) : . k(o-k-F0 Wl = - Ne 21' N H 6k+1[N]Sk Ainsi, appliquant un critère du maximum de vraisemblance, le rang du symbole 20 modulant le signal chirp de base et correspondant au signal reçu, correspond à l'indice k maximisant la densité de probabilité du symbole observé en réception, ou, s'agissant d'une densité Gaussienne, à l'indice k minimisant l'argument de la fonction Gaussienne, i.e. la quantité N-1 2 k(o-k+n) Yn - Ne PT N H6k+n[N]Skf n=o De manière équivalente, après développement du module au carré et 25 changement de variable de n vers N-n, il apparaît que le rang du symbole correspondant au signal reçu s'exprime en fonction de l'indice k maximisant la quantité (N 2 ic(o-k-n) -N* -ne -te N lio-k-n[N]Sk n=1 3052615 22 où 9(.) désigne la partie réelle. De manière équivalente, le complexe conjugué de l'argument de la partie réelle ci-dessus pourrait être pris. En d'autres termes, N composantes de décision DI, I de 0 à N-1, permettant l'estimation du rang du symbole porté par le signal reçu peuvent être déterminées sur la 5 base de cette expression prise pour les différentes hypothèses de rang de symbole possibles (i.e. les N hypothèses correspondant à k allant de 0 à N-1 dans l'expression ci dessus). Chacune des N composantes de décision DI correspond alors à la quantité ci-dessus prise pour l'hypothèse de rang de symbole correspondant, et l'estimé Ti du rang du symbole porté par le signal reçu s'exprime alors en fonction de la composante de 10 décision, d'indice k, notée composante Dk, présentant un extremum de valeur parmi les N composantes de décision DI ainsi déterminées. Dans le cas général où le signal chirp de référence correspond à un signal chirp de base modulé par un symbole de référence de rang r dans la constellation de symboles, un calcul équivalent permet de définir les N composantes de décision DI 15 obtenues en sortie du module de génération 305, la composante de décision d'indice k, Dk, s'exprimant comme: dans le premier mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y(nTc.).sr* (nTc.) et sur w(nTc.).sr* (nTc.)) : rn . k(crk-n) Dk = 9(EnN =1YN*-nejlT N e 4IT N I I o-k-n[N]Sk) (Eq-6a) 20 dans le deuxième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y(nTc.).sr* (nTc.) et sur w(nTc.).sr* (nTc.)) : ,k(crk-n) Dk = 9(EnN =P-Yn*ejlT N e 4IT N H o-k-n[N]Sk) (Eq-6b) dans le troisième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y* (nTc.)sr(nTc.) et sur w* (nTc.)sr(nTc.)) : . rn . k(ak-n) 25 Dk = 9(EnN =1:11- Yn*e N ejlT N lio-*k-n[N]Sk*) (Eq-6c) dans le quatrième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y* (nTc.)sr(nTc.) et sur w* (nTc.)sr(nTc.)) : rn . k(ak-n) Dk = (EnN =1YN* -ne 4IT N ejlT N Hu* k-n[N]Sk*) (Eq-6d) 3052615 23 Comme discuté ci-dessus, dans des variantes, c'est le complexe conjugué de l'argument de la partie réelle définissant Dk qui est pris dans les équations (Eq-6a) à (Eq6d). Dans une variante, le canal de propagation radioélectrique se réduit à un trajet 5 unique (e.g. dans le cas d'une liaison point à point en vue directe). Dans ce cas, la réponse impulsionnelle donnée par l'équation (Eq-3) se réduit à un seul terme d'amplitude Ao. De même, dans l'hypothèse d'une synchronisation parfaite du récepteur, on a ro = O. Il apparaît alors à partir des équations (Eq-4a) et (Eq-4b) que tous les termes H/ sont nuls pour I allant de 1 à N-1, et que seul Ho est non nul. 10 Ainsi, dans ce cas particulier où le canal de propagation se réduit à un canal AWGN (pour « Additif White Gaussian Noise » en anglais), les N composantes de décision DI obtenues en sortie du module de génération 305 et données dans le cas général par les équations (Eq-6a) à (Eq-6d) se simplifient et la composante de décision d'indice k, Dk, s'exprime comme : 15 dans le premier mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y (nT,) sr* (nT,) et sur w(nTc)sr* (nT,)) : ( Dk = 9-t YN* _o_k[N]e2inerT i Ho S k k (ni-7a) dans le deuxième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y(nTc)sr* (nT,) et sur w(nTc)sr* (nT,)) : k 20 Dk = 9-t(Y0.*k[N]e2iner7HOSk (Eq-7b) dans le troisième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y* (nT,)sr(nT,) et sur w* (nTc)sr(nT,)) : Dk = 91 (r;k[N] e 217r0-7 I; f ms . i; (Eq_70 dans le quatrième mode de réalisation précité (correspondant à l'application d'une 25 transformée de Fourier inverse sur y* (nTc)sr(nT,) et sur w* (nTc)sr(nT,)) : k Dk = 91 YN*o_ _k[N]e ( 2irc6Tr ims I*, (Eq_7d) Comme discuté ci-dessus, dans des variantes, c'est le complexe conjugué de l'argument de la partie réelle définissant Dk qui est pris dans les équations (Eq-7a) à (Eq7d). 3052615 24 Il apparaît ainsi sur les équations (Eq-7a) à (Eq-7d) que le récepteur optimal en canal AWGN au sens du maximum de vraisemblance appliqué aux échantillons pris en sortie de transformée de Fourier, directe ou inverse, fait intervenir un terme Sk (dont l'expression est donnée par l'équation (Eq-4c)) dont la phase varie quadratiquement en 5 fonction de l'indice de l'échantillon considéré dans les composantes de décision Dk permettant d'estimer le symbole reçu. Cette variation quadratique est directement liée à la variation quadratique de la phase instantanée du signal reçu. La prise en compte de la loi de variation particulière de cette phase instantanée permet ainsi d'implémenter le récepteur optimal au sens du 10 maximum de vraisemblance à un coût analytique comparable à celui lié au récepteur de l'état de la technique qui base la décision uniquement sur le module des échantillons en sortie de transformée de Fourier comme décrit dans le document de brevet EP 2 449 690 B1. Il apparaît également dans ce cas que le seul coefficient lié au canal de 15 propagation présent dans les équations (Eq-7a) à (Eq-7d), i.e. le coefficient Ho, se réduit à une constante de normalisation indépendante de l'indice k. Cependant, il apparaît que la phase de ce terme Ho (phase liée au temps de propagation subi par le signal reçu depuis son émission) se retrouve sommée avec la phase d'autres termes dépendant de k dans l'argument de la fonction partie réelle apparaissant dans les équations (Eq-7a) à 20 (Eq-7d). Ainsi, bien qu'indépendant de k, le terme Ho impacte tout de même l'indice k correspondant à la composante de décision Dk présentant un extremum de valeur parmi les N composantes de décision. Par ailleurs, reconsidérant les équations (Eq-6a) à (Eq-6d), il apparaît maintenant pour un canal présentant des trajets multiples que les termes de couplage de Dk 25 pondérant les échantillons Yn, pour n différent de k, sont proportionnels à un coefficient de canal Hvk-n[N] dépendant uniquement de l'écart entre les indices des échantillons de signal considérés en sortie de transformée de Fourier, directe ou inverse. En effet, l'invariance dans le temps de la réponse impulsionnelle du canal conduit à des termes représentatifs de l'interférence entre symboles dépendant uniquement de l'écart entre 30 les indices des échantillons de signal considérés. 3052615 25 Cependant, la variation quadratique de la phase du signal reçu impose que le couplage entre échantillons ne soit pas invariant dans le temps pour un écart donné entre indices d'échantillons considérés. Plus particulièrement, le terme Sk dont la phase varie quadratiquement en fonction de l'indice de l'échantillon considéré, et qui est 5 intrinsèquement lié à la structure même de la forme d'onde utilisée, est ici également présent. Ainsi, la prise en compte de ces deux effets dans la structure même des N composantes de décision utilisées pour estimer le symbole reçu permet d'implémenter un récepteur au sens du maximum de vraisemblance en présence d'un canal de 10 propagation présentant des trajets multiples, tout en permettant de travailler dans le domaine fréquentiel, i.e. en travaillant sur les échantillons en sortie de transformée de Fourier, directe ou inverse. Lors d'une étape E44, un estimé k du rang k du symbole porté par le signal reçu est décidé à partir de l'indice de la composante de décision Dk qui présente un 15 extremum de valeur parmi les N composantes déterminées lors de l'étape E43. Plus particulièrement, l'estimé k correspond à k = r + arg miax{Dk} [N] La combinaison des étapes E43 et E44 permet alors de mettre en oeuvre une étape E46 d'estimation du symbole reçu. Il apparaît à la lumière des expressions des composantes de décision Dk 20 données par les (Eq-6a) à (Eq-6d) ou (Eq-7a) à (Eq-7d), que dans certains modes de réalisation, les coefficients de canal H1, I allant de 0 à N-1, doivent être connus pour la mise en oeuvre de l'étape de décision E44 . Dans un mode de réalisation, les coefficients de canal H/ sont initialisés à une valeur par défaut, e.g. Ho est mis à 1 et les coefficients de canal H1, I allant de 1 à N-1, 25 sont mis à 0 pour permettre l'initialisation de la réception. Ainsi, la réception de premiers symboles peut avoir lieu et une obtention des coefficients de canal H1, I allant de 0 à N-1, peut alors être réalisée, comme décrit ci-dessous en lien avec l'étape E45, pour une mise en oeuvre ultérieure de l'étape de décision E44. 3052615 26 Lors d'une étape E45, les N coefficients de canal H1, I allant de 0 à N-1 sont ainsi obtenus. Dans un mode de réalisation, les caractéristiques du canal de propagation sont connues (e.g. dans une configuration statique) et les N coefficients de canal obtenus 5 correspondent alors à N coefficients de canal prédéterminés qui peuvent être directement chargés à l'initialisation dans le module de décision 306. Dans un autre mode de réalisation, les caractéristiques du canal de propagation sont inconnues à l'avance (e.g. en cas de mobilité du récepteur et/ou de l'émetteur) et les N coefficients de canal obtenus correspondent à N coefficients de canal R1 estimés 10 lors d'une sous-étape E451. Plus particulièrement, la méthode décrite base cette estimation sur les échantillons délivrés par le module de transformée de Fourier discrète 304 lors d'une mise en oeuvre préalable des étapes E40 à E42 ainsi que sur le rang d'au moins un symbole prédéterminé correspondant. 15 Dans une variante, les symboles prédéterminés en question sont des symboles d'une séquence d'apprentissage (e.g. un préambule ou une séquence d'apprentissage d'une trame radio), permettant par là-même une estimation robuste des coefficients de canal. Dans le cas d'une transmission LoRa®, il s'agit alors d'une pluralité de signaux chirp de base, i.e. correspondant à un symbole de rang 0 dans la constellation, de pente 20 positive ou négative (i.e. la valeur de a varie entre +1 et -1 d'un chirp à l'autre). Dans une autre variante, les symboles prédéterminés en question sont des symboles de données dont le rang a été préalablement déterminé lors de l'exécution d'une étape E44 précédente, permettant par là-même d'affiner l'estimation des coefficients de canal au cours de la réception. 25 Dans un mode de réalisation, cette estimation est effectuée sur un seul symbole reçu afin de simplifier cette étape d'estimation et de réduire la consommation globale de l'objet connecté embarquant la technique décrite. Dans un autre mode de réalisation, cette estimation est effectuée sur la base d'une pluralité de symboles reçus, permettant par là-même de moyenner l'estimation 30 afin de réduire sa variance. 3052615 27 De manière générale, si l'on considère Ns symboles pour estimer les N coefficients de canal H1, I allant de 0 à N-1, en notant ki le rang du i-ème de ces Ns symboles dans la constellation de N symboles, et ri le rang du symbole de référence utilisé lors de la réception de ce i-ème symbole, les équations (Eq-5a) à (Eq-5d) nous 5 donnent l'expression des N échantillons du signal transformé Yi(`), I allant de 0 à N-1, obtenus en sortie du module de transformée de Fourier 304 dans les quatre modes de réalisation précités lors de la réception de ce i-ème symbole. Par manipulation algébrique, il est possible d'isoler les N coefficients de canal H/ dans ces équations. Ainsi, adoptant une notation vectorielle pour plus de clarté et 10 notant H le vecteur dont les composantes sont les N coefficients du canal H1, on peut écrire à partir des équations (Eq-5a) à (Eq-5d) que Y'(i) = H + W'(i) (Eq-8) avec H0 H = 111 HN-1 et avec les composantes du vecteur 0i) données par : 15 Dans le premier mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y (nTc)sr* (nT,) et sur w(nTc)s* (nT,)) par : +2jrcri(o-o-i(ki-ri)) (i) e Ne N 170- 0 e+2jrc N e2jrc N 1-o-i(lci-ri)[N] rip-t-cri(ki-ri)) 2 .7r(ki-ri)-(N-1) (i) e +2 jrc N e YN-1-o-i(lci-ri)[N]_ 0i) = 1S,* r N (Eq-9a) Dans le deuxième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier inverse sur y(nTc)si.K (nT,) et sur w(nT,)4 (nT,)) par : 2i7Tri(0-i(kri)-0) y(i) N e N ui(ki-ri)-0[N] 2inrderi(ki-ri)-1) (i) 20 0i) = S* N e N yei(ki-ri)-1[N] (Eq-9b) 2jrcri(o-i(lci-r)-N+1) e2j _ .7r(1ci-ri)(N-1) (i) Yo-i(lci-ri)-N+1[N] 3052615 28 Dans le troisième mode de réalisation précité (correspondant à l'application d'une transformée de Fourier directe sur y*(nTc)sr(nT,) et sur w*(nTc)sr(nT,)) par : 2i7Tri(0-i(kri)-0) 2i7,(Ici-ri).0 y(0* N e ei(ki-ri)-0[N] 2inrderi(ki-ri)-1) 2i7r(ki-riI1 (0* 0i) = 1S* r N e N Yo-i(ki-ri)-1[N] (Eq-9c) N ki-riLN1 2jrcri(o-i(ki-ri)-N+1) zIn ._(;ci-ri)-(N-1) e Y") cri(ki-ri)-N+1[N]- Dans le quatrième mode de réalisation précité (correspondant à l'application d'une 5 transformée de Fourier inverse sur y*(nTc)sr(nT,) et sur w*(nTc)sr(nTc)) par : ri(0-0-i(ki-ri)) 2i7,(ki-ri).o (0* e+2jrc N e N YO-6i(ki-ri)[1\1] ri(t-o-i(ki-ri)) 2 ,n(ki-r).1 (0* 0i) = S* e+217r N e i N y (Eq-9d) kt-ri[N] 1-o-i(ki-ri)[N] ri(N-1-cri(ki-ri)) rr Oci-ri).(N-1)y(i)* e+2jrc e2j _ et avec W'(i) un vecteur dont la 1-ème composantes est proportionnelle à l'échantillon W1 obtenus en sortie du module de transformée de Fourier 304 lors de la réception du ième symbole utilisé pour l'estimation du canal. Il apparaît ainsi que W'(i) est un vecteur 10 Gaussien blanc et centré. Le vecteur H peut alors être estimé sur la base d'un critère de maximum de vraisemblance. La densité de probabilité du vecteur W'(i) étant Gaussienne, le vecteur Ro R = fit estimé de H maximisant la densité de probabilité du symbole observé en [ fiN-1 réception, sachant qu'un symbole de rang k a été émis, correspond au vecteur R 15 minimisant l'argument de la fonction Gaussienne, i.e. la quantité Ns-1 Dr° - i=o où I1.11 désigne la norme Hermitienne. Après développement du carré de cette norme, il apparaît que R s'exprime comme la moyenne sur les Ns symboles considérés des vecteurs 0i), i.e. =Ns 1= 0 Y'(`) (Eq-9e) 3052615 29 avec le vecteur donné par les équations (Eq-9a) à (Eq-9d) suivant le mode de réalisation précité considéré. On décrit maintenant, en relation avec la figure 5 une simplification dans l'estimation des paramètres de canal selon un mode de réalisation de l'invention. 5 Plus particulièrement, reconsidérant les équations (Eq-4a) et (Eq-4b), il apparaît que les variations de l'argument de la fonction ON(.), i.e. 6 Tc + 1, restent faibles autour de I entier. En effet, dans la technologie LoRa®, Tc est choisit égal à 8i.ts, ce qui reste faible devant la dispersion observée dans la plupart des canaux de propagation radioélectrique connus (i.e. devant les écarts entre les délais Tp, p non nul, associés à 10 chaque trajet au-delà du délai de trajet principal, qui est souvent le trajet direct, et le délai de ce trajet principal). A titre d'exemple, les modèles de canaux de propagation en milieu urbain donnés dans le document de standardisation « 3GPP TS 45.005 V8.8.0: 3rd Generation Partnership Project; Technical Specification Group GSM/EDGE Radio Access Network; Radio transmission and reception », publié par l'ETSI en avril 2010, 15 donnent des écarts de délais, i.e. correspondant à Tij TO, inférieurs à 5i.ts. De la sorte, supposant une synchronisation parfaite du récepteur, ce qui revient à considérer ro = 0 dans les équations précédentes, un développement limité de ON(. ) autour des valeurs de son argument multiple entier conduit à pouvoir exprimer les coefficients de canal H/ pour / # 0 comme 20 re (Eq-10) 1 Ill Ioue N N sin-N avec Io un paramètre s'exprimant en fonction des paramètres du canal de propagation comme P-1 tP0.(fTp \ 2 Tp) eiN U'c) Tc P=0 Il apparaît ainsi que l'ensemble des paramètres H1, I allant de 0 à N-1, peut être déterminé sur la base de seulement deux paramètres, ce qui simplifie drastiquement 25 l'étape d'estimation de canal. Dans une variante, les deux paramètres en question sont Ho et un autre des H/ avec I différent de zéro. En effet, l'équation (Eq-10) nous montre que les paramètres H/ avec I différent de zéro peuvent être déduits de l'un d'eux. Dans cette variante, le Io = o- - P Tc 3052615 30 paramètre Ho et le paramètre H/ considéré peuvent être estimés à partir des équations (Eq-9e) et (Eq-9a) à (Eq-9d) suivant le mode de réalisation précité considéré. En effet, ces paramètres Ho et H/ sont respectivement la première et la 1-ème composante du vecteur H défini ci-dessus et peuvent ainsi être estimés suivant la technique décrite 5 pour estimer ce vecteur. Dans une autre variante, les deux paramètres en question sont Ho et le paramètre Io introduit dans l'équation (Eq-10). Le paramètre /0 peut ainsi être alternativement estimé en injectant l'équation (Eq-10) dans l'équation (Eq-8), conduisant à 01) = H +W'(1) = IoC +W'(1) 10 avec - le N N sin 2 e N - 2e N sin N Ir N -1 e N N sin 71-(N - 1) N Reprenant alors un critère de maximum de vraisemblance appliqué à cette équation pour déterminer le paramètre Io, un calcul similaire à celui décrit ci-dessus en relation avec l'obtention de l'équation (Eq-9e) donne 1 TEN.-s Cil YI") 1.0 C* Ns. LI1 avec le vecteur 0i) donné par l'une des équations (Eq-9a) à (Eq-9d) suivant le mode de 15 réalisation précité considéré, avec C*T le vecteur transposé du vecteur C*, lui-même étant obtenu en conjuguant chaque composante de C. '10 représente dans cette formule l'estimation de /0. Dans le mode de réalisation où l'estimation de canal est effectuée sur un seul symbole reçu, les équations ci-dessus restent valables en considérant Ns=1. 20 Par ailleurs, dans la variante où l'ensemble des paramètres H1, I allant de 0 à N- 1, est déterminé sur la base des deux paramètres Ho et Io, les expressions des N composantes de décision D/ obtenues en sortie du module de génération 305 et C= 3052615 31 permettant l'estimation du symbole émis, données dans le cas général par les équations (Eq-6a) à (Eq-6d), se simplifient sur la base de l'équation (Eq-10), et la composante de décision d'indice k, Dk, s'exprime comme : dans le premier mode de réalisation précité (correspondant à l'application d'une 5 transformée de Fourier directe sur y(nTc)sr* (nT,) et sur w(nTc)sr* (nT,)) : Dk = ff . rk , Io - SkHo YN* _o_k[Nie2in e r N + IN ( Ho =1 YN- n=o-k[N] n N sinE(0-k-n) e2inrNn e N j. n(2 k+1)(o-k-n)i) N (Eq-11a) dans le deuxième mode de réalisation précité (correspondant à l'application d'une 10 transformée de Fourier inverse sur y(nTc)sr* (nT,) et sur w(nTc)sr* (nT,)) : Dk = S kilo ( rk 17:-Ic[Nle2j"7 + El \ I -ni=o YrK H n=o-k[N] n N sinL(o-k-n) e2j.TcrNn e j N .n(2k+1)(o-k-n)i) N (Eq-11b) dans le troisième mode de réalisation précité (correspondant à l'application d'une 15 transformée de Fourier directe sur y* (nTc)sr(nT,) et sur w*(nTc)sr(nT,)) : Dk = rk 2jnrn jn(2k+1)(o-k-n)i) -2jrco-- Io v,1\1-1 n (ski A; ko-*k[N]e N + 7 L n=0 Yn* NsinE(o-k-n)e N e N ° n=o-k[N] N (Eq-11c) dans le quatrième mode de réalisation précité (correspondant à l'application d'une 20 transformée de Fourier inverse sur y* (nTc)sr(nT,) et sur w*(nTc)sr(nT,)) : Dk = rk 2jnrn jn(2k-n)(o-k-n)l) -2jrco-- Io v,N rc Sk* HO* [YN* -o-k[N]e N -1.. n=1 YN*-1/ N e Hi; - N sinI(o-k-n) n=o-k[N] (Eq-11d) Comme indiqué précédemment, dans des variantes, c'est le complexe conjugué 25 de l'argument de la partie réelle définissant Dk qui est pris dans les équations (Eq-11a) à (Eq-11d). Par ailleurs, il apparaît à la lumière de l'équation (Eq-10) (et donc des équations (Eq-11a) à (Eq-11d) dérivées de cette équation (Eq-10)) que l'approximation de la 3052615 32 fonction ON(.) (approximation autorisée par le choix d'une valeur du temps chip Tc. grande devant les écarts de délai associés à chaque trajet au-delà du délai principal, e.g. comme dans la technologie LoRal dans l'expression des coefficients de canal H1, pour I entier non nul, fait apparaître une variation de l'amplitude de ces termes H/ comme la 5 fonction f (l) = H représentée sur la figure 5. Il apparaît ainsi une décroissance N sin-N exponentielle de l'amplitude des coefficients H/ en fonction de l'indice I, l'amplitude du coefficient H10 étant divisée par 10 relativement à celle de H1. Par conséquent, l'effet du canal peut être correctement modélisé en ne tenant compte que d'un nombre restreint de paramètres H1, par exemple les N' premiers 10 coefficients de canal d'indice I allant de 0 à N'-1, simplifiant par là-même les traitements embarqués dans le récepteur pour la décision des symboles reçus en présence d'un canal de propagation présentant des multi-trajets. Dans une variante, les N' coefficients de canal (N' ) sont obtenus par application de la méthode générale décrite ci-dessus en relation avec les équations (Eq- 1-10 15 8) et (Eq-9a) à (Eq-9e) appliquées au vecteur H = Hl . Les vecteurs 0i) à considérer pour la mise en oeuvre de cette méthode, suivant le mode de réalisation parmi les quatre modes de réalisation précités considérés, sont alors ceux donnés par les équations (Eq-9a) à (Eq-9d), mais restreint à leurs N' premiers termes. Dans une autre variante, les N' coefficients de canal considérés sont déterminés 20 à partir de seulement deux paramètres comme décrit ci-dessus en relation avec les équations (Eq-10) et suivantes (e.g. H0 et un autre des H/ avec I différent de zéro, ou bien H0 et In). Là encore, les vecteurs considérés doivent être restreints aux N' premiers termes. Dans encore une autre variante, seuls N' coefficients de canal sont pris en 25 compte parmi les N possibles, mais ce ne sont pas les N' premiers coefficients de canal, i.e. les coefficients de canal d'indice I inférieur à N'. Dans ce cas, la méthode générale décrite ci-dessus en relation avec les équations (Eq-8) et (Eq-9a) à (Eq-9e) s'applique, mais N-N' coefficients de canal correspondant sont présupposés nuls. Il en est de même lorsque les N' coefficients de canal considérés sont déterminés à partir de seulement 3052615 33 deux paramètres comme décrit ci-dessus en relation avec les équations (Eq-10) et suivantes. Ceci permet de simplifier la structure de l'organe d'estimation du symbole reçu lorsqu'une caractéristique du canal de propagation peut être présupposée. On décrit maintenant, en relation avec la figure 6 les performances obtenues 5 lorsque la technique décrite est utilisée dans le cas d'un canal de propagation du type AWGN comparativement à celles obtenues lorsque la technique de l'art antérieur est utilisée. Dans cette situation, l'obtention des coefficients de canal selon la technique décrite effectuée à l'étape E45 se résume à l'obtention d'un seul paramètre Ho, les 10 autres termes H/ étant nuls pour I allant de 1 à N-1 comme décrit ci-dessus en relation avec la figure 4. Par ailleurs, la décision du rang des symboles reçus lors de l'étape E44 est basée dans ce cas sur l'utilisation de composantes de décision Dk données par l'équation parmi les équations (Eq-7a) à (Eq-7d) correspondant au mode de réalisation précité considéré, et déterminées lors de l'étape E43. 15 Selon la technique de l'art antérieur décrite dans le document de brevet EP 2 449 690 B1, le rang du symbole reçu est déterminé uniquement sur la base de l'échantillon en sortie de transformée de Fourier présentant l'amplitude maximale, indépendamment de toute information de phase. Il apparaît alors que l'utilisation de la technique décrite (courbe 600b), permet 20 un gain de l'ordre de 1 décibel sur le rapport Eb/NO (i.e. sur le rapport de l'énergie par bit reçu ramené à la densité spectrale de puissance de bruit) nécessaire pour obtenir un taux d'erreur binaire, ou TEB, donné par rapport à la technique connue (courbe 600a). A TEB donné, un tel gain sur le rapport Eb/NO se traduit directement sur le rapport signal à bruit requis à l'entrée du récepteur. Il en découle un gain correspondant 25 sur la portée du système global et donc sur la couverture des cellules du réseau considéré. En pratique, 1 décibel de gain sur le rapport signal à bruit à l'entrée du récepteur correspond à un accroissement conséquent de la portée de 12%. Les gains escomptés lorsque le canal de propagation présente des évanouissements sont encore supérieurs, la technique décrite permettant en effet de 3052615 34 corriger l'interférence entre symboles résultant des multi trajets et donc d'améliorer la discrimination entre le symbole émis et ses adjacents. Les figures 7a et 7b présentent des exemples de structures de dispositif de démodulation 300, 300' des symboles reçus permettant la mise en oeuvre d'un procédé 5 de démodulation décrit en relation avec la figure 4 selon différents modes de réalisation de l'invention. Le dispositif de démodulation 300, 300' comprend une mémoire vive 703, 713 (par exemple une mémoire RAM), une unité de traitement 702, 712 équipée par exemple d'un processeur, et pilotée par un programme d'ordinateur stocké dans une 10 mémoire morte 701, 711 (par 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 703, 713 avant d'être exécutées par le processeur de l'unité de traitement 702, 712. Ces figures 7a et 7b illustrent seulement une manière particulière, parmi 15 plusieurs possibles, de réaliser le dispositif de démodulation 300, 300' afin qu'il effectue certaines étapes du procédé détaillé ci-dessus en relation avec la figure 4 (dans l'un quelconque des différents modes de réalisation). 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 20 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 matériel). Dans le cas où le dispositif de démodulation 300, 300' 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 (tel 25 que par exemple une disquette, un CD-ROM ou un DVD-ROM) ou non, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. REVENDICATIONS1. Procédé de démodulation d'un signal reçu, ledit signal reçu résultant de la modulation d'un signal chirp de base dont une fréquence instantanée (102, 102', 102", 102"') varie linéairement entre une première fréquence instantanée f0 et une deuxième fréquence instantanée f1 pendant un temps symbole Ts, ladite modulation correspondant, pour un symbole de rang s d'une constellation de N symboles, s étant un entier de 0 à N-1, à une permutation circulaire du motif de variation de ladite fréquence instantanée sur ledit temps symbole Ts, obtenue par un décalage temporel de s fois une durée temporelle élémentaire Tc, telle que N*Tc=Ts, et de la transmission du signal chirp modulé dans un canal de transmission, caractérisé en ce qu'il comprend une étape d'estimation (E46) d'un symbole porté par ledit signal reçu, mettant en oeuvre les sous-étapes suivantes : détermination (E43) de N composantes de décision, à partir dudit signal reçu et d'un signal chirp de référence obtenu en modulant ledit signal chirp de base par un symbole de référence correspondant à un symbole de rang r dans ladite constellation, une composante de décision d'indice I, notée composante D1, étant fonction d'un terme dont la phase dépend quadratiquement de I, avec I un entier de 0 à N-1 ; décision (E44) du rang Ti du symbole porté par ledit signal reçu, à partir de la composante de décision, d'indice k, notée composante Dk, présentant un extremum de valeur parmi lesdites N composantes de décision. 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape d'estimation d'un symbole comprend en outre les étapes suivantes, pour N échantillons dudit signal reçu et pour N échantillons dudit signal chirp de référence, pris aux mêmes instants multiples de Tc : conjugaison (E40) desdits N échantillons dudit signal chirp de référence, respectivement desdits N échantillons dudit signal reçu, délivrant N échantillons d'un signal chirp conjugué ; multiplication (E41) terme à terme desdits N échantillons dudit signal chirp conjugué par lesdits N échantillons dudit signal reçu, respectivement dudit signal chirp de 3052615 36 référence, délivrant N échantillons d'un signal multiplié ; transformation de Fourier (E42) directe ou inverse dudit signal multiplié, délivrant N échantillons Y1 d'un signal transformé, avec I un entier allant de 0 à N-1 ; et en ce que ladite composante Dk est en outre fonction d'un terme proportionnel à 5 une amplitude de l'échantillon dudit indice k, Yk, parmi lesdits N échantillons Y1 dudit signal transformé, ainsi que de la phase dudit échantillon Yk. 3. Procédé selon la revendication 2, caractérisé en ce que ladite composante Dk est en outre fonction d'un sous-ensemble de N' échantillons Yn parmi lesdits N échantillons Y1 dudit signal transformé, avec n 10 différent de ak, avec N' N, et avec a un paramètre appartenant à {-1,1}. 4. Procédé selon la revendication 3, caractérisé en ce que le procédé comprend une étape d'obtention (E45) de N coefficients de canal, et en ce qu'un échantillon d'indice n dudit sous-ensemble d'échantillons Yn est pondéré 15 par un coefficient de couplage proportionnel au coefficient de canal Hak-nLN r dépendant de la différence entre les indices ak et n, et à un terme dont l'argument dépend quadratiquement dudit indice k, et en ce que ledit terme proportionnel à une amplitude de l'échantillon Yk est un coefficient de canal Ho indépendant de k. 20 5. Procédé selon la revendication 4, caractérisé en ce que ladite composante Dk est fonction d'un terme proportionnel à : k(crk-n) la partie réelle de la somme VN V* `.2jITT. r e-2-"T N H 0-k_n[N]Sk, ou du complexe conjugué de ladite somme, lorsque ladite transformation de Fourier est une transformée de Fourier directe et lorsque ledit signal chirp conjugué correspond 25 à la conjugaison dudit signal chirp de référence ; ou rn k(crk-n) la partie réelle de la somme EnNZ(1, Yn*e2iniTre-2-in N o-k-n[N]Sic, ou du complexe conjugué de ladite somme, lorsque ladite transformation de Fourier est une transformée de Fourier inverse et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal chirp de référence ; ou rn . k(ak-n) 30 2jrciTre2jrc N u la partie réelle de la somme EnNZ(1, Yn*e - "o-* k-n[N]Sk* ou du complexe 3052615 37 conjugué de ladite somme, lorsque ladite transformation de Fourier est une transformée de Fourier directe et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal reçu ; ou rn . k(ak-n) la partie réelle de la somme EN N- Y* e-21n-N e21re n=1 n N H*cric-n[N] ou du Sk* 5 complexe conjugué de ladite somme, lorsque ladite transformation de Fourier est une transformée de Fourier inverse et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal reçu ; k2 avec Sk = (-1)k ejne N et avec a un paramètre appartenant à {-1,1}. 6. Procédé selon la revendication 5, 10 caractérisé en ce que lesdits coefficients de canal Hak-nLN r sont nuls pour n différent de ak. 7. Procédé selon l'une quelconque des revendications 4 à 6, caractérisé en ce que l'étape d'obtention (E45) comprend en outre une estimation (E451) desdits coefficients de canal à partir desdits N échantillons Yn dudit signal 15 transformé et d'au moins un symbole k, prédéterminé. 8. Procédé selon la revendication 7, lesdits coefficients de canal estimés formant 110 un vecteur R = H1 , ladite estimation desdits coefficients étant effectuée sur la base FIN-1 de Ns symboles reçus, ki désignant le rang du i-ème desdits Ns symboles dans la constellation de N symboles, ri désignant le rang d'un symbole de référence utilisé lors 20 de la réception dudit i-ème symbole, Yi(i) désignant N échantillons dudit signal transformé obtenu lors de ladite réception dudit i-ème symbole ) caractérisé en ce que ledit vecteur fi s'exprime comme R = -N1 Ni s( 1 y f (i) avec ri(o-cri(ki-ri)) 2 ,n(ki-ri)-o (i) e+2jn Ne N Y 0-0-i(ki-ri)[N] 1+, ri(t-cri(ki-ri)) 2jrc ,,(I) e+2jn N r ri(N-1-cri(ki-ri)) (ki-ri)-(N 2jrc -1 )17(i) e +2j7[ N e N I lorsque ladite transformation de Fourier correspond à une transformée de Fourier directe et 25 lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal chirp de 3052615 38 e N e lcNl Y r do- dici- r i)-o) 2 .n(i -rio (i) o-i(ki-ri)-0[N] 2 inr do- i(k i)-1) y( uki_ri)_1[N] e 2inri(cri(ki-Nr)-N+1)e2-in(1:i-riNI(N-1)17(:)(ki-ri)-N+1[N1- lorsque ladite référence ; ou - 0i) Ski-ri[N] transformation de Fourier correspond à une transformée de Fourier inverse et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal chirp de 5 référence ; ou e 2j1Tri(cri(kiN-ri)-())e2j7r(ki-Nri).°Y(iiik o-i-ri)-0[N] 2in ri(o-i(kiN-ri)-1) e 2 j7r(ki-N eri ri).1y(ii e ki_ri)_1[N] 2inri(o-i(ki-ri)-N+1) 2 .7r (1ci-ri).(N-1-) ,,,* I N Y") e N e o-i(k-ri)-N+1[N] -17'(1) = 1+S K lorsque l'ci-ri[N] ladite transformation de Fourier correspond à une transformée de Fourier directe et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal reçu ; ou +2ri(o-o-i(ki-ri)) .7r(ki-ri).0 (0* ePT N e 21 N YO-Cri(ki-ri)[N] +2 i(1- o- i-r i)) 2i (Ic i-r i). 1 (0* e,in N enN Yl-cri(ki-ri)[N] - 0i) ski-ri[N] (ki-r%, (0* e e N I 2j+2jn rc 10 lorsque ladite transformation de Fourier correspond à une transformée de Fourier inverse et lorsque ledit signal chirp conjugué correspond à la conjugaison dudit signal reçu ; k2 avec Sk = (-1)k ejne N et avec a un paramètre appartenant à {-1,1}. 9. Procédé selon l'une quelconque des revendications 7 ou 8, 15 caractérisé en ce que ladite étape d'estimation desdits coefficients de canal comprend les sous-étapes suivantes : - calcul de paramètres représentatifs dudit coefficient de canal Ho et d'un autre desdits coefficients de canal ; - obtention de paramètres représentatifs des coefficients de canal restants à partir 20 desdits paramètres calculés. 10. Procédé selon l'une quelconque des revendications 4 à 9, 3052615 39 caractérisé en ce que ledit coefficient de canal d'indice I non nul est inversement rrl proportionnel à sin-. 11. Procédé selon l'une quelconque des revendications 7 à 10, caractérisé en ce que ledit symbole prédéterminé est un symbole d'une séquence 5 d'apprentissage ou un symbole reçu dont le rang k a été décidé lors d'une exécution précédente de ladite étape d'estimation de symbole. 12. Produit programme d'ordinateur, comprenant des instructions de code de programme pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 11, lorsque ledit programme est exécuté sur un ordinateur. 10 13. Dispositif de démodulation (300, 300') d'un signal reçu, ledit signal reçu résultant de la modulation d'un signal chirp de base dont une fréquence instantanée (102, 102', 102", 102"') varie linéairement entre une première fréquence instantanée f0 et une deuxième fréquence instantanée f1 pendant un temps symbole Ts, ladite modulation correspondant, pour un symbole de rang s d'une constellation de N 15 symboles, s étant un entier de 0 à N-1, à une permutation circulaire du motif de variation de ladite fréquence instantanée sur ledit temps symbole Ts, obtenue par un décalage temporel de s fois une durée temporelle élémentaire Tc, telle que N*Tc=Ts, et de la transmission du signal chirp modulé dans un canal de transmission, caractérisé en ce qu'il comprend une machine (702, 712) de calcul reprogrammable ou 20 une machine (702, 712) de calcul dédiée, apte à et configurée pour : déterminer de N composantes de décision, à partir dudit signal reçu et d'un signal chirp de référence obtenu en modulant ledit signal chirp de base par un symbole de référence correspondant à un symbole de rang r dans ladite constellation, une composante de décision d'indice I, notée composante D1, étant fonction d'un 25 terme dont la phase dépend quadratiquement de I, avec I un entier de 0 à N-1 ; décider du rang k du symbole porté par ledit signal reçu, à partir de la composante de décision, d'indice k, notée composante Dk, présentant un extremum de valeur parmi lesdites N composantes de décision.