Свойства спектральной плотности мощности. Спектральная плотность мощности случайного процесса. Методы вычисления СПМ. Функции спектральной плотности Максимальная спектральная плотность мощности

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

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

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

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

Разделив эту энергию на получим среднюю мощность k-й реализации на отрезке Т

При увеличении Т энергия возрастает, однако отношение стремится к некоторому пределу. Совершив предельный переход получим

представляет собой спектральную плотность средней мощности рассматриваемой реализации.

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

Опуская индекс k, получаем окончательное выражение для средней мощности случайного процесса

Если рассматривается случайный процесс с ненулевым средним значением то спектральную плотность следует представить в форме

Под энергией сигнала иЦ) понимают величину

Если сигнал имеет конечную длительность Т, т.е. не равен нулю на отрезке времени [-Т/ 2, Т/ 2], то его энергия

Запишем выражение для энергии сигнала, используя формулу (2.15):

где

Полученное равенство называется равенством Парсеваля. Оно определяет энергию сигнала через временную функцию или спектральную плотность энергии, которая равна |5(/0))| 2 . Спектральная плотность энергии называется также энергетическим спектром.

Рассмотрим сигнал, существующий на ограниченном интервале времени. К такому сигналу применимо равенство Парсеваля. Следовательно,

Разделим левую и правую части равенства на интервал времени, равный Г, и устремим этот интервал к бесконечности:

С увеличением Т энергия незатухающих сигналов возрастает,

однако отношение может стремиться к определенному пределу. Этот предел называется спектральной плотностью мощности С(со). Размерность спектральной плотности мощности: [В 2 Дц].

Автокорреляционная функция

Автокорреляционная функция сигнала и (?) определяется следующим интегральным выражением:

где т - аргумент, определяющий функцию Я(х) и имеющий размерность времени; и(? + т) - исходный сигнал, сдвинутый во времени на величину -т.

Автокорреляционная функция имеет следующие свойства.

1. Значение автокорреляционной функции при сдвиге т = О равно энергии сигнала Е:

2. Автокорреляционная функция при сдвигах т Ф 0 меньше энергии сигнала:

3. Автокорреляционная функция является четной функцией, т.е.

В справедливости свойств 2 и 3 убедимся на примере.

Пример 2.6. Вычислить автокорреляционные функции сигналов: видеосигнала, представленного на рис. 2.7, я, и радиосигнала с теми же амплитудой и длительностью. Несущая частота радиосигнала равна щ, а начальная фаза равна 0.

Решение. Первую задачу решим графическим способом. Автокорреляционная функция определяется интегралом от произведения функции и (?) и ее смещенной во времени копии. Смещение видеосигнала найдем из уравнения? + т = 0. График функции м(? + т) приведен на рис. 2.7, б. Площадь, определяемая графиком произведения м(?)м(? + т) (рис. 2.7, в), равна

Функция Д(т) определяется уравнением прямой (рис. 2.7, г). Функция имеет максимум, если значение аргумента т = 0, и равна 0, если т = т и. Для других значений аргумента /?(т)

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

Рис. 2.7.

видеоимпульса:

а - прямоугольный видеоимпульс; б - задержанный во времени прямоугольный импульс; в - произведение импульсов; г - автокорреляционная функция

Окончательное выражение для автокорреляционной функции

Функция приведена на рис. 2.7, г и имеет треугольный вид.

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

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

Выражение для автокорреляционной функции радиоимпульса состоит из двух слагаемых. Первое из них определяется произведением треугольной функции и гармонического сигнала. На выходе согласованного фильтра это слагаемое реализуется в виде ромбовидного радиоимпульса. Второе слагаемое определяется произведением треугольной функции и функций (втд^/лг, расположенных в точках т = +т и. Значения функций (втх)/:*:, которые оказывают заметное влияние на второе слагаемое автокорреляционной функции, весьма быстро убывают при изменении аргумента т от -т и до оо и от т и до -°о. Решив уравнение

можно найти интервалы задержки, в пределах которых значения функций (втлс)/;*; еще влияют на поведение функции /?(т). Для положительных значений задержки

где 7о - период гармонического сигнала.

Аналогично находится интервал для отрицательных значений задержки.

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

Выявим связь автокорреляционной функции #(т) со спектральной плотностью энергии сигнала |5(/со)| 2 . Для этого выразим сдвинутый во времени сигнал и(1ь + т) через его спектральную плотность 5(/со):

Подставим данное выражение в выражение (2.21). В результате получим

Нетрудно убедиться также в справедливости равенства

Разделим обе части равенства (2.23) на интервал времени Т и устремим величину Т к бесконечности:

С учетом формулы (2.20) перепишем полученное выражение:

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

Обобщением понятия «автокорреляционная функция» является взаимно корреляционная функция, которая представляет собой скалярное произведение двух сигналов:

Рассмотрим основные свойства взаимно корреляционной функции.

1. Перестановка сомножителей под знаком интеграла изменяет знак аргумента взаимно корреляционной функции:

В приведенных преобразованиях использована замена t + т = х.

  • 2. Взаимно корреляционная функция, в отличие от автокорреляционной функции, не является четной относительно аргумента т.
  • 3. Взаимно корреляционная функция определяется обратным преобразованием Фурье от произведения спектральных плотностей сигналов u(t), v(t) :

Эта формула может быть выведена аналогично формуле (2.22).

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

сигналом v(t ) = Uq(?)

где R(t) - автокорреляционная функция непериодического сигнала u 0 (t).

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

Взаимно корреляционная функция Я ил> (т) между сигналом u(t ) и сигналом

равна - длительность сигнала v(t).

Действительно, вследствие того что период сигнала u(t ) равен Т и

взаимно корреляционная функция где

Вычисляя предел функции (2п + 1)7? м Мо (т) при п -> определим выражение для автокорреляционной функции периодического сигнала:

Размерность функции: [В 2 /Гц].

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

Разделим последнее выражение на интервал, равный (2п + 1 )Т. В результате получим функцию


так как вследствие периодичности функции - т + Т) = - т).

Полученная формула определяет функцию В(т) как предел отношения автокорреляционной функции сигнала, существующего в интервале времени (2п + 1 )Т, к этому интервалу и стремлении его к бесконечности. Этот предел для периодически повторяющегося сигнала называется автокорреляционной функцией периодического сигнала. Размерность этой функции: [В 2 ].

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

где 0)1 = 2п/Т.

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

Пример 2.7. Вычислить периодическую автокорреляционную функцию сигнала и(ф) = А бш СИ. По найденной функции, ограниченной одним периодом, определить спектральную плотность мощности.

Решение. Подставляя в выражение (2.26) заданный сигнал, получим выражение для периодической автокорреляционной функции:

Полученное выражение подставим в формулу (2.24) и найдем спектральную плотность мощности:

Пример 2.8. Для периодической нормированной автокорреляционной функции шумоподобного сигнала (М-последовательности с периодом N = 1023) вычислить спектральную плотность мощности. (Периодическая функция для последовательности меньшей длины (IV= 15) приведена на рис. 3.39.)

Решение. Для сравнительно большого периода ЛГ = 1023 значения автокорреляционной функции в интервале Т - То > т > То, где То - длительность импульса шумоподобной последовательности, примем равными нулю. В этом случае автокорреляционная функция определяется периодически повторяющейся с периодом Т последовательностью треугольных импульсов. Основание каждого треугольника равно 2то, а его высота равна 1. Уравнение, определяющее автокорреляционную функцию в пределах одного периода, равно В(т) = 1 - |т|/хо- Учитывая четность этой функции, определим коэффициенты ряда Фурье:

При вычислении интеграла использована формула

Подставляя вычисленные коэффициенты в формулу (2.27), ползшим

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

Корреляционные функции цифровых сигналов связаны с корреляционными функциями последовательностей символов. Для кодовой последовательности (см. § 1.3) конечного числа N

двоичных символов автокорреляционная функция записывается в виде

где - двоичные символы, равные 0 или 1, или символы, равные -1, 1; д = О, 1, 2, ..., N - .

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

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

Сигналы, построенные с использованием кодов или их кодовых последовательностей, называются кодированными сигналами.

Большинство свойств автокорреляционной функции кодовой последовательности совпадает с рассмотренными выше свойствами автокорреляционной функции сигнала.

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

Если символы равны -1, 1, то г(0) = N.

Значения автокорреляционной функции при других сдвигах меньше г(0).

Автокорреляционная функция кодовой последовательности является четной функцией.

Обобщением автокорреляционной функции является взаимно корреляционная функция. Для кодовых последовательностей одинаковой длины эта функция

где 2 } 0 6/, - символы соответственно первой и второй последовательности.

Многие свойства функции г 12 (д) совпадают со свойствами взаимно корреляционной функции рассмотренных выше сигналов. Если функция г^(д), I Ф для любой пары кода при сдвиге д = О равна нулю, то такие коды называются ортогональными. Краткое описание некоторых используемых в системах связи кодов приведено в приложениях 2-4.

Взаимно корреляционная функция между кодовой последовательностью и периодически повторяющейся той же последовательностью называется периодической автокорреляционной функцией кодовой последовательности. Выражение для функции следует из выражений (2.25), (2.26):

где г(д) - непериодическая автокорреляционная функция кодовой последовательности; д - значение сдвига между последовательностями.

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

где а/г, а^+ц - элементы кодовой последовательности.

Периодическая автокорреляционная функция кодовой последовательности равна взаимно корреляционной функции, вычисленной для кодовой последовательности и циклически сдвинутых символов этой последовательности. Циклически сдвинутые кодовые последовательности, полученные по исходной последовательности а 0 = а 0 ,а { ,а 2 , ..., а м _ ь приведены ниже. Кодовая последовательность а { получена в результате сдвига исходной последовательности а 0 па один символ вправо и переноса последнего символа а дм в начало сдвинутой последовательности. Остальные последовательности получены аналогично:

Пример 2.9. Вычислить автокорреляционную и периодическую автокорреляционную функцию кодированного сигнала (рис. 2.8, а)

где и 0 (О - прямоугольный импульс с амплитудой А и длительностью т и.

Этот сигнал построен из прямоугольных импульсов, знак которых определяется весовыми коэффициентами: а 0 = ,а. = 1, а 2 = -1, а их число N = 3. Длительность сигнала равна Зт и.

Решение. Подставляя выражение для сигнала в формулу (2.21), получим

Произведем замену переменной t - кт н на х:

Обозначим: & - т = - и заменим дискретные переменные &, т на переменные к, ц. В результате получим

График автокорреляционной функции для заданного сигнала показан на рис. 2.8, б. Эта функция зависит от автокорреляционной функции /? 0 (т) прямоугольного импульса и значений автокорреляционной функции г(

Рис. 2.8. Автокорреляционная функция кодированного сигнала: а - кодированный сигнал; 6 - автокорреляционная функция сигнала; в - автокорреляционная функция периодического сигнала

Вычислим периодическую автокорреляционную функцию, используя рассчитанную выше автокорреляционную функцию, полученные значения автокорреляционной функции кодовой последовательности и формулу (2.28).

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

Подставим заданное значение N = 3 в полученную формулу:

С учетом значений автокорреляционной функции кодовой последовательности К+З) = 0, г(+ 2) = -1, г(+1) = О, КО) = 3 запишем окончательное выражение для одного периода периодической автокорреляционной функции сигнала:

График функции приведен на рис. 2.8, в.

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

Для незатухающих во времени сигналов энергия бесконечно велика и интеграл (1.54) расходится. Задание спектра амплитуд невозможно. Однако средняя мощность Рср, определяемая соотношением

оказывается конечной. Поэтому применяется более широкое понятие "спектральная плотность мощности". Определим ее как производную средней мощности сигнала по частоте и обозначим Сk(щ):

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

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

Для установления связи между спектральной плотностью Сk(щ) и спектром амплитуд воспользуемся сигналом u(t), существующим на ограниченном интервале времени (-T<. t

где - спектральная плотность мощности сигнала, ограниченного во времени.

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

Функция автокорреляции детерминированного сигнала

Теперь в частотной области имеется две характеристики: спектральная характеристика и спектральная плотность мощности. Спектральной характеристике, содержащей полную информацию о сигнале u(t), соответствует преобразование Фурье в виде временной функции. Выясним, чему соответствует во временной области спектральная плотность мощности, лишенная фазовой информации.

Следует предположить, что одной и той же спектральной плотности мощности соответствует множество временных функций, различающихся фазами. Советским ученым Л.Я. Хинчиным и американским ученым Н. Винером практически одновременно было найдено обратное преобразование Фурье от спектральной плотности мощности:


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

Справедливо и второе интегральное соотношение для пары преобразования Фурье:

Пример 1.6 Определить временную функцию· автокорреляции гармонического сигнала u(t) = u0 cos(t-ц). В соответствии с (1.64)

Проведя несложные преобразования


окончательно имеем

Как и следовало ожидать, ru() не зависит от ц и, следовательно, (1.66) справедливо для целого множества гармоник, различающихся фазами.

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

Спектральные плотности ряда сигналов приведены в табл. 4.2.

1) Импульс прямоугольной формы (табл. 4.2, поз. 4). Колебание, изобра­женное на рис. (4.28, а), можно записать в виде

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

График спектральной плотности (рис. 4.28, а) построен на основе прове­данного ранее анализа спектра периодической последовательности однополярных, прямоугольных импульсов (4.14). Как видно из (рис. 4.28, б), функция обра­щается в нуль при значениях аргумента = n , где п - 1, 2, 3, ... - лю­бое целое число. При этом угловые частоты равны = .

Рис. 4.28. Импульс прямоугольной формы (а) и его спектральная плотность (б)

Спектральная плотность импульса при численно равна его площади, т.еG (0)=A . Это положение справедливо для импульса s (t ) произвольной формы. Действительно, полагая в общем выражении (4.41) = 0, получим

т. е. площадь импульса s (t ).

Таблица 4.3.

Сигнал s (t )

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

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

Рис. 4.29. Графики амплитудного (а) Рис. 4.30. Импульс прямоугольной формы, и фазового (б) спектров сдвинутый на время

При сдвиге импульса вправо (за­паздывание) на время (рис. 4.30) фазовый спектр изменяется на величи­ну, определяемую аргументом множителяexp() (табл. 4.2, поз. 9). Результирующий фазовый спектр запаздывающего импульса изо­бражен на рис. 4.29, б пунктирной ли­нией.

2) Дельта-функция (табл. 4.3, поз. 9). Спектральную плотность – функции находим по формуле (4.41), используя фильтрующее свойствоδ -функции:

Таким образом, амплитудный спектр равномерный и определяется пло­щадьюδ -функции [= 1], а фазовый спектр равен нулю [= 0].

Обратным преобразованием Фурье от функции = 1 пользуются как одним из определенийδ -функции:

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

Амплитудный и фазовый спектры функции показаны в табл. 4.3, поз. 10. Обратное преобразование Фурье от функции имеет вид

3) Гармоническое колебание (табл. 4.3, поз. 12). Гармони­ческое колебание не является абсолютно интегрируемым сигналом. Тем не ме­нее для определения его спектральной плотностиприменяют прямое пре­образование Фурье, записывая формулу (4.41) в виде:

Тогда с учетом (4.47) получаем

δ(ω) – дельта-функции, смещенные по оси частот на частоту , соответственно вправо и влево относительно. Как видно из (4.48), спектральная плотность гармонического колебания с конечной амплитудой принимает бесконечно боль­шое значение на дискретных частотахи.

Выполняя аналогичные преобразования, можно получить спектральную плотность колебания (табл. 4.3, поз. 13)

4) Функция вида (табл. 4.3, поз. 11)

Спектральная плотность сигнала в виде постоянного уровня А определяется по формуле (4.48), положив = 0:

5) Единичная функция (или единичный скачок) (табл. 4.3, поз. 8). Функция не является абсолютно интегрируемой. Если представить как предел экспоненциального импульса , т. е.

то спектральную плотность функцииможно определить как предел спектральной плотности экспоненциального импульса (табл. 4.3, поз. 1) при :

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

Поэтому пределом первого слагаемого можно считать функцию . Преде­лом второго слагаемого является функция. Окончательно получим

Наличие двух слагаемых в выражении (4.51) согласуется с представлением функции в виде 1/2+1/2sign(t ). Постоянной составляющей 1/2 со­гласно (4.50) соответствует спектральная плотность , а нечетной функции - мнимое значение спектральной плотности .

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

6) Комплексный экспоненциальный сигнал (табл. 4.3, поз. 16). Если представить функциюв виде

то на основании линейности преобразования Фурье и с учетом выражений (4.48) и (4.49) спектральная плотность комплексного экспоненциального сигнала

Следовательно, комплексный сигнал обладает несимметричным спект­ром, представленным одной дельта-функцией, смещенной на частотувправо относительно.

7) Произвольная периодическая функция. Представим произвольную перио­дическую функцию (рис. 4.31, а) комплексным рядом Фурье

где - частота следования импульсов.

Коэффициенты ряда Фурье

выражаются через значения спектральной плотности одиночного импуль­са s (t ) на частотах (n =0, ±1, ±2, ...). Подставляя (4.55) в (4.54) и поль­зуясь соотношением (4.53), определяем спектральную плотность перио­дической функции:

Согласно (4.56) спектральная плотность произвольной периодической функции имеет вид последовательности-функций, смещенных друг от­носительно друга, на частоту (рис. 4.31, б). Коэффициенты при δ -функциях изменяются в соответствии со спектральной плотностьюодиночного им­пульсаs (t ) (пунктирная кривая на рис. 4.31,б).

8) Периодическая последовательность δ-функций (табл. 4.3, поз. 17). Спект­ральная плотность периодической последовательности –функций

определяется по формуле (4.56) как частный случай спектральной плотности периодической функции при = 1:

Рис.4.31. Произвольная последовательность импульсов (а) и её спектральная плотность (б)

Рис. 4.32. Радиосигнал (а), спектральные плотности радиосигнала (в) и его огибающей (б)

и имеет вид периодической последовательности δ -функций, умноженных на ко­эффициент .

9) Радиосигнал с прямоугольной огибающей. Радиосигнал, представленный на (рис. 4.32,а), можно записать как

Согласно поз. 11 табл.4.2 спектральная плотность радиосигнала полу­чается путем сдвига спектральной плотностипрямоугольной огибающей по оси частот на вправо и влево с уменьшением ординат в два раза, т. е.

Это выражение получается из (4.42) путем замены частоты на частоты– сдвиг вправо и- сдвиг влево. Преобразование спектра огибающейпоказано на (рис. 4.32, б, в).

Примеры расчета спектров непериодических сигналов приведены так же в .

Подразумевая под случайным процессом множество (ансамбль) функций времени, необходимо иметь в виду, что функциям, имеющим различную фор­му, соответствуют различные спектральные характеристики. Усреднение комплексной спектральной плотности, определяемойой (1.47), по всем функциям приводит к нулевому спектру процесса (при М[х (t )]=0 ) из-за случайности и независимости фаз спектральных составляющих в различных реализациях.

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

Спектральная плотность средней мощности представляет со­бой среднюю мощность, приходящуюся на 1 Гц при заданной частоте ω . Размерность функции W (ω) , являющейся отношением мощности к полосе частот, есть

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

Выделив из ансамбля какую-либо реализацию x k (t ) и ограничив ее дли­тельность конечным интервалом Т , можно применить к ней обычное преоб­разование Фурье и найти спектральную плотность X kT (ω). Тогда энергию рассматриваемого отрезка реализации можно вычислить с помощью форму­лы:

(1.152)

Разделив эту энергию на T , получим среднюю мощность k-й реализации на отрезке Т

(1.153)

При увеличении Т энергия Э кТ возрастает, однако отношение стремится к некоторому пределу. Совершив предельный переход , получим:

г
де

представляет собой спектральную плотность средней мощности рассматри­ваемой k-й реализации.

В общем случае величина W k (ω) должна быть усреднена по множеству реализаций. Ограничиваясь в данном случае рассмотрением стационарного и эргодического процесса, можно считать, что найденная усреднением по одной реализации функция W k (ω) характеризует весь процесс в целом. Опуcкая индекс k, получаем окончательное выражение для средней мощности случайного процесса

Для процесса с нулевым средним

(1.156)

Из определения спектральной плотности (1.155) очевидно, что W х (ω) является четной и неотрицательной функцией ω.

1.5.3 Соотношение между спектральной плотностью и ковариационной функцией случайного процесса

С одной стороны, скорость изменения х(t ) во времени определяет шири­ну спектра. С другой стороны, скорость изменения х (t) определяет ход ковариационной функции. Очевидно, что между W х (ω) и К х (τ) имеется тес­ная связь.

Теорема Винера - Хинчина утверждает, что К х (τ) и W x (ω) связаны между собой преобразованиями Фурье:

(1.157)

(1.158)

Для случайных процессов с нулевым средним аналогичные выражения имеют вид:

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

Рис.1.20. Широкополосный и узкополосный спектры случайного процесса; границы центральной полосы: ±F 1

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

Если в выражение 1.158 подставить W x (ω) = W 0 = const, то получим

где δ(τ) - дельта-функция.

Для белого шума с бесконечным и равномерным спектром корреляцион­ная функция равна нулю для всех значений τ, кроме τ = 0 , при котором R x (0) обращается в бесконечность. Подобный шум, имеющий игольчатую структуру с бесконечно тонкими случайными выбросами, иногда называют дельта-коррелированным процессом. Дисперсия белого шума бесконечно велика.

Вопросы для самопроверки

    Назовите основные характеристики случайного сигнала.

    Как связаны математически корреляционная функция и энергетический спектр случайного сигнала.

    Какой случайный процесс называется стационарным.

    Какой случайный процесс называется эргодическим.

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

    Какой сигнал называется аналитическим.