How to get the phase in function of time of the hht (Hilbert-H​uang-trans​form?

9 visualizzazioni (ultimi 30 giorni)
Hi,
we try to decompensate a signal with matlab with a emd (empirical mode decomposition) and with hht(Hilbert-Huang transform). From the IMFs, we try to compute the instantaneous phase IP. It is already possible with python:
https://emd.readthedocs.io/en/stable/emd_tutorials/00_quick_start/emd_tutorial_00_start_01_quicksift.html
IP, IF, IA = emd.spectra.frequency_transform(imf, sample_rate, 'hilbert')
Is it also possible with matlab? I have only found this article about emd: https://de.mathworks.com/help/signal/ref/emd.html
Thanks for your help.

Risposte (1)

AH
AH il 27 Ott 2022
Instantaneous phase can be obtained by taking the integral from the instantaneous frequency, where is the initial phase at time and is the instantaneous frequency. To compute the integral numerically, cumtrapz can be used.
load('sinusoidalSignalExampleData.mat','X','fs')
t = (0:length(X)-1)/fs;
[imf,residual,info] = emd(X,'Interpolation','pchip');
iph = cumtrapz(t,imf);
figure
subplot(211),plot(t,imf(:,3)),xlabel("Time [sec]"),ylabel("Inst. Freq."),grid("minor")
subplot(212),plot(t,iph(:,3)),,xlabel("Time [sec]"),ylabel("Inst. Phase."),grid("minor")
  3 Commenti

Accedi per commentare.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by