Théorie de la détection du signal

La théorie de la détection du signal est une application de la théorie de la décision concernant la détection d'une information ou d'un motif dans un signal mesuré en présence de bruit. Cette théorie décrit comment choisir de manière optimale un critère permettant de discerner l'information jugée utile - le signal - du bruit environnant, en tenant compte de la nature probabiliste de la source du signal, des différentes sources de bruit et du récepteur, à savoir l'outil de mesure du signal. Elle définit en outre des critères de qualité de la détection obtenue et plus généralement de la performance du récepteur.

Les débuts de la théorie de la détection du signal remontent au développement des récepteurs radar.

Les principales applications de cette théorie sont, outre le radar, les télécommunications, le contrôle qualité, et l'astronomie d'observation.

Cas de la détection binaire

C'est le cas le plus simple, dans lequel on cherche à différencier deux états possibles, généralement la présence ou l'absence d'une information dans un ou plusieurs signaux déterministes en présence de bruit aléatoire. Cela revient aussi à un problème de classification automatique supervisée dans le cas où il y a deux classes. On considère un ensemble de n signaux pouvant être émis aléatoirement par l'une de deux sources H₁ et H₂, avec des probabilités p₁ et p₂ (telles que $p_1+p_2=1$), appelées probabilités a priori. Ceci produit une variable aléatoire continue X sur un espace d'observation $\chi \subset {ℝ}^n$ dont les probabilités conditionnelles sont définies par les fonctions densité de probabilité f_X|H1 et f_X|H2.

On cherche une partition de l'espace d'observation $\chi ={\chi }_1\cup {\chi }_2$ telle que ${\chi }_1\cap {\chi }_2=\mathrm{\varnothing }$. Etant donnée une réalisation x de X, c'est à dire une mesure du signal, on lui associe les décisions :

• D₁ si x est dans ${\chi }_1$
• D₂ si x est dans ${\chi }_2$

Relations entre les hypothèses, l'espace d'observation et les décisions

Il y a donc quatre types d'événement possibles :

• Décision D₁ sous l'hypothèse H₁, appelé vrai négatif
• Décision D₁ sous l'hypothèse H₂, appelé faux négatif
• Décision D₂ sous l'hypothèse H₁, appelé faux positif
• Décision D₂ sous l'hypothèse H₂, appelé vrai positif

auxquels sont associés les probabilités conditionnelles $P(D_i|H_j)=\underset{{\chi }_i}{\int }{f}_{X|H_j}(x)dx$.

La règle de décision se base sur le rapport de vraisemblance $\mathrm{\Lambda }(x)=\frac{f_{X|H_2}(x)}{f_{X|H_1}(x)}$ et une valeur de seuil ${\lambda }_s$ tels que

$\mathrm{\Lambda }(x)\ge {\lambda }_s$ conduit à prendre la décision D₂

et

$\mathrm{\Lambda }(x)<{\lambda }_s$ conduit à prendre la décision D₁

.

Cette règle permet d'exprimer les probabilités conditionnelles de décision en fonction du rapport de vraisemblance :

• $P(D_1|H_j)=\underset{0}{\overset{{\lambda }_s}{\int }}{f}_{\mathrm{\Lambda }|H_j}(\lambda )d\lambda$
• $P(D_2|H_j)=\underset{{\lambda }_s}{\overset{\mathrm{\infty }}{\int }}{f}_{\mathrm{\Lambda }|H_j}(\lambda )d\lambda$

avec $f_{\mathrm{\Lambda }|H_j}$ la fonction densité de probabilité associée à la probabilité conditionnelle de la variable aléatoire $\mathrm{\Lambda }$ sous hypothèse H_j.

La performance de la détection est quantitifée au moyen :

• des probabilités totales de bonne ou mauvaise détection $P(BD)=P(D_1|H_1)p_1+P(D_2|H_2)p_2$ et $P(MD)=P(D_2|H_1)p_1+P(D_1|H_2)p_2$, qui dépendent du niveau de seuil ${\lambda }_s$ choisi,
• de la caractéristique opérationnelle du récepteur représentant graphiquement l'évolution du taux de vrai positifs $P(D_2|H_2)$ et du taux de faux positifs $P(D_2|H1)$ en fonction du seuil ${\lambda }_s$.

Critère de Bayes

Toutes les décisions ne sont pas équivalentes quant à leur utilité. On associe donc un coût à chaque événement $C_{i,j}$ correspondant à la décision D_i sous hypothèse H_j. En général les mauvaises décisions sont plus coûteuses et on a $C_{2,1}>C_{1,1}$ et $C_{1,2}>C_{2,2}$. On cherche ensuite le niveau de seuil ${\lambda }_s$ qui minimise le coût moyen

$C=\sum _{i=1}^2\sum _{j=1}^2{C}_{i,j}P(D_i|H_j)p_j$

.

Celui ci est obtenu pour

${\lambda }_s=\frac{(C_{1,1}-C_{2,1})p_1}{(C_{2,2}-C_{1,2})p_2}$

c'est le seuil obtenu selon le critère de Bayes.

Maximum a posteriori

Si on choisit de pénaliser uniquement et uniformément les mauvaises détections $C_{1,2}=C_{2,1}=1$ et $C_{1,1}=C_{2,2}=0$, le coût moyen devient $C=P(D_2|H_1)p_1+P(D_1|H_2)p_2=P(MD)$. Ce choix particulier revient à minimiser la probabilité totale de mauvaise détection. Le niveau de seuil devient alors

${\lambda }_s=\frac{p_1}{p_2}$

On peut montrer que dans ce cas la règle de décision peut se reformuler sous la forme suivante :

$P(H_2|X)\ge P(H_1|X)$ conduit à la décision D₂ $P(H_1|X)>P(H_2|X)$ conduit à la décision D₁

$P(H_j|X)$ étant la probabilité conditionnelle a posteriori de l'hypothèse H_j.

Cette règle revient donc à choisir à chaque fois le maximum de la probabilité a posteriori.

Critère minimax

Le critère de Bayes suppose la connaissance des probabilités a priori, ce qui n'est pas toujours le cas. On va donc essayer de se placer dans le pire cas du critère de Bayes, c'est à dire rechercher une des probabilités a priori, par exemple $p_1$, de sorte à avoir

$\underset{p_1}{\max }\underset{{\lambda }_s(p_1)}{\min }C(p_1,{\lambda }_s)$

$C(p_1,{\lambda }_s)$ étant le coût moyen précédemment défini.

La valeur de seuil reste celle définie par le critère de Bayes

${\lambda }_s(p_1)=\frac{(C_{1,1}-C_{2,1})p_1}{(C_{2,2}-C_{1,2})(1-p_1)}$

et la probabilité qui maximise le coût $C(p_1,\lambda (p_1))$ est solution de l'équation caractéristique

$C_{2,2}-C_{1,1}+(C_{1,2}-C_{2,2})P(D_1|H_2)-(C_{2,1}-C_{1,1})P(D_2|H_1)=0$

Les deux probabilités conditionnelles $P(D_i|H_j),i\ne j$ dépendent implicitement de $p_1$ au travers de ${\lambda }_s$.

On peut représenter graphiquement le critère minimax en remarquant que puisque $P(D_1|H_2)+P(D_2|H_2)=1$ l'équation caractéristique est aussi l'équation d'une droite dans le plan $(P(D_2|H_1),P(D_2|H_2))$. Son intersection avec la courbe caractéristique opérationnelle du récepteur donne directement les probabilités conditionnelles satisfaisant le critère minimax.

Intersection de la droite représentative de l'équation caractéristique avec la courbe caractéristique opérationnelle du récepteur

Critère de Neyman-Pearson

Les critères précédents supposent tous qu'il est possible de définir le coût associé aux différentes décisions. Toutefois lorsque ceci n'est pas possible, il existe une autre approche consistant à fixer a priori le taux de faux positif $P(D_2|H_1)$ à une valeur $\alpha$ et à rendre maximum le taux de bonne détection $P(D_2|H_2)$.

Ceci revient à chercher un niveau de seuil ${\lambda }_s$ intervenant dans la règle de décision

$\mathrm{\Lambda }(x)\gtrless {\lambda }_s$

tel que

$P(D_2|H_1)=\underset{{\lambda }_s}{\overset{\mathrm{\infty }}{\int }}{f}_{\mathrm{\Lambda }|H_1}(\lambda )d\lambda =\alpha$

On peut interpréter graphiquement le critère de Neyman-Pearson à l'aide de la courbe caractéristique opérationnelle du récepteur. Le taux de faux positif en abscisse est directement donné par le niveau désiré $\alpha$. D'autre part, une propriété de la caractéristique opérationnelle du récepteur est que sa pente est égale au niveau de seuil ${\lambda }_s$. Le critère de Neyman-Pearson revient donc à mesurer la pente de la courbe au point d'abscisse $\alpha$ considéré.

Application au cas d'un signal en présence de bruit gaussien

On se place dans le cas où une seule mesure (n = 1) est produite par deux hypothèses

• H₁ : la variable aléatoire X est distribuée suivant une loi normale d'espérance mathématique ${\mu }_1$ et de variance ${\sigma }^2$,
• H₂ : la variable aléatoire X est distribuée suivant une loi normale d'espérance ${\mu }_2$ (que l'on suppose plus grande que ${\mu }_1$) et de même variance ${\sigma }^2$.

Autrement dit, le signal peut prendre deux valeurs déterministes ${\mu }_1$ et ${\mu }_2$ auxquelles un bruit gaussien est ajouté.

Le critère du maximum a posteriori fournit la règle de décision suivante sur la valeur mesurée x :

$x\ge x_s$ conduit à la décision D₂ $x<x_s$ conduit à la décision D₁

avec une valeur de seuil $x_s=\frac{{\mu }_1+{\mu }_2}{2}+\frac{{\sigma }^2}{{\mu }_2-{\mu }_1}\log \left(\frac{p_1}{p_2}\right)$

Fonctions densité de probabilité conditionnelle pour les deux hypothèses, distribuées suivant deux lois normales de même variance et d'espérances ${\mu }_1$ et ${\mu }_2$ et probabilités conditionnelles de mauvaise détection $P(D_i|H_j),i\ne j$ pour une valeur $x_s$ du seuil de détection

Les probabilités conditionnelles sont données ici par :

• $P(D_1|H_j)=\mathrm{\Phi }(\frac{x_s-{\mu }_j}{\sigma })$
• $P(D_2|H_j)=1-\mathrm{\Phi }(\frac{x_s-{\mu }_j}{\sigma })$

dans lesquelles $\mathrm{\Phi }$ est la fonction de répartition de la loi normale centrée réduite, soit encore $\mathrm{\Phi }(x)=\frac{1+erf(x/\sqrt{2})}{2}$, erf étant la fonction d'erreur.

La probabilité totale de mauvaise détection résultant de ce choix est donc

$\begin{array}{cc}P(MD)& =P(D_2|H_1)p_1+P(D_1|H_2)p_2\\ & =(1-\mathrm{\Phi }(\frac{x_s-{\mu }_1}{\sigma }))p_1+\mathrm{\Phi }(\frac{x_s-{\mu }_2}{\sigma })p_2\end{array}$

Une autre facçon d'évaluer la détection est de tracer la courbe de la caractéristique opérationnelle du récepteur, c'est à dire $P(D_2|H_2)$ en fonction de $P(D_2|H_1)$ lorsqu'on fait varier $x_s$ de $-\mathrm{\infty }$ à $+\mathrm{\infty }$.

Caractéristique opérationnelle du récepteur correspondant à différentes valeurs du ratio $\frac{\sigma }{{\mu }_2-{\mu }_1}$

Notes et références

• Modèle:Ouvrage
• Modèle:Ouvrage

Voir aussi

• Courbe ROC
• Lemme de Neyman-Pearson