hht

Преобразование Гильберта-Хуанга

свернуть все на странице

Синтаксис

hs = hht(imf)

hs = hht(imf,fs)

[hs,f,t] = hht(___)

[hs,f,t,imfinsf,imfinse] = hht(___)

[___] = hht(___,Name,Value)

hht(___)

hht(___,freqlocation)

Описание

пример

hs = hht(imf) возвращает Гильбертов спектр hs из сигнала, заданного внутренним режимом, функционирует imf. hs полезно для анализа сигналов, которые включают смесь сигналов, спектральное содержимое которых изменяется вовремя. Используйте hht выполнять Гильбертов спектральный анализ на сигналах идентифицировать локализованные функции.

пример

hs = hht(imf,fs) возвращает Гильбертов спектр hs из сигнала, произведенного на уровне fs.

пример

[hs,f,t] = hht(___) возвращает вектор частоты f и временной вектор t в дополнение к hs. Эти выходные аргументы могут использоваться с любым из предыдущих входных синтаксисов.

пример

[hs,f,t,imfinsf,imfinse] = hht(___) также возвращает мгновенные частоты imfinsf и мгновенные энергии imfinse из внутреннего режима функционирует для диагностики сигнала.

пример

[___] = hht(___,Name,Value) оценки параметры спектра Гильберта с дополнительными опциями, заданными одним или несколькими Name,Value парные аргументы.

пример

hht(___) без выходных аргументов строит Гильбертов спектр в окне текущей фигуры. Можно использовать этот синтаксис с любым из входных параметров в предыдущих синтаксисах.

hht(___,freqlocation) строит Гильбертов спектр с дополнительным freqlocation аргумент, чтобы задать местоположение оси частоты. Частота представлена на y - ось по умолчанию.

Примеры

свернуть все

Гильбертов спектр квадратичного щебета

Скрипт Open Live Script

Сгенерируйте Модулируемый гауссовым образом квадратичный щебет. Задайте частоту дискретизации 2 кГц и длительность сигнала 2 секунд.

fs = 2000;
t = 0:1/fs:2-1/fs;
q = chirp(t-2,4,1/2,6,'quadratic',100,'convex').*exp(-4*(t-1).^2);
plot(t,q)

Figure contains an axes object. The axes object contains an object of type line.

Используйте emd визуализировать внутренние функции режима (IMFs) и невязку.

emd(q)

Вычислите IMFs сигнала. Используйте 'Display' пара "имя-значение", чтобы вывести таблицу, показывающую количество отсеивания итераций, относительной погрешности и отсеивания, останавливает критерий каждого МВФ.

imf = emd(q,'Display',1);

Current IMF  |  #Sift Iter  |  Relative Tol  |  Stop Criterion Hit  
      1      |        2     |    0.0063952   |  SiftMaxRelativeTolerance
      2      |        2     |       0.1007   |  SiftMaxRelativeTolerance
      3      |        2     |      0.01189   |  SiftMaxRelativeTolerance
      4      |        2     |    0.0075124   |  SiftMaxRelativeTolerance
Decomposition stopped because the number of extrema in the residual signal is less than the 'MaxNumExtrema' value.

Используйте вычисленный IMFs, чтобы построить Гильбертов спектр квадратичного щебета. Ограничьте частотный диапазон от 0 Гц до 20 Гц.

hht(imf,fs,'FrequencyLimits',[0 20])

Figure contains an axes object. The axes object with title Hilbert Spectrum contains 4 objects of type patch.

Выполните эмпирическое разложение моды и визуализируйте гильбертов спектр сигнала

Скрипт Open Live Script

Загрузите и визуализируйте неустановившийся непрерывный сигнал, состоявший из синусоидальных волн с отличным изменением в частоте. Вибрация отбойного молотка и звук фейерверков являются примерами неустановившихся непрерывных сигналов. Сигнал производится на уровне fs.

load('sinusoidalSignalExampleData.mat','X','fs')
t = (0:length(X)-1)/fs;
plot(t,X)
xlabel('Time(s)')

Figure contains an axes object. The axes object contains an object of type line.

Смешанный сигнал содержит синусоидальные волны с различной амплитудой и значениями частоты.

Чтобы создать Гильбертов график спектра, вам нужны внутренние функции режима (IMFs) сигнала. Выполните эмпирическое разложение моды, чтобы вычислить IMFs и остаточные значения сигнала. Поскольку сигнал не является гладким, задайте 'pchip'как метод интерполяции.

[imf,residual,info] = emd(X,'Interpolation','pchip');

Таблица, сгенерированная в командном окне, показывает, что количество отсеивает итерации, относительную погрешность, и отсеять критерий остановки каждого сгенерировал МВФ. Эта информация также содержится в info. Можно скрыть таблицу путем добавления 'Display',0 пара значение-имя.

Создайте Гильбертов график спектра с помощью imf компоненты, полученные с помощью эмпирического разложения моды.

hht(imf,fs)

Figure contains an axes object. The axes object with title Hilbert Spectrum contains 9 objects of type patch.

Частота по сравнению с графиком временной зависимости является разреженным графиком с вертикальной цветной полосой, указывающей на мгновенную энергию в каждой точке в МВФ. График представляет мгновенный спектр частоты каждого компонента, анализируемого от исходного смешанного сигнала. Три IMFs появляются в графике с отличным изменением в частоте в 1 секунду.

Гильбертов спектр песни кита

Скрипт Open Live Script

Загрузите файл, который содержит аудиоданные от Тихоокеанского голубого кита, произведенного на уровне 4 кГц. Файл от библиотеки вокализаций животных, обеспеченных Программой исследований Биоакустики Корнелльского университета. Масштаб времени в данных сжат на коэффициент 10, чтобы повысить тангаж и выполнить более слышимые вызовы. Преобразуйте сигнал в расписание MATLAB® и постройте его. Четыре функции выделяются от шума в сигнале. Первое известно как трель, и другие три известны как стоны.

[w,fs] = audioread('bluewhale.wav');
whale = timetable(w,'SampleRate',fs);
stackedplot(whale);

Figure contains an object of type stackedplot.

Используйте emd визуализировать первые три внутренних функции режима (IMFs) и невязку.

emd(whale,'MaxNumIMF',3)

Вычислите первые три IMFs сигнала. Используйте 'Display' пара "имя-значение", чтобы вывести таблицу, показывающую количество отсеивания итераций, относительной погрешности и отсеивания, останавливает критерий каждого МВФ.

imf = emd(whale,'MaxNumIMF',3,'Display',1);

Current IMF  |  #Sift Iter  |  Relative Tol  |  Stop Criterion Hit  
      1      |        1     |      0.13523   |  SiftMaxRelativeTolerance
      2      |        2     |     0.030198   |  SiftMaxRelativeTolerance
      3      |        2     |      0.01908   |  SiftMaxRelativeTolerance
Decomposition stopped because maximum number of intrinsic mode functions was extracted.

Используйте вычисленный IMFs, чтобы построить Гильбертов спектр сигнала. Ограничьте частотный диапазон от 0 Гц до 1 400 Гц.

hht(imf,'FrequencyLimits',[0 1400])

Figure contains an axes object. The axes object with title Hilbert Spectrum contains 3 objects of type patch.

Вычислите Гильбертов спектр для той же области значений частот. Визуализируйте Гильбертовы спектры трели и стонов как сетчатый график.

[hs,f,t] = hht(imf,'FrequencyLimits',[0 1400]);

mesh(seconds(t),f,hs,'EdgeColor','none','FaceColor','interp')
xlabel('Time (s)')
ylabel('Frequency (Hz)')
zlabel('Instantaneous Energy')

Figure contains an axes object. The axes object contains an object of type surface.

Вычислите гильбертовы параметры спектра сигнала

Скрипт Open Live Script

load('sinusoidalSignalExampleData.mat','X','fs')
t = (0:length(X)-1)/fs;
plot(t,X)
xlabel('Time(s)')

Figure contains an axes object. The axes object contains an object of type line.

Смешанный сигнал содержит синусоидальные волны с различной амплитудой и значениями частоты.

Чтобы вычислить Гильбертовы параметры спектра, вам нужен IMFs сигнала. Выполните эмпирическое разложение моды, чтобы вычислить внутренние функции режима и остаточные значения сигнала. Поскольку сигнал не является гладким, задайте 'pchip' как метод интерполяции.

[imf,residual,info] = emd(X,'Interpolation','pchip');

Вычислите Гильбертовы параметры спектра: Гильбертов спектр hs, вектор частоты f, временной вектор t, мгновенная частота imfinsf, и мгновенная энергия imfinse.

[hs,f,t,imfinsf,imfinse] = hht(imf,fs);

Используйте вычисленные Гильбертовы параметры спектра для частотно-временного анализа и диагностики сигнала.

VMD многокомпонентного сигнала

Скрипт Open Live Script

Сгенерируйте многокомпонентный сигнал, состоящий из трех синусоид частот 2 Гц, 10 Гц и 30 Гц. Синусоиды производятся на уровне 1 кГц в течение 2 секунд. Встройте сигнал в белый Гауссов шум отклонения 0,01 ².

fs = 1e3;
t = 1:1/fs:2-1/fs;
x = cos(2*pi*2*t) + 2*cos(2*pi*10*t) + 4*cos(2*pi*30*t) + 0.01*randn(1,length(t));

Вычислите IMFs сигнала с шумом и визуализируйте их в 3-D графике.

imf = vmd(x);
[p,q] = ndgrid(t,1:size(imf,2));
plot3(p,q,imf)
grid on
xlabel('Time Values')
ylabel('Mode Number')
zlabel('Mode Amplitude')

Figure contains an axes object. The axes object contains 5 objects of type line.

Используйте вычисленный IMFs, чтобы построить Гильбертов спектр многокомпонентного сигнала. Ограничьте частотный диапазон [0, 40] Гц.

hht(imf,fs,'FrequencyLimits',[0,40])

Figure contains an axes object. The axes object with title Hilbert Spectrum contains 5 objects of type patch.

Вычислите гильбертов спектр сигнала вибрации

Скрипт Open Live Script

Симулируйте сигнал вибрации от поврежденного подшипника. Вычислите Гильбертов спектр этого сигнала и ищите дефекты.

Средний диаметр подшипника 12 см, в подшипнике восемь тел качения. Каждое тело качения имеет диаметр 2 см. Внешнее кольцо остается стационарным, внутренне кольцо вращается со скоростью 25 оборотов в секунду. Акселерометр производит измерения колебаний подшипника с частотой дискретизации 10 кГц.

fs = 10000;
f0 = 25;
n = 8;
d = 0.02;
p = 0.12;

Сигнал вибрации от здорового подшипника включает несколько порядков ведущей частоты.

t = 0:1/fs:10-1/fs;
yHealthy = [1 0.5 0.2 0.1 0.05]*sin(2*pi*f0*[1 2 3 4 5]'.*t)/5;

Резонанс взволнован в вибрации подшипника на полпути посредством процесса измерения.

yHealthy = (1+1./(1+linspace(-10,10,length(yHealthy)).^4)).*yHealthy;

Резонанс вводит дефект во внешней гонке подшипника, который приводит к прогрессивному износу. Дефект вызывает ряд ударов, которые повторяются на частоте передачи мяча внешняя гонка (BPFO) подшипника:

$B P F O = \frac{1}{2} n f_{0} [1 - \frac{d}{p} \cos θ],$

где $f_{0}$ пропускная способность, $n$ количество прокручивающихся элементов, $d$ диаметр прокручивающихся элементов, $p$ диаметр тангажа подшипника, и $θ$ угол контакта подшипника. Примите угол контакта 15 ° и вычислите BPFO.

ca = 15;
bpfo = n*f0/2*(1-d/p*cosd(ca));

Используйте pulstran (Signal Processing Toolbox) функция, чтобы смоделировать удары как периодическое обучается синусоид с 5 миллисекундами. Каждая синусоида на 3 кГц является оконной окном с плоской вершиной. Используйте закон о степени, чтобы ввести прогрессивный износ в сигнале вибрации подшипника.

fImpact = 3000;
tImpact = 0:1/fs:5e-3-1/fs;
wImpact = flattopwin(length(tImpact))'/10;
xImpact = sin(2*pi*fImpact*tImpact).*wImpact;

tx = 0:1/bpfo:t(end);
tx = [tx; 1.3.^tx-2];

nWear = 49000;
nSamples = 100000;
yImpact = pulstran(t,tx',xImpact,fs)/5;
yImpact = [zeros(1,nWear) yImpact(1,(nWear+1):nSamples)];

Сгенерируйте сигнал вибрации BPFO путем добавления ударов на здоровый сигнал подшипника. Постройте сигнал и выберите 0,3 вторых интервала, запускающиеся в 5,0 секунд.

yBPFO = yImpact + yHealthy;

xLimLeft = 5.0;
xLimRight = 5.3;
yMin = -0.6;
yMax = 0.6;

plot(t,yBPFO)

hold on
[limLeft,limRight] = meshgrid([xLimLeft xLimRight],[yMin yMax]);
plot(limLeft,limRight,'--')
hold off

Figure contains an axes object. The axes object contains 3 objects of type line.

Увеличьте масштаб выбранного интервала, чтобы визуализировать эффект ударов.

xlim([xLimLeft xLimRight])

Figure contains an axes object. The axes object contains 3 objects of type line.

Добавьте белый Гауссов шум в сигналы. Задайте шумовое отклонение $1 / 15 0^{2}$ .

rn = 150;
yGood = yHealthy + randn(size(yHealthy))/rn;
yBad = yBPFO + randn(size(yHealthy))/rn;

plot(t,yGood,t,yBad)
xlim([xLimLeft xLimRight])
legend('Healthy','Damaged')

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Healthy, Damaged.

Используйте emd (Signal Processing Toolbox), чтобы выполнить эмпирическое разложение моды здорового сигнала подшипника. Вычислите первые пять внутренних функций режима (IMFs). Используйте 'Display' аргумент значения имени, чтобы вывести таблицу, показывающую количество отсеивания итераций, относительной погрешности и отсеивания, останавливает критерий каждого МВФ.

imfGood = emd(yGood,'MaxNumIMF',5,'Display',1);

Current IMF  |  #Sift Iter  |  Relative Tol  |  Stop Criterion Hit  
      1      |        3     |     0.017132   |  SiftMaxRelativeTolerance
      2      |        3     |      0.12694   |  SiftMaxRelativeTolerance
      3      |        6     |      0.14582   |  SiftMaxRelativeTolerance
      4      |        1     |     0.011082   |  SiftMaxRelativeTolerance
      5      |        2     |      0.03463   |  SiftMaxRelativeTolerance
Decomposition stopped because maximum number of intrinsic mode functions was extracted.

Используйте emd без выходных аргументов, чтобы визуализировать первые три IMFs и невязку.

emd(yGood,'MaxNumIMF',5)

Вычислите и визуализируйте IMFs дефектного сигнала подшипника. Первый эмпирический режим показывает высокочастотные удары. Этот высокочастотный режим увеличения энергии как износ прогрессирует.

imfBad = emd(yBad,'MaxNumIMF',5,'Display',1);

Current IMF  |  #Sift Iter  |  Relative Tol  |  Stop Criterion Hit  
      1      |        2     |     0.041274   |  SiftMaxRelativeTolerance
      2      |        3     |      0.16695   |  SiftMaxRelativeTolerance
      3      |        3     |      0.18428   |  SiftMaxRelativeTolerance
      4      |        1     |     0.037177   |  SiftMaxRelativeTolerance
      5      |        2     |     0.095861   |  SiftMaxRelativeTolerance
Decomposition stopped because maximum number of intrinsic mode functions was extracted.

emd(yBad,'MaxNumIMF',5)

Постройте Гильбертов спектр первого эмпирического режима дефектного сигнала подшипника. Первый режим получает эффект высокочастотных ударов. Энергия увеличений ударов как износ подшипников прогрессирует.

figure
hht(imfBad(:,1),fs)

Figure contains an axes object. The axes object with title Hilbert Spectrum contains an object of type patch.

Гильбертов спектр третьего режима показывает резонанс в сигнале вибрации. Ограничьте частотный диапазон от 0 Гц до 100 Гц.

hht(imfBad(:,3),fs,'FrequencyLimits',[0 100])

Figure contains an axes object. The axes object with title Hilbert Spectrum contains an object of type patch.

Для сравнения постройте Гильбертовы спектры первых и третьих режимов здорового сигнала подшипника.

subplot(2,1,1)
hht(imfGood(:,1),fs)
subplot(2,1,2)
hht(imfGood(:,3),fs,'FrequencyLimits',[0 100])

Figure contains 2 axes objects. Axes object 1 with title Hilbert Spectrum contains an object of type patch. Axes object 2 with title Hilbert Spectrum contains an object of type patch.

Входные параметры

свернуть все

`imf` — Внутренняя функция режима
матрица | расписание

Внутренняя функция режима в виде матрицы или расписания. imf любой сигнал, конверт которого симметричен относительно нуля и чьи количества экстремальных и нулевых пересечений отличаются самое большее один. emd используется, чтобы анализировать и упростить сложные сигналы в конечное число внутренних функций режима, требуемых выполнять Гильбертов спектральный анализ.

hht обработки каждый столбец в imf как внутренняя функция режима. Для получения дополнительной информации о вычислении imf, смотрите emd.

`fs` — Частота дискретизации
2π (значение по умолчанию) | положительная скалярная величина

Частота дискретизации в виде положительной скалярной величины. Если fs не предоставляется, нормированная частота 2π используется для расчета Гильбертов спектр. Если imf задан как расписание, частота дискретизации выведена из него.

`freqlocation` — Местоположение оси частоты на графике
`'yaxis'` (значение по умолчанию) | `'xaxis'`

Местоположение оси частоты на графике в виде 'yaxis' или 'xaxis'. Чтобы отобразить данные о частоте по y - ось или x - ось графика, задайте freqlocation как 'yaxis' или 'xaxis' соответственно.

Аргументы name-value

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'FrequencyResolution',1

`FrequencyLimits` — Пределы частоты, чтобы вычислить Гильбертов спектр
`[0, fs/2]` (значение по умолчанию) | 1 2 вектор с целочисленным знаком

Пределы частоты, чтобы вычислить Гильбертов спектр в виде разделенной запятой пары, состоящей из 'FrequencyLimits'и 1 2 вектор с целочисленным знаком. FrequencyLimits задан в Гц.

`FrequencyResolution` — Разрешение частоты, чтобы дискретизировать частотный диапазон
(f_high-f_low)/100 (значение по умолчанию) | положительная скалярная величина

Разрешение частоты, чтобы дискретизировать пределы частоты в виде разделенной запятой пары, состоящей из 'FrequencyResolution'и положительная скалярная величина.

Задайте FrequencyResolution в Гц. Если 'FrequencyResolution' не задан, значение (_fhigh-_flow)/100 выведено из FrequencyLimits. Здесь, _fhigh является верхним пределом FrequencyLimits и _flow является нижним пределом.

`MinThreshold` — Минимальное пороговое значение Гильбертова спектра
`-inf` (значение по умолчанию) | скаляр

Минимальное пороговое значение Гильбертова спектра в виде разделенной запятой пары, состоящей из 'MinThreshold'и скаляр.

MinThreshold элементы наборов hs к 0, когда соответствующие элементы $10 \log_{10} (h s)$ меньше MinThreshold.

Выходные аргументы

свернуть все

`hs` — Гильбертов спектр сигнала
разреженная матрица

Гильбертов спектр сигнала, возвращенного как разреженная матрица. Используйте hs для частотно-временного анализа и идентифицировать локализованные функции в сигнале.

`f` — Значения частоты
вектор

Значения частоты сигнала, возвращенного как вектор. hht использует вектор частоты f и временной вектор t создать Гильбертов график спектра.

Математически, f обозначается как: f = _flow: _fres: _fhigh, где _fres является разрешением частоты.

`t` — Временные стоимости
вектор | `duration` массив

Временные стоимости сигнала, возвращенного как вектор или a duration массив. hht использует временной вектор t и вектор частоты f создать Гильбертов график спектра.

t возвращен как:

Массив, если imf задан как массив.
A duration массив, если imf задан как однородно произведенное расписание.

`imfinsf` — Мгновенная частота каждого МВФ
вектор | матрица | расписание

Мгновенная частота каждого МВФ, возвращенного как вектор, матрица или расписание.

imfinsf имеет одинаковое число столбцов как imf и возвращен как:

Вектор, если imf задан как вектор.
Матрица A, если imf задан как матрица.
Расписание, если imf задан как однородно произведенное расписание.

`imfinse` — Мгновенная энергия каждого МВФ
вектор | матрица | расписание

Мгновенная энергия каждого МВФ, возвращенного как вектор, матрица или расписание.

imfinse имеет одинаковое число столбцов как imf и возвращен как:

Вектор, если imf задан как вектор.
Матрица A, если imf задан как матрица.
Расписание, если imf задан как однородно произведенное расписание.

Алгоритмы

Преобразование Гильберта-Хуанга полезно для выполнения частотно-временного анализа неустановившихся и нелинейных данных. Процедура Хилбэрт-Хуана состоит из следующих шагов:

emd или vmd разлагает набор данных x на конечное число внутренних функций режима.
Для каждой внутренней функции режима, _xi, функции hht:
1. Использование hilbert (Signal Processing Toolbox), чтобы вычислить аналитический сигнал, $z_{i} (t) = x_{i} (t) + j H {x_{i} (t)}$ , где H {_xi} является преобразованием Гильберта _xi.
2. Экспрессы _zi как $z_{i} (t) = a_{i} (t) e^{j θ_{i} (t)}$ , где _ai (t) является мгновенной амплитудой и $θ_{i} (t)$ мгновенная фаза.
3. Вычисляет мгновенную энергию, $| a_{i} (t) |^{2}$ , и мгновенная частота, $ω_{i} (t) \equiv d θ_{i} (t) / d t$ . Если дали частота дискретизации, hht преобразует $ω_{i} (t)$ к частоте в Гц.
4. Выводит мгновенную энергию в imfinse и мгновенная частота в imfinsf.
Когда названо без выходных аргументов, hht строит энергию сигнала в зависимости от времени и частоты, с цветом, пропорциональным амплитуде.

Ссылки

[1] Хуан, Норден Э и Сэмюэль С П Шен. Преобразование Гильберта-Хуанга и Его Приложения. 2-е издание 16 редактора. Междисциплинарные Математические Науки. НАУЧНЫЙ МИР, 2014. https://doi.org/10.1142/8804.

[2] Хуан, Норден Э., Жаохуа Ву, Стивен Р. Лонг, Кеннет К. Арнольд, Ксиэньяо Чен и Карин Блэнк. “НА МГНОВЕННОЙ ЧАСТОТЕ”. Усовершенствования в Адаптивном Анализе данных 01, № 02 (апрель 2009): 177–229. https://doi.org/10.1142/S1793536909000096.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Указания и ограничения по применению:

Аргументы, заданные с помощью пар "имя-значение", должны быть константами времени компиляции.
Расписания не поддерживаются для генерации кода.

Смотрите также

emd | vmd

Темы

Галерея частоты времени

Введенный в R2018a

Документация

hht

Синтаксис

Описание

Примеры

Гильбертов спектр квадратичного щебета

Выполните эмпирическое разложение моды и визуализируйте гильбертов спектр сигнала

Гильбертов спектр песни кита

Вычислите гильбертовы параметры спектра сигнала

VMD многокомпонентного сигнала

Вычислите гильбертов спектр сигнала вибрации

Входные параметры

`imf` — Внутренняя функция режима
матрица | расписание

`fs` — Частота дискретизации
2π (значение по умолчанию) | положительная скалярная величина

`freqlocation` — Местоположение оси частоты на графике
`'yaxis'` (значение по умолчанию) | `'xaxis'`

Аргументы name-value

`FrequencyLimits` — Пределы частоты, чтобы вычислить Гильбертов спектр
`[0, fs/2]` (значение по умолчанию) | 1 2 вектор с целочисленным знаком

`FrequencyResolution` — Разрешение частоты, чтобы дискретизировать частотный диапазон
(f_high-f_low)/100 (значение по умолчанию) | положительная скалярная величина

`MinThreshold` — Минимальное пороговое значение Гильбертова спектра
`-inf` (значение по умолчанию) | скаляр

Выходные аргументы

`hs` — Гильбертов спектр сигнала
разреженная матрица

`f` — Значения частоты
вектор

`t` — Временные стоимости
вектор | `duration` массив

`imfinsf` — Мгновенная частота каждого МВФ
вектор | матрица | расписание

`imfinse` — Мгновенная энергия каждого МВФ
вектор | матрица | расписание

Алгоритмы

Ссылки

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

Темы

Документация Wavelet Toolbox

Поддержка

Документация

hht

Синтаксис

Описание

Примеры

Гильбертов спектр квадратичного щебета

Выполните эмпирическое разложение моды и визуализируйте гильбертов спектр сигнала

Гильбертов спектр песни кита

Вычислите гильбертовы параметры спектра сигнала

VMD многокомпонентного сигнала

Вычислите гильбертов спектр сигнала вибрации

Входные параметры

imf — Внутренняя функция режима матрица | расписание

fs — Частота дискретизации2π (значение по умолчанию) | положительная скалярная величина

freqlocation — Местоположение оси частоты на графике 'yaxis' (значение по умолчанию) | 'xaxis'

Аргументы name-value

FrequencyLimits — Пределы частоты, чтобы вычислить Гильбертов спектр [0, fs/2] (значение по умолчанию) | 1 2 вектор с целочисленным знаком

FrequencyResolution — Разрешение частоты, чтобы дискретизировать частотный диапазон (f_high-f_low)/100 (значение по умолчанию) | положительная скалярная величина

MinThreshold — Минимальное пороговое значение Гильбертова спектра -inf (значение по умолчанию) | скаляр

Выходные аргументы

hs — Гильбертов спектр сигнала разреженная матрица

f — Значения частоты вектор

t — Временные стоимости вектор | duration массив

imfinsf — Мгновенная частота каждого МВФ вектор | матрица | расписание

imfinse — Мгновенная энергия каждого МВФ вектор | матрица | расписание

Алгоритмы

Ссылки

Расширенные возможности

Генерация кода C/C++ Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

Темы

Документация Wavelet Toolbox

Поддержка

`imf` — Внутренняя функция режима
матрица | расписание

`fs` — Частота дискретизации
2π (значение по умолчанию) | положительная скалярная величина

`freqlocation` — Местоположение оси частоты на графике
`'yaxis'` (значение по умолчанию) | `'xaxis'`

`FrequencyLimits` — Пределы частоты, чтобы вычислить Гильбертов спектр
`[0, fs/2]` (значение по умолчанию) | 1 2 вектор с целочисленным знаком

`FrequencyResolution` — Разрешение частоты, чтобы дискретизировать частотный диапазон
(f_high-f_low)/100 (значение по умолчанию) | положительная скалярная величина

`MinThreshold` — Минимальное пороговое значение Гильбертова спектра
`-inf` (значение по умолчанию) | скаляр

`hs` — Гильбертов спектр сигнала
разреженная матрица

`f` — Значения частоты
вектор

`t` — Временные стоимости
вектор | `duration` массив

`imfinsf` — Мгновенная частота каждого МВФ
вектор | матрица | расписание

`imfinse` — Мгновенная энергия каждого МВФ
вектор | матрица | расписание

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.