How to get the phase in function of time of the hht (Hilbert-Huang-transform?

16 次查看（过去 30 天）

thomas baranowski 2022-10-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1836248-how-to-get-the-phase-in-function-of-time-of-the-hht-hilbert-huang-transform

评论： thomas baranowski 2022-11-2

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.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

AH 2022-10-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1836248-how-to-get-the-phase-in-function-of-time-of-the-hht-hilbert-huang-transform#answer_1085398

在 MATLAB Online 中打开

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 个评论
显示 1更早的评论隐藏 1更早的评论

thomas baranowski 2022-11-2

thomas baranowski 2022-11-2

Here are is the signal and time files.

Best Regards

请先登录，再进行评论。

请先登录，再回答此问题。

类别

Signal Processing Signal Processing Toolbox Transforms, Correlation, and Modeling Transforms Hilbert and Walsh-Hadamard Transforms

在 Help Center 和 File Exchange 中查找有关 Hilbert and Walsh-Hadamard Transforms 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by