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

wcoh = wcoherence(x,y) возвращает когерентность вейвлета в квадрате величиной, которая является мерой корреляции между сигналами x и y в плоскости частоты времени. Когерентность вейвлета полезна для анализа неустановившихся сигналов. Входные параметры x и y должна быть равная длина, 1D, сигналы с действительным знаком. Когерентность вычисляется с помощью аналитического вейвлета Morlet.

[wcoh,wcs] = wcoherence(x,y) возвращает перекрестный спектр вейвлета x и y. Можно использовать фазу значений перекрестного спектра вейвлета, чтобы идентифицировать относительную задержку между входными сигналами.

[wcoh,wcs,period] = wcoherence(x,y,ts) использует положительный duration ts как интервал выборки. Длительность ts используется для расчета преобразование шкалы к периоду, period. Массив длительности period имеет тот же формат, как задано в ts.

[wcoh,wcs,f] = wcoherence(x,y,fs) использует положительную частоту дискретизации, fs, вычислить преобразование шкалы к частоте, f. Частота дискретизации fs находится в Гц.

[wcoh,wcs,f,coi] = wcoherence(___) возвращает конус влияния, coi, для когерентности вейвлета в циклах на выборку. Если вы задаете частоту дискретизации, fs, конус влияния находится в Гц.

[wcoh,wcs,period,coi] = wcoherence(___,ts) возвращает конус влияния, coi, в циклах в единицу времени.

[___,coi,wtx,wty] = wcoherence(___) возвращает непрерывный вейвлет преобразовывает (CWT) x и y в wtx, wty, соответственно. wtx и wty используются в формировании перекрестного спектра вейвлета и оценок когерентности.

[___] = wcoherence(___,Name,Value) задает дополнительные опции с помощью одного или нескольких аргументов пары "имя-значение". Этот синтаксис может использоваться в любом из предыдущих синтаксисов.

wcoherence(___) без выходных аргументов строит когерентность вейвлета и конус влияния в текущей фигуре. Из-за обратной связи между частотой и периодом, график, который использует интервал выборки, является инверсией графика использование частота дискретизации. Для областей, где когерентность превышает 0.5, графики, которые используют стрелы отображения частоты дискретизации, чтобы показать задержку фазы y относительно x. Стрелы расположены с интервалами вовремя и шкала. Направление стрел соответствует задержке фазы на модульном круге. Например, вертикальная стрелка указывает на π/2 или задержку фазы четверти цикла. Соответствующая задержка вовремя зависит от длительности цикла.

Примеры

Когерентность вейвлета двух Sine wave

Используйте wcoherence по умолчанию настройки, чтобы получать когерентность вейвлета между синусоидой со случайным шумом и модулируемым частотой сигналом с уменьшающейся частотой в зависимости от времени.

t  = linspace(0,1,1024);
x = -sin(8*pi*t) + 0.4*randn(1,1024);
x = x/max(abs(x));
y = wnoise('doppler',10);
wcoh = wcoherence(x,y);

Расчет когерентности по умолчанию использует аналитический вейвлет Morlet, 12 речи на октаву и сглаживает 12 шкал. Количество по умолчанию октав равно floor(log2(numel(x)))-1, который в этом случае равняется 9.

Эффект интервала выборки на когерентности вейвлета

Получите данные о когерентности вейвлета для двух сигналов, задав интервал выборки 0,001 секунд. Оба сигнала состоят из двух синусоид (10 Гц и 50 Гц) в белом шуме. Синусоиды имеют различные поддержки времени.

Установите генератор случайных чисел на его настройки по умолчанию для воспроизводимости. Затем создайте два сигнала.

rng default;
t = 0:0.001:2;
x = cos(2*pi*10*t).*(t>=0.5 & t<1.1)+ ...
cos(2*pi*50*t).*(t>= 0.2 & t< 1.4)+0.25*randn(size(t));
y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+...
sin(2*pi*50*t).*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t));
subplot(2,1,1)
plot(t,x)
title('X')
subplot(2,1,2)
plot(t,y)
title('Y')
xlabel('Time (seconds)')

Получите когерентность двух сигналов.

[wcoh,~,period,coi] = wcoherence(x,y,seconds(0.001));

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

figure
period = seconds(period);
coi = seconds(coi);
h = pcolor(t,log2(period),wcoh);
h.EdgeColor = 'none';
ax = gca;
ytick=round(pow2(ax.YTick),3);
ax.YTickLabel=ytick;
ax.XLabel.String='Time';
ax.YLabel.String='Period';
ax.Title.String = 'Wavelet Coherence';
hcol = colorbar;
hcol.Label.String = 'Magnitude-Squared Coherence';
hold on;
plot(ax,t,log2(coi),'w--','linewidth',2)

Используйте wcoherence(x,y,seconds(0.001)) без любых выходных аргументов. Этот график включает стрелы фазы и конус влияния.

wcoherence(x,y,seconds(0.001));

Воздействие частоты дискретизации на когерентности вейвлета

Получите когерентность вейвлета для двух сигналов, задав частоту дискретизации 1 000 Гц. Оба сигнала состоят из двух синусоид (10 Гц и 50 Гц) в белом шуме. Синусоиды имеют различные поддержки времени.

Установите генератор случайных чисел на его настройки по умолчанию для воспроизводимости и создайте два сигнала.

rng default
t = 0:0.001:2;
x = cos(2*pi*10*t).*(t>=0.5 & t<1.1)+...
    cos(2*pi*50*t).*(t>= 0.2 & t< 1.4)+0.25*randn(size(t));
y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+...
    sin(2*pi*50*t).*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t));

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

wcoherence(x,y,1000)

Получите преобразование шкалы к частоте выход в f.

[wcoh,wcs,f] = wcoherence(x,y,1000);

Эффект количества сглаживавших шкал на когерентности вейвлета

Получите когерентность вейвлета для двух сигналов. Оба сигнала состоят из двух синусоид (10 Гц и 50 Гц) в белом шуме. Используйте количество по умолчанию шкал, чтобы сглаживать. Это значение эквивалентно количеству речи на октаву. Оба значения по умолчанию значений к 12.

Установите генератор случайных чисел на его настройки по умолчанию для воспроизводимости. Затем создайте два сигнала и получите когерентность.

rng default;
t = 0:0.001:2;
x = cos(2*pi*10*t).*(t>=0.5 & t<1.1)+ ...
cos(2*pi*50*t).*(t>= 0.2 & t< 1.4)+0.25*randn(size(t));
y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+...
sin(2*pi*50*t).*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t));
wcoherence(x,y)

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

wcoherence(x,y,'NumScalesToSmooth',18)

Эффект порога отображения фазы на когерентности вейвлета данных о погоде

Сравните эффекты использования различных порогов отображения фазы на когерентности вейвлета.

Постройте когерентность вейвлета между временными рядами Эль-Ниньо и Всем Средним индексом Ливня Индии. Данные производятся ежемесячно. Задайте интервал выборки как 1/12 года, чтобы отобразить периоды в годах. Используйте порог отображения фазы по умолчанию 0,5, который показывает стрелы фазы только там, где когерентность больше или равна 0,5.

load ninoairdata;
wcoherence(nino,air,years(1/12));

Установите порог отображения фазы к 0,7. Количество уменьшений стрел фазы.

wcoherence(nino,air,years(1/12),'PhaseDisplayThreshold',0.7);

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

`x` — Входной сигнал
вектор действительных значений

Входной сигнал в виде вектора действительных значений. x должен быть 1D, сигнал с действительным знаком. Эти два входных сигнала, x и y, должен быть та же длина и должна иметь по крайней мере четыре выборки.

`y` — Входной сигнал
вектор действительных значений

Входной сигнал в виде вектора действительных значений. y должен быть 1D, сигнал с действительным знаком. Эти два входных сигнала, x и y, должен быть та же длина и должна иметь по крайней мере четыре выборки.

`ts` — Интервал выборки
длительность с положительной скалярной величиной вводится

Интервал выборки, также известный как период выборки в виде duration с входом положительной скалярной величины. Допустимой длительностью является yearsдничасысекунды, и minutes. Можно также использовать duration функция, чтобы задать ts. Вы не можете использовать календарную длительность (caldayscalweekscalmonthscalquarters, или calyears).

Вы не можете задать обоих частота дискретизации fs и период выборки ts.

`fs` — Частота дискретизации
положительная скалярная величина | `[]`

Частота дискретизации в виде положительной скалярной величины.

Если вы задаете fs как пустой, wcoherence использование нормировало частоту в циклах/выборке. Частота Найквиста ½.

Вы не можете задать обоих частота дискретизации fs и период выборки ts.

Аргументы в виде пар имя-значение

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

Пример: 'PhaseDisplayThreshold',0.7; задает порог для отображения векторов фазы.

`'FrequencyLimits'` — Пределы частоты
двухэлементный скалярный вектор

Пределы частоты, чтобы использовать в wcoherenceВ виде двухэлементного вектора с положительными строго увеличивающимися элементами. Первый элемент задает самую низкую пиковую частоту полосы пропускания и должен быть больше или быть равен продукту частоты пика вейвлета в герц и в два раза стандартных отклонениях, разделенных на длину сигнала. Второй элемент задает самую высокую пиковую частоту полосы пропускания и должен быть меньше чем или равен частоте Найквиста. Основа 2 логарифма отношения максимальной частоты к минимальной частоте должны быть больше или быть равны 1/NV, где NV является количеством речи на октаву.

Если вы указываете, что частота ограничивает вне допустимой области значений, wcoherence обрезает пределы минимальным и максимальным допустимым значениям. Используйте cwtfreqbounds с набором вейвлета к 'amor' определить частоту ограничивает для различной параметризации когерентности вейвлета.

Пример: 'FrequencyLimits',[0.1 0.3]

`'PeriodLimits'` — Пределы периода
двухэлементный массив длительности

Пределы периода, чтобы использовать в wcoherenceВ виде двухэлементного массива длительности со строгим увеличением положительных элементов. Первый элемент должен быть больше или быть равен 2×ts где ts период выборки. Основа 2 логарифма отношения минимального периода к максимальному периоду должны быть меньше чем или равны-1/NV, где NV является количеством речи на октаву. Максимальный период не может превысить длину сигнала, разделенную на продукт двух раз стандартные отклонения вейвлета и частоты пика вейвлета.

Если вы указываете, что период ограничивает вне допустимой области значений, wcoherence обрезает пределы минимальным и максимальным допустимым значениям. Используйте cwtfreqbounds с набором вейвлета к 'amor' определить период ограничивает для различной параметризации когерентности вейвлета.

Пример: 'PeriodLimits',[seconds(0.2) seconds(1)]

Типы данных: duration

`'VoicesPerOctave'` — Количество речи на октаву
12 (значений по умолчанию) | даже целое число от 10 до 32

Количество речи на октаву, чтобы использовать в когерентности вейвлета в виде ровного целого числа от 10 до 32.

`'NumScalesToSmooth'` — Количество шкал, чтобы сглаживать
положительное целое число

Количество шкал, чтобы сглаживать вовремя и масштабироваться в виде положительного целого числа, меньше чем или равного одной половине N, где N является количеством шкал в вейвлете, преобразовывает. Если незаданный, NumScalesToSmooth значения по умолчанию к минимуму floor(N/2) и VoicesPerOctave. Функция использует фильтр скользящего среднего значения, чтобы сглаживать через шкалу. Если ваша когерентность является шумной, можно задать больший NumScalesToSmooth значение, чтобы сглаживать когерентность больше.

`'NumOctaves'` — Количество октав
положительное целое число

Количество октав, чтобы использовать в когерентности вейвлета в виде положительного целого числа между 1 и floor(log2(numel(x)))-1. Если вы не должны исследовать более низкие значения частоты, используйте меньший NumOctaves значение.

'NumOctaves' пара "имя-значение" не рекомендуется и будет удалена в будущем релизе. Рекомендуемый способ изменить частоту или область значений периода когерентности вейвлета с 'FrequencyLimits' или 'PeriodLimits' пары "имя-значение". Вы не можете задать обоих 'NumOctaves' и 'FrequencyLimits' или 'PeriodLimits' пары "имя-значение". Смотрите cwtfreqbounds.

`'PhaseDisplayThreshold'` — Порог для отображения векторов фазы
0.5 (значение по умолчанию) | действительный скаляр между 0 и 1

Порог для отображения векторов фазы в виде действительного скаляра между 0 и 1. Эта функция отображает векторы фазы для областей с когерентностью, больше, чем или равный заданному пороговому значению. Понижение порогового значения отображает больше векторов фазы. Если вы используете wcoherence с любыми выходными аргументами, PhaseDisplayThreshold значение проигнорировано.

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

`wcoh` — Когерентность вейвлета
матрица

Когерентность вейвлета, возвращенная как матрица. Когерентность вычисляется с помощью аналитического вейвлета Morlet по логарифмическим масштабам со значением по умолчанию 12 речи на октаву. Количество по умолчанию октав равно floor(log2(numel(x)))-1. Если вы не задаете интервал выборки, частота дискретизации принята.

`wcs` — Перекрестный спектр вейвлета
матрица комплексных чисел

Перекрестный спектр вейвлета, возвращенный как матрица комплексных чисел. Можно использовать фазу значений перекрестного спектра вейвлета, чтобы идентифицировать относительную задержку между входными сигналами.

`period` — Преобразование шкалы к периоду
массив длительности

Преобразование шкалы к периоду, возвращенное как массив длительности. Значения преобразования вычисляются из периода выборки, заданного в ts. Каждый period элемент имеет тот же формат как ts.

`f` — Преобразование шкалы к частоте
вектор

Преобразование шкалы к частоте, возвращенное как вектор. Вектор содержит пиковые значения частоты для вейвлетов, использовался для расчета когерентности. Если вы хотите к выходу f, но не указывайте, что частота дискретизации ввела, fs, возвращенная когерентность вейвлета находится в циклах на выборку.

`coi` — Конус влияния
массив типа double | массив длительности

Конус влияния за когерентность вейвлета, возвращенную или как массив типа double или как массив длительности. Конус влияния указывает, где краевые эффекты происходят в данных о когерентности. Если вы задаете частоту дискретизации, fs, конус влияния находится в Гц. Если вы задаете интервал выборки или период, ts, конус влияния находится в периоды. Из-за краевых эффектов, верьте меньше к областям очевидной высокой когерентности, которые находятся вне или перекрывают конус влияния. Конус влияния обозначается пунктирной линией.

Для получения дополнительной информации смотрите Граничные эффекты и Конус Влияния.

`wtx` — Непрерывное преобразование вейвлета `x`
матрица

Непрерывное преобразование вейвлета x, возвращенный как матрица.

`wty` — Непрерывное преобразование вейвлета `y`
матрица

Непрерывное преобразование вейвлета y, возвращенный как матрица.

Больше о