FR 3030964 A1 20160624 FR 1402936 A 20141219 INDENTIFICATION CONJOINTE DE SIGNAUX CONFONDUS EN TELECOMMUNICATIONS NUMERIQUES NON COOPERATIVES La présente invention concerne généralement les méthodes et systèmes de traitement du signal en télécommunications numériques, et plus particulièrement ceux de séparation aveugle de tels signaux ou de signaux de caractéristiques similaires. En télécommunications numériques, un émetteur cherche à transmettre une suite de bits d'information à un (ou plusieurs) récepteur(s). L'émetteur effectue une suite d'opérations de codage afin de rendre possible la transmission de l'information sur un support physique que l'on appelle aussi canal de propagation (tel qu'une fibre optique, ou tel qu'est plus particulièrement le cas ici, par propagation non contrainte dans l'espace). Ce canal de propagation n'est, généralement, pas parfait dans le sens où le signal reçu n'est pas une réplique exacte du signal émis (perturbations thermiques, réverbérations, dispersion). Au niveau de l'émetteur, on se donne une constellation (ASK, M-PSK, MQAM) qui est un ensemble fini de symboles dans le plan complexe et une fonction qui à nombre défini de bits consécutifs associe un symbole de la constellation. Chaque symbole de la suite est multiplié par une forme d'onde continue (ou filtre de mise en forme), décalée dans le temps pour former le signal en bande de base. Ce signal en bande de base est transposé à une certaine fréquence, dite fréquence porteuse, puis émis par une interface radio (une antenne dans le cas d'une transmission sans fil). A la réception, on effectue les opérations inverses pour reconstituer les bits d'information à partir du signal mesuré sur l'interface radio. Notamment, la démodulation est l'opération qui identifie à la réception les symboles effectivement émis par l'émetteur à partir du signal reçu. On parle de démodulation aveugle dès qu'un des paramètres de la transmission est inconnu de l'émetteur. En télécommunications coopératives, la plupart des paramètres sont définis par contrat entre l'émetteur et le récepteur (fréquence de porteuse, type de constellation, décalage temporel entre les formes d'onde, puissance d'émission). Les seules inconnues qui subsistent sont d'éventuelles erreurs de synchronisation et l'effet du canal de propagation. Ces paramètres sont, en général, estimés par émission régulière de séquences, dites pilotes connues des deux parties. L'opération compensant les effets du canal avant démodulation est appelée égalisation. En télécommunications non coopératives en revanche, l'intégralité des paramètres est inconnue - ou connue en ordre de grandeur seulement - et il n'y a pas de séquences pilotes pour faciliter l'égalisation. La séparation de sources est un problème générique où l'on cherche à extraire d'un signal composite les signaux source qui le constituent. La séparation aveugle de sources désigne la sous-classe des problèmes de séparation quand les caractéristiques des signaux source et de leur mélange sont inconnues ou très partiellement connues a priori. Des hypothèses statistiques sur les signaux source et/ou sur la nature de leur mélange s'avèrent nécessaires pour que ce problème soit soluble. Il est, par exemple, couramment supposé : que les sources sont mutuellement indépendantes ; que le nombre maximal de sources constituantes est connu ; que le mélange est linéaire ; qu'il existe des statistiques discriminantes entre les sources constituantes. Nous nous intéressons ici à la séparation et à la démodulation aveugle de signaux de télécommunication numérique dans un cadre non coopératif. Nous nous intéressons plus spécifiquement à un cas où il n'existe pas de statistiques discriminantes claires entre les sources. A titre d'exemple, on cite l'identification de signaux émis selon le protocole de transmission satellitaire « DoubleTalk Carrier-in-Carrier®». Cette technique de transmission satellitaire utilise une même bande de fréquence pour l'échange d'information en continu entre deux pairs. Un signal DoubleTalk intercepté consiste, donc, en la superposition en temps et en fréquence de deux signaux de mêmes caractéristiques. L'identification de ces signaux est un problème de séparation et de démodulation aveugle de sources. Parmi les méthodes de démodulation aveugle en télécommunications numériques, on s'intéresse à celles qui reposent : sur l'algorithme de classification EM (pour « Expectation- Maximization » en anglais, soit « Espérance-Maximisation ») pour l'estimation des paramètres du signal ; et sur l'algorithme de Viterbi pour la démodulation. Les documents (S. Barembruch, « Méthodes approchées de maximum de vraisemblances pour la classification et identification aveugles en communications numériques », rapport de thèse, Telecom ParisTech, 7 mars 2011) et (E. Punskaya, « Sequential Monte Carlo methods for digital communications », rapport de thèse, université de Cambridge, 2003) divulguent deux méthodes proches alliant le filtrage particulaire aux algorithmes EM et Viterbi. Dans ces deux documents, le modèle du signal à démoduler s'écrit comme suit: y(i) = + (i) j=0 ou encore : y (i) = +77(i) où S(i)est la version vectorisée des symboles s(i-J) à s(i) qui appartiennent à une constellation C, h est un filtre linéaire modélisant le canal de propagation, et n est un bruit blanc gaussien de variance 0-2. L'algorithme EM permet d'estimer le filtre h et la variance 0-2 selon un paradigme probabiliste et l'algorithme de Viterbi permet ensuite de reconstituer la séquence g de symboles vectorisés qui maximise la probabilité des observations du signal composite y. L'algorithme EM permet d'approcher itérativement la valeur du vecteur de paramètres 61= (h, o-2) qui maximise la vraisemblance du signal observé. A chaque itération, on applique successivement une étape E (pour espérance) et une étape M (pour maximisation). A l'itération k, on a un estimé 19(k) = (h(k), 0_2(k), ) des paramètres du système. - Etape E: cette étape cherche à évaluer l'espérance conditionnelle de la log-Vraisemblance des observations : Le(k)(0) = Etlog P0( (y(1: 0,5(1: I) = X(1: I)) I X(1: I) - 0(k)(. ly(1: /))) qui se réécrit Lek(0) =,o(k)(x) - L' e(y(01X(i) = X) (+cste(6)) i xeci Avec pour tout temps i, X(i) est la variable vectorielle dans CI qui désigne une valeur possible pour le vecteur de symboles S(i) effectivement émis, l'écriture y(1:1) désigne la suite d'observations (y(1), , y(I)) du signal composite y, f3i,e(k)(X) = Pek)(X(i) =X ly(1: /)) est la distribution de lissage et L' e(y(i)IX = X) = log 13 6,(y(i)IX = X) = IY - (111X)I2 /2o-2 . Le terme cste ne dépend pas de O. Il n'a donc pas besoin d'être calculé, et est omis dans la définition de l'espérance conditionnelle de la log-Vraisemblance Le(k). L'enjeu de l'étape E est de calculer les distributions de lissage /30(k)(X), ce qui se fait en calculant successivement les distributions dites de filtrage (Pek)(X(i) = X ly(1: i))) puis de lissage (igi,e(k)(X)), au moyen par exemple de l'algorithme Forward-Backard de Baum et Welch. - Etape M: une fois les quantités /30(k)(X) calculées, l'étape M met à jour les paramètres 61(k) en 19(k+1) selon 6,0c+1) = argmaxe Le(k)(0) En théorie, la suite de vecteurs 0(k) converge vers le vecteur 61 qui maximise la vraisemblance des observations. Cependant, appliqué à la démodulation d'un signal de communication numérique, l'algorithme EM est mis en défaut en termes de complexité calculatoire assez rapidement avec la taille du filtre h et/ou la taille de la constellation C. La contribution des deux documents cités ci-dessus est de dépasser cette limitation de complexité avec les notions de filtrage et de lissage particulaire qui permettent d'approximer le calcul de la distribution de lissage /3900(X)en limitant le nombre d'états X à évaluer. Ainsi, ils considèrent que igi,e(k)(X).----- Er [0,1] qui dépend éventuellement de l'échantillon y(i) ; le calcul de la distribution de filtrage Po(k)(X(i)= X ly(1),-,y(0) est réalisé récursivement, selon les étapes : - En entrée, soit t((zr,i_,,er,i_i))1, une approximation particulaire de P 9(k) - = X ly(1),---,y(i - 1)) ; - s rj1,r seuil o poser pour tout r, ii-1,r = Zi_ix et ei-1,r = - 51 Er > seuil o tirer une variable multinomiale: (Y1,-,YR) (Les Yr sont des variables aléatoires entières dont la somme vaut R et l'espérance de chaque variable vaut R - ) o soit E: -> e(t) = inf fr, '.21 Les documents cités proposent et comparent différentes méthodes pour déduire des coefficients de filtrage les coefficients de lissage. Nous retenons en particulier la méthode dite « Fixed Lag Smoothing ». Après l'estimation des paramètres du modèle du signal composite par cet algorithme EM-particulaire, les états de la technique cités ci-dessus utilisent une adaptation de l'algorithme de Viterbi pour réaliser la démodulation. Cette adaptation - dite Viterbi Stochastique - limite l'exploration du treillis du signal aux états des particules obtenues lors de la dernière itération de l'algorithme EM. Cette approche permet ainsi de résoudre le problème de la démodulation numérique avec une complexité applicable à des filtres et/ou des constellations de grande dimension. Ces états de la technique présentent plusieurs limitations pour être directement adaptés au problème de séparation et démodulation aveugle de signaux numériques confondus parce qu'entre autres, leur formalisation se limite à une seule source ; ils supposent qu'une synchronisation est acquise entre l'émetteur et le récepteur. Cependant, il est impossible pour un intercepteur de se synchroniser sur plusieurs signaux en même temps. Il s'ensuit que dans le cas d'un mélange de signaux, on ne peut supposer qu'une synchronisation est acquise, compromettant ainsi l'exploitation du modèle purement séquentiel des méthodes citées ; ils considèrent un canal de propagation causal. Cependant, le défaut de synchronisation induit sur chaque voie (c.à.d. sur chaque signal du mélange) un filtre non causal équivalent. En outre, un problème spécifique à l'exemple des communications numériques « DoubleTalk Carrier-in-Carriere» cité ci-dessus est que les signaux à séparer sont, théoriquement, émis à la même fréquence, aux mêmes instants et avec le même débit. Toutefois, dans la pratique, les caractéristiques temps-fréquence de ces signaux peuvent légèrement différer en raison d'une synchronisation non-parfaite des émetteurs de ces signaux. Le point de fonctionnement de ce système s'appuie sur des paramètres de synchronisation très proches (mais pas rigoureusement identiques) qui sont, en l'espèce, inconnus ou au mieux connus par leur ordre de grandeur. Il s'ensuit qu'il est impossible de prétendre séparer dans tous les cas de figure les signaux par les techniques classiques de filtrage temps-fréquence ou en négligeant un des deux signaux en faveur de l'autre. Un objet de la présente invention est de proposer des méthodes de séparation et de démodulation aveugles en télécommunications numériques de signaux compris dans un signal composite. Un autre objet de la présente invention est de proposer des méthodes permettant de séparer et de démoduler en aveugle des signaux numériques fonctionnant aussi bien avec des signaux de paramètres différents qu'avec des signaux de paramètres de modulation très proches, voire identiques. Un autre objet de la présente invention est de proposer des méthodes permettant de séparer et de démoduler en aveugle des signaux numériques en temps réel, c'est-à-dire avec des complexités en temps de traitement et en stockage de données linéaires selon le nombre d'échantillons du signal composite prélevés. Un autre objet de la présente invention est de proposer des méthodes permettant de séparer et de démoduler en aveugle deux signaux numériques stables numériquement. Un autre objet de la présente invention est de proposer des méthodes permettant de séparer et de démoduler en aveugle deux signaux numériques avec un canal et des paramètres évoluant dans le temps. Un autre objet de la présente invention est la démodulation aveugle et conjointe des différentes voies d'un signal « DoubleTalk 0 Carrier-in- Carriere». A ces fins, l'invention se rapporte, selon un premier aspect, à une méthode temps réel de séparation et de démodulation aveugle de signaux de télécommunication numérique, appelés voies, à partir de l'observation au moyen d'un unique capteur d'un signal composite comprenant ces signaux, les paramètres de ces voies incluant leur type de modulation, leur amplification, leur déphasage, leur temps de retard au niveau du capteur, leur fréquence et leur temps symbole, ces paramètres pour les différentes voies pouvant être différents, sensiblement ou parfaitement égaux, cette méthode comprenant les étapes suivantes : acquisition d'une première pluralité d'observations du signal composite réalisées au moyen de l'unique capteur ; estimation, à partir de la première pluralité d'observations acquise, des paramètres des voies au sens du Maximum de Vraisemblance par un algorithme Espérance-Maximisation, l'espérance conditionnelle de la log-vraisemblance étant calculée, dans cet algorithme, récursivement par une méthode de filtrage-lissage particulaire ; démodulation conjointe des voies selon un algorithme de Viterbi stochastique. La méthode de filtrage-lissage particulaire comprend une approximation particulaire de la distribution de filtrage, cette approximation particulaire assignant, en plus des poids, des multiplicités aux particules de cette approximation, une multiplicité étant le nombre de particules représentant un même état. La méthode de filtrage-lissage particulaire comprend une approximation particulaire de la distribution de lissage, cette approximation particulaire assignant, en plus des poids, des multiplicités aux particules de cette approximation, une multiplicité étant le nombre de particules représentant un même état. Cette méthode comprend, en outre, une étape d'acquisition d'une deuxième pluralité d'observations du signal composite réalisées au moyen de l'unique capteur ; une étape d'estimation, à partir de la deuxième pluralité d'observations acquise, des paramètres des voies au sens du Maximum de Vraisemblance par un algorithme Espérance-Maximisation, l'espérance conditionnelle de la log-vraisemblance étant calculée, dans cet algorithme, récursivement par une méthode de filtrage-lissage particulaire ; une étape d'estimation à long terme des paramètres, cette estimation à long terme des paramètres étant associée à la deuxième pluralité d'observations et étant une combinaison linaire des paramètres estimés à partir de la première pluralité d'observations et des paramètres estimés de la deuxième pluralité d'observations. L'invention se rapporte, selon un deuxième aspect, à un produit programme d'ordinateur implémenté sur un support mémoire, susceptible d'être mis en oeuvre au sein d'une unité de traitement informatique et comprenant des instructions pour la mise en oeuvre d'une méthode telle que présentée ci- dessus. D'autres objets et avantages de l'invention apparaîtront à la lumière de la description de modes de réalisation, faite ci-après en référence aux dessins annexés dans lesquels : les figures 1-6 illustrent schématiquement différentes étapes d'une méthode de séparation aveugle de sources selon divers modes de réalisation ; les figures 7-15 illustrent schématiquement les résultats obtenus à différentes étapes de la méthode de séparation appliquée à un signal composite « DoubleTalk Carrier-in-Carrier®». Dans un mode de réalisation, on enregistre au moyen d'un capteur unique un signal composite y(t) comprenant un mélange de N signaux source (sn(t))'=1..N. Les signaux (.97,(0).',..L.N, appelés dans la suite voies, contiennent les symboles d'information (sn(i))i=1..oe issus respectivement des constellations C n = (CO g=1..Gn et portés respectivement par des fonctions d'onde hn(.). Le signal composite y(t) s'écrit en théorie: i=+c° y(t) = An. el2' 5 flit hn(t - jT, - ± (t) n=1 où: An est le gain complexe du signal source sn (c.-à-d. An = lAnle'21r9n avec lAI l'amplification (ou gain) et cp' le déphasage (ou l'erreur de phase) sur la voie n) ; hn est la forme analytique du filtre de mise en forme de la voie n (par exemple, un filtre en Cosinus Surélevé) ; sn est la suite de symboles de la constellation Cn émis sur la voie n et s(i) est le i-ème symbole émis sur la voie n ; 8fn est le résidu de porteuse (ou plus généralement, la fréquence) sur la voie n. En particulier, ces résidus peuvent être sensiblement égaux ; Tn est le temps symbole sur la voie n. En particulier, ces temps symbole peuvent être sensiblement égaux (T1 = T2 = - - - TN) ; rn et le temps de retard sur la voie nau niveau de l'unique capteur. Ces temps sont choisis, sans perte de généralité, entre -Tn et 0; no est une perturbation aléatoire. Le signal composite y(t) est échantillonné à la période Te : j=+09 Y(i) Y(iTe) =1An. e1271-8fnTei hn(iTe iTn - Tn)sn(i) n(i) n=1 Te est choisie inférieure au double des temps symbole Tn, c'est-à-dire que pour tout 1 On note tn(i,j)(= tn,1i(Tn,T71)) l'écart temporel entre le &érine échantillon intercepté et la j-ème composante du vecteur Xn(i). Le signal y s'approxime en tronquant les fonctions d'onde selon l'équation y(i) rie,=1 An. ei27n5fnTei (Eii==21j+1 hn(tn(i, n, j (0) + n(i) On adjoint à chaque constellation Cn une relation d'ordre notée (Ci x x Cn)21+1, [Xn]' ((01,1([4]),...,OLN([4])),...,(02j+1,1([Xn1),...,021+1,n([4]))), avec pour tout n, Xn =[Xnili=1...2 j+i et OLni([Xn]n) = Xnrj ... On définit la relation d'ordre sur C1.21 x x c.N2J+1 par +1 pour U1 c Ci2j+1 r2J+1 " r2J+1 x x c.2J+1 . UN u2 N ' U2 (te, no), V./ Le module 2 implémente l'algorithme EM par filtrage et lissage particulaire. Ce module 2 fonctionne sur K itérations par trame ym et donne en sortie un vecteur paramètre OZ) = ([Anacm), Tn(mk), 6fn(m", Tn(km)in=1..1^1, 0-72//")) estimé sur la trame ym du signal composite y. Ce module 2 requière la définition de deux paramètres additionnels : Om pour initialiser l'estimation récursive du vecteur paramètre 0 ; et Pa(m,0) qui correspond à une « approximation particulaire » du dernier échantillon du signal précédant directement le premier échantillon de la trame yin (c'est-à-dire, l'échantillon ym_i(I), I étant le nombre d'échantillons par trame). Dans le cas où m= 1, cette approximation peut, par exemple, résulter d'un tirage uniforme des états possibles. Le module 3 réalise le suivi récursif à long terme des paramètres du signal. Ce module 3 prend en entrée le vecteur Om(K) délivré par le module 2 pour la trame m et un vecteur 0m-1, donnant une estimation à long terme des paramètres du signal jusqu'à la trame m-1. Ce paramètre Om_i provient en 30 général de ce même module 3 dans la trame m-1. Le module 3 produit l'estimation à long terme des paramètres Om en combinant linéairement Om(K) et Om_i. Un module 4 de démodulation met en oeuvre une variante de l'algorithme de Viterbi Stochastique à partir de la trame ym courante du signal composite et du vecteur paramètre Om estimé sur le long terme par le module 3. Ce module 4 de démodulation fournit en sortie une estimation .§1,m,-.-,'N,m des symboles d'information si,m,..-,sN,m émis sur la durée de la trame ym par chacune des voies .51,-,sN. Les symboles d'information estimés .§1,m,---,:ssN,m au fil des trames yi,..-,ym, respectivement, sur les voies, sont, respectivement, chargés dans des mémoires tampon (voir « buffer » sur la figure 2) pour former en sortie du système de séparation aveugle 10 des flux de symboles d'information parallèles correspondants aux signaux si,..-,sN. La figure 3 détaille le module 2 pour la mise en oeuvre de l'algorithme EM par filtrage et lissage particulaire et, en particulier, les interactions entre l'étape E mise en oeuvre par le module 2.1 et l'étape M mise en oeuvre par le module 2.2 au fil des itérations k=1,-,K sur une trame ym . Le module 2 encapsulant l'algorithme EM fournit en sortie une estimation du vecteur paramètre Or) sur la trame ym au bout de K itérations à partir de : la trame courante yn ; un jeu de particules noté Pa(m,O)approximant le signal à l'échantillon précédant le premier échantillon de la trame ym ; et un vecteur paramètre d'initialisation OZ initialisant l'algorithme EM. A l'itération k, le module 2.1 pour la mise en oeuvre de l'étape E calcule à partir de Pa(m,0) et du vecteur paramètre Om(k-1) (c.à.d. celui obtenu à l'itération précédente) un jeu de particules noté P-6,-(ma:/) dérivant de la distribution de lissage sur l'ensemble de la trame ym. Dans le cas où k= 1 et m> 1, 07(nic-1) = om-1 (donné par le module 3 sur la trame ym-i) ; si k= 1 et m = em(k-i.) =Oc, est un vecteur paramètre d'initialisation de l'algorithme EM. Selon les modes de réalisations, il peut être fourni par un opérateur, être choisi aléatoirement ou être un sous-produit du module de prétraitement 1. A l'itération k, le module 2.2 pour la mise en oeuvre de l'étape M raffine l'estimation de oTn(k) en Om(k+1) à partir de la trame ym, de am(k) et des particules P(m,1:/) obtenues par le module 2.1 pour la mise en oeuvre de l'étape E (voir figure 3). A cet égard, à l'itération k, l'étape E mise en oeuvre par le module 2.1 évalue l'espérance conditionnelle de la log-vraisemblance des observations : L9(-i) (6) = Esp (log P 0m( (ym (1: I), [Sn(1: I) = X na: pin) I Exn(1: Peoc-i)(. n)} qui se réécrit Le-i)(0) = > ,Gie-i)([4(i)].) - L'emcym(i) I [sn(i) = xn(oin) (+r(Om)) i n xn(DEei où : ,gie_i>axnumn) = Pee-,asn(i) = xnuAn iym(i),---,ymcm est la distribution de lissage et 2 x-11^1 j=21 L' 0m() I [X(i) = Mn) = 1Ym(i) Ane2n1 fnTel. T n)) /20-2 n=1 j=0 L" est un terme non-nécessaire à la variante de l'étape M mise en oeuvre par le module 2.2. Ce terme L" n'a donc pas besoin d'être calculé. Dans la suite, l'indice m de la trame ym est sous-entendu pour alléger les notations. Le(k)(0) se calcule au moyen d'une approximation particulaire. Une approximation particulaire est une écriture parcimonieuse de l'espace d'état adaptée à l'évaluation de Le(k)(0) dans le sens où le nombre de points d'évaluation est très petit par rapport aux dimensions de l'espace d'état. On considère deux types d'approximations particulaires. Selon une première approximation particulaire, dite simple, une approximation particulaire de la distribution de filtrage Pek)([4(0]' ly(1),-, y(0) s'écrit x-1R P e(k)([Xn(O]fl 13(1) ,y (j)) 4,i saxn(i)k-[zn,r(oin) où : R est un entier appelé nombre maximal de particules donné à l'initialisation de cette méthode particulaire; pour tout i et pour tout r De la même façon, une approximation particulaire simple de la distribution de lissage fli,o(k)([Xn(O]n) = PemaXn(OLIY(1),--,y(n) est une écriture de la forme: (k)([Xn(i)k)Wr,j (Vn(Din [Zi n,r(Oln) où : pour tout i et pour tout r Cette approximation peut, notamment, s'obtenir à partir de l'approximation particulaire de la distribution de filtrage P ek)([Xn(i)]n ly(1), , y(0) ci-dessus. Selon une seconde approximation particulaire, dite multiple, l'approximation particulaire de la distribution de filtrage PemaXn(Ok IY(1),-,Y(0) est une écriture de la forme Reff(i) P e(k)([Xn(O] IY (1) - - - y(0)-= - 11r,i 8aXn(i)in [Zn,r(017,) où : - R est un entier appelé nombre maximal de particules donné à l'initialisation de cette méthode particulaire ; - pour tout i,Reff(i) R est un entier appelé nombre effectif de particules à la génération i ; pour tout i et pour tout r Reff(i), er,i est un réel entre 0 et 1 appelé poids de filtrage ; pour tout i et pour tout r Reff(i), iirj est un entier appelé multiplicité de la particule ; - pour tout i, Eritr,i4,ti = 1 ; - pour tout n, tout i et tout r Reff(i),Vn,r(01n=Kzn,r,i(0) E Cn2J+1 est un n vecteur de symboles appelé support de l'approximation particulaire. Le filtrage particulaire selon cette première méthode particulaire est le calcul récursif de Po(k)([Xn(O] IY(1),-,Y(0) pour tout i selon l'approximation ci-20 dessus. De la même façon, une approximation particulaire multiple de la distribution de lissage 13i,e(k)([X,i(i)]) = P0ck»n(OlnlY(1),-,y(I)) est une écriture de la forme: Rieff(i) 130(k)([Xn(01n) .""' (13r,i - Vr,i 6([Xn(i)] [Zin,r(O] où : pour tout i et pour tout r R, cor,i est un réel entre 0 et 1 appelé poids de lissage ; pour tout i et pour tout r Reff(i), vri est un entier appelé multiplicité de la particule ; pour tout i, Er cor,ivr,i =1 ; pour tout n, tout i et tout r R, [Z' ,(0] (0) E Cn"" est un j n vecteur de symboles appelé support de l'approximation particulaire. Deux modes de réalisation peuvent être suivis: selon le paradigme de la première approximation particulaire ou selon le paradigme de la deuxième (dite approximation particulaire multiple) présentés ci-dessus. Avantageusement, l'approximation particulaire multiple ci-dessus de la distribution de filtrage, ainsi que celle de la distribution de lissage permettent d'éliminer la redondance des particules. Cette approximation multiple permet, en effet, de passer outre la redondance des particules en rassemblant les particules représentant un même état. Cette approximation particulaire multiple, qui permet de limiter le nombre de particules effectivement calculées et donc de diviser le temps de calcul et l'espace mémoire nécessaire à l'exécution de l'algorithme, est détaillée dans la suite. Les variables auxiliaires suivantes sont introduites: pour tout i, on note Vi:KReff(i)]] 41[1,Reff(i - 1)]1 l'application qui à une particule de la génération i désignée par son indice r associe une particule de la génération i - 1 appelée son ancêtre Vi(r); pour tout i et toute voie n, on se donne une variable Trn,i e {0,1} appelée indicatrice de saut ; pour tout i et tout r E j{1,Reff(011, on introduit un poids supplémentaire ani e [0,1] appelé poids auxiliaire, notamment tel que pour tout EReff(i) r ar,i - = 1 La définition précise de ces quantités auxiliaires se fera par récurrence dans la spécification complète du déroulement de la méthode pour la séparation aveugle des signaux. Pour assurer la stabilité des calculs, les poids auxiliaires et les poids de filtrage sont remplacés, respectivement, par leurs logarithmes notés Lani et Les écarts de temps tn(i,j) entre l'échantillon i du signal composite y et le j - ème échantillon du vecteur de symboles 4(0 contribuant à ym(i) sont calculés en même temps que les particules. Les indicatrices de saut,n-i e {0,1} vaut 1 s'il est nécessaire de considérer l'arrivée d'un nouveau symbole entre les vecteurs Sn(i - 1) et S(i) et 0 sinon. On note les ensembles : Pc(,k) 0n, = (Reff(0, [Un [(t( i, , (Larj,Lri 4tr,i,[Zn,r(01., Vi(r))r) (appelé approximation auxiliaire de l'échantillon i) ; (m, i) = (12 ef f(i) , [Tr n,dn, Rt , (gr,i , (01 n, Vi(r)),..) (appelé approximation particulaire de la distribution de filtrage en 0; Pco(k)(111, = (11 eff(0, RtnO,D) , (cor"! ,vr [Z n r(01n (appelé approximation n r particulaire de la distribution de lissage en 0. On définit sous le même format des ensembles de particules multiples supplémentaires : Pé(ck) 011, = (1keff0), [Ttn,d. Rtn (Lécr,i, [2n,r (r)) produit par le module 2.1.1 de la figure 4 et utilisé par le module 2.1.2 (appelé approximation auxiliaire rééchantillonnée) ; 1),-,k)(rn, = (rZeff(i), Rtn 0, D)i n ((T)r,i [2n,r W]jr) produit par le module 2.1.4 de la figure 4 et utilisé par le module 2.2 (appelé approximation mélangée de la distribution de lissage en j). Ces ensembles sont les unités d'échange de données entre les sous modules constituant l'algorithme EM particulaire mis en oeuvre par les modules 2.1.1- 2.1.4 et 2.2 de la figure 4. Notons que dans le détail de ces ensembles, les indices de trame yni et d'itération k ont été sous-entendus pour alléger les notations. On utilise la notation Pe(m,1:0 (respectivement Pe(m,1:I) ) pour désigner l'ensemble constitué de Pe(m,1),---,PC)(m,l) (respectivement p (k) (m, p(k) (m, /)) ). Les ensembles Pa(m,i) et P(m, i) ont de plus la propriété d'être classés selon l'ordre lexicographique sur Ci2J+1 x ---xCn2J+1 --= (C1 x x Cri)21+1 : pour tout r, i E [[1, Reff0A, r 5_ É [Z,,,(01n [Z',t(O]n. Dans la méthode présentée, un calcul de la forme logEiciebi,bi Les méthodes d'estimation particulaire mettent en jeu un noyau de transition (fonction de l'espace d'état au carré dans l'intervalle [0,1]). Comme les fonctions d'onde sont susceptibles de ne pas être causaux, les noyaux de transitions étudiés et utilisés dans les documents cités dans l'état de la technique ne sont donc pas appropriés au problème qui nous intéresse ici. On se donne un entier 8de l'ordre de J- max Tn/Te. On définit pour tout i w1([4(0]7,[4(i- 1)]n) = [x,i(i+6)]nPe(k)(Yi,-* ,Yi+8 ,[Xn(i ,[Xn(i + (Mn), wi([4(i = E[xnco]nwi([4()1.,[4(i-1)1n) et wi([X.,(0],.,, [X, - = WiaXn(Oin, [Xn(i - 1)])/ WiaXn(i - ) qui est un noyau de transition de C121+1 x ---x Cn2J+1dans lui-même. Des approximations de type Monte-Carlo sont utilisées en pratique pour évaluer ce noyau. Le module 2, et en particulier le module 2.1 réalisant l'étape E, se décompose en sous modules comme détaillé dans la figure 4 - pour chaque itération k o pour chaque échantillon i - le module 2.1.1 correspond à une étape de rééchantillonnage des particules qui prend en entrée l'ensemble de particules Pe)(m,i -1) et donne en sortie l'ensemble de particules rééchantillonnées Pe((k)(m, i - 1) ; - le module 2.1.2 correspond à une étape de filtrage particulaire qui prend en entrée Pe(ck)(m,i) , la trame yin , l'estimation du vecteur paramètre 07(0 obtenue à l'itération précédente et qui donne en sortie la distribution de filtrage Pe(k)(m,i) pour l'échantillon i ainsi que l'approximation auxiliaire Pc(ck(m,i) pour l'échantillons i; - le module 2.1.3 correspond à une étape de lissage particulaire. Dans un mode de réalisation, cette étape de lissage particulaire met en oeuvre un algorithme dit « fixed-lagg smoothing » qui prend en entrée la distribution de filtragePe(k)(m,i) de l'échantillon i et donne en sortie la distribution de lissage Pe)(m,i--.d) de l'échantillon i , A étant un paramètre entier défini à l'instanciation du système ; o les distributions de lissage PC) (m, i - A) pour i = 1,-,/ sont agrégés dans un buffer (voir figure 4) afin de produire en sortie la distribution de lissage de la trame entière P 0) définir: eo: fteff( - 1)11 -> [[1, Reff(1 - 1)1 par eo(i-) = inffr,Er, 0) définir la bijection E1,r: [ri Reff (r, -> K = Ki E K, Edoctixdoe)Pd(cri) 1-1 - noter pour tout n,tout 1[1 , Reff(r, 011 (0 = = LWi ([Zn,r(i - 1)]n) 0 = O 1.7. = + log Pe(() (y(i)1 [Z',(rj.),i1n) - o V; 0-, = Poser Reff(i) = Eràee1-1)Reff(r,i) - définir la bijection (01,02): 111 ,Reff()1 -> Urr_ftieff(i-1){r} x D,Reff(r,)1 une permutation telle que r1 Pour l'étape de lissage mise en oeuvre par le module 2.1.3 (figure 4), la méthode connue sous le nom de « fixed-lag smoothing » est par exemple utilisée. On se donne un entier A, - en entrée, la distribution de filtrage de l'échantillon i Pt(m,i) = (Reff(i),[un,dri,(gr,i 41-r,i,[Zn,r(01n,Vi(r)) poser V1,0(r) -= V1(r) pour tout 1 Le module 2.1.4 est configuré pour mettre en oeuvre une étape de mélange stochastique des particules obtenues par le module 2.1.3. Il prend, donc, en entrée l'ensemble des particules lissées P,(m,1:/) en renvoie un nouvel ensemble de particules noté Pia(m,1:/). Ce module 2.1.4 permet de s'immuniser contre les risques de convergences vers un optimum local de l'algorithme EM. Il consiste, pour tout i, à ne retenir de P'(m, 0 qu'une seule particule tirée au sort. Il est à noter que, si cette étape est appliquée à chaque itération k de chaque trame yin, l'algorithme devient instable. On ne l'applique, donc, qu'à certains couples (m, k), ceux qui appartiennent à un ensemble It, cet ensemble étant pour être de moins en moins dense avec l'évolution de (m,k).Par exemple, on peut considérer It tel que, (m, k) E It ssi -1)- K + k est le carré d'un entier. La procédure mise en oeuvre dans ce module est : en entrée: P(k)(m,1: 1) = ((R ef f ,{(tri(i, l)) .1 , (cor i ,vr [Z )) est n r l'approximation de la distribution de lissage ; si (m, k) It o pour tout i e o poser peff(i) =1 o tirer un indice ro dans [[1,Re'ff(i)]] avec une probabilité proportionnelle à coi,r, o poser tà,,i = 1, = 1 et [2a(i)]n = [Znixo(O] - si (m, k) It o pour tout i E ft1,I, o poser r?eff(i)----Re'ff(i) o pour tout r E fieff (0]] o poser ior,i = 1, -er,i = 1 et [2n,r(0] = [4,r(0]7.' en sortie : P (m, 1:1) = ((P'ef f ,l(tn(i, fin, fr)r,i [2n,r (Oljr)) est l'approximation mélangée de la distribution de lissage. Le module 2.2, mis en évidence sur les figures 3 et 4, correspond à l'étape M de l'algorithme EM. A l'itération k, il prend en entrée la trame de signal ym, le vecteur paramètre em (k-1) estimé à l'itération précédente et les particules mélangées P(m,1:0 obtenues en fin de l'étape E (c.à.d. par le module 2.1 et en particulier le sous module 2.1.4). Il donne en sortie une estimation raffinée (k) em du vecteur paramètre. Cette estimation est choisie pour faire décroitre l'espérance conditionnelle de la log-vraisemblance des observations selon le vecteur paramètre 61: (k) n(k-1) (k-i) (uni )e(k-1) (Uni ). m De nombreuses méthodes peuvent être envisagées à cet effet. Dans un mode de réalisation, la procédure suivante est mise en place: On note le vecteur paramètre Om(k) estimé à la k-ème itération sur la m-ème oc) = (mn(krn),Tn(rnk), fn(ink) Tem)]n=1..N, 0_711(k)) trame em sous la forme 61m(k) = (em(9) où L ' est le nombre total de paramètres à estimer. La permutation entre les paramètres et les indices I peut être quelconque, voire aléatoire et différente à chaque itération. on note e® = pour tout 1,1 1 ou éventuellement par le module d'amorçage 1.3 si m = 1. A la convergence de l'algorithme EM sur la trame yin, tous les estimateurs de paramètres ont une variance proportionnelle à 0-2m(K). On pose donc : Pm = 2 (K) -2 Y - a m -2 1 - 2 = cro = +00 Y1 ' -2 , 2 (K) am-1 0" m Puis on actualise les paramètres selon : em = + (1 - fm) - (Or) - cm-1) Le paramètre de suivi y fixé par l'utilisateur entre 0 et 1 permet de donner de l'élasticité pour prendre en compte l'évolution du canal. Les quantités fm et eirl, sont respectivement un paramètre de suivi et une variance modifiés afin de combiner linéairement de façon optimale l'estimation sur la trame m de l'algorithme EM Om(k) et l'estimation à long terme jusqu'à la trame m-1,8m_i. Le module 4 de démodulation, détaillé dans la figure 5, démodule les voies. Ce module 4 implémente un algorithme de Viterbi stochastique. Ce module 4 de démodulation produit essentiellement les trames de symboles estimés .§1,m,§2,m,-..,N,in pour chacune des voies. 2 (K) Y - a m Le module 4 de démodulation comprend deux blocs fonctionnels. Le sous-module 4.1 est une instanciation de l'étape E de l'algorithme EM, limitée à la reproduction des modules 2.1.1 et 2.1.2 correspondant à l'étape de filtrage particulaire. Ce sous-module 4.1 II prend en entrée : - la trameym ; - le vecteur paramètre Om estimé à long terme sur la trame ym par le module 3 pour le suivi à long terme du vecteur paramètre O. Le sous-module 4.1 réalise avec ce vecteur paramètre Om une nouvelle étape de filtrage et produit en sortie : un ensemble de particules fi(m,1:/) à destination du sous-module 4.2. Cet ensemble peut, selon l'approximation particulaire utilisée, être un ensemble de particules ou un ensemble de particules multiples. Il est à noter que dans cette sortie, les poids de filtrages obtenus ne sont pas utilisés dans la suite ; une approximation particulaire auxiliaire Pa(m,I) du dernier échantillon de la trame ym. C'est cette approximation qui sert à initialiser le module 2 et en particulier le module 2.1.1 de la trame ym+i via leur paramètre d'entrée Pa(m +1,0) Le sous-module 4.2 réalise la démodulation des voies par un algorithme de Viterbi stochastique. Il prend en entrée : le vecteur paramètre Om estimés à long terme jusqu'à la trame ym par le module 3 pour le suivi à long terme du vecteur paramètre e; le support de l'approximation particulaire /5(m,1:0 de la trame donné par le sous-module 4.1. On obtient par l'algorithme de Viterbi stochastique la séquence vectorisée de maximum a posteriori (P,(0]). Les symboles émis sur chaque voie sont estimés par la séquence ((in)) selon la procédure: 15in5In - en entrée : les séquences vectorielles de maximum a posteriori : ([]n) et les indicatrices de saut [n] pour tout 1 i I et tout 1 Comme montré sur la figure 6, le module 1 comprend trois sous-modules. - le sous-module 1.1 assure la « bufferisation » du signal y en observation à son entrée pour former les trames ym ; le sous-module 1.2 procède conjointement au sous-échantillonnage et au filtrage de la trame ym pour donner en sortie la trame ym qui est le signal directement manipulé par les autres modules du système de séparation aveugle 10. Ce sous-module 1.2 utilise la connaissance plus ou moins précise des temps symboles Tn disponible dans le système à travers le vecteur 0,-1 délivré par le module 3 pour le suivi à long terme du vecteur paramètre O. Pour la première trame, le vecteur paramètre Bo est estimé par exemple en aveugle par le sous-module 1.3; le sous-module 1.3 permet de réaliser une première estimation grossière du vecteur paramètre O . Cette étape est utilisée uniquement sur la première trame du signal et est appelé étape d'amorçage. Par exemple les fréquences peuvent être amorcées à partir de la propriété de circularité des constellations considérées, les temps symbole par la méthode de la raie spectrale, les amplifications par un partage équitable de l'énergie du signal intercepté, les temps de retard et les erreurs de phase qui sont a priori bornés peuvent être pris égaux à 0 ou à des valeurs tirées aléatoirement. La procédure de filtrage et décimation conjointes du signal mise en oeuvre dans le sous-module d'adaptation 1.2 est la suivante : en entrée: la trame ym(i)=Einv=iAn. ei2n8 f N (Eiii+cecoitn(iTe. - j Tn - inbn(n) (i) échantillonnée à la cadence Teo ; -lin est un filtre en racine de cosinus surélevé de temps symbole Tn ; les Tn sont sensiblement égaux ; une estimation des temps d'échantillonnage de chaque voie Ti,-,TN via le paramètre Om_1 ; on pose T, min Tn et fi le filtre en racine de cosinus surélevé de temps symbole 2T, ; on filtre et sous-échantillonne conjointement ym en ym à la cadence T, selon l'équation Ym(i) = Ijirn(j) - ri UT e0 iTe) e0 - en sortie, la trame ym est donnée en entrée de la chaine de traitement aval. L'algorithme de séparation proposé est une généralisation de l'algorithme EM calculé sur une approximation particulaire. Cet algorithme est appliqué à un modèle analytique du signal composite pour en extraire les signaux sources. Le décodage est assuré par un algorithme de Viterbi Stochastique. Le terme temps réel est justifié, ici, dans la mesure où le signal est « bufferisé » en trames de I (100 à quelques milliers) symboles et que la complexité temporelle de l'algorithme de séparation et ses besoins en ressources mémoire sont linéaires selon la taille de la trame. Ainsi il est possible de l'implémenter de telle façon que le temps de traitement d'une trame soit égale au temps d'émission de ses I échantillons quelle que soit la valeur de la taille. Ainsi la taille des trames peut être adaptée à des signaux dont les paramètres sont constants dans le temps (trames longues) ou à des signaux dont les paramètres varient rapidement dans le temps (trames courtes), le module de suivi à long terme se charge dans ce dernier cas de composer une estimation des paramètres de variance moindre. Avantageusement, les systèmes et méthodes décrits ci-dessus introduisent la notion de particules avec multiplicité qui permet de limiter grandement le nombre de particules à traiter et donc d'accélérer largement l'exécution de l'algorithme. En effet, avec ce paradigme, le nombre de particules effectif utilisé pour une trame lors d'une itération EM s'adapte automatiquement à la difficulté de ce traitement. Ainsi, sur les premières itérations des premières trames, beaucoup de particules sont effectivement utilisées ; quand la convergence des paramètres est amorcée, ce nombre diminue progressivement avant de se stabiliser à un nombre de particules largement plus petit que celui qui serait nécessaire dans le paradigme classique de filtrage particulaire quand la convergence est acquise. Avantageusement, les méthodes et systèmes décrits ci-dessus permettent le traitement d'un large éventail de configuration : de différents types de modulations (par exemples M-ASK, M-PSK, QPSK, M-QAM), le type de modulation peut être différent sur chacune des voies ; des bandes de fréquences variées pour les différentes voies, qu'elles soient différentes ou égales (ou sensiblement égales à cause d'erreurs de synchronisation entre les voies à envisager par exemple pour la séparation de signaux Double Talk) ; des débits symboles variés pour les différentes voies, qu'ils soient différents ou égaux (ou sensiblement égaux à cause d'erreurs de synchronisation entre les voies, à envisager par exemple pour la séparation de signaux Double Talk) ; des gains de transmission au niveau des différentes sources qu'ils soient différents ou sensiblement égaux ; de temps de retard et les erreurs de phase sur chaque voie peuvent être quelconques. Les systèmes décrits sont capables de détecter automatiquement, de combiner de façon optimale et d'exploiter à son avantage toute asymétrie entre les signaux source. Dans une mise en oeuvre illustrative de la méthode de séparation aveugle 20 décrite ci-dessus, un signal composite comprenant deux signaux « DoubleTalk Carrier-in-CarrierO» est simulé selon les paramètres suivants : le signal y est composé de deux voies : N = 2 ; la période d'échantillonnage à l'interception est prise comme unité 25 temporelle de référence : Te = 1 ; les filtres de mise en forme i et 1-1.2 sont des filtres en racine de cosinus surélevé de roll-off 0.35 ; les temps de symboles sont T1 = 2 et T2 = 2 en unités réduites ; les résidus de porteuse sont de 8f1 = i.i0, 8f2 = 3.10-6 en unités 30 réduites ; les décalages de phase sont de (pi =1T/10 et de (p2 = Tr/4; les retards sont -c1 = -0.3 - Ti et T2 = -0.8 - T2; - les deux modulations sont des QPSK; les amplifications sont telles que 1A21= 0.95 -1A11 le bruit est blanc, gaussien, le rapport entre la puissance moyenne des signaux et celle du bruit est de 20dB. Le signal composite simulé est analysé sur 20 trames de taille équivalente à 500 symboles par voie. Après une étape d'adaptation du signal, la méthode de séparation aveugle est instanciée avec les réglages suivants : le nombre maximal de particules est fixé à 200 ; le paramètre de lissage est A = 40 ; le canal est supposé quasi-invariant : y = 0.9 ; il y a K = 10 itérations EM par trame. Le déroulement de l'algorithme de séparation sur ce signal composite est montré sur les figures 7 à 14 et sur la figure 15. Les figures 7 à 14 montrent visuellement les résultats de quelques itérations typiques de l'algorithme, les sous figures représentent respectivement de gauche à droite et de haut en bas : une reconstruction de la voie 1 dans le plan complexe (cette reconstruction est obtenue en soustrayant au signal composite l'estimation de la voie 2 et les interférences inter-symboles sur la voie 1. Les niveaux de gris représentent la densité locale des mesures sur la voie 1 reconstruite). La valeur de SINR est une estimation en décibels du rapport entre la puissance du signal sur la voie 1 et toutes les perturbations qui s'y superposent (interférences entre symboles, voie 2 et bruit). une reconstruction de la voie 2 (cette reconstruction est obtenue en soustrayant au signal composite l'estimation de la voie 1 et les interférences inter-symboles sur la voie 2. Les niveaux de gris représentent la densité locale des mesures sur la voie 2 reconstruite). La valeur de SINR est une estimation en décibels du rapport entre la puissance du signal sur la voie 2 et toutes les perturbations qui s'y superposent (interférences entre symboles, voie 1 et bruit). une reconstruction du signal composite (cette reconstruction est obtenue en soustrayant au signal composite les interférences inter- symboles sur la voie 1 et sur la voie 2. Les niveaux de gris représentent une densité locale de mesures sur le signal composite reconstruit. Les points noirs correspondent à ce que serait le signal composite sans bruit, les couples de chiffres au-dessus de chaque point désignent les couples d'états des constellations des voies 1 et 2 qui correspondent à cet état de signal composite). les données quantitatives dans cette sous figure sont :Tr : le numéro de la trame traitée, It : l'itération courante, MP : le nombre effectif moyen de particules utilisées pour traiter l'itération courante et SI : l'écart type du bruit estimé à cette itération. une reconstruction dans le domaine fréquentiel de la voie 1, de la voie 2 et du bruit. une reconstruction dans le domaine temporel des filtres de mise en forme sur chacune des voies Sur la figure 7 on observe les signaux après la phase d'adaptation (0 itérations). Ils sont indiscernables. Les SINR pour les voies 1 et 2 ne peuvent pas encore être estimés ; leur valeur inconnue est notée « ? ». Sur la figure 8 on observe les signaux après reconstruction avec les paramètres d'initialisation de l'algorithme (c'est-à-dire, après une seule itération) ; le niveau de bruit est élevé et on ne relève pas à l'ceil de structures pertinentes sur les différentes reconstructions. Sur la figure 9, après 4 itérations, on remarque que le niveau de bruit a diminué, on commence à prévoir une QPSK sur la reconstruction de la voie 2. La figure 10 montre la dernière itération de l'algorithme sur la première trame ; sur chacune des voies, on voit une QPSK se former. Le niveau de bruit est encore cependant élevé, l'algorithme n'a pas fini de converger à la fin des itérations allouées à la première trame. Cependant on voit sur la figure 11 que la convergence se poursuit sur la 30 deuxième trame : le changement de trame n'a pas créé de bouleversement dans l'algorithme, le mode de convergence est continu. C'est à la quatrième itération de la deuxième trame (figure 12) que la convergence de l'algorithme est acquise. Les QPSK sont bien dessinées sur chacune des voies reconstruites et le niveau de bruit estimé est inférieur au niveau de bruit du signal intercepté (l'étape d'adaptation a aussi fait gagner quelques dB de SNR). Les figure 13 et 14 sont typiques de la suite des résultats de l'algorithme : la convergence des paramètres reste stable avec le temps, l'algorithme suit l'évolution du canal et corrige l'accumulation des erreurs d'estimation des résidus de porteuses en faisant légèrement évoluer les erreurs de phases au fil des trames. Il est à noter que le paradigme des particules avec multiplicité se montre, dans cet exemple, efficace pour adapter le nombre de particules effectives à la difficulté du cas à traiter. Si dans les premières itérations 88 à 14 particules effectives sont utilisées (sur les 200 disponibles), ce nombre se stabilise à 4-5 dans les itérations suivantes, ce qui est un gain en besoin de ressources temps et matérielles important par rapport à l'état de l'art en démodulation aveugle. La figure 15 présente l'évolution de l'estimation des paramètres au fil des trames (estimation « globale ») et des itérations (estimation « locale »). On y observe notamment que l'estimation des erreurs de phase évolue légèrement pour compenser la variance sur l'estimation des résidus de porteuses. REVENDICATIONS1. Méthode temps réel de séparation et de démodulation aveugle de signaux de télécommunication numérique, appelés voies, à partir de l'observation au moyen d'un unique capteur d'un signal composite comprenant ces signaux, les paramètres de ces voies incluant leur type de modulation, leur amplification, leur déphasage, leur temps de retard au niveau du capteur, leur fréquence et leur temps symbole, ces paramètres pour les différentes voies pouvant être différents, sensiblement ou parfaitement égaux, cette méthode comprenant les étapes suivantes : acquisition d'une première pluralité d'observations du signal composite réalisées au moyen de l'unique capteur ; estimation, à partir de la première pluralité d'observations acquise, des paramètres des voies au sens du Maximum de Vraisemblance par un algorithme Espérance-Maximisation, l'espérance conditionnelle de la log-vraisemblance étant calculée, dans cet algorithme, récursivement par une méthode de filtrage-lissage particulaire ; démodulation conjointe des voies selon un algorithme de Viterbi stochastique. 2. La méthode de la revendication précédente, dans laquelle la méthode de filtrage-lissage particulaire comprend une approximation particulaire de la distribution de filtrage, cette approximation particulaire assignant, en plus des poids, des multiplicités aux particules de cette approximation, une multiplicité étant le nombre de particules représentant un même état. 3. La méthode de la revendication 1 ou 2, dans laquelle la méthode de filtrage-lissage particulaire comprend une approximation particulaire de la distribution de lissage, cette approximation particulaire assignant, en plus des poids, des multiplicités aux particules de cette approximation, une multiplicité étant le nombre de particules représentant un même état. 4. La méthode de l'une quelconque des revendications précédentes, comprenant, en outre, une étape d'acquisition d'une deuxième pluralité d'observations du signal composite réalisées au moyen de l'unique capteur ;une étape d'estimation, à partir de la deuxième pluralité d'observations acquise, des paramètres des voies au sens du Maximum de Vraisemblance par un algorithme Espérance-Maximisation, l'espérance conditionnelle de la log-vraisemblance étant calculée, dans cet algorithme, récursivement par une méthode de filtrage-lissage particulaire ; une étape d'estimation à long terme des paramètres, cette estimation à long terme des paramètres étant associée à la deuxième pluralité d'observations et étant une combinaison linaire des paramètres estimés à partir de la première pluralité d'observations et des paramètres estimés de la deuxième pluralité d'observations. 5. Produit programme d'ordinateur implémenté sur un support mémoire, susceptible d'être mis en oeuvre au sein d'une unité de traitement informatique et comprenant des instructions pour la mise en oeuvre d'une méthode selon l'une quelconque des revendications précédentes.