FR2678738A1 - Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions. - Google Patents

Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions. Download PDF

Info

Publication number
FR2678738A1
FR2678738A1 FR8603234A FR8603234A FR2678738A1 FR 2678738 A1 FR2678738 A1 FR 2678738A1 FR 8603234 A FR8603234 A FR 8603234A FR 8603234 A FR8603234 A FR 8603234A FR 2678738 A1 FR2678738 A1 FR 2678738A1
Authority
FR
France
Prior art keywords
cell
signal
noise
noise reference
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR8603234A
Other languages
English (en)
Other versions
FR2678738B1 (fr
Inventor
Siefert Eric Pierre
Zolesio Jean-Luc Emmanuel
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Laboratoire Central de Telecommunications SA
Original Assignee
Laboratoire Central de Telecommunications SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Laboratoire Central de Telecommunications SA filed Critical Laboratoire Central de Telecommunications SA
Priority to FR8603234A priority Critical patent/FR2678738B1/fr
Priority to DE19873707594 priority patent/DE3707594A1/de
Publication of FR2678738A1 publication Critical patent/FR2678738A1/fr
Application granted granted Critical
Publication of FR2678738B1 publication Critical patent/FR2678738B1/fr
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • G01S7/2921Extracting wanted echo-signals based on data belonging to one radar period
    • G01S7/2922Extracting wanted echo-signals based on data belonging to one radar period by using a controlled threshold
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
    • G01S13/534Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi based upon amplitude or phase shift resulting from movement of objects, with reference to the surrounding clutter echo signal, e.g. non coherent MTi, clutter referenced MTi, externally coherent MTi

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

L'invention concerne un procédé de détection autogressive d'un signal sinusoidal complexe dans du bruit et d'estimation de sa fréquence et un complexe dans du bruit et d'estimation de sa fréquence et un dispositif en faisant application, notamment pour la détection d'une cible dans un signal de radar à impulsions. Après filtrage préalable (CCr1, CCrp, FT1, FTp) pour obtenir un signal utile dans du bruit blanc, un calculateur de coefficient (Cr1) détermine le coefficient de réflexion d'ordre 1 dans l'analyse autorégressive du signal filtré ; ce coefficient de réflexion est envoyé à un circuit de détection (CD) qui le compare à une référence de bruit calculée pour en déduire l'information X(i), de la présence ou l'absence de cible dans la cellule de résolution i. Le coefficient de réflexion est aussi envoyé à un circuit de calcul d'argument (CCA) qui extrait l'argument dudit coefficient complexe et en déduit la fréquence du signal utile. L'invention s'applique en particulier à la détection de cibles radar.

Description

La présente invention se rapporte à un procédé de détection autorégressive d'un signal sinusoïdal complexe dans du bruit et d'estimation de sa fréquence, notamment pour la détection d'une cible monofrêquentielle dans un signal de radar a impulsions.
Elle se rapporte également à un dispositif pour la mise en oeuvre de ce procédé.
Dans le domaine du radar, il est connu d'effectuer une détection des cibles par une intégration, cohérente ou non, des échos reçus.
L'intégration cohérente, qui est la méthode présentant la plus grande sensibilité de détection, s 'effectue généralement par l'intermédiaire d'un calcul de transformée de Fourier discrete qui, en théorie, est le traitement optimal si la fréquence du signal analysé correspond a une des fréquences d'analyse de la transformée de Fourier. Par contre si ce n'est plus le cas, c' est-a-dire quand la fréquence du signal est comprise entre deux fréquences d'analyse de la transformée de Fourier discrete, il y a une perte de sensibilité à la détection. Cette perte est encore accentuée par le fait qu'on utilise une pondération du signal temporel.
Cette technique de traitement est une technique appliquée dans le domaine temporel.
Par ailleurs, dans l'analyse spectrale, on applique depuis un certain nombre d'années des méthodes d'analyse autorégressive. En particulier, ces méthodes ont été utilisées pour éliminer de façon adaptative le fouillis dans les signaux radar.
Cependant, ces méthodes étaient ainsi toujours utilisées comme un outil de traitement appliqué dans le domaine temporel.
Un objet de l'invention est donc de définir un nouveau procédé de détection ayant de meilleures performances que les techniques connues et reposant sur la constatation que l'analyse autorégressive d'un signal radar peut être considérée comme définissant un nouvel espace de représentation des réponses radar. En effet, dans le cas d'une analyse autorégressive du signal radar, le premier coefficient de réflexion contient l'information relative à la présence d'une cible.
Un autre objet de l'invention est donc un procédé de détection mettant en oeuvre des criteres de caractérisation autorégressive des cibles dans un signal radar, ou plus généralement, d'un signal sinusoïdal complexe dans du bruit.
Selon l'invention, il est prévu un procédé de détection autorégressive d'un signal sinusoidal complexe dans du bruit et d'estimation de sa fréquence, caractérisé en ce qu'il consiste à - effectuer une analyse autorégressive du signal pour obtenir le coeffi
cient de réflexion d'ordre 1 - déterminer une référence de bruit fixe ou calculée; - comparer les valeurs dudit coefficient de réflexion à ladite référence
de bruit pour détecter la présence dudit signal sinusoidal, et - calculer l'argument dudit coefficient de réflexion pour en déduire la
fréquence dudit signal sinusoïdal.
Un tel procédé, meme s'il a été présenté dans le cadre du problème de la détection de cibles dans un signal radar, peut etre appliqué dans d'autres domaines (sonar, etc...) et concerner aussi bien l'analyse de signaux impulsionnels que continus, numériques ou analogiques.
Par rapport à une détection via pondération et transformée de
Fourier discrète, le procédé selon l'invention présente une meilleure sensibilité et permet une estimation de la fréquence du signal utile plus précise et non biaisée et une mise en oeuvre simple.
Par ailleurs, la détection peut se faire soit par rapport à une référence de bruit fixe, soit par rapport à une référence de bruit calculée ce qui conduit alors à une détection à taux de fausse alarme constant.
L'invention sera mieux comprise et d'autres caractéristiques et avantages apparaîtront à l'aide de la description ci-après et des dessins joints où - la figure 1 est un diagramme spectral explicatif - la figure 2 représente un schéma de principe d'un dispositif selon
l'invention dans le cas d'un radar à impulsions ; - la figure 3 est le schéma d'un mode de réalisation de détecteur auto
régressif selon l'invention ; - la figure 4 est un diagramme explicatif - la figure 5 est un schéma d'un circuit de calcul de référence de bruit
pour le détecteur de la figure 3, et - la figure 6 représente le schéma d'une variante de circuit de calcul de
référence de bruit.
Bien que cela ne doive pas etre considéré comme limitatif de l'invention, on se placera, pour une plus grande simplicité de l'explication, dans le cas du traitement d'un signal radar contenant des échos de cibles utiles dans du bruit blanc ou sensiblement blanc. On considèrera en outre qu'on est dans le cas d'un radar à impulsions fournissant un signal complexe sous forme d'échantillons numériques successifs.
il faut noter que l'on peut toujours se ramener à l'hypothèse envisagée d'un signal utile dans du bruit blanc en effectuant un filtrage préalable du signal radar. Une méthode particulièrement efficace a été développée par la demanderesse pour l'élimination adaptative du fouillis par modélisation autorégressive et est décrite en particulier dans les demandes de brevet français nO 83 14673 déposée le 15 septembre 1983, intitulée "Procédé et dispositif pour la détermination adaptative des coefficients d'un filtre transversal d'ordre n permettant d'éliminer les signaux de fouillis dans un radar Doppler à impulsions", et nO 84 04722 déposée le 27 mars 1984, intitulée "Procédé et dispositif d'élimination adaptative de fouillis dans un radar Doppler à impulsions".
Comme on l'a déjà mentionné plus haut, en effectuant une analyse autorégressive du signal, l'information relative à la cible est contenue dans le coefficient de réflexion d'ordre 1. Ce coefficient est, on le sait, donné par la relation
Figure img00030001

où xn est la valeur de l'échantillon complexe de rang n reçu d'une cellule de résolution i, n variant de 1 à N si N est le nombre d'échantillons reçus d'une cellule et correspondant aux N impulsions radar émises pendant le passage du faisceau radar sur cette cellule.
Cependant, la simple utilisation de ce coefficient de réflexion d'ordre 1 ne permet pas une détection optimale comme on va le voir cidessous. En effet, soient deux types de réalisations différentes notées Ci et V correspondant respectivement à des réalisations de bruit blanc et à des réalisations de signal sinusoïdal dans du bruit blanc.
il est clair qu'aucune détection n'a lieu lorsque les coefficients de réflexion estimés relatifs à deux réalisations r1( et rl(y) ont le meme module t
Figure img00030002
Cependant, comme on peut considérer que le coefficient rl est R(-1) pratiquement égal au rapport R(0) ' où R(x) désigne la fonction d'autocorrélation au point x du signal reçu de la cellule en cours d'examen, il est évident qu'en moyenne, la fonction d'autocorrélation R' > (0) pour la réalisation # est supérieure à R#(0). On a donc :
Figure img00040001

où ss a une valeur supérieure à 0 et où r1 désigne l'espérance mathématique de r1.
D'autre part, la figure 1 illustre un autre inconvénient que l'on rencontre lorsque l'on utilise simplement le module du coefficient r1 pour la détection. Sur cette figure, on a représenté le spectre de deux réalisations # différentes mais pour lesquelles la différence A entre pic et plancher du spectre est la même. Dans le cas d'une détection par transformation de Fourier discrète qui définit un seuil de détection S fixe, seule la réalisation V2 donne lieu à une détection, conséquence de sa puissance par rapport à une référence fixe. Au contraire, dans le cas de la détection par le coefficient de réflexion r1, les deux réalisations donneront lieu à la même conclusion puisque le coefficient de réflexion ne tient compte que de la différence A.
Une modification de principe apportée selon l'invention pour remédier à ces inconvénients a donc été d'introduire dans la détection autorégressive une référence fixe de puissance. Selon une caractéristique de l'invention, on considère donc un coefficient r1 modifié tel que
R(-1)
r1 =
Max (R(0),Kb) où K est un coefficient prédéterminé supérieur à 1.
Avec ce nouveau coefficient, on a, pour la réalisation # :
Figure img00040002

puisque R(O) est égale à b, puissance moyenne du bruit.
Pour la réalisation 9, on a suivant le cas 1) Si b+ss # Kb
Figure img00040003

2) Si b+ss > Kb
Figure img00040004
Donc, dans les deux cas, il y aura détection de la réalisation # car
Figure img00050001

alors que cela n'était pas possible avec le coefficient de réflexion r non modifié.
Comme on le voit, le cas le plus defavorable pour le coefficient modifié est le premier cas où b+ss # Kb. On doit donc se placer dans ce cas pour choisir un coefficient K optimum qui est fonction, par ailleurs, de l'écart type du bruit radar, de la probabilité de fausse alarme (Pfa) imposée et du nombre d'échantillons N. Cette optimisation est obtenue en prenant pour K la valeur qui donne, dans ce cas aussi
Figure img00050002

c'est-à-dire:
K=l+f3
b
La valeur de K optimum peut etre prédéterminée en fonction des paramètres choisis (Pfa, N, ...).
Le coefficient de détection, dérivé du coefficient de réflexion, qui a été défini ci-dessus sera désigné par ^ ,ARC pour le distinguer d'une variante de coefficient de détection selon l'invention qui sera définie ci-dessous. Pour ce coefficient de réflexion 21 ARC on a pris comme référence de bruit
RARC = Max (R(0),Kb)
On peut définir un autre coefficient de détection r1,ARM, non normé, en prenant, comme référence de bruit RARM,la puissance moyenne de bruit.On a
f - R(-1)
1,ARM b
On voit que l'on a alors dans les memes conditions que précédemment
Figure img00050003

et, dans tous les cas
Figure img00050004
Là aussi, on obtient une détection de la réalisation # puisque :
Figure img00050005
A partir de ces principes ainsi exposés, un dispositif de détection autorégressive selon l'invention est représenté sur la figure 2.
Le signal radar sous forme d'échantillons xnm (échantillon complexe numéro n reçu de la cellule de résolution numéro m) est d'abord appliqué à un dispositif pour l'élimination adaptative de fouillis du type décrit plus en détail dans la demande de brevet nO 84 04722 citée ci-dessus.
Ce dispositif comporte des cellules de filtre en treillis d'ordre 1 à p, FT1 à FTp, et des calculateurs de coefficient de réflexion associés CCr1 à
CCrp fournissant la valeur moyenne du coefficient de réflexion R1 à Rp sur la zone perturbée considérée.
Chaque cellule de filtre en treillis, telle la cellule FT1, comprend de manière connue un registre 10 introduisant un retard R, égal à une période d'échantillonnage, dans la voie retardée, des mémoires 11, 12 conservant en mémoire les signaux d'entrée pendant le temps kR, nécessaire au calcul du coefficient de réflexion, et des multiplieurs 13, 14 et des additionneurs-soustracteurs 15, 16 exécutant les calculs nécessaires à partir du coefficient de réflexion et de son conjugué désigné par le signe *.
Le signal filtré obtenu sur la voie directe ou prédictive en sortie de la cellule FTp comprend le signal utile de cible dans du bruit blanc. Ce signal est envoyé d'une part à un calculateur de coefficient de réflexion d'ordre 1, Cri, et d'autre part à un circuit de détection CD qui reçoit ce coefficient et une valeur de seuil aJustable et qui fournit le cas échéant un signal de présence de cible X(i) dans la cellule de numéro i après avoir effectué le calcul de la référence de bruit.
Enfin, un circuit de calcul d'argument CCA extrait l'argument du coefficient de réflexion qui permet d'obtenir la fréquence FD(i) de la cible éventuellement présente dans la cellule de résolution i.
La figure 3 représente le schéma d'un mode de réalisation du dispositif de détection autorégressive selon l'invention. Comme on l'a vu dans l'exposé des principes de l'invention, on utilise aussi bien pour le coefficient de détection rl,ARC que pour le coefficient le numérateur du coefficient de réflexion r1, sensiblement égal à R(-1).
Le dispositif de la figure 3 comprend donc des moyens de calcul pour, à partir des N échantillons d'une cellule de résolution i, calculer l'expression
Figure img00060001
Ces moyens de calcul reçoivent les parties réelle w (xmn) et imaginaire f (xn) correspondant aux N échantillons successifs reçus de la cellule m. Ils comprennent deux registres 20, 21 permettant de retarder ces parties réelle et imaginaire d'un retard R, égal à une période d'échantillonnage, quatre multiplieurs accumulateurs 22 à 25, un additionneur 26 et un soustracteur 27 fournissant respectivement la partie réelle et la partie imaginaire de l'expression
Figure img00070001
Ces valeurs sont envoyées à deux mémoires 28, 29 qui les retardent d'un temps k'R.
Le traitement effectué par le dispositif de la figure 3 sera mieux compris en se reportant au diagramne de la figure 4. Sur ce diagramme est représenté le panorama radar. Une cellule de résolution est délimitée par les lignes de gisement espacées entre elles d'une largeur du faisceau radar et par les cercles de distance espacés entre eux de la valeur de la résolution en distance du radar. La référence de bruit nécessaire pour la détection est calculée sur une zone limitée Z entourant la cellule Ci en cours de traitement. Dans l'exemple de réalisation choisi, on s' est limité, dans un but de simplification, à une zone s'étendant sur une seule largeur de faisceau et comprenant
M cellules C1 à CM, la cellule Ci se trouvant pratiquement au centre de la zone. Mais on pourrait, bien sûr, prendre une zone s'étalant sur plusieurs largeurs de faisceau de part et d'autre.
Ceci étant, il est clair que le temps de stockage k > R dans les mémoires 28 et 29 est tel que, lorsque l'on inscrit les valeurs corres pondant à la cellule m, on lit les valeurs correspondant à la
M M-i cellule i = m - 2 ou i = m - 2 selon que le nombre M de cellules de la zone Z est pair ou impair. On obtient donc aux sorties des mémoires 28 et 29 les parties réelle et imaginaire de l'expression
Figure img00070002

qui constitue le numérateur du coefficient de détection pour la cellule i en cours de traitement. Ces parties réelle et imaginaire sont envoyées à un calculateur de module 30.
D'autre part, les parties réelle et imaginaire des échantillons
m reçus xn sont envoyéesàun ensemble de calcul de la fonction d'autocorrélation au point zéro du signal émanant de la cellule m. Cet ensemble de calcul comprend des multiplieurs accumulateurs 31 et 32 et un additionneur 33 et fournit la valeur
Figure img00080001
Cette valeur est envoyée à un circuit de calcul de référence de bruit 34 qui calcule soit la référence RARC soit la référence RARE.
Cette référence de bruit est multipliée, dans un multiplieur 35, par une valeur de seuil déterminée, en particulier, en fonction de la probabilité de fausse alarme souhaitée. La sortie du multiplieur 35 est envoyée à un comparateur 36 qui reçoit d'autre part le module fourni par le circuit 30. Ce comparateur fournit un signal de présence de cible X(i) dans la cellule i si
Figure img00080002

R ou ARM Seuil soit
II 1,ARC ou ARMAI > Seuil
Enfin, il est prévu des circuits d'extraction de l'argument du coefficient de réflexion, ctest- -dire ici de l'argument de la quantité complexe
Figure img00080003

qui est égal à la pulsation #D(i) du signal utile de cible.
Pour cela, on utilise une table contenue dans une mémoire morte 40 du type PROM. Cette mémoire est adressée par les bits de plus fort poids des parties réelle et imaginaire fournies par les mémoires 28 et 29. Si on utilise une mémoire PROM courante à treize bits d'adresse, on peut par exemple l'adresser par les six bits de plus fort poids et le bit de signe de la partie imaginaire et par les six bits de plus fort poids de la partie réelle.On obtient ainsi, en mettant de côté le bit de signe de la partie réelle, un angle compris entre - 2 et + 2 qui vérifie avec l'argument vrai r les relations indiquées dans le tableau ci-après.
Figure img00090001
<tb> SIGNE <SEP> PARTIE <SEP> REELLE <SEP> SIGNE <SEP> PARTIE <SEP> IMAGINAIRE <SEP> VALEUR <SEP> DE <SEP> <SEP>
<tb> <SEP> + <SEP> +
<tb> <SEP> + <SEP> n-p <SEP>
<tb> <SEP> + <SEP> - <SEP> #
<tb>
- - # - #
On corrige donc la valeur #, fournie par la mémoire morte 40, à l'aide d'un élément de mémoire 41 et d'un additionneur ou soustracteur 42 commandés par le bit de signe de la partie réelle venant de . la mémoire 28. Si ce signe est positif, l'élément 41 fournit la valeur O et le circuit 42 est additionneur. Si ce signe est négatif, l'élément 41 fournit la valeur Tr et le circuit 42 est soustracteur.
La figure 5 représente le schéma détaillé du circuit 34 de calcul de référence de bruit lorsque l'on utilise la référence RARE.
Ce circuit de calcul de RARM qu'on désigne par CRB, reçoit la valeur
Figure img00090002

correspondant à la cellule de résolution m dont les
N échantillons ont été déjà regroupés par un redistributeur utilisé pour le filtrage adaptatif, comme décrit dans la demande de brevet nO 84 04722 déjà mentionnée. Cette valeur est envoyée à une première mémoire 340 qui les stocke et les restitue avec un retard k'R identique à celui des mémoires 28 et 29 de la figure 3. En sortie de la mémoire 344 on a donc la valeur
Figure img00090003

correspondant à la cellule i en cours de traitement. Il s'agit pour le circuit CRB de fournir la référence de bruit RARM pour la cellule i :
Figure img00090004

où M est, comme on l'a vu, le nombre de cellules de la zone Z (figure 4) entourant la cellule i.
Le circuit CRB assure une remise à jour continue de la référence RARM. Pour cela, lors de la réception des signaux de la cellule m, il faut ajouter les contributions de cette cellule et de la cellule i-l précédemment traitée, qui n'avait pas été prise en compte dans le calcul de la référence, et retrancher les contributions de la plus ancienne cellule de la zone, c'est-a-dire ici la cellule m-M, et de la cellule i en cours de traitement qui ne doit pas intervenir dans le calcul de la référence de bruit. Ceci est réalisé à l'aide de la mémoire 341 introduisant un retard k'R, du registre 342 introduisant un retard R, des additionneurs 343 à 345 et du soustracteur accumulateur 346. Un multiplieur 347 assure ensuite la multiplication par la constante M 1 .
La figure 6 représente le schéma du circuit 34 de calcul de référence de bruit RARC. Ce circuit doit calculer la référence de bruit selon la relation
Figure img00100001
On utilise pour cela le circuit CRB de la figure 5 qui fournit d'une part le second terme
Figure img00100002

de la comparaison, prélevé à la sortie de la mémoire 340, et d'autre part le premier terme de la comparaison, après multiplication de RARM par le coefficient K prédéterminé à l'aide du multiplieur 348. Ces deux termes sont envoyés à un circuit de comparaison et de sélection 349 qui fournit la valeur maximum comme référence de bruit RARC.
Comme on peut le constater, la charge de calcul, et donc le nombre de circuits nécessaires, n' est pas très élevée. Elle est moins importante que celle exigée pour un traitement par transformée de Fourier discrète. Ainsi qu'on l'a déjà mentionné, la sensibilité de détection est nettement améliorée. D'autre part, il est clair qu'on peut très facilement reconfigurer le système lorsque l'on modifie par exemple le nombre d'échantillons des réponses radar.
Bien entendu, l'exemple de réalisation décrit n'est nullement limitatif de l'invention. En particulier, celle-ci n'est pas limitée au domaine des signaux radar et peut s'appliquer à la détection de tout signal sinusoidal dans du bruit.

Claims (11)

  1. REVENDICATIONS
    fréquence dudit signal sinusoidal.
    de bruit pour détecter la présence dudit signal sinusoïdal, et - calculer l'argument dudit coefficient de réflexion pour en déduire la
    cient de réflexion d'ordre 1 - déterminer une référence de bruit fixe ou calculée - comparer les valeurs dudit coefficient de réflexion à ladite référence
    I. Procédé de détection autorégressive d'un signal sinusoïdal complexe dans du bruit et d'estimation de sa fréquence, caractérisé en ce qu'il consiste à : - effectuer une analyse autorégressive du signal pour obtenir le coeffi
  2. 2. Procédé selon la revendication 1, caractérisé en ce qu'avant d'effectuer ladite analyse autorégressive, on filtre le signal pour obtenir un bruit sensiblement blanc.
  3. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que, pour déterminer ladite référence de bruit, on calcule la puissance moyenne de bruit dans une zone temporelle entourant la partie de signal examinée.
  4. 4. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que, pour déterminer ladite référence de bruît, on calcule d'une part le produit par un coefficient K prédéterminé, supérieur à 1, de la puissance moyenne de bruit dans une zone temporelle entourant la partie de signal examinée et, d'autre part, la valeur de la fonction d'autocorrélation au point zéro du signal dans ladite partie examinée et on prend le maximum de ces deux valeurs calculées.
  5. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que, pour comparer les valeurs du coefficient de réflexion à ladite référence de bruit, on calcule, à partir de ces valeurs et de cette référence de bruit, un coefficient de détection égal au rapport entre le produit du coefficient de réflexion par la valeur de la fonction d'autocorrélation au point zéro du signal dans ladite partie examinée d'une part et ladite référence de bruit d'autre part, et en ce que l'on compare ledit coefficient de détection à une valeur de seuil prédéterminée pour en déduire s'il existe ou non un signal sinusoidal dans ladite partie de signal examinée.
  6. 6. Procédé selon la revendication 5, pour détecter à partir d'un signal radar échantillonné complexe la présence d'une cible monofréquentielle dans un environnement de bruit blanc et estimer sa fréquence, ledit signal radar étant constitué pour chaque cellule de résolution de
    N échantillons successifs correspondant aux N impulsions radar émises pendant le passage du faisceau radar sur ladite cellule, caractérisé en ce que ledit produit du coefficient de réflexion par la valeur de ladite fonction d'autocorrélation au point zéro pour une cellule i donnée est obtenu à partir des valeurs des N échantillons complexes xi, reçus de ladite cellule i, (où n est le numéro d'ordre de l'échantillon et varie entre 1 et N) selon l'expression
    Figure img00120001
    et en ce que la comparaison du coefficient de détection de ladite cellule au seuil prédéterminé est réalisée en comparant ledit produit au produit de la référence de bruit par ledit seuil.
  7. 7. Procédé selon la revendication 6, caractérisé en ce que ladite référence de bruit est déterminée en faisant la moyenne, sur M-l cellules entourant ladite cellule i, des valeurs de la fonction d'autocorrélation au point zéro du signal reçu de chacune de ces
    M-1 cellules selon la relation
    Figure img00120002
    où RARM représente la référence de bruit, égale à la puissance moyenne de bruit dans ladite zone de M-1 cellules entourant la cellule i.
  8. 8. Procédé selon la revendication 6, caractérisé en ce que ladite référence de bruit est déterminée en faisant le produit par un coefficient K prédéterminé, supérieur à 1, de la moyennne, sur M-1cellules entourant ladite cellule i, des valeurs de la fonction d'autocorrélation au point zéro du signal reçu de chacune de ces M-l cellules selon l'expression
    Figure img00120003
    en calculant la valeur de la fonction d'autocorrélation au point zéro du signal reçu de la cellule i selon l'expression
    Figure img00120004
    et en prenant le maximum de ces deux expressions selon la relation
    Figure img00120005
    où RARC représente la référence de bruit.
  9. 9. Dispositif de détection autorégressive d'un signal de cible monofréquentielle dans un environnement de bruit blanc et d'estimation de sa fréquence à partir du signal échantillonné complexe d'un radar à impulsions, ledit signal radar étant constitué pour chaque cellule de résolution de N échantillons successifs correspondant aux N impulsions radar émises pendant le passage du faisceau radar sur ladite cellule, ledit dispositif mettant en oeuvre le procédé selon l'une des revendications 7 ou 8 et étant caractérisé en ce qu'il comprend - des premiers moyens de calcul (20 à 27) pour calculer, à partir des
    N échantillons reçus relatifs à chaque cellule de résolution de rang m, les parties réelle et imaginaire de l'expression
    Figure img00130001
    éventuellement détecté dans la cellule i.
    mémoires (28, 29), pour fournir la pulsation wb(i) du signal de cible
    dont les parties réelle et imaginaire sont délivrées par lesdites
    présence de cible dans la cellule i en cours de traitement, et - des moyens d'extraction (40 à 42) de l'argument de la quantité complexe
    multiplieur (35) et pour délivrer, le cas échéant, un signal X(i) de
    lesdits moyens de calcul de module (30) au produit fourni par ledit
    RARM) par ledit seuil prédéterminé ;; - un circuit de comparaison (36) pour comparer le module fourni par
    résolution, comprenant la cellule i et les M-I cellules l'entourant ; - un multiplieur (35) pour faire le produit de ladite référence (RARC ;
    lesdits seconds moyens de calcul (31 à 33) pour les M cellules de
    référence de bruit (RARC ;RARM) à partir des valeurs fournies par
    - un circuit de calcul de référence de bruit (34) pour calculer ladite
    Figure img00130002
    fonction d'autocorrélation au point zéro du signal reçu pour chaque cellule de résolution m selon l'expression
    cellule i - des seconds moyens de calcul (31 à 33) pour calculer la valeur de la
    mémoires (28, 29), le module dudit produit correspondant à la
    réelle et imaginaire relatives à la cellule i fournies par lesdites
    traitement - des moyens (30) de calcul de module pour fournir, à partir des parties
    naire et restituer celles correspondant à la cellule i en cours de
    cette cellule - des mémoires (28, 29) pour conserver lesdites parties réelle et imagi
    par la valeur de ladite fonction d'autocorrélation au point zéro pour
    produit du coefficient de réflexion d'ordre 1 pour ladite cellule m
  10. 10. Dispositif de détection autorégressive selon la revendication 9, caractérisé en ce que ledit circuit de calcul de référence de bruit (34) comprend - une mémoire (340) recevant les valeurs de la fonction d'autocorrélation
    au point zéro, fournies par lesdits seconds moyens de calcul (31 à 33)
    et correspondant aux signaux des cellules de résolution successivement
    reçus, et restituant la valeur correspondant à la cellule i en cours
    de traitement ;; - des troisièmes moyens de calcul (341 à 346) pour calculer et mettre à
    jour la somme, sur M-1 cellules entourant ladite cellule i, des valeurs
    de la fonction d'autocorrélation au point zéro selon l'expression
    Figure img00140001
    - un multiplieur (347) pour effectuer le produit de la valeur fournie
    par lesdits troisièmes moyens de calcul par la constante M 1 et fournir
    M-i
    la référence de bruit (RARM) basée sur la puissance moyenne de bruit, selon la relation
    Figure img00140002
    Figure img00150001
    de ces deux valeurs comme référence de bruit (RARC) selon la relation
    la cellule i fournie par ladite mémoire (340) et qui délivre le maximum
    la valeur de la fonction d'autocorrélation au point zéro relative à
    la valeur délivrée par lesdits moyens multiplieurs et d'autre part
    et - un circuit de comparaison et de sélection (349) qui reçoit d'une part
    constante K , s où K est un coefficient prédéterminé, supérieur à 1,
    valeur fournie par lesdits troisièmes moyens de calcul par la
    - des moyens multiplieurs (347, 348) pour effectuer le produit de la
    Figure img00140003
    valeurs de la fonction d'autocorrélation au point zéro selon l'expression
    jour la somme, sur M-1 cellules entourant ladite cellule i, des
    de traitement ; ; - des troisièmes moyens de calcul (341 à 346) pour calculer et mettre à
    reçus, et restituant la valeur correspondant à la cellule i en cours
    et correspondant aux signaux des cellules de résolution successivement
    au point zéro, fournies par lesdits seconds moyens de calcul (31 à 33)
  11. 11.Dispositif de détection autorégressive selon la revendication 9, caractérisé en ce que ledit circuit de calcul de référence de bruit (34) comprend - une mémoire (340) recevant les valeurs de la fonction d'autocorrélation
FR8603234A 1986-03-07 1986-03-07 Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions. Expired - Fee Related FR2678738B1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
FR8603234A FR2678738B1 (fr) 1986-03-07 1986-03-07 Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions.
DE19873707594 DE3707594A1 (de) 1986-03-07 1987-03-09 Verfahren fuer die autoregressive erfassung eines komplexen sinusfoermigen signals in rauschen und die schaetzung seiner frequenz und vorrichtung fuer die durchfuehrung des verfahrens, insbesondere fuer die erfassung eines zieles in einem impulsradarsignal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR8603234A FR2678738B1 (fr) 1986-03-07 1986-03-07 Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions.

Publications (2)

Publication Number Publication Date
FR2678738A1 true FR2678738A1 (fr) 1993-01-08
FR2678738B1 FR2678738B1 (fr) 1994-06-03

Family

ID=9332871

Family Applications (1)

Application Number Title Priority Date Filing Date
FR8603234A Expired - Fee Related FR2678738B1 (fr) 1986-03-07 1986-03-07 Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions.

Country Status (2)

Country Link
DE (1) DE3707594A1 (fr)
FR (1) FR2678738B1 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105640505B (zh) * 2016-01-13 2018-06-19 昆明理工大学 一种基于ar模型谱估计的脉搏信号随机噪声去噪方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
L' ONDE ELECTRIQUE vol. 64, no. 4, Août 1984, PARIS , FR pages 38 - 43 J. L. ZOLESIO 'Filtrage autoregressif d' elimination de fouillis' *

Also Published As

Publication number Publication date
DE3707594A1 (de) 1992-09-10
FR2678738B1 (fr) 1994-06-03

Similar Documents

Publication Publication Date Title
EP0261736B1 (fr) Dispositif d&#39;exploration par echographie ultrasonore d&#39;un milieu en mouvement
EP0749083B1 (fr) Procédé et dispositif de détermination du spectre de fréquence d&#39;un signal
EP0225667B1 (fr) Appareil d&#39;exploration par échographie ultrasonore de milieux en mouvement et notamment d&#39;écoulements sanguins ou d&#39;organes tels que le coeur
FR2696011A1 (fr) Procédé et dispositif de réglage du seuil de détection d&#39;un radar.
FR2678125A1 (fr) Procede et appareillage de commande adaptative d&#39;une camera d&#39;imagerie electronique.
CH625053A5 (fr)
EP0014619B1 (fr) Dispositif de filtrage dynamique non linéaire du bruit de mesure angulaire dans un radar, et système radar le comportant
EP0166836A1 (fr) Procédé de caractérisation de la structure d&#39;un milieu, et dispositif mettant en oeuvre ce procédé
CA1226654A (fr) Dispositif de discrimination d&#39;echos radar
FR2662348A1 (fr) Dispositif de mesure et de visualisation par echographie ultrasonore de debit d&#39;un ecoulement sanguin et de dilatation du vaisseau associe.
FR2678738A1 (fr) Procede de detection autoregressive d&#39;un signal sinusouidal complexe dans du bruit et d&#39;estimation de sa frequence pour un radar a impulsions.
EP0464115A1 (fr) Procede et dispositif d&#39;analyse spectrale en temps reel de signaux instationnaires complexes
EP0821322A1 (fr) Procédé de filtrage temporel du bruit dans une image d&#39;une séquence d&#39;images numérisées et dispositif mettant en oeuvre ce procédé
FR2549957A1 (fr) Procede et appareil pour synthetiser un signal continu estime a partir de segments d&#39;un signal gaussien resultant de la mesure doppler par ultra-sons de l&#39;ecoulement d&#39;un fluide
EP0334711B1 (fr) Dispositif d&#39;élimination du fouillis mobile dans un radar
FR2817973A1 (fr) Methode de detection et de positionnement d&#39;objets basee sur deux etapes de formation numerique de faisceaux d&#39;un reseau phase de capteurs
EP0651263B1 (fr) Procédé de détermination du rang d&#39;ambiguité en distance d&#39;échos radar
EP1081505A1 (fr) Dispositif numérique de corrélation et procédé de corrélation pour systèmes télémétriques
FR2672395A1 (fr) Procede et dispositif de reduction des effets des bruits parasites sur la detection d&#39;une cible par un systeme comprenant une pluralite de capteurs elementaires.
FR3038391B1 (fr) Procede et dispositif de traitement d&#39;un signal spectral
FR3051959A1 (fr) Procede et dispositif pour estimer un signal dereverbere
FR2724028A1 (fr) Procede d&#39;estimation aveugle de retards differentiels entre deux signaux
EP0977051B1 (fr) Procédé de restitution de la sensibilité d&#39;un radar en présence d&#39;une pollution électromagnétique impulsionnelle
FR2689645A1 (fr) Procédé et dispositif de limitation de fausses alarmes par tri de vitesse pour radar doppler cohérent à impulsions à diversité de fréquences de répétition.
FR3030721A1 (fr) Dispositif et procede de detection de poissons dans un dispositif de franchissement

Legal Events

Date Code Title Description
ST Notification of lapse