Discussion :

[najbet]

salut
comment peut-on calculer l'énergie d'un signal sous matlab
merci

[sango85]

Quelle est ta formule à coder (plusieurs formules d'énergies et de puissances) ? Et où se situent tes difficultés ?

[najbet]

En fait mon signal est un signal que j'ai acquis par un carte son, donc je recupere des donnée. pour l'energie je n'arrive pas a le faire sur matlab, la formule est bien

Energie c'est bien l'integral du module du signal au carré sur -infini à + l'infini

J'ai fait ça:

Code :

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

1
2
3
4
5
6
7
 m = abs(data1);
 n = square(m);
 
 subplot(223);
 c = 0:0.1:2.2E5;
 i = trapz(c,n);
 plot(i);

et il me dis que:

Error using ==> trapz at 59
LENGTH(X) must equal the length of the first non-singleton dimension of Y.

Error in ==> t at 18
i = trapz(c,n);

j'ai joins la figure du signal que je recupère
merci

[sango85]

Energie c'est bien l'integral du module du signal au carré sur -infini à + l'infini

Oui si tu es en continu, mais en échantillonné je ne pense pas.
Je peux simplement t'aiguiller en disant que tu peux calculer la puissance RMS de ton signal en faisant :

Code :

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

>> Prms = sqrt(mean(signal.^2));

Bon courage...

[najbet]

comment puis-je calculer la puissance RMS par bande de frequence?

[sango85]

Là, c'est un peu plus délicat. Il est nécessaire de calculer la FFT, mais le principe est le même (sauf qu'au lieu de moyenner, on somme...).

Code :

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

1
2
3
4
5
6
signal = ...;
Fech = ...;
dF = Fech/size(signal);
%% Calcul de ta FFT
FFT_signal = fft(signal,NFFT)/NFFT;
Prms_totale = sqrt(sum(FFT_signal.^2))

Et pour calculer la puissance dans une bande de fréquence, avec ton vecteur-fréquences, tu sélectionnes les points correspondant à ta bande puis tu réalises le même calcul sur ta sélection de points..

[najbet]

salut
meci

En fait je voudrais afficher sur un meme graphe la distribution de puissance par bande de frequence de mon signal
je sais pas trop si ce que je dis est comprehensible

[ancrou]

Parles tu de DSP (densité spectrale de puissance)?

lien Wiki