Discussion :

[Caroline_1]

Bonjour,

Je fais un analyse fréquentielle à partir d'un signal de houle tel que :

function ana_freq(n,a,b)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Analyse fréquentielle
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 
% clear x
% clear t
%filename=(['result' num2str(n)]);
matrix=load(['result' num2str(n) '.lvm'],'-ascii');
x=0.11329.*matrix(:,1)-0.007166;    % Relation d'etalonnage
x=x-mean(x);
 
for i=1:length(x)-1
t(i)=(i-1)*1/32;
end
%%%%%%%%%% Spectre d'énergie
xb=x(a*32:b*32);
tb=t(a*32:b*32);
Y=fft(xb);
N=length(xb);
df=1/(max(tb)-min(tb));
f=(0:1:N/2)*df;
%%Affichage du spectre
figure(n)
hold on
loglog(f(2:N/2+1),abs(Y(2:N/2+1)))
title(['S(f) -', num2str(n)])
xlabel('f (Hz)')
ylabel('S(f)')
hold off

Et je voulais savoir comment récupérer mon amplitude à partir du spectre que je viens de tracer. Ce serait sympa si vous pouviez m'aider ! Merci

Caro

[phryte]

Bonjour.

comment récupérer mon amplitude à partir du spectre

Peut-être avec ifft.

[Caroline_1]

Non ifft ne serverait qu'à récupérer mon signal initial que j'ai déjà.

[biomanred]

Je comprends pas trop si tu cherches l'amplitude du spectre (le module) ou l'amplitude du signal temportel?

[Caroline_1]

Je veux récupérer l'amplitude de mon signal grâce à mon spectre, je sais que c'est possible via une formule que je ne connais pas.
Il y a une autre méthode que j'ai déjà faite pour obtenir l'amplitude de mon signal en prenant les passages par zéro en up crossing mais ce n'est pas celle là qui m'intéresse.