Преобразование Гильберта в среде LabVIEW

№ 2’2012
На протяжении многих десятилетий в рамках радиоэлектроники разрабатываются и совершенствуются методы обработки сигналов, а также методы их анализа и выделения полезной информации. Среди них особое место занимает преобразование Гильберта.

Преобразование Гильберта

Преобразованием Гильберта сигнала s(t) называется интегральное преобразование вида:

где несобственный интеграл понимается в смысле главного значения по Коши, то есть:

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

Аналитический сигнал

Преобразование Гильберта позволяет сформировать из анализируемого сигнала — действительной функции времени s(t) — аналитический сигнал:

В результате в распоряжении исследователя оказываются три функции времени:

  • мгновенная амплитуда (огибающая сигнала):
  • мгновенная частота:
  • мгновенная фаза:

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

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

Все это наглядно свидетельствует о полезности применения преобразования Гильберта для детального анализа сигналов.

Алгоритмическая реализация

Эффективная алгоритмическая реализация преобразования Гильберта на практике [1] опирается на соотношения:

Поэтому, если сигнал представлен в форме разложения в ряд Фурье, то у сигнала, преобразованного по Гильберту, меняются местами действительная и мнимая части спектра, и спектр сигнала, преобразованного по Гильберту S(t), можно записать так:

С учетом (6) алгоритм вычисления преобразования Гильберта, используя дискретное (быстрое) преобразование Фурье (БПФ), можно представить в виде, показанном на рис. 1.

Эффективный алгоритм преобразования Гильберта

Рис. 1. Эффективный алгоритм преобразования Гильберта

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

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

В статье описана процедура использования пакета прикладных программ LabVIEW [3] и приведен ряд примеров его применения для решения задач, имеющих практическую значимость.

В пакете прикладных программ LabVIEW реализован алгоритм построения преобразования Гильберта входного сигнала s(f). Он отображается пиктограммой

и включает в себя последовательное выполнение следующих операций:

  • Преобразование Фурье входной последовательности X: Y = F(X).
  • Обнуление постоянной составляющей: YQ = 0.0.
  • Выполнение логической операции: если длина последовательности четна, то компонента на частоте Найквиста приравнивается к нулю.
  • Выполнение операции формирования последовательности Hk: Hk = –jsign(k)Yk.

Полученное преобразование Гильберта ŝ(t) сигнала s(t) используется для формирования огибающей сигнала ω(t), его мгновенной частоты f(t) и мгновенной фазы φ(t). Так, для определения огибающей и мгновенной фазы входного сигнала найденный выше массив данных необходимо преобразовать согласно (3) и (5). Утилита, реализующая алгоритм указанного преобразования, изображена на рис. 2.

Структурная схема формирователя огибающей ω(t) и мгновенной фазы φ(t) входного сигнала (t)

Рис. 2. Структурная схема формирователя огибающей ω(t) и мгновенной фазы φ(t) входного сигнала (t)

Для определения мгновенной частоты f(t) выполняется последовательность действий согласно (4).

Структурная схема формирователя огибающей ω(t), мгновенной фазы φ(t) и мгновенной частоты f(t) входного сигнала s(t) представлена на рис. 3.

Схема формирователя огибающей ω(t), мгновенной фазы φ(t) и мгновенной частоты f(t) входного сигнала s(t)

Рис. 3. Структурная схема формирователя огибающей ω(t), мгновенной фазы φ(t) и мгновенной частоты f(t) входного сигнала s(t)

Примеры

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

  1. ЛЧМ-сигнал формируется в соответствии с формулой:
  2. где A — амплитуда; a = 2π(f2–f1)/n; b = 2πf1; f1 — начальная частота; f2 — конечная частота; n — число отсчетов.

    Огибающая мгновенная фаза и мгновенная частота аналитического ЛЧМ-сигнала представлены на рис. 4.

    Огибающая аналитического ЛЧМ-сигнала: а) мгновенная фаза; б) мгновенная частота

    Рис. 4. Огибающая аналитического ЛЧМ-сигнала: а) мгновенная фаза; б) мгновенная частота

  3. Затухающий сигнал формируется в соответствии с формулой:
  4. где α — величина затухания; Ω — несущая частот.

    Огибающая затухающего аналитического сигнала представлена на рис. 5.

    Сигнал и огибающая затухающего аналитического сигнала

    Рис. 5. Сигнал и огибающая затухающего аналитического сигнала

  5. АМ-сигнал формируется в соответствии с формулой:
  6. где ω и Ω — несущая и модулирующая частоты; m — индекс глубины модуляции.

    Огибающая мгновенная фаза и мгновенная частота аналитического АМ-сигнала представлены на рис. 6.

    АМ-сигнал и его огибающая

    Рис. 6. АМ-сигнал и его огибающая

  7. Фронт радиоимпульса формируется в соответствии с формулой:
  8. где ω — несущая частота; Heaviside(x) — функция Хэвисайда; Т — величина времени задержки относительно t = 0.

Фронт радиоимпульса и его огибающая представлены на рис. 7.

Фронт радиоимпульса и его огибающая

Рис. 7. Фронт радиоимпульса и его огибающая

Таким образом, с помощью пакета LabVIEW можно формировать и анализировать различные типы аналитических сигналов.

Литература

  1. Вопросы подповерхностной радиолокации. Коллективная монография / Под ред. А. Ю. Гринева. М.: Радиотехника, 2005.
  2. Френкc Л. Теория сигналов / Под ред. В. Е. Вакмана. М.: Советское радио, 1974.
  3. Суранов А. Я. LabVIEW 8.20. Справочник по функциям. М.: ДМК Пресс, 2007.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *