CEEMDAN (англ. Complete Ensemble Empirical Mode Decomposition with Adaptive Noise) -- метод разложения на эмпирические функции, основанный на множественном добавлении белого шума к различным копиям исходного сигнала. Сами моды же вычисляются, как среднее значение от всех реализаций.
Определим оператор , который возвращает i-ю моду входного сигнала S, подсчитанную при помощи EMD. Также обозначим белый шум с нулевым средним и дисперсией, равной единице. Тогда если x[n] - входные данные, то разложение можно описать следующим алгоритмом:
1. Используя EMD, найдем первую моду I реализаций и подсчитаем
;
2. На первой итерации (k=1) найдем первый остаток ;
3. Раскладываем до первой моды реализации используя EMD, и затем определяем вторую моду:
;
4. Для k=2,...,K считаем k-е остатки:
;
5. Раскладываем до первой моды реализации используя EMD, и затем определяем (k+1) моду:
,
Шаги с 4 по 6 повторяются до тех пор, пока из получившегося остатка уже невозможно выделить новую моду (у остатка нет как минимум двух экстремумов). Конечный остаток удовлетворяет соотношению:
,
где K -- итоговое количество мод. В итоге входной сигнал x[n] удовлетворяет разложению:
.
Данный подход, основанный на многочисленном добавлении шума, а, как следствие, и локальных экстремумов, позволяет точнее разложить зашумленные процессы. Также такой подход эффективен при работе с импульсными источниками. Благодаря CEEMDAN, удается обойти такие проблемы EMD, как перемешивание мод: моды с близкими частотами теперь действительно две разные моды, и не может быть такого, что в одной моде присутствуют различные частоты. Кроме того, удается избежать размазывания энергии по модам.
Фильтрация шумов
Рис. 14: Зависимость энергии модовой функции белого шума от ее периода. По обоим осям логарифмические координаты.
Одним из основных преимуществ эмпирической декомпозиции является то, что базисные функции подбираются в процессе разложения, то есть для каждого конкретного случая они уникальные. Однако, несмотря на простоту нахождения модовых функций, у них есть ряд преимуществ. Во-первых, каждая мода несет определенную частоту, что сильно упрощает спектральный анализ сигнала. Во-вторых, выделяются моды различной амплитуды, и это позволяет фильтровать шумовые моды, то есть такие, которые несут малую долю энергии. В-третьих, анализируя белый шум, можно найти зависимость энергии моды шума от ее периода. Это будет прямая (см Рис. 14). Если при анализе сигнала получились моды, которые ложатся на данную прямую, можно прийти к выводу, что данные моды -- шумовые, и рассматривать сигнал уже с вычетом таких мод.
2.2 Экспериментальные установки
Экспериментальная база на горе Арагац
Первая часть моей работы состоит в реализации метода Эмпирической Модовой Декомпозиции и разработке системы определения локации микроразрядов, как это было описано в [14]. Необходимы были данные радиоизлучения в ОВЧ и УВЧ диапазонах с трех антенн, расположенных на перпендикулярных базовых линиях. Данные предоставляла научная база, расположенная в Армении на горе Арагац. Поскольку научная станция расположена на горе, влияние атмосферы на получаемые данные меньше, а грозовые облака находятся ближе к станции, чем при измерениях с земли.
В работе использовались данные с трех плоских антенн, отфильтрованные в полосе частот от 24 до 82 МГц. Также на установке присутствует четвертая триггерная антенна, при срабатывании которой велась запись последующих 400 мс, к которым добавлялись предшествующие срабатыванию триггера 100 мс. Таким образом, длительность записи в один файл составляла 0.5 секунды.
Отфильтрованные данные с антенн попадали в 4-х канальный осциллограф с частотой дискретизации 5000 сэмплов в секунду, то есть 1250 сэмплов в секунду на канал, что соответствует длительность сэмпла в 6.4 наносекунды. На Рис. 15 и Рис. 16 показана конфигурация антенн и фотография здания, где они находились.
Рис. 15: Три плоские антенны, расположенные на научной станции на горе Арагац. Антенны расположены на двух перпендикулярных базовых линиях.
Рис. 16: Три плоские антенны УВЧ интерферометра на научной базе на горе Арагац
На Рис. 17 показан пример отображения данных с интерферометра. Для хранения и отображения используется приложение PicoSoft.
Рис. 17: Данные, получаемые интерферометром. Верхние три временных ряда - запись УВЧ радиоизлучения с трех плоских антенн. Нижний временной ряд - данные с триггерной антенны. Данные хранятся в формате приложения PicoSoft. Изображенные данные получены 25 сентября 2017 года, в 09:57:02.834
Экспериментальная база в Тарусе
В Тарусе также, как и на Арагаце, находится четырехканальный интерферометр, антенны которого формируют прямоугольный треугольник. На Рис. 18 представлена схема расположения антенн в г. Таруса. Данные того же формата, что и на Арагаце.
2.3 Используемое программное обеспечение
Как уже упоминалось выше, данные конвертируются в формат приложения PicoSoft, в котором хранятся и отображаются в этом же приложении. Моя разработка велась на языке Python 3.6.3. Дополнительно использовались библиотеки numpy, scipy, pandas для вычислений и matplotlib для отображения результатов программирования.
3. Эксперимент
3.1 Реализация методов EMD и CEEMDAN
Рис. 18: Схема расположения антенн интерферометра на научной базе в Тарусе
Здесь представлены результаты работы реализованного метода модовой декомпозиции.
Рис. 19: Сигнал, состоящий из суммы трех синусов и шума
Рис. 20: Результат разложения сигнала с Рис. 18 методом EMD
Теперь рассмотрим, как справляются с разложением импульсных сигналов оба метода. В качестве входных сигналов использовался дельта импульс и зашумленный Гаусс.
Результаты представленны на рисунках Рис. 22 - Рис. 27.
Рис. 21: Результат разложения сигнала с Рис. 18 методом CEEMDAN
Рис. 22: Дельта импульс
Рис. 23: Разложение дельта-импульса методом EMD
Рис. 24: Разложение дельта импульса методом CEEMDAN
Рис. 25: биполярный гауссовый импульс и белый шум
Рис. 26: Разложение биполярного гаусса методом EMD
3.2 Построение карты эволюции молнии
Рис. 27: Разложение биполярного гаусса методом CEEMDAN
Здесь представлены результаты расчета временной задержки, углов азимута и элевации, а также полученная картина молнии.
Рис. 28: Два временных ряда, полученные плоскими антеннами интерферометра на научной станции на горе Арагац
Москва
17 июня 2019
Рис. 29: Графики функций кросс-корреляции: построенный методом перевода данных в частотное пространство (сверху) и обычный, реализованный во временном пространстве (снизу).
Слева от каждого графика представлена увеличенная область вблизи максимума функций. В этой области также отображены данные о значении временной задержки между сигналами (Рис. 27) и время вычисления временной задержки.
Из результатов видно, что время вычисления временной задержки при помощи метода перевода в частотный домен и обратно (GCC-PHAT) больше, чем на порядок ниже, чем при построении обычной функции кросс-корреляции. С учетом того, что в одном файле очень много подобных событий, выбор метода вычисления временной задержки играет ключевую роль в производительности вычислений.
3.3 Корреляция между количеством событий и их масштабом
Благодаря высокому временному разрешению интерферометров и хорошему соотношению сигнал-шум, в полученных данных удалось найти масштабную зависимость.
Ниже представлены исследуемые временные ряды с гистограммами. На всех анализируемых событиях видна одна и та же зависимость количества событий в зависимости от величины поля, которая, как известно, зависит от пространственного масштаба разряда.
Рис. 30: Местоположения внутриоблачных микроразрядов на небесной сфере.
Рис. 31: Здесь и далее:
Важно отметить, что, исходя из представленных выше данных, в грозовом облаке присутствует два механизма формирования зарядов, причем каждый удовлетворяет масштабной зависимости: доминируют мелкомасштабные события, в то время как крупномасштабные составляют малую часть.
Также интересно посмотреть за эволюцией наклона прямой: будет ли она становится более пологой со временем? Это, соответственно, означало бы, что со временем все мелкомасштабные события перерастают в крупномасштабные, что соответствует обсуждаемому выше механизму слияния стримеров (раздел 1.5.3 Компьютерное моделирование). Для этого необходимо проводить измерения с боМльшим временным окном.
4. Результаты
Как известно, физику, как науку, можно поделить на экспериментальную и теоретическую. Это две составляющие, которые тесно связаны между собой. Любая идея должна подтверждаться экспериментами, чтобы стать самодостаточной теорией, поэтому любой физик рано или поздно встречается с экспериментом. А где эксперимент, там измерения, и где измерения, там помехи, шумы и анализ данных.
Поскольку физика молниевых разрядов - одна из самых больших загадок физики атмосферы, в этой области существует множество теорий о том или ином процессе, и все они требуют экспериментального подтверждения. А так как экспериментальные измерения внутренних параметров грозовых облаков очень сильно затруднены (1.3 Трудности проведения измерений), интерпретация и анализ любых данных о внутриоблачных разрядах крайне важны и сложны. Поэтому меня, как физика и программиста, и заинтересовала эта актуальная и нетривиальная задача анализа сильно зашумленных нестационарных данных.
В работе рассмотрены методы разложения нестационарного и нелинейного сигнала, альтернативные спектральному анализу Фурье, который дает нефизические результаты для подобных данных. В качестве таких методов выступали алгоритмы Эмпирической Модовой Декомпозиции (Empirical Mode Decomposition) и его улучшения CEEMDAN (Complete Ensemble Empirical Mode Decomposition with Adaptive Noise). В этих методах базис разложения задается на ходу и уникален для каждого сигнала, что хорошо подходит при рассмотрении молниевых разрядов (см. Рис. 27: Разложение биполярного гаусса методом CEEMDAN).
В ходе тестирования реализованного метода CEEMDAN было обнаружено, что при разложении белого шума по собственным модам их энергия и частота связаны линейно в двойном логарифмическом масштабе (см. Рис. 14), что позволяет фильтровать шумовые моды при разложении сигнала.
По аналогии с [14], был разработан алгоритм определения положения разрядов на небесной сфере, который после тестировался на данных, полученных от научных баз в Тарусе и Армении (Рис. 30).
Также, в результате более подробного анализа полученных данных, была найдена закономерность в количестве разрядов определенной силы: число разрядов уменьшается с увеличением их пространственного масштаба. Поскольку большой интерес представляет поведение наклона подобной гистограммы с течением времени, имеется потребность в измерениях с увеличенным временным окном. Ожидается, что со временем мелкомасштабных разрядов будет становиться меньше, в то время как крупномасштабные разряды, сформированные механизмом слияния стримеров, будут встречаться чаще.
Несмотря на то, что эта работа посвящена обработке данных интерферометра средствами программирования, физические нити непрерывно пронизывают её страницы, преображая пустой поток нулей и единиц в физически осознанные данные, воссоздавая удивительные и загадочные природные явления. В попытках разгадать тайны молний ученые находят все новые вопросы, ответы на которые требуют все новых экспериментов. В результате новых измерений будут получаться новые данные.
И именно в них таятся все разгадки физики молний.
Список литературы
1. Д. И. Иудин, С. С. Давыденко, В. М. Голиб, М. С. Долгоносов и Л. М. Зеленый, «Физика молний: новые подходы к моделированию и перспективы спутниковых наблюдений,» Успехи физических наук, т. 188, № 8, pp. 850-864, Август 2018.
2. J. R. Dwyer и M. A. Uman, «The physics of lightning,» Physics Reports, т. 534, № 4, pp. 147-241, 2014.
3. V. A. Rakov и M. A. Uman, Lightning: Physics and Effects, ред., т. , , : Cambridge University Press, , p. .
4. V. A. Rakov, «New insights into lightning processes gained from triggered-lightning experiments in Florida and Alabama,» Geophysical Research Letters, 1998.
5. V. A. Rakov, M. A. A. Uman и K. J. J. Rambo, «A review of ten years of triggered-lightning experiments at Camp Blanding, Florida,» Atmospheric Research, pp. 503-517, 28 11 2004.
6. J. R. e. a. Dwyer, «A ground level gamma-ray burst observed in association with rocket-triggered lightning,» Geophysical Research Letters, 2004.
7. V. A. Rakov и M. A. Uman, Lightning: physics and effects, ред., т. , , : Cambridge University Press, 2007, p. .
8. R. H. Golde, Lightning, Vol. 1, Physics of Lightning, Academic Press, N.Y., 1977, p. 368.
I. B. Cohen, Benjamin Franklin's Experiments: a new edition of Franklin's Experiments and observations on electricity, ред., т. , , : Harvard University Press, 1941, p. .
9. C. V. Boys, «Progressive Lightning: A New Stereoscope,» Nature, т. 131, № 3310, pp. 492-494, 1933.
10. S. . Flьgge и S. A. Korff, «Gas Discharges II. Vol. 22 of Handbuch der Physik,» Physics Today, т. 10, № 1, pp. 45-46, 1957.
11. F. . Pockels, «Ueber das magnetische Verhalten einiger basaltischer Gesteine,» Annalen der Physik, т. 299, № 13, pp. 195-201, 1897.
12. K. . Berger, «Novel observations on lightning discharges: Results of research on Mount San Salvatore,» Journal of The Franklin Institute-engineering and Applied Mathematics, т. 283, № 6, pp. 478-525, 1967.