0% ont trouvé ce document utile (0 vote)
28 vues1 page

Clear All

Le document présente un script MATLAB pour analyser des données vibratoires extraites d'un fichier Excel. Il calcule des paramètres statistiques tels que la moyenne, l'écart-type, la valeur de crête et la kurtosis du signal. Une visualisation du signal en fonction du temps est également incluse.

Transféré par

zakiben305
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
28 vues1 page

Clear All

Le document présente un script MATLAB pour analyser des données vibratoires extraites d'un fichier Excel. Il calcule des paramètres statistiques tels que la moyenne, l'écart-type, la valeur de crête et la kurtosis du signal. Une visualisation du signal en fonction du temps est également incluse.

Transféré par

zakiben305
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

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);

Vous aimerez peut-être aussi