FR 2873876 A1 20060203 FR 0408296 A 20040727 2873876 1 PROCEDE D'ETALEMENT DE SPECTRE AVEC EFFACEMENT DE SOUS-BANDES. La présente invention concerne un procédé et un dispositif d'étalement de spectre. Elle s'applique notamment, mais non exclusivement aux transmissions par satellite et aux télécommunications mobiles, utilisant des techniques d'étalement 10 de spectre. Lors de transmissions brouillées par des interféreurs à bande de fréquences limitée, les récepteurs mettent en général en place des stratégies de détection et d'élimination des brouilleurs. Il est également possible d'utiliser une technique d'étalement spectral qui limite les influences du brouilleur. Néanmoins, il n'existe aucune méthode simple pour réaliser un étalement spectral excluant certaines bandes de fréquences prédéfinies. Par ailleurs, on a déjà proposé d'utiliser des filtres linéaires périostiques numériques LPTV (Linear Periodic Time Varying) pour généraliser les techniques de transmission DS-CDMA (Direct Sequence Code Division Multiple Access) classiques, dans lesquelles la séquence pseudo-aléatoire est remplacée par un filtre de ce type utilisé par exemple comme un entrelaceur. La présente invention a pour but de proposer un procédé d'étalement de spectre excluant certaines bandes de fréquences prédéfinies. Cet objectif est atteint par la prévision d'un procédé d'étalement de spectre d'un signal d'entrée numérique. Selon l'invention, ce procédé comprend des étapes consistant à : choisir des coefficients d'un filtre linéaire périodique numérique réversible, de manière à ce qu'il étale le spectre du signal d'entrée pour occuper une bande de fréquences étalée à l'exception d'au moins une sous-bande de fréquences effacée, incluse dans la bande de fréquences étalée, et traiter le signal d'entrée à l'aide du filtre linéaire périodique numérique de manière à obtenir un signal de sortie à spectre étalé. De préférence, la sortie du filtre linéaire périodique est apte à conserver l'enveloppe du signal d'entrée. Selon un mode de réalisation préféré de l'invention, le filtre linéaire périodique numérique présente une transformée en z dans une représentation fréquentielle 5 de la forme suivante: N-1 Y(z) = Tp É X(zWp) p=0 dans laquelle: N-1 / 21TC \ T = 1 E Wg(k)+pk w = exp P Nk=O N i N est la période du filtre, et g(k) pour les nombres entiers k compris entre 0 et N 1 sont des coefficients à 10 déterminer. De préférence, le choix des coefficients g(k) comprend des étapes consistant à annuler les termes Tp du filtre pour toutes les valeurs de p telles que: E/fmax bmin N p Er 'min borax N\ -1 B N - B 1 N expression dans laquelle [fmin, fmax] représente la bande de fréquences du signal numérique d'entrée, B représente la bande de fréquence étalée, [bmin, bmax] représente la sous-bande de fréquences effacée, incluse dans la bande de fréquences B, E(a) représente la fonction partie entière et [a]N représente la fonction "modulo N". Avantageusement, l'annulation d'un terme Tp du filtre linéaire périodique consiste à regrouper tous les entiers compris entre 0 et N-1 en couples (k:0, kl) tels que: Wg(kO)+pkO + Wg(kl)+pkl = 0. Selon un mode de réalisation préféré de l'invention, le filtre linéaire périodique est applicable à un signal d'entrée à enveloppe constante. Selon un autre mode de réalisation préféré de l'invention, le signal d'entrée est 30 préalablement suréchantillonné. De préférence, avant d'être suréchantillonné, le signal d'entrée est traité par un modulateur numérique présentant une constellation à m points. 2873876 -3- Selon un mode de réalisation préféré de l'invention, le signal de sortie du filtre linéaire périodique est suréchantillonné. Selon un mode de réalisation préféré de l'invention, le signal de sortie 5 suréchantillonné est traité par un filtre passe-bas en racine de cosinus surélevé. Selon un mode de réalisation préféré de l'invention, le filtre linéaire périodique présente une période de la forme 2n, n étant un nombre entier positif. Un mode de réalisation préféré de l'invention sera décrit ci-après, à titre d'exemple non limitatif, avec référence aux dessins annexés dans lesquels: La figure 1 représente un exemple de chaîne de transmission numérique mettant en oeuvre le procédé d'étalement de spectre 15 selon l'invention; La figure 2 montre plus en détail le filtre linéaire périodique figurant dans la chaîne de transmission représentée sur la figure 1; La figure 3 est une représentation bifréquentielle d'un exemple de 20 filtre du type représenté sur la figure 2; La figure 4 montre plus en détail un mode de réalisation du filtre représenté sur la figure 2; Les figures 5 à 13 montrent des exemples de spectres susceptibles d'être obtenus grâce au filtre représenté sur la figure 4. La figure 1 représente une chaîne de transmission de signaux numériques e(n) par étalement de spectre. Cette chaîne comprend tout d'abord un modulateur numérique 2, qui effectue un étiquettage des signaux d'entrée e(n) se présentant sous la forme de trains de bits. Le signal de sortie du modulateur numérique 2 est appliqué à l'entrée d'un premier suréchantillonneur 3 qui suréchantillonne le signal e(n) modulé en répétant L fois chaque échantillon de ce signal. Le signal suréchantillonné x(n) en sortie du suréchantillonneur 3 est appliqué à l'entrée d'un filtre linéaire périodique numérique 1 selon l'invention. Les échantillons y(n) obtenus en sortie du filtre 1 sont ensuite traités par un second 2873876 -4- suréchantillonneur 4 qui multiplie par un facteur K le nombre d'échantillons du signal qu'il traite, par exemple en insérant K échantillons nuls entre chaque échantillon du signal appliqué en entrée. Le signal de sortie du second suréchantillonneur 4 est enfin traité par un filtre numérique passe-bas 5 assurant une mise en forme spectrale du signal et délivrant un signal s(n) de sortie de la chaîne de transmission. Le modulateur numérique est par exemple un modulateur à déplacement de phase (MDP-M, M-PSK, M-MSK) et/ou d'amplitude (MAQ-M), à m points de 10 constellation. Selon l'invention, le filtre 1 est choisi de manière à obtenir un signal de sortie dont le spectre est étalé avec des sous-bandes effacées. Ce filtre est un filtre linéaire périodique variant en fonction du temps (LPTV) illustré schématiquement sur la figure 2 sous la forme d'une représentation commutée en entrée et en sortie (MIMO: Multiple Input Multiple Output). Le filtre 1 comporte dans une représentation polyphase une fonction de transfert en z matricielle H(z) = [H; i (z)] , comportant N voies (N entrées et N sorties), et un sélecteur de voie d'entrée synchrone avec un sélecteur de voie de sortie. Les sélecteurs de voies d'entrée et de sortie sont commandés à une cadence égale à Fe, Fe étant la fréquence d'échantillonnage du signal d'entrée x(n), et N étant la période du filtre. La fréquence du signal dans les voies d'entrée et de sortie est donc égale à Fe/N. Ce filtre transforme le signal d'entrée x(n) compris dans la bande de fréquences 25 [fmin, fm.] en un signal de sortie y(n) occupant une bande de fréquences B avec au moins une sous-bande B' effacée. Ce filtre présente une réponse impulsionnelle variant en fonction du temps avec une périodicité N et peut être représenté par la relation suivante: + o y(n) _ E h(n, k) É x(n k) (1) k= co h(n, k) étant la réponse impulsionnelle du filtre à l'instant n. Comme le filtre présente une période égale à N, h vérifie l'égalité suivante: h(n+N, k) = h(n, k) (2) Dans la représentation polyphase à N composantes des signaux d'entrée x(n) et 2873876 -5- de sortie y(n) du filtre 1, les signaux x(n) et y(n) peuvent se décomposer de la manière suivante: x(n) = [xo(n), ... , xN-1(n)] et y(n) _ [yo(n), ... , yN-1(n)I (3) avec xk(n) = x(nN + k) et yk(n) = y(nN + k), k étant compris entre 0 et N-1. La transformée en z de la composante xk(n) peut être obtenue à l'aide de la formule suivante: +co Xk (Z) _ E xk (n)Z n (4) n = oo On montre aisément que la transformée en z Yk(z) de la composante polyphase 10 k du signal y(n) s'écrit en fonction de la transformée en z du signal x(n) de la manière suivante: N 1 Yk (Z) = E, Xi (z) É Hk i (Z) (5) j=o avec: X(Z) = [Xo(z) X1(z) ... XN_I (z)] (6) Y(z) = [Yo(z) Y1(z) ... YN-1(z)], (7) H00(z) Ho 1(z) Ho,N-1(z) H1,o (z) H1,1 (z) H1,N I (z) H(z) _ (8) HN I,O(Z) HN_I(Z) .. HN 1,N I(Z) et (9) Y(z) = H(z) X(z) Un tel filtre peut également être représenté sous forme fréquentielle en transformée en z de la manière suivante: N 1 Y(z) = Tp(zWp) X(zW") (10) p=o dans laquelle W est la racine Nième de l'unité : W = exp(-2iit/N), et les termes 25 Tp(z) représentent des filtres modulateurs tels que: N 1 N 1 Tp(z) = z '1, Hjk(ZN)ÉZk Wp' (11) N i=o k=o La figure 3 est une représentation bifréquentielle d'un exemple de filtre linéaire périodique de période N = 4 qui permet d'exprimer le spectre du signal d'entrée x(n). Cette figure montre en traits épais les zones où les réponses fréquentielles des filtres modulateurs Tp(z) sont non nuls (supérieures à un seuil donné). Pour pouvoir déterminer en réception le signal x(n) en fonction du signal y(n) reçu, le filtre 1 doit être inversible. En outre, dans des applications avec transmission par satellite, les signaux sont amplifiés de manière non linéaire. Dans ce type d'applications, le signal après étalement doit être à enveloppe complexe constante. Ces conditions sont remplies si l'on choisit un modulateur à déplacement de phase et un filtre linéaire périodique numérique ayant en représentation polyphase une matrice dont seules les composantes Hk,k(z) = Hk(z) sont non nulles, et choisies égales à Wg(k), g(k) (k variant de 0 à N-1) étant des coefficients réels. Il en résulte que: 1 N 1 T (z) (z) = Tp = E Wg(k)+pk N k=o Un tel filtre est illustré par la figure 4. Le filtre inverse utilisable en réception H'k(z) donnant x(n) en fonction de y(n) pour la composante k est alors de la forme: H'k(z) = W g(k) (13) Par ailleurs, en transformée en z, le signal d'origine X(z) peut être obtenu à partir du signal reçu Y(z) par la relation suivante: N 1 X(z) = ETp(zWp) É Y(zWp) (14) p=0 Il en résulte que: 1 N 1 Tp(z) = Tp = W-g(k)+ pk (15) N k=O Selon l'invention, il est possible de déterminer les paramètres g(k) du filtre de manière à obtenir au moins une sous-bande B' effacée. Or, pour annuler dans la sous-bande de fréquences B' = [bmin, bmax] un signal de bande [fmin, fmax], on démontre qu'il faut annuler les termes Tp pour toutes les valeurs de p comprises dans le segment [purin, pmax] défini modulo N tel que: _E/fmax b bmin min N (16) B N (12) 2873876 -7- _ E/Emin bmax N\ -1 B, N Dans ces formules, l'expression "[a]N" signifie "a modulo N", c'est--à-dire représente le reste de la division entière de a par N, et "E(a)" représente la partie entière de a. Inversement, lorsque la bande de fréquences B' est occupée par un signal de brouillage, il s'agit d'éviter que le signal récupéré par le filtre LTPV inverse du filtre 1 ne contienne des résidus de brouillage dans la bande [fmin, furax]. Il s'avère que cette condition est vérifiée si les filtres modulateurs inverses T'p sont nuls pour toutes les valeurs de p comprises dans l'intervalle [p'min, p'max] On démontre que: fmin bmax purin = B N +11 N Lorsque le filtre linéaire périodique numérique choisi est directement utilisé pour ré-étiquetter les symboles complexes en bribes à un rythme Rc, on choisit la bande de fréquences [fmin, furax] égale à [-Rc/2, Rc/2] et la bande B égale à Fe. Pour annuler les termes Tp du filtre 1 pour une valeur pl de p (nombre entier), on se propose d'annuler deux à deux les quantités Wg(k)+pik Cette opération consiste à regrouper les éléments de l'ensemble EN des entiers compris entre 0 et N-1 en couples (ka, k,) tels que: Wg(ko)+Ptko + wg(ki)+Piki = 0 (20) ce qui permet d'obtenir N/2 conditions portant sur les N coefficients g(k) du 25 filtre: g(ko(l)) g(ki(l))= [pl (kl (j) ko (j + 2 (21) IN pour toutes les valeurs de j comprises entre 0 et 2 -1. Les valeurs de g(k) sont donc liées deux par deux par 2 équations linéaires. L'algorithme qui suit est applicable au cas particulier où la période N du filtre est choisie de la forme N = 2 , n étant un nombre entier. P max (17) P max E/ffax bmin N B N (18) (19) 2873876 -8- Pour trouver un premier ensemble de couples (ko, k1) d'éléments de l'ensemble EN des entiers compris entre 0 et N 1, on choisit un entier al compris entre 0 et n 1, tel que k1(j) = ko(j) + 2a' pour toutes les valeurs de j comprises entre 0 et N 1 = 2n-' 1. 2a' représente donc la "distance" entre les éléments ko, kl 2 formant les couples (ko, k1). Il s'avère que le fait d'annuler un terme Tp pour une certaine valeur de p entraîne l'annulation d'autres termes Tp. En effet, si les termes Tp sont nuls pour deux valeurs pa et Pb, on a la relation suivante: Pa(k1G) ko(1))+ = Pb(kia) ko(i))+ (22) 2 N 2 N pour toutes les valeurs de j comprises entre 0 et 2 1. (23) Si k1(j) = ko(j) + 2a', il résulte de l'équation (22) que: [(Pb Pa)2a']N = 0, soit Pb = [Pa]2n-a, si N=2n. Par conséquent, si l'on annule un terme Tp du filtre pour une valeur de p, on annule les termes Tp, du filtre pour toutes les valeurs p' = [p+jÉ2na']N, ce qui réduit la qualité d'étalement si l'effacement des sousbandes ainsi effacées est indésirable. Si par exemple on annule le terme To, on annule également tous les termes TU2n-a,]N. Si on choisit a1 = n 1, tous les termes T2p sont aussi annulés. Pour effacer une bande de fréquences B' _ [b,,,;n, b,,,ax], il faut que les termes Tp du filtre s'annulent pour toutes les valeurs de p consécutives situées dans l'intervalle [p,n,p, Pmax] É Il s'agit donc trouver de nouveaux couples d'éléments de l'ensemble EN qui doivent être tous distincts des précédents pour conserver l'annulation de T. A cet effet, on regroupe deux par deux les couples utilisés à l'itération précédente en quadruplets, en choisissant une nouvelle valeur a2 différente de al telle que la "distance" entre les couples regroupés soit égale à 2a2, c'est-à-dire que l'on regroupe en quadruplets les couples (ko, k1) et (k'o, k',) tels que: k'o=ko+2aZ. (24) Cette condition permet d'obtenir par quadruplet (ko, k1, k'o, k'1) deux nouvelles équations de la forme (21), qui sont en fait équivalentes: g(ko (i)) g(ko G)) = [Pooi) ko G)) + 2 (25) N g(ki (j)) g(ki (l)) = [p0(k; (j) ki (l)) + 2 (26) N soit une seule équation par quadruplet, c'est-à-dire 4 = 2n 2 équations linéaires 5 à N inconnues. A l'itération suivante, on regroupe deux à deux en groupes de 8 entiers les quadruplets (groupes de 4 entiers) tels que le premier élément du second groupe est distant de 2a3 du premier élément du premier groupe, chaque groupe ainsi obtenu fournit une équation supplémentaire sur les coefficients g(k) à déterminer. Cette itération peut être répétée autant de fois que nécessaire dans une certaine limite, pour effacer totalement une bande B' de la largeur souhaitée, en prenant soin à chaque itération j de choisir toutes différentes des valeurs a; de distance entre les groupes regroupés, telles que k'o = ko + tai. A l'itération j, le degré de liberté pour le choix des coefficients g(k) est alors égal à 2 -i. On démontre que la largeur maximale de la bande susceptible d'être effacée correspond à un nombre maximum de filtres consécutifs nuls égal à 4(n 3) + 3. Bien entendu d'autres formes de valeurs de la période N du filtre peuvent être choisies. Ainsi, si l'on choisit N de la forme 2 + 4, on peut constater que l'étalement produit par le filtre est de meilleure qualité (moins d'annulations indésirables de termes Tp). Les figures 5 à 13 illustrent les performances du dispositif représenté sur la figure 1 sous la forme de courbes spectrales des signaux x(n) et y(n) en fonction de la fréquence normalisée (f/Fe en abscisse), dans le cas où le modulateur numérique d'entrée 2 présente une constellation à deux points (M=2), le suréchantillonneur 4 présente un facteur K égal à 4, et le filtre de sortie 5 est en racine de cosinus surélevé. Les spectres du signal y(n) sur les figures 5 à 10 sont obtenus avec un filtre linéaire périodique 1 de période N égale à 210. Les spectres des figures 5 à 7 sont obtenus avec un faible étalement (le suréchantillonneur 3 présente un 2873876 -10- facteur d'étalement L égal à 8). Sur ces figures, les courbes spectrales du signal de sortie y(n) sont obtenues lorsque les termes Tp du filtre 1 sont annulés au voisinage de la fréquence normalisée -0,07, pour respectivement trois, quatre et neuf valeurs de p consécutives, ce qui correspond à un degré de liberté D pour le choix des coefficients g(i) du filtre 1 respectivement égal à 27, 26 et 2. Les spectres des figures 8 à 10 sont obtenus avec un fort étalement (le suréchantillonneur 3 présente un facteur L égal à 128). Sur ces figures, les courbes spectrales du signal de sortie y(n) sont également obtenues lorsque les termes Tp du filtre 1 sont annulés au voisinage de la fréquence normalisée -0,07 pour respectivement trois, quatre et neuf valeurs de p consécutives, ce qui correspond à un degré de liberté D pour le choix des coefficients g(i) du filtre 1 respectivement égal à 27, 26 et 2. Sur les figures 11 à 13, les spectres sont obtenus avec un filtre linéaire périodique 1 dont la période N est égale à 2' et avec un faible étalement (L = 8). Sur ces figures, les courbes spectrales du signal de sortie y(n) sont également obtenues lorsque les termes Tp du filtre 1 sont annulés au voisinage de la fréquence normalisée -0,07 pour respectivement trois, quatre et six valeurs de p consécutives, ce qui correspond à un degré de liberté D pour le choix des coefficients g(i) du filtre 1 respectivement égal à 24, 23 et 2. -11-REVENDICATIONS 1. Procédé d'étalement de spectre d'un signal d'entrée numérique (x(n)), caractérisé en ce qu'il comprend des étapes consistant à : choisir des coefficients d'un filtre linéaire périodique numérique (1) réversible, de manière à ce qu'il étale le spectre du signal d'entrée (x(n)) pour occuper une bande de fréquences étalée (B) à l'exception d'au moins une sous-bande (B') de fréquences effacée, incluse dans la bande de fréquences étalée, et traiter le signal d'entrée à l'aide du filtre linéaire périodique numérique de manière à obtenir un signal de sortie (y(n)) à spectre étalé. 2. Procédé selon la revendication 1, caractérisé en ce que la sortie du filtre linéaire périodique (1) est apte à 15 conserver l'enveloppe du signal d'entrée (x(n)). 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que le filtre linéaire périodique numérique présente une transformée en z dans une représentation fréquentielle de la forme suivante: N-1 Y(z) = Tp É X(zW p) p=0 dans laquelle: 1 N-1 Tp = Wg(k)+pk W = exp N k=0 les nombres entiers k compris entre 0 et N 1 sont des coefficients à déterminer. 4. Procédé selon la revendication 3, caractérisé en ce que le choix des coefficients g(k) comprend des étapes consistant à annuler les termes Tp du filtre pour toutes les valeurs de p telles que _E/fmax Bbmin N p _E /f min bmax N _l 1N B N, expression dans laquelle [fmin, fmax] représente la bande de fréquences du signal numérique d'entrée (x(n)), B représente la bande de fréquence étalée, [b,,,i,,, bmax] représente la sous-bande de fréquences (B') effacée, incluse dans la bande de fréquences B, E(a) représente la fonction partie entière et [a]N représente la 2i n N, N est la période du filtre, et g(k) pour -12-fonction "modulo N". 5. Procédé selon la revendication 4, caractérisé en ce que l'annulation d'un terme Tp du filtre linéaire périodique (1) 5 consiste à regrouper tous les entiers compris entre 0 et N-1 en couples (ka, kt) tels que: Wg(ko)+pko + wg(ki)+pk, = 0. 6. Procédé selon l'une des revendications 1 à 5, 10 caractérisé en ce que le filtre linéaire périodique (1) est applicable à un signal d'entrée (x(n)) à enveloppe constante. 7. Procédé selon l'une des revendications 1 à 6, caractérisé en ce que le signal d'entrée (x(n)) est préalablement 15 suréchantillonné. 8. Procédé selon la revendication 7, caractérisé en ce qu'avant d'être suréchantillonné, le signal d'entrée (x(n)) est traité par un modulateur numérique présentant une constellation à m points. 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce que le signal de sortie (y(n)) du filtre linéaire périodique (1) est suréchantillonné. 10. Procédé selon la revendication 9, caractérisé en ce que le signal de sortie suréchantillonné est traité par un filtre passe-bas (5) en racine de cosinus surélevé. 11. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que le filtre linéaire périodique (1) présente une période de la forme 2n, n étant un nombre entier positif.