clear all
clc
close all
data_bon =xlsread('bon_VHA_gvd_15Hz_3.xlsx'); % Charger les données
%t = data_bon.data(:,1); % Extraction de la colonne temps
M =size(data_bon); %taille de matrice
t = data_bon(:,1); %lire la 1ere colonne
Dt= t(2)-t(1);
fe = 1/Dt ; % frequance d'echan
N = length(data_bon(:, 2)); %nbr d'echant
ta = N*Dt;
figure;
plot(t, data_bon(:, 2)); % Exemple pour une direction
xlim([0.5 1.5]); % Ajustement temporel
title('Signal Vibratoire en Déplacement');
xlabel('Temps (s)');
ylabel('Déplacement (mm)');
grid on;
signal = data_bon(:, 2);
RMS = sqrt(sum(signal.^2) / length(signal));
Peak = max(abs(signal));
CF = Peak / RMS;
mu = sum(signal) / length(signal); % Moyenne du signal
sigma = sqrt(sum((signal - mu).^2) / length(signal)); % Écart-type
Kurtosis = sum((signal - mu).^4) / (length(signal) * sigma^4);