Brevet FR2669489A1 — publié le 19920522

Procédé et dispositif de reconnaissance de modulations La présente invention concerne un procédé et un dispositif de reconnaissance de modulations à partir de spectres instantanés. Pour reconnaître la classe de modulation de plusieurs émissions radioélectriques, il existe plusieurs techniques telles que celles décrites par exemple dans les articles suivants intitulés: "Automatic Classification of High Frequency Signals", Friedrich JONDIRAL, member EURASIP, AEG TELEFUNKEN, Signal processing 9 ( 1985) p 177-190; "Computer Simulation of an Automatic Classification Procedure for Digitally Modeled Communication Signals with Unknown Parameters", F.F LIEDTKE (FGAN), Signal processing 6 ( 1984) p 311-323; "Automatic Modulation Recognition Using Domain Parameters", Janet AISBETT, Electronics Research Laboratory, Defense Science and Technology Organization, Department of Defense Australia, Signal processing 13 ( 1987) p 323-328; "Sampling and Algorithms and Modulation Recognition", T G. CALLAGHAN WATKINS-JOHNSON, Microwaves and RF September 1985; "Modulation Clasification Based on Statistical Moments", Jackie E. HIPP, ph D, Southwest Research Institute, San Antonio, Texas, 1986 IEEE; "Receiver Recognition of Analogue Modulation Types", P M. FABRIZA, L B LOPES and G B LOCKHART. Celles-ci ne permettent cependant pas la reconnaissance d'émissions de façon simultanée lorsque la bande de fréquence surveillée est très étendue, car les signaux traités sont toujours issus de démodulateurs qui ont des bandes passantes limitées. D'autre part, l'acquisition précise des fréquences porteuses qui est effectuée généralement au moyen de synthétiseurs de fréquence prend un temps qui est préjudiciable à la rapidité des traitements. Egalement, suivant d'autres méthodes connues, la classification des émissions a lieu en mettant en oeuvre des algorithmes de maximum de vraisemblance appliqués sur un nombre déterminé de paramètres, ces paramètres résultant d'une synthèse plus ou moins parfaite de l'information au niveau de l'émission des informations portées par chaque canal Cette synthèse qui a pour but de réduire la redondance et l'interdépendance des paramètres est en général nécessaire pour permettre l'application de méthodes de classification particulières comme celle par exemple connue sous le nom de méthode de classification de BAYES Cependant, les algorithmes mis en oeuvre et nécessaires pour réduire la redondance et l'interdépendance des paramètres nécessite de longues opérations de traitement qui prennent du temps de calcul et qui sont préjudiciables à la rapidité des détections mises en oeuvre. Le but de l'invention est de pallier les inconvénients pré-cités. Le procédé de reconnaissance de la modulation d'émissions radio-électriques à partir de spectres instantanés d'émission relevés dans une bande de fréquence déterminée par un analyseur -de spectre à transformée de Fourier rapide du type consistant à calculer les paramètres suivants pour chaque raie de spectre d'émission relevé dans la bande de fréquence déterminée, une amplitude moyenne de toutes les raies k des spectres contenus dans la bande de fréquence déterminée, un rapport signal à bruit RS Bk, un écart type E Tk en amplitude de chaque raie du spectre et un coefficient de corrélation COR(k,k') d'amplitude de chaque raie k avec les raies homologues des spectres des émissions contenues dans la bande de fréquence déterminée caractérisé en ce qu'il consiste à comparer au travers d'un réseau neurones les paramètres de chaque spectre d'émission à des paramètres d'émission attendue, et à déclarer qu'une émission correspond à une émission attendue lorsque l'écart détecté par la comparaison entre les paramètres de l'émission et ceux attendus est minimum. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard des dessins annexés qui représentent: la figure 1 un mode de réalisation d'un dispositif selon l'invention; la figure 2 un mode d'organisation du réseau de neurones de la figure 1; la figure 3 une représentation schématique d'un neurone la figure 4 une courbe de réponse d'un neurone. Le dispositif représenté à la figure 1 comprend un récepteur 1, un analyseur de spectre 2, un bloc d'extraction 3, représenté à l'intérieur d'une ligne fermée en pointillés, un réseau de neurones 4 et un module d'apprentissage 5 Le bloc d'extraction 3 comprend un bloc de calcul de paramètres primaires 6 représenté à l'intérieur d'une ligne fermée en pointillés et un bloc de détection 7. Le bloc de calcul des paramètres primaires 6 est formé d'un bloc de calcul intermédiaire 8 et d'un bloc de filtres 9. Le récepteur 1, est un récepteur radio à large bande, calé sur la fréquence centrale de la bande de fréquence à analyser Le signal capté par le récepteur est filtré de façon connue en sortie des étages de fréquence intermédiaire de celui-ci Ce récepteur est réalisé de façon connue pour capter des signaux radio fréquence transmis dans la bande HF normalisée par la commission des communications fédérale des Etats-Unis connue encore sous l'abréviation FCC, ladite bande normalisée s'étendant dans la gamme 3 à 30 M Hz correspondant à des longueurs d'onde de 100 à 10 mètres c'est-à-dire des ondes décamétriques Ces signaux peuvent être modulés suivant des modulations normalisées du type AO, Ai, A 3, A 3 H+, A 3 H-, A 3 J+, A 3 J-, F 3, Fl L'analyseur de spectre 2 comprend de façon connue un convertisseur analogique numérique et un bloc de calcul non représentés Le convertisseur analogique fournit 2 N point résultant de l'échantillonnage à la vitesse 2 B de la bande du signal à analyser. Le signal fourni par le convertisseur analogique numérique est multiplié également de façon connue par une fenêtre de pondération du type "Blackman-Harris à quatre termes". Le bloc de calcul intermédiaire 8 effectue le calcul d'une raie sur deux du spectre du signal en appliquant aux échantillons du signal fourni par la fenêtre de pondération une transformation de FOURIER rapide encore connue sous l'abréviation FFT, suivant les méthodes décrites par exemple dans le livre de E ORAN BRIGHAM ayant pour titre "the Fast Fourier Transform" Printice Hall Inc. Englewood Cliffs New Jersey 1974, ou dans le livre de P M. BEAUFILS et RAMI ayant pour titre "le filtrage numérique" p 117 à édité par SYBEX en 1986. La résolution obtenue est égale à deux B/N. Les résultats de calcul de la FFT sont fournis suivant une période Delta T et une résolution Delta F fixe (exemples Delta T = 2 à 8 ms, Delta F = 1 K Hz à 250 Hz). Le bloc d'extraction est composé des blocs de calcul de paramètre 6 et de bloc de détection 7 pour permettre de traiter l'ensemble des émissions de la bande B en temps réel Le bloc de calcul intermédiaire 8 anticipe le plus possible le calcul des paramètres primaires par des calculs répétitifs identiques pour tous les canaux et indépendants des émissions présentes dans la bande B du récepteur 1. Le bloc de détection 7 calcule à partir des résultats des calculs intermédiaires, la position dans la bande B des différents émetteurs, (début, fin) et il estime le niveau de bruit NVB. Le bloc de filtres 9 calcule à partir des informations de calcul intermédiaire et de position des émetteurs dans la bande B, les vecteurs de forme de chaque émission constituée par les valeurs des paramètres primaires. Le bloc d'extraction 3 exécute ces tâches de façon périodique, la période étant un multiple M de Delta T En fixant par exemple M à 300 spectres et Delta T à 4 ms, la période d'extraction est alors de 1,2 secondes. Pendant chaque période d'extraction, le bloc d'extraction 3 calcule pour chaque canal k de la bande B le rapport signal à bruit RSB, l'écart type en amplitude E Tk et un coefficient de corrélation COR k, k' inter-raie du canal k avec un nombre déterminé de canaux fixés, par exemple égal à 9 tel que k' = k + 1 à k + 9 qui fixe également la largeur de bande LB de chaque canal. Ces calculs sont effectués conjointement par les blocs de calcul intermédiaires 8, de détection 7 et le filtre 9 qui sont formés de façon connue par des processeurs de traitement du signal convenablement programmés. Ces processeurs sont programmés pour effectuer les opérations de calcul de rapport signal à bruit RS Bk de chaque canal k, les calculs d'écart type E Tk sur I'amplitude de chaque canal k et des calculs de coefficient de corrélation COR(k,k') entre les variations d'amplitude des canaux k et k'. Le calcul du rapport signal à bruit RS Bk pour un canal k a lieu en considérant l'amplitude A Mk de la raie k de chaque spectre N et en calculant l'amplitude moyenne Ank de la raie k durant la période d'extraction suivant les relations. M A Mk = 1 Amk en d Bm ( 1) M m= 1 et RS Bk = A Mk NVB en d Bm ( 2) o NVB correspond au niveau moyen de bruit. Les calculs de l'écart type E Tk des paramètres COR(k,k') ont lieu en application des relations M E Tk = l (Amk A Mk) ( 3) m=m 1 M Z (Am,k A Mk)(Amk' A Mk') COR(k,k') m=l M 4) (Amk A Mk)2,_ (Am' A Mk)22 rnml m=l Le bloc de filtre 9 effectue ensuite la sélection de Ni valeur de rapport signal à bruit RS Bk, de N 2 valeurs d'écart type E Tk et de N 3 valeurs des coefficients de corrélation COR(k,k'). Les canaux sélectionnés sont pris arbitrairement à partir du centre de l'émission Si la largeur de bande de l'émission est trop faible les valeurs manquantes sont complétées par des quantités moyennes A titre d'exemple, un choix possible peut être déterminé comme suit avec Ni = 7, N 2 = 5 et N 3 = 10. Un mode de réalisation du réseau de neurones 4 est représenté à la figure 2 Ce réseau comporte deux sous-réseaux 10 et 11 Le sous-réseau 10 comprend également NC 1 neurones cachés ainsi que 9 neurones de sortie, un par type de modulation Le sous-réseau 11 est couplé au premier sous-réseau 10 et comprend 9 entrées correspondant aux 9 sorties du sous-réseau 10, NC 2 neurones cachés et une sortie correspondant à l'information de validation Dans l'architecture représentée à la figure 2 chaque neurone utilisé du type connu des neurones sygmoides dont un schéma de principe est représenté à la figure 3 Suivant ce type d'organisation l'état électrique de la sortie d'un neurone est défini en fonction de l'état de ces entrées par une relation de la forme O ( 5) N. 1 + Exp ( L Omega X; Tétai) j=i ou O désigne la sortie du neurone i, Teta désigne le seuil associé au neurone i, X désigne l'entrée j du neurone i, -N est le nombre d'entrée du neurone i, et Omega i désigne le poids attribué à l'entrée j. A partir des paramètres issus du bloc de filtre 9 le réseau calcule les sorties de tous les neurones par application de la relation ( 5). Le type de modulation trouvé est celui pour lequel le neurone associé présente un niveau de signal de sortie maximal A partir des sorties du sous-réseau 10, le sous-réseau 11 calcule une information de validation Si cette information de validation est supérieure à un seuil S la reconnaissance est validée sinon il y a rejet L'apprentissage effectué par le bloc d'apprentissage 5 est de type de celui connu sous la désignation "rétro-propagation"' qui est décrite dans l'article parallèle "Distributed Processing" de Rummel Hartington WILLIAMS de 1986 Cet apprentissage consiste pour chaque entrée présentée, à comparer la sortie du réseau avec la sortie attendue La différence constitue un vecteur d'erreur utilisé pour corriger les poids Omegaji des neurones de sortie. L'erreur est ensuite rétro-propagée suivant un algorithme de type gradiant Les corrections s'effectuent alors de la sortie vers l'entrée du réseau L'apprentissage a lieu d'abord en corrigeant les coefficients Omegaji du sous-réseau 10 puis ceux du sous-réseau 11. REVENDICATIONS 1 Procédé de reconnaissance de la modulation d'émissions radioélectriques à partir de spectres instantanés d'émissions relevés dans une bande de fréquence déterminée par un analyseur de spectre ( 2) à transformée de Fourier rapide du type consistant à calculer ( 3) les paramètres suivants pour chaque raie de spectre d'émission relevé dans la bande de fréquence déterminée, une amplitude moyenne de toutes les raies k des spectres contenus dans la bande de fréquence déterminée, d'un rapport signal à bruit RS Bk, un écart type E Tk en amplitude de chaque raie du spectre et un coefficient de corrélation COR(k, k') d'amplitude de chaque raie k avec les raies homologues des spectres des émissions contenues dans la bande de fréquence déterminée caractérisé en ce qu'il consiste à comparer au travers d'un réseau neurones les paramètres de chaque spectre d'émission à des paramètres d'émission attendue, et à déclarer qu'une émission correspond à une émission attendue lorsque l'écart détecté par la comparaison entre les paramètres de l'émission et ceux attendus est minimum. 2 Procédé selon la revendication 1 caractérisé en ce qu'il consiste à associer dans le réseau de neurones, un neurone à chaque type de modulation. 3 Procédé selon la revendication 2 caractérisé en ce qu'il consiste à déclarer qu'un type de modulation est trouvé lorsque le niveau de sortie du neurone associé est maximum. 4 Procédé selon l'une quelconque des revendications 2 et 3 caractérisé en ce que le niveau de sortie Qi d'un neurone i est déterminé par une relation de la forme Oi = Ni 1 + Exp(-Z Omega i X Teta) j=l o: Teta désigne un seuil associé au neurone i X désigne une entrée j Ni est le nombre d'entrées du neurone i et Omega 1 i désigne le poids attribué à l'entrée j. Procédé selon l'une quelconque des revendications 2 et 3 caractérisé en ce que le niveau de sortie Oi d'un neurone est déterminé par une relation de la forme Ni O = SIGN ( 5 Omega X; Teta) j=-1 o: Tetai désigne un seuil associé au neurone i X désigne une entrée j Ni est le nombre d'entrées du neurone i et Omegaii désigne le poids attribué à l'entrée j. 6 Procédé selon l'une quelconque des revendications 1 à 4 caractérisé en ce qu'il consiste à ajuster dans une phase d'apprentissage les poids Omega ii de chaque neurone correspondant à un type d'émission attendu en comparant pour chaque émission attendue la sortie du réseau avec la sortie attendue en rétropageant l'erreur obtenue suivant un algorithme de type gradiant qui minimise cette erreur. 7 Dispositif pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 6 caractérisé en ce qu'il comprend un récepteur ( 1) couplé à un réseau de neurones ( 4) par l'intermédiaire d'un analyseur de spectres ( 2) et d'un dispositif d'extraction de paramètres d'émission ( 3). 8 Dispositif selon la revendication 7 caractérisé en ce que le dispositif d'extraction comprend un bloc de filtres ( 9) de détermination de paramètre primaire d'émissions couplé à un bloc de calcul intermédiaire ( 8) et à un bloc de détection d'émission ( 7). 9 Dispositif selon la revendication 8 caractérisé en ce que le bloc de filtres ( 9), le bloc de calcul intermédiaire ( 8) et le bloc de détection d'émission ( 7), sont fournis par des processeurs de traitement du signal. Dispositif selon la revendication 9 caractérisé en ce qui -les processus de traitement du signal sont programmés pour calculer pour chaque raie k de spectre d'émission, une amplitude moyenne des raies de rang k des spectres contenus dans la bande de fréquence, un rapport signal à bruit RS Bk, un écart type E Tk en amplitude et un coefficient de corrélation COR(k, k') d'amplitude avec les raies homologues dans les spectres d'émission contenues dans la bande de fréquence.