Brevet FR2678738A1 — publié le 19930108

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 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. 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 # : 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 2) Si b+ss > Kb 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é. 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, ...). 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 Là aussi, on obtient une détection de la réalisation # puisque : 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 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. 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. 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 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. 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. SIGNE PARTIE REELLE SIGNE PARTIE IMAGINAIRE VALEUR DE + + + n-p + - # - - # - # 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. 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. 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. REVENDICATIONS 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 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 sinusoïdal, et - calculer l'argument dudit coefficient de réflexion pour en déduire la fréquence dudit signal sinusoidal. 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 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 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. 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 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. 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 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 et en prenant le maximum de ces deux expressions selon la relation où 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 des N échantillons reçus relatifs à chaque cellule de résolution de rang m, les parties réelle et imaginaire de l'expression produit du coefficient de réflexion d'ordre 1 pour ladite cellule m par la valeur de ladite fonction d'autocorrélation au point zéro pour cette cellule - des mémoires (28, 29) pour conserver lesdites parties réelle et imagi naire et restituer celles correspondant à la cellule i en cours de traitement - des moyens (30) de calcul de module pour fournir, à partir des parties réelle et imaginaire relatives à la cellule i fournies par lesdites mémoires (28, 29), le module dudit produit correspondant à la cellule i - des seconds moyens de calcul (31 à 33) pour calculer la valeur de la fonction d'autocorrélation au point zéro du signal reçu pour chaque cellule de résolution m selon l'expression - un circuit de calcul de référence de bruit (34) pour calculer ladite référence de bruit (RARC ;RARM) à partir des valeurs fournies par lesdits seconds moyens de calcul (31 à 33) pour les M cellules de 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 ; RARM) par ledit seuil prédéterminé ;; - un circuit de comparaison (36) pour comparer le module fourni par lesdits moyens de calcul de module (30) au produit fourni par ledit multiplieur (35) et pour délivrer, le cas échéant, un signal X(i) de 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 dont les parties réelle et imaginaire sont délivrées par lesdites mémoires (28, 29), pour fournir la pulsation wb(i) du signal de cible éventuellement détecté dans la cellule i. 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 - 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 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 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 - des moyens multiplieurs (347, 348) pour effectuer le produit de la valeur fournie par lesdits troisièmes moyens de calcul par la constante K , s où K est un coefficient prédéterminé, supérieur à 1, et - un circuit de comparaison et de sélection (349) qui reçoit d'une part la valeur délivrée par lesdits moyens multiplieurs et d'autre part la valeur de la fonction d'autocorrélation au point zéro relative à la cellule i fournie par ladite mémoire (340) et qui délivre le maximum de ces deux valeurs comme référence de bruit (RARC) selon la relation