Estimation spectrale Méthodes non-paramétriques

Finalement, l’estimateur de Bartlett est:
ˆSB(ω) = 1
N
K−1
i=0
L−1
x(n + iL)exp(−jωn)
n=0
2
. (32)
Puisque les périodogrammes dans ˆ SB(ω) sont calculés
en utilisant des séquences de longueur L, larésolution
est:
rés
ˆSB(ω)
=0.89 2π
L
2π
=0.89K , (33)
N
qui est K fois plus grande (ou pire) que celle du
périodogramme.
Code en MATLAB
function Px = bart(x,K);
N = length(x);
L = floor(N/K);
Px = 0;
n1 = 1;
for i=1:K
Px = Px + periodogram(x(n1:n1+L-1))/K;
n1 = n1 + L;
end
INRS-EMT J. Benesty 15