HS matlab code:

y=load ('signal');

N=length(y);

fs=; %Sampling Frequency;

%---- Filtring ------------- %----------------------------

%Input fmin et fmax du filter

yf=fitlter(y,fmin,fmax,fs);

%---- Hilbert ------- %---------------------

H=hilbert(yf);

yh=abs(H);

nfft=1024;

p=abs(fft(yh,nfft));

Good luck

Thanks for your reply Hocine Bendjama

I have decomposed the main signal into some IMFs using EMD and now as I understand, I shoud write the code after filtering part? Also, write imf instead of yf?

Is it correct?

H=hilbert(imf);

yh=abs(H);

nfft=1024;

p=abs(fft(yh,nfft));

Yes, it is correct, but nfft=1024 depends on the representation and the length of the signal, in your case (imf)