Discussion :

[Mc geek]

Bonjour tout le monde,

J'ai reçu un TP de Matlab à faire et je bloque sur un point :
On me demande de calculer la puissance d'un signal échantillonné non-périodique. On nous a fournit la formule suivante :

La puissance d'un signal échantillonné x(n) est estimée en chaque instant nTe par :
P(n) = (1/(2*K+1)) * SOMME(x(k)^2, k=n-K, n+K)
Avec (2K+1)*Te la durée de la fenêtre temporelle pour l'estimation.

Pour pouvoir faire des essais j'ai pris 1ms comme fenêtre temporelle afin d'avoir un signal quasi-stationnaire durant cette période, ce qui donne :
K=(1*10^(-3)-Te)/(2*Te) avec Te = 1/44100.
Mon problème, c'est que je ne vois pas comment utiliser la fonction sum() pour pouvoir lui donner l'intervalle demandé.
Je vous fournis mon code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
[y, Fs] = audioread('bonjour.wav');
Ts = 1/Fs;
K = ((1*10.^(-3)-Ts)/Ts);
puissanceStereo = (1/(2*K+1))*sum((y.^2));
puissanceMono = (puissanceStereo(1) + puissanceStereo(2))/2;
disp(puissanceStereo);
disp(puissanceMono);

Auriez-vous des idées ?
Merci d'avance,

NB : Je précise que je sais que l'on peut utiliser la fonction sum() avec des dimensions en paramètres, mais quand j'ai tenté de mettre n+K et n-K en paramètres, il est marqué qu'il faut que ce soit des valeurs entières.

[phryte]

Bonjour,

Peut être :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
[y, Fs] = audioread('Bonjour.wav');
Ts = 1/Fs;
K = fix(((1*10.^(-3)-Ts)/Ts));
N=length(y);
for n=K+1:N-K
    k=n-K:n+K;
    P(n)=(1/(2*K+1))*sum(y(k).^2);
end
plot(1:N-K,P)
grid

[Mc geek]

Oui effectivement ça marche, merci !
Juste pour savoir si j'ai bien compris, en fait le vecteur y est crée à chaque étape de la boucle for, c'est ça ?

[phryte]

Oui, sur un horizon glissant de 2k+1.