%%% CALCUL DE LA FFT DU SIGNAL %%%
clc
clear all
load ('D:\TD\Maintenance industrielle\TP_Diagnostic des machines électriques\TP1\TP1\sain20.mat');
tail_fft=      100000;         %modéfier cette valeur selon le nombre de points du signal
              
tail_fenetre=tail_fft;  %taille de la fenetre.

fe=10000 ;        %fréquence d'échantillonnage.

 %Le spectre du courant:
[B,F,T]=specgram(ia,tail_fft,fe,hanning(tail_fenetre),1);
 
ab=abs(B);            %L'amplitude de la FFT

ab_nor=ab/(max(max(ab))); %Normalisation par rapport au maximum:

 ab_nor_lo=20*log10(ab_nor);     % echelle logarithmique.
plot(ia);
figure(2)
 plot(F, ab_nor_lo,'b');hold on