FR2819963A1 - Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee - Google Patents
Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee Download PDFInfo
- Publication number
- FR2819963A1 FR2819963A1 FR0100849A FR0100849A FR2819963A1 FR 2819963 A1 FR2819963 A1 FR 2819963A1 FR 0100849 A FR0100849 A FR 0100849A FR 0100849 A FR0100849 A FR 0100849A FR 2819963 A1 FR2819963 A1 FR 2819963A1
- Authority
- FR
- France
- Prior art keywords
- state
- receiver
- continuous
- particulate
- receiver according
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/005—Control of transmission; Equalising
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0204—Channel estimation of multiple channels
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03178—Arrangements involving sequence estimation techniques
- H04L25/03331—Arrangements for the joint estimation of multiple sequences
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Récepteuroptimal particulaire délivrant de façon récurrente l'estimation conjointe à maximum de vraisemblance de l'information-source (synchronisationet symboles discret, paramètres continus) transmise par modulationpulsée, avec codes possibles d'accès et de canal, et captéesous forme numérique bruitée après le canal global incluantmulti-trajets et effets Doppler ainsi que filtre d'entrée et corrélateur par le code d'accès.Il est caractérisé par KX L X M X N processeurs numériques initiaux, réels ou virtuels,dits particulaires, recevant en parallèle les données numériqueset calculant la mesure de probabilité globale, conditionnelleà ces données jusqu'à l'instant courant, pour M séquences distinctes m parmi la combinatoire croissante des symboles libres de la source,de l'état x des paramètres continus probabilisés ensomme de N mesures conditionnelles localement multinormales indicéespar n, dites particules étendues gaussiennes. Un opérateur deredistribution à maximum de vraisemblance gère leurs morts etleur distribution conformément à leurs masses de probabilitéconditionnelle, éliminant rapidement les états peu vraisemblablesde synchronisation (k, l) durant la période d'acquisition. Le résultatdélivré comporte notamment les maxima de vraisemblance marginauxcorrespondant au message le plus vraisemblable utile en télécommunications,et les paramètres continus les plus vraisemblables utiles en télépositionnement,avec possibilité de modularité des deux fonctions et, pour ladernière nommée, discrimination possible des multi-trajetsinférieure à la période de (sur) échantillonnage.
Description
<Desc/Clms Page number 1>
Récepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux à modulation pulsée
----------------
La présente invention a pour objet de restituer conjointement de manière probabiliste optimale la totalité de l'information à la fois digitale (symboles et synchronisations) et continue (paramètres d'onde et de mouvement) transmise par signaux à modulation pulsée, avec codes possibles d'accès et de canal, reçus sous forme de mesures échantillonnées bruitées à travers un canal de transmission comportant des multi-trajets et effets Doppler.
----------------
La présente invention a pour objet de restituer conjointement de manière probabiliste optimale la totalité de l'information à la fois digitale (symboles et synchronisations) et continue (paramètres d'onde et de mouvement) transmise par signaux à modulation pulsée, avec codes possibles d'accès et de canal, reçus sous forme de mesures échantillonnées bruitées à travers un canal de transmission comportant des multi-trajets et effets Doppler.
Le récepteur suivant l'invention est fondé sur l'intégration cohérente longue que permet l'exploration et la mémorisation pertinentes, par processeurs numériques virtuels ou réels dits particulaires, de l'espace d'état global et de ses probabilités. Cela concerne plus particulièrement la probabilisation des paramètres continus sous forme de sommes de densités locales multinormales représentatives de leur espace tangent, dites particules étendues gaussiennes et sujettes à redistribution dynamique, ainsi que l'arborescence des séquences à symboles libres présentes dans le signal source, sous forme de d'états discrets en nombre réduit. Ce récepteur a la particularité de pouvoir discriminer par le calcul des effets Doppler différents ainsi que des multi-trajets dont le décalage temporel est inférieur à la période d'échantillonnage des mesures d'entrée. Sa mise en oeuvre, récurrente par son principe même, permet le fonctionnement et l'application en temps réel et assure, pour une capacité de calcul donnée (nombre de processeurs réels parallèles, ou équivalents virtuels séquentiels), l'estimation optimale de l'information-source au sens du maximum de vraisemblance.
Le secteur technique de l'invention est celui de la réception numérique de l'information transmise sous modulation pulsée, à travers le canal global de transmission incluant les filtres et les multi-trajets possibles entre émetteur et récepteur dotés de mobilité relative.
<Desc/Clms Page number 2>
Parmi les applications de cette invention figurent : - la réception optimale des messages de télécommunication, notam- ment en téléphonie cellulaire pour usage mobile tels les systèmes UMTS, EDGE, ou GSM.
- l'estimation optimale des signaux de télépositionnement, soit par voie passive satellitaire tels GPS, GNSS Galileo, ou terrestre tels LORAN-C, soit par voie aussi bien active que passive tel le RADAR, ou encore selon le mode acoustique du SONAR.
Le contexte précis est celui des signaux de transmission comportant les éléments suivants, inclus dans la figure 2 : - une pulsation porteuse No qui peut être annulée ou déplacée au moyen d'un démodulateur 2, par descente de la fréquence en bande utile dite de base.
- une enveloppe modulatrice, pouvant faire l'objet de codes emboîtés (code d'accès, code de canal) dont les impulsions respectives de largeur 8 et A sont les symboles d'un alphabet de nombres complexes ou constellation (de phase-amplitude) et forment des séquences comportant notamment les symboles libres du message utile dans le code canal.
- un bruit perturbateur additif supposé être du bruit blanc gaussien dans la bande fréquentielle de base.
- un convertisseur analogique-numérique (3/3'sur la figure 2) produisant des données indicées par leur ordre temporel t, dont le filtre d'entrée assure soit l'anti-repliement (télépositionnement) soit l'adaptation à l'impulsion élémentaire (télécommunications) en préservant la blancheur du bruit additif en temps discret.
- une déformation linéaire par superposition de multi-trajets constituée, après échantillonnage, d'une somme de retards indicés de 0 à R, ordre maximum, et pondérés chacun par leur amplitude à valeur complexe.
- une déformation non-linéaire de décalage fréquentiel due aux vitesses relatives possibles entre émetteur et récepteur ainsi qu'à l'erreur résiduelle dans la descente en fréquence, le tout constituant l'oscillation à la fréquence dite Doppler.
<Desc/Clms Page number 3>
Après passage par le corrélateur 4 avec le possible code d'entrée, les données à traiter 5, notées y, ont pour modèle :
où U désigne l'union exclusive, E la sommation numérique, j le nombre imaginaire pur (j2 -1), C, (m) est la valeur complexe du symbole à l'instant t, avant convolution par le filtre d'entrée de réponse f (. ), du message produit au moyen du code canal modulant le code d'accès Ctl, dans leurs états respectifs (k, l) de synchronisation, par une séquencesource m parmi toutes celles possibles, sous forme de sortie :
Ctk (m) =h (#tk (m))
d'une machine séquentielle à coder dont l'état interne # obéit à l'équation dynamique de transition à valeurs discrètes notée : < Pi= F ( (t, bt), d'état initial (po=k, bt étant le symbole-source à l'instant t, actualisant la séquencesource correspondante m. r est la partie entière des retards en unités d'échantillonnage #, #r,vd reste réel, vr étant l'indice de multiplicité de ces retards pour une même partie entière.
où U désigne l'union exclusive, E la sommation numérique, j le nombre imaginaire pur (j2 -1), C, (m) est la valeur complexe du symbole à l'instant t, avant convolution par le filtre d'entrée de réponse f (. ), du message produit au moyen du code canal modulant le code d'accès Ctl, dans leurs états respectifs (k, l) de synchronisation, par une séquencesource m parmi toutes celles possibles, sous forme de sortie :
Ctk (m) =h (#tk (m))
d'une machine séquentielle à coder dont l'état interne # obéit à l'équation dynamique de transition à valeurs discrètes notée : < Pi= F ( (t, bt), d'état initial (po=k, bt étant le symbole-source à l'instant t, actualisant la séquencesource correspondante m. r est la partie entière des retards en unités d'échantillonnage #, #r,vd reste réel, vr étant l'indice de multiplicité de ces retards pour une même partie entière.
Les coefficients Ard constituent les amplitudes complexes des trajets indexés de r=0 (trajet direct) à R correspondant au retard
maximum, pour une pulsation Doppler mdnumérotée d de 0 à l'entier maximum D. La pulsation Doppler md est liée à la vitesse relative vd entre émetteur et récepteur par la relation connue md = Movd/c. Les retards Tr, vd sont liés à la vitesse Vd, aux instants échantillonnés, par l'équation dynamique : Tr,vd(t) = #r,vd (t-1)+#vd (t)/c
x=(#d, At,vd, #r,vd) résume le vecteur du paramétrage continu, tel qu'utilisé par la suite et symbolisé en figure 1. v est le bruit additif à temps discret après conversion analogique-numérique, comprenant interférences et colorations possibles.
maximum, pour une pulsation Doppler mdnumérotée d de 0 à l'entier maximum D. La pulsation Doppler md est liée à la vitesse relative vd entre émetteur et récepteur par la relation connue md = Movd/c. Les retards Tr, vd sont liés à la vitesse Vd, aux instants échantillonnés, par l'équation dynamique : Tr,vd(t) = #r,vd (t-1)+#vd (t)/c
x=(#d, At,vd, #r,vd) résume le vecteur du paramétrage continu, tel qu'utilisé par la suite et symbolisé en figure 1. v est le bruit additif à temps discret après conversion analogique-numérique, comprenant interférences et colorations possibles.
<Desc/Clms Page number 4>
Trois remarques sont à noter : i) Lorsque les variables de retard Tr, v ne sont pas recherchées, comme c'est le cas en télécommunications, un changement de variable dans l'équation cohérente de y les absorbe de deux façons : soit sous la forme
(Xi, rdétant les amplitudes des réflexions retardées de r unités, soit sous la forme
aid étant les amplitudes du canal global des réflexions. ii) Lorsque le filtre d'entrée approche le passe-bas idéal, l'influence des symboles passés s'éteint lentement (en 1/t). Il en résulte un canal réverbérant de mémoire longue, rendant observables de faibles écarts entre retards, et facilitant la discrimination des multi-trajets proches du trajet direct, notamment de moins d'une période d'échantillonnage. Cet avantage, utile en télépositionnement, se manifeste dans le Jacobien de la paramétrisation non-linéaire en amplitudes et retards, qu'utilise l'estimation particulaire ci-après au moyen de particules locales Gaussiennes explorant le plan tangent associé à ce Jacobien. iii) On note que les données échantillonnées brutes, après filtre d'entrée, sont à bruit additif gaussien malgré l'interférence multi-codes d'accès en considérant ceux-ci comme relevant du théorème central limite, mais qui s'avère coloré pour ces mêmes raisons (en communications) ainsi qu'en raison du suréchantilonnage (en télépositionnement). Cette coloration du bruit de réception Vu, connue pour chacun des deux cas concernés évoqués ci-dessus, s'inclue aisément dans le modèle précédent sous la forme de la matrice de corrélation At, t, exprimant l'espérance mathématique E [vt, vt.]. Le caractère Gaussien du bruit coloré de réception, ainsi caractérisé par une augmentation paramétrique selon des procé-
(Xi, rdétant les amplitudes des réflexions retardées de r unités, soit sous la forme
aid étant les amplitudes du canal global des réflexions. ii) Lorsque le filtre d'entrée approche le passe-bas idéal, l'influence des symboles passés s'éteint lentement (en 1/t). Il en résulte un canal réverbérant de mémoire longue, rendant observables de faibles écarts entre retards, et facilitant la discrimination des multi-trajets proches du trajet direct, notamment de moins d'une période d'échantillonnage. Cet avantage, utile en télépositionnement, se manifeste dans le Jacobien de la paramétrisation non-linéaire en amplitudes et retards, qu'utilise l'estimation particulaire ci-après au moyen de particules locales Gaussiennes explorant le plan tangent associé à ce Jacobien. iii) On note que les données échantillonnées brutes, après filtre d'entrée, sont à bruit additif gaussien malgré l'interférence multi-codes d'accès en considérant ceux-ci comme relevant du théorème central limite, mais qui s'avère coloré pour ces mêmes raisons (en communications) ainsi qu'en raison du suréchantilonnage (en télépositionnement). Cette coloration du bruit de réception Vu, connue pour chacun des deux cas concernés évoqués ci-dessus, s'inclue aisément dans le modèle précédent sous la forme de la matrice de corrélation At, t, exprimant l'espérance mathématique E [vt, vt.]. Le caractère Gaussien du bruit coloré de réception, ainsi caractérisé par une augmentation paramétrique selon des procé-
<Desc/Clms Page number 5>
dés connus, permet le fonctionnement du récepteur particulaire décrit ci-après, sans autre procédure que celle d'une augmentation de dimension de son état particulaire.
Rappelons le problème posé par la réception :
Si l'on connaît l'information-source à transmettre, ainsi que les éléments constitutifs du canal de transmission (réflexions et effets Doppler inclus), on sait en déduire au moyen d'un modèle cohérent le signal qui sera capté à la réception possiblement mobile. La réception consiste à résoudre le problème inverse, qui est précisément celui que la solution de la présente invention a pour propos : Le contenu de l'information-source étant inconnu, ainsi que celui des perturbations de canal et de réception, il s'agit d'en faire à chaque instant une estimation qui rende compte au mieux des données reçues jusque là, en un sens optimal du point de vue probabiliste. Cette estimation concerne la valeur d'un jeu de variables d'état ou paramètres, de nature discrète et/ou continue, nécessaire et suffisant pour restituer l'information-source utile. Plus précisément, pour obtenir une telle solution optimale, il s'agit de construire la distribution de probabilité de ces variables, conditionnellement à l'ensemble des observations accumulées au cours du temps écoulé depuis un certain instant initial.
Si l'on connaît l'information-source à transmettre, ainsi que les éléments constitutifs du canal de transmission (réflexions et effets Doppler inclus), on sait en déduire au moyen d'un modèle cohérent le signal qui sera capté à la réception possiblement mobile. La réception consiste à résoudre le problème inverse, qui est précisément celui que la solution de la présente invention a pour propos : Le contenu de l'information-source étant inconnu, ainsi que celui des perturbations de canal et de réception, il s'agit d'en faire à chaque instant une estimation qui rende compte au mieux des données reçues jusque là, en un sens optimal du point de vue probabiliste. Cette estimation concerne la valeur d'un jeu de variables d'état ou paramètres, de nature discrète et/ou continue, nécessaire et suffisant pour restituer l'information-source utile. Plus précisément, pour obtenir une telle solution optimale, il s'agit de construire la distribution de probabilité de ces variables, conditionnellement à l'ensemble des observations accumulées au cours du temps écoulé depuis un certain instant initial.
Le problème posé, pour la situation générale que décrit le modèle ci-dessus, n'a pas de solution exacte réalisable par les procédés de réception actuellement connus, ni de solution approchée satisfaisante par ceux actuellement mis en oeuvre lorsque l'on aborde des situations critiques courantes telles que multi-trajets rapprochés et multi-Doppler. La difficulté réside à la fois dans la nature conjointe discrète-continue du problème dynamique d'estimation et dans la présence de non-linéarités sur les variables continues, conduisant à une double explosion combinatoire et dimensionnelle hors de portée des structures calculatoires envisagées jusqu'à présent.
Ainsi, pour un jeu de variables discrètes fixées, il existe bien une estimation constructive des paramètres continus au moyen de l'algorithme connu de Kalman, en utilisant une approximation linéaire, mais celle-ci présente deux écueils : d'une part, elle ne
<Desc/Clms Page number 6>
rend pas compte des non-linéarités et multi-determinations créant des optima locaux parasites qui font obstacle à l'estimation optimale ; d'autre part, son caractère optimal en dimension finie n'est plus préservé avec l'évolution conjointe des variables discrètes libres que présente le problème réel.
Inversement, pour des paramètres continus fixés, il existe bien un estimateur optimal de séquences à symboles libres du message, au moyen de l'algorithme connu de Viterbi, en utilisant le modèle à états discrets correspondant. Mais celui-ci perd sa réalisabilité en dimension finie lorsque les variables continues deviennent inconnues et doivent elles-mêmes être estimées conjointement.
Le procédé le plus proche du présent propos, et visant à résoudre le problème d'optimisation conjointe évoqué est celui connu sous le nom de PSP ("Per Survivor Processing"), faisant l'objet de publications et d'un brevet déposé aux Etats-Unis d'Amérique en 1995. Ce procédé aborde la difficulté d'estimation optimale conjointe discrète-continue décrite ci-dessus, en appliquant conjointement l'algorithme de Viterbi pour les états discrets du problème d'une part, et l'algorithme d'estimation linéaire des paramètres continus restreint aux séquences optimales retenues pour chaque état discret, d'autre part. Or, cette séparation arbitraire par juxtaposition des deux aspects du problème ne rend pas compte de leur interdépendance créant l'explosion combinatoire et dimensionnelle sus-dites. L'impossibilité d'atteindre l'optimalité par cette voie a été confirmée par preuves mathématiques dans des articles scientifiques antécédents tels que celui de Z. S. Roth et K. A. Loparo"Nonlinear filtering problems with finite dimensional algebras", paru dans la revue Systems & Control Letters, N. 7, 1986.
La solution originale au problème technique exposé ci-dessus, qui est l'objet de la présente invention, étend et complète les principes généraux d'estimation particulaire introduits par une invention antécédente du même auteur intitulée"Procédé et système d'estimation optimale non-linéaire des processus dynamiques en temps réel", dont le brevet a été déposé en France sous le numéro 94/07274, en Europe sous le numéro 9595256.5-2206, et aux Etats-
<Desc/Clms Page number 7>
Unis d'Amérique comme US Patent 5933352. Cette invention antécédente fait usage de particules ponctuelles de Dirac dont les sup- ports reproduisent l'évolution stochastique a priori des variables d'état par tirages aléatoires, pour être ensuite pondérées conditionnellement aux données suivant la formule de Bayes, et éventuellement redistribuées en conséquence.
La présente invention de récepteur particulaire comporte en premier lieu deux innovations générales, absentes de la précédente, essentielles dans les performances du dit récepteur tant en précision qu'en temps de calcul : - elle introduit pour estimer les variables continues des particules diffuses dites gaussiennes, c'est à dire des distributions à support non ponctuel, localisées par leurs moyennes et covariances. L'intérêt de sommes de telles particules, outre une capacité d'approximation universelle des mesures de probabilité, tout comme les mesures de Dirac, est de permettre cela avec un nombre réduit d'entre elles dû à leur caractère étendu et leur déplacement automatique conditionnellement aux données.
elle introduit dans l'exploration-sélection combinatoire des états discrets une exploration arborescente déterministe ainsi qu'une redistribution, également déterministe. Cette dernière est fondée sur la mémorisation d'un nombre fixé d'entre eux, ceux aux probabilités les plus élevées conditionnellement aux données recueillies, réalisant ainsi l'affectation probabiliste optimale des capacités calculatoires imposées.
elle introduit dans l'exploration-sélection combinatoire des états discrets une exploration arborescente déterministe ainsi qu'une redistribution, également déterministe. Cette dernière est fondée sur la mémorisation d'un nombre fixé d'entre eux, ceux aux probabilités les plus élevées conditionnellement aux données recueillies, réalisant ainsi l'affectation probabiliste optimale des capacités calculatoires imposées.
D'autres innovations, plus spécifiques au problème de réception, apparaissent dans l'exposé technique ci-après, et sont rappelées en revendications.
Le récepteur particulaire suivant l'invention, pour l'estimation optimale conjointe de l'information digitale et continue en modulation pulsée, intervient selon la figure 2 sur les données 5 issues du signal d'antenne 1 après la descente en fréquence de base par un démodulateur 2 et la conversion numérique assurée par l'échantillonneur 3'muni de son filtre d'entrée 3, soit adapté à la fréquence de (sur) échantillonnage (télépositionnement), soit adapté à l'impulsion élémentaire (télécommunications) et après le
<Desc/Clms Page number 8>
corrélateur 4, au code d'accès ; il comporte une unité de calcul 6 constituée de KXLXMXN processeurs initiaux 600,,,,"réel ou vir- tuels, dits particulaires, recevant en parallèle les dites données ; chacun de ces processeurs est porteur, dans les mémoires associées, d'une représentation particulaire locale nommée klmn de l'état composite (k, l, m, x) du signal à estimer, selon la structure suivante qu'illustre la figure 1 : - le couple (k, l) représente l'état possible des synchronisations de code d'accès et de code canal parmi les K XL possibles.
- l'indice entier m identifie une des séquences de symboles retenues parmi les M plus vraisemblables, capacité de mémorisation finie imposée par construction.
- l'indice entier n identifie, pour chaque triplet (k, l, m) précédent, la distribution multinormale locale Gmnkl (x) que comporte en chaque point de son espace tangent la probabilisation de l'état x des paramètres continus précédemment définis (md, Ar,vd, Tr,vd ou md, ai, rd suivant les deux cas) sous forme d'une somme globale de la forme :
avec
cette somme ayant la propriété de représentation probabiliste universellement convergente pour N
Chaque particule klmn représente donc une densité locale gaussienne mémorisée numériquement par l'état composite (k, l, m, \I), dont les deux dernières composantes sont respectivement la moyenne et la covariance de la densité Gaussienne associée et le réel scalaire p (k, l, m, n) qui représente la probabilité ou poids de la particule klmn, le tout formant l'état particulaire.
avec
cette somme ayant la propriété de représentation probabiliste universellement convergente pour N
Chaque particule klmn représente donc une densité locale gaussienne mémorisée numériquement par l'état composite (k, l, m, \I), dont les deux dernières composantes sont respectivement la moyenne et la covariance de la densité Gaussienne associée et le réel scalaire p (k, l, m, n) qui représente la probabilité ou poids de la particule klmn, le tout formant l'état particulaire.
Chacun des processeurs 600mnkl comporte un générateur 604mnkl de transition d'état des dites particules entre deux instants d'échantillonnage t-1 et t qui, tenant compte de la dynamique des symboles aussi bien imposés par le code d'accès que possibles dans
<Desc/Clms Page number 9>
le message du code canal, ainsi que de la cinématique du mouvement relatif émetteur-recepteur, délivre les nouvelles valeurs 612 de kl l'état particulaire à l'instant t. Dans chaque processeur 600mnkl, un correcteur 605mnkl calcule, à partir des dites valeurs 612mnkl et des données échantillonnées 5 les nouvelles valeurs 613 de cet état particulaire. Chaque processeur particulaire est donc caractérisé par des informations qu'il est chargé de calculer de façon récurrente et de garder en mémoire en permanence : l'état (k, 1, m, n, Pn) de la particule gaussienne klmn, ainsi que son poids p (k, l, m, n). Les composantes k, l représentent, rappelons-le, les états discrets de synchronisation de code canal et code d'accès ; la composante m désigne la branche considérée dans l'arborescence des symboles libres du message ;'les composantes et pn sont la moyenne et la covariance de la n-ième gaussienne nécessaire à la représentation de la densité de probabilité du paramètre continu x, sachant k, l, m, n. Le poids de chaque particule Gaussienne et sa localisation dans l'espace d'état reconstituent ainsi collectivement la probabilité sur tout l'espace d'état. Cette probabilité, comme annoncé au paragraphe précédent, évolue par prédictioncorrection conditionnelle aux données, ce qui est l'objet de l'estimation particulaire introduite.
La description détaillée du fonctionnement séquentiel d'un processeur particulaire 600mink1 est la suivante : - L'initialisation : A l'instant initial t=0, le commutateur 607mnkl délivre au générateur d'évolution 604mnkl les composantes 610mnkl de l'état initial de la particule klmn, fournies par le générateur d'état initial 603mnkl selon une loi de probabilité a priori représentative de l'incertitude sur tout l'espace d'état.
Concernant (k, l), les K X L états sont équiprobables et tous initialement retenus. Concernant m, l'état initial se réduit à l'alphabet des symboles libres. Concernant x, les n gaussiennes initiales doivent au moins couvrir les possibilités multi-modales dues à la non linéarité (multiplicité des vitesses Doppler aveugles, par exemple).
<Desc/Clms Page number 10>
Aux instants t suivants, le commutateur 607maki délivre au génékl rateur d'évolution 604maki les composantes de l'état particulaire 615 après redistribution à l'instant t-1. kl Le générateur d'évolution 604maki calcule alors les nouvelles composantes d'état particulaires 612maki à l'instant t conformément aux transitions d'état entre t et t-1. L'état discret (k, l) de la particule notée klmn reste invariant tant qu'il n'est pas concerné par l'organe de redistribution 606. L'état discret m de la particule klmn, qui représente la séquence m retenue depuis l'instant précédent, est momentanément augmenté des branches représentant chacune des possibilités de symbole suivant dans l'alphabet considéré, avant qu'il ne soit redistribué par 606 après la mesure.
L'état continu x évolue sur son espace tangent suivant les contraintes cinématiques linéarisées du mouvement, à faible dérive entre t et t-1, et fournit les nouvelles moyennes un et covariances
Pn locales de la particule Gaussienne klmn.
Pn locales de la particule Gaussienne klmn.
Si l'on désigne par Gmnki (x, tint-1) la distribution localement Gaussienne obtenue après la prédiction et par p (k, l, m, n, t|t-1) son poids, on peut suivre ci-dessous sa transformation après chaque nouvelle donnée 5, sous la notation conventionnelle t|t-1#t : - La correction :
A partir des composantes de l'état particulaire 612. kl et de la mesure 5 à l'instant t, le pondérateur fournit la nouvelle particule pondérée conditionnellement à cette nouvelle donnée par application de la formule de Bayes, sous la forme :
A partir des composantes de l'état particulaire 612. kl et de la mesure 5 à l'instant t, le pondérateur fournit la nouvelle particule pondérée conditionnellement à cette nouvelle donnée par application de la formule de Bayes, sous la forme :
Le caractère local des particules Gaussiennes considérées permettant la linéarisation par rapport à x de l'expression des mesures y, on obtient ainsi la multiplication de deux Gaussiennes locales qui fournit après normalisation une nouvelle particule Gaussienne locale représentée par :
<Desc/Clms Page number 11>
- La Redistribution :
Après pondération par la donnée à l'instant t, les particules Gaussiennes se trouvent toutes dotées de nouveaux poids, fonction de leurs lieux de localisation. La procédure de redistribution les utilise pour redistribuer sur ces lieux de localisation les particules renormalisées en les affectant en nombre conforme au maximum de vraisemblance.
Après pondération par la donnée à l'instant t, les particules Gaussiennes se trouvent toutes dotées de nouveaux poids, fonction de leurs lieux de localisation. La procédure de redistribution les utilise pour redistribuer sur ces lieux de localisation les particules renormalisées en les affectant en nombre conforme au maximum de vraisemblance.
Cette procédure de redistribution introduit un couplage entre tous les processeurs 600 de l'unité 6 via l'organe de redistribution. Celui-ci étant activé, il redistribue tous les états par-
kl kl ticulaires 612mu1 en 615.
kl kl ticulaires 612mu1 en 615.
Ainsi disparaissent très rapidement les états discrets de synchronisation (k, l) peu probables, abaissant leur cardinalité KXL durant l'acquisition à celle des états rémanents de probabilité p (k, l) non négligeable, comme indiqué en figure 1, libérant la capacité de calcul correspondante pour la poursuite en temps réel.
Concernant la variable d'état de séquence m, sa cardinalité momentanément augmentée durant l'étape de prédiction reprend sa valeur permanente M après redistribution, ce nombre de séquences retenues étant limité par principe à une capacité calculatoire imposée.
- L'estimation :
L'organe terminal 609 a pour fonction de délivrer l'estimation récurrente 7 à maximum de vraisemblance. Le maximum de vraisemblance conjoint correspond à l'argument (k*, 1*, m*, x*) du maximum des densités p (k, l, m, x). Les maxima de vraisemblance marginaux correspondent à : - l'argument mt du maximum de la marginale p (m), désignant le message le plus vraisemblable, à usage des télécommunications.
L'organe terminal 609 a pour fonction de délivrer l'estimation récurrente 7 à maximum de vraisemblance. Le maximum de vraisemblance conjoint correspond à l'argument (k*, 1*, m*, x*) du maximum des densités p (k, l, m, x). Les maxima de vraisemblance marginaux correspondent à : - l'argument mt du maximum de la marginale p (m), désignant le message le plus vraisemblable, à usage des télécommunications.
- l'argument xt du maximum de la marginale p (x), désignant la valeur la plus vraisemblable des variables dynamiques x, à usage du télépositionnement.
Le résultat est un nouveau type de récepteur numérique. Son procédé de calcul particulaire est adapté aux difficultés nonlinéaires et combinatoires d'estimation conjointe des variables continues et discrètes et permet une intégration cohérente préser-
<Desc/Clms Page number 12>
vant le plus longtemps possible, pour une capacité de calcul donnée, tous les états de probabilité suffisante pour contribuer à l'estimation finale, ceci pour une extraction optimale de l'information dans les situations critiques rencontrées couramment (multi-trajets rapprochés, Doppler à vitesses aveugles, faibles rapports signal/bruit).
Claims (7)
- REVENDICATIONS ------------- 1. Récepteur particulaire pour l'estimation optimale conjointe de l'information digitale (symboles, synchronisations) et continue (paramètres d'onde et de mouvement) dans les signaux (1) à modulation pulsée avec codes possibles d'accès et de canal, après démodulation (2), à travers le canal global incluant multi-trajets et effets Doppler, filtre d'entrée (3), échantillonneur (3') ainsi que corrélateur (4) par le code d'accès, et caractérisé en ce que le récepteur comporte : - K X L X M X N processeurs numériques initiaux (600mn), réels ou virtuels, dits particulaires représentent chacun un état (k, l) de la synchronisation des codes, une séquence m de la combinatoire des symboles libres, un état vectoriel réel x des paramètres continus probabilisé en somme de N mesures locales multinormales sur leur espace tangent par des particules gaussiennes klmn à distributions de probabilité respectives Gmn (x, tlt) de poids associé p (k, l, m, n, t|t), nouveau type de particules de portée générale.instant t la densité multinormale d'évolution a priori Gmnkl (x, tilt-1) de chaque particule klmn, en tenant compte de la combinatoire des séquences augmentées du symbole libre possible à l'instant t courant et de la dynamique de la variable d'état continue x.- Un générateur (604mil), initialisé par (603mnkl), calcule à chaque- Un correcteur (605mnkl) calcule à chaque instant t l'état particulaire (615=kil) de klmn à partir de sa valeur à l'instant précédent t-1 et des nouvelles données Yt (5), au moyen de la formule de Bayes utilisant la densité de probabilité du bruit additif.(600nml) redistribue tous les états particulaires (613mol) en (615nmkl), en les affectant de façon conforme au maximum de vraisem- blance.- L'organe de redistribution (606) couplé à tous les processeursL'organe terminal (609) délivre l'estimation récurrente (7) du maximum de vraisemblance conjoint p (k*, l*, m*, x*) des variables discrètes et continues, ainsi que les maxima de vraisemblance marginaux respectifs p (mit) et p (xt) du message m utile en télécommunications et des paramètres x utiles en télépositionnement.<Desc/Clms Page number 14>
- 2. Récepteur selon la revendication 1, caractérisé en ce que l'organe de redistribution (606) réduit automatiquement la cardinalité K X L des états de synchronisation (k, l) après la période initiale d'acquisition, ne laissant subsister que les seules particules de poids non négligeable.
- 3. Récepteur selon les revendications 1 et 2, caractérisé e ce que le générateur (604mol) augmente après chaque donnée la cardinalité de la variable m par les possibilités de symbole nouveau à cet instant, cette cardinalité reprenant sa valeur de consigne M après redistribution par l'organe (606), pour mémoriser un nombre constant de séquences les plus vraisemblables dans la combinatoire croissante des symboles libres.
- 4. Récepteur selon les revendications 1,2, 3, caractérisé en ce que le traitement optimal effectué repose sur la probabilité conditionnelle aux données (5) pouvant comporter un bruit d'interférence coloré (télécommunications) ou un suréchantillonnage du code d'accès qui colore le bruit additif résultant (télépositionnement).
- 5. Récepteur selon les revendications 1,2, 3,4, caracyérisé en ce qu'en télépositionnement, il utilise la réponse impulsionnelle longue du filtre d'entrée approchant le passe-bas idéal pour discriminer par le calcul les multi-trajets proches du trajet direct de moins qu'une période de (sur) échantillonnage.
- 6. Récepteur selon les revendications 1,2, 3,4, caractérisé en ce qu'en l'absence d'estimation explicite du retard (cas des télécommunications), un changement de variable fait apparaître les amplitudes de réponse échantillonnée du canal, absorbant la multiplicité des retards entiers ou non tout en laissant inchangés les autres éléments du récepteur, et permettant ainsi la modularité des deux fonctions de télécommunication et de télépositionnement.
- 7. Récepteur selon les revendications 1,2, 4, caractérisé en ce qu'en l'absence de code canal (cas du RADAR), les indices k et 1 disparaissent en donnant des particules gaussiennes bi-indicées mn.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0100849A FR2819963B1 (fr) | 2001-01-23 | 2001-01-23 | Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee |
PCT/FR2002/000284 WO2002060090A1 (fr) | 2001-01-23 | 2002-01-23 | Recepteur particulaire pour l"estimation optimale conjointe de l"information discrete et continue dans les signaux a modulation pulsee |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0100849A FR2819963B1 (fr) | 2001-01-23 | 2001-01-23 | Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2819963A1 true FR2819963A1 (fr) | 2002-07-26 |
FR2819963B1 FR2819963B1 (fr) | 2004-07-16 |
Family
ID=8859112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0100849A Expired - Fee Related FR2819963B1 (fr) | 2001-01-23 | 2001-01-23 | Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee |
Country Status (2)
Country | Link |
---|---|
FR (1) | FR2819963B1 (fr) |
WO (1) | WO2002060090A1 (fr) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110135291B (zh) * | 2019-04-29 | 2023-03-24 | 西北工业大学 | 一种低信噪比信号的参数估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5142551A (en) * | 1991-02-28 | 1992-08-25 | Motorola, Inc. | Signal weighting system for digital receiver |
US5231648A (en) * | 1991-03-21 | 1993-07-27 | Northern Telecom Limited | Adaptive equalizer for digital cellular radio |
WO1999012290A2 (fr) * | 1997-08-29 | 1999-03-11 | Nokia Networks Oy | Procede pour l'estimation de parametres et recepteur |
US5933352A (en) * | 1994-06-08 | 1999-08-03 | Salut; Gerard | Method and a system for non-linear optimal estimation of dynamic processes in real time |
-
2001
- 2001-01-23 FR FR0100849A patent/FR2819963B1/fr not_active Expired - Fee Related
-
2002
- 2002-01-23 WO PCT/FR2002/000284 patent/WO2002060090A1/fr not_active Application Discontinuation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5142551A (en) * | 1991-02-28 | 1992-08-25 | Motorola, Inc. | Signal weighting system for digital receiver |
US5231648A (en) * | 1991-03-21 | 1993-07-27 | Northern Telecom Limited | Adaptive equalizer for digital cellular radio |
US5933352A (en) * | 1994-06-08 | 1999-08-03 | Salut; Gerard | Method and a system for non-linear optimal estimation of dynamic processes in real time |
WO1999012290A2 (fr) * | 1997-08-29 | 1999-03-11 | Nokia Networks Oy | Procede pour l'estimation de parametres et recepteur |
Also Published As
Publication number | Publication date |
---|---|
FR2819963B1 (fr) | 2004-07-16 |
WO2002060090A1 (fr) | 2002-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0665665B1 (fr) | Procédé et dispositif permettant à un modem de se synchroniser sur un transmetteur de données numériques par voie hertzienne en présence de brouilleurs | |
EP0291979B1 (fr) | Procédé de démodulation de signaux modulés numériquement et dispositif de mise en oeuvre d'un tel procédé | |
EP1259833A1 (fr) | Procede de localisation de sources radioelectriques au moyen d'un radiogoniometre haute resolution deux voies | |
EP0991952B1 (fr) | Procede et dispositif de radiogoniometrie cooperative en transmission | |
FR2975193A1 (fr) | Procede et systeme de localisation d'interferences affectant un signal de radionavigation par satellite | |
FR2620583A1 (fr) | Systeme de transmission par modems comportant des canaux principal et secondaire | |
WO2012025306A1 (fr) | Dispositif spatio temporel multi-antennes multi-correlateurs pour la rejection des multi-trajets des systemes de navigation | |
EP2095150B1 (fr) | Procédé et dispositif de réception d'un signal de radionavigation à modulation boc | |
EP3543745B1 (fr) | Dispositif multi-antennes pour la rejection de multi-trajets dans un systeme de navigation par satellite et procede associe | |
EP1482648B1 (fr) | Récepteur de signal ultra large bande et procédé de réception associé | |
FR2971653A1 (fr) | Procede d'acquisition avec multi-correlation frequentielle amelioree | |
EP1961138B1 (fr) | Procede et dispositif permettant le suivi de doppler pour modem a large bande | |
FR2882479A1 (fr) | Procede et dispositif de synchronisation de liaisons rectilignes ou quasi-rectilignes en presence d'interferences | |
EP1235399B1 (fr) | Procédé et dispositif d'estimation de la matrice de covariance du canal de transmission d'un système UMTS | |
FR2819963A1 (fr) | Recepteur particulaire pour l'estimation optimale conjointe de l'information digitale et continue dans les signaux a modulation pulsee | |
FR2798794A1 (fr) | Procede de traitement de signaux de navigation a trajets multiples dans un recepteur presentant une pluralite d'antennes | |
FR2973615A1 (fr) | Procede d'integration de signaux emis depuis un emetteur vers au moins un recepteur ultra large bande (ulb) ainsi que dispositif pour la mise en oeuvre du procede | |
EP0010481A1 (fr) | Procédé et dispositif de visualisation de cibles mobiles, notamment dans un équipement radar à période de répétition variable | |
EP1204884A1 (fr) | Systeme de reperage acoustique par bouees sous-marines | |
EP1477004B1 (fr) | Modulateur de frequence pour transmission numeriques | |
EP3400457A1 (fr) | Procédé d'élimination d'un signal en provenance d'un radar de bord | |
EP1949548A1 (fr) | Procede de detection de trajets en transmission impulsionnelle et dispositif correspondant | |
EP1221792A1 (fr) | Procédé et récepteur autodidacte pour déterminer les paramètres spatio-temporels d'un canal de propagation | |
WO2023208957A1 (fr) | Procédé de détection d'une pluralité de signaux de géolocalisation | |
CN112764064A (zh) | 一种实现预处理的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TP | Transmission of property | ||
ST | Notification of lapse |
Effective date: 20091030 |
|
RN | Application for restoration | ||
FC | Decision of inpi director general to approve request for restoration | ||
PLFP | Fee payment |
Year of fee payment: 16 |
|
PLFP | Fee payment |
Year of fee payment: 17 |
|
PLFP | Fee payment |
Year of fee payment: 18 |
|
PLFP | Fee payment |
Year of fee payment: 20 |