ELE103 PartI v23
ELE103 PartI v23
ELE103 PartI v23
D. ROVIRAS
[email protected]
[email protected]
Tel : 01 40 27 25 67
Accs 11, 2me tage, Bureau 11B2.37
2015-2016 version 23
1. Introduction
2. Le traitement du signal
Station de radio n2
Station de radio n1
t t
Traitement Traitement
Modulation Modulation
Amplification Amplification
Transmission Transmission
mon poste
de radio
CNAM ELE 103 D. Roviras 7
Introduction
A quoi sert le cours de
Bases de Traitement du Signal ?
Pr-amplification
Traitement
Dmodulation
Amplification
(5) Dmodulation
Entre Sortie
Systme
Relations entre/sortie ?
Codeur Codeur
CAN
source canal
Source
ELE102-103-112 Modulateur
Dcodeur Dcodeur
CNA
source canal
Signal ELE102-103 ELE102-112 ELE102-112
ELE103 CNAM ELE 103 D. Roviras 11
Introduction
Rappels sur le filtrage et
les distributions
1. Filtrage par un SLIT
Linarit :
x1 (t ) SLIT y1 (t ) x2 (t ) SLIT y2 (t )
Invariance temporelle :
x1 (t ) SLIT y1 (t ) x1 (t ) SLIT y1 (t )
y : amplitude de la sortie
A
0 x: amplitude de lentre
x1(t)=constante = 1 y1(t)=constante=A
x2(t)=constante = 2 y2(t)=constante=A
x3(t)= x1(t)+x2(t) y3 (t ) = A y1 (t ) + y2 (t )
y : amplitude de la sortie
Pente de 1
0 x: amplitude de lentre
x1(t)=constante = -1 y1(t)=constante=0
x2(t)=constante = 1 y2(t)=constante=1
x3(t)= x1(t)+x2(t) y3 (t ) = 0 y1 (t ) + y2 (t )
x(t) y(t)
g(t)
x1 (t ) y1 (t ) = g (t ) x1 (t )
x2 (t ) = x1 (t ) y2 (t ) = g (t ) x1 (t ) y1 (t )
x(t) y(t)
g(t)
g(t)
2
1
t
x1(t) x2(t)=x1(t-)
y1(t) y2(t)
t
CNAM ELE 103 D. Roviras 17
Filtrage
Systme Linaire Invariant Temporellement (SLIT)
Caractrisation dun SLIT : Rponse frquentielle H(f)
H( f )
f f
Phase de(H ( f 0 ) ) = 1
A1
H ( f0 ) =
A0
CNAM ELE 103 D. Roviras 18
Filtrage
Systme Linaire Invariant Temporellement (SLIT)
Caractrisation dun SLIT : Rponse impulsionnelle h(t)
h(t)
(t ) SLIT h(t )
1/
(t ) = lim
0
t
0
H ( f ) = TF [h(t )]
CNAM ELE 103 D. Roviras 19
Filtrage
Relation entre/sortie dun SLIT
x(t ) SLIT y (t )
h(t)
x(t )
Rponse au rectangle de
largeur D et de hauteur
t
x(k.D): V(t-kD)
Invariance temporelle
et linarit
xe (t ) ye (t )
SLIT
h(t)
t t
k.D
Rectangle de largeur D et de hauteur x(k.D)
1 +
D (t ) tend vers (t)
D ye (t ) x(kD).D.h(t kD)
k =
1
V (t ) tend vers h(t) +
D
y (t ) =
x( ).h(t )d
CNAM ELE 103 D. Roviras 21
Filtrage
Relation entre/sortie dun SLIT
+
y (t ) = x( ).h(t )d = [x * h](t )
De faon ne pas alourdir les notations on crira plus simplement :
y (t ) = x(t ) * h(t )
On verra dans le chapitre suivant que lon a :
Y ( f ) = X ( f ).H ( f )
Transforme de
TF [x(t ) * h(t )] = X ( f ).H ( f )
Fourier
Transforme de TF [x(t ).h(t )] = X ( f ) * H ( f )
Fourier
h( )
h(t0 )
+
t0
x( ).h(t
0 )d
y (t )
t
CNAM ELE 103 D. Roviras 24
Filtrage
Notion de causalit
h( )
Remarque : dans certains calculs thoriques nous utiliserons des filtres non
causaux. Cest une libert qui permet de simplifier les calculs mais il est bien
clair que ces systmes non causaux sont non ralisables physiquement.
Soit x(t) un rectangle centr sur 0 de largeur D et de hauteur 1/D. Laire de x(t) vaut 1.
En faisant tendre D vers 0 on a x(t) qui tend vers une impulsion baptise impulsion de
Dirac
1
D (t ) (t )
D
1/D
D tend vers 0
-D/2 D/2 t 0 t
Dirac .Fonction
x(t ). (t ) = x(0). (t )
Dirac.Fonction x(t ). (t t0 ) = x(t0 ). (t t0 )
c
b
Re
a
b
phase de c = arctg = Angle du vecteur c
a
*
An des + +
An = An*
complexes n = n =
An des *
+ +
An = ( An )
*
complexes
n = n =
exp( j ) = cos( ) + j sin( ) exp( 2jft ) = cos( 2ft ) + j. sin( 2ft )
exp( 2jft ) = 1
(exp(2jft ) )* = exp(2jft )
exp( a ). exp(b) = exp( a + b)
d
(exp( K .t ) ) = K . exp( K .t )
dt
CNAM ELE 103 D. Roviras 30
Filtrage
Rappel sur les fonctions trigonomtriques
sin( a ) sin(b) = 1 cos(a b) 1 cos(a + b)
2 2
cos(a ) cos(b) = 1 cos(a b) + 1 cos(a + b)
2 2
sin( a ) cos(b) = 1 sin( a b) + 1 sin( a + b)
2 2
cos 2 (a ) = (1 + cos(2a )) 2
sin 2 ( a ) = (1 cos( 2a )) 2
1
sin(u)du = cos(u) sin( k .u ) du =
k
. cos( k .u)
2. Transforme de Fourier
x(to)=x(to+T)
On a un spectre de raies
+
n
x(t ) = X n exp 2j t
n = T
to +T
1 n
Xn = x(t ) exp 2j t dt
T to T
Rappel : ( j ) = 1
2
a0 + n n
x(t ) = + an cos 2 t + bn sin 2 t
2 n =1 T T
to +T
2 n
an = x(t ) cos 2 t dt
T to T
to +T
2 n
bn = x(t ) sin 2 t dt
T to T
CNAM ELE 103 D. Roviras 36
SF et TF
Quelques proprits des coefficients de Fourier
Rappel : (a + jb ) = a jb
*
x(t ) relle X n = X *
n
Xn Phase( X n )
n n
5
3 1
1
3
9
7 Signal sinusodal redress
en simple alternance:
21 1, 3 et 9 harmoniques
Signal carr: 1, 3, 5, 7 et 21 harmoniques
Identit de Parseval :
to +T + 2
1
Puissance = x(t ) dt = X
2
n
T to n =
T to
T to n = T n = T
to +T to +T to +T +
1 1 + 2 1 n p
x(t ) dt = n + n p
2 *
X dt X X exp 2 j t dt
T to
T to n = T to nn ,pp= T
to +T + Intgrale gale 0
1
x(t ) dt =
2 2
Xn
T to n =
+
n
x(t ) = X n exp 2j t
n = T
p + n p
x(t ) exp 2j t = X n exp 2j t exp 2j t
T n = T T
to +T to +T +
p n p
x(t ) exp 2j t dt =
T
n =
X n exp 2j t exp 2j t dt
T T
to to
to +T to +T +
n p
= X p dt +
n =
X n exp 2j t exp 2j t dt
T T
to to
Intgrale gale 0
-D/2 D/2 t
-T 0
T
nD
sin
AD T AD nD sin(x)
Xn = = sin c Dfinition : sin c( x) =
T nD T T x
AD/T
0 f
-2/D -1/D -1/T 1/T 2/T 3/T 1/D 2/D
+
X(f ) = x (t ) exp ( 2jft )dt
+
x(t ) = X ( f ) exp (2jft )df
+ + + +
= x(t ) exp( 2jft )dt X ( f )* df = X ( f )X ( f )* df = X ( f ) 2 df
dt n
Symtrie X(t) x(-f)
Multiplication x(t)y(t) X ( f )* Y ( f )
Convolution x(t)*y(t) X ( f ).Y ( f )
x(t ) cos(2f ot )
1 1
Modulation 2
X ( f fo ) + X ( f + fo )
2
SF et TF
CNAM ELE 103 D. Roviras 48
Proprits de la Transforme de Fourier :
Dmonstrations :
Conjugaison : changement de variable
Fonction relle : changement de variable
Dcalage temporel : changement de variable
Dcalage frquentiel : changement de variable
Multiplication et convolution
+
TF ( x(t ) * y (t ) ) = (x(t ) * y(t ))exp( 2jft )dt
+ + +
= ( x(t ) * y (t ) ) exp( 2jft )dt = (x( ) y (t ) )d exp( 2jft )dt
+
+ +
= x( ) ( y (t ) ) exp( 2jft )dt d = x( )(TF ( y (t ) ))d
+ +
= x( )(TF ( y(t )))d = x( )Y ( f ) exp(2jf )d = X ( f )Y ( f )
SF et TF
CNAM ELE 103 D. Roviras 49
Proprits de la Transforme de Fourier :
Dmonstrations :
Conjugaison : changement de variable
Fonction relle : changement de variable
Dcalage temporel : changement de variable
Dcalage frquentiel : changement de variable
Multiplication et convolution
+
TF ( x(t ) ) = x(t ) exp( 2jft )dt posons u = t
+
TF ( x(t ) ) = x(u ) exp( 2jf [u + ])du
+
=
x(u ) exp( 2jfu )du. exp(2jf )
= TF (x(t ) ). exp(2jf )
SF et TF
CNAM ELE 103 D. Roviras 50
Transforme de Fourier usuelles (1/5)
x(t)= Rectangle = T (t ) T.sinc(f.T)= T.[sin(.f.T)/.f.T]
amplitude=1, largeur=T, centr sur 0
cos( 2. . fo. t ) 1
.[ ( f + fo ) + ( f fo ) ]
2
1 ( f )
(t ) 1
sin( 2. . fo.t ) j
.[ ( f + fo ) ( f fo ) ]
2
1 1
cos( 2 . . fo .t + 0 ) ( f + fo ). exp( j . 0 ) + ( f fo ). exp( j . 0 )
2 2
SF et TF
CNAM ELE 103 D. Roviras 51
Transforme de Fourier usuelles (2/5)
sin ( Bt ) X(f) = Rectangle de largeur B
x (t ) = B
( Bt ) entre - B/2 et + B/2 de hauteur 1
sin ( Bt )
2
X(f) = Triangle de largeur 2B
x ( t ) = B
( Bt ) entre - B et + B de hauteur 1
SF et TF
CNAM ELE 103 D. Roviras 52
Transforme de Fourier usuelles (3/5)
U(t) = Echelon unit 1 1
. ( f ) +
2 j .2. . f
Signe(t) 1
j. . f
exp(-a.t).U(t) 1
(U(t) = Echelon unit de temps) a + j .2. . f
2. a
exp(-a.|t|) a 2 + ( 2. . f ) 2
exp( . t 2 ) exp( . f 2 )
2. . fo
exp( a . t ).sin( 2. . fo. t ).U ( t )
( a + j.2. . f ) 2 + ( 2. . fo) 2
SF et TF
CNAM ELE 103 D. Roviras 53
Transforme de Fourier usuelles (4/5)
SF et TF
CNAM ELE 103 D. Roviras 54
Transforme de Fourier usuelles (5/5)
+
n
n =
X n . ( f
T
)
+
avec X n = A.sin c ( n / 2 )
2. A. T ( t n. T ) A
n =
Signal carr de largeur T/2, entre +A et -A et X n = 0 pour n pair ou nul
priodique de priode T
2. A
avec X n = pour n = 1, 5,....
n.
2. A
X n = pour n = 3, 7, ....
n.
+ +
n
A.
n =
( t n. T )
n =
X n . ( f )
T
Signal carr de largeur , entre 0 et +A et A.
avec X n = .sin c( n / T )
priodique de priode T T
+ 1 + n
. ( f )
(t n.T )
n =
T n = T
SF et TF
CNAM ELE 103 D. Roviras 55
Exemple de calcul de TF usuelles :
Calcul classique
Rectangle = T (t )
x(t)= Triangle Convolution de deux rectangles
damplitude=1, largeur=2.T,
centr sur 0
Calcul par TF-1 de (f-fo)
exp(2. . j. fo. t )
Dcomposition en
cos(2. . fo. t ) exponentielles
Dcomposition en
sin( 2. . fo.t ) exponentielles
Calcul classique ou limite du
(t )
rectangle
1 Limite du rectangle
t
T(t)
t
sin(fT )
x(t ) = x(t ). T (t ) X ( f ) = X ( f ) * T
fT
La fonction sin(fT)/fT stendant de moins linfini plus linfini, le support
spectral de X(f) est donc infini
SF et TF
CNAM ELE 103 D. Roviras 57
Allure gnrale dun signal et de sa TF
Si x(t) est trs pointu, sa TF sera trs tale
x(t)
t
X(f)
x(t)
t
X(f)
SF et TF
CNAM ELE 103 D. Roviras f 58
TF de signaux priodiques
Soit x(t) un signal priodique de priode T
x(t) est dcomposable en srie de Fourier
+
n
x(t ) = X n exp 2j t
n = T
+
n
X ( f ) = X n ( f )
n = T
Signal Bits
x(t) Numrisation
mettre
Signal Bits
Quantification
x(t) sur n bits mettre
Echantillonnage
Fe
Dmonstration
Echantillonnage Fe=1/Te
+ 1 + k
X e ( f ) = X ( f ) * TF (t k .Te) = X ( f ) * f
k = Te k = Te
1 + k
Xe( f ) = X f
Te k = Te
CNAM ELE 103 D. Roviras 62
SF et TF
Peigne de Dirac et chantillonnage
1 +
Xe( f ) = X ( f k .Fe )
Te k =
Fe>2.B k=0
X(f )
f
-B 0 +B
-Fe Fe
-B +B f
k=-1 X e ( f ) k=1 k=2
k=-2
Fe<2.B
f
k=0
-2Fe -Fe Fe 2Fe
Fe>2.B
Passage du signal chantillonn au signal temporel : voir TD
Echantillonnage
Fe > 2. fmax Fe
Choix de Fe ?
Signal
x(t) x(k.Te)
Echantillonnage
Fe > 2. fmax Fe
Choix de n ?
Li au bruit de quantification
Notion de QoS
Pas de
quantification
n bits par
x(k.Te) Quantification
chantillon
Signal Filtre Db
x(t) Anti CAN
Bits
repliement mettre
n Bits par
chantillon
Echantillonnage
Fe
Db = Fe . n
CNAM ELE 103 D. Roviras 70
SF et TF
Signaux nergie finie
1. Notion de corrlation
Signaux borns temporellement
Signaux de carr intgrable
Objectifs :
Dualit corrlation/spectre
Introduction de la notion de corrlation
Comment mesurer la ressemblance dun signal avec une version dcale de lui-mme ?
3
-1
-2
Signal x1(t)
-3
-4
0 100 200 300 400 500 600 700 800 900 1000
1
Signal x2(t)
0
-1
-2
-3
0 100 200 300 400 500 600 700 800 900 1000
2 2
0 0
-2
-2
-4
0 100 200 300 400 500 600 700 800 900 1000 630 632 634 636 638 640 642 644 646 648 650
4
4
x2(n) et x2(n-10)
2
2
0
0
-2
-2
-4
0 100 200 300 400 500 600 700 800 900 1000
632 634 636 638 640 642 644 646 648
+
Mesure de ressemblance entre x(t) et x(t + ) : x(t ) x(t + )dt
t t
Transforme de Fourier :
Densit
Spectrale dnergie
Sx2(f)
Sx1(f)
f f
+
K xx( ) = (t)x(t + )dt
*
Autocorrlation x
+
K xy( ) = x (t)y(t + )dt
*
Intercorrlation
Symtrie Hermitienne
K xx( ) = K xx* ( ) K xy( ) = K *yx ( )
+ 2
K xx( 0 ) = x(t)
dt = Energie de x(t)
Ingalit de Schwartz 2
K xx() K xx( 0) K xy() K xx( 0) K yy( 0)
Si x(t) est rel, Kxx() est une fonction relle et Si x(t) et y(t) sont rels, Kxx()
paire est relle
K xx( ) = x * ( ) * x( ) K xy( ) = x * ( ) * y ( )
(K xx( ) )' = K x x( ) = K xx ( )
' '
(K xy( ) ) = K x ' y( ) = K xy ' ( )
'
Rappel de + 2 + +
[ ]
X ( f ) = X ( f ) X ( f ) = TF x * ( t ) .TF [x(t )] = TF x * ( t ) * x(t )
2 *
[ ]
[ ]
X ( f ) = TF x* (t ) * x(t ) = TF [K xx (t )]
2
+
S xx ( f ) = S x ( f ) = K xx (t ) exp(2jft ) dt
f2 +f2
Energie entre f1 et f2 = S
f1
xx ( f )df + S
+ f1
xx ( f )df
Units:
Si x(t) est en Volts, Kxx(t) est en V2s et Sx(f) en V2s/Hz soit en V2s2
exp(+ a ) exp(a )
K xx( ) = si < 0 K xx( ) =
2a 2a
1 1 1 1 2a
= + =
[
2a [a - 2jf ] [ a - 2jf ] 2a a 2 + (2f )2 ] Sx ( f ) = 2
1
a + (2f )
2
0 1/a t
Faire varier a,
1/2a Kxx()
Consquences ?
0 1/a
Sx(f)
1/a2
1/(2a2)
0 a/(2.) f
T + T
T / 2
Signaux physiquement ralisables
Signaux priodiques
+T / 2
1
K xx( ) = lim (t)x(t + )dt
*
Autocorrlation x
T + T
T / 2
+T / 2
1
Intercorrlation K xy( ) = lim (t)y(t + )dt
*
x
T + T
T / 2
Units:
Si x(t) est en Volts, Kxx(t) est en V2
Si x(t) est en Unit-arbitraire, Kxx(t) est en (Unit-arbitraire)2
Autocorrlation Intercorrlation
Symtrie Hermitienne
K xx( ) = K xx* ( ) K xy( ) = K *yx ( )
+T / 2 2
1
K xx( 0 ) = lim
T + T x(t)
T / 2
dt = Puissance de x(t)
Ingalit de Schwartz 2
K xx() K xx( 0) K xy() K xx( 0) K yy( 0)
Si x(t) est rel, Kxx() est relle et paire Si x(t) et y(t) sont rels, Kxy()
est relle
(K xx( ) )' = K x x( ) = K xx ( )
' ' (K xy( ) ) = K x ' y( ) = K xy ' ( )
'
Puissance entre f1 et f2 = S
f1
xx ( f )df + S
+f1
xx ( f )df
Units:
Si x(t) est en Volts, Sx(f) est en V2/Hz
Si x(t) est en Unit-arbitraire, Sx(f) est en (Unit-arbitraire)2/Hz
CNAM ELE 103 D. Roviras 85
SPF
Proprits de la Densit Spectrale de Puissance
+
K xx ( ) = TF 1
[S x ( f )] = S x ( f ) exp(+2jf )df
+
K xx (0) = S
x ( f )df
x(t)=A.cos(2fot+)
+T / 2 + To / 2
1 1
K xx( ) = lim x (t)x(t + )dt = (t)x(t + )dt avec To = 1/f o
* *
x
T + T To
T / 2 To / 2
+ To / 2
A2
K xx( ) =
To cos(2f t + ) cos(2f
To / 2
o o (t + ) + )dt
1 1
Rappel : cos(a) cos(b) = cos(a + b) + cos(a b)
2 2
+ To / 2
A2
K xx( ) =
2To {cos(4f t + 2 ) + cos(2f )}dt
To / 2
o o
+ To / 2 + To / 2
A2 A2 A2
K xx( ) = cos( 4f o t + 2 )dt + cos( 2f o )dt = cos( 2f o )
2To To / 2
2To To / 2
2
A2
K xx( ) = cos(2f o ) Vrifier les proprits de la fonction
dautocorrlation de signaux puissance finie
2
CNAM ELE 103 D. Roviras 87
SPF
Exemple de calcul de fonction dautocorrlation et de DSP
A2
K xx( ) = cos(2f o )
2
A2
K xx( 0) = : Puissance d' un cosinus d' amplitude A et de phase quelconque
2
A2 A2 A2
S x(f ) = TF cos(2f o ) = ( f + fo ) + ( f fo )
2 4 4
x(t)=A.exp(2jfot)
+ To / 2 2 + To / 2
1 A
K xx( ) = (t)x(t + )dt = exp(2jf t ) exp(2jf (t + ))dt
*
x o o
To To / 2
To To / 2
2 + To / 2
A
K xx( ) = exp(2jf )dt = A exp( 2jf o )
2
o
To To / 2
K xx( ) = A exp(2jf o ) S x (f ) = A ( f f o )
2 2
n =
n
T n = T
Vrifier quen intgrant Sx(f) on retrouve Parseval, ou bien que la fonction dautocorrlation
en zro est bien la puissance
+T / 2 +T / 2 +
1 1 1
K xx ( ) = x (t)x(t + )dt = + = T (t)x(t + )dt
* * *
x T (t)x(t )dt x
T T / 2
T T / 2
T
+ +
1 1
K xx ( ) = xT (t)x(t + )dt = xT ( u)x( u)du avec u = -t
* *
T T
[ ] *
k = +
1 * 1
K xx ( ) = xT ( u) * x(u ) (u = ) = xT ( u) * xT (u ) * (u kT ) (u = )
T T k =
{ }
k = + k = +
1 * 1
K xx ( ) =
T
x
T
( u) * xT (u ) *
k =
(u kT )
(u = ) = K
T xT
(u ) *
k =
(u kT ) (u = )
k = + k = +
1 1 k
K xx ( ) = K xT ( ) * ( kT ) S x ( f ) = 2 S xT ( f ). ( f )
T k = T k = T
CNAM ELE 103 D. Roviras 91
SPF
Exemple de fonction dautocorrlation et DSP de signaux priodiques
k = + k = +
1 1 k
K xx ( ) = K xT ( ) * ( kT ) S x ( f ) = 2 S xT ( f ). ( f )
T k = T k = T
TF
KxT() SxT(f)
xT(t)
0 T t -T 0 T 0 f
x(t) Kx()
TF Sx(f)
0 t -T 0 T 0 f
-T T
-1/T 1/T
Notion de probabilit:
Nombre tirages avec 4
Tirage dun d six faces : Proba(4) = lim
nombre de tirages Nombre total de tirages
Probabilit :
Ensemble des rsultats de lexprience alatoire :
= partition de lensemble
Probabilit : Application de vers lensemble [0,1]
={1,2,3,4,5,6}
[0,1]
Calculer P(1), P(2), P(3), P(4), P(5), P(6)
Calculer P(1,2,3) = P(tirage = 1 ou 2 ou 3)
, P()=1
AI B
B, P(B)
A, P(A)
P(A I B) = P( A / B )P(B) = P( B / A) P( A)
B vrifi
A
Surface ( A B )
P(B / A) =
Surface ( A)
Surface ( A B ) Surface ( A B ) Surface ( A)
P(A B ) = = .
Surface () Surface ( A) Surface ()
P(A B ) = P(B / A) / P ( A)
CNAM ELE 103 D. Roviras 97
Proprits gnrales : formule des probabilits totales
Probabilits et vnements disjoints et complets:
A
A
B
AIB AIB
P(B) = P( A I B ) + P( A I B)
Ai I Aj = vide si i j N
N P(B) = P(B / Ai )P ( Ai )
U Ai =
i =1
i =1
Probabilits conditionnelles
P( A / B ) =
rgle de Bayes : P( B / A) P ( A)
P( B)
Exemple : Transmission de 0 et de 1 travers un canal de transmission
Emission : P(0)=0.9 P(1)=0.1
p
0 0
Canal binaire symtrique
p = Probabilit de ne pas faire une erreur
1-p = Probabilit derreur = 0.1 1-p
1 p 1
Calculer P(0 reu/1 mis)= P(dtecter 0/1mis)
Calculer P(1 reu/ 1 mis)
Calculer P(0 mis/ 1 reu)
Calculer P(1 mis/ 1 reu)
Calculer P(0 mis/ 0 reu) et P(1 mis/ 0 reu)
CNAM ELE 103 D. Roviras 99
Probabilits
Proprits gnrales : formule de Bayes
0.5
(50%) TR2 CD2
0.1 0.05
0.2 0.2
(20%) TR3 CD3
0.7
P(TR 2 / CD1) =
P(CD1 / TR 2) P(TR 2) 0.3 . 0.5
=
( Bayes ) P (CD1) P(CD1)
3
P(CD1) =
(Probabilits totales)
P(CD1 / TRi) P(TRi) = 0.8.0.3 + 0.3.0.5 + 0.1.0.2 = 0.41
i =1
P(TR 2 / CD1) =
0.3 . 0.5
= 0.3659
0.41
P(TR 2 / CD 2 ) =
P(CD 2 / TR 2) P(TR 2) 0.5 . 0.5
=
( Bayes ) P(CD 2) P(CD 2)
3
P(CD2) =
(Probabilits totales)
P(CD2 / TRi) P(TRi) = 0.15.0.3 + 0.5.0.5 + 0.2.0.2 = 0.335
i =1
P(TR 2 / CD 2 ) =
0.5 . 0.5
= 0.7463
0.335
P(TR 2 / CD3) =
P(CD3 / TR 2) P(TR 2) 0.2 . 0.5
=
( Bayes ) P(CD3) P(CD3)
3
P(CD3) =
(Probabilits totales)
P(CD3 / TRi) P(TRi) = 0.05.0.3 + 0.2.0.5 + 0.7.0.2 = 0.255
i =1
P(TR 2 / CD3) =
0.2 . 0.5
= 0.3922
0.255
On peut vrif ier que :
3
P(TR2) =
(Probabili ts totales)
P (TR 2 / CDi ) P (CDi ) =0.3659.0.41 + 0.7463.0.335 + 0.3922.0.255 = 0.5
i =1
P(TR1 / CD3) =
P (CD3 / TR1) P (TR1) 0.05 . 0.3
= = 0.0588
( Bayes ) P(CD3) 0.255
Une variable alatoire sera note dans le cours par une majuscule.
x(t)
Vo + Suite de valeurs
VA : X
n(t)
VA continue :
Proba(X=1.025083947)=0
La notion de probabilit dapparition dune valeur na pas de sens
Fonction de rpartition de la VA X :
Cest une fonction note FX(u) FX (uo ) = P( X < uo )
FX(u) est non dcroissante et comprise entre 0 et 1
FX ( ) = 0 FX ( + ) = 1
Sur les deux exemples prcdents, vrifier les proprits de la fonction de rpartition
Programme Matlab : ELE103_VA_1_et_2dim.m
Nombre denfants
Pas 4 enfants
denfant 1 enfant 2 enfants 3 enfants ou plus Total
Un seul adulte
Un homme
- actif
occup 18 262 104 367 45 888 13 878 4 278 186 673
- autre 58 693 32 621 9 414 3 273 2 024 106 025
Une femme
- active
occupe 81 176 499 239 259 611 70 433 17 578 928 037
- autre 332 807 223 046 118 586 56 359 33 066 763 864
Un homme et une femme
Deux actifs
occups 1 911 611 1 754 773 1 856 785 569 551 114 582 6 207 302
Un seul actif
occup
- lhomme 803 818 617 879 737 205 437 348 204 825 2 801 075
- la femme 577 933 187 951 114 714 43 240 17 790 941 628
Pas dactifs
occups 3 708 032 195 983 113 056 73 897 71 210 4 162 178
Total 7 492 332 3 615 859 3 255 259 1 267 979 465 353 16 096 782
FX (uo ) = P( X < uo )
FX(u) est non dcroissante et comprise entre 0 et 1
FX ( ) = 0 FX ( + ) = 1
P( x0 X < x1 ) = FX ( x1 ) FX ( x0 )
1
FX(u)
X : VA discrte
x1 x2 x3 x4 x5 u
1
FX(u)
X : VA continue
FX ( x)
uo
FX (uo ) = P ( X < uo ) = dx
x
uo
= [FX ( x)] = FX (uo ) FX () = FX (uo )
FX (x)
p X (x) =
x
FX (x)
p X (x) = pX(u)
X : VA continue
x
pX(x) est une fonction positive ou nulle u
+
pX(u) X : VA discrte
p
X (x)dx = 1
uo x1 x2 x3 u
P ( X < uo ) =
p
X (x)dx = FX (uo )
x1
P( x0 X < x1 ) = p X (x)dx = FX ( x1 ) FX ( x0 )
x0
Probabilits
CNAM ELE 103 D. Roviras 110
Variables alatoires discrtes et continues
Exemple de densit de probabilit de VA :
X est uniformment rpartie entre -1 et +1. Cela veut dire que la probabilit est
identique quelque soit x sur -1/+1
Calculer la densit de probabilit de cette VA continue
Tracer la densit de probabilit
Vrifier les proprits de la densit de probabilit
Calculer et tracer la fonction de rpartition de la VA X
Calculer P(0<X<0.5)
Calculer P(X=0.5)
pX(u)
uo
u
1 (u uo )
2 uo=moyenne
p X (u) = exp
= cart type
2 2 2
2 = variance
0.35
0.25
Signal 1 : Gaussienne de moyenne 0 et
dcart type sigma = 1
0.2
0.05
0
-5 0 5 10
10
-5
0 500 1000 1500 2000 2500
+ + 1 (u uo )2
P( X > S ) = p X (u ) du = exp du
S S
2 2 2
u uo
Changement de variable : v =
+ 1 v2
P( X > S ) = S uo exp dv. Pas de primitive:
2 2 Table
v S uo
2
+ 1
P ( X > S ) = S u o exp dv = Q
2 2
CNAM ELE 103 D. Roviras 115
Probabilits
Variables alatoires discrtes et continues
Loi Normale ou Gaussienne
Q(z)
p( z ) =
1 1 2
exp z
3.17 e-5
2. 2
+ 9.87 e-10
Q( z ) = p(u). du
z
1.28 e-12
0 2 4 6 7 z
CNAM ELE 103 D. Roviras 116
Probabilits
Variables alatoires discrtes et continues
+
Loi Normale rduite 1 1
Q(z) = z 2 .
exp u 2 .du
2
0 z u
z<0
Q(z) (non tabule)
z 0 u
+ z
1 1 1 1
Q(z) =
z 2 .
exp u 2 .du = 1
2
2 .
exp u 2 .du =
2
+
1 1
= 1
z 2 .
exp u 2 .du = 1 Q ( z )
2
CNAM ELE 103 D. Roviras 117
Probabilits
Variables alatoires discrtes et continues
Loi Normale ou Gaussienne
1 z
Q ( z ) = .erfc
2 2
Xmax
C3
C2
C1
Xmin
1 2 3 4 5 . N
Histogramme d'une Gaussienne avec 1000 points et 5 classes Histogramme d'une Gaussienne avec 1000 points et 100 classes
500 35
450
30
400
350 25
300
20
250
200 15
150
10
100
50 5
0
-4 -3 -2 -1 0 1 2 3 4 0
-4 -3 -2 -1 0 1 2 3 4
3500
6
3000
5
2500
4
2000
3
1500
2
1000
1
500
0 0
-5 -4 -3 -2 -1 0 1 2 3 4 5 -5 -4 -3 -2 -1 0 1 2 3 4
Petit Jeu :
Pendant 50% du temps je gagne 1
Pendant 50% du temps je gagne -1.5 (en fait cest une perte de 1.5)
Combien vais-je gagner en moyenne ?
+
Moment d ' ordre k = E ( X k ) = u k . p X (u )du
+
E (Y ) = E ( f ( X ) ) = f (u ). p X (u )du
+
E (Y ) = E ( f ( X ) ) = f (u ). p X (u )du
1 carat 1/carat = 1
100 carats 100/carat = 10000
Je gagne en moyenne pour chaque recherche : 0.9 x 1 + 0.1 x 10000 = 1000.9
+
E (Y ) = E ( f ( X ) ) = f (u ). p X (u )du
Le Taux dErreur de Bit (TEB) lorsque le gain du canal est de 1 est donn par
lexpression suivante:
( )
TEB ( gain = 1) = Q 1 .SNR
2
SNR est le rapport Signal sur Bruit (Ps/Pn) en mission,
+
1 u2
Q( z ) =
2 z exp 2 du
Le Taux dErreur de Bit (TEB) lorsque le gain du canal est de g est donn par
lexpression suivante:
Fonction Q(z)
0
10
-1
10
-2
10
-3
10
-4
10
-5
10
-6
10
-7
10
1 1.5 2 2.5 3 3.5 4 4.5 5
CNAM ELE 103 D. Roviras 129
Probabilits
Exemple
Exemple :
g est une VA discrte qui prend trois valeurs : 0.5, 1 et 2
P(g=0.5)=0.25, P(g=1)=0.5, P(g=2)=0.25
Calculer le TEB pour g gal une constante de 1 et un SNR de 5dB
Calculer le TEB pour g qui suit la loi de la VA discrte et un SNR de 5dB
SNRdB=5dB SNRlin=3.162
( )
TEB ( g = 1) = Q (1) .SNRlin = Q 3.162 = Q(1.77 ) = 0.0377
(
TEB( g = 0.5) = Q (0.5) .SNRlin = Q 0.25 . 3.162 )
= Q(0.8891) = 0.187
SNRdB=5dB SNRlin=3.162 g=2
( ) (
TEB( g = 2) = Q 2 .SNRlin = Q 4 . 3.162
2
)
= Q(3.55) = 1.87 10 4
SNRdB=5dB SNRlin=3.162 g = VA discrte
3
TEB ( g variable) = E (TEB ( g ) ) = E ( fonction( g ) ) = P( g i ).TEB ( g i )
i =1
-1
10
-2
10
-3
10
-4
10
0 1 2 3 4 5 6 7 8 9 10
2
X (
Variance = = E ( X E ( X ) ) = E ( X X )
2
) ( 2
)
Ecart type = X = Variance
2
X ( )
= E ( X E ( X ) ) = E (X ) (E ( X ) )
2 2 2
Densit de probabilit
n
p X 1,..., Xn (u1 ,...., un ) = FX 1,..., Xn (u1 ,...., u n )
u1u2 ...un
Proprits
FX1,..,Xn(u1,,un) : fonction non dcroissante et positive ou
nulle
FX 1,..., Xn ( ,...., ) = 0 FX 1,..., Xn (+ ,....,+ ) = 1
pX1,,Xn(u1,,un) : fonction positive ou nulle
1
Fonction de rpartition :
1/2
1/2
1 1/4
0
0 1
X2
Densit de probabilit :
1
0 X1
0 1
CNAM ELE 103 D. Roviras 137
Probabilits
Variables alatoires multidimensionnelles
Exemple de VA multidimensionnelle:
Tracer la fonction de rpartition et la densit de
probabilit de la VA suivante:
VA X deux dimensions X1 et X2
X1 : tirage dun d six faces
X2 : tirage dune pice en pile ou face (0 ou 1)
Exemple de VA multidimensionnelle:
Tracer la fonction de rpartition et la densit de
probabilit de la VA suivante:
VA X deux dimensions X1 et X2
X1 : tirage dun d six faces
X2 : somme du tirage du premier d et dun second d
Exemple de VA multidimensionnelle:
Tracer la densit de probabilit de la VA suivante:
VA X deux dimensions X1 et X2
X1 : Gaussienne de moyenne nulle et de variance un
X2 : Gaussienne de moyenne nulle et de variance un
Exemple de VA multidimensionnelle:
Tracer la densit de probabilit de la VA suivante:
VA X deux dimensions X1 et X2
X1 : Gaussienne de moyenne nulle et de variance un
X2 : Somme de X1 et dune autre Gaussienne de moyenne nulle et de variance un
600
400
100
200
80
0 60
100
80 40
60
40 20
Programme Matlab : 20
0 0
ELE103_Illustration_VA_continues_2dim.m
Histogramme VA multidimensionnelle
2000
1500
1000
500
0
100
100
80
50 60
Programme Matlab : 40
ELE103_Illustration_VA_continues_2dim.m 20
0 0
Programme Matlab :
ELE103_Illustration_VA_continues_2dim.m
2500
2000
1500
1000
500
0
100
50 100
80
60
40
0 20
0
CNAM ELE 103 D. Roviras 142
Probabilits
Variables alatoires multidimensionnelles
Exemple de VA multidimensionnelle:
X1 : Gaussienne de moyenne nulle et de variance un
X2 : Somme de X1 et dune VA uniformment rpartie entre +1 et -1
Programme Matlab :
Histogramme VA multidimensionnelle ELE103_Illustration_VA_continues_2dim.m
2500
2000
1500
1000
500
0
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100
1000
500
0
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100
VA continues VA discrtes
+ n
p X 1 (u1 ) = p X 1, X 2 (u1 , u2 )du2 p X 1 (u1 ) = PX 1, X 2 (u1 , ui )
i =1
+ m
p X 2 (u2 ) = p X 1, X 2 (u1 , u2 )du1 p X 2 (u2 ) = PX 1, X 2 (u j , u2 )
j =1
Cxy
Coefficient de corrlation : XY =
XY
CNAM ELE 103 D. Roviras 146
Probabilits
Variables alatoires multidimensionnelles
Proprits :
RXY = C XY + X Y
2
X +Y = + + 2.C XY
2
X
2
Y
X
pY (v) = J .[ p X (u )] u = f 1 ( v ) = .[ p X (u )] u = f 1 ( v )
Y
Soit (X1,X2,,Xn) un vecteur de VA continue caractris par sa densit de
probabilit pX1X2Xn(u1,u2,un), alors le vecteur de VA continue
(Y1,Y2,,Yn)=f(X1,X2,,Xn) a une densit de probabilit donne par :
Exemples :
Soit X une VA continue uniformment rpartie entre [-1, +1]. Calculer la
densit de probabilit de 2.X+3
X
pY (v) = J .[ p X (u )] u = f 1 ( v ) = .[ p X (u )] u = f 1 ( v )
Y
X 1 1 1
Y = X 2 X = Y = Y 1/ 2 = 0.5 1/ 2 J = 0.5 1/ 2 = 0.5
Y Y v v
u v = u 2 u = + v et v
[ p X (u )] u = f 1 (v ) = [ p X (u )] u =+ v et v = 1 / 2 + 1 / 2 = 1
1
0.5 pour 0 v 1
pY (v) = v
0 ailleurs
X 1 Y1 X 1 X 1
X 2 Y 2 = f =
X 2 X 1 + X 2
X 1 X 1
1 2 1 2
Y 1
pY 1Y 2 (v1 , v2 ) = J .[ p X 1 X 2 (u1 , u 2 )] (u ,u ) = f 1 ( v ,v ) = det
X 2
[
Y 2 . p
X 2 X 1 X 2 1 2
(u , u ) ] ( u1 ,u 2 ) = f 1 ( v1 ,v2 )
Y 1 Y 2
u1 v1 u1 u1 v1
u v = u + u u = v v
2 2 1 2 2 2 1
X 1 X 1
Y 1 Y 2 = det 1 1 = 1 = 1
J = det
X 2 X 2 0 1
Y 1 Y 2
4 Histogramme VA X et Y
Histogramme VA X et Y x 10
8
6000 X VA unif rpartie entre [-1,+1]
X VA unif rpartie entre [-1,+1]
Y=2.X+3 Y=X.X
7
5000
6
4000
5
3000 4
3
2000
1000
1
0 0
-6 -4 -2 0 2 4 6 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
4 Histogramme VA X et Y 4
x 10 x 10 Histogramme VA X et Z=X.Y
18
X VA unif rpartie entre [-1,+1] 3.5 X VA unif rpartie entre [-1,+1]
Y=X.X.X Z=X.Y
16
3
14
2.5
12
10 2
8
1.5
6
1
4
2 0.5
0 0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
5000 5000
4000 4000
3000 3000
2000 2000
1000 1000
0 0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -3 -2 -1 0 1 2 3
Bruit rsultant de la somme des N signaux indpendants (b) et signal individuel (r) Histogramme de la somme des N signaux indpendants
3 12000
2 10000
1 8000
0 6000
-1 4000
-2 2000
-3 0
0 100 200 300 400 500 600 700 800 900 1000 -6 -4 -2 0 2 4 6
Bruit rsultant de la somme des N signaux indpendants et signal individuel (r) Histogramme de la somme des N signaux indpendants
6 8000
4 7000
6000
2
5000
0
4000
-2
3000
-4
2000
-6
1000
-8 0
0 100 200 300 400 500 600 700 800 900 1000 -10 -8 -6 -4 -2 0 2 4 6 8 10
6
5000
4
4000
2
0 3000
-2
2000
-4
1000
-6
-8 0
0 100 200 300 400 500 600 700 800 900 1000 -10 -8 -6 -4 -2 0 2 4 6 8 10