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 PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2921—Extracting wanted echo-signals based on data belonging to one radar period
- G01S7/2922—Extracting wanted echo-signals based on data belonging to one radar period by using a controlled threshold
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
- G01S13/524—Discriminating 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/534—Discriminating 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.
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.
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.
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
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.
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
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 :
où ss a une valeur supérieure à 0 et où r1 désigne l'espérance mathématique de r1.
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.
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 # :
puisque R(O) est égale à b, puissance moyenne du bruit.
puisque R(O) est égale à b, puissance moyenne du bruit.
Donc, dans les deux cas, il y aura détection de la réalisation # car
alors que cela n'était pas possible avec le coefficient de réflexion r non modifié.
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
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, ...).
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
et, dans tous les cas
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
et, dans tous les cas
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.
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
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
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.
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
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.
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
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
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
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
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
qui est égal à la pulsation #D(i) du signal utile de cible.
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
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.
<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>
<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.
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
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
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 :
où M est, comme on l'a vu, le nombre de cellules de la zone Z (figure 4) entourant la cellule i.
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
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 :
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
On utilise pour cela le circuit CRB de la figure 5 qui fournit d'une part le second terme
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.
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)
- REVENDICATIONSfré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 lacient 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érenceI. 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. 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. 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. 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. 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. 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 deN é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'expressionet 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. 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 cesM-1 cellules selon la relationoù 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. 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'expressionen calculant la valeur de la fonction d'autocorrélation au point zéro du signal reçu de la cellule i selon l'expressionet en prenant le maximum de ces deux expressions selon la relationoù RARC représente la référence de bruit.
- 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 desN échantillons reçus relatifs à chaque cellule de résolution de rang m, les parties réelle et imaginaire de l'expressionéventuellement détecté dans la cellule i.mémoires (28, 29), pour fournir la pulsation wb(i) du signal de cibledont les parties réelle et imaginaire sont délivrées par lesditesprésence de cible dans la cellule i en cours de traitement, et - des moyens d'extraction (40 à 42) de l'argument de la quantité complexemultiplieur (35) et pour délivrer, le cas échéant, un signal X(i) delesdits moyens de calcul de module (30) au produit fourni par leditRARM) par ledit seuil prédéterminé ;; - un circuit de comparaison (36) pour comparer le module fourni parré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 deréférence de bruit (RARC ;RARM) à partir des valeurs fournies par- un circuit de calcul de référence de bruit (34) pour calculer laditefonction d'autocorrélation au point zéro du signal reçu pour chaque cellule de résolution m selon l'expressioncellule i - des seconds moyens de calcul (31 à 33) pour calculer la valeur de lamémoires (28, 29), le module dudit produit correspondant à laréelle et imaginaire relatives à la cellule i fournies par lesditestraitement - des moyens (30) de calcul de module pour fournir, à partir des partiesnaire et restituer celles correspondant à la cellule i en cours decette cellule - des mémoires (28, 29) pour conserver lesdites parties réelle et imagipar la valeur de ladite fonction d'autocorrélation au point zéro pourproduit du coefficient de réflexion d'ordre 1 pour ladite cellule m
- 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élationau point zéro, fournies par lesdits seconds moyens de calcul (31 à 33)et correspondant aux signaux des cellules de résolution successivementreçus, et restituant la valeur correspondant à la cellule i en coursde 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 valeursde la fonction d'autocorrélation au point zéro selon l'expression- un multiplieur (347) pour effectuer le produit de la valeur fourniepar lesdits troisièmes moyens de calcul par la constante M 1 et fournirM-ila référence de bruit (RARM) basée sur la puissance moyenne de bruit, selon la relationde ces deux valeurs comme référence de bruit (RARC) selon la relationla cellule i fournie par ladite mémoire (340) et qui délivre le maximumla valeur de la fonction d'autocorrélation au point zéro relative àla valeur délivrée par lesdits moyens multiplieurs et d'autre partet - un circuit de comparaison et de sélection (349) qui reçoit d'une partconstante 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 lavaleurs de la fonction d'autocorrélation au point zéro selon l'expressionjour la somme, sur M-1 cellules entourant ladite cellule i, desde traitement ; ; - des troisièmes moyens de calcul (341 à 346) pour calculer et mettre àreçus, et restituant la valeur correspondant à la cellule i en courset correspondant aux signaux des cellules de résolution successivementau point zéro, fournies par lesdits seconds moyens de calcul (31 à 33)
- 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
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105640505B (zh) * | 2016-01-13 | 2018-06-19 | 昆明理工大学 | 一种基于ar模型谱估计的脉搏信号随机噪声去噪方法 |
-
1986
- 1986-03-07 FR FR8603234A patent/FR2678738B1/fr not_active Expired - Fee Related
-
1987
- 1987-03-09 DE DE19873707594 patent/DE3707594A1/de not_active Withdrawn
Non-Patent Citations (1)
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'exploration par echographie ultrasonore d'un milieu en mouvement | |
EP0749083B1 (fr) | Procédé et dispositif de détermination du spectre de fréquence d'un signal | |
EP0225667B1 (fr) | Appareil d'exploration par échographie ultrasonore de milieux en mouvement et notamment d'écoulements sanguins ou d'organes tels que le coeur | |
FR2696011A1 (fr) | Procédé et dispositif de réglage du seuil de détection d'un radar. | |
FR2678125A1 (fr) | Procede et appareillage de commande adaptative d'une camera d'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'un milieu, et dispositif mettant en oeuvre ce procédé | |
CA1226654A (fr) | Dispositif de discrimination d'echos radar | |
FR2662348A1 (fr) | Dispositif de mesure et de visualisation par echographie ultrasonore de debit d'un ecoulement sanguin et de dilatation du vaisseau associe. | |
FR2678738A1 (fr) | Procede de detection autoregressive d'un signal sinusouidal complexe dans du bruit et d'estimation de sa frequence pour un radar a impulsions. | |
EP0464115A1 (fr) | Procede et dispositif d'analyse spectrale en temps reel de signaux instationnaires complexes | |
EP0821322A1 (fr) | Procédé de filtrage temporel du bruit dans une image d'une séquence d'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'un signal gaussien resultant de la mesure doppler par ultra-sons de l'ecoulement d'un fluide | |
EP0334711B1 (fr) | Dispositif d'élimination du fouillis mobile dans un radar | |
FR2817973A1 (fr) | Methode de detection et de positionnement d'objets basee sur deux etapes de formation numerique de faisceaux d'un reseau phase de capteurs | |
EP0651263B1 (fr) | Procédé de détermination du rang d'ambiguité en distance d'é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'une cible par un systeme comprenant une pluralite de capteurs elementaires. | |
FR3038391B1 (fr) | Procede et dispositif de traitement d'un signal spectral | |
FR3051959A1 (fr) | Procede et dispositif pour estimer un signal dereverbere | |
FR2724028A1 (fr) | Procede d'estimation aveugle de retards differentiels entre deux signaux | |
EP0977051B1 (fr) | Procédé de restitution de la sensibilité d'un radar en présence d'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 |