FR 2730881 A1 19960823 FR 9502051 A 19950222 Description L'invention concerne un système de séparation de sources pour traiter des signaux d'entrée formés par des mélanges linéaires instantanés de signaux primaires issus des sources et pour délivrer au moins un signal primaire estimé, les mélanges faisant intervenir des coefficients de mélange, le système comprenant des moyens de séparation de sources ayant une pluralité d'entrées reliées aux signaux d'entrée et au moins une sortie pour délivrer le signal primaire estimé, les moyens de séparation de sources déterminant adaptativement des coefficients de séparation servant à extraire les signaux primaires estimés. Elle concerne également l'application d'un tel système à la réception de signaux électriques, sonores, électromagnétiques. Il peut s'agir par exemple d'une antenne, d'un autoradio, d'un téléphone mains libres. On connaît la technique de séparation de sources de signaux primaires qui consiste à traiter des mélanges de signaux primaires pour fournir une estimation de chaque signal primaire. Ladite technique s'applique à des signaux primaires issus de sources indépendantes qui ne sont disponibles que sous forme desdits mélanges. Il peut s'agir de mélanges linéaires convolutifs ou de mélanges linéaires instantanés. Ils peuvent être engendrés par des mécanismes de propagation de signaux primaires et/ou par des mécanismes de superposition de signaux de plusieurs sources ou d'autres causes. Dans le cas général, la technique de séparation de sources opère en aveugle, c'est-à-dire que les sources sont supposées inconnues, indépendantes, avec des mélanges inconnus. Pour cela on détecte plusieurs échantillons desdits mélanges à partir desquels la mise en oeuvre d'algorithmes de séparation permet de restituer une ou plusieurs estimations des signaux primaires d'origine. Une telle technique, appliquée à la séparation de mélanges linéaires instantanés de signaux, est par exemple divulguée dans le document intitulé "Separation of Independent Sources from Correlated Inputs" de J.L. Lacoume et P.Ruiz, IEEE Transactions on Signal Processing, vol. 40, n 12, décembre 1992, pages 3074 à 3078. Pour effectuer la séparation de sources, c'est-à-dire obtenir en sortie une estimation de chaque source formant le mélange, ce document révèle une procédure de calcul de cumulants. Pour cela il enseigne d'adapter des paramètres du système de séparation de sources de façon à annuler des cumulants de signaux de sortie qui s'expriment en fonction de cumulants déterminés sur les signaux d'entrée issus des mélanges, les cumulants étant d'ordre supérieur à deux. En annulant ces cumulants, on déduit indirectement des coefficients de mélange inverse pour effectuer une transformée inverse de la transformée effectuée par l'application des opérations de mélange aux signaux primaires. L'enseignement de ce document conduit à donner au système de séparation de sources une structure directe. De plus, comme les signaux d'entrée sont statistiquement dépendants car issus de mélanges, les équations qui relient les cumulants des signaux de sortie aux cumulants des signaux d'entrée sont très complexes. Une telle technique s'avère donc être complexe à mettre en oeuvre et ne fournit pas de solution simple dans le cas de signaux réels issus généralement de mélanges linéaires instantanés. Un but de l'invention est de simplifier le traitement conduisant à distinguer les sources en délivrant des signaux primaires estimés. Ce but est atteint avec un système de séparation de sources qui comprend des moyens de caractérisation ayant une pluralité d'entrées reliée aux signaux d'entrée pour estimer, à partir de ceux-ci, les coefficients de mélange, et ayant au moins une sortie pour délivrer des coefficients estimés de mélange, les moyens de séparation ayant au moins une autre entrée pour recevoir les coefficients estimés de mélange, les moyens de séparation comprenant des moyens pour transformer les coefficients estimés de mélange en coefficients de séparation. Avantageusement, on obtient une estimation des coefficients de mélange eux-mêmes et non pas celle des coefficients de mélange inverses. Cette propriété assure une séparation de sources beaucoup plus précise, les signaux primaires estimés étant plus proches des signaux primaires. De plus, les moyens de séparation de sources peuvent alors avoir une structure récursive qui est plus compacte. Préférentiellement, les moyens de caractérisation comprennent des moyens d'estimation pour calculer, sur les signaux d'entrée, des estimations de cumulants d'ordre supérieur à deux et des moyens de calcul de coefficients estimés pour transformer les estimations desdits cumulants en coefficients estimés de mélange. Avantageusement, on simplifie le traitement en établissant une relation entre les cumulants des signaux d'entrée et les cumulants des signaux de sources. Comme les sources sont statistiquement indépendantes, il s'ensuit que les équations sont bien plus simples, donc plus faciles à résoudre que dans l'art antérieur, ce qui nécessite des moyens de traitement moins importants. Plus particulièrement, les moyens de calcul de coefficients estimés comprennent un réseau de neurones pour transformer, selon une fonction de transformation prédéterminée, les estimations desdits cumulants en coefficients estimés de mélange. Un domaine intéressant d'utilisation de l'invention est celui dans lequel les sources de signaux primaires et/ou les récepteurs sont mobiles entraînant des conditions variables de transmission. Ces différents aspects de l'invention et d'autres encore seront apparents et élucidés à partir des modes de réalisation décrits ci-après. L'invention sera mieux comprise à l'aide des figures suivantes données à titre d'exemples non limitatifs qui représentent: Figure 1: un schéma d'un système de séparation de sources selon l'invention. Figure 2: un schéma d'une unité de caractérisation selon l'invention. Figure 3: un schéma d'une unité de séparation de sources selon l'invention comprenant une unité de transformation des coefficients estimés de mélange en coefficients de séparation. Figure 4: un schéma partiel d'un premier exemple de réalisation de l'unité de séparation dans le cas d'une structure directe. Figure 5: un schéma partiel d'un second exemple de réalisation de l'unité de séparation dans le cas d'une structure récursive. Sur la figure 1 sont représentées les sources primaires , Sl à Sn, constituées par exemple par différents émetteurs radio- fréquences. Pour identifier et séparer ces émetteurs, on place des capteurs Cl - Cn en divers points de réception. Chacun de ces capteurs est par exemple une antenne suivie d'un amplificateur et d'un démodulateur. Les capteurs fournissent des signaux El(t) à En(t) issus de mélanges de signaux primaires Xl(t) à Xn(t) qui, soit sont délivrés directement par les sources Sl à Sn, soit servent à moduler des signaux diffusés par les sources S1 à Sn. Les mélanges qui se produisent entre les signaux primaires sont directement liés à la propagation de ces signaux dans l'espace. Ces mélanges peuvent être rendus complètement ou approximativement linéaires et instantanés en utilisant des émetteurs et des capteurs adéquats notamment en choisissant les types de modulateurs et démodulateurs utilisés, ainsi que les fréquences de fonctionnement du système, par exemple en choisissant des émetteurs utilisant une modulation d'amplitude. Sur la figure 1, les signaux détectés El(t) à En(t) entrent dans une unité 10 de séparation de sources qui délivre des signaux primaires estimés Xl(t) Xn(t), la variable t étant par exemple le temps. Les signaux mélangés d'entrée El(t) à En(t) sont des mélanges linéaires instantanés issus, dans le cas général, de p signaux primaires inconnus Xl(t) à Xp(t) tels que: (1) Ei(t)= a.j Xj(t) avec lsi n et l1jsp, i et j étant des variables courantes. Sur la figure 1 on a n=p. Ces signaux primaires sont supposés statistiquement indépendants préférentiellement jusqu'à au moins l'ordre 4. Cette préférence peut être déplacée à au moins l'ordre 3 dans certains cas. Les coefficients aij sont des coefficients de mélange qui définissent les contributions du signal Xj(t) vers le signal mélangé Ei(t). Les coefficient de mélanges aij, qui sont inconnus, sont constants ou lentement variables. Pour la suite de la description on définit des signaux normalisés primaires tels que: (2) Yj(t) = aj.Xj(t) pour l1j p et des coefficients de mélange normalisés tels que: (3) aij=aij/ajj pour iej. Pour i=j, on a aii= ajj=1. Après cette normalisation, l'équation (1) devient: (4) Ei(t) = Yi(t) + Ea,,.Yj(t) pour l1i n, l1j p et ji. -7 Les coefficients normalisés aij de mélange seront par la suite appelés simplement coefficients de mélange. Pour fournir des signaux primaires estimés Xk(t) à partir des signaux d'entrée Ei(t), l'unité de séparation 10, dans le cas o elle possède une structure directe, met en oeuvre des coefficients de séparation Cki. Elle délivre les signaux primaires estimés en effectuant la transformation inverse (k indice courant de sortie): (5) Xk(t) =. Cki. E(t) avec l i n et l1k n. Les coefficients de séparation Cki sont déterminés dans l'unité de séparation 10 à partir de coefficients estimés &ij de mélange qui constituent une estimation des coefficients aij de mélange qui sont inconnus. Cette estimation est effectuée dans une unité de caractérisation 15 qui, à partir des signaux d'entrée Ei(t), délivre les coefficients estimés de mélange (figure 1). D'une manière analogue, l'unité de séparation (10) peut avoir une structure récursive mettant en oeuvre des coefficients de séparation dki. La figure 2 représente un schéma de l'unité de caractérisation 15. Elle comprend une unité d'estimation 16 suivie d'une unité 17 de calcul de coefficients estimés. L'unité d'estimation 16 reçoit l'ensemble des signaux d'entrée Ei(t) et délivre une estimation d'un ensemble A de grandeurs caractéristiques des signaux d'entrée Ei(t). Ces grandeurs caractéristiques de l'ensemble A sont choisies pour qu'il soit possible premièrement de les calculer à partir des signaux d'entrée Ei(t), et deuxièmement de les exprimer en faisant intervenir les coefficients de mélange ai, et un ensemble B de grandeurs caractéristiques se rapportant aux signaux primaires normalisés Yj(t) bien que ceux-ci soient inconnus, les grandeurs caractéristiques de l'ensemble B étant de même nature que les grandeurs caractéristiques de l'ensemble A. Dans la suite de la description les grandeurs caractéristiques choisies vont être des cumulants calculés à partir des signaux mélangés Ei(t). Les relations qui relient l'ensemble A et l'ensemble B par l'intermédiaire des coefficients de mélange constituent une propriété intrinsèque du système décrit. Selon l'invention, l'unité d'estimation 16 détermine une estimation de l'ensemble A, estimation suffisante pour qu'il soit possible d'en déduire des coefficients estimés &ij de mélange dans l'unité 17 de calcul de coefficients estimés. Selon un mode préférentiel de réalisation, les grandeurs caractéristiques constituant l'ensemble A sont des auto-cumulants associés à chaque signal mélangé Ei(t) préalablement centré, ainsi que des cumulants croisés associés à chaque couple de signaux mélangés Ef(t) et E g(t) préalablement centrés (f et g étant des indices entiers courants, l1fsn, l1g n). De manière analogue, les grandeurs caractéristiques constituant l'ensemble B sont des auto- cumulants associés à chaque signal primaire normalisé Yi(t) préalablement centré. La définition de ces grandeurs est en soi bien connue de l'homme du métier. Selon l'invention, on considère les cumulants (auto et croisés) d'ordre supérieur à 2. Les cumulants les plus faciles à calculer sont ceux qui ont l'ordre le plus faible. Néanmoins les cumulants d'ordre trois doivent être suffisamment grands pour que le calcul se limite à ces seuls cumulants. Cette condition ne se rencontre que pour un nombre limité de types de sources. Préférentiellement, selon l'invention, on calcule les cumulants d'ordre 4 ce qui permet de couvrir une plus large gamme de signaux rencontrés expérimentalement. L'homme du métier peut choisir des cumulants ayant un ordre différent sans sortir du cadre de l'invention. Dans un souci de clarté, un rappel des définitions des cumulants est fourni ci- après en se limitant aux seuls cumulants nécessaires à l'exemple traité, c'est-à-dire aux auto-cumulants et aux cumulants croisés d'ordre 4, en se limitant à deux signaux, et sans faire intervenir de décalage temporel entre les deux signaux considérés. Dans le cas général, considérons deux signaux U(t) et V(t) qui varient en fonction d'une variable t, par exemple le temps. On définit des signaux centrés par: (6) u(t) = U(t) - E{ U(t)} (7) v(t) = V(t) - E{ V(t)} Dans les équations ci-dessus la notation E{.} représente l'espérance mathématique de l'expression entre parenthèses. Les cinq cumulants d'ordre 4 associés à ces signaux centrés sont (en abrégeant par u et v respectivement les signaux u(t) et v(t)): (8) cum4(u) = cum(u,u,u,u) = E{u4} - 3[E{u2}]2 (9) cum31(u,v) = cum(u,u,u, v)= E{u3v} - 3E{u2} E{uv} (10) cum22(u,v) = cum(u,u,v,v)= E{u2v2} - E{u2}E{v2} -2[E{uv}]2 (11) cum31(v,u) = cum(u,v,v,v)= E{v3u} - 3E{v2}E{uv} (12) cum4(v) = cum(v,v,v,v)= E{v4} -3[E{v2}]2 Dans la pratique on ne calcule pas les cumulants exacts, mais seulement des cumulants estimés. Différentes techniques sont connues de l'homme de métier pour calculer des cumulants. A titre d'exemple, deux méthodes seulement sont fournies ci-après. Une première méthode, simple pour le calcul de cumulants, consiste à effectuer les opérations suivantes: 1. Dans une fenêtre temporelle, prélever M échantillons pour chaque signal U(t) et V(t). 2. A partir de ces échantillons, effectuer des estimations des espérances mathématiques des signaux U(t) et V(t) en calculant deux estimateurs et . Chaque estimateur est par exemple la moyenne arithmétique du signal considéré calculée sur les M échantillons prélevés. Ces estimateurs et sont des nombres réels. 3. A partir des équations 6 et 7 et des estimateurs et , on déduit des estimateurs et se rapportant aux signaux centrés u(t) et v(t). Ces seconds estimateurs contiennent M échantillons qui sont, par exemple, calculés de la manière suivante. On calcule l'estimateur respectif ou des signaux centrés u(t) ou v(t) en soustrayant de chaque échantillon du signal initial U(t), respectivement V(t), l'estimateur de l'espérance mathématique E{U(t)}, respectivement l'estimateur de l'espérance mathématique E{V(t)}. [équations 6 et 7 appliquées aux estimateurs]. 4. On applique la même méthode de calcul pour les estimateurs de toutes les espérances mathématiques qui figurent dans les équations 8 à 12. Par exemple, l'estimateur de E{u4} sera obtenu en calculant la moyenne arithmétique de [ ]4 sur M échantillons. De même pour les autres espérances mathématiques à calculer. 5. On en déduit une estimation pour chaque cumulant défini par les équations 8 à 12 en remplaçant les espérances mathématiques par leurs valeurs estimées respectives. On peut obtenir une plus grande précision dans l'estimation des cumulants en appliquant une seconde méthode qui consiste à subdiviser la fenêtre contenant les M échantillons en plusieurs sous-fenêtres. La première méthode est appliquée pour chaque sous-fenêtre pour délivrer des cumulants partiels. On calcule ensuite les cumulants pour l'ensemble des M échantillons (équations 8 à 12) en calculant les moyennes arithmétiques des cumulants partiels. D'autres méthodes de calcul des cumulants peuvent être mises en oeuvre. Afin de bien faire comprendre l'invention, considérons, à titre d'exemple, le cas simple de la séparation de 2 sources en utilisant 2 capteurs. Dans cet exemple, seuls seront considérés les auto-cumulants et les cumulants croisés d'ordre 4. On détecte deux signaux El(t) et E2(t) qui, après avoir été centrés, deviennent les signaux el(t) et e2(t). De même, yl(t) et y2(t) sont les signaux centrés correspondant aux signaux normalisés Yl(t) et Y2(t) supposés indépendants. En appliquant ce qui a été exposé précédemment on obtient les 5 équations suivantes: (13) cum4(el) = cum4(yl) + a412 cum4(y2) (14) cum31(el,e2) = a21 cum4(yl) + a312 cum4(y2) (15) cum22(el,e2) = a221 cum4(yl) + a212 cum4(y2) (16) cum31(e2,el) = a321 cum4(yl) + a12 cum4(y2) (17) cum4(e2) = a421 cum4(yl) + cum4(y2) Ces équations renferment: i) les éléments de l'ensemble A déjà défini, c'est-à- dire les éléments caractérisant les signaux fournis par les capteurs, à savoir: cum4(el), cum31(el,e2), cum22(el,e2), cum31(e2,el), cum4(e2). ii) les éléments de l'ensemble B déjà défini, c'est-à- dire les éléments caractérisant les signaux primaires normalisés, à savoir: cum4(yl), cum4(y2) qui sont inconnus. iii) deux coefficients de mélange inconnus: a12, a21. Disposant de cinq équations, il est donc théoriquement possible de déterminer les éléments inconnus. Dans le cas de plus de deux sources, avec un nombre de capteurs égal au nombre de sources, les cumulants cum4 (équations 13 et 17) et les cumulants cum31 (équations 14 et 16) par exemple se généralisent respectivement tels que: n (18) cum4(er) = cum4(Yr) + a 4j cum4(yj) j-l, jr (19) cum31(e", e,) = a,zcum4(y,) + a35cum4(y,) + a3. ajaLcum4 (yj) avec l1jÉn, jer et jes, et l1r6n, l1s n, et r*s,r,s étant des variables entières courantes. Le but de ce calcul (cas à deux sources et deux capteurs) est de déterminer les coefficients de mélanges a12, a21. En réalité, le calcul est conduit pour déterminer uniquement des estimations des coefficients de mélange a12, a21 sous forme de coefficients estimés 12l, a21' Le calcul des coefficients estimés a12, &21 est effectué dans l'unité de caractérisation (15) représentée sur la figure 2. Premièrement, l'unité d'estimation (16) détermine les estimations des cumulants des signaux mélangés, puis l'unité (17) de calcul de coefficients estimés délivre les coefficients estimés &12, &21- L'unité d'estimation (16) calcule les cumulants cum4(el) à cum4(e2) selon les équations 8 à 12. L'unité d'estimation DSP peut être constituée d'un calculateur, d'un microprocesseur ou d'une unité de traitement numérique du signal. Lorsque tous les cumulants précédents ont été calculés, ils sont traités par l'unité de calcul de coefficients 17 qui calcule les coefficients de mélange a12 et a21 à partir des équations 13 à 17. En réalité, la résolution de ces équations fournit non pas les coefficients vrais a12 et a21 mais des valeurs proches constituant les coefficients estimés a12 et a21 du fait que ce sont des cumulants estimés et non pas des cumulants exacts qui sont utilisés pour résoudre les équations 13 à 17. L'unité de calcul de coefficients 17 peut elle aussi être constituée d'un calculateur, d'un microprocesseur ou d'un dispositif de traitement numérique du signal appliquant une méthode de résolution d'équations. Pour résoudre ces équations, il est possible d'utiliser des méthodes de résolution d'équations non-linéaires qui sont connues de l'homme de métier. Il peut s'agir des méthodes de Gauss- Newton, de Levenberg-Marquardt, de Powell-Fletcher ou des méthodes de descente de gradient classiques. Une autre approche consiste à réduire le nombre d'équations en combinant entre elles les équations 13 à 17. On obtient ainsi les équations suivantes: (20) "a1cum4 (el) +a12cum4 (e2) -a"l2cum31 (el, e) -cum31(e2, e1) =0 (21) î21cum4(el) +oe2cunm4(e2) -cum31 (e, e2) -al2a2lcum31 (e2, el) =0 L'unité de calcul des coefficients 17 est alors programmée pour résoudre les équations 20 et 21 par exemple à l'aide d'une des méthodes connues indiquées préalablement. Il est possible de combiner ces deux équations pour ne conserver que des équations renfermant un seul coefficient de mélange a12 ou a21 plus faciles à résoudre. Une autre approche peut consister non pas à résoudre exactement le système d'équations 13 à 17 comme cela vient d'être dit, mais à se limiter à réaliser une approximation de la fonction qui relie les cumulants estimés calculés aux coefficients estimés dij de mélange qui sont à déterminer. En effet, d'une manière globale, les équations 13 à 17 peuvent s'écrire sous la forme: (22) aij = F(A). L'équation 22 n'a pour objet que de relier les coefficients inconnus aij avec l'ensemble A sans prendre en considération l'ensemble B qui, bien qu'inconnu également, n'est pas à déterminer. Si G(.) est une fonction qui constitue une approximation de F(.) et si G(.) est plus facile à mettre en oeuvre expérimentalement, on peut alors se limiter à calculer: (23) aJi = G(A). Cette fonction F(.) est caractéristique du nombre de signaux primaires et de capteurs nécessaires à la mise en oeuvre de l'application considérée. Ceci constitue un avantage car le système est par lui-même indépendant des autres paramètres de l'application. On conçoit aisément que l'unité de calcul de coefficients 17 peut subir préalablement dans une première phase un apprentissage pour résoudre la fonction F(.). Dans une seconde phase, les données mesurées expérimentalement peuvent être traitées selon la méthode apprise. Or, il est connu qu'il est possible d'obtenir une approximation d'une fonction F(.) par une fonction G(.) en utilisant un réseau de neurones multicouches ou un perceptron multicouches ou un réseau de neurones en arbre. L'utilisation d'un réseau de neurones est particulièrement intéressante car elle ne nécessite pas de résoudre, à chaque fois, pour chaque lot de cumulants, les équations 13 à 17 dans le système de séparation lui-même. Il suffit que le réseau de neurones apprenne préalablement à résoudre la fonction G(.) pour ensuite l'utiliser pour chaque lot de cumulants, le système de séparation de sources devenant ainsi dédié au type de fonction G(.) appris. Pour adapter le système de séparation de sources à un nombre de capteurs ou à un nombre de signaux primaires à séparer différent, il suffit de faire réapprendre au réseau de neurones la nouvelle fonction G(.) relative à ce nombre nouveau de capteurs et/ou à ce nombre nouveau de signaux primaires. Pour générer la base d'exemples nécessaires à l'apprentissage, on choisit des valeurs arbitraires d'autocumulants de signaux primaires et des valeurs arbitraires de coefficients de mélanges normalisés aij (membres de droite des équations 13 à 17). A l'aide de ces équations, on calcule les autocumulants et les cumulants croisés des signaux mélangés (membre de gauche desdites équations). Puis on introduit ces cumulants des signaux mélangés sur les entrées du réseau de neurones et on impose que ses sorties délivrent les coefficients de mélanges normalisés aij préalablement choisis. Ces coefficients délivrés a.ij constituent les résultats à obtenir en sortie du réseau de neurones. L'apprentissage du réseau de neurones va donc consister à déterminer les coefficients synaptiques du réseau de neurones pour approximer la correspondance signaux d'entrée/signaux de sortie du réseau de neurones conformément aux résultats des équations 13 à 17 précalculés. Ceci présente l'avantage de ne pas avoir à résoudre les équations complexes 13 à 17. Les coefficients estimés iij de mélange ainsi calculés sont entrés dans l'unité 10 de séparation de sources pour démêler les signaux d'entrée Ei(t) et fournir des signaux primaires estimés Xk(t). La séparation de sources effectuée par l'unité 10 de séparation de sources peut s'écrire sous la forme: n ú ( t) =n (24) Xkt) = 3 Cki. E(t) i-l pour une structure directe, et n (25) Xkt = Ek(t) d. i ( t) i-1 isk pour une structure récursive. Les coefficients Cki' dki sont des coefficients de séparation qui sont déduits des coefficients estimés &ij de mélange qui entrent dans l'unité 10 de séparation. La transformation des coefficients estimés de mélanges aij en coefficients de séparation Cki, dki dépend de la structure directe ou récursive de l'unité 10 de séparation. Lorsque l'unité 10 de séparation possède une structure récursive, on peut être conduit à effectuer l'une ou l'autre des transformations suivantes: (26) dki = &ki ou dki = ak, (i)/aio,(i) La première de ces transformations correspond au cas o le signal primaire estimé correspondant au signal primaire de kième rang se retrouve sur la kième sortie de l'unité de séparation. Mais il peut apparaître des cas o des problèmes de stabilité ne permettent pas d'assurer cette correspondance. Il peut alors être nécessaire de faire que le signal primaire estimé correspondant au signal primaire de kième rang soit délivré par une autre sortie de l'unité de séparation. Il faut alors déterminer des coefficients de séparation dki appropriés permettant d'effectuer cette permutation. C'est l'objet de la seconde transformation qui fait intervenir sur les indices une loi de permutation a(i). Dans le cas de deux sources, avec deux coefficients estimés &12 et "21, les coefficients de séparation dki sont choisis tels que: (27) - si &12.&21 S1, on a d12 = &12 et d21 = a2. pour la première transformation, et (28) - si a&.22l >1, on a d,2 = 1/a21 et d21 = 1/&2 pour la seconde transformation. Lorsque l'unité 10 de séparation de sources possède une structure directe, il faut faire une inversion de matrice telle que: 2 + C 1a -lksn (29) C i ki.a ij = 0 ljn, jk, i',ioj avec Ckk = 1. Les transformations des équations 26 à 28 dans le cas de la structure récursive ou l'inversion de matrice de l'équation 29 dans le cas de la structure directe sont réalisées dans l'unité DSP de transformation 11, par exemple un calculateur, un microprocesseur ou une unité de traitement numérique du signal. Dans ce cas (figure 3), l'unité 10 de séparation de sources comprend l'unité de transformation 11 qui transforme les coefficients estimés &ij de mélanges en coefficients de séparation Cki, dki et une sous-unité 12 de séparation proprement dite qui reçoit, sur une pluralité d'entrées, les signaux mélangés Ei(t) et qui délivre au moins un signal primaire estimé Xk(t). Les coefficients de séparation Cki, dki arrivant sur une autre pluralité d'entrées. La figure 4 représente une partie d'une sous-unité de séparation 12k à structure directe. Elle comprend une pluralité d'entrées I1 à In chacune recevant un signal mélangé El(t) à En(t). Chacune de ces entrées est reliée à un moyen de multiplication 131 à 13n dudit signal par le coefficient de séparation Ckl à Ckn affecté à l'entrée. Les sorties de tous les moyens de multiplication 131 à 13n sont reliés à un sommateur 125 pour additionner tous les signaux et délivrer le signal primaire estimé Xk(t). La sous-unité de séparation 12 (figure 3) contient autant de sous-unités partielles 12k que de signaux primaires estimés Xk(t). La figure 5 représente une sous-unité de séparation 12 ayant une structurerécursive pour un exemple destiné à délivrer deux signaux primaires estimés X1(t) et X2(t) à partir de deux signaux mélangés E1(t) et E2(t). Cette structure comprend un premier sommateur 112 ayant une entrée reliée au signal E1(t) et une sortie qui délivre le signal primaire estimé X1(t). Un second sommateur 212 a une entrée reliée au signal E2(t) et une sortie qui délivre le signal primaire estimé X2(t). Une autre entrée du premier sommateur 112 est reliée à la sortie du second sommateur 212 à travers un moyen de multiplication 111 qui pondère le signal de sortie du second sommateur par un coefficient -d12. De la même manière, une autre entrée du second sommateur 212 est reliée à la sortie du premier sommateur 112 à travers un autre moyen de multiplication 211 qui pondère le signal de sortie du premier sommateur par un coefficient -d21. Les sommateurs 112, 212 et les moyens de multiplication 111, 211 peuvent faire partie d'un calculateur, d'un microprocesseur ou d'une unité de traitement numérique du signal correctement programmé pour réaliser les fonctions décrites. REVENDICATIONS 1. Système de séparation de sources pour traiter des signaux d'entrée formés par des mélanges linéaires instantanés de signaux primaires issus des sources et pour délivrer au moins un signal primaire estimé, les mélanges faisant intervenir des coefficients de mélange, le système comprenant des moyens de séparation de sources ayant une pluralité d'entrées reliées aux signaux d'entrée et au moins une sortie pour délivrer le signal primaire estimé, les moyens de séparation de sources déterminant adaptativement des coefficients de séparation servant à extraire les signaux primaires estimés caractérisé en ce que le système comprend des moyens de caractérisation ayant une pluralité d'entrées reliée aux signaux d'entrée pour estimer, à partir de ceux-ci, les coefficients de mélange, et ayant au moins une sortie pour délivrer des coefficients estimés de mélange, les moyens de séparation ayant au moins une autre entrée pour recevoir les coefficients estimés de mélange, les moyens de séparation comprenant des moyens pour transformer les coefficients estimés de mélange en coefficients de séparation. 2. Système selon la revendication 1 caractérisé en ce que les moyens de caractérisation comprennent des moyens d'estimation pour calculer, sur les signaux d'entrée, des estimations de cumulants d'ordre supérieur à deux et des moyens de calcul de coefficients estimés pour transformer les estimations desdits cumulants en coefficients estimés de mélange. 3. Système selon la revendication 2 caractérisé en ce que les moyens de calcul de coefficients estimés comprennent un réseau de neurones pour transformer, selon une fonction de transformation prédéterminée, les estimations desdits cumulants en coefficients estimés de mélange. 4. Système selon les revendications 2 ou 3 caractérisé en ce que les moyens de séparation sont des moyens à structure directe, les moyens de séparation comprenant une unité de transformation des coefficients estimés de mélange en coefficients de séparation adaptés à la structure directe. 5. Système selon les revendications 2 ou 3 caractérisé en ce que les moyens de séparation sont des moyens à structure récursive, les moyens de séparation comprenant une unité de transformation des coefficients estimés de mélange en coefficients de séparation adaptés à la structure récursive.