Corrections
Rappels de probabilité ¶ Exercice 1 ¶ E [ x ] = 1 2 et σ = 1 12 \mathbb{E}[x] = \frac{1}{2}
\qquad\text{et}\qquad
\sigma = \frac{1}{\sqrt{12}} E [ x ] = 2 1 et σ = 12 1 Exercice 2 ¶ E ^ [ x ] = 14 , 5 et σ ^ [ x ] = 21 3 ≈ 2 , 646 \hat{\mathbb{E}}[x] = 14,5
\qquad\text{et}\qquad
\hat{\sigma}[x] = \sqrt{\frac{21}{3}} \approx 2,646 E ^ [ x ] = 14 , 5 et σ ^ [ x ] = 3 21 ≈ 2 , 646 Le bruit dans les signaux ¶ Exercice 3 ¶ Pour rappel, l’espérance d’un bruit blanc est nulle.
Par définition de l’estimateur de la variance, l’écart-type σ b \sigma_b σ b du bruit est égal à la racine carrée de la puissance du bruit multipliée par N N − 1 \frac{N}{N-1} N − 1 N puisque
σ b 2 = 1 N − 1 ∑ n = 1 N b n 2 = N N − 1 P b . \sigma_b^2 = \frac{1}{N-1} \sum_{n=1}^N b_n^2 = \frac{N}{N-1}P_b. σ b 2 = N − 1 1 n = 1 ∑ N b n 2 = N − 1 N P b . Expression de l’écart-type du bruit en fonction du RSB et des échantillons du signal x x x :
σ b = 1 0 − RSB / 20 ∑ n x n 2 N − 1 . \sigma_b = 10^{-\text{RSB}/20} \sqrt{ \frac{\sum_n x_n^2}{N-1} }. σ b = 1 0 − RSB /20 N − 1 ∑ n x n 2 . Traitement des signaux bruités ¶ Exercice 4 ¶ L’intercorrélation R y x R_{yx} R y x est définie par :
R y x [ n ] = ∑ m = − ∞ + ∞ y [ n + m ] x [ m ] . R_{yx}[n] = \sum_{m=-\infty}^{+\infty} y[n+m] x[m]. R y x [ n ] = m = − ∞ ∑ + ∞ y [ n + m ] x [ m ] . Par changement de variable k = n + m ⇔ m = k − n k=n+m \Leftrightarrow m=k-n k = n + m ⇔ m = k − n , l’intercorrélation peut s’écrire :
R y x [ n ] = ∑ n = − ∞ + ∞ y [ k ] x [ k − n ] . R_{yx}[n] = \sum_{n=-\infty}^{+\infty} y[k] x[k-n]. R y x [ n ] = n = − ∞ ∑ + ∞ y [ k ] x [ k − n ] . Le produit de convolution a pour définition :
( y ∗ h ) [ n ] = ∑ k = − ∞ + ∞ y [ k ] h [ n − k ] (y*h)[n] = \sum_{k=-\infty}^{+\infty} y[k] h[n-k] ( y ∗ h ) [ n ] = k = − ∞ ∑ + ∞ y [ k ] h [ n − k ] Il y a égalité entre le produit de convolution et l’intercorrélation lorsque h [ n − k ] = x [ k − n ] h[n-k]=x[k-n] h [ n − k ] = x [ k − n ] ,
ce qui signifie que la réponse impulsionnelle h [ n ] h[n] h [ n ] du filtre est égale à x [ − n ] x[-n] x [ − n ] .
Le filtre adapté est donc bien un filtre dont la réponse impulsionnelle est le signal x x x retourné dans le temps.
Exercice 5 ¶ La sortie du filtre moyenneur est :
z [ n ] = 1 M ∑ m = n − M − 1 2 n + M − 1 2 y [ m ] = ∑ m = n − M − 1 2 n + M − 1 2 y [ m ] 1 M = ∑ m = n − M − 1 2 n + M − 1 2 y [ m ] h [ n − m ] z[n] = \frac{1}{M} \sum_{m=n-\frac{M-1}{2}}^{n+\frac{M-1}{2}} y[m]
= \sum_{m=n-\frac{M-1}{2}}^{n+\frac{M-1}{2}} y[m] \frac{1}{M}
= \sum_{m=n-\frac{M-1}{2}}^{n+\frac{M-1}{2}} y[m] h[n-m] z [ n ] = M 1 m = n − 2 M − 1 ∑ n + 2 M − 1 y [ m ] = m = n − 2 M − 1 ∑ n + 2 M − 1 y [ m ] M 1 = m = n − 2 M − 1 ∑ n + 2 M − 1 y [ m ] h [ n − m ] où
h [ n − m ] = { 1 / M si m ∈ { n − M − 1 2 , … , n + M − 1 2 } 0 sinon = { 1 / M si n − m ∈ { − M − 1 2 , … , M − 1 2 } 0 sinon \begin{align*}
h[n-m]
&=
\begin{cases}
1/M \text{ si } m \in \{n-\frac{M-1}{2},\dots,n+\frac{M-1}{2}\} \\
0 \text{ sinon}
\end{cases} \\
&=
\begin{cases}
1/M \text{ si } n-m \in \{-\frac{M-1}{2},\dots,\frac{M-1}{2}\} \\
0 \text{ sinon}
\end{cases}
\end{align*} h [ n − m ] = { 1/ M si m ∈ { n − 2 M − 1 , … , n + 2 M − 1 } 0 sinon = { 1/ M si n − m ∈ { − 2 M − 1 , … , 2 M − 1 } 0 sinon soit
h [ n ] = { 1 / M si n ∈ { − M − 1 2 , … , M − 1 2 } 0 sinon h[n] =
\begin{cases}
1/M \text{ si } n \in \{-\frac{M-1}{2},\dots,\frac{M-1}{2}\} \\
0 \text{ sinon}
\end{cases} h [ n ] = { 1/ M si n ∈ { − 2 M − 1 , … , 2 M − 1 } 0 sinon Le filtre moyenneur s’écrit donc bien comme la convolution entre le signal observé y y y et une porte h h h .
Exercice 6 ¶ La moyenne des mesures y i y_i y i est :
∀ n , 1 I ∑ i = 1 I y i [ n ] = 1 I ∑ i = 1 I x [ n ] ⏟ x ~ [ n ] + 1 I ∑ i = 1 I b i [ n ] ⏟ b ~ [ n ] \forall n,\qquad
\frac{1}{I} \sum_{i=1}^I y_i[n]
= \underbrace{ \frac{1}{I} \sum_{i=1}^I x[n] }_{\tilde{x}[n]} + \underbrace{ \frac{1}{I} \sum_{i=1}^I b_i[n] }_{\tilde{b}[n]} ∀ n , I 1 i = 1 ∑ I y i [ n ] = x ~ [ n ] I 1 i = 1 ∑ I x [ n ] + b ~ [ n ] I 1 i = 1 ∑ I b i [ n ] où x ~ [ n ] = x [ n ] \tilde{x}[n] = x[n] x ~ [ n ] = x [ n ] .
On montre déjà que l’espérance de b ~ \tilde{b} b ~ est nulle (rappelons que les b i b_i b i sont blanc et gaussien, donc d’espérance nulle) :
E [ b ~ ] = E [ 1 I ∑ i = 1 I b i [ n ] ] = 1 I ∑ i = 1 I E [ b i [ n ] ] = 1 I ∑ i = 1 I 0 = 0 \mathbb{E}\left[\tilde{b}\right]
= \mathbb{E}\left[\frac{1}{I} \sum_{i=1}^I b_i[n]\right]
= \frac{1}{I} \sum_{i=1}^I \mathbb{E}\left[b_i[n]\right]
= \frac{1}{I} \sum_{i=1}^I 0
= 0 E [ b ~ ] = E [ I 1 i = 1 ∑ I b i [ n ] ] = I 1 i = 1 ∑ I E [ b i [ n ] ] = I 1 i = 1 ∑ I 0 = 0 On peut donc calculer la variance de b ~ \tilde{b} b ~ :
V [ b ~ ] = E [ b ~ 2 ] − E [ b ~ ] 2 = E [ b ~ 2 ] − 0 = E [ 1 I 2 ( ∑ i = 1 I b i [ n ] ) 2 ] = 1 I 2 ∑ i = 1 I E [ b i [ n ] 2 ] + 1 I 2 ∑ i = 1 I ∑ j ≠ i E [ b i [ n ] b j [ n ] ] \begin{align*}
\mathbb{V}\left[\tilde{b}\right]
&= \mathbb{E}\left[\tilde{b}^2\right] - \mathbb{E}\left[\tilde{b}\right]^2 \\
&= \mathbb{E}\left[\tilde{b}^2\right] - 0 \\
&= \mathbb{E}\left[\frac{1}{I^2} \left( \sum_{i=1}^I b_i[n] \right)^2 \right] \\
&= \frac{1}{I^2} \sum_{i=1}^I \mathbb{E}\left[b_i[n]^2\right] + \frac{1}{I^2} \sum_{i=1}^I \sum_{j \neq i} \mathbb{E}\left[b_i[n]b_j[n]\right]
\end{align*} V [ b ~ ] = E [ b ~ 2 ] − E [ b ~ ] 2 = E [ b ~ 2 ] − 0 = E ⎣ ⎡ I 2 1 ( i = 1 ∑ I b i [ n ] ) 2 ⎦ ⎤ = I 2 1 i = 1 ∑ I E [ b i [ n ] 2 ] + I 2 1 i = 1 ∑ I j = i ∑ E [ b i [ n ] b j [ n ] ] or, puisque les échantillons b i b_i b i du bruit sont indépendants, alors E [ b i [ n ] b j [ n ] ] = 0 \mathbb{E}\left[b_i[n]b_j[n]\right]=0 E [ b i [ n ] b j [ n ] ] = 0 quels que soient i i i et j j j .
Il reste donc :
V [ b ~ ] = 1 I 2 ∑ i = 1 I E [ b i [ n ] 2 ] = 1 I V [ b i ] \mathbb{V}\left[\tilde{b}\right]
= \frac{1}{I^2} \sum_{i=1}^I \mathbb{E}\left[b_i[n]^2\right]
= \frac{1}{I} \mathbb{V}\left[b_i\right] V [ b ~ ] = I 2 1 i = 1 ∑ I E [ b i [ n ] 2 ] = I 1 V [ b i ] Ce derner résultat implique que le RSB du signal moyenné RSB ~ \widetilde{\text{RSB}} RSB
est égal au RSB d’un signal y i y_i y i plus un facteur 10 log ( I ) 10\log(I) 10 log ( I ) :
RSB ~ = 10 log 10 ( P x ~ P b ~ ) = 10 log 10 ( P x ~ 1 I P b i ) = 10 log 10 ( I ) + 10 log 10 ( P x ~ P b i ) = 10 log 10 ( I ) + RSB \begin{align*}
\widetilde{\text{RSB}}
&= 10 \log_{10} \left(\frac{P_{\tilde{x}}}{P_{\tilde{b}}}\right) \\
&= 10 \log_{10} \left(\frac{P_{\tilde{x}}}{\frac{1}{I}P_{b_i}}\right) \\
&= 10 \log_{10}(I) + 10 \log_{10} \left(\frac{P_{\tilde{x}}}{P_{b_i}}\right) \\
&= 10 \log_{10}(I) + \text{RSB}
\end{align*} RSB = 10 log 10 ( P b ~ P x ~ ) = 10 log 10 ( I 1 P b i P x ~ ) = 10 log 10 ( I ) + 10 log 10 ( P b i P x ~ ) = 10 log 10 ( I ) + RSB Exercice 7 ¶ θ = 2 , 5 \boldsymbol{\theta}=2,5 θ = 2 , 5 est la valeur de la constante qui approxime au mieux les données.
On peut calculer J = 0 , 5 \mathcal{J}=0,5 J = 0 , 5 .
θ = ( 7 3 1 ) \boldsymbol{\theta}=\begin{pmatrix}\frac{7}{3} \\ 1\end{pmatrix} θ = ( 3 7 1 ) regroupe les coefficients de la droite qui approxime au mieux les données.
Celle-ci est donc d’équation z = 7 3 + n z = \frac{7}{3} + n z = 3 7 + n .
On peut calculer J = 2 3 \mathcal{J}=\frac{2}{3} J = 3 2 .