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