Этап 1. Выделение особенностей и выбор длины участков шума. Формально для нахождения длины участка разбиения требуется построить детектор, который позволяет выбрать одну из двух гипотез: H0 — выбранный участок соответствует фону (наиболее часто повторяющимся участкам с одинаковыми признаками); H1 – существует момент времени τ, такой, что исследуемый участок шума соответствует фону при t<τ и нехарактерным (маргинальным) особенностям, вызванным звуками отдельных пчел при t≥τ.
В анализируемом шуме выбираем пробный интервал [tk,tk+τ], при этом tk+τ < T, где tk = kΔτ, k=0, 1, 2, 3, … . Временные интервалы τ и Δτ определяются характерными временами, наиболее важными для динамики изменения исследуемого шума. Значение Δτ соответствует интервалу между отчетами оцифрованного сигнала, определяемого по теореме Котельникова, и зависит от длины сигнала T и частоты дискретизации. Для нестационарного шума имеется набор характерных времен τ, определяемых локальной структурой исследуемого шумового процесса.
Для выбранного участка шума длиной τ определяем случайную последовательность nj (j=1,... N), которая представляет собой разность значений исходных данных и их среднего арифметического. Таким образом, случайная последовательность nj представляет собой временной ряд в равноотстоящих точках хj. Далее проведем ее рекуррентное интегрирование по формуле
Jj = Jj–1 + (хj – хj–1)Rj, (1)
где
Rj = 0,5(nj + nj–1). (2)
Полученную последовательность Jj сглаживаем с помощью оптимального метода линейного сглаживания (ОМЛС) с гауссовым ядром по минимальной величине относительной ошибки (В.А. Тобоев, 2008). Такая процедура позволяет отделить относительные флуктуации от тренда*, содержащего только низкочастотные составляющие.
* Здесь и ниже под трендами мы понимаем сглаженные участки временной акустической последовательности, которые определяют, по сути, усредненную функцию средних значений, отражающую поведение этой функции на больших временах.
Длину участка τопт выбираем по максимальному значению обобщенной функции корреляции Пирсона (ОФКП), вычисленному между выделенными трендами для участков одинаковой длины (В.А. Тобоев, 2008).
Этап 2. Кластеризация выделенных участков по статистически однородному признаку (признаковое описание). На данном этапе анализируем свойства сигнала на каждом выбранном участке. Пусть Sk – участок шума длиной τопт, соответствующий k-му фрагменту нестационарного акустического шума. Единый алгоритм обработки данных каждого отдельного участка включает в себя следующие основные пункты. Для выбранного участка шума вновь находим интегрированную последовательность и ее сглаженное среднее значение. Полученный тренд вычитаем из проинтегрированной последовательности, что позволяет определить искомую случайную последовательность распределения относительных флуктуаций Yk. У нее вновь вычисляем тренд с помощью ОЛМС, который дает горизонтальную линию, практически совпадающую с осью абсцисс. Последовательность Yk упорядочиваем по амплитудам. Когда они подчиняются условию y1k > y2k >...>yNk, упорядоченный набор величин {yik} образует последовательность ранжированных амплитуд (ПРА).
Для количественного сравнения относительных флуктуаций различных участков сигнала, огибающие ПРА заменяем подгоночными функциями F(t) вида
F(t)=Citαj, (3)
что позволяет разделить амплитуды нестационарного сигнала yjk (j=1,... N) на некие оптимальные статистические группы (кластеры) с параметрами (Cm, αm), что соответствует редуцированному описанию рассматриваемого сигнала.
В.А.ТОБОЕВ
Чувашский государственный университет

Синтетические адаптогены и биостимулятор…
апр 1, 2015
Горчица — дар людям и пчелам…
мая 6, 2016
Геометрия ячеек пчел
окт 7, 2017
Стандартные ульи выгоднее …
март 31, 2019
Температурный стресс и пчелиный яд…
июль 11, 2015
Cs-137 в продукции пчел и медоносах на ф…
мая 5, 2015
Шкафы-секции для хранения сотов…
янв 25, 2020
Интенсивная технология производства мато…
нояб 13, 2015
Роевня — какой она должна быть?…
сен 5, 2022
Формирование отводков-резервов…
март 6, 2017
Противороевой метод c тихой сменой маток…
июнь 20, 2017
Как шмели определяют качество пыльцы?…
фев 2, 2017
Масса капель меда и скорость его течения…
июнь 5, 2022
Технология использования пчел на главном…
авг 13, 2021
Пчеловодство без роев или «чересполосица…
июль 17, 2015