Спектральная плотность
В статистической радиотехнике и физике при изучении детерминированных сигналов и случайных процессов широко используется их спектральное представление в виде спектральной плотности, которая базируется на преобразовании Фурье.
Если процесс имеет конечную энергию и квадратично интегрируем (а это нестационарный процесс), то для одной реализации процесса можно определить преобразование Фурье как случайную комплексную функцию частоты:
((1)) |
Однако она оказывается почти бесполезной для описания ансамбля. Выходом из этой ситуации является отбрасывание некоторых параметров спектра, а именно спектра фаз, и построении функции, характеризующей распределение энергии процесса по оси частот. Тогда согласно теореме Парсеваля энергия
((2)) |
Функция 
Перейдем теперь к стационарному в широком смысле центрированному случайному процессу , реализации которого с вероятностью 1 имеют бесконечную энергию и, следовательно, не имеют преобразования Фурье. Спектральная плотность такого процесса может быть найдена на основании теоремы Винера-Хинчина как преобразование Фурье от корреляционной функции:
((3)) |
Если существует прямое преобразование, то существует и обратное преобразование Фурье, которое по известной
определяет
:
((4)) |
Если полагать в формулах (3) и (4) соответственно
и
, имеем
((5)) |
((6)) |
Формула (6) с учетом (2) показывает, что дисперсия определяет полную энергию стационарного случайного процесса, которая равна площади под кривой спектральной плотности. Размерную величину
можно трактовать как долю энергии, сосредоточенную в малом интервале частот от
до
. Если понимать под
случайный (флуктуационный) ток или напряжение, то величина
будет иметь размерность энергии [В 2 /Гц] = [В 2 с]. Поэтому
иногда называют энергетическим спектром. В литературе часто можно встретить другую интерпретацию:
– рассматривается как средняя мощность, выделяемая током или напряжением на сопротивлении 1 Ом. При этом величину
называют спектром мощности случайного процесса.
Свойства спектральной плотности
- Энергетический спектр стационарного процесса (вещественного или комплексного) – неотрицательная величина:
- Энергетический спектр вещественного стационарного в широком смысле случайного процесса есть действительная и четная функция частоты:
- Корреляционная функция
и энергетический спектр
стационарного в широком смысле случайного процесса обладают всеми свойствами, характерными для пары взаимных преобразований Фурье. В частности, чем «шире» спектр
тем «уже» корреляционная функция
, и наоборот. Этот результат количественно выражается в виде принципа или соотношения неопределенности.
Спектральная плотность сигнала
Пусть интервал разложения сигнала (см. рис. 2.1) стремится к бесконечности. При его увеличении частота = 2п/Т уменьшается до бесконечно малой величины:

Расстояние между спектральными компонентами при этом уменьшается до бесконечно малой величины, а значения превращаются в текущие значения частоты со (см. рис. 2.2). Интервал разложения стремится к бесконечной величине. Это позволяет при вычислении предела ряда Фурье в комплексной форме заменить знак суммы знаком интеграла, основную частоту О)! = 2п/Т — на ?/со, а кратную частоту к(о< заменить текущей частотой со:

Интеграл, который записан в скобках выражения (2.13), обозначим

Тогда выражение (2.13) запишется более компактно:

Выражения (2.14) и (2.15) называются соответственно прямым и обратным преобразованиями Фурье. Функция 5(/со) называется
спектральной плотностью. Она является комплексной и имеет размерность [В/Гц], если размерность сигнала и <Р)[В].
Преобразование Фурье (2.14) может быть вычислено на основе общих правил интегрирования, если сигнал удовлетворяет условию абсолютной интегрируемости:

Это условие означает, что преобразование (2.14) существует для тех сигналов, площадь под кривой |м(?)| которых ограничена.
К этому классу не относятся, например, периодические сигналы, которые не удовлетворяют условию абсолютной интегрируемости. Однако это не означает, что для периодических сигналов спектральная плотность не может быть вычислена. Методы вычислений, специально разработанные для этих целей, используют так называемые обобщенные функции. Примером обобщенной функции является дельта-функция. Некоторые свойства дельта-функции приведены в приложении 1.
Преобразуем спектральную плотность сигналов, которые удовлетворяют условию абсолютной интегрируемости. Такие сигналы ограничены во времени.
С учетом формулы Эйлера перепишем выражение (2.14):
где 
Модуль |5(/со)| называется спектральной плотностью амплитуд сигнала или амплитудно-частотной характеристикой
(АЧХ) спектральной плотности сигнала. Функция ср(со) определяет фазо-частотную характеристику (ФЧХ) спектральной плотности сигнала. АЧХ и ФЧХ спектральной плотности являются непрерывными функциями частоты.
Перейдем к анализу спектральной плотности сигналов, не удовлетворяющих условию абсолютной интегрируемости. Такие сигналы не ограничены во времени и имеют бесконечно большую энергию.
На основе сигнала Ц)(?), удовлетворяющего условию абсолютной интегрируемости, построим периодически повторяющийся сигнал

и вычислим его спектральную плотность:
где 
Размерность спектральной плотности периодически повторяющегося сигнала определяется размерностью спектральной плотности непериодического сигнала, из которого формируется периодически повторяющийся сигнал, т.е. [В/Гц].
Первый сомножитель полученного выражения в равенстве (2.16) определяет спектральную плотность ограниченного во времени сигнала и0(?), второй — спектральную плотность периодически повторяющейся дельта-функции 
Убедимся в этом, вычислив указанную плотность:

При вычислении интеграла использовано фильтрующее свойство дельта-функции (см. приложение 1).
Если периодически повторяющуюся дельта-функцию разложить в ряд Фурье в комплексной форме, то се спектральную плотность можно выразить иначе:


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

Эта функция равна нулю, если со Ф к(х)ь и равна если со = к(о<. Подставим в (2.16) новое выражение 5ф(/со):

Спектральная плотность периодически повторяющегося сигнала определяется значениями спектральной плотности ограниченного во времени сигнала г/0(?), отсчитанными через интервал, равный со^ = 2л/Т.
Вычислим значение спектральной плотности ограниченного отрезком времени Т сигнала:

Умножим левую и правую части равенства на коэффициент 2/Т:

где а(/&а>1) — спектр ограниченного во времени сигнала в базисе экспоненциальных функций.
С учетом последней формулы спектральную плотность периодически повторяющегося сигнала запишем в виде

где модуль спектра определяется в базисе экспоненциальных функций формулой (2.9), а спектр фаз — формулой (2.10).
Значения АЧХ и ФЧХ спектральной плотности ограниченного во времени сигнала г/о(0> отсчитанные через интервал (щ = 2п/Т в точках частотной оси кщ, к = 0, ±1, ±2. определяют АЧХ и ФЧХ спектральной плотности этого периодического сигнала.
Рассмотрим некоторые свойства спектральной плотности сигнала, удовлетворяющие условию абсолютной интегрируемости.
- 1. Спектральная плотность (2.14) — это комплексная и непрерывная функция частоты со, определенная в бесконечном интервале частот.
- 2. АЧХ и ФЧХ спектральной плотности удовлетворяют уравнениям

где +(л)? — выбранные значения частот.
3. Преобразования Фурье (2.14), (2.15) являются линейными преобразованиями. Поэтому спектральная плотность суммы сигналов равна сумме спектральных плотностей этих сигналов, а сумма сигналов определяется обратным преобразованием Фурье от суммы их спектральных плотностей:

где Uj(t) — i-й сигнал; б’/О’оз) — спектральная плотность г-го сигнала.
4. Спектральная плотность сигнала, ограниченная бесконечно малыми интервалами 2лА/(рис. 2.3) вблизи, например, частот -со0, со(), определяет гармонический сигнал с бесконечно малой амплитудой.
Убедимся в этом, считая, что из-за малости А/ значения спектральной плотности около частот -ю(), (н() равны соответственно S(—jco0) = |А(70)0)| _ — / “ о) и5(усо0) = |5(y’oe0)|e ;(p( — tO(,) . Найдем сигнал по спектральной плотности, приведенной на рис. 2.3:

Рис. 2.3. Спектральная плотность гармонического сигнала с бесконечно малой амплитудой
Поскольку в бесконечно малых интервалах спектральная плотность остается постоянной, можно вынести за знак интегралов выражения |50’со0)|е ;ф(10о) и |50’м0)|е — — ,ф(а)о) :

Как следует из полученной формулы, амплитуда полученного сигнала определяется значением спектральной плотности, функцией (бшл — )/^ и весьма малым диапазоном частот А/. При стремлении Д/ к нулю функция (81пх)/х стремится к единице, а амплитуда становится равной нулю.
5. Если все составляющие спектральной плотности ограниченного во времени сигнала сдвигаются по фазе на +(л)?о> то этот сиг- нал сдвигается во времени на величину +?0. Действительно:

6. При передаче ограниченного во времени сигнала через линейный четырехполюсник, АЧХ которого в полосе пропускания равна постоянной величине К0, а фазовая характеристика ср(со) = = -а)?0> форма этого сигнала остается неизменной, а сигнал запаздывает во времени на величину ?0: 
Пример 2.2. Рассчитать спектральную плотность задержанного во времени одиночного прямоугольного импульса с длительностью ти и амплитудой А, ее АЧХ и ФЧХ.
Решение. Спектральная плотность задержанного на время ?0 импульса равна

где м(?) — импульс, который расположен в начале координат;

Вычисления дают следующий результат:

Запишем эту плотность в виде
где

Последнее выражение определяет спектральную плотность сигнала и(?). В диапазоне частот [2я/ти] спектральная плотность является положительной величиной, д(со) = = 1. Поэтому в этом диапазоне фазовая характеристика ф(со) = 0, так как ,ф(со) = -1, и т.д.
АЧХ спектральной плотности задержанного импульса совпадает с АЧХ спектральной плотности сигнала «(?), а ФЧХ определяется выражением

Спектральная плотность прямоугольного импульса г/(?), АЧХ и ФЧХ этой плотности изображены на рис. 2.4.

Рис. 2.4. Спектральная плотность, АЧХ и ФЧХ спектральной плотности прямоугольного импульса
Пример 2.3. Вычислить спектральную плотность кодированного сигнала

где ак — элементы кодового слова, равные -1 или 1, т.е. = +1, и0(0 — прямоугольный импульс с амплитудой А и длительностью ти.
Решение. Применим формулу (2.14):

После замены переменной
, получим 
Пример 2.4. Вычислить спектральную плотность периодического сигнала, записанного в виде ряда Фурье в тригонометрической форме [см. формулу (2.11)]. Записать выражения АЧХ и ФЧХ постоянной, синусной и косинусной составляющих этого ряда.
Решение. Функции, определяющие формулу (2.11), — периодические, за исключением постоянной составляющей. Эту составляющую аппроксимируем периодической косинусной функцией с частотой, которая стремится к нулю:

Вычислим спектральную плотность периодического сигнала u(t) = = a cos fit, записав его в виде

где 

Спектральная плотность сигнала щ(():

Значение первого слагаемого, стоящего в скобках выражения, равно 1, если со = -Q, и равно 0 для других дискретных значений частоты со = kfl, k = 0, 1, ±2, ±3, ±4, . Значение второго слагаемого равно 1, если со = Q, и равно 0 для других дискретных значений частоты to = kQ, k = 0, -1, ±2, ±3, ±4, . Учитывая это, найдем спектральную плотность, АЧХ и ФЧХ спектральной плотности периодического сигнала u(t) = a cos Q?:

Значения АЧХ спектральной плотности в точках частотной оси со = +?2 равны паТ/(2п) = аТ/2.
Значения ФЧХ спектральной плотности гармонического сигнала в точках частотной оси со = равны 0.
По формуле спектральной плотности косинусоидального сигнала можно найти спектральную плотность постоянной составляющей:

АЧХ спектральной плотности постоянной составляющей определяется значением

Вычисление спектральной плотности синусоидального сигнала аналогично вычислению спектральной плотности косинусоидального сигнала.
Запишем периодический сигнал u(t) = bsinQ? в виде

где 

Спектральная плотность сигнала и0(О :

По найденному выражению найдем спектральную плотность периодического сигнала u(t) = b sin Qt:

АЧХ спектральной плотности этого сигнала в точках частотной оси со = +П:

Значения ФЧХ спектральной плотности сигнала в точках частотной оси со = +П равны -я/2, п/2.
Полученные формулы для спектральных плотностей гармонических сигналов позволяют найти спектральную плотность суммы этих сигналов:

где
— модуль спектра, равный амплитуде гармонического
сигнала; ф(П) = -экЛ%(Ь/а) — значение фазы спектра, равное значению начальной фазы этого сигнала.
Ряд Фурье в тригонометрической форме (2.11) содержит бесконечно большое число сумм гармонических сигналов:

Спектральная плотность этой суммы находится по последнему выражению спектральной плотности заменой П = ко)^. Используя эту формулу и формулу спектральной плотности постоянной составляющей, получим выражение спектральной плотности сигнала, записанного в виде ряда Фурье в тригонометрической форме:

где
— модуль спектра; ф^о^) = — значение фазы спектра, равное значению начальной фазы гармонического сигнала.
Для периодической последовательности импульсов, приведенной на рис. 2.1,

Спектральная плотность 

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

Рис. 2.5. Спектральная плотность периодической последовательности
График содержит информацию о постоянной составляющей и гармонических сигналах, входящих в ряд Фурье в тригонометрической форме.
Пример 2.5. По спектральной плотности, вид которой приведен на рис. 2.6, вычислить выражение для сигнала «(?)•

Рис. 2.6. Спектральная плотность сигнала
Решение. Спектральная плотность сигнала ограничена значениями частоты, равными -сов, сов. Найдем сигнал:
2 / UMKD_TMFIR_Anufriev / Konspekt_lekcij / Lekciya_5
Для практических приложений является важным установление связи между преобразованием сигнала и соответствующим этому преобразованию изменением спектральных характеристик.
Спектральная плотность сигнала, смещенного во времени.
Предположим, что сигнал
произвольной формы, существующий на интервале от
до
, имеет спектральную плотность
. Найдем спектральную плотность этого же сигнала при условии его задержки на интервал
, например преобразователем, называемым линией задержки. Функция времени задержанного сигнала при сохранении его формы запишется в виде:
(5.1)
Спектральная плотность задержанного сигнала
очевидно имеет вид:
(5.2)
Вводя новую переменную интегрирования
, получим:
=
. (5.3)
Из этого соотношения видно, что задержка во времени сигнала
на интервал
приводит к изменению фазовой характеристики спектра
(спектра фаз) на величину
. Очевидно, что в общем случае при сдвиге сигнала во времени на величину
его фазовый спектр изменится на величину
. Спектр амплитуд этого сигнала (модуль спектральной плотности) от положения сигнала на временной оси не зависит.
Спектральная плотность сигнала, сжатого во времени.
Пусть сигнал
длительностью
подвергся сжатию во времени в соответствии с рисунком:

Новый сжатый сигнал
связан с исходным сигналом соотношением:
,
(5.4)
Длительность сжатого сигнала очевидно равна
. Определим спектральную плотность сжатого сигнала
:
(5.5)
Вводя новую переменную интегрирования
, получаем:
(5.6)
Интеграл в правой части выражения (5.6) есть не что иное, как спектральная плотность исходного сигнала
при частоте
, т.е.:
. (5.7)
Итак, при сжатии сигнала в
раз на временной оси имеем:
уменьшение модуля спектральной плотности в
раз;
расширение во столько же раз его спектральных составляющих на оси частот.
Очевидно, при расширении исходного сигнала во времени (т.е. при
) имеют место обратные процессы: сужение спектра и увеличение модуля спектральной плотности.
Можно также показать, что длительность сигнала и ширина его спектра амплитуд не могут быть одновременно ограничены конечными интервалами: если длительность сигнала ограничена, то спектр его неограничен, и, наоборот, сигнал с ограниченным спектром длится бесконечно долго. Говорят, что ширина спектра и длительность импульса связаны соотношением неопределенности:
, где
– длительность импульса, а
– ширина спектра (практическая ширина),
– постоянная, зависящая от формы импульса (в первом приближении принимают
).
Спектральная плотность на выходе сумматора сигналов.
Преобразование Фурье, определяющее спектральную плотность заданного сигнала, является линейным преобразованием. Если на вход сумматора подать некоторую совокупность сигналов

…, обладающих спектральными плотностями соответственно
,
,
, …, то взвешенной сумме сигналов на выходе сумматора
будет соответствовать спектральная плотность:
, (5.8)
где
– произвольные числовые коэффициенты. Для доказательства этого утверждения следует подставить сумму сигналов в формулу для прямого преобразования Фурье.
Спектральная плотность продифференцированного сигнала.
Пусть сигнал
и его спектральная плотность
заданы. Подадим сигнал на вход линейного устройства, осуществляющего дифференцирование сигнала. Сигнал на выходе дифференцирующего устройства будет иметь вид:
, (5.9)
где
– константа преобразования.
Используя свойство преобразования Фурье, часто записываемое в виде:
, (5.10)
где символом
обозначено прямое преобразование Фурье над функцией
,
(5.11)
Итак, дифференцирование сигнала по времени эквивалентно алгебраической операции умножения спектральной плотности на множитель
, часто называемым оператором дифференцирования в частотной области. Отметим, что при дифференцировании скорость изменения сигнала во времени возрастает, а сам сигнал укорачивается. Следовательно, модуль спектра производной имеет большие значения в области высоких частот по сравнению с модулем исходного сигнала.
Спектральная плотность сигнала на выходе интегратора. Во многих измерительных системах находят применение так называемые интеграторы – преобразователи, выходной сигнал которых пропорционален интегралу от входного воздействия
:
,
– константа преобразования. (5.12)
По аналогии с операцией дифференцирования нетрудно найти формулу связи спектральных плотностей сигналов на входе и выходе интегратора:
. (5.13)
Так как модуль знаменателя растет линейно с увеличением частоты, это свидетельствует о том, что рассматриваемый интегратор действует подобно фильтру нижних частот, ослабляя высокочастотные спектральные составляющие входного сигнала.
Практическая ширина спектра сигнала.
Реальные сигналы всегда ограничены во времени, следовательно, их амплитудный спектр теоретически неограничен. Однако реальные сигналы генерируются и передаются устройствами, содержащими инерционные элементы (например, емкости и индуктивности в электрических цепях и прочих преобразователях). Поэтому они не могут содержать гармонических составляющих сколь угодно высоких частот.
В связи с этим возникает необходимость ввести в рассмотрение модели сигналов, обладающих как конечной длительностью, так и ограниченным спектром. При этом в соответствии с каким-либо критерием дополнительно ограничивается либо ширина спектра, либо длительность сигнала, либо оба параметра одновременно.
Чаще всего в качестве такого критерия используется энергетический критерий, согласно которому практическую ширину амплитудного спектра
выбирают так, чтобы в нем была сосредоточена подавляющая часть энергии сигнала. Для этого используют равенство Парсеваля, позволяющее определить энергию сигнала
либо через функцию
, описывающую реальный сигнал длительностью
, либо через модуль ее спектральной плотности
:
. (5.14)
Практическая ширина спектра сигнала, сосредоточенная в диапазоне частот от 0 до некоторого значения
, определяется из соотношения:
. (5.15)
Здесь
– граничная частота, определяющая верхнее значение спектра сигнала;
– коэффициент, достаточно близкий к 1 (на практике его значение выбирают в интервале от 0.9 до 0,998 в зависимости от требований к качеству воспроизведения сигнала). Значение
означает, что в полосе частот от
до
содержится 99 % энергии сигнала. Значение граничной частоты находят, решая трансцендентное уравнение (5.15).
Практическая ширина спектра экспоненциального импульса.
Задача: определить граничную частоту спектра сигнала вида
, (5.16)
ориентируясь на практическую ширину спектра сигнала с 
Принять следующие значения:
мВ,
1/сек.
Сигналы такого типа, называемые экспоненциальными видеоимпульсами, генерируются многими датчиками (например, полупроводниковыми детекторами при регистрации параметров ионизирующих излучений) . Условие
обеспечивает достаточно быстрое (экспоненциальное) уменьшение значений сигнала с ростом времени.
Спектральная плотность такого сигнала равна:
(5.17)
Подставляя пределы, получаем:
. (5.18)
Спектр амплитуд такого импульса равен:
(5.19)
Трансцендентное уравнение, решение которого позволяет определить
, имеет вид:
. (5.20)
Так как
,
получаем:
. Отсюда
.
Практическую длительность подобных импульсов в измерительной технике обычно определяют из условия десятикратного уменьшения уровня сигнала:

. Отсюда
.
Спектральная плотность и ее свойства. Теоремы о спектрах
Спектральная плотность и сигнал связаны между собой парой преобразований Фурье:


Все свойства спектральной плотности объединены в основных теоремах о спектрах.
I. Свойство линейности.

Если имеется некоторая совокупность сигналов причём ,…, то взвешенная сумма сигналов преобразуется по Фурье следующим образом:


Здесь — произвольные числовые коэффициенты.
II. Теорема о сдвигах.

Предположим, что для сигнала известно соответствие . Рассмотрим такой же сигнал, но возникающий на секунд позднее. Принимая точку за новое начало отсчёта времени, обозначим этот смещённый сигнал как . Введём замену переменной: . Тогда ,

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

Предположим, что исходный сигнал подвергнут изменению масштаба времени. Это означает, что роль времени играет новая независимая переменная (- некоторое вещественное число.) Если > 1, то происходит “ сжатие” исходного сигнала; если же 0<<1, то сигнал “растягивается” во времени. Если , то :


Произведём замену переменной , тогда , откуда следует:

При сжатии сигнала в раз на временной оси во столько же раз расширяется его спектр на оси частот. Модуль спектральной плотности при этом уменьшается в раз.
Очевидно, что при растягивании сигнала во времени ( т.е. при <1) имеет место сужение спектра и увеличение модуля спектральной плотности.
IV. Теорема о спектре производной и неопределённого интеграла.



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

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

Представляя экспоненциальную функцию рядом Тейлора:

подставляя этот ряд в (2.6) и ограничиваясь первыми двумя членами ряда, находим


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

Вторая часть теоремы. Рассмотренная функция является неопределённым интегралом по отношению к функции. Интеграл это есть, значит — его спектральная плотность, а из формулы (2.7) равна:


Таким образом, множитель служит оператором интегрирования в частотной области.
V. Теорема о свёртке.
При суммировании сигналов их спектры складываются. Однако спектр произведения сигналов не равен произведению спектров, а выражается некоторым специальным интегральным соотношением между спектрами сомножителей.

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


Применив обратное преобразование Фурье, выразим сигнал через его спектральную плотность и подставим результат в (2.9):


Изменив порядок интегрирования, будем иметь:



Интеграл, стоящий в правой части называют свёрткой функций и . Символически операция свёртки обозначается как *


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

Операция свёртки коммутативна, т.е. допускает изменения порядка следования преобразуемых функций:

Теорема о свёртке может быть обращена: если спектральная плотность некоторого сигнала представляется в виде произведения , причём


и , то сигнал является свёрткой сигналов и , но уже не в частотной, а во временной области:

Пусть два сигнала и , в общем случае комплексные, определены своими обратными преобразованиями Фурье:


Найдём скалярное произведение этих сигналов, выразив один из них, например , через его спектральную плотность:


Здесь внутренний интеграл представляет собой спектральную плотность сигнала поэтому:

Скалярное произведение двух сигналов с точностью до коэффициента пропорционально скалярному произведению их спектральных плотностей.
((1))
((2))
((3))
((4))
((5))
((6))