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

Лекция № 2

Метрологические
основы методов экологических исследований

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

В
случае необходимости обнаружения
какого-либо компонента обычно фиксируют
появление
аналитического сигнала — появление
осадка, окраски, линии в спектре и т. д.
Появление аналитического сигнала должно
быть надежно зафиксировано. При
определении количества компонента
измеряется величина
аналитического
сигнала: масса осадка, сила тока,
интенсивность линии спектра и т. д. Затем
рассчитывают содержание компонента с
использованием функциональной зависимости
аналитический сигнал — содержание:
y=f(c),
которая устанавливается расчетным или
опытным путем и может быть представлена
в виде формулы, таблицы или графика.
Содержание при этом может быть выражено
абсолютным количеством определяемого
компонента в молях, в единицах массы
или через соответствующие концентрации.

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

фона,
обусловленного примесями определяемого
компонента и мешающими компонентами в
растворах, растворителях и матрице
образца, а также «шумами», возникающими
в измерительных приборах, усилителях
и другой аппаратуре. Эти шумы не имеют
отношения к определяемому компоненту,
но накладываются на его собственный
аналитический сигнал. Задача аналитика
состоит в том, чтобы максимально снизить
величину аналитического сигнала фона
и, главное, сделать минимальными его
колебания. Обычно аналитический сигнал
фона учитывают при проведении контрольного
(холостого) опыта
,
когда через все стадии химического
анализа проводится проба, не содержащая
определяемого компонента. Полезным
сигналом при этом будет аналитический
сигнал, равный разности измеренного
аналитического сигнала и аналитического
сигнала фона

.

Метод градуировочного
графика ,

метод добавок…Метод
стандартов.

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

Наиболее
распространен метод
градуировочного графика
.
При этом в координатах аналитический
сигнал — содержание компонента строят
график с использованием образцов
сравнения с различным и точно известным
со­держанием определяемого компонента.
Затем, измерив величину аналитиче­ского
сигнала в анализируемой пробе, находят
содержание определяемого компонента
по градуировочному графику (рис. 2.2).

В
методе
стандартов

измеряют аналитический сигнал в образце
сравнения (эталонном образце) с известным
содержанием компонента и в анализируемой
пробе: уэт
= Scэт
и ух
= Scх,
где S
— коэффициент пропорциональности. Если
определенное в идентичных условиях
значение S заранее известно, то ‘можно
провести расчет по формуле cх
= yх/S.
Обычно же применяют соотношение уэтх
= сэтх,
откуда

Сх=YxCэт/Yэт

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

Cxэт,1+(Сэт,2
– С
эт,1)(Yx
Yэт,1)
/
Yэт,2
Yэт,1

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


Рис.
2.2. Метод градуировочного
Рис. 2.3. Метод добавок
графика

При
определении содержания расчетным
методом берут две аликвоты раствора
анализируемой пробы. В одну из них вводят
добавку определяемого компонента
известного содержания. В обеих пробах
измеряют аналитический сигнал — Ух и
Ух+цоб- Неизвестную концентрацию
определяемого компонента рассчитывают
по формуле

Cx=YxVдобСдоб
/
Yx+добVдоб
+ (
Yx+доб
Yx)V

При
определении содержания компонента
графическим методом берут аликвоты
анализируемой пробы: 1, 2, 3, …, п. В аликвоты
2, 3, …, п вводят известные, возрастающие
количества определяемого компонента.
Во всех аликвотах измеряют аналитический
сигнал и строят график в координатах
аналитический сигнал — содержание
определяемого компонента, приняв за
условный нуль содержание определяемого
компонента в аликвоте без добавки
(аликвота 1). Экстраполяция полученной
прямой до пересечения с осью абсцисс
дает отрезок, расположенный влево от
условного нуля координат, величина
которого в выбранном масштабе и единицах
измерения соответст­вует искомому,
содержанию (с^) определяемого компонента
(рис. 2.3).

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

Для
построения градуировочного графика,
наилучшим образом удовлетворяющего
экспериментальным данным, обычно
используют метод наименьших квадратов
(МНК). Основное положение МНК утверждает,
что если для каждой из т экспериментальных
точек провести на оптимальную кривую
прямые, параллельные оси ординат, то
для Оптимальной (теоретической) кривой
сумма квадратов отклонений точек от
кривой ( д,) должна быть минимальной, т.
е.

В
химическом анализе чаще всего используют
прямолинейные градуировочные графики,
построенные для определенного диапазона
определяемых содержаний, т. е. в области
значений, предусмотренных данной
методикой. Уравнение прямой можно
записать в виде у = а + Ъх. Если имеется
т экспериментальных точек (xi>>),
x-iYi,
…, х^Уу,), то, используя постулат МНК
,можем найти параметры а и Ъ прямой,
наилучшим образом удовлетворяющей
экспериментальным данным

Заметим,
что параметры а и Ъ тем точнее и,
следовательно, построенный градуировочный
график тем ближе к теоретическому, чем
в более широком диапазоне определяемых
содержаний строится график и чем большее
число образцов сравнения (эталонов) т
было взято для его построения.

Во
всех методах определения неизвестного
содержания компонента используют
функциональную зависимость у = Sx.
Коэффициент чувствительности S
(иногда его называют просто чувствительность)
характеризует отклик аналитического
сигнала на содержание компонента.
Коэффициент чувствительности — это
значение первой производной градуировочной
функции при данном определенном
содержании. Для прямолинейный
градуировочных графиков – это тангенс
угла наклона прямой (см. рис.2,2)

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

Во
всех рассмотренных способах используют
образцы сравнения (эталоны), т. е. образцы,
пробы, растворы с точно установленным
содержанием компонента. Методы анализа,
использующие образцы сравнения — это
так называемые относительные методы
химического анализа. Абсолютных
(безэталонных) методов в аналитической
химии немного — например, методы
гравиметрии, прямой кулонометрии,
некоторые варианты радиохимических
методов.

Образцы
сравнения для относительных методов
анализа могут быть приготовлены из
химически чистых, устойчивых веществ
известного состава (стандартные
вещества). В этом случае содержание
определяемого компонента вычисляют по
химической формуле стандартного
вещества. Возможно приготовление
образцов сравнения в отдельной
лаборатории, учреждении, отрасли, когда
содержание компонента устанавливают
разными методами, на разных приборах
многие аналитики. Наиболее надежные
результаты получают, когда в качестве
образцов сравнения используют стандартные
образцы (СО) — специально приготовленные
материалы, состав и свойства которых
достоверно установлены и официально
аттестованы специальными го­сударственными
метрологическими учреждениями.

В
химических методах анализа

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

Аналитическая реакция представляет
собой химическую реакцию аналитически
активной формы определяемого компонента
с другим веществом – аналитическим
реагентом, в результате чего происходят
видимые или доступные для регистрации
изменения в анализируемой системе. Для
химических методов анализа характерно
использование специальных реагентов,
способных к практически полному
необратимому
взаимодействию с определяемым компонентом.
Само измерение
носит косвенный характер
,
поскольку непосредственно определяется
не число определенных частиц, а их
совокупная масса
(в гравиметрии
),
расход
реагента (в титриметрии
)
или другие параметры, пропорциональные
массе
.

Чувствительность
реакции

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

Предельная
концентрация


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

Предельная
концентрация может быть выражена в
единицах молярной концентрации
определяемого вещества или как отношение
1 массы определяемого вещества к массе
растворителя. Например, для реакции
определения ионов железа по реакции с
роданид-ионом:

Предельная
концентрация равна 610-6
моль/л или 1: 3106.
Это значит, что можно обнаружить ионы
железа по реакции образования окрашенного
роданида железа при концентрации железа
610-6
моль/л или при содержании 1 г ионов
железа

в 3106
г растворителя (наблюдается едва заметная
розовая окраска) . При меньшей концентрации
железа окрашивания раствора при
добавлении роданида не обнаруживается,
при большей концентрации железа
добавление роданида вызывает появление
хорошо заметной красной окраски раствора.

Открываемый
минимум


минимальное количество вещества, которое
может быть обнаружено в данных условиях.
Величина открываемого минимума зависит
от предельной концентрации и объема
раствора. По существу, открываемый
минимум – это масса вещества, которая
содержится в данном объеме раствора
при концентрации раствора равной
предельной концентрации.

Например,
для реакции ионов железа с роданид-ионом
предельная концентрация равна 610-6
моль/л. Тогда открываемый минимум m
может быть определен, используя
определение молярной концентрации:


;


;

, г

Так,
если для анализа взяты разные объемы
исследуемых растворов, то значения
открываемого минимума также будут
разными:

  • Объем исследуемого
    раствора 1 мл

  • Объем исследуемого
    раствора 10 мл


.

Чувствительность
реакции зависит от многих факторов
:
объема раствора, концентрации, порядка
добавления реактивов, рН и температуры
раствора, условий наблюдения за
результатами реакции и др. Данные
о чувствительности какой-либо реакции,
приведенные без указания условий ее
выполнения, являются величинами весьма
неопределенными и невоспроизводимыми
.
Чем чувствительнее реакция, тем больше
возникает трудностей и препятствий при
анализе малых концентраций веществ,
тем больше следует обращать внимания
на самые разные факторы. Так, при
определении чувствительности реакции
или метода приходится считаться с тем
обстоятельством, что в контрольном
опыте с одними только реактивами и
чистой водой также может наблюдаться
слабый результат реакции – бледное
окрашивание в случае цветных реакций,
слабое помутнение в реакциях осаждения
и др. Поэтому при анализе используют
так называемую «холостую
пробу
».

Физические
методы анализа – это методы, в которых
аналитический сигнал представляет
собой регистрируемую величину какого-либо
физического свойства вещества либо
результат взаимодействия излучения с
веществом.
В
физических методах анализа для получения
аналитического сигнала используется
не аналитическая реакция, но аналитический
процесс
.
Например, в ИК-спектроскопии под
аналитическим процессом следует понимать
процесс взаимодействия анализируемого
вещества с внешним по отношению к нему
источником электромагнитного излучения,
сопровождающийся поглощением
характеристических порций энергии –
квантов энергии. Обычно в физических
методах анализа для определения наличия
и величины аналитического сигнала
необходимо определенное устройство –
детектор, входящий в состав соответствующего
аналитического прибора. Физические
методы анализа являются инструментальными
методами.

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

Физико-химические
методы анализа включают элементы как
физических, так и химических методов
анализа.

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

Х
– аналитический
сигнал.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #

Макеты страниц

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

Аналитический сигнал. Формула Эйлера

представляющая гармоническое колебание в виде суммы даух комплексно-сопряженных функций, наводит на мысль о том, что произвольный сигнал s(t) с известной спектральной плотностью можно записать как сумму двух составляющих, каждая из которых содержит или только положительные, или только отрицательные частоты

Назовем функцию

аналитическим сигналом, отвечающим вещественному колебанию s(t). Первый из интегралов в правой части формулы (5,37) путем замены переменной преобразуется к виду

Поэтому формула (5.37) устанавливает связь между сигналами или

Мнимая часть аналитического сигнала

называется сопряженным сигналом по отношению к исходному колебанию s(t). Итак, аналитический сигнал

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

Введение аналитического и сопряженного сигналов, безусловно, не позволяет подучить каких-либо новых сведений, которые не содержались бы в математической модели сигнала s(t). Однако эти новые понятия открывают прямой путь к созданию систематических методов исследования узкополосных колебаний.

На конкретном примере покажем способ вычисления аналитического сигнала по известному спектру исходного сигнала.


Пример 5.6. Пусть – идеальный низкочастотный сигнал с известными параметрами (см. § 5.1).

В этом случае аналитический сигнал

Выделяя вещественную и мнимую части, получаем

Графики этих двух сигналов приведены на рис. 6,3.

Рис. 5.3. Исходный и сопряженный сигналы: 1 — идеальный низкочастотный сигнал; 2 — сопряженный с ним сигнал


Спектральная плотность аналитического сигнала.

Исследуем спектральную плотность аналитического сигнала, т. е. функцию связанную с -преобразованием Фурье:

На основании формулы (5.38) можно утверждать, что эта функция отлична от нуля лишь в области положительных частот:

Если — спектральная плотность сопряженного сигнала, то в силу линейности преобразования Фурье

Поэтому равенство (5.42) будет выполняться только в случае, когда спектральные плотности исходного и сопряженного сигналов связаны между собой следующим образом:

Абстрактно можно представить себе такой способ получения сопряженного сигнала: исходное колебание подается на вход некоторой системы, которая осуществляет поворот фаз всех спектральных составляющих на угол —90° в области положительных частот и на угол 90° в области отрицательных частот, не изменяя их по амплитуде. Систему, обладающую подобными свойствами, называют квадратурным фильтром.

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

Формула (5.44) показывает, что спектральная плотность сопряженного сигнала есть произведение спектра исходного сигнала и функции — . Поэтому сопряженный сигнал представляет собой свертку даух функций: , которая является обратным преобразованием Фурье по отношению к функции .

Для удобства вычислений представим эту функцию в виде предела:

Тогда

Таким образом, сопряженный сигнал связан с исходным сигналом соотношением

Можно поступить и по-иному, выразив сигнал через который полагается известным. Для этого достаточно заметить, что из (5.44) вытекает следующая связь между спектральными плотностями:

Поэтому соответствующая формула будет отличаться от (5.45) лишь знаком:

Формулы (5.45) и (5.46) известны в математике под названием прямого и обратного преобразований Гильберта.

Символическая запись их такова:

Поскольку функция называемая ядром этих преобразований, имеет разрыв при интегралы (5.45) и (5.46) следует понимать в смысле главного значения. Например:

Некоторые свойства преобразований Гильберта.

Простейшее свойство этих интегральных преобразований — их линейность:

при любых постоянных в чем можно убедиться непосредственно.

Ядро преобразования Гильберта есть нечетная функций аргумента относительно точки а, значит, сигнал, сопряженный к константе, тождественно равен нулю:

Важное свойство преобразования Гильберта состоит в следующем: если при каком-либо t исходный сигнал s(t) достигает экстремума (максимума или минимума), то в окрестности этой точки сопряженный сигнал проходит через нуль. Чтобы убедиться в этом, нужно на одном чертеже совместить графики s(t) и ядра . Пусть значение t близко к тому при котором функция экстремальна. Поскольку сигнал является здесь четной функцией, а ядро нечетной, будет наблюдаться компенсация площадей фигур, ограниченных горизонтальной осью и кривой, которая описывает подынтегральную функцию преобразования Гильберта. Образно говоря, если исходный сигнал изменяется во времени «подобно косинусу», то сопряженный с ним сигнал будет изменяться «подобно синусу».

Отметим, что преобразования Гильберта имеют нелокальный характер: подведение сопряженного сигнала в окрестности какой-либо точки зависит от свойств исходного сигнала на всей оси времени, хотя наибольший вклад дает, конечно, достаточно близкая окрестность рассматриваемой точки.

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

Вычислим сигналы, сопряженные с гармоническими колебаниями и Результаты можно получить непосредственно из формулы (5.45). Однако проще поступить таким образом. Пусть некоторый произвольный сигнал задан своим Фурье-представлением:

На основании соотношения (5.44) находим аналогичное представление сопряженного сигнала:

Рассматривая формулы (5.48) и (5.49) совместно, находим следующие законы преобразования Гильберта:

Преобразование Гильберта для узкополосного сигнала

Пусть известна функция — спектральная плотность комплексной огибающей узкополосного сигнала s(t) с опорной частотой . Согласно формуле (5.36), спектр данного сигнала

Первое слагаемое в правой части соответствует области частот второе — Тогда на основании формулы (5.44) спектр сопряженного сигнала

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

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

то в соответствии с равенством (5.53) комплексная огибающая сопряженного сигнала

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

Отсюда следует, что узкополосному сигналу

соответствует сопряженный по Гильберту сигнал

Вычисление огибающей, полной фазы и мгновенной частоты.

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

Целесообразность такого определения можно проверить на примере узкополосного сигнала. Используя формулы (5.54) и (5.55), находим, что огибающая такого сигнала

В § 5.3 данная формула была получена из других соображений.

По определению, полная фаза любого сигнала равна аргументу аналитического сигнала

(5-57)

Наконец, мгновенная частота сигнала есть производная полной фазы по времени:

Рассмотрим примеры, иллюстрирующие вычисление указанных характеристик узкополосных сигналов.


Пример 5.7. Дано простое гармоническое колебание

В этом случае сопряженный сигнал Огибающая исходного сигнала

естественно, не зависит от времени и равна его амплитуде.

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

Пример 5.8. Колебание s(t) является суммой двух гармонических составляющих с различными амплитудами и частотами:

Поскольку

огибающая такого сигнала изменяется во времени по закону

Полная фаза сигнала

Для вычисления мгновенной частоты сяедует воспользоваться формулой (5.58), которая приводит к следующему результату:

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

Пример 5.9. Рассмотрим идеальный полосовой сигнал s(t), спектр которого при отличен от нуля лишь на отрезке

Соответствующий аналитический сигнал

Огибающая исходного полосового сигнала

Наконец, мгновенная частота сигнала

Выполнив несложные преобразования, находим, что в данном случае не зависит от времени и равна центральной частоте интервала, в котором сосредоточен спектр.


Итак, зная аналитический сигнал, можно однозначно определять огибающую и мгновенную частоту узкополосного колебания, не применяя несколько искусственное понятие опорной частоты. Более того, формулы (5.56)-(5.58) сохраняют смысл применительно к сигналам произвольного вида, не обязательно удовлетворяющим условиям квазигармоничности (узкополосности).

Заключительные замечания.

Теория аналитического сигнала применительно к задачам теории колебаний и волн была развита в 40-х годах в работах Габора [30]. Однако преобразования Гильберта появились в математике еще в начале XX в. в связи с так называемой краевой задачей теории аналитических функций [10]. Сущность этой задачи состоит в следующем.

Пусть — комплексная переменная, — функция, аналитическая в верхней полуплоскости, т. е. при На вещественной оси, являющейся границей области аналитичности, функция имеет как вещественную, так и мнимую части:

Требуется найти закон, связывающий между собой функции .

Решение задачи дается преобразованиями Гильберта:

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

В последнее время методы, основанные на понятиях аналитического сигнала и преобразований Гильберта, прочно вошли в арсенал теоретической радиотехники. Некоторые интересные проблемы в этой области описаны в [261.

1

Оглавление

  • Предисловие ко второму изданию
  • Введение
  • 1. Радиотехнические сигналы
  • 1.1. Классификация радиотехнических сигналов
  • 1.2. Динамическое представление сигналов
  • 1.3. Геометрические методы в теории сигналов
  • 1.4. Теория ортогональных сигналов
  • Результаты
  • Глава 2. Спектральные представления сигналов
  • 2.1. Периодические сигналы и ряды Фурье
  • 2.2. Спектральный анализ непериодических сигналов. Преобразование Фурье
  • 2.3. Основные свойства преобразования Фурье
  • 2.4. Спектральные плотности неинтегрируемых сигналов
  • 2.5. Преобразование Лапласа
  • Результаты
  • Глава 3. Энергетические спектры сигналов. Принципы корреляционного анализа
  • 3.1. Взаимная спектральная плотность сигналов. Энергетический спектр
  • 3.2. Корреляционный анализ сигналов
  • 3.3. Автокорреляционная функция дискретного сигнала
  • 3.4. Взаимокорреляционная функция двух сигналов
  • Результаты
  • Глава 4. Модулированные сигналы
  • 4.1. Сигналы с амплитудной модуляцией
  • 4.2. Сигналы с угловой модуляцией
  • 4.3. Сигналы с внутриимпульсной частотной модуляцией
  • Результаты
  • Глава 5. Сигналы с ограниченным спектром
  • 5.1. Некоторые математические модели сигналов с ограниченным спектром
  • 5.2. Теорема Котельникова
  • 5.3. Узкополосные сигналы
  • 5.4. Аналитический сигнал и преобразование Гильберта
  • Результаты
  • Глава 6. Основы теории случайных сигналов
  • 6.1. Случайные величины и их характеристики
  • 6.2. Статистические характеристики систем случайных величин
  • 6.3. Случайные процессы
  • Результаты
  • Глава 7. Корреляционная теория случайных процессов
  • 7.1. Спектральные представления стационарных случайных процессов
  • 7.2. Дифференцирование и интегрирование случайных процессов
  • § 7.3. Узкополосные случайные процессы
  • Результаты
  • 2. Радиотехнические цепи, устройства и системы
  • 8.1. Физические системы и их математические модели
  • 8.2. Импульсные, переходные и частотные характеристики линейных стационарных систем
  • 8.3. Линейные динамические системы
  • 8.4. Спектральный метод
  • 8.5. Операторный метод
  • Результаты
  • Глава 9. Воздействие детерминированных сигналов на частотно-избирательные системы
  • 9.1. Некоторые модели частотно-избирательных цепей
  • 9.2. Частотно-избирательные цепи при иирокополосных входных воздействиях
  • 9.3. Частотно-избирательные цепи при узкополосных входных воздействиях
  • Результаты
  • Глава 10. Воздействие случайных сигналов на линейные стационарные цепи
  • 10.1. Спектральный метод анализа воздействия случайных сигналов на линейные стационарные цепи
  • 10.2. Источники флуктуационных шумов в радиотехнических устройствах
  • Результаты
  • Глава 11. Преобразования сигналов в нелинейных радиотехнических цепях
  • 11.1. Безынерционные нелинейные преобразования
  • 11.2. Спектральный состав тока в безынерционном нелинейном элементе при гармоническом внешнем воздействии
  • 11.3. Нелинейные резонансные усилители и умножители частоты
  • 11.4. Безынерционные нелинейные преобразования суммы нескольких гармонических сигналов
  • 11.5. Получение модулированных радиосигналов
  • 11.6. Амплитудное, фазовое и частотное детектирование
  • 11.7. Воздействие стационарных случайных сигналов на безынерционные нелинейные цепи
  • Результаты
  • Глава 12. Преобразования сигналов в линейных параметрических цепях
  • 12.1. Прохождение сигналов через резистивные параметрические цепи
  • 12.2. Энергетические соотношения в параметрических реактивных элементах цепи
  • 12.3. Принципы параметрического усиления
  • 12.4. Воздействие гармонических сигналов на параметрические системы со случайными характеристиками
  • Результаты
  • Глава 13. Элементы теории синтеза линейных частотных фильтров
  • 13.1. Частотные характеристики четырехполюсников
  • 13.2. Фильтры нижних частот
  • 13.3. Реализация фильтров
  • Результаты
  • Глава 14. Активные цепи с обратной связью и автоколебательные системы
  • 14.1. Передаточная функция линейной системы с обратной связью
  • 14.2. Устойчивость цепей с обратной связью
  • 14.3. Активные RС-фильтры
  • 14.4. Автогенераторы гармонических колебаний. Режим малого сигнала
  • 14.5. Автогенераторы гармонических колебаний. Режим большого сигнала
  • Результаты
  • Глава 15. Дискретные сигналы. Принципы цифровой фильтрации
  • 15.1. Модели дискретных сигналов
  • 15.2. Дискретизация периодических сигналов
  • 15.3. Теория z-преобразования
  • 15.4. Цифровые фильтры
  • 15.5. Реализация алгоритмов цифровой фильтрации
  • 15.6. Синтез линейных цифровых фильтров
  • Результаты
  • Глава 16. Некоторые вопросы теории помехоустойчивости радиоприема
  • 16.1. Выделение полезного сигнала с помощью линейного частотного фильтра
  • 16.2. Оптимальная линейная фильтрация сигналов известной формы
  • 16.3. Реализация согласованных фильтров
  • 16.4. Оптимальная фильтрация случайных сигналов
  • 16.5. Сравнение помехоустойчивости радиосистем с амплитудной и частотной модуляцией
  • Результаты
  • Заключение
  • Приложения
  • 1. Функции Уолша и их некоторые свойства
  • Список рекомендуемой литературы

Аналитический сигнал (аналитическое представление сигнала) — используемое в теории обработки сигналов математическое представление аналогового сигнала в виде комплекснозначной аналитической функции времени. Обычный, действительный сигнал x является при этом действительной частью аналитического представления xa.

Идея преобразования — оставить лишь неотрицательные частоты в спектре сигнала, достаточные для его восстановления в силу эрмитовой симметрии: X(-f) = overline{X(f)}.

Аналитический сигнал является обобщением понятия комплексной амплитуды на случай сигналов, отличных от гармонического.

Определение[править | править код]

Пусть x(t) — представляющая сигнал действительнозначная функция, преобразование Фурье (т.е. спектр) которой обозначим X(f),[1], а u(f) — функция Хевисайда.

Тогда:


begin{align}
X_mathrm{a}(f) &  stackrel{mathrm{def}}{=}  
begin{cases}
  2 X(f), & mbox{for } f > 0, \
  X(f), & mbox{for } f = 0, \
  0, & mbox{for } f < 0,
end{cases} \
&= X(f)cdot 2 mathrm{u}(f)
end{align}

содержит лишь неотрицательную часть спектра X(f).

Подвергая полученный спектр {displaystyle X_{mathrm {a} }(f)} обратному преобразованию Фурье, мы и получаем аналитический сигнал:


begin{align}
x_mathrm{a}(t) &= underbrace{mathcal{F}^{-1}{X(f)}}_{x(t)}  *  underbrace{mathcal{F}^{-1}{2 mathrm{u}(f)}}_{delta(t) + jcdot {1 over pi t}} quad quad mbox{ } \
&= x(t) + junderbrace{left[x(t) * {1 over pi t}right]}_{hat{x}(t)},
end{align}

где * — свёртка, hat{x}(t) — преобразование Гильберта функции x(t), а j означает мнимую единицу.

Примеры[править | править код]

Пусть  {displaystyle x(t)=cos(omega _{0}t)} для некоторой частоты {displaystyle omega _{0}>0}

Тогда:

{displaystyle {hat {x}}(t)=cos(omega _{0}t-{begin{matrix}{frac {pi }{2}}end{matrix}})=sin(omega _{0}t)}
{displaystyle x_{mathrm {a} }(t)=cos(omega _{0}t)+jcdot sin(omega _{0}t)=e^{jomega _{0}t}}

Это комплексная функция с возрастающим по времени аргументом.

Практические применения[править | править код]

Устранение «отрицательных частот» используется в аналоговой передаче звука (АМ-радиовещание, аналоговая телефонная связь) для экономии полосы частот.

Примечания[править | править код]

  1. Обращаем внимание на то, что fпеременная частоты, а не функция.

См. также[править | править код]

  • Комплексное сопротивление

Внешние ссылки[править | править код]

  • Преобразование Гильберта. Аналитический сигнал.
  • Преобразование Гильберта. Аналитический сигнал.

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

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

Пока управляемая функция не имеет отрицательных частотных составляющих (то есть она все еще аналитическая), преобразование из комплексного обратно в реальное – это просто вопрос отбрасывания мнимой части. Аналитическое представление является обобщением концепции фазора : в то время как вектор ограничен неизменными во времени амплитудой, фазой и частотой, аналитический сигнал допускает параметры, изменяющиеся во времени.

Содержание

  • 1 Определение
    • 1.1 Отрицательные частотные составляющие
  • 2 Примеры
    • 2.1 Пример 1
    • 2.2 Пример 2
    • 2.3 Пример 3
  • 3 Свойства
    • 3.1 Мгновенная амплитуда и фаза
    • 3.2 Комплексная огибающая / основная полоса
  • 4 Расширения аналитического сигнала до сигналов нескольких переменных
    • 4.1 Многомерный аналитический сигнал на основе специального направления
    • 4.2 Моногенный сигнал
  • 5 См. Также
    • 5.1 Приложения
  • 6 Примечания
  • 7 Ссылки
  • 8 Дополнительная литература
  • 9 Внешние ссылки

Определение

Передаточная функция для создания аналитического сигнала

Если s (t) { displaystyle s (t)}s (t) – функция с действительным знаком с преобразованием Фурье S (f) { displaystyle S (f)}S (f) , тогда преобразование имеет эрмитову симметрию относительно оси f = 0 { displaystyle f = 0}f = 0 :

S (- f) = S (f) ∗, { Displaystyle S (-f) = S (f) ^ {*},}S (-f) = S (f) ^ {*},

где S (f) ∗ { displaystyle S (f) ^ {*}}S (е) ^ {*} – это комплексно сопряженное из S (f) { di splaystyle S (f)}S (f) . Функция:

S a (f) ≜ {2 S (f), для f>0, S (f), для f = 0, 0, для f < 0 = 2 u ⁡ ( f) ⏟ 1 + sgn ⁡ ( f) S ( f) = S ( f) + sgn ⁡ ( f) S ( f), {displaystyle {begin{aligned}S_{mathrm {a} }(f)triangleq {begin{cases}2S(f),{text{for}} f>0, \ S (f), { text {for}} f = 0, \ 0, { text {for}} f <0end{cases}}\=underbrace {2operatorname {u} (f)} _{1+operatorname {sgn}(f)}S(f)=S(f)+operatorname {sgn}(f)S(f),end{aligned}}}{displaystyle {begin{aligned}S_{mathrm {a} }(f)triangleq {begin{cases}2S(f),{text{for}} f>0, \ S (f), { text {for}}  f = 0, \ 0, { text {for}}  f <0end{cases}}\=underbrace {2operatorname {u} (f)} _{1+operatorname {sgn}(f)}S(f)=S(f)+operatorname {sgn}(f)S(f),end{aligned}}}

, где

содержит только неотрицательные частотные компоненты S (е) { Displaystyle S (f)}S (f) . И операция обратима из-за эрмитовой симметрии S (f) { displaystyle S (f)}S (f) :

S (f) = {1 2 S a (f), для f>0, S a (f), для f = 0, 1 2 S a (- f) ∗, для f < 0 (Hermitian symmetry) = 1 2 [ S a ( f) + S a ( − f) ∗ ]. {displaystyle {begin{aligned}S(f)={begin{cases}{frac {1}{2}}S_{mathrm {a} }(f),{text{for}} f>0, \ S _ { mathrm {a}} (f), { text {for}} f = 0, \ { frac {1} {2}} S _ { mathrm {a}} (- f) ^ {*}, { text {for}} f <0 {text{(Hermitian symmetry)}}end{cases}}\={frac {1}{2}}[S_{mathrm {a} }(f)+S_{mathrm {a} }(-f)^{*}].end{aligned}}}{begin{aligned}S(f)={begin{cases}{frac {1}{2}}S_{{mathrm {a}}}(f),{text{for}} f>0, \ S _ {{ mathrm {a}}} (f), { text {for}}  f = 0, \ { frac {1} {2}} S _ {{ mathrm {a}}} ( -f) ^ {*}, { text {for}}  f <0 {text{(Hermitian symmetry)}}end{cases}}\={frac {1}{2}}[S_{{mathrm {a}}}(f)+S_{{mathrm {a}}}(-f)^{*}].end{aligned}}

. Аналитический сигнал из s (t) { displaystyle s (t)}s (t) – это обратное преобразование Фурье для S a (f) { displaystyle S _ { mathrm {a}} (f)}S _ {{ mathrm {a}}} (е) :

sa (t) ≜ F – 1 [S a (f)] = F – 1 [S (f) + sign ⁡ (f) ⋅ S (f)] = F – 1 {S (f)} ⏟ s (t) + F – 1 {sign ⁡ (f)} ⏟ j 1 π t ∗ F – 1 {S (f)} ⏟ s (t) ⏞ свертка = s (t) + j [1 π t ∗ s (t)] ⏟ H ⁡ [s (t)] = s (t) + js ^ (t), { displaystyle { begin {align} s _ { mathrm {a}} (t) треугольникq { mathcal {F}} ^ {- 1} [S _ { mathrm {a}} (f)] \ = { mathcal {F}} ^ {- 1} [S (f) + operatorname {sgn} (f) cdot S ( f)] \ = underbrace {{ mathcal {F}} ^ {- 1} {S (f) }} _ {s (t)} + overbrace { underbrace {{ mathcal {F }} ^ {- 1} { operatorname {sgn} (f) }} _ {j { frac {1} { pi t}}} * underbrace {{ mathcal {F}} ^ {- 1} {S (f) }} _ {s (t)}} ^ { text {convolution}} \ = s (t) + j underbrace { left [{1 over pi t } * s (t) right]} _ { operatorname { mathcal {H}} [s (t)]} \ = s (t) + j { hat {s}} (t), конец {выровнен}}}{ displaystyle { begin {align} s _ { mathrm {a}} (t)  треугольникq { mathcal {F}} ^ {- 1} [ S _ { mathrm {a}} (f)] \ = { mathcal {F}} ^ {- 1} [S (f) +  operatorname {sgn} (f)  cdot S (f)]   =  underbrace {{ mathcal {F}} ^ {- 1}  {S (f) }} _ {s (t)} +  overbrace { underbrace {{ mathcal {F}} ^ { -1}  { operatorname {sgn} (f) }} _ {j { frac {1} { pi t}}} *  underbrace {{ mathcal {F}} ^ {- 1}  {S (f) }} _ {s (t)}} ^ { text {convolution}} \ = s (t) + j  underbrace { left [{1  over  pi t} * s (t)  right]} _ { operatorname { mathcal {H}} [s (t)]} \ = s (t) + j { hat {s}} (t),  end {align}}}

где

. Отметим, что s (t) = s (t) ∗ δ (t), { displaystyle s (t) = s (t) * delta (t),}{ displaystyle s (t) = s (t) *  delta (t),} это также можно выразить как операцию фильтрации, которая непосредственно удаляет отрицательные частотные составляющие :

sa (t) = s (t) ∗ [δ (t) + j 1 π t] ⏟ F – 1 {2 u (f)}. { Displaystyle s _ { mathrm {a}} (t) = s (t) * underbrace { left [ delta (t) + j {1 over pi t} right]} _ {{ mathcal {F}} ^ {- 1} {2u (f) }}.}{ displaystyle s _ { mathrm {a}} (t) = s (t) *  underbrace { left [ delta (t) + j {1  over  pi t}  right]} _ {{ mathcal {F }} ^ {- 1}  {2u (f) }}.}

Отрицательные частотные компоненты

Так как s (t) = Re ⁡ [sa (t)] { displaystyle s (t) = operatorname {Re} [s _ { mathrm {a}} (t)]}s (t) =  operatorname {Re} [s _ {{ mathrm {a}}} (t)] , восстановление отрицательных частотных составляющих – это простой вопрос отбрасывания Im ⁡ [ sa (t)] { displaystyle operatorname {Im} [s _ { mathrm {a}} (t)]} operatorname {Im} [s _ {{ mathrm {a}}} (t)] , что может показаться нелогичным. Мы также можем отметить, что комплексное сопряжение s a ∗ (t) { displaystyle s _ { mathrm {a}} ^ {*} (t)}s _ {{ mathrm {a}}} ^ {*} (t) содержит только отрицательные частотные компоненты. И поэтому s (t) = Re ⁡ [sa ∗ (t)] { displaystyle s (t) = operatorname {Re} [s _ { mathrm {a}} ^ {*} (t)]}{ displaystyle s (t) =  operatorname {Re} [s _ { mathrm {a}} ^ {*} (t)]} восстанавливает подавленные положительные частотные составляющие. Другая точка зрения состоит в том, что мнимая составляющая в любом случае – это член, который вычитает частотные составляющие из s (t). Оператор Re { displaystyle operatorname {Re}} operatorname {Re} удаляет вычитание, создавая впечатление добавления новых компонентов.

Примеры

Пример 1

s (t) = cos ⁡ (ω t), { displaystyle s (t) = cos ( omega t),}s (t) =  cos ( omega t), где ω>0. { displaystyle omega>0.}omega>0.

Тогда:

s ^ (t) = cos ⁡ (ω t – π 2) = sin ⁡ (ω t), { displaystyle { hat {s}} ( t) = cos left ( omega t – { frac { pi} {2}} right) = sin ( omega t),}{ displaystyle { hat {s}} (t) =  cos  left ( omega t - { frac { pi} {2}}  right) =  sin ( omega t),}
sa (t) = s (t) + js ^ (T) знак равно соз ⁡ (ω T) + J грех ⁡ (ω T) = Ej ω T. { Displaystyle s _ { mathrm {a}} (t) = s (t) + j { hat {s }} (t) = cos ( omega t) + j sin ( omega t) = e ^ {j omega t}.}s _ {{ mathrm {a}}} (t) = s (t) + j { hat {s}} (t) =  cos ( omega t) + j  sin ( omega t) = e ^ {{j  omega t}}. Третье равенство – это формула Эйлера.

A следствие из формулы Эйлера : cos ⁡ (ω t) = 1 2 (ej ω t + ej (- ω) t). { displaystyle cos ( omega t) = { tfrac {1} {2}} left (e ^ {j omega t} + e ^ {j (- omega) t} right).}{ displaystyle  cos ( omega t) = { tfrac {1} {2}}  left (e ^ {j  omega t} + e ^ {j (-  omega) t}  right).} В общем, аналитический представление простой синусоиды получается выражением ее в терминах комплексных экспонент, отбрасывая отрицательную частотную составляющую и удваивая положительную частотную составляющую энт. А аналитическое представление суммы синусоид – это сумма аналитических представлений отдельных синусоид.

Пример 2

Здесь мы используем формулу Эйлера, чтобы идентифицировать и отбросить отрицательную частоту.

s (t) знак равно соз ⁡ (ω t + θ) = 1 2 (ej (ω t + θ) + e – j (ω t + θ)) { displaystyle s (t) = cos ( омега t ​​+ theta) = { frac {1} {2}} left (e ^ {j ( omega t + theta)} + e ^ {- j ( omega t + theta)} right)}{ displaystyle s ( t) =  cos ( omega t +  theta) = { frac {1} {2}}  left (e ^ {j ( omega t +  theta)} + e ^ {- j ( omega t +  theta)}  right)}

Тогда:

sa (t) = {ej (ω t + θ) = ej | ω | t ⋅ e j θ, если ω>0, e – j (ω t + θ) = e j | ω | t ⋅ e – j θ, если ω < 0. {displaystyle s_{mathrm {a} }(t)={begin{cases}e^{j(omega t+theta)} = e^{j|omega |t}cdot e^{jtheta },{text{if}} omega>0, \ e ^ {- j ( omega t + theta)} = e ^ {j | omega | t} cdot e ^ {- j theta }, { text {if}} omega <0.end{cases}}}s_{{mathrm {a}}}(t)={begin{cases}e^{{j(omega t+theta)}}  = e^{{j|omega |t}}cdot e^{{jtheta }},{text{if}} omega>0, \ e ^ {{- j ( omega t +  theta)}} =  e ^ {{j |  omega | t}}  cdot e ^ {{- j  theta}}, { text {if}}   omega <0.end{cases}}

Пример 3

Это еще один пример использования метода преобразования Гильберта для удаления отрицательных частотных составляющих. Заметим, что ничто не мешает нам вычислить sa (t) { displaystyle s _ { mathrm {a}} (t)}s _ {{ mathrm {a}}} (t) для комплексного значения s (t) { displaystyle s (t)}s (t) . Но это может быть не обратимое представление, потому что исходный спектр в целом несимметричен. Таким образом, за исключением этого примера, в общем обсуждении предполагается, что вещественное значение s (t) { displaystyle s (t)}s (t) .

s (t) = e – j ω t { displaystyle s (t) = e ^ {- j omega t}}s (t) = e ^ {{- j  omega t} } , где ω>0 { displaystyle omega>0}omega>0 .

Затем:

s ^ (t) = je – j ω t, { displaystyle { hat {s}} (t) = je ^ {- j omega t},}{ hat {s}} (t) = je ^ {{ -j  omega t}},
sa (t) = e – j ω t + j 2 e – j ω t = e – j ω T – E – J ω T = 0. { Displaystyle s _ { mathrm {a}} (t) = e ^ {- j omega t} + j ^ {2} e ^ {- j omega t} = e ^ {- j omega t} -e ^ {- j omega t} = 0.}s _ {{ mathrm {a}}} (t) = e ^ {{- j  omega t}} + j ^ {2} e ^ {{- j  omega t}} = e ^ {{- j  omega t}} - e ^ {{- j  omega t}} = 0.

Свойства

Мгновенные амплитуда и фаза

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

Аналитический сигнал также может быть выражен в полярных координатах :

sa (t) = sm (t) ej ϕ (t), { displaystyle s _ { mathrm {a}} (t) = s _ { mathrm {m}} (t) e ^ {j phi (t)},}s _ {{ mathrm {a}}} (t) = s _ {{ mathrm {m}}} (t) e ^ {{j  phi (t)}},

где вводятся следующие изменяющиеся во времени величины:

На прилагаемой диаграмме синяя кривая изображает s (t) { displaystyle s (t)}s (t) и красная кривая изображает соответствующую см (t) { displaystyle s _ { mathrm {m}} (t)}s _ {{ mathrm {m}}} (t) .

Производная по времени от развернутой мгновенной фазы имеет единицы радиан / секунду, и называется мгновенной угловой частотой:

ω (t) ≜ d ϕ dt (t). { displaystyle omega (t) треугольникq { frac {d phi} {dt}} (t).}{ displaystyle  omega (t)  треугольникq { frac {d  phi} {dt}} (t).}

мгновенная частота (в герцах ) следовательно:

f (t) ≜ 1 2 π ω (t). { displaystyle f (t) Triangleq { frac {1} {2 pi}} omega (t).}{ displaystyle f (t)  треугольникq { frac {1} {2  pi}}  omega (t).}

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

Комплексная огибающая / основная полоса

Аналитические сигналы часто смещены по частоте (преобразованы с понижением частоты) в сторону 0 Гц, что может создавать [несимметричные] отрицательные частотные компоненты:

sa ↓ (t) ≜ sa (t) е – j ω 0 T знак равно sm (t) ej (ϕ (t) – ω 0 t), { displaystyle {s _ { mathrm {a}}} _ { downarrow} (t) треугольник q s _ { mathrm {a}} (t) e ^ {- j omega _ {0} t} = s _ { mathrm {m}} (t) e ^ {j ( phi (t) – omega _ {0} t)},}{ displaystyle { s _ { mathrm {a}}} _ { downarrow} (t)  треугольник s _ { mathrm {a}} (t) e ^ {- j  omega _ {0} t} = s _ { mathrm {m }} (т) е ^ {j ( phi (t) -  omega _ {0} t)},}

где ω 0 { displaystyle omega _ {0}} omega _ {0} – произвольная опорная угловая частота.

Эта функция имеет различные названия, такие как сложный конверт и сложный baseband. Сложный конверт не уникален; он определяется выбором ω 0 { displaystyle omega _ {0}} omega _ {0} . Эта концепция часто используется при работе с сигналами полосы пропускания. Если s (t) { displaystyle s (t)}s (t) является модулированным сигналом, ω 0 { displaystyle omega _ {0}} omega _ {0} может быть приравнивается к его несущей частоте.

. В других случаях ω 0 { displaystyle omega _ {0}} omega _ {0} выбирается где-то посередине желаемой полосы пропускания. Тогда простой фильтр нижних частот с действительными коэффициентами может вырезать интересующую часть. Другой мотив – снизить максимальную частоту, что снижает минимальную скорость выборки без псевдонимов. Сдвиг частоты не подрывает математическую управляемость представления комплексного сигнала. Таким образом, в этом смысле сигнал, преобразованный с понижением частоты, остается аналитическим. Однако восстановление представления с действительным знаком – это уже не просто извлечение реального компонента. Может потребоваться преобразование с повышением частоты, и если сигнал был дискретизирован (дискретное время), также может потребоваться интерполяция (повышающая дискретизация ), чтобы избежать алиасинг.

Если ω 0 { displaystyle omega _ {0}} omega _ {0} выбрано больше максимальной частоты sa (t), { displaystyle s _ { mathrm {a}} (t),}s _ {{ mathrm {a}}} (t), , затем sa ↓ (t) { displaystyle {s _ { mathrm {a}}} _ { downarrow} (t)}{ displaystyle {s _ { mathrm {a}}} _ { downarrow} (t)} не имеет положительных частот. В этом случае извлечение реального компонента восстанавливает их, но в обратном порядке; низкочастотные компоненты теперь высокие, и наоборот. Это может использоваться для демодуляции типа сигнала одной боковой полосы, называемого нижней боковой полосой или инвертированной боковой полосой.

Другие варианты задания частоты

Иногда ω 0 { displaystyle омега _ {0}} omega _ {0} выбрано, чтобы минимизировать

∫ 0 + ∞ (ω – ω 0) 2 | S a (ω) | 2 d ω. { displaystyle int _ {0} ^ {+ infty} ( omega – omega _ {0}) ^ {2} | S _ { mathrm {a}} ( omega) | ^ {2} , d omega.}{ displaystyle  int _ {0} ^ {+  infty} ( omega -  omega _ {0}) ^ {2} | S _ { mathrm {a}} ( omega) | ^ {2} , d  omega.}

В качестве альтернативы можно выбрать ω 0 { displaystyle omega _ {0}} omega _ {0} , чтобы минимизировать среднеквадратичную ошибку при линейной аппроксимации развернутой мгновенной фазы ϕ (t) { displaystyle phi (t)} phi (t) :

∫ – ∞ + ∞ [ω (t) – ω 0] 2 | s a (t) | 2 dt { displaystyle int _ {- infty} ^ {+ infty} [ omega (t) – omega _ {0}] ^ {2} | s _ { mathrm {a}} (t) | ^ {2} , dt} int _ {{-  infty}} ^ {{+  infty}} [ omega (t) -  omega _ {0}] ^ {2} | s _ {{ mathrm {a} }} (t) | ^ {2} , dt

или другая альтернатива (для некоторого оптимального θ { displaystyle theta} theta ):

∫ – ∞ + ∞ [ϕ (t) – (ω 0 t + θ)] 2 dt. { displaystyle int _ {- infty} ^ {+ infty} [ phi (t) – ( omega _ {0} t + theta)] ^ {2} , dt.} int _ {{-  infty}} ^ {{+  infty}} [ phi (t) - ( om ega _ {0} t +  theta)] ^ {2} , dt.

В В области частотно-временной обработки сигналов было показано, что аналитический сигнал был необходим для определения распределения Вигнера – Вилля, чтобы метод мог иметь желаемые свойства, необходимые для практических приложений.

Иногда фраза «комплексная огибающая» имеет более простой смысл комплексной амплитуды вектора (постоянной частоты); в других случаях комплексная огибающая s m (t) { displaystyle s_ {m} (t)}{ displaystyle s_ { m} (t)} , как определено выше, интерпретируется как зависящее от времени обобщение комплексной амплитуды. Их соотношение мало чем отличается от реального случая: изменяющаяся огибающая обобщающая константа амплитуда.

Расширение аналитического сигнала на сигналы нескольких переменных

Концепция аналитического signal четко определен для сигналов одной переменной, которая обычно является временем. Для сигналов двух или более переменных аналитический сигнал может быть определен по-разному, и ниже представлены два подхода.

Многомерный аналитический сигнал на основе специального направления

Прямое обобщение аналитического сигнала может быть выполнено для многомерного сигнала после того, как установлено, что подразумевается под отрицательными частотами для Это дело. Это можно сделать, введя единичный вектор u ^ { displaystyle { boldsymbol { hat {u}}}}{ boldsymbol { hat {u}}} в область Фурье и пометив любой частотный вектор ξ { displaystyle { boldsymbol { xi}}}{ boldsymbol  xi} как отрицательное, если ξ ⋅ u ^ < 0 {displaystyle {boldsymbol {xi }}cdot {boldsymbol {hat {u}}}<0}{ boldsymbol  xi}  cdot { boldsymbol { hat {u}}} <0 . Затем получают аналитический сигнал путем удаления всех отрицательных частот и умножения результата на 2 в соответствии с процедурой, описанной для случая сигналов с одной переменной. Однако нет конкретного направления для u ^ { displaystyle { boldsymbol { hat {u}}}}{ boldsymbol { hat {u}}} , которое следует выбрать, если нет дополнительных ограничений. Следовательно, выбор u ^ { displaystyle { boldsymbol { hat {u}}}}{ boldsymbol { hat {u}}} является специальным или зависит от приложения.

Моногенный сигнал

Реальная и мнимая части аналитического сигнала соответствуют двум элементам векторно-значного сигнала, как это определено для сигналов с одной переменной. Однако моногенный сигнал может быть расширен до произвольного числа переменных прямым способом, создавая (n + 1) -мерную векторнозначную функцию для случая сигналов с n переменными.

См. Также

Приложения

Примечания

Ссылки

Дополнительная литература

  • Леон Коэн, Частотно-временной анализ, Prentice Hall, Upper Saddle River, 1995.
  • Frederick W. King, Hilbert Transforms, vol. II, Cambridge University Press, Кембридж, 2009.
  • B. Боашаш, Частотно-временной анализ и обработка сигналов: исчерпывающий справочник, Elsevier Science, Oxford, 2003.

Внешние ссылки

From Wikipedia, the free encyclopedia

In mathematics and signal processing, an analytic signal is a complex-valued function that has no negative frequency components.[1]  The real and imaginary parts of an analytic signal are real-valued functions related to each other by the Hilbert transform.

The analytic representation of a real-valued function is an analytic signal, comprising the original function and its Hilbert transform. This representation facilitates many mathematical manipulations. The basic idea is that the negative frequency components of the Fourier transform (or spectrum) of a real-valued function are superfluous, due to the Hermitian symmetry of such a spectrum. These negative frequency components can be discarded with no loss of information, provided one is willing to deal with a complex-valued function instead. That makes certain attributes of the function more accessible and facilitates the derivation of modulation and demodulation techniques, such as single-sideband.

As long as the manipulated function has no negative frequency components (that is, it is still analytic), the conversion from complex back to real is just a matter of discarding the imaginary part. The analytic representation is a generalization of the phasor concept:[2] while the phasor is restricted to time-invariant amplitude, phase, and frequency, the analytic signal allows for time-variable parameters.

Definition[edit]

Transfer function to create an analytic signal

If s(t) is a real-valued function with Fourier transform S(f), then the transform has Hermitian symmetry about the f=0 axis:

S(-f)=S(f)^{*},

where S(f)^{*} is the complex conjugate of S(f).
The function:

{displaystyle {begin{aligned}S_{mathrm {a} }(f)&triangleq {begin{cases}2S(f),&{text{for}} f>0,\S(f),&{text{for}} f=0,\0,&{text{for}} f<0end{cases}}\&=underbrace {2operatorname {u} (f)} _{1+operatorname {sgn}(f)}S(f)=S(f)+operatorname {sgn}(f)S(f),end{aligned}}}

where

contains only the non-negative frequency components of S(f). And the operation is reversible, due to the Hermitian symmetry of S(f):

{begin{aligned}S(f)&={begin{cases}{frac  {1}{2}}S_{{mathrm  {a}}}(f),&{text{for}} f>0,\S_{{mathrm  {a}}}(f),&{text{for}} f=0,\{frac  {1}{2}}S_{{mathrm  {a}}}(-f)^{*},&{text{for}} f<0 {text{(Hermitian symmetry)}}end{cases}}\&={frac  {1}{2}}[S_{{mathrm  {a}}}(f)+S_{{mathrm  {a}}}(-f)^{*}].end{aligned}}

The analytic signal of s(t) is the inverse Fourier transform of S_{{mathrm  {a}}}(f):

{displaystyle {begin{aligned}s_{mathrm {a} }(t)&triangleq {mathcal {F}}^{-1}[S_{mathrm {a} }(f)]\&={mathcal {F}}^{-1}[S(f)+operatorname {sgn}(f)cdot S(f)]\&=underbrace {{mathcal {F}}^{-1}{S(f)}} _{s(t)}+overbrace {underbrace {{mathcal {F}}^{-1}{operatorname {sgn}(f)}} _{j{frac {1}{pi t}}}*underbrace {{mathcal {F}}^{-1}{S(f)}} _{s(t)}} ^{text{convolution}}\&=s(t)+junderbrace {left[{1 over pi t}*s(t)right]} _{operatorname {mathcal {H}} [s(t)]}\&=s(t)+j{hat {s}}(t),end{aligned}}}

where

Noting that {displaystyle s(t)=s(t)*delta (t),} this can also be expressed as a filtering operation that directly removes negative frequency components:

{displaystyle s_{mathrm {a} }(t)=s(t)*underbrace {left[delta (t)+j{1 over pi t}right]} _{{mathcal {F}}^{-1}{2u(f)}}.}

Negative frequency components[edit]

Since s(t)=operatorname {Re}[s_{{mathrm  {a}}}(t)], restoring the negative frequency components is a simple matter of discarding operatorname {Im}[s_{{mathrm  {a}}}(t)] which may seem counter-intuitive. We can also note that the complex conjugate s_{{mathrm  {a}}}^{*}(t) comprises only the negative frequency components. And therefore {displaystyle s(t)=operatorname {Re} [s_{mathrm {a} }^{*}(t)]} restores the suppressed positive frequency components. Another viewpoint is that the imaginary component in either case is a term that subtracts frequency components from s(t). The operatorname {Re} operator removes the subtraction, giving the appearance of adding new components.

Examples[edit]

Example 1[edit]

s(t)=cos(omega t),   where  omega >0.

Then:

{displaystyle {begin{aligned}{hat {s}}(t)&=cos left(omega t-{frac {pi }{2}}right)=sin(omega t),\s_{mathrm {a} }(t)&=s(t)+j{hat {s}}(t)=cos(omega t)+jsin(omega t)=e^{jomega t}.end{aligned}}}

The last equality is Euler’s formula, of which a corollary is {textstyle cos(omega t)={frac {1}{2}}left(e^{jomega t}+e^{j(-omega )t}right).} In general, the analytic representation of a simple sinusoid is obtained by expressing it in terms of complex-exponentials, discarding the negative frequency component, and doubling the positive frequency component. And the analytic representation of a sum of sinusoids is the sum of the analytic representations of the individual sinusoids.

Example 2[edit]

Here we use Euler’s formula to identify and discard the negative frequency.

{displaystyle s(t)=cos(omega t+theta )={frac {1}{2}}left(e^{j(omega t+theta )}+e^{-j(omega t+theta )}right)}

Then:

s_{{mathrm  {a}}}(t)={begin{cases}e^{{j(omega t+theta )}}  = e^{{j|omega |t}}cdot e^{{jtheta }},&{text{if}} omega >0,\e^{{-j(omega t+theta )}}= e^{{j|omega |t}}cdot e^{{-jtheta }},&{text{if}} omega <0.end{cases}}

Example 3[edit]

This is another example of using the Hilbert transform method to remove negative frequency components. We note that nothing prevents us from computing s_{{mathrm  {a}}}(t) for a complex-valued s(t). But it might not be a reversible representation, because the original spectrum is not symmetrical in general. So except for this example, the general discussion assumes real-valued s(t).

s(t)=e^{{-jomega t}}, where omega >0.

Then:

{displaystyle {begin{aligned}{hat {s}}(t)&=je^{-jomega t},\s_{mathrm {a} }(t)&=e^{-jomega t}+j^{2}e^{-jomega t}=e^{-jomega t}-e^{-jomega t}=0.end{aligned}}}

Properties[edit]

Instantaneous amplitude and phase[edit]

A function in blue and the magnitude of its analytic representation in red, showing the envelope effect.

An analytic signal can also be expressed in polar coordinates:

s_{{mathrm  {a}}}(t)=s_{{mathrm  {m}}}(t)e^{{jphi (t)}},

where the following time-variant quantities are introduced:

In the accompanying diagram, the blue curve depicts s(t) and the red curve depicts the corresponding s_{{mathrm  {m}}}(t).

The time derivative of the unwrapped instantaneous phase has units of radians/second, and is called the instantaneous angular frequency:

{displaystyle omega (t)triangleq {frac {dphi }{dt}}(t).}

The instantaneous frequency (in hertz) is therefore:

{displaystyle f(t)triangleq {frac {1}{2pi }}omega (t).}  [3]

The instantaneous amplitude, and the instantaneous phase and frequency are in some applications used to measure and detect local features of the signal. Another application of the analytic representation of a signal relates to demodulation of modulated signals. The polar coordinates conveniently separate the effects of amplitude modulation and phase (or frequency) modulation, and effectively demodulates certain kinds of signals.

Complex envelope/baseband[edit]

Analytic signals are often shifted in frequency (down-converted) toward 0 Hz, possibly creating [non-symmetrical] negative frequency components:

{displaystyle {s_{mathrm {a} }}_{downarrow }(t)triangleq s_{mathrm {a} }(t)e^{-jomega _{0}t}=s_{mathrm {m} }(t)e^{j(phi (t)-omega _{0}t)},}

where omega _{0} is an arbitrary reference angular frequency.[2]

This function goes by various names, such as complex envelope and complex baseband. The complex envelope is not unique; it is determined by the choice of omega _{0}. This concept is often used when dealing with passband signals. If s(t) is a modulated signal, omega _{0} might be equated to its carrier frequency.

In other cases, omega _{0} is selected to be somewhere in the middle of the desired passband. Then a simple low-pass filter with real coefficients can excise the portion of interest. Another motive is to reduce the highest frequency, which reduces the minimum rate for alias-free sampling. A frequency shift does not undermine the mathematical tractability of the complex signal representation. So in that sense, the down-converted signal is still analytic. However, restoring the real-valued representation is no longer a simple matter of just extracting the real component. Up-conversion may be required, and if the signal has been sampled (discrete-time), interpolation (upsampling) might also be necessary to avoid aliasing.

If omega _{0} is chosen larger than the highest frequency of s_{{mathrm  {a}}}(t), then {displaystyle {s_{mathrm {a} }}_{downarrow }(t)} has no positive frequencies. In that case, extracting the real component restores them, but in reverse order; the low-frequency components are now high ones and vice versa. This can be used to demodulate a type of single-sideband signal called lower sideband or inverted sideband.

Other choices of reference frequency are sometimes considered:

In the field of time-frequency signal processing, it was shown that the analytic signal was needed in the definition of the Wigner–Ville distribution so that the method can have the desirable properties needed for practical applications.[5]

Sometimes the phrase “complex envelope” is given the simpler meaning of the complex amplitude of a (constant-frequency) phasor;[a][b]
other times the complex envelope {displaystyle s_{m}(t)} as defined above is interpreted as a time-dependent generalization of the complex amplitude.[c] Their relationship is not unlike that in the real-valued case: varying envelope generalizing constant amplitude.

Extensions of the analytic signal to signals of multiple variables[edit]

The concept of analytic signal is well-defined for signals of a single variable which typically is time. For signals of two or more variables, an analytic signal can be defined in different ways, and two approaches are presented below.

Multi-dimensional analytic signal based on an ad hoc direction[edit]

A straightforward generalization of the analytic signal can be done for a multi-dimensional signal once it is established what is meant by negative frequencies for this case. This can be done by introducing a unit vector {boldsymbol  {hat  {u}}} in the Fourier domain and label any frequency vector {boldsymbol  xi } as negative if {boldsymbol  xi }cdot {boldsymbol  {hat  {u}}}<0. The analytic signal is then produced by removing all negative frequencies and multiply the result by 2, in accordance to the procedure described for the case of one-variable signals. However, there is no particular direction for {boldsymbol  {hat  {u}}} which must be chosen unless there are some additional constraints. Therefore, the choice of {boldsymbol  {hat  {u}}} is ad hoc, or application specific.

The monogenic signal[edit]

The real and imaginary parts of the analytic signal correspond to the two elements of the vector-valued monogenic signal, as it is defined for one-variable signals. However, the monogenic signal can be extended to arbitrary number of variables in a straightforward manner, producing an (n + 1)-dimensional vector-valued function for the case of n-variable signals.

See also[edit]

  • Practical considerations for computing Hilbert transforms
  • Negative frequency

Applications[edit]

  • Single-sideband modulation
  • Quadrature filter
  • Causal filter

Notes[edit]

  1. ^ “the complex envelope (or complex amplitude)”[6]
  2. ^ “the complex envelope (or complex amplitude)”, p. 586 [7]
  3. ^ “Complex envelope is an extended interpretation of complex amplitude as a function of time.” p. 85[8]

References[edit]

  1. ^ Smith, J.O. “Analytic Signals and Hilbert Transform Filters”, in Mathematics of the Discrete Fourier Transform (DFT) with Audio Applications, Second Edition, https://ccrma.stanford.edu/~jos/r320/Analytic_Signals_Hilbert_Transform.html, or https://www.dsprelated.com/freebooks/mdft/Analytic_Signals_Hilbert_Transform.html, online book, 2007 edition, accessed 2021-04-29.
  2. ^ a b Bracewell, Ron. The Fourier Transform and Its Applications. McGraw-Hill, 2000. pp. 361-362
  3. ^ B. Boashash, “Estimating and Interpreting the Instantaneous Frequency of a Signal-Part I: Fundamentals”, Proceedings of the IEEE, Vol. 80, No. 4, pp. 519–538, April 1992
  4. ^ Justice, J. (1979-12-01). “Analytic signal processing in music computation”. IEEE Transactions on Acoustics, Speech, and Signal Processing. 27 (6): 670–684. doi:10.1109/TASSP.1979.1163321. ISSN 0096-3518.
  5. ^ B. Boashash, “Notes on the use of the Wigner distribution for time frequency signal analysis”, IEEE Trans. on Acoustics, Speech, and Signal Processing , vol. 26, no. 9, 1987
  6. ^ Hlawatsch, Franz; Auger, François (2013-03-01). Time-Frequency Analysis. John Wiley & Sons. ISBN 9781118623831.
  7. ^ Driggers, Ronald G. (2003-01-01). Encyclopedia of Optical Engineering: Abe-Las, pages 1-1024. CRC Press. ISBN 9780824742508.
  8. ^ Okamoto, Kenʼichi (2001-01-01). Global Environment Remote Sensing. IOS Press. ISBN 9781586031015.

Further reading[edit]

  • Leon Cohen, Time-frequency analysis, Prentice Hall, Upper Saddle River, 1995.
  • Frederick W. King, Hilbert Transforms, vol. II, Cambridge University Press, Cambridge, 2009.
  • B. Boashash, Time-Frequency Signal Analysis and Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003.

External links[edit]

  • Analytic Signals and Hilbert Transform Filters

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