FR 2838582 A1 20031017 FR 0204616 A 20020412 mémoires pour alimenter autant de codeurs. d PROCEDE ET DISPOSITIF DE DETECTION DE DONNEES TRANSMISES PAR ETALEMENT DE SPECTRE Domaine technique et art antérieur La présente invention concerne un prc édé de détection de donnces transmises par étalement de spectre par séquence directe ainsi qu'un dispositif pour mettre en _uvre le procédé. Le domaine d'application de l' invention est le domaine des transmissions numériques à étalement de spectre par séquence directe plus communément appelées transmissions numériques SSD (ESS3 pour " Etalemet de Spectre par Séquence Directe "j. Les transmissions numériques ESSD sont utilisces dans les systèmes à accès muliples ar répartition de coces cor,.muném_rt appelés systèmes 7'IRC (MRC pou " Accès Multiple à Répartition par les Codes ") cu encore systèmes CDMA (CDMA pour " Code Division Multiple Access). Dans un système CDMA, les utilisateurs partagent la même bande de fréquences. La discrimination entre utilisateurs est rendue possible par l' assignation de codes d'étalement différents à chaque utilisateur. Le récepteur d'un utilisateur particulier reçoit non seulement le signal qui lui est destiné mais également des signaux destinés à dautres utilisateurs. I1 est alors nécessaire de minimiser les interférences qui résultent des signaux requs destinés aux autres utilisateurs, interférences communément appelées bruit d'accès multiples ou bruit MAI (MAI pour " Multiple Access InterferenCe "). Par ailleurs, il est également B 13857.3/PR nécessaire de diminuer le bruit inhérent au canal de transmission. Le signal rec,u en entrée d'un récepteur numérique est représenté par le vecteur complexe e tel que: e=A.d+n (1) o - A est une matrice complexe qui dépend de la réponse impulsionnelle du canal et des codes d'étalement relatifs au sysième CDMA, - d est un vecteur complexe constitué des donnéss émises par les utilisateurs et que l' on cherche à c;c_, - - n es' un vecteur de bruit aléatoire. 1, Le ré-epteur calaule un vecteur d, qui es le vecteur estimé de d, à partir de la connaissance cu vecteur e et de la matrice A préalablement évaluce ou calaulée. Le vecteur d doit être aussi proche que possible de d. Pour cela, on cherche à annuler les distorsions et les perturLations introduites par la chaîne de transmission entre l'émetteur et le récepteur. Cette opération d'annulation des distorsions et des perturbations s'appelle égalisation. Une égalisation parfaite conduit à d = d. Une difficulté essentielle réside dans le fait que le récepteur ne dispose pas de suffisamment d'informations pour résoudre l'équation (1). La résolution de l'équation (1) est alors basée sur l'utilisation de critères d'optimisation dans le but d'obtenir une égalisation parfaite. Diverses techniques sont connues pour résoudre B 13857.3/PR 3 2838582 l'équation (1). Parmi ces techniques, la technique de détection conjointe consiste à estimer la globalité du vecteur d, c'est-à-dire la totalité des donnces émises par l' ensemble des utilisateurs. Dans ce cas, le vecteur estimé d est calculé par une équation du type: d L-l AH e (2) dans laquelle:. - Ll est la matrice inverse d'une matrice L, .. communément appelée matrice d'égalisation, qui est lG une matrice carrée de taille NxN et dont l'expression dépend de critères d'optimisation choisis pour résoudre l'équation (1), et - A représente la matrice symétrique herr.itierne de la matrice A (c'est-à-dire la matrice transposée conjuguce de A). D frérents exemp'es ce détection conjointe peuvert être ctés. ur premier ex-mple concerne une détection selon laquelle la matrice L stécrit AH.A et qui vise à annuler le bruit MAI (égalisation ZF, ZF pour " Zero Forcing "). Un deuxTème exemple concerne une détection selon laquelle la matrice L s'écrit AH.A+62I, o 62 est la variance du bruit n et I la matrice identité (égalisation MMSE, MMSE pour " Minimum Mean Square Error "). Dans ce dernier cas, la distance entre le vecteur estimé d et le vecteur d est évalué sans distinguer l'effet du bruit n de l'effet du bruit MAI. De façon générale, pour résoudre l'équation (2), le récepteur commence par calculer le vecteur z=AH.e à partir de la matrice A qu'il connaît et du vecteur des données reques e. L'équation (2) peut alors B 13857.3/PR se mettre sous la forme: z = L. d (3) Une technique classique de résolution de l'équation (3) consiste à inverser la matrice L et à multiplier à droite les 2 termes de l'équation (2) par la matrice L-l, inverse de la matrice L. Cette technique doit cependant être écartée dans les cas o la matrice L est de grande taille et o le temps disponible pour résoudre l'équation (3) est très limité. C'est la cas, par exemple, du mode TDD (TDD pour " Time Division Duplex ") du système UMTS (UMTS pour " Universal Mobile Telecommunication System "). La résolution de l'équation (3) coit alors en e==et être ef_ectuée en moins de 0,666 ms pour une matrice L dont la taille p-ut ê_r_ "'evirn OOx_CC c_- -c erts complexes. a technqe classique par -nversion de la matrice L mentionnée ci-dessus ne permet pas la résolution de léquation (3). Devant la complexité des problèmes posés, des architectures vectorielles ont été propososs (cf. " A vector multiprocessor for real-time multi-user detection in spread spectrum communication ", IEEE Internationa1 Conference on Application Specific Systems, ArchitectureS and Processors, Boston, Juillet 2000). Ces architectures n'ont cependant été conçues que pour résoudre des systèmes linéaires de petite taille (matrice 32 x 32). La transpos.tion diecte de telles architectures à des systèmes de détection du domaine de l' invention dont la taille peut atteindre, par exemple, 500 x 500 (cf. ci-dessus) n'est pas envisageable car cela conduirait également à mettre en B 13857.3/PR _uvre des systèmes matériels, par exemple des processeurs, d'une extrême complexité. Une autre technique, qui nécessite moins de calculs, consiste à tirer partie des propriétés de la matrice L pour en faire une décompos it ion de Cholesky donnée par la formule L=U.U o U est une matrice triangulaire supérieure de taille NxN dont les éléments notés uij sont déterminés à partir des éléments de la matrice L notés li,. La résolution de l'équation (3) se déroule alors en 2 étapes. Une première étape est une décomposition de Cholesky de la matrice L pour obtenir la matrice U. Il vient alors: A,. Z = L.d = U.U.d (4) Une deuxième étape consiste résoudre successivement ltéauation: Z = UH.y () qui permet d'obtenir la matrice y à partir des matrices UH et z connues, et l'équation: y=U.d (6) qui permet d'obtenir d à partir des matrices U et y caleuléss précédemment. Ces opérations sont programmées sur un ordinateur classique, voire un processeur de traitement de signal (DSP). Cependant, les performances de ceux-ci ne permettent également pas d'atteindre les performances temps éel pour les matr ces de grande taille mentionnses ci-dessus. L' invention ne présente pas les inconvénients mentionnés ci-dessus. B 13857.3/PR Exposé de l' invention En effet, l' invention concerne un dispositif de détection coujointe de données transmises par étalement de spectre par séquence directe, le dispositif S comprenant des moyens de mise en _uvre d'une décompos it ion de Choleeky pour cal cul er une matrice U telle que: L = UH x U. o L est une matrice d'égalisation et U une matrice lO triangulaire supérieure. Les moyens de mise en _uvre d'une décomposition de Choleeky comprennent: - des moyens pour calculer les coefficients diagonaux ui,i de la matrice U selon la formule: i-1 Ui i = li, i - luk/ l ( 1) OU k=1 1 lii est le coefficient diagoral de la lgne de rar.g i de la matrice L et uki est le coefficient de la ligne de rang k et de la colonne de rang i de la matrice U. et - des moyens pour calculer les coefficients non diagonaux de la matrice sous la forme: 1-1 * 1 i j - uk 1.uk j ui, j = ui'i les moyens pour calauler les coefficients non diagonaux comprenant: un réseau de N cellules l4C montées en série, chaque cellule MAC comprenant une entrée d'accumulation, une sortie d' accumulation, une entrée de retour, une sortie de retour et une entrce de coefficient, la B 13857.3/PR sortie d'accumulation de la cellule de rang p (p=1 à N) étant reliée à l'entrce d' accumulation de la cellule de rang p+1, - une cellule diagonale comprenant une entrée daccumulation, une entrce de coefficient de soustraction, une entrée de coefficient de division, une sortie de résultat et une sortie de retour, l' entrée d' accumulation de la cellule diagonale étant reliée à la sortie d'accumulation de la cellule MAC de rang N. - des moyens de chargement des coefficients u k,i sur les entrées de retour des cellules MAC, - des moyens de chargement des coefficients ux,, sur les entrces de coefficient des cellules MAC, - des moyens de chargement des coe-ic-ents ui,i sur l'entr'e de coefic ent "e Qivision ce la cellule diagona'e, et - des moyens de chargement des coefficients li,j sur l'entrée de coefficient de soustaction de la cellule diagonale. Selon un premier mode de éalisation de l' invention, les moyens de chargement des coefficients U k,i sur les entrées de retour des cellules MAC sont constitués d'un ensemble de N multiplexeurs (Ml, M2, MN) commandés par un premier état de commande, chaque multiplexeur ayant une première entrce, une deuxième entrce et une sortie, 'a sortie d'un multiplexcur de rang p (p=1, 2,..., N) étant reliée à l'entrée de retour de la cellule MC de même rang, la première entrce du multiplexeur de rang p étant relice à un coefficient U k,i différent d'un multiplexeur à l'autre, la deuxième B 13857.3/PR 8 2838582 entrce du multiplexcur de rang p étant reliée à la sortie de retour de la cellule MAC de rang p+1, la deuxième entrée du multiplexeur de rang N étant relice à la sortie de retour de la cellule diagonale, le premier état de commande établissant une connexion électrique directe entre la première entrée et la sortie de chaque multiplexeur. Selon un deuxième mode de réalisation de l' invention, la sortie de retour de la cellule MAC de rang p (p=1, 2,..., N) est relise à l'éntrée de retour de la cellule MAC de rang p-l, les moyens de chargement des coefficients u k,i sur les entrées de retour des cellules MAC sont constitués d'un multiplexcur commardé par un premier état de commande, le multiplexeur ayant -m--^s _n:^-e, u-e deuxime entr_e et une sc--ie, la sortie du mu' lexeur étant reliée G l'entrce de retour de la cellule C de rang N. la première entrée du multiplexeur étant reliée, successivement, à un parmi les coefficients u k,i et la deuxième entrce du multiplexeur étant relice à la sortie de retour de la cellule diagonale, le premier état de commande établissant une connexion électrique directe entre la première entrée du multiplexeur et la sortie du multiplexeur. Selon encore une caractéristique supplémentaire du dispositif de l' invention, les cellules MAC, la cellule diagonale et l' ensemble des multiplexeurs comman"és par un deuxième état de commarde constituent un réseau permettant de calculer le vecteur estimé d qui correspond au vecteur d des donnces émises par les utilisateurs. B 13857.3/PR L' invention concerne également un procédé de détection conjointe de données transmises par étalement de spectre par séquence directe, le procédé comprenant une étape de décompos it ion de Choleeky pour cal aule r une matrice U telle que: L = U x U. o L est une matrice d'égalisation et U une matrice triangulaire supérieure. L'étape de décomposition de Cholesky comprend: - - ure étape pour ca'culer les coeficients diagonaux uii de la matrice U selon la formule: r i-1 i,i =;11..,.' - '^,i 6;) k=1 . _- le c ieni_ diagonal ce 1.e de rang i de la matrice L et uk,i est le coefficient de la ligne 1- ce rang k et de la colonne de rang i de la matrice U. et - une étape pour calculer les coefficients non diagonaux u;j de la matrice sous la forme: i-1 * - i k=1 ui,i = ui,i l'étape pour calauler les coefficients non diagonaux comprenant - une étape de chargement des coefficients u k,i sur des ertrées de retour de cellules MAC montéss en série, chaque cellule MAC comprenant une entrée d'accumulation, une sortie d'accumulation, une entrce de retour, une sortie de retour et une entrée de coefficient, la sortie d' accumulation de la cellule B 13857.3/PR de rang p (p=1 à N) étant reliée à l'entrée d' accumulation de la cellule de rang p+1, - une étape de chargement des coefficients Uk,j sur les entrées de coefficient des cellules MAC, - une étape de chargement des coefficients uii sur une entrée de coefficient de division dune cellule diagonale, la cellule diagonale (DIAG) comprenant une entrée de coefficient de soustraction, une sortie de résultat, une entrée de retour, une sortie de retour et une entrée de coefficient de division, l'entrée d' accumulation de la cellule diagonale étant reliée à la sortie d'accumulatior de la cellule MAC de rang N. - _ - une étape de chargement des coefic ents li su - _te--r_e G_ ^_CtUr de la celTul- c ac-nae Selon le premier mode de réalisarion de -nvention, es coeicients u,i sort chargés en parallèle sur les entrées de retour des cellules MAC. Selon le deuxième mode de réalisation de l' invention, les coefficients u k,i sont chargés successivement sur les entrces de retour des cellules MAC, à partir de la cellule MAC de rang N. Selon une caractéristique supplémentaire de l' invention, le procédé comprend une étape de calcul du vecteur estimé d qui correspond au vecteur d des données émises par les utilisateurs. L'étape de calcul du vecteur estimé d comprend: - une première étape ae rcaolution de système linéaire triangulaire par la réaurrence B 13857.3/PR Yi = H Zi - UlkYk k=1 pour la résolution de l'équation z=UH.y, avec y=U.d, et - une deuxTème étape de résolution de système linéaire triangulaire par la réaurrence N di = 1 Yi - Uikdk k=i+1 pour la résolution de l'équation y=U.d. Avantageusement, le procédé de détection conjointe selon l'invençicr me. en oeuvre les opérations de décomposition de Co'-sky et de réso ution ae s-ys_è;-.,e _-' angua -e 'a e des m_,.es circuits, permettant ainsi de diminuer la complex té matérielle du dispositif. Des moyens de contrôle permettent de contrôler le chargement séquentiel des différents circuits (cellules MAC, cellule diagonale, multiplexcurs) par les coefficients appropriés. Dans les équations ci-dessus, les coefficients uii, uij, a, aa, b, accu sont généralement des nombres complexes. Il faut cependant noter que l' invention concerne également le cas o ces coefficients sont des nombres récls. Par ailleurs, il faut noter que la matrice U est une matrice bande. Le calcul des coefficients ui,i et ui,j mentionné ci-dessus n'est donc effectué que pour des coefficients uk,i et uk,j non nuls. Dans la suite de la description, les coefficients uk,i et Uk,j non nuls B 13857.3/PR seront préférentiellement notés ukp,i et Ukpj avec p=l, 2,...., N. Brèves descrintion des fiqures D'autres caractéristiques et avantages de l' invention apparaîtront à la lecture d'un mode de réalisation préférentiel fait en référence aux figures jointes parmi lesquelles - les figures 1A et 1B représentent un premier type de cellule utilisce dans le dispositif de détection selon l' invention; - la figure 2 illustre un exemple de réalisation du premier type de cellule représentée aux figures A et 1B; - les figures 3A et 3B représentent ur deuxTème type de cellule utilisée dans le dispositif de actectior selon l'invention, - la figure 4 représente un premier mode de réalisation du dispositif de détection selon l' invention, - la figure 5 représente un deuxTème mode de réalisation du dispositif de détection selon l'invertion, - la figure 6 représente la mise en _uvre de l'étape de décomposition de Choleeky à l' aide d'un dispositif de détection selon le premier mode de réalisation de l' invention, - la figure 7 représente la mise en uvre de l'étape de décomposition de Cholesky à l' aide d'un dispositif de détection selon le deuxième mode de réalisation de l' invention, B 13857.3/PR 13 2838582 - la figure 8 représente un exemple de mise en _uvre de l'étape de décomposition de Cholesky à l' aide d'un dispositif de détection selon le premier mode de réalisation de l' invention, - la figure 9 représente la mise en _uvre d'une première étape de résolution de système triansulaire à l' aide d'un dispositif de détection selon le premier mode de réalisation de l'invention, - la figure 10 représente la mise en _uvre d'une deuxième étape de résolution de système triangulaire qui succède à la première étape de résolution de système triangulaire à l' aide d'un dispositif de détection selon le premier mode de réalisation de l' invention. Sur toues les figures, les mêmes références désignent les mêmes él^m-rts. DescriPtion détaillee de modes de mise en _uvre de l' invention Les figures 1A et 1B représentent un premier type de cellule utilisée dans le dispositif de détection selon l' invention. Ce premier type de cellule, appelée cellule MAC (MAC pour " Multiplication ACcumulation "), est munie de trois entrées et de deux sorties. Les figures 1A et 1B représentent respectivement les signaux entrant et sortant d'une cellule MAC aux instants t et t+l. Si à l' instant t, les trois entrées reçoivent respectivement les grandeurs accu, x et a (cf. figure A), alors, à l' instant t+l, les deux sorties délivrent les grandeurs x et accu+ax (cf. figure 1B). Pour des raisons de B 13857.3/PR commodité, dans la suite de la description, les entrées/sorties d'une cellule MAC seront nommées comme suit: - l'entrée qui reçoit à l' instant t la grandeur accu sera nommoe entrée d' accumulation de cellule MAC, - l'entrée qui reçoit à l' instant t la grandeur a sera nommée entrée de chargement de coefficient de cellule MAC, - ltentrée qui reçoit à l' instant t la grandeur x sera nommée entrée de retour de cellule MAC, - la sortie qui délivre à l' instant t+1 la grandeur accu+ax sera nommée sortie daccumulation de cellule îê'C, - la sortie qui délivre à l' instant t+1 la grandeur x sera nommée sorie de retour de cellule FLC. La tigu--e -e"résene u-. exem;c'e de celule MAC utilisée dans le dispositif de détection selon l' invention. La cellule comprend deux opérateurs de retard 1, 2, un multiplicateur 3 et un additionneur 4. A l' instant t, un coefficient a est appliqué sur l'entrée de chargement de coefficient, une donnce x(t) entre dans l'opérateur de retard 1 et une donnée accu(t) entre dans l'additionneur 4. Le multiplicateur 3 multiplie le coefficient a et la donnce x(t). L'additionneur 4 somme la donnée accu(t) et la donnée ax(t) délivrce par le multiplicateur 3. La donnée accu(t) + ax(t) qui entre dans l'opérateur de retard 2 à l' instant t en sort à l' instant t+1 et la donnée x(t) qui entre dans l'opérateur de retard 1 à l' instant t en ressort à l' instant t+1. B 13857.3/PR Les figures 3A et 3B représentent un deuxTème type de cellule utilisse dans le dispositif de détection selon l' invention; Ce deuxième type de cellule, appelée cellule diagonale, est également munie de trois entrées et de deux sorties. Les figures 3A et 3B représentent respectivement les signaux entrant et sortant d'une cellule diagonale aux instants t et t+1. Si à l' instant t, les trois entrées reçoivent respectivement les grandeurs accu, b et ad (cf. figure 3A), alors, à l' instant t+1, les deux sorties délivrent la grandeur (b-accu)/ad (cf. figure 3B). Pour des raisons de commodité, dans la suite de la description, les entrées/sorties d'une cellule diagonale seront nommoes _ c^rme s_it: l'entrée qui reçoit à l' instant t la grandeur accu sera nommée entrce d'accumulation de cellule diagonale, - l'entrée qui reçoit à l' instant t la grandeur ad sera nommée entrée de chargement du coefficient ad ou entrce de coefficient de division de cellule diagonale, - l'entrée qui reçoit à l' instant t la grandeur b sera nommoe entrce de coefficient de soustraction de cellule diagonale, - les sorties qui délivrent à l' instant t+1 la grandeur (b-accu)/ad seront nommoes, l'une, sortie de résultat de cellule diagonale et, l'autre, sortie de retour de cellule diagonale, selon le sens de cTrculation des données, comme cela sera précisé ultérieurement. B 13857.3/PR La figure 4 représente un premier mode de réalisation de dispositif de détection selon l' invention. Le dispositif de détection selon le premier mode de réalisation de l' invention comprend un circuit pour mettre en _uvre l'équation (1), un ensemble de N cellules MAC, MAC1, MAC2, MACN, une cellule diagonale DIAG et un ensemble de N multiplexcurs M1, M2 À., MN. Aux fins du calcul de l'équation (1), le circuit 5 comprend, par exemple, un opérateur d'élévation au carré, un additionneur, un soustracteur et un opérateur de calaul de racine carrée. L' entrée d' accumulation de la cellule MAC e rang p (p=1, 2,..., N) est relice à la sortie d' accumulation de la cellule rc de rans p-1 et la sortie d'accumulation de la cellule MAC de rang N est reliée à l'entrée d' accumulation de la cellule diagonale. Les entrces et sorties de retour des cellules MAC sont reliées par l'intermédiaire de multiplexeurs. Chaque multiplexeur a deux entrces et une sortie. L'entrée de retour de la cellule MAC de rang p (p=1, 2,..., N-1) est alors reliée à la sortie du multiplexcur de rang p dont une des deux entrées est reliée à la sortie de retour de la cellule MAC de rang p+l. L'entrce de retour de la cellule MAC de rang N est relice à la sortie du multiplexcur de rang N dont une des deux entrées est relice à la sortie de retour de la cellule diagonale. La figure 5 eprésente un deuxième mode de réalisation de dispositif de détection selon l' invention. B 13857.3/PR Le dispositif de détection selon le deuxtème mode de réalisation de l' invention comprend un circuit pour mettre en _uvre l'équation (1), un ensemble de N cellules MAC, MAC1, MAC2, MACN, une cellule diagonale DIAG et un multiplexeur MN ayant chacun deux entrces et une sortie. Les entrées/sorties d' accumulation des cellules MAC sont reliées entre elles comme dans le cas de la figure 4. La sortie de retour de la cellule MAC de rang p (p=1, 2,..., N) est ici reliée à l'entrce de retour de la cellule MAC de rang p-1. Le multiplexeur MN a sa sor_ie re' ée à l'entrce de retour de la cellule MAC de rang N et une première entrse reli=^ à 'a sorti de retour de la cellule diagonale. -- Cor me cela a été mentionné pré_édemment, La décomposition de Choleeky consiste à calculer les coefficients uij d'une matrice U telle que: L UH x U La matrice U est une matrice triangulaire supérieure. Le coefficient diagonal de la ligne i s'écrit: I i-1 Ui, i = ili, i - luk, il ( 1) et les coefficients non diagonaux de la ligne i s'écrivent: li j - uk i.uk j uii uki,=i1 (2) Selon l'invention, les coefficients de la matrice U sont calaulés ligne par ligne, de la première 13857.3/PR ligne à la dernière ligne, le coefficient diagonal d'une ligne étant calaulé avant les coefficients non diagonaux de la ligne. A titre d'exemple non limitatif, pour une matrice U de taille 3x3, les coefficients peuvent successivement être calaulés comme suit: - ull, puis ul2, puis ul3 pour la première ligne, puis - u22, puis u23 pour la deuxième ligne, puis - U33 pour la troisième ligne. Les coefficients diagonaux uii sont calculés à l' aide du circuit 5 qui met en _uvre l'équation (1). Le principe du calaul des coefficients non diagonaux ui,j de la matrice U va maintenant être décrit en référence aux figures 6 et 7. Lors du calaul de l'élément diagonal ui,j, les ^5 entrées de re_ou:' des cellu'es MAC sont cha--éss avec les coefficients uk,i. Le flot de retour est figé pendant le calcul de l' ensemble des coefficients uij d'une même ligne puisque les mêmes uk,i sont nécessaires au calcul de tous ces coefficients. Deux modes de réalisation sont possibles pour * le chargement des coefficients uk i. La figure 6 représente un premier mode de chargement des coefficients uk, i. Le premier mode de chargement des coefficients uk,i est le mode de chagement associé au premier mode de réalisation du dispositif de l' invention. Les multiplexeurs commandés par un premier état de commande permettent un chargement en parallèle des B 13857.3/PR Uk i sur les entrées de retour des cellules MAC lors du calcul des coefficients non diagonaux de la matrice U. La figure 7 représente un deuxTème mode de chargement des coefficients uk,i. Le deuxième mode de chargement des coefficients uk,i est le mode de chargement associé au deuxième mode de réalisation du dispositif de l' invention. Selon le deuxième mode de chargement, les uk i sont chargés l'un après l'autre. Le multiplexeur MN commandé par un premier état de commande permet d' appliquer les uk,i sur l' entrse de retour de la cellule MAC de rang N. Se' cn ce deuxième mode de chargemert, le ci;-c-it 5 pour ica'cul des el=.i,ens diagonaux peut être relié au multiplexeur MN comme cela est représenté en figure 7. Le coefficient uki qui est appliqué sur une entrée du circuit 5 est alors également appliqué en entrce d'un opérateur de conjugaison C qui délivre le coefficient u k,i qui est appliqué en entrée du i-1 2 multiplexeur MN. Le calcul de la quantité k |uk,i| est alors réalisé pendant le chagement des registres du flot de retour. Une fois ce chargement terminé, le flot est figé pour calauler des éléments non diagonaux. De facon plus générale, l'ordre dans lequel sont appliqués les coefficients u ki sur les différentes entrées de retour des cellules MAC n'a pas de réelle importance. En effet les produits ukiukj peuvent être calculés dans n'importe quel ordre avant d'être sommés. B 13857.3/PR La seule condition à réaliser pour le chargement des coefficients u ki sur les entrées de retour des cellules MAC est le maintien du même indice i jusqu'à temps que i-1 * li, i uk, i Uk, i la quantité uij = k=1 soit calaulée. ui,i Une condition à réaliser pour le calcul du coefficient ui,j est d'appliquer successivement, sur 1'entrée de coefficient de la cellule MAC qui rec,oit le coefficient u ki sur son entrée de retour, l' ensemble des coefficients uk,j, j pouvant prendre toutes les valeurs possibles. Les coefficierts diagonaux et non diagonaux de la matrice U calaulés sont stockés au fur et à mesure d _U' d=_a--r -t i C A c tte firl,;e dspositif de détection coujointe selon l' invention comprend des circuits de mémorisation (non représentés sur les figures). Un exemple non limitatif de calaul de coefficient ui,j de matrice U de taille supérieure à 4 va maintenant être donné en description de la figure 8. Cet exemple correspond à une mise en oeuvre de l' invention selon le mode de chargement de coefficients représenté en figure 6. Le coefficient U4s est calaulé dans l'exemple décrit. Le dispositif de détection comprend trois cellules MAC (MACl, MAC2, MAC3) et une cellule diagonale (DIAG). On suppose les entrces de retour des cellules MAC1, MAC2, MAC3 chargées par les coefficients respectifs u 14, U 24, U*34. B 13857.3/PR A l' instant t=1, le coefficient ul5 est appliqué sur l'entrée de coefficient de la cellule MAC1. Cette cellule voit alors sa sortie d' accumulation passer à u*l4ul5 à l' instant t=2. A l ' instant t=2, c' est le coefficient u25 qui est appliqué sur l'entrée de coefficient de la cellule MAC2. Cette cellule voit alors sa sortie d'accumulation passer à u24.u25 + accu o accu est la valeur appliquce sur l' entrée d' accumulation de la cellule MAC2 qui vaut u l4ul5. La sort e d' accumulation de la cellule MAC2 vaut donc u14ul5+u24u2s. A l' instant t=3, l'entrée de coefficient de la cellule MAC3 est alimentée pa U3, La sortie daccumulation de MAC3 passe alors à la valeur p1=U 34.U3s+u 24U-5+u 14-U15- La cellule dia"cnz.le voit donc -'a VG' eu p sur scr en_-e c'a_cumul ation. En z.pp2iquant respect veme-t u e. 145 sur les entrées d'accumulation et de retour de la cellule diagonale, la sortie de résultat de la cellule diagonale passe à U45= 45 P1 qui est la valeur que l'on cherche à u44 calculer. Liaccumulation des termes de Pl progresse vers la droite de la figure (sens des données accumulées). Les cellules MAC laissses sur la gauche ne sont alors plus utilisées. Il est alors avantageusement possible de réaffecter ces cellules au calcul de termes produits destinés aux autres éléments de la ligne courante suivante. A titre d'exemple non limitatif, le tableau ci après décrit la façon dont les coefficients Ukj sont B 13857. 3/PR affectés aux cellules MAC pour le calcul des U45, U47 de lexemple précédent. TABLEAU TEMPS MAC1 MAC2 MAC3 DIAG I DIAG Res | a a a a b t=1 uls t=2 ul6 u2s. t=3 U17 U25 U3; t=4 U27 U36 U44 145 t=5 U37 U44 146 U45 t=6 U44 147 U46 Une fois la matrice U détermince, le procédé selon l'-nvc-icr. cempfend les étapes successives suivantes: - une première étape de résolution de sys.ème linéaire triangulaire par la réaurrence: Yi = H Zi - UikYk pour la résolution de l'équation z=U.y, avec y=U.d, et - une deuxième étape de résolution de système linéaire triangulaire par la récurrence: N di = Uii Yi Uikdk k=i+1 pour la résolution de l'équation y=U.d. B 13857.3/PR Pour les première et deuxième étapes de résolution de système linéaire triangulaire mentionnées ci-dessus, les multiplexeurs sont positionnés de façon que la sortie de retour de la cellule MAC de rang p (p=1, 2,..., N) soit directement reliée à l'entrée de retour de cellule MAC de rang p-l, la sortie de retour de la cellule diagonale étant directement reliée à l'entrée de retour de la cellule MAC de rang N. La cellule diagonale et les cellules MAC constituent alors un réseau tel que célui décrit, par exemple, dans le document intitulé " Introduction to VLSI systems " (Kung, H.T. and Leiserson, C.E., chapter Systolic arrays for VLSI (chap 8.3!' Addison-Wesley, 1980). De fac,on connue en soi, un tel réssau permet de résoudre une équation du type Ax=b o A est une matrice triangulaire connue, b un vecteur connu et x un vecteur inconnu à déterminer. Chaque composante du vecteu solution x est calculé selon la récurrence suivante: x(io) = o k) = xik 1) + aik.xk pour 0 B 13857.3/PR dans le réscau, il accumule les produits calaulés dans chaque cellule MAC. La valeur finale est calculée dans la cellule diagonale. Elle est alors réinjectée dans le réseau pour cTrculer inchangée en sens inverse dans le réseau, permettant ainsi aux éléments suivants x(.) (j > i) d'accumuler les produits ajixi. Plusieurs éléments x(k) sont ainsi calculés en parallèle de manière pipeline. Les entrces et les sorties des cellules MAC entre lesquelles circulent les éléments x(k) calaulés par accumulat ions success ives sont respect ivement l es entrées d' accumulation et les sorties d' accumulation de cellul s > rie même, les e-tréss et les sorties des cellules MAC entre lesquelles circulent les éléments x(k-) réinlectés à pa-rtir de la cellule diagonale sont respectivement les entrces de retour et les sorties de retour des cellules MAC et les entrées des cellules MAC et de la cellule diagonale sur lesquelles sont appliquées les coefficients aik sont les entrces de chargement de coefficient. Les coefficients de la matrice A sont chargés dans le réseau diagonale par diagonale. De ce fait, si A est une matrice bande de largeur Lg. le nombre de cellules MAC du réseau peut être réduit à Lg-1. Ainsi, la taille du réseau n'est-elle pas liée à la taille du système, mais à la largeur de bande de la matrice A. Les figures 9 et 10 représentent, respectivement, le calcul des coefficients Yi de la matrice y et le calcul des coefficients di du vecteur B 13857.3/PR estimé d à l' aide d'un dispositif selon le premier mode de réalisation de l'invention. Pour la première étape de résolution de système linéaire triangulaire, ce sont les coefficients de la matrice transposée de la matrice coujuguce de U qui chargent les entrées de coefficient des cellules MAC et de la cellule diagonale (cf. figure 9). Pour la deuxième étape de résolution de système linéaire triangulaire, ce sont les coefficients de la matrice U qui chargent les entrées de coefficient des cellules MAC et de la cellule diagonale (cf. figure 10). Les figures 9 et 10 concernent un dispositif de cétect_ion se2o.i le premie-' rmcce ce rc.'is=.iori d- l'irventior. Bien aue non représenté sur les figures, il est clair que l' invention concerne également le calcul de la matrice y et le calcul du vecteur d à l' aide d'un dispositif de détection selon le deuxième mode de réalisation de l' invention. Dans le cadre d'une application de l' invention aux téléphones mobiles, le procédé selon l' invention est. par exemple, capable de traiter des donnces transmises selon les modes de communicaticn FDD et TD3 (FDD pour " Frequency Division Duplex " et TDD pour " Time Division Duplex "). Le traitement des données peut avantageusement être effectué en temps réel. B 13857.3/PR REVENDICATIONS 1. Dispositif de détection conjointe de données transmises par étalement de spectre par séquence directe, le dispositif comprenant des moyens de mise en _uvre d'une décomposition de Cholesky pour calculer une matrice U telle que: L = UH x U.. OU L est une matrice d'égalisation et U une matrice triangulaire supérieure, caractérisé en ce que les moyens de mise en _uvre d'une décomposition de Cholesky comprennent: - des moyens pour calauler les coefficients diagonaux ui,i de la matrice U selon la formule: I i-1 uii = gli,i - lukil (1) o k=1 lii est le coefficient diagonal de la ligne de rang i (i=1 à Xi) de la matrice L et uki est le coefficient de la ligne de rang k (k=1 à Xk) et de la colonne de rang i de la matrice U. - des moyens pour calauler les coefficients non diagonaux de la matrice sous la forme: li j - uk, i.Uk, j U1, - ui, i les moyens pour calauler les coefficients non diagonaux comprenant - un réseau de N cellules MAC (MAC1, MAC2,..., MACN) montées en série, chaque cellule MAC comprenant une entrée d'accumulation, une sortie d'accumulation, une B 13857.3/PR entrée de retour, une sortie de retour et une entrée de coefficient, la sortie d'accumulation de la cellule de rang p (p=1 à N) étant reliée à l'entrée d'accumulation de la cellule de rang p+1, - une cellule diagonale (DIAG) comprenant une entrée d' accumulation, une sortie de résultat, une entrée de retour, une sortie de retour et une entrée de coefficient, l'entrce d'accumulation de la cellule diagonale étant reliée à la sortie d' accumulation de la cellule MAC de rang N. - des moyens de chargement des coefficients u k,i sur les entrées de retour des cellules MAC, - des moyens de chargement des coefficients Ukj sur les entrées de coefficient des cellules MAC, - des moyens de chargement des coefficients uii sur l'entrée de coefficient de la cellule diagonale, et - des moyens de chargement des coefficients lij sur l'entrée de retour de la cellule diagonale. 2. Dispositif selon la revendication 1, caractérisé en ce que les moyens de chargement des coefficients u k,i sur les entrées de retour des cellules MAC sont constitués d'un ensemble de N multiplexeurs (M1, M2,..., MN) commandés par un premier état de commande, chaque multiplexeur ayant une première entrée, une deuxième entrée et une sortie, la sortie d'un multiplexeur de rang p (p=1, 2,..., N) étant relice à l'entrée de retour de la cellule MAC de même rang, la première entrée du multiplexeur de rang p étant reliée à un coefficient u k,i différent d'un multiplexeur à l'autre, la deuxième entrée d'un B 13857.3/PR multiplexeur de rang p étant reliée à la sortie de retour de la cellule MAC de rang p+l, la deuxième entrée du multiplexeur de rang N étant reliée à la sortie de retour de la cellule diagonale, le premier état de commande établissant une connexion électrique directe entre la première entrée et la sortie de chaque multiplexeur. 3. Dispositif selon la revendication 1, caractérisé en ce que, la sortie de retour de la cellule MAC de rang p (p=1, 2,..., N) étant reliée à l'entrée de retour de la cellule MAC de rang p-l, les moyens de chargement des coefficients u k,i sur les entrées de retour des cellules MAC sont constitués d'un multiplexeur (MN) commandé par un premier état de commande, le multiplexeur ayant une première entrée, une deuxième entrée et une sortie, la sortie du multiplexeur étant reliée à l'entrée de retour de la cellule MAC de rang N. la première entrée du multiplexcur étant relice, successivement, à un parmi les coefficients u k,i et la deuxième entrée du multiplexeur étant reliée à la sortie de retour de la cellule diagonale, le premier état de commande établissant une connexion électrique directe entre la première entrée du multiplexeur et la sortie du multiplexeur. 4.Dispositif selon la revendication 2 ou 3, caractérisé en ce que les cellules MAC, la cellule diagonale et l' ensemble des multiplexcurs commandés par B 13857.3/PR un deuxième état de commande constituent un réseau permettant de calculer un vecteur estimé d tel que: d = L1 AH e, o - L est la matrice d'égalisation, - AH est une matrice symétrique hermitienne d'une matrice A qui est la matrice produit de convolution entre une matrice C de codes d'étalement utilisés pour coder les données à transmettre et une matrice h qui représente la réponse impulsionnelle du canal de transmission, et - e est une matrice élaborce à partir des données reques. 5. Procédé de détection conjointe de données transmises par étalement de spectre par séquence directe, le procédé comprenant une étape de décompos it ion de Choleeky pour calauler une matrice U telle que: L UH x U o L est une matrice d'égalisation et U une matrice triangulaire supérieure, caractérisé en ce que I'étape de décompos it ion de Choleeky comprend: une étape pour calculer les coefficients diagonaux uii de la matrice U selon la formule: I i-1 Uii = li,i - Uk,il (1) o k=1 li,i est le coefficient diagonal de la ligne de rang i de la matrice L et uk,i est le coefficient de la ligne de rang k et de la colonne de rang i de la matrice U. - une étape pour calauler les coefficients non B 13857.3/PR diagonaux de la matrice sous la forme: i-1 * j k=1 ui,i = ui,i l'étape pour calculer les coefficients non diagonaux comprenant - une étape de chargement des coefficients u k,i sur des entrées de retour de cellules MAC montées en série, chaque cellule MAC comprenant une entrée d'accumulation, une sortie d'accumulation, une entrée de retour, une sortie de retour et une entrée de coefficient, la sortie d' accumulation de la cellule de rang p (p=1 à N) étant reliée à l'entrée d'accumulation de la cellule de rang p+1, - une étape de chargement des coefficients uk, jl (j=1, 2,...,N+1) sur les entrées de coefficient des cellules MAC, - une étape de chargement du coefficient ui,i sur l'entrée de coefficient d'une cellule diagonale, la cellule diagonale (DIAG) comprenant une entrée d' accumulation, une sortie de résultat, une entrée de retour, une sortie de retour et une entrée de coefficient, l'entrce d' accumulation de la cellule diagonale étant reliée à la sortie d' accumulation de la cellule MAC de rang N. et - une étape de chargement du coefficient li,j sur l'entrée de retour de la cellule diagonale, 6. Procédé selon la revendication 5, caractérisé en ce que les coefficients u k,i sont B 13857.3/PR chargés en parallèle sur les entrées de retour des cellules MAC (MAC1, MAC2,..., MACN). 7. Procédé selon la revendication 5, caractérisé en ce que les coefficients u k,i sont chargés successivement sur les entrces de retour des cellules MAC (MAC1, MAC2,..., MACN). 8. Procédé selon l'une quelconque des revendications 5 à 7, caractérisé en ce qu'il comprend une étape supplémentaire de calaul d'un vecteur de données estimées. 9. Procédé selon la revendication 8, caractérisé en ce que l'étape supplémentaire comprend: - une première étape de résolution de système linéaire triangulaire par la récurrence Yi = H Zi - UikYk1 pour la résolution de l'équation z=UH.y, avec y=U.d o d est le vecteur de données estimées défini par l'équation: d = L1 AH e, o - L est la matrice d'égalisation, - AH est une matrice symétrique hermitienne d'une matrice A qui est la matrice produit de convolution entre une matrice C de codes d'étalement utilisés pour coder le s données à transmet t re et une mat ri ce h qui représente la réponse impusionnelle du canal de transmission, B 13857.3/PR - e est une matrice élaborée à partir des données re,cues, et - une deuxième étape de résolution de système linéaire triangulaire par la réaurrence di = - pi - Uikdk k=i+1 pour la resolution de l'équation y=U.d.