FR 2894416 A1 20070608 FR 0512311 A 20051205 PROCEDE ET DISPOSITIF DE SELECTION DES PARAMETRES D'ETALEMENT D'UN SYSTEME OFDM-CDMA La présente invention se rapporte à un procédé et à un dispositif de sélection des paramètres d'étalement d'un système OFDM-CDMA (multiplexage par division de fréquences orthogonales - accès multiple à répartition par codes, en anglais "Orthogonal Frequency Division Multiplexing - Code Division Multiple Access"). Plus précisément, l'invention a pour objet la sélection des paramètres de l'étalement à deux dimensions pour le lien descendant d'un système multi-porteuses avec un accès multiple de type CDMA. Il s'agit des paramètres d'étalement temporel et fréquentiel, notés respectivement ST et SF. L'invention trouve en particulier une application dans le domaine des télécommunications, pour des systèmes de communication exploitant la technique d'accès multiple CDMA combinée à une transmission OFDM, tels que par exemple les systèmes MC-CDMA (accès multiple à répartition par codes - multi-porteuses, en anglais "Multi-Carrier - Code Division Multiple Access"), MC-DS-CDMA (accès multiple à répartition par codes - multi-porteuses - séquence directe, en anglais "Multi-Carrier - Direct Sequence - Code Division Multiple Access") ou les systèmes à étalement temps/fréquence à deux dimensions (2D). Dans la suite, le terme générique OFDM-CDMA désigne un système utilisant un étalement temps/fréquence. La figure 1 illustre l'architecture générale d'un émetteur numérique 25 OFDM-CDMA classique pour une liaison descendante d'une station de base vers un terminal mobile. Les bits ao(n), ... , aK_l (n) des K utilisateurs sont d'abord envoyés à un module de codage de canal respectif 10o, ..., 10K_1. Le module de codage de canal comprend un encodeur, qui peut être par exemple un codeur de type 30 convolutif, un turbocodeur ou un codeur du type LDPC ("Low-Density Parity- Check Code"). Le module de codage de canal comprend en outre un dispositif de poinçonnage et un entrelaceur bit. Il fournit en sortie des données bo(n), ..., bK-1(n). Les données binaires bo(n), ..., bK_1(n) sont ensuite envoyées à une unité de modulation I/Q respective 120, ..., 12K_1, qui réalise la modulation des bits encodés, par exemple une modulation d'amplitude en quadrature à 4 états MAQ-4 et fournit en sortie des données do(n), ..., dK_1(n). Les données do(n), ..., dK_l(n) des différents utilisateurs sont ensuite traitées par un module d'étalement 14. Puis les signaux étalés sont traités par un module d'allocation des chips 16, qui les place sur une grille temps/fréquence. Le signal résultant est alors transmis à un modulateur OFDM, qui comprend successivement une unité de conversion série-parallèle 18, qui fournit en sortie des données xo(n), ..., xN_l(n), N étant le nombre de sous-porteuses, une unité de transformation de Fourier rapide inverse 20, une unité de conversion parallèle-série 22 et une unité d'insertion de préfixe cyclique 24. Les symboles obtenus sont émis sur un canal de transmission. Pour une description des techniques classiques de I'OFDM, on se reportera utilement à l'article de W. ZHENDAO et G. B. GIANNAKIS intitulé "Wireless Multicarrier Communications - Where Fourier meets Shannon", IEEE Signal Processing Magazine, vol. 17, n 3, pages 29 à 48, mai 2000. La figure 2 détaille le module d'étalement 14 de l'émetteur. A chaque symbole dk(n) de l'utilisateur k (k = 0, ..., K-1) est d'abord affectée une amplitude 'Pk par un module 26k. Puis on procède dans un module 28k à une élévation de cadence d'un facteur L et enfin à un filtrage numérique ck(z) dans un module 30k, les coefficients de filtrage étant égaux aux chips de la séquence d'étalement de l'utilisateur k. Chaque symbole dk(n) de l'utilisateur k est étalé par une séquence de L chips ck(z). Les signaux étalés de tous les utilisateurs sont ensuite additionnés par un module additionneur 32. Comme décrit plus haut, le module d'allocation des chips 16 répartit les échantillons provenant du module d'étalement 14 sur une grille temps/fréquence. On suppose que le facteur d'étalement L est égal à SFXST, où SF est le paramètre d'étalement dans le domaine fréquentiel et ST est le paramètre d'étalement dans le domaine temporel. Les figures 3 à 5 représentent des grilles temps-fréquence respectivement dans des cas où ST = 1 (figure 3), SF = 1 (figure 4) et SF et ST sont quelconques (figure 5). La figure 3 illustre le cas où ST = 1. Dans ce cas, on retrouve les caractéristiques d'un système MC-CDMA classique. On transmet F = N/SF symboles de données par code dans un symbole OFDM. La figure 4 illustre le cas où SF = 1. Dans ce cas, on retrouve les 10 caractéristiques d'un système MC-DS-CDMA classique. On transmet F = N symboles de données par code sur ST symboles OFDM. La figure 5 illustre le cas où SF et ST sont quelconques. Dans ce cas, on transmet F = N/SF symboles de données par code sur ST symboles OFDM. 15 Le module d'allocation des chips 16 représenté sur la figure 1 fournit en sortie un vecteur de taille N qui peut s'exprimer suivant l'équation (1), la taille de ce vecteur correspondant à la taille de la transformation de Fourier rapide : K-1 xi[gSF + p] = 1VPkdk[q]ck[pST + i] (1) k=0 où : 20 • i = 0, ..., ST-1 correspond à l'indice du symbole OFDM, • q = 0, ..., F-1 correspond à l'indice de la sous-bande, • p = 0, ..., SF-1 correspond à l'indice de la sous-porteuse, • SF est le facteur d'étalement dans le domaine fréquentiel, • ST est le facteur d'étalement dans le domaine temporel, 25 • K est le nombre total d'utilisateurs, Pk est la puissance associée au k1eme utilisateur, • • les dk sont les symboles associés au kième utilisateur, et • ck est le code d'étalement pour le kième utilisateur. Après transposition dans le domaine temporel par l'intermédiaire de 30 l'unité de transformation de Fourier rapide inverse 20, un préfixe cyclique est ajouté par l'unité d'insertion de préfixe cyclique 24. Il contient NG > W-1 échantillons, où W est la durée maximale de la réponse impulsionnelle du canal global. La figure 6 présente la structure d'un récepteur numérique OFDMCDMA classique correspondant à l'émetteur illustré sur la figure 1. Un module de synchronisation grossière 60 effectue, d'une part, la détection du début du symbole OFDM et, d'autre part, l'estimation initiale du décalage LF entre les fréquences porteuses de l'émetteur et du récepteur et du décalage AT entre les horloges d'échantillonnage de l'émetteur et du récepteur. Cette estimation est dite grossière dans le sens où la variance de l'estimation est élevée. Cela permet d'entrer dans la plage de fonctionnement du récepteur, mais nécessite une correction dite fine, pour atteindre les performances désirées. Ainsi, la synchronisation est réalisée en deux phases. Après une synchronisation grossière, le préfixe cyclique est éliminé par un module 62 de suppression de préfixe cyclique et le signal est conditionné en vecteurs de N échantillons r;(m), m = 0, ..., N-1, l'indice i indiquant le numéro du symbole OFDM reçu. Dans la suite, on présente les traitements effectués pour décoder les symboles de l'utilisateur k = O. Les résultats obtenus sont supposés identiques pour les autres utilisateurs. Les échantillons ri(m) sont ensuite fournis à une unité de conversion série-parallèle 64 puis à une unité de transformation de Fourier rapide 66 et une unité d'estimation de canal 68 elle-même reliée à une unité d'estimation des conditions de propagation du canal 70. L'unité d'estimation de canal 68 est également reliée à un égaliseur linéaire 72 qui applique un coefficient g;[gSF+p] par sous-porteuse. L'égaliseur 72 est relié à un module de corrélation 74, qui fournit en sortie les symboles estimés do (n),...,dK_, (n) . On supposera que les coefficients de I'égaliseur linéaire sont calculés indépendamment des codes d'étalement. Il peut s'agir des égaliseurs linéaires MRC (combinaison à rapport maximum, en anglais "Maximum Ratio Combining"), EGC (combinaison à gain identique, en anglais "Equal Gain Combining"), ZF (forçage à zéro, en anglais "Zero Forcing") ou MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Error"), bien connus de l'homme du métier. Les symboles estimés obtenus en sortie du module de corrélation 74 peuvent s'exprimer suivant l'équation (2) : SF-1ST-1 do[q] _ Igi[gSF + p]Co[pST + i]xi[gSF + p]hi[gS p=0 i=0 SF-1ST-1 + lgi[gSF + p]Co[pST + i]wi[gSF + p] p=0 i=0 où: • gi[gsF+p] est le coefficient appliqué par l'égaliseur à la pième sous-porteuse de la gième sous-bande du iième symbole OFDM, • co est le code d'étalement de l'utilisateur 0, 10 • le signe * désigne le complexe conjugué, • xi[gSF+p] correspond aux symboles fournis en sortie par l'unité de transformation de Fourier rapide 66 sur la pième sous-porteuse de la gième sous-bande du i1ème symbole OFDM, • hi[gSF+p] correspond à l'atténuation du canal pour la pième sous-porteuse 15 de la gième sous-bande du i1ème symbole OFDM, et • wi[gSF+p] correspond à l'échantillon de bruit blanc additif gaussien de variance Q-2 sur la pième sous-porteuse de la gième sous-bande du iième symbole OFDM. Les données sont ensuite envoyées à une unité de démodulation I/Q 20 souple 76. L'opération de démodulation I/Q dite dure consiste à retrouver les valeurs binaires émises à partir des symboles complexes issus du détecteur linéaire. Lorsqu'on utilise un décodeur de canal à entrées souples, les valeurs optimales à injecter dans le décodeur de canal sont des valeurs souples. On parle de valeurs souples dans le sens où ce ne sont pas directement les 25 valeurs dures "0" ou "1". Ainsi, si on utilise une modulation MAQ-16 (4 bits), la démodulation I/Q souple consiste à calculer 4 valeurs souples correspondant aux 4 bits de la modulation MAQ-16. Les valeurs souples (ou métriques) optimales à injecter dans le décodeur de canal à entrées souples correspondent au Logarithme du Rapport de Vraisemblance (LRV). + p] (2) Après l'opération de démodulation I/Q souple vient le processus de décodage de canal, effectué par des modules 780, ..., 78K_1. Ces modules réalisent les opérations de désentrelacement bit, de dépoinçonnage et de décodage (par exemple, par un décodeur de Viterbi dans le cas d'un codeur convolutif à l'émission, par l'algorithme Max Log-MAP dans le cas d'un turbocodeur ou par l'algorithme Min-Sum dans le cas d'un codage LDPC). Cela permet de retrouver les données binaires émises. Dans l'état de la technique, l'unité d'estimation des conditions de propagation du canal 70 permet d'évaluer le Rapport Signal A Interférence + Bruit (RSIB) instantané et/ou l'étalement maximum des retards et/ou la fréquence maximum du Doppler. Ces trois paramètres mesurent les conditions de propagation du canal. Selon les valeurs de ces paramètres renvoyées à l'émetteur par boucle de retour, l'émetteur peut décider d'adapter les paramètres d'étalement SF et ST afin d'améliorer la qualité de transmission. Néanmoins, les méthodes de calcul de l'étalement maximum des retards et de la fréquence Doppler ne sont pas faciles à mettre en oeuvre. Par ailleurs, le RSIB instantané n'est pas un critère optimal dans la sélection des paramètres d'étalement lorsqu'on utilise un processus de codage/décodage de canal. On présente maintenant la notion de RSIB asymptotique. Lorsque les dimensions du système deviennent grandes en terme de nombre et taille des codes d'étalement, la théorie des matrices aléatoires offre des outils d'analyse très puissants permettant d'obtenir des valeurs explicites du RSIB à la sortie de l'égaliseur, tout en tenant compte de la propriété d'orthogonalité des codes d'étalement. Ce régime asymptotique est néanmoins obtenu pour des valeurs courantes des facteurs d'étalement (par exemple, taille 32). L'expression du RSIB asymptotique pour l'utilisateur de rang 0 d'un système OFDM-CDMA pour la sous-bande de rang q est donnée par la formule suivante : RSIBo[g] = SO[g] I[q] + N[q] (3) où So[.] correspond à la puissance du signal utile après égalisation et corrélation, I[.] correspond à la puissance du bruit d'accès multiple créé par les autres codes d'étalement et N[.] correspond à la puissance du bruit gaussien filtré par l'égaliseur et la séquence d'étalement, avec : 2 So[q] = po 1 SF-1ST-1 ù 1gi[gSF + p]hi[gSF + p] L p=o i=0 (4) où: • po est la puissance appliquée à l'utilisateur de rang 0, • L = SFXST, • gi[gSF+p] est le coefficient d'égalisation pour la pième sous-porteuse de la 10 gième sous-bande pour le iième symbole OFDM, ce coefficient d'égalisation dépendant du canal de transmission, et • hi[gSF+p] correspond à l'atténuation du canal pour la pième sous-porteuse de la gième sous-bande pour le i1ème symbole OFDM, ( SF-1ST-1 ù,+p] p=0 i=0 1 SF-1ST-1 ù E Ig~[gSF + p]hi[gSF + p] ~ L p=0 i=0 I[q] = ap hi[qSF + p] 2 2 2" (5) 15 où : a = K/L est la charge du système (K étant le nombre total d'utilisateurs), et K-1 p = ~pk est la puissance moyenne des codes interférents, et K -1 k=1 2 1 SF-1ST-1 2 N[q] = a ù g [qSF + p] (6) L p=o i=o 20 où 62 est la variance du bruit. Si le système considéré est un système multicellulaire, la variance du bruit a2 contient la puissance du bruit thermique No ainsi que la puissance reçue par l'utilisateur de rang 0 provenant des autres stations de base interférentes. 25 Ces différentes valeurs sont obtenues indépendamment de la valeur des codes d'étalement, tout en tenant compte de la propriété d'orthogonalité des codes. On montre que le RSIB asymptotique permet d'estimer précisément et facilement le Taux d'Erreur Binaire (TEB) avant codage canal. Cependant, un système réel emploie un codeur de canal et il convient de prendre en compte son effet sur les performances du système. L'article de N. MAEDA, Y. KISHIYAMA, H. ATARASHI et M. SAWAHASHI intitulé "Variable Spreading Factor OFCDM with two dimensional spreading that prioritizes time domain spreading for forward link broadband wireless access", VTC Spring 2003, pages 127 à 132, Jeju Island, Corée, avril 2003, met en évidence l'impact du schéma de modulation et de codage sur les performances (TEB) d'un système OFDM-CDMA, en fonction des paramètres du canal de propagation (bande et temps de cohérence du canal) et des paramètres d'étalement SF et ST. De plus, dans l'article de N. MAEDA, Y. KISHIYAMA, K. HIGUCHI, H. ATARASHI et M. SAWAHASHI intitulé "Experimental Evaluation of Throughput Performance in Broadband Packet Wireless Access Based on VSFOFCDM and VSD-CDMA", IEEE PIMRC 2003, pages 6 à 11, l'utilisation du RSIB instantané mesuré par la formule (3) ci-dessus pour sélectionner les paramètres d'étalement SF et ST est peu avantageuse, car ce RSIB ne dépend pas des paramètres d'étalement, puisqu'il est obtenu en mesurant la puissance du signal utile et la puissance des interférences et du bruit thermique. Il ne permet donc pas de sélectionner directement les paramètres d'étalement SF et ST optimaux lorsque le système utilise un module de codage/décodage de canal. En effet, il existe plusieurs valeurs de RSIB représentant un même bloc de données codées entrant dans le décodeur de canal. Ces multiples valeurs de RSIB sont dues au fait que le canal de propagation varie en temps et en fréquence à l'intérieur d'un bloc d'information codée. II est par conséquent inefficace d'optimiser le RSIB donné par la formule (3) car il faut tenir compte du codage canal et de la modulation. L'invention a pour but de remédier aux inconvénients de l'art antérieur. Dans ce but, la présente invention propose un procédé de sélection des paramètres d'étalement fréquentiel (SF) et/ou temporel (ST) pour un système de communication incluant un émetteur et un récepteur et utilisant une modulation d'amplitude en quadrature et mettant en oeuvre une technique d'accès multiple à répartition par codes ou CDMA et un multiplexage par division de fréquences orthogonales ou OFDM, dans lequel le domaine des fréquences est découpé en F = N/SF sous-bandes de SF sous-porteuses, N étant la taille de la transformation de Fourier rapide utilisée par le récepteur et SF étant le paramètre d'étalement fréquentiel, et le domaine temporel est découpé en T = Ns/ST blocs de ST symboles OFDM, Ns étant le nombre de symboles OFDM utilisés pour coder une trame de données et ST étant le paramètre d'étalement temporel, ce procédé étant remarquable en ce qu'il comporte des étapes suivant lesquelles : on estime les coefficients de l'atténuation du canal de propagation ; - on détermine, pour l'ensemble des F sous-bandes de SF sous-porteuses et des T blocs de ST symboles OFDM, au moins un Rapport Signal à Interférence + Bruit (RSIB) asymptotique, à partir des coefficients de l'atténuation du canal de propagation et du type d'égaliseur utilisé par le récepteur ; on en déduit, pour chaque bloc d'information codée obtenu à l'issue de la modulation, le Rapport Signal à Interférence + Bruit asymptotique compressé, qui s'exprime par une formule analytique dépendant des paramètres d'étalement fréquentiel (SF) et temporel (ST) ; et - on sélectionne les paramètres d'étalement fréquentiel et temporel qui maximisent la valeur du Rapport Signal à Interférence + Bruit asymptotique compressé. Ainsi, l'optimisation des paramètres d'étalement fréquentiel SF et temporel ST se fait par l'intermédiaire du calcul du RSIB asymptotique compressé. Cette optimisation permet d'améliorer la qualité de la transmission entre l'émetteur et le récepteur selon les caractéristiques du canal de propagation, à savoir, la bande de cohérence B,oh du canal et le temps de cohérence Tmh du canal. Le RSIB compressé peut par ailleurs être utilisé pour estimer le taux d'erreur binaire ou le taux d'erreur paquet au niveau du récepteur, ce qui indique la qualité de la communication. Avantageusement, on détermine plusieurs RSIB asymptotiques et en particulier, on détermine tous les RSIB asymptotiques y(i, j) qui figurent dans la formule donnée ci-après du RSIB asymptotique compressé, noté RSIBcompressé. Cela permet d'obtenir une sélection plus précise des paramètres d'étalement et d'améliorer de ce fait la transmission. En particulier, lorsque les valeurs du RSIB asymptotique y(i,j) varient beaucoup à l'intérieur du bloc d'information codée obtenu à l'issue de la modulation, il est avantageux d'utiliser toutes ces valeurs (qui sont donc au nombre de FxT) pour calculer le RSIB asymptotique compressé. En revanche, lorsque les valeurs du RSIB asymptotique y(i,j) varient peu à l'intérieur du bloc d'information codée, il est avantageux de n'utiliser qu'un faible nombre de ces valeurs. Dans un mode particulier de réalisation, à l'étape de calcul du RSIB asymptotique compressé, on calcule le RSIB asymptotique compressé au moyen de la formule suivante : ~ L F-1 T-1 _Y(i.j) RSIBcompressé = ù%`'In NN • g i=0 j=0 ou : • a, est un scalaire connu à l'avance par le récepteur, dépendant de la modulation et du codage, dont le mode d'obtention sera décrit en détail plus loin, • L est le facteur d'étalement : L = SFXST, • i et j sont des entiers positifs, et • y(i, j) est le RSIB asymptotique de la i'eme sous-bande du j'erre bloc de ST symboles OFDM. Dans un mode particulier de réalisation, à l'étape de calcul du RSIB asymptotique compressé, on calcule le RSIB asymptotique compressé de façon itérative en utilisant l'approximation suivante : In(ea + eb ) , max(a, b) + fc û b ) où fcù b ) est une table pré-calculée ne dépendant que de a ù b ) et approchant la fonction In(1 + e-la-bl )Cela permet de simplifier le calcul du RSIB asymptotique compressé, en évitant d'effectuer des opérations d'exponentiation et de calcul de logarithmes. Dans un mode particulier de réalisation, l'étape de sélection consiste à procéder à une recherche exhaustive du couple de paramètres d'étalement fréquentiel et temporel qui maximise la valeur du RSIB asymptotique compressé. Ce mode de réalisation est particulièrement simple à mettre en oeuvre lorsque le nombre de valeurs du couple de paramètres d'étalement est faible, par exemple, de l'ordre de 10. En variante, l'étape de sélection consiste à appliquer une méthode de recherche par la section Golden. Cette variante est particulièrement avantageuse lorsque l'utilisateur se déplace lentement, par exemple à une vitesse de l'ordre de 3 km/h, ou dans le cas où le canal varie peu en fréquence. Selon une autre variante, l'étape de sélection consiste à appliquer une méthode d'optimisation sous contrainte utilisant la technique du multiplicateur de Lagrange, où la contrainte s'écrit g(SF,ST) = SFST - L, SF étant le paramètre d'étalement fréquentiel, ST étant le paramètre d'étalement temporel et L étant le facteur d'étalement. Cette variante permet d'obtenir d'excellents résultats en termes de qualité. II est à noter que l'ensemble des modes de réalisation présentés ci-dessus sont applicables quel que soit le type d'égaliseur utilisé. Néanmoins, l'invention propose en outre divers modes particuliers de réalisation dans des cas particuliers d'égaliseurs. Dans un mode particulier de réalisation, lorsque le procédé conforme à l'invention est mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Error') mono-utilisateur, le RSIB asymptotique y(i, j) de la i1ème sous- bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante : po ' n2(i,j) hjST+k [i SF + p] hjST+k[i'SF +p] 2 SF-1ST-1 L p=o k=o hjST+k[j'SF +p] 2 \ 2 hjsT+k[i'SF +p] + ~ 4 2 2 -T12(i,j) ap où 2 1 SF-1ST-1 III, j) = ùIl Il p=0 k=0hjsT+k[j cF +p] hj.sT+k[i' SF + p] 2 + où: • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = KIL est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents et les pi sont les puissances des codes d'étalement interférents, avec i # 0, • p et k sont des entiers positifs, • hjsT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i SF+p)ième sous-porteuse du (j.ST+k)ième symbole OFDM, et • a2 est la variance du bruit. Pour un égaliseur MMSE mono-utilisateur et dans le cas particulier où K-1 1 K-1 K lp; = ù~p; , le RSIB asymptotique est défini par la relation suivante : Y(i, j) = po ' r1(i, j) ap(1û rl(i, j)) SF-1ST-1 n(i, _ ù L p=0 k=0 13 h;.sT+k [i . SF + p] 2 2 + h;.sT+k [i SF + p] K-1 P KZ,=1 Pi avec les mêmes notations que précédemment. Ces modes de réalisation avec égaliseur MMSE permettent de réaliser un compromis entre la réduction de l'interférence d'accès multiple et le rehaussement du niveau de bruit. Dans un autre mode particulier de réalisation, lorsque le procédé conforme à l'invention est mis en oeuvre en liaison avec un récepteur comportant un égaliseur MRC (combinaison à rapport maximum, en anglais "Maximum Ratio Combining"), le RSIB asymptotique y(i, j) de la i1ème sous-bande du j1ème bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante : Y(i,J) = F po r12(i,J) ( S -1ST-1 ~ ap L hisT+k[i SF +p]4 ûr~2(i,j) +62r1(i,j) p=0 k=0 1 SF-1ST-1 2 avec fl(i, J) = L. h,sT+k [i SF + p] p=0 k=0 où: • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, ou • hj.sT+k [i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)ième sous-porteuse du (j.ST+k)ième symbole OFDM, et • a2 est la variance du bruit. Ce mode de réalisation avec égaliseur MRC est particulièrement avantageux dans un contexte mono-utilisateur pour optimiser le rapport signal à bruit en sortie de ('égaliseur. Dans un autre mode particulier de réalisation, lorsque le procédé conforme à l'invention est mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Erro?') multi-utilisateur, le RSIB asymptotique y(i, j) de la iième sous-bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est solution de l'équation implicite suivante : 1 SF-1ST-1 1=-, û p=0 k=0 hj.sT+k [i ' SF + p] 1 A(a,Pk,Y(i, j)) + a' B(a, Pk ,Y(i, j)) 2 a2 où: • L est le facteur d'étalement : L = SFXST, • a = K / L , K étant le nombre de codes d'étalement transmis, • p et k sont des entiers positifs, • hj sT+k [i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)ième sous-porteuse du (j.ST+k)'ème symbole OFDM, et 20 • a2 est la variance du bruit, • pk est la puissance appliquée au kième utilisateur, et : A(a, Pk, Y(i, j)) = (1û a) Y(i, j) + a • Y(i, j) . m(pk, Y(i, j)) Pk B(a,Pk,Y(i,J)) = 1 - Pk m(Pk,Y(i,j)) m(pk, Y(i, J)) = f~ power (u) du Pk + Y(i, J) • u où power (p) est la distribution limite de la puissance lorsque L et K tendent vers l'infini et a = K / L reste constant. Ce mode de réalisation est particulièrement avantageux dans un contexte multi-utilisateurs. Dans un autre mode particulier de réalisation, lorsque le procédé conforme à l'invention est mis en oeuvre en liaison avec un récepteur comportant un égaliseur EGC (combinaison à gain identique, en anglais "Equal Gain Combining"), le RSIB asymptotique y(i, j) de la Mme sous-bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante : Po . T12(i,j) SF 2 hjST+k[i . SF +p] T12(bJ) i + a2 ap 1 sF-IsT-1 avec rl(i,j)=ùE L- p=0 k=0 OU : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • h;.sT+k [i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)ième sous-porteuse du (j.ST+k)ième symbole OFDM, et • 62 est la variance du bruit. Sachant que I'égaliseur EGC ne corrige que la phase du signal et permet de conserver l'amplitude de chaque sous-porteuse après égalisation, ce mode de réalisation permet de conserver l'orthogonalité des codes dans un canal à bruit non limité (contrairement à un égaliseur MRC) et ainsi de diminuer le terme d'interférence. Dans le même but que celui indiqué plus haut, la présente invention propose également un dispositif de sélection des paramètres d'étalement h;.sT+k [i • SF + p] fréquentiel et/ou temporel pour un système de communication incluant un émetteur et un récepteur et utilisant une modulation d'amplitude en quadrature et mettant en oeuvre une technique d'accès multiple à répartition par codes ou CDMA et un multiplexage par division de fréquences orthogonales ou OFDM, dans lequel le domaine des fréquences est découpé en F = N/SF sous-bandes de SF sous-porteuses, N étant la taille de la transformation de Fourier rapide utilisée par le récepteur et SF étant le paramètre d'étalement fréquentiel, et le domaine temporel est découpé en T = NS/ST blocs de ST symboles OFDM, Ns étant le nombre de symboles OFDM utilisés pour coder unetrame de données et ST étant le paramètre d'étalement temporel, ce dispositif étant remarquable en ce qu'il comporte : - un module pour estimer les coefficients de l'atténuation du canal de propagation ; - un module pour déterminer, pour l'ensemble des F sous-bandes de SF sous-porteuses et des T blocs de ST symboles OFDM, au moins un Rapport Signal à Interférence + Bruit (RSIB) asymptotique, à partir des coefficients de l'atténuation du canal de propagation et du type d'égaliseur utilisé par le récepteur ; - un module de calcul pour en déduire, pour chaque bloc d'information codée obtenu à l'issue de la modulation, le Rapport Signal à Interférence + Bruit asymptotique compressé, qui s'exprime par une formule analytique dépendant des paramètres d'étalement fréquentiel et temporel ; et - un module pour sélectionner les paramètres d'étalement fréquentiel et temporel qui maximisent la valeur du Rapport Signal à Interférence + Bruit asymptotique compressé. Toujours dans le même but, la présente invention propose en outre une unité de simulation d'un système radio-mobile adaptée à mettre en oeuvre un procédé de sélection des paramètres d'étalement tel que ci-dessus. Toujours dans le même but, la présente invention propose en outre une unité de simulation d'un système radio-mobile comprenant un dispositif de sélection des paramètres d'étalement tel que ci-dessus. Toujours dans le même but, la présente invention propose également : - un dispositif de téléphonie mobile sans fil adapté à mettre en oeuvre un procédé de sélection des paramètres d'étalement tel que ci-dessus, 5 et -un dispositif de téléphonie mobile sans fil comprenant un dispositif de sélection des paramètres d'étalement tel que ci-dessus. En effet, l'invention peut également s'appliquer à un dispositif de réception tel qu'un téléphone mobile, dans le cadre d'un procédé d'émission et 10 de réception de signaux où des signaux sont transmis par un dispositif transmetteur vers des dispositifs de téléphonie sans fil et où au moins un de ces dispositifs de téléphonie sans fil calcule de la façon décrite plus haut le RSIB asymptotique compressé qui le concerne. Ce dispositif de téléphonie mobile peut ainsi sélectionner les 15 paramètres d'étalement aptes à optimiser le RSIB asymptotique compressé. Les paramètres d'étalement optimaux sont alors renvoyés au dispositif transmetteur afin d'optimiser la transmission pour ce mobile. Cette optimisation des paramètres d'étalement peut être faite directement au niveau du dispositif transmetteur si le mobile renvoie au transmetteur l'information d'atténuation du 20 canal et le type d'égaliseur utilisé. Les caractéristiques particulières et les avantages du dispositif de sélection, des unités de simulation d'un système radio-mobile et des dispositifs de téléphonie mobile étant similaires à ceux du procédé de sélection, ils ne sont pas répétés ici. 25 D'autres aspects et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit de modes particuliers de réalisation, donnés à titre d'exemples non limitatifs. La description se réfère aux dessins qui l'accompagnent, dans lesquels : - la figure 1, déjà décrite, représente de façon schématique un 30 émetteur OFDM-CDMA classique ; - la figure 2, déjà décrite, représente de façon schématique l'architecture du module d'étalement de la figure 1 ; - les figures 3 à 5, déjà décrites, représentent des grilles temps-fréquence pour diverses valeurs des facteurs d'étalement dans les domaines fréquentiel et temporel ; - la figure 6, déjà décrite, représente de façon schématique un 5 récepteur OFDM-CDMA classique ; - la figure 7 est un schéma de principe de la technique EESM ; - la figure 8 est un organigramme illustrant les principales étapes d'un procédé de sélection des paramètres d'étalement conforme à la présente invention, dans un mode particulier de réalisation ; 10 - la figure 9 illustre le découpage en blocs et en sous-bandes d'un "slot", c'est-à-dire du nombre de symboles OFDM utilisés pour coder une trame de données ; - la figure 10 représente de façon schématique dans un mode particulier de réalisation des dispositifs de téléphonie mobile mettant en oeuvre 15 la présente invention ; et - la figure 11 représente de façon schématique dans un mode particulier de réalisation une unité susceptible de comporter un simulateur de système radio-mobile mettant en oeuvre la présente invention. Lors des travaux du 3GPP, une technique appelée "Effective 20 Exponential SIR Mapping" (EESM) a été proposée pour trouver un RSIB compressé qui permette d'estimer le TEB en sortie du décodeur de canal. On se reportera utilement à ce sujet aux documents de la société Ericsson intitulés "Effective SNR mapping for modelling frame error rates in multiple-state channels", 3GPP2-C30-20030429-010 et "System-level evaluation of OFDM - 25 furtherconsiderations", TSG-RAN WG1 #35, Lisbonne, Portugal, 17-21 novembre 2003, R1-031303. Le schéma de principe de I'EESM est illustré sur la figure 7. On suppose qu'une trame codée transporte N symboles MAQ (par exemple MAQ-2, dite aussi BPSK, ou MAQ-4, dite aussi QPSK, ou MAQ-16) ayant chacun un RSIB 30 (non asymptotique), noté RSIBn, n = 1, ..., N. La technique EESM consiste à trouver un unique RSIB, appelé RSIB compressé, qui permette d'estimer le taux d'erreur binaire en sortie du décodeur de canal. Pour cela, on construit une table de conversion qui donne le TEB en sortie du décodeur de canal en fonction du RSIB entrant pour un canal de propagation gaussien. La construction de cette table est réalisée soit par simulation, soit analytiquement lorsque la théorie le permet. Le RSIB compressé est donné par la formule suivante : R 1 "états p/ RSIB (7) RSIBcompressé - - In N ex états 1=1 (~ \ r // où [3 est un paramètre qui dépend du schéma de modulation et de codage et "états représente le nombre d'états du canal sur le bloc d'information codée. Il est à noter que le paramètre 13 peut être obtenu par simulation informatique de la même façon que le paramètre A qui figure dans l'équation (13) donnée plus loin. On se reportera donc à la description détaillée de la méthode d'obtention de A. Conformément à la présente invention, on applique cette technique de compression du RSIB au RSIB asymptotique obtenu dans le cas d'un système OFDM-CDMA par la formule (3) donnée en introduction, de façon à obtenir une valeur de RSIB asymptotique par bloc d'information codée, qui sera appelée RSIB asymptotique compressé. On suppose que le domaine des fréquences est découpé en F = N/SF sous-bandes de SF sous-porteuses, N étant la taille de la transformation de Fourier rapide utilisée par le récepteur et SF étant le paramètre d'étalement fréquentiel, et le domaine temporel est découpé en T = Ns/ST blocs de ST symboles OFDM, NS étant le nombre de symboles OFDM utilisés pour coder une trame de données et ST étant le paramètre d'étalement temporel. Comme le montre l'organigramme de la figure 8, dans un mode particulier de réalisation, une première étape 80 du procédé de sélection des paramètres d'étalement consiste à estimer les coefficients hp[n] représentant l'atténuation du canal de propagation pour la nième sous-porteuse du p1eme symbole OFDM, pour p= 1, ..., Nsetn= 1, ..., N. Puis ces estimées du canal permettent, lors d'une étape 82, de déterminer pour l'ensemble des F sous-bandes de SF sous-porteuses et des T blocs de ST symboles OFDM au moins un RSIB asymptotique selon la formule (3). On suppose que les paramètres suivants sont connus préalablement à la détermination du RSIB asymptotique : • p : puissance moyenne des codes d'étalement interférents, • po : puissance du signal étalé par le code d'étalement dont on souhaite calculer le RSIB, • 62 : variance du bruit gaussien additif (contenant également la puissance des cellules interférentes dans le cas d'un environnement multi-cellulaire), • hp[n] : atténuation du canal pour la nième sous-porteuse du pième 10 symbole OFDM ; • L : facteur d'étalement, c'est-à-dire longueur des codes d'étalement, • K : nombre de codes d'étalement transmis, • N : taille de la transformation de Fourier rapide, • Ns : taille d'un "slot", c'est-à-dire nombre de symboles OFDM utilisés 15 pour coder une trame de données. Certains des paramètres ci-dessus sont connus par l'intermédiaire de la signalisation et d'autres peuvent être obtenus avec un mécanisme d'estimation de canal. Un slot contenant NS symboles OFDM de N sous-porteuses, il y a 20 NNs/L RSIB asymptotiques différents, notés y(i,j) , à calculer selon l'équation (3). Comme le montre la figure 9, le domaine des fréquences peut être découpé en F = N/SF sous-bandes de SF sous-porteuses, tandis que le domaine temporel est découpé en T = Ns/ST blocs de ST symboles OFDM. Le RSIB asymptotique y(i, j) est le RSIB de la i1ème sous-bande du j'ème bloc, i = 0, ..., F-1, j = 0, ..., T-1. 25 Dans le cas où le procédé de sélection des paramètres d'étalement conforme à la présente invention est mis en oeuvre dans un système d'émission-réception où le récepteur comporte un égaliseur MMSE mono-utilisateur, le RSIB asymptotique y(i, j) est défini par la formule suivante : i 21 Po ' 112(i,J) 4 2 + 2 -r12(i,J) 1 SF-1ST-1 i(i,J) _ ~ ~ p=0 k=0 OU = 1 K-1 a-lpi K ;=o K-1 P K-p' où : h,sT+k[iSF +pl 2 h,ST+kUSF +pl + h,sT+k [i SF + pl h;sT+k [i SF + pl 2 SF 1ST 1 hi.sT+k [i SF + pl /1 1z L p-0 k =O hJ.ST+k[i . SF + p1~2 + I L SE F-ISIT-1 p=0 k=0 ap 62 2 2 • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = KIL est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents et les p; sont les puissances des codes d'étalement interférents, avec i ~ 0, • p et k sont des entiers positifs, • hi.sT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i SF+p)1ème sous-porteuse du (j.ST+k)ième symbole OFDM, et • a2 est la variance du bruit. Toujours pour un égaliseur MMSE mono-utilisateur, dans le cas K-1 1 K-1 particulier où K -1 1p' = K , où les p;, pour i # 0, sont les puissances des i=o codes d'étalement interférents, le RSIB asymptotique est défini par la relation suivante : 22 Y(i,j) = Po Tl(i, aP(1ù J)) SF-1ST-1 11(i, j) = ~ 1 1 p=0 k=0 hj.sT+k [i SF + P] 2 hj S T [i • SF + p] + l,l 2 (8) ou ap K-1 P KùP' avec les mêmes notations que précédemment. Pour un égaliseur MRC, le RSIB asymptotique Y(i, j) est défini par la formule suivante : Po • Tl2(i,j) 4 h;.sT+k[i • SF + P] ù rl2(i,j) + a2Tl(i,j) ) Y(i, J) = (9) ap 1 SF-1ST-1 où rl(i, j) = ù E hi ST+k [i • SF + p] , avec les mêmes notations que I- p=O k=0 précédemment. Pour un égaliseur MMSE multi-utilisateur, le RSIB asymptotique y(i, j) est solution de l'équation implicite suivante : 1 SF-1ST-1 1 = 1 11 2 (10) p=0 k=0 h;sT+k [i • SF + p] avec les mêmes notations que précédemment et : A(a,Pk= (1ù a) 70' i) + a • Y(i,j) • m(Pk,Y(i,j)) Pk B(a,Pkj)) =1ù Pk • m(Pk,Y(i,j)) (11) m(Pk, Y(i, J)) _ Power (u) du ~Pk +Y(i,j) u où power (p) est la distribution limite de la puissance lorsque la taille L des codes d'étalement et le nombre K de codes d'étalement transmis tendent vers l'infini et 15 a = K / L reste constant. Par exemple, si le système comporte K0 classes de 2 2 A(a,Pk+ a • B(a,Pk,Y(i,j)) Kc-1 puissance pi, ..., PKc, alors ['power (u) = I yi8(p ù pi), où y; est le ratio de codes i=o appartenant à la classe i de puissance pi, et 8(.) est la distribution de Dirac. Pour un égaliseur EGC, le RSIB asymptotique est défini par la relation suivante : y(I,1) = 2 h;.sT+k[i SF +pl ù g2(1,1) i Po '~12(i,l) ap (12) SF-1ST-1 avec r1(i, = L ~ p=0 k=0hjST+k Ii .SF + pl et les mêmes notations que précédemment. Après l'étape 82 de détermination des RSIB asymptotiques y(i, j) , l'étape 84 consiste à en déduire le RSIB asymptotique compressé, au moyen de la formule suivante : F-1 T-1 Y(i l) RSIBcompressé = ùXi In - L lle NNs i=o i=o où le scalaire a, ne dépend que de la modulation et du schéma de codage (rendement du codeur et entrelaceur bit). En effet, ? est indépendant du canal de propagation, de la charge du système et de la combinaison (SF,ST). Les différentes valeurs de ?Li sont connues à l'avance par le récepteur. Le scalaire peut être obtenu par simulation informatique. Les simulations requises sont réalisées suivant la procédure suivante. Tout d'abord, pour une trame i d'information codée, on tire un ensemble de réalisations de canal correspondant au modèle de canal considéré. Cela permet, pour la trame i d'information codée, d'obtenir un RSIB asymptotique compressé qu'on notera RSIBicompressé(2), dépendant de ?. Ensuite, à l'aide d'une table (obtenue également par simulation) donnant le Taux d'Erreur Binaire (TEB) à la sortie du décodeur de canal en fonction du RSIB à son entrée pour un canal gaussien, on obtient un Taux d'Erreur Binaire théorique dépendant de ?, qu'on notera TEB'théorique(%). Ce TEB théorique dépendant de a, est ensuite comparé au TEB réel, noté TEBiréel, mesuré par simulations à la sortie du décodeur de canal. (13) Afin d'obtenir plusieurs valeurs de TEB'théorigue(2) pour une même valeur de 2, ces opérations sont répétées sur un grand nombre de trames (par exemple 1000 trames). Chaque TEB'théorique(?) est ensuite comparé au TEB'réei. Le scalaire a, est alors obtenu par la minimisation, sur l'ensemble des trames d'information codées, de la différence entre TEB'rée, et TEB'théor;gue(2). Le scalaire a, dépendant du schéma de modulation et de codage, le processus est répété pour l'ensemble des schémas de modulation et de codage considérés. L'unité d'estimation des conditions de propagation du canal 70 (voir 10 figure 6) peut donc calculer le RSIB asymptotique compressé pour différents couples (SF,ST) et déterminer la combinaison optimale. Ainsi, à l'étape 86, on sélectionne les valeurs des paramètres d'étalement fréquentiel SF et temporel ST qui maximisent le RSIB asymptotique compressé calculé précédemment. 15 Ce processus d'optimisation du RSIB compressé peut être réalisé de diverses façons. Pour sélectionner, à l'étape 86, le meilleur couple (SF,ST), on peut procéder à une recherche de type exhaustive, en calculant le RSIB asymptotique compressé pour différentes configurations de SF et ST. Par exemple, si le facteur 20 d'étalement L = SFXST vaut 32, on peut calculer le RSIB asymptotique compressé pour (SF,ST) E {(1,32), (4,8), (8,4), (32,1)}. La valeur la plus élevée parmi les quatre valeurs du RSIB ainsi obtenues fournira ainsi la configuration (SF,ST) à retenir, pour un couple modulation/codage donné. En variante, on peut appliquer une méthode de recherche par la 25 section Golden. On se reportera utilement à ce sujet à l'ouvrage de W. H. PRESS, B. P. FLANNERY, S. A. TEUKOLSKY et W. T. VETTERLING publié en 1988 aux éditions Cambridge Press University, intitulé "Numerical Recipes in C" et plus particulièrement à son chapitre 10.1, pages 397 à 402. Il s'agit d'une méthode permettant de réduire l'intervalle de recherche de façon optimale lorsqu'on 30 considère un système à une inconnue, soit le paramètre SF, soit le paramètre ST. Cette méthode sera avantageusement utilisée dans le cas d'un utilisateur qui se déplace lentement, par exemple à une vitesse de l'ordre de 3 km/h, ou dans le cas où le canal varie peu en fréquence. Dans le premier cas, où l'utilisateur se déplace lentement, on fixe la valeur du paramètre ST (par exemple ST = 16) et on recherche SF par la méthode de Golden. Dans le cas où le canal varie peu en fréquence (par exemple lorsque la bande de cohérence est supérieure à SF), on fixe la valeur de SF (par exemple SF = 16) et on recherche ST par la méthode de Golden. Selon une autre variante, pour sélectionner le meilleur couple (SF,ST), on peut appliquer une méthode d'optimisation sous contrainte à l'aide de la technique du multiplicateur de Lagrange, avec pour contrainte g(SF,ST) = SFST - L. La technique du multiplicateur de Lagrange est bien connue de l'homme du métier ; on pourra en trouver une description sur Internet à l'adresse suivante : http://ocw.mit.edu/NR/rdonlyres/Mathematics/l 8-466Mathematical-StatisticsSpring2003/F6138C74-AA64-4B63-9A6038A4833AECC B/0/appf.pdf La présente invention propose également une simplification du calcul du RSIB asymptotique compressé donné par la formule (13), évitant de calculer les logarithmes et exponentielles des quantités mises en jeu. Pour cela, on utilise l'approximation suivante : In(ea + eb) , max(a,b) + fc (la ù b,) (14) où 0a ù bl) est une table pré-calculée ne dépendant que de (la ù b) et approchant de façon très précise la fonction ln(1 + e-la- bl ). On se reportera utilement à ce sujet à l'article de P. ROBERTSON, P. HOEHER et E. VILLEBRUN intitulé "Optimal and sub-optimal maximum a posteriori algorithms suitable for turbo-decoding", European Transactions on Telecommunications, pages 119 à 125, 1997. Le calcul de l'approximation du RSIB asymptotique compressé se fait donc de façon itérative : en supposant qu'on en connaisse n-1 valeurs, on en déduit la n1ème de la façon suivante : In(ea' +•••+ea )=ln(A+ea ) avec A=ea' +...+ea0 =ea In(ea' + • • + en ) = max(In(A), an) + fc In(A) ù an I) = max(a, ann) + fc (a ù an) (15) L'invention peut être appliquée à un dispositif récepteur tel qu'un téléphone mobile. Comme le montre la figure 10, on considère un système comportant un réseau 200 de diffusion de téléphonie mobile composé d'un serveur de réseau, d'une infrastructure de transmission et d'un parc d'appareils de réception sans fil mobiles ou portables 202, 204, par exemple des téléphones mobiles associés au réseau. Des messages 206, 208 sont envoyés aux appareils portables 202, 204 tandis que ceux-ci réémettent des informations 210, 212 en retour, par exemple l'information d'atténuation du canal ou les paramètres d'étalement aptes à optimiser la transmission. Les appareils portables peuvent donc réaliser le calcul du RSIB asymptotique compressé suivant la formule (13) afin de sélectionner les paramètres d'étalement aptes à maximiser ce RSIB asymptotique compressé. Ces paramètres d'étalement sont ensuite renvoyés à une station de base émettrice 214, qui va pouvoir ajuster divers paramètres d'émission, comme par exemple les paramètres d'étalement ou le schéma de modulation-codage. L'optimisation des paramètres d'étalement peut être faite directement au niveau de l'émetteur si le dispositif mobile renvoie à l'émetteur l'information d'atténuation du canal et le type d'égaliseur utilisé. La présente invention peut également être appliquée aux simulateurs de systèmes radio-mobiles. Un tel simulateur système est un logiciel permettant de modéliser un environnement multi-cellulaire afin d'obtenir les performances d'un réseau en termes de débit à l'intérieur d'une cellule, de distribution des débits moyens des utilisateurs à l'intérieur d'une cellule ou de distribution des délais moyens dans la transmission de paquets. Le simulateur système est implanté dans une unité du type de celle représentée sur la figure 11 comportant au moins un processeur 100, un ensemble de mémoires RAM 102, pour le stockage des données, et de mémoires ROM 104, par exemple pour le stockage d'instructions du programme. Ces divers éléments sont reliés par un bus 106. Un élément périphérique tel qu'un clavier 108 permet à l'utilisateur du simulateur système d'entrer des données en réponse à un message affiché sur un écran de visualisation 110. L'utilisateur peut ainsi sélectionner par exemple le modèle de mobilité (position et vitesse du mobile) ou un modèle pour le canal de propagation. L'intégration du calcul du RSIB asymptotique compressé suivant la formule (13) dans un simulateur système permet non seulement la sélection des paramètres d'étalement optimaux pour un système OFDM-CDMA, mais aussi l'obtention du taux d'erreur paquet après décodage de canal, ce qui est un indicateur de la qualité de la communication. REVENDICATIONS 1. Procédé de sélection des paramètres d'étalement fréquentiel (SF) et/ou temporel (ST) pour un système de communication incluant un émetteur et un récepteur et utilisant une modulation d'amplitude en quadrature et mettant en oeuvre une technique d'accès multiple à répartition par codes ou CDMA et un multiplexage par division de fréquences orthogonales ou OFDM, dans lequel le domaine des fréquences est découpé en F = N/SF sous-bandes de SF sous-porteuses, N étant la taille de la transformation de Fourier rapide utilisée par le récepteur et SF étant le paramètre d'étalement fréquentiel, et le domaine temporel est découpé en T = Ns/ST blocs de ST symboles OFDM, Ns étant le nombre de symboles OFDM utilisés pour coder une trame de données et ST étant le paramètre d'étalement temporel, ledit procédé étant caractérisé en ce qu'il comporte des étapes suivant lesquelles : - on estime (80) les coefficients de l'atténuation du canal de propagation ; - on détermine (82), pour l'ensemble des F sous-bandes de SF sous-porteuses et des T blocs de ST symboles OFDM, au moins un Rapport Signal à Interférence + Bruit (RSIB) asymptotique, à partir des coefficients de l'atténuation du canal de propagation et du type d'égaliseur utilisé par le récepteur ; - on en déduit (84), pour chaque bloc d'information codée obtenu à l'issue de la modulation, le Rapport Signal à Interférence + Bruit asymptotique compressé, qui s'exprime par une formule analytique dépendant desdits paramètres d'étalement fréquentiel (SF) et temporel (ST) ; et - on sélectionne (86) lesdits paramètres d'étalement fréquentiel (SF) et temporel (ST) qui maximisent la valeur du Rapport Signal A Interférence + Bruit asymptotique compressé. 2. Procédé selon la revendication 1, caractérisé en ce qu'à l'étape (84) de calcul du RSIB asymptotique compressé, on calcule le RSIB asymptotique compressé au moyen de la formule suivante :( L F-1 T-1 _Y(i.i) NN Rslgcompressé = ZE ù%In e S i-o )-o OU : • 2 est un scalaire connu à l'avance par le récepteur, dépendant de la modulation et du codage, • L est le facteur d'étalement : L = SFXST, • i et j sont des entiers positifs, et • y(i, j) est le RSIB asymptotique de la i'ème sous-bande du jième bloc de ST symboles OFDM. 3. Procédé selon la revendication 2, caractérisé en ce qu'à l'étape (84) de calcul du RSIB asymptotique compressé, on calcule le RSIB asymptotique compressé de façon itérative en utilisant l'approximation suivante : In(ea +eb) max(a,b)+fcùb ) où 0a --b) est une table pré-calculée ne dépendant que de (la ùb) et 15 approchant la fonction In(1 + e-la-bl ) 4. Procédé selon la revendication 1, 2 ou 3, caractérisé en ce que l'étape de sélection (86) consiste à procéder à une recherche exhaustive du couple de paramètres d'étalement fréquentiel et temporel (SF,ST) qui maximise la valeur du RSIB asymptotique compressé. 20 5. Procédé selon la revendication 1, 2 ou 3, caractérisé en ce que l'étape de sélection (86) consiste à appliquer une méthode de recherche par la section Golden. 6. Procédé selon la revendication 1, 2 ou 3, caractérisé en ce que l'étape de sélection (86) consiste à appliquer une méthode d'optimisation sous 25 contrainte utilisant la technique du multiplicateur de Lagrange, où la contrainte s'écrit g(SF,ST) = SFST - L, SF étant le paramètre d'étalement fréquentiel, ST étant le paramètre d'étalement temporel et L étant le facteur d'étalement. 7. Procédé selon l'une quelconque des revendications précédentes, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE 30 (minimisation de l'erreur quadratique moyenne, en anglais "Minimum MeanSquare Error") mono-utilisateur, caractérisé en ce que le RSIB asymptotique y(i,j) de la iième sous-bande du jlème bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante : po'~12(i,1) Y(i, _ 4 2 ap )2 i2(i,i) 1 sF-1sT-1 IhjsT+k [i . SF + ~ ~~( 2 hj.sT+k [i SF + pl + 62 S1S~1 hj.ST+k L. SF + pi Lr 2 / + L p-0 k 0 hjST+k Ij SF + p~12 + ) hjST+k [j ' SF + p] 2 h;sT+k [i ' S F+ p] + 2 T1(i, i) _ ~ ~ p=0 k=0 1 SF-1ST-1 où : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes 10 d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents et les pi sont les puissances des codes d'étalement interférents, avec i ~ 0, • p et k sont des entiers positifs, • hi sT+k [i . SF + p] est l'atténuation du canal de propagation pour la 15 (i.SF+p)ième sous-porteuse du (j.ST+k)jème symbole OFDM, et • 62 est la variance du bruit. 8. Procédé selon la revendication précédente, dans lequel on a 1 K-1 1 K-1 1p; = -ER , où K est le nombre de codes d'étalement transmis et les pi, K -1 ;=1 K ;=0 pour i # 0, sont les puissances des codes d'étalement interférents, caractérisé en 20 ce que le RSIB asymptotique est défini par la relation suivante : 2 31y(i,1) = po ' î1(i, j) ap(1û ri(i,j)) 1 SF-1sT-1 ~(i, j) _ L ~ p=0 k=0 2 6 = CCP K1 P Kû1~p' jST+kli SF +p] 2 jST+k[j SF +p] + Ou 2 où : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K / L est la charge du système, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • hj.sT+k [i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)ième sous-porteuse du (j.ST+k)ième symbole OFDM, et • 62 est la variance du bruit. 9. Procédé selon l'une quelconque des revendications 1 à 6, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MRC (combinaison à rapport maximum, en anglais "Maximum Ratio Combining"), caractérisé en ce que le RSIB asymptotique y(i,j) de la iième sous-bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 etj = 0, ..., T-1, est défini par la formule suivante : _ po . 12(i,j) 7(',J) ( 1 SF-1ST 1 4 2 a p hjST+k Ii . SF + p] û (i, j) + 6` 1(i, l) P=0 k=0 / 1 SF-1ST-1 avec r)0,1) _ -, L p=0 k=0 hjsT+k[i . SF + p] 2 où: • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB,• L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • hi.sT+k [i SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)'eme sous-porteuse du (j.ST+k)ieme symbole OFDM, et • 6-2 est la variance du bruit. 10. Procédé selon l'une quelconque des revendications 1 à 6, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Error") multi- utilisateur, caractérisé en ce que le RSIB asymptotique y(i, j) de la Mme sous- bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est solution de l'équation implicite suivante : 1 SF-1ST-1 1 1=ù~~ 2 L P=o k_p a h; sT+k [i SF + P] où : • L est le facteur d'étalement : L = SFXST, • a = K / L , K étant le nombre de codes d'étalement transmis, • p et k sont des entiers positifs, 20 • hisT+k[i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)'ème sous-porteuse du (j.ST+k)ième symbole OFDM, et • 62 est la variance du bruit, • pk est la puissance appliquée au k'ème utilisateur, et : A(a,Pk,Y(i,j)) = (1û a) Y(i,J) + a • Y(i,j) m(Pk,Y(i,j)) Pk B(a,Pk,Y(i,j)) = 1- Pk m(Pk,Y(i,j)) m(Pk, YO, J)) _ Power (u) du 'Pk + Y(i,j) • u 2 A(a,Pk, Y(i,j))+a B(a,Pk,Y(i,j))où power (p) est la distribution limite de la puissance lorsque L et K tendent vers l'infini et a = K / L reste constant. 11. Procédé selon l'une quelconque des revendications 1 à 6, mis en oeuvre en liaison avec un récepteur comportant un égaliseur EGC (combinaison à gain identique, en anglais "Equal Gain Combining"), caractérisé en ce que le RSIB asymptotique y(i,j) de la hème sous-bande du jieme bloc de ST symboles OFDM, pour i = 0, ..., F-1 etj = 0, ..., T-1, est défini par la formule suivante : Y(i,l) = Po . T12(i,J) ( 1 SF-1ST_1 2 ap ~1 1 h;sT+k[i SF+p] -T12(i,j) p= k=0 +a` sr-1ST-1 avec j) = ~ 1 1 p=0 k=0h;.sT+k[i.SF +p] où : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes 15 d'étalement transmis, • est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • hi .sT+k [i . SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)ième sous-porteuse du (j.ST+k)ième symbole OFDM, et 20 • 62 est la variance du bruit. 12. Dispositif de sélection des paramètres d'étalement fréquentiel (SF) et/ou temporel (ST) pour un système de communication incluant un émetteur et un récepteur et utilisant une modulation d'amplitude en quadrature et mettant en oeuvre une technique d'accès multiple à répartition par codes ou 25 CDMA et un multiplexage par division de fréquences orthogonales ou OFDM, dans lequel le domaine des fréquences est découpé en F = N/SF sous-bandes de SF sous-porteuses, N étant la taille de la transformation de Fourier rapide utilisée par le récepteur et SF étant le paramètre d'étalement fréquentiel, et ledomaine temporel est découpé en T = Ns/ST blocs de ST symboles OFDM, Ns étant le nombre de symboles OFDM utilisés pour coder une trame de données et ST étant le paramètre d'étalement temporel, ledit dispositif étant caractérisé en ce qu'il comporte : - des moyens pour estimer les coefficients de l'atténuation du canal de propagation ; - des moyens pour déterminer, pour l'ensemble des F sous-bandes de SF sous-porteuses et des T blocs de ST symboles OFDM, au moins un Rapport Signal à Interférence + Bruit (RSIB) asymptotique, à partir des coefficients de l'atténuation du canal de propagation et du type d'égaliseur utilisé par le récepteur ; - des moyens (70) de calcul pour en déduire, pour chaque bloc d'information codée obtenu à l'issue de la modulation, le Rapport Signal à Interférence + Bruit asymptotique compressé, qui s'exprime par une formule analytique dépendant desdits paramètres d'étalement fréquentiel (SF) et temporel (ST) ; et - des moyens pour sélectionner lesdits paramètres d'étalement fréquentiel (SF) et temporel (ST) qui maximisent la valeur du Rapport Signal à Interférence + Bruit asymptotique compressé. 13. Dispositif selon la revendication 12, caractérisé en ce que les moyens (70) de calcul du RSIB asymptotique compressé sont adaptés à calculer le RSIB asymptotique compressé au moyen de la formule suivante : L F-1 T-1 Y(iJ) RSIBcompressé = -k ln 1le NN S;=o J=O où : • a, est un scalaire connu à l'avance par le récepteur, dépendant de la modulation et du codage, • L est le facteur d'étalement : L = SFXST, • i et j sont des entiers positifs, et • y(i, j) est le RSIB asymptotique de la i1ème sous-bande du jième bloc de ST symboles OFDM. 14. Dispositif selon la revendication 13, caractérisé en ce que les moyens (70) de calcul du RSIB asymptotique compressé sont adaptés à calculer le RSIB asymptotique compressé de façon itérative en utilisant l'approximation suivante : In(ea + eb) max(a,b) + fc (a û b ) où fc(a ûb) est une table pré-calculée ne dépendant que de (la û b) et approchant la fonction In(1 + e* )15. Dispositif selon la revendication 12, 13 ou 14, caractérisé en ce que les moyens de sélection sont adaptés à procéder à une recherche exhaustive du couple de paramètres d'étalement fréquentiel et temporel (SF,ST) qui maximise la valeur du RSIB asymptotique compressé. 16. Dispositif selon la revendication 12, 13 ou 14, caractérisé en ce que les moyens de sélection sont adaptés à appliquer une méthode de recherche par la section Golden. 17. Dispositif selon la revendication 12, 13 ou 14, caractérisé en ce que les moyens de sélection sont adaptés à appliquer une méthode d'optimisation sous contrainte utilisant la technique du multiplicateur de Lagrange, où la contrainte s'écrit g(SF,ST) = SFST - L, SF étant le paramètre d'étalement fréquentiel, ST étant le paramètre d'étalement temporel et L étant le facteur d'étalement. 18. Dispositif selon l'une quelconque des revendications 12 à 17, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Error') mono-utilisateur, caractérisé en ce que le RSIB asymptotique y(i, j) de la i'ème sous-bande du jieme bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante :( 36 Po ' ~12(i,1) 2 h;.ST+k U SF + pll z 12 hjsT+k[i'SF +p]i + ) 2 SF-1ST-1 6 + L p=0 k=0 7 4 1 SF-1ST-1 L~~7 1 SF-1ST-1 J) _ ~ Il p=0 k=0h;.sT+k [i SF + p] 2 lII2 ST+k [j . SF + p]l + / ap h h;ST+k[G-SF +p] 2 h;.sT+k[i • SF + p) + Il 2 -T12(i,j) où 62 - 1 K-1 aùlpi K ;=o K-1 p K-11,p' où : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents et les p; sont les puissances des codes d'étalement interférents, avec i ~ 0, • p et k sont des entiers positifs, • hj.ST+k [i • SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)Ième sous-porteuse du (j.ST+k)ième symbole OFDM, et • 62 est la variance du bruit. 19. Dispositif selon la revendication précédente, dans lequel on a K-1 K-1 = -Ip; , où K est le nombre de codes d'étalement transmis et les pi, K -1 i=1 K ;=o pour i # 0, sont les puissances des codes d'étalement interférents, caractérisé en ce que le RSIB asymptotique est défini par la relation suivante : 37Y(i,J) = po ' fl(i,j) apO ù 71(i3O 2 h;sT+k[i . Sr + p] hj.ST+k [i ' SF + p] 2+ 1 SF-1ST-1 11(i, j) = ù, 1 1 p=0 k=0 où p= où : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K I L est la charge du système, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • h; sT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)lème sous-porteuse du (j.ST+k)ième symbole OFDM, et • 62 est la variance du bruit. 20. Dispositif selon l'une quelconque des revendications 12 à 17, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MRC (combinaison à rapport maximum, en anglais "Maximum Ratio Combining"), caractérisé en ce que le RSIB asymptotique y(i,j) de la iième sous-bande du j'ème bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, T-1, est défini par la formule suivante : Po '112(i,j) 4 h;.sT+k[i'SF +p] -2(i,j) +62r1(i,j) / 70, J) ( 1 SF-1ST-1 a,p -, p=0 k=0 1 SF-1ST-1 avec n(i, j) = ù 1 1 L p=0 k=0 h;sT+k [i ' Sr + p] 2 où :• po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K/L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, • h. sT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i SF+p)ième sous-porteuse du (j.ST+k)ieme symbole OFDM, et • 6z est la variance du bruit. 21. Dispositif selon l'une quelconque des revendications 12 à 17, mis en oeuvre en liaison avec un récepteur comportant un égaliseur MMSE (minimisation de l'erreur quadratique moyenne, en anglais "Minimum Mean Square Erro(') multi-utilisateur, caractérisé en ce que le RSIB asymptotique y(i, j) de la i1ème sous-bande du jieme bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est solution de l'équation implicite suivante : 1 $F-13T-1 1=~ ~ ~ 6z p=0 k=0 2 1 A(a,pk,y(i,j)) + ' B(a,pk,y(i,j)) hj sT+k [) .SF +p] OU : • L est le facteur d'étalement : L = SFXST, • a = K / L , K étant le nombre de codes d'étalement transmis, • p et k sont des entiers positifs, • h;.sT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i.SF+p)'ème sous-porteuse du (j.ST+k)ième symbole OFDM, et • az est la variance du bruit, 25 • pk est la puissance appliquée au kieme utilisateur, et :A(a,pk,Y(i,j)) = (1 ù a) Y0,1) + a. Y(i,j) .m(pk,Y(i,j)) Pk B(a,pk,Y(i,j)) = 1 ù Pk m(pk,Y(i,j)) m(pk,Y(i,j)) _ f- power(u) du pk + Y(i, j) u où power (p) est la distribution limite de la puissance lorsque L et K tendent vers l'infini et a = K / L reste constant. 22. Dispositif selon l'une quelconque des revendications 12 à 17, mis en oeuvre en liaison avec un récepteur comportant un égaliseur EGC (combinaison à gain identique, en anglais "Equal Gain Combining"), caractérisé en ce que le RSIB asymptotique y(i,j) de la i1ème sous-bande du jième bloc de ST symboles OFDM, pour i = 0, ..., F-1 et j = 0, ..., T-1, est défini par la formule suivante : Y(i, j) _ Po .1120, j) ~ 2 h;sT+k[i ~ SF + p] ù T12(i,j) ) ap 1 SF-1ST-1 avec rt(i, j) _ L -1 1 h;.sT+k [i SF + p=0 k=0 OU : • po est la puissance du signal étalé par le code d'étalement dont on souhaite déterminer le RSIB, • L est le facteur d'étalement : L = SFXST, • a = K /L est la charge du système, K étant le nombre de codes d'étalement transmis, • p est la puissance moyenne des codes d'étalement interférents, • p et k sont des entiers positifs, 20 • hj.sT+k [i • SF + p] est l'atténuation du canal de propagation pour la (i SF+p)1ème sous-porteuse du (j.ST+k)ième symbole OFDM, et • 6-2 est la variance du bruit.23. Unité de simulation d'un système radio-mobile, caractérisée en ce qu'elle est adaptée à mettre en oeuvre un procédé de sélection des paramètres d'étalement selon l'une quelconque des revendications 1 à 11. 24. Unité de simulation d'un système radio-mobile, caractérisée en ce qu'elle comporte un dispositif de sélection des paramètres d'étalement selon l'une quelconque des revendications 12 à 22. 25. Dispositif de téléphonie mobile sans fil, caractérisé en ce qu'il est adapté à mettre en oeuvre un procédé de sélection des paramètres d'étalement selon l'une quelconque des revendications 1 à 11. 26. Dispositif de téléphonie mobile sans fil, caractérisé en ce qu'il comporte un dispositif de sélection des paramètres d'étalement selon l'une quelconque des revendications 12 à 22.