FR 2919135 A1 20090123 FR 0756530 A 20070717 MÉTHODE DE DÉTECTION DE SIGNAUX CYCLOSTATIONNAIRES DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne la détection de signaux cyclostationnaires. Elle trouve application dans les domaines aussi variés que celui des systèmes radio opportunistes, celui du traitement d'antennes ou du tatouage radio. ÉTAT DE LA TECHNIQUE ANTÉRIEURE L'encombrement de plus en plus important du spectre a conduit à envisager des systèmes de télécommunication capables de coexister avec les systèmes à bandes de fréquence déjà allouées, dits primaires. Deux stratégies de coexistence font actuellement l'objet d'importantes recherches. La première consiste à utiliser un niveau de signal très faible grâce à un très fort étalement de spectre, c'est la voie suivie par les systèmes ultra-large bande encore dénommés UWB (Ultra Wide Band). La seconde consiste à utiliser de manière opportuniste une partie du spectre ponctuellement ou temporairement inoccupée, c'est la voie dite de radio opportuniste (ou Cognitive Radio). On trouvera une description de la radio opportuniste dans la thèse fondatrice de J. Mitola intitulée Cognitive radio : an integrated agent architecture for software defined radio , Royal Institute of Technology, Stockholm, PhD Dissertation, 8 Mai 2000. Avant de pouvoir émettre dans une bande donnée, il est nécessaire qu'un émetteur puisse déterminer si un signal radio, notamment un signal TDMA (Time Division Multiple Access) comme un signal GSM, ou un signal étalé spectralement de type CDMA (Code Division Multiple Access) ou encore un signal OFDM (Orthogonal Frequency Division Multiplexing) est présent dans cette bande. Etant donné que le signal étalé peut être noyé dans le bruit, un simple détecteur de puissance ne permet pas d'obtenir cette information de manière fiable. De manière générale, on distingue les méthodes de détection à l'aveugle pour lesquelles on ne dispose d'aucune information a priori sur le signal dont on veut détecter la présence et les méthodes d'estimation de type semi-aveugle pour lesquelles on connaît certains paramètres du signal à détecter. Une méthode de détection semi-aveugle de la présence d'un signal TDMA, CDMA ou OFDM à partir d'un critère de cyclostationnarité a été décrite dans l'article de M. Oner et al. intitulé Cyclostationnarity based air interface recognition for software radio systems publié dans Proc. of 15th IEEE International Symposium on Personal, Indoor and Mobile Radio Communications, Sept. 2004, pp. 1947-1951. Plus précisément, cette méthode permet de tester à partir de la fonction d'autocorrélation du signal reçu si ce dernier présente une fréquence cyclique caractéristique d'une interface radio prédéterminée. Cette méthode ne fonctionne cependant pas dans un contexte aveugle. En outre, même dans un contexte semi-aveugle, lorsque le signal à détecter peut prendre différentes valeurs de paramètres de transmission, par exemple un signal CDMA pouvant avoir différents facteurs d'étalement spectral, cette méthode ne permet pas d'obtenir un taux de fausse alarme ou CFAR (Constant False Alarm Rate) prédéterminé. La détection de signaux cyclostationnaires est utilisée dans des domaines très variés, par exemple celui du tatouage radio. Le tatouage est un signal de très faible puissance que l'on insère dans un signal audio pour y introduire des données cachées, aux fins notamment de traquer les copies illicites. L'article de L. de C.T. Gomes et al. intitulé Tatouage audio exploitant des propriétés de cyclostationnarité publié dans la revue Traitement du Signal, 2001, Vo. 19, N 1, pp. 1-9, propose un signal de tatouage cyclostationnaire codant des données cachées à l'aide d'une pluralité de fréquences cycliques. La détection publique du tatouage audio est réalisée au moyen d'une détection desdites fréquences cycliques dans le signal audio. Toutefois, cette méthode de détection ne permet pas de déterminer, pour une grande variété de signaux audio, si un tatouage radio est présent avec un taux de fausse alarme prédéterminée. Le but de la présente invention est de proposer une méthode robuste permettant de déterminer avec un degré de fiabilité prédéterminé, dans un contexte aveugle ou semi-aveugle, si un signal cyclostationnaire est présent dans un signal reçu ou lu, ce pour des conditions de faible rapport signal à bruit. EXPOSÉ DE L'INVENTION La présente invention est définie, selon un premier mode de réalisation, par une méthode de détection d'un signal cyclostationnaire dans un signal à analyser, reçu d'un émetteur par un récepteur ou lu d'un support d'enregistrement par un dispositif de lecture, dans laquelle on estime la valeur d'une fonction de discrimination J s'exprimant comme une forme quadratique des coefficients de corrélation cyclique du signal à analyser pour un ensemble de décalages temporels non nuls et un ensemble de 4 fréquences U V,4J 6 ~ déterminer cycliques, puis l'on compare la valeur à un taux de mauvaise détection ri pour si ledit signal cyclostationnaire est 1-y présent dans le signal à analyser, où *,x) est la fonction gamma inférieure incomplète normalisée, v est le produit du cardinal dudit ensemble de décalages temporels non nuls et du cardinal dudit ensemble de fréquences cycliques, U est la largeur de la fenêtre d'autocorrélation sur laquelle sont calculés les coefficients de corrélation cyclique, et 64 est le carré de la variance du bruit. Selon une variante, ledit ensemble de fréquences cycliques se réduit à la fréquence nulle et la valeur de la fonction de discrimination est calculée par : où R (m) est la valeur de la fonction d'autocorrélation du signal à analyser pour un décalage temporel non nul m EIM, et lm étant ledit ensemble de décalages temporels non nuls. Avantageusement, la valeur de la fonction d'autocorrélation pourra être estimée à partir des échantillons s(u) du signal à analyser par : R (m)=ùEs U u=o où U est la largeur de la fenêtre d'autocorrélation et m est un décalage temporel exprimé en nombre de périodes d'échantillonnage. Selon une seconde variante, la fonction de discrimination est calculée par : 1 2 J= A 2R où R est le vecteur de taille v constitué des coefficients de corrélation cyclique R Ck(m) pour lesdits 20 décalages temporels non nuls et lesdites fréquences cycliques, et A est une matrice proportionnelle à la matrice de corrélation de ces coefficients. Selon une troisième variante, la fonction de discrimination pourra être calculée par : 1 U-1 25 où R Ck(m), me lm sont les coefficients de corrélation cyclique, lm est ledit ensemble de décalages temporels non nuls et IK est ledit ensemble de fréquences cycliques. Avantageusement, les coefficients de corrélation cyclique pourront être estimés à partir des échantillons s(u) du signal à analyser par : 1 U-1 Rak (m) = Ls(u+m)s*(u)eahmku U où îk est une fréquence cyclique appartenant à IK et m est un décalage temporel, exprimé en nombre de périodes d'échantillonnage. L'invention est également définie, selon un second mode de réalisation, par une méthode de détection d'un signal cyclostationnaire dans une pluralité Q de signaux à analyser, reçus d'un émetteur par une pluralité de récepteurs ou lus d'un support d'enregistrement par une pluralité de dispositifs de lecture, dans laquelle, pour chacun desdits signaux à analyser, on estime la valeur d'une fonction de discrimination Jq s'exprimant comme une forme quadratique des coefficients de corrélation cyclique de ce signal pour un ensemble de décalages temporels non nuls et un ensemble de fréquences cycliques, et l'on iQ u compare la valeur 1ùy Lvq,L 4 Jq à un taux de mauvaise 0=1 q=16q / détection ri pour déterminer si ledit signal cyclostationnaire est présent dans lesdits signaux à analyser, où y(v,x) est la fonction gamma inférieure incomplète normalisée, vq est le produit du cardinal dudit ensemble de décalages temporels non nuls et du cardinal dudit ensemble de fréquences cycliques pour le signal à analyser sq(t), Uq est la largeur de la fenêtre d'autocorrélation sur laquelle sont calculés les coefficients de corrélation cyclique et 6q est le carré de la variance du bruit, relatifs au signal à analyser. Selon une première variante, pour chacun desdits signaux à analyser, ledit ensemble de fréquences cycliques se réduit à la fréquence nulle et la valeur de la fonction de discrimination est calculée, pour ce signal, par : a Jq = R (m ) mEIMq où R (m) est la valeur de la fonction d'autocorrélation du signal à analyser sq(t), pour un décalage temporel non 20 nul me 1mq, IMq étant ledit ensemble de décalages temporels non nuls pour ce signal. Avantageusement, pour chacun desdits signaux à analyser, ladite valeur de la fonction d'autocorrélation pourra être estimée à partir des échantillons s(u) de ce signal par : 0 1 U-1 Rsq (m) = ù ~sq (u + m)s (u) U ,i=0 25 où U est la largeur de la fenêtre d'autocorrélation et m est un décalage temporel exprimé en nombre de périodes d'échantillonnage. Selon une seconde variante, pour chaque signal sq(t) à analyser, la fonction de discrimination est calculée par : 1 2 Jq = Aq 2Rq où Rq est le vecteur de taille vq constitué des coefficients de corrélation cyclique Rak(m) pour lesdits décalages temporels non nuls et lesdites fréquences cycliques, et Aq est une matrice proportionnelle à la matrice de corrélation de ces coefficients. Selon une troisième variante, pour chaque signal sq(t) à analyser, la fonction de discrimination est calculée par : Jq = L mEIMq CCkEIK, Rak (m) s, 2 où Rak(m), m EIMq sont les coefficients de corrélation cyclique du signal sq(t), IMq est ledit ensemble de décalages temporels non nuls et IKq est ledit ensemble de fréquences cycliques relatifs à ce signal. Avantageusement, pour chaque signal sq(t) à analyser, les coefficients de corrélation cyclique pourront être estimés à partir des échantillons sq(u) de ce signal par : 1 U-1 _ ù Ls(u + m)s* lu\e2Z7L0Gku U u=o où îk est une fréquence cyclique appartenant à IK et m est un décalage temporel, exprimé en nombre de périodes d'échantillonnage. BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de 15 réalisation préférentiel de l'invention fait en référence aux figures jointes parmi lesquelles : La Fig. 1 représente un organigramme de la méthode de détection d'un signal cyclostationnaire, selon un premier mode de réalisation de l'invention ; 20 La Fig. 2 représente un organigramme de la méthode de détection d'un signal cyclostationnaire, selon un second mode de réalisation de l'invention ; La Fig. 3 représente un organigramme de la méthode de détection d'un signal cyclostationnaire, 25 selon un troisième mode de réalisation de l'invention ; La Fig. 4 représente un organigramme de la méthode de détection d'un signal cyclostationnaire, selon un quatrième mode de réalisation de l'invention. 10 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Nous considèrerons dans la suite le cas d'un signal cyclostationnaire noyé dans du bruit. Nous désignerons ici par signal cyclostationnaire (au sens large) un signal susceptible d'être représenté par une fonction aléatoire s(t) telle que sa fonction d'autocorrélation définie par : RS (t,ti) = E{s(t)s* (t ùt )} (1) où E{} désigne l'espérance mathématique et est une fonction périodique. Dans ce cas, la fonction d'autocorrélation peut être développée en série de Fourier . 15 RS (t ,t ) _ LRâk (t 2iIakt k l~ où les coefficients complexes R Ck(ti) sont appelés coefficients de corrélation cyclique et les valeurs 20 réelles âk sont appelées fréquences cycliques du signal. Si l'on se place dans le cas d'un signal cyclostationnaire à temps discret, l'expression (2) devient . 25 RS(n,m) _ LRak (m2i7cakn (3) k (2) où t = nTe = mTe et âk = îk où Te est la période e d'échantillonnage. Dans la suite, nous nous référerons indifféremment à l'une ou l'autre formulation sans perte de généralité. Nous considérons dans la suite un signal à analyser, reçu d'un émetteur par un ou plusieurs récepteurs, ou bien un signal lu d'un support d'enregistrement, par un ou plusieurs dispositifs de lecture. La méthode de détection selon l'invention permet de déterminer avec un degré de fiabilité donné, c'est-à-dire, de manière équivalente, avec un taux de mauvaise détection donnée, si un signal cyclostationnaire est présent dans le signal à analyser. Nous supposerons dans la suite que le signal cyclostationnaire est noyé dans un bruit blanc gaussien centré, avec un rapport signal à bruit très faible, autrement dit très négatif lorsqu'il est exprimé en dB. Nous noterons conventionnellement Ho l'hypothèse selon laquelle le signal à analyser ne contient que du bruit et H1 l'hypothèse selon laquelle le signal à analyser contient un signal cyclostationnaire noyé dans ce bruit. La méthode de détection selon l'invention utilise de manière avantageuse une fonction de discrimination définie comme une forme quadratique des coefficients de corrélations cycliques du signal à analyser, les coefficients étant relatifs à une pluralité de décalage temporels non nuls et une pluralité de fréquences cycliques. La fonction de discrimination pourra prendre la forme : J= L L R `k (m)2 (4) mEIM îkEIK où IM est un ensemble discret de valeurs non nulles de décalage temporel de la fonction d'autocorrélation et 10 IK est un ensemble discret de valeurs de fréquences cycliques. Le cas échéant, l'ensemble IM et/ou l'ensemble IK peu(ven)t être réduit (s) à un singleton. On notera M = card(IM) et K = card(IK) . 15 En théorie, lorsque le signal à analyser ne présente pas de cyclostationnarité, c'est-à-dire lorsque l'on se trouve dans l'hypothèse Ho, la valeur de J est nulle. En pratique, on estime J en estimant les 20 coefficients de corrélation cyclique dans le domaine temporel de la manière suivante : U- 1 1 * 8 ( m ) ) _ ù Es(u+m)s (u)e2z"ku U u=o 25 où U est le temps d'observation, aussi appelé largeur de fenêtre d'autocorrélation. L'estimée J de J étant alors obtenue simplement par : (5) J = L L F?., (m)2 13 (6) mEIM îkEIK Lorsque l'estimée J est supérieure à un seuil 5 prédéterminé Io, on pourra en conclure qu'un signal cyclostationnaire est présent dans le signal à analyser. En revanche, lorsque l'estimée J est inférieure à ce seuil, on pourra en conclure que le signal à analyser ne contient que du bruit. 10 Nous donnerons ci-après quelques exemples non limitatifs de fonction de discrimination lorsque le signal cyclostationnaire est un signal étalé spectralement ou un signal OFDM. 15 En général, un signal étalé spectralement, reçu d'un émetteur par un récepteur peut s'exprimer de la manière suivante : L-1 20 s(t) _ ~apEc,eg~tû~T~ pTs)+b(t) (7) p ou c,e, = 0,..,Lû1 est la séquence d'étalement spectral, Tc est la période bribe (chip), T, = LT', est la période symbole, ap sont les symboles d'information transmis 25 par l'émetteur, g(t) est la convolution de l'impulsion de mise en forme d'un symbole par la réponse impulsionnelle du canal entre l'émetteur et le récepteur, b(t) est une fonction aléatoire décrivant le bruit. De manière similaire, un signal OFDM peut s'écrire, après démodulation en bande de base dans le récepteur, sous la forme suivante : y L t Nù1 2in ' (tùDTT ùk(N+D )Tc ) s(t)= N ~g(t-k(N+D)T~)e2in4f.Ean+kNe NT, k n=0 (8) où E est la puissance du signal, N est le nombre de porteuses du multiplex OFDM, an sont les symboles d'information appartenant à un alphabet de modulation, typiquement BPSK, QPSK ou QAM, 1/T, est le débit des symboles d'information, l'intervalle inter-porteuse valant 1/NT,, D est la taille du préfixe cyclique exprimée en fraction de la durée utile ]=NT, g(t) est une impulsion de mise en forme des symboles OFDM ayant un support temporel [0,(N+D)TA destiné à apodiser le spectre du signal, 4f est un résidu de porteuse dû à l'imperfection de la démodulation RF. Lorsque la méthode de détection selon l'invention opère dans un contexte aveugle, c'est-à-dire lorsque l'on ne dispose d'aucune information a priori sur le signal cyclostationnaire qui pourrait se trouver noyé dans le bruit, on se contente alors de calculer la fonction de discrimination pour la fréquence cyclique nulle. Autrement dit, l'ensemble des fréquences cycliques se réduit à un singleton, IK ={0} Pour détecter à l'aveugle la présence d'un signal étalé spectralement noyé dans du bruit, on pourra prendre comme ensemble de valeurs de décalage temporel, l'ensemble défini par : max , Mmin,}u{Mmin,..maxJ où et M,,,ax sont des entiers positifs non nuls. Typiquement, =1 et Mmax =10 . 10 Pour détecter à l'aveugle un signal OFDM noyé dans du bruit, on prendra comme ensemble de valeurs de décalage temporel, le singleton : lm= lm= 15 où Tu est la durée utile, et Te est la période d'échantillonnage. Dans un contexte semi-aveugle on échantillonnera 20 le signal à analyser, CDMA ou OFDM à sa fréquence chip, 1 soit ù . Pour un signal CDMA dont on connaît le facteur d'étalement spectral L, on pourra prendre respectivement comme ensemble de fréquences 25 cycliques et ensemble de décalages temporels : max , Mmin,}u{Mmin,..maxJ De manière similaire, dans le cas de la détection d'un signal OFDM, on pourra prendre respectivement comme ensemble de fréquences cycliques et ensemble de décalages temporels : N lm = et IK = D Revenant au cas général, la méthode de détection selon l'invention permet avantageusement de décider entre les hypothèses H1 et Ho, autrement dit de décider si le signal à analyser contient ou ne contient pas un signal cyclostationnaire, ce avec un taux de mauvaise détection prédéterminé. On peut montrer (voir annexe I) que, si l'hypothèse Ho est réalisée, les estimés Rak des coefficients de corrélation cyclique suivent une loi 64 gaussienne centrée de variance . On peut montrer en U outre que (voir annexe II), sous la même hypothèse, les estimés respectifs de deux coefficients de corrélation cyclique s et Rak2(m2) sont décorrélés, c'est-à- dire sont tels que E Ra (ml)(Rak2 (m2))* Ho}=O si m1 ≠m2 . En revanche, si m1= m2, les estimés de ces coefficients, sas (ml) et Raki (m2), sont corrélés. Nous envisagerons successivement le cas d'un contexte aveugle et celui d'un contexte semi-aveugle.25 Dans un contexte aveugle, la fonction de discrimination proposée se réduit à : car l'on ne connaît pas a priori les fréquences cycliques du signal cyclostationnaire que l'on recherche. L'estimée J s'exprime comme un somme de M 10 variables gaussiennes centrées, sa densité de probabilité conditionnelle suit une loi du x2, soit : U 1 ( U`Mùlù3 Pr(JHo)= 4 `14 e 6 (Mù1)(2 6 On en déduit que sa fonction de répartition est donnée par : U Pr(J Jo Ho) = 7 M, Jo 4 6 / et, par conséquent, U Pr(J>IoHo)=1ù7 M,Io 4 25 où y est la fonction gamma incomplète inférieure normalisée, c'est-à-dire : 6 72 (10) 15 20 (12) 18 7(M,x) _ (Ml 1)i Jo tM-ie-tdt D'après l'expression (12), on voit que si : (13) (14) 1ùy où ri est un taux prédéterminé de mauvaise détection, on peut conclure que, si J > Jp, l'hypothèse Ho est 10 réalisée avec une probabilité inférieure à r. En d'autres termes, la probabilité d'absence de cyclostationnarité sachant que J > Jp est égale à t~. La Fig. 1 représente de manière schématique la 15 méthode de détection d'un signal cyclostationnaire selon un premier mode de réalisation de l'invention. A l'étape 110 le signal à analyser s(t) est filtré puis démodulé en bande base, si nécessaire. Le signal ainsi filtré est échantillonné en 120 à 1 20 la fréquence de Nyquist ù . Te A l'étape 130, on estime les coefficients de coefficients cycliques du R (m) grâce à l'expression (5), c'est-à-dire : 25 ) (15) A l'étape 140 on estime la valeur de la fonction de discrimination selon l'expression (9). On note J cette estimée. On calcule, en 150, la valeur de l'expression U 1ùy M,J 4 et on la compare en 160 à un taux de mauvaise (ou fausse) détection 11 prédéterminé. Si cette valeur est inférieure à 1 , on conclut en 172 à la présence d'un signal cyclostationnaire. Dans la négative, on conclut en 171 à l'absence de cyclostationnarité du signal à analyser. Selon une seconde variante de réalisation, on dispose de plusieurs versions sq(t), q = 1,..,Q d'un signal susceptible de contenir un même signal cyclostationnaire. Ces différentes versions peuvent être analysées conjointement pour déterminer si un signal cyclostationnaire est effectivement présent ou absent. Les signaux à analyser peuvent être, selon le 20 cas, les signaux reçus par les différentes antennes de réception d'un système SIMO (Single Input Multiple Output), les signaux reçus par une pluralité de capteurs ou encore les signaux lus d'un support d'enregistrement au moyen d'une pluralité de 25 dispositifs de lecture. La Fig. 2 illustre de manière schématique la méthode de détection d'un signal cyclostationnaire selon un second mode de réalisation de l'invention. Ce 30 mode de réalisation utilise une pluralité Q de signaux 6 à analyser et comprend Q voies de traitement, chaque voie q opérant sur un signal sq(t), et comprenant des étapes 210q à 240q identiques aux étapes 110 à 140 de la Fig. 1. Autrement dit, chaque voie permet d'obtenir une estimée Jq de la fonction de discrimination Jq relative au signal sq(t). A la différence du premier mode de réalisation, on utilise une fonction de discrimination synthétique : QU J=E 4Jq q=16q où Mq, Uq et 6q sont respectivement le nombre de décalage temporels non nuls, la largeur de la fenêtre d'autocorrélation, et le carré de la puissance de bruit relatifs au signal sq(t). La fonction de répartition conditionnelle est alors la suivante : (16) Pr(J>JI) ,q=1,..QHO)=1ùy (Q Mq,Jo \q 1 / (17) Une fois les estimées Jq calculées aux étapes Q QU 240q, on calcule en 250, la valeur de 1ùy EMq,E 4 Jq q=l q=16q et on la compare en 260 à un taux de mauvaise détection ri prédéterminé. Selon que cette valeur est inférieure ou supérieure à 1 , on conclut respectivement, en 272 ou en 271, à la présence ou à l'absence de cyclostationnarité dans les signaux à analyser. Nous supposerons dans la suite que la méthode de détection est utilisée dans un contexte semi-aveugle, c'est-à-dire que les fréquences cycliques du signal cyclostationnaire que l'on recherche sont connues a priori. Dans ce cas, la fonction de discrimination est définie par (6) où IK est constitué par l'ensemble de ces fréquences connues. Comme on l'a vu plus haut les estimés des coefficients de corrélation cyclique Rai (ml) et Raki (m2) sont décorrélés si ml #m2 et corrélés pour ml = m2 . Dans ce dernier cas, le coefficient de corrélation de ces coefficients vaut (cf. Annexe II) . {Ra (ml )(1\ak2 (m2 ))*64 Ho} _ (ak k2) (18) avec (îk,,îk2)= -ak2)(U-1) sin(~(ak ù îkz)U) (18' ) sin (7c k, -ak2 )) On introduit le vecteur R de taille v =MK des estimés des coefficients de corrélation cycliques : où (.)T désigne l'opération de transposition. On a alors . (19) où R est la norme euclidienne de R. La matrice de corrélation du vecteur R sous l'hypothèse H0 n'est autre que : 4 r=6 A U avec : (20) A= (A 0 0" 0 A 0 . 0 0 0 A i et (îI,î2) (al,aK )' 1 (î2,îK) (a2,aK) 1 i (21) A= On remarquera que la matrice de corrélation t, de taille vxv est diagonale par blocs puisque les estimés de coefficients cycliques sont décorrélés pour des décalages temporels distincts. 5 La méthode de détection en contexte semi-aveugle utilise une fonction de discrimination modifiée dont les coefficients de corrélation cyclique sont blanchis, autrement dit décorrélés, soit : 2 (22) Jwhit = 1 A 2R On est ainsi ramené à une somme quadratique de variables aléatoires gaussiennes centrées. De manière 10 similaire au cas aveugle on considère la fonction de répartition de Jwhit Pr(Jwhit JoHo)=1ùy (24) U~ La méthode de détection compare 1ûy v,Jp 4 avec 6 20 un taux de mauvaise détection prédéterminé 1 et en déduit si l'hypothèse Ho ou l'hypothèse H1 est vérifiée avec un taux d'erreur maximal rl. La Fig. 3 illustre de manière schématique la méthode de détection de stationnarité selon un troisième mode de réalisation de l'invention. Après avoir été éventuellement filtré et démodulé en bande de base à l'étape 310, le signal à analyser s(t) est échantillonné à une fréquence 1 , dépendant des Te caractéristiques spectrales du signal. Par exemple, si le signal est un signal étalé spectralement ou un signal OFDM, la fréquence d'échantillonnage ne sera 1 Tc On estime ensuite, à l'étape 320, les coefficients de corrélation cyclique à partir de l'expression (5), ce pour un ensemble IK de fréquences cycliques connues et un ensemble de décalages temporels IM. 15 On calcule ensuite, à l'étape 330, les coefficients (ak1,ak2) de la matrice A. Etant donné que A est hermitienne et que sa diagonale est constituée K(Kù1) de 1, seuls 2 coefficients sont à calculer. A l'étape 340, on calcule les coefficients de 1 20 corrélation blanchis A 2R et on obtient en 350 la valeur estimée de la fonction de discrimination modifiée Jwhit au moyen de l'expression (22) On en déduit à l'étape suivante, 360, la valeur de v Jwhit , valeur qui est ensuite comparée en 370 à 25 un taux d'erreur prédéterminé r~. autre que la fréquence chip 1-y Selon que ladite valeur est inférieure ou supérieure à 1 , on conclut, en 382 à la présence ou, en 381, à l'absence d'un signal cyclostationnaire. On remarquera d'après l'expression (18') que si la condition : 1 max ki ak ù ki ≠k2 2 U est vérifiée, autrement dit si la fenêtre d'autocorrélation est suffisamment large, les estimés des coefficients de corrélation cyclique sont décorrélés et la matrice de corrélation se réduit à (25) 4 r = I~xV Oë U est la matrice unité de taille vxv . Ivxv Dans ce cas, les étapes de calcul 330 et 340 sont inutiles. La Fig. 4 illustre de manière schématique la méthode de détection d'un signal cyclostationnaire selon un quatrième mode de réalisation de l'invention. Ce mode de réalisation utilise une pluralité Q de signaux à analyser et comprend Q voies de traitement. Les signaux à analyser peuvent être, comme pour le second mode de réalisation, les signaux reçus par les différentes antennes de réception d'un système SIMO (Single Input Multiple Output), les signaux reçus par une pluralité de capteurs ou encore les signaux lus d'un support d'enregistrement au moyen d'une pluralité de dispositifs de lecture. Chaque voie de traitement q comprend une pluralité d'étapes 410q à 450q identiques aux étapes 310 à 350 de la Fig. 3. En particulier, chacune des étapes 450q fournit une estimée Jwhit q de la valeur de la fonction de discrimination pour le signal sq(t), soit Jqwhit = Aq 2Rq où Rq est le vecteur de taille v 2 q constitué des coefficients de corrélation cyclique Rak(m) pour les décalages temporels non nuls et les fréquences cycliques du signal U sq(t), et Aq = 4 I'q (5q où Tq est la matrice de corrélation de ces coefficients. A la différence du troisième mode de réalisation, on calcule à l'étape 460 la valeur de : Q Q Uq Jwhit 1ùy Evq, E 4 q \q 1 q 16q vq= MgKq, où Mq et Kq sont respectivement le nombre de décalages temporels et le nombre de fréquences cycliques relatifs au signal sq(t) ; Uq et 64 la largeur de la fenêtre d'autocorrélation, et le carré de la 20 puissance de bruit relatifs à ce même signal. (Q u La valeur de 1ùy Evq,E 4 Jghit est ensuite 0=1 q=16q / comparée en 470 à un taux de mauvaise détection prédéterminér~. (26)15 Selon que cette valeur est inférieure ou supérieure à 1 , on conclut comme précédemment à la présence, en 482, ou à l'absence, en 481, de cyclostationnarité dans les signaux à analyser. La remarque précédente concernant la décorrélation des coefficients de corrélation cyclique s'applique ici également. En choisissant U suffisamment grand, les matrices de corrélation (ou certaines d'entre elles seulement) se réduisent à 4 hg = U Iv xvg où I~ xvg est la matrice unité de taille 9 9 vgxvq. Les étapes de calcul 330q et 340q sont alors inutiles. ANNEXE 2 Les coefficients de corrélation cyclique peuvent être estimés par : Rak( U-1 ) = 1 Ls(u+m)s*(u)e'Z7L0Gku U u=0 On se place dans l'hypothèse H0, autrement dit les s(u) sont des variables aléatoires de bruit i.i.d. D'après le théorème de la limite centrale, tend en loi vers une densité de probabilité gaussienne quand U tend vers l'infini. Sa moyenne E{Rak (m)} tend vers Rak (m) =0 pour m =0. Le moment d'ordre 4 de s(u) s'exprime en fonction 15 du cumulant d'ordre 4 : E{s(ul + m )s* (ul )s* (u2 + m)s(u2)}= cum (s (ul + m )s* (u1 )s*(u2 + m)s(u2 )) + t {s(ul + m )s* (ul )I{s* (u2 + m)s(u2)}+ t {s(ul + m)s(u2 )Jt {s* (ul )s* (u2 + m )l + t {s(ul + m )s* (u2 + m))E{s* (ul) s(u2 )} Le bruit étant gaussien : cum(s(ul +m)s*(ul)s*(u2 +m)s(u2))= 0 Le bruit étant circulaire, c'est-à-dire sa phase étant équirépartie sur [0,27L], on a : t{s(ul + m)s(u2 )}t{s* (ul )s* (u2 + m)}= 0 20 25 Par conséquent, la variance asymptotique s'exprime à partir des 2'd et 4ième termes . U-1 Uù1 lim /t{Rak(m)2}= lira L Et{s(ul+m)s*(ul) -2itaku1 t{s* (u2 + m)s(u u1=0 u2=O U-1 U-1 + lim 12 L L t {s(ul + m )s* (u2 + m)}E{s* (ui )s(u2) é Zi~ak (u uz ) U~~ U ul=0uZ=0 2 Le premier terme tend vers Rak(m) =0 quand U tend vers l'infini. Le second terme est non nul si et seulement si 10 ul =u2 et dans ce cas, puisque EIs(u)2}=62, le second 4 terme tend vers 6 (seuls U facteurs sont non nuls U sous le signe somme).5 ANNEXE II Nous montrerons dans cette annexe que deux 5 coefficients de corrélation cyclique Rak (m1) et Rak2 (m2) sont décorrélés si ml #m2 . Sinon, si ml = m2 , ces deux coefficients sont corrélés. 10 lim U-*oo 1t U-1 U-1 lim Lt{s(ul+m)s*(u1)s*(u2+m)s(u2) UU2 ,. =0u2=0 On a : /E}Raki (ml )(Ra k2 (m2 ))*1 -2in (aklul -ak2u2 e De même que dans l'annexe I, on peut exprimer ce moment d'ordre 4 en fonction du cumulant et de produits 15 de moments d'ordre 2. Le cumulant s'annule, un terme est nul du fait de la circularité du bruit, un terme tend vers t} Rak (m1)(Rak2 (m2 1=0. En définitive seul reste le terme : lim ~~{Rakl (mi )(Rak2 (m2~)* u c. U-1 U-1 lim 12 t L L t {s(ul + mi )s* (u2 + m2))E{S u2 =0 u2 =0 20 * 1)s(u2)e -2in (aklul -ak2u2 Etant donné que le bruit est décorrélé, E{s(ul +ml)s(u2 +m2)}= 0 , Si ml ≠ m2 . Par conséquent, dans ce cas lim (El ek' (m1)(lek2 (m2 ))* =O. En revanche si ml =m2, on a : ~~ Ra \ml )(Rak2 (m2 ))* _ 4 L e-2in (akiui -ak2u2 a U-1 U~~\ r) / U u=0 = 6 4 e-i~ (aklul -ak2u2 )(u-1) sin(n (aki ùîk2) ) sin &[ (akl ùîk2 )) Les coefficients Rak' (ml) et Rak2 (m2) sont donc bien corrélés. REVENDICATIONS 1. Méthode de détection d'un signal cyclostationnaire dans un signal à analyser, reçu d'un émetteur par un récepteur ou lu d'un support d'enregistrement par un dispositif de lecture, caractérisée en ce que l'on estime la valeur d'une fonction de discrimination J s'exprimant comme une forme quadratique des coefficients de corrélation cyclique du signal à analyser pour un ensemble (IM) de décalages temporels non nuls et un ensemble (IK) de fréquences cycliques, et que l'on compare la valeur U 1ùy v, 4J à un taux de mauvaise détection ri pour 6 ~ déterminer si ledit signal cyclostationnaire est présent dans le signal à analyser, où *,x) est la fonction gamma inférieure incomplète normalisée, v est le produit du cardinal (M) dudit ensemble de décalages temporels non nuls et du cardinal (K) dudit ensemble de fréquences cycliques, U est la largeur de la fenêtre d'autocorrélation sur laquelle sont calculés les coefficients de corrélation cyclique, et 64 est le carré de la variance du bruit. 2. Méthode de détection selon la revendication 1, caractérisée en ce que ledit ensemble de fréquences cycliques se réduit à la fréquence nulle et que la valeur de la fonction de discrimination est calculée par : 33J= où R (m) est la valeur de la fonction d'autocorrélation du signal à analyser pour un décalage temporel non nul me lm* , et lm étant ledit ensemble de décalages temporels non nuls. 3. Méthode de détection selon la revendication 2, caractérisée en ce que ladite valeur de la fonction d'autocorrélation est estimée à partir des échantillons s(u) du signal à analyser par : Ro s(m)=ù Es(u+m)s*(u) U u=o où U est la largeur de la fenêtre d'autocorrélation et m est un décalage temporel exprimé en nombre de périodes d'échantillonnage. 4. Méthode de détection selon la revendication 1, 20 caractérisée en ce que la fonction de discrimination est calculée par : 1 2 J= A 2R 25 où R est le vecteur de taille v constitué des coefficients de corrélation cyclique R Ck(m) pour lesdits 1 U-1 34 décalages temporels non nuls et lesdites fréquences cycliques, et A est une matrice proportionnelle à la matrice de corrélation de ces coefficients. 5. Méthode de détection selon la revendication 1, caractérisée en ce que la fonction de discrimination est calculée par : J = L L R `k ( `IZ mEIM îkEIK où R Ck(m), me lm sont les coefficients de corrélation cyclique, IM est ledit ensemble de décalages temporels non nuls et IK est ledit ensemble de fréquences cycliques. 6. Méthode de détection selon la revendication 4 ou 5, caractérisée en ce que les coefficients de corrélation cyclique sont estimés à partir des échantillons s(u) du signal à analyser par : 1 U-1 = ù Ls(u + m)s* (u\e2Z7L0Gku U u=o où îk est une fréquence cyclique appartenant à IK et m est un décalage temporel, exprimé en nombre de périodes 25 d'échantillonnage. 7. Méthode de détection d'un signal cyclostationnaire dans une pluralité Q de signaux à 15 20analyser (sq(t)), reçus d'un émetteur par une pluralité de récepteurs ou lus d'un support d'enregistrement par une pluralité de dispositifs de lecture, caractérisée en ce que pour chacun desdits signaux à analyser, on estime la valeur d'une fonction de discrimination Jq s'exprimant comme une forme quadratique des coefficients de corrélation cyclique de ce signal pour un ensemble (IMq) de décalages temporels non nuls et un ensemble (IKq) de fréquences cycliques, et que l'on (Q u compare la valeur 1-y Evq,E 4Jq à un taux de mauvaise 0=1 q=16q / détection Il pour déterminer si ledit signal cyclostationnaire est présent dans lesdits signaux à analyser, où y(v,x) est la fonction gamma inférieure incomplète normalisée, vq est le produit du cardinal (Mq) dudit ensemble de décalages temporels non nuls et du cardinal (Kq) dudit ensemble de fréquences cycliques pour le signal à analyser sq(t), Uq est la largeur de la fenêtre d'autocorrélation sur laquelle sont calculés les coefficients de corrélation cyclique et 6q est le carré de la variance du bruit, relatifs au signal à analyser. 8. Méthode de détection selon la revendication 7, caractérisée en ce que, pour chacun desdits signaux à analyser, ledit ensemble de fréquences cycliques se réduit à la fréquence nulle etque la valeur de la fonction de discrimination est calculée, pour ce signal, par : 2 Jq = R (m) mEIMq où R (m) est la valeur de la fonction d'autocorrélation du signal à analyser sq(t), pour un décalage temporel non nul me IMq , IMq étant ledit ensemble de décalages temporels non nuls pour ce signal. 10 9. Méthode de détection selon la revendication 8, caractérisée en ce que, pour chacun desdits signaux à analyser, ladite valeur de la fonction d'autocorrélation est estimée à partir des 15 échantillons s(u) de ce signal par : U-i k0 q(m)=ùLsq(u+m)s(u) U u=o où U est la largeur de la fenêtre d'autocorrélation et 20 m est un décalage temporel exprimé en nombre de périodes d'échantillonnage. 10. Méthode de détection selon la revendication 7, caractérisée en ce que, pour chaque 25 signal sq(t) à analyser, la fonction de discrimination est calculée par :5Jq = Aq 2Rq où Rq est le vecteur de taille vq constitué des coefficients de corrélation cyclique Rak(m) pour lesdits décalages temporels non nuls et lesdites fréquences cycliques, et Aq est une matrice proportionnelle à la matrice de corrélation de ces coefficients. 11. Méthode de détection selon la revendication 7, caractérisée en ce que, pour chaque signal sq(t) à analyser, la fonction de discrimination est calculée par : 37 2 Jq= L mEIMq CCkEIK, Rak (m) s, 2 où Rak(m), m EIMq sont les coefficients de corrélation cyclique du signal sq(t), IMq est ledit ensemble de décalages temporels non nuls et IKq est ledit ensemble de fréquences cycliques relatifs à ce signal. 12. Méthode de détection selon la revendication 10 ou 11, caractérisée en ce que, pour chaque signal sq(t) à analyser, les coefficients de corrélation cyclique sont estimés à partir des 25 échantillons sq(u) de ce signal par :20 38 U-1 s( m) _ ù s(u+m)s*(u)e'l- ku U u=0 où îk est une fréquence cyclique appartenant à IK et m est un décalage temporel, exprimé en nombre de périodes d'échantillonnage.