Unité de la densité spectrale des énergies sur MATLAB

[invite0819340a]

Bonjour,

Je vous prie de m'aider afin de résoudre ce problème qui est le suivant.

Je veux effectuer pour mon signal temporel une Power Spectral Density sur Matlab pour connaître les énergies contenues dans mon signal (signal temporel); x: temps (s) et y: des élévations en cm.


Je vous prie de me dire si mon code est bon: le nombre total de mes données est 28798 dont la fréquence échantillonnage est de 0.016 Hz.


Aussi, en lisant la doc en faisant open psd j'ai vu que la PSD est en watt/Hz, alors que je veux l'avoir en cm2/Hz car mes (y) sont en cm en focntion de (t).

Je vous remercie d'avance.
code:

load tempor.txt
y=tempor(:,1);
x=tempor(:,2);
Y = fft(y,28798);
Pyy = Y.*conj(Y)/28798;
f = 0.016/28798*(0:14398);
plot(f,Pyy(1:14399))
title('Power spectral density')
xlabel('Frequency (Hz)')

-----