7. Спектральный анализ периодических

СИГНАЛОВ

7.1. Разложение периодических сигналов

по ортогональным функциям

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

Рис. 7.1

При расчете разнообразных сигналов удобно представить их взвешенной суммой заданных функций времени вида

, (7.1)

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

Чтобы разложение в ряд (7.1) было взаимно однозначным

144

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

(7.2)

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

Для ортогонального базиса коэффициенты разложения определяются выражением

. (7.3)

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

, (7.4)

полиномы Чебышева, Лагранжа, Эрмита и др. В цифровой технике применяют ортогональные дискретные функции Уолша, Радамахера.

7.2. Ряд Фурье

Ряд Фурье для действительной периодической функции времени является ее разложением по ортогональному базису (7.4) и имеет вид

145

(7.5)

Компоненту ряда Фурье вида

(7.6)

называют -й гармоникойсигнала,

(7.7)

- частота первой гармоники,-постоянная составляющая сигнала,

, (7.8)

-амплитуда -й гармоникисигнала,

, (7.9)

, (7.10)

, (7.11)

-начальная фаза -й гармоникисигнала,

146

(7.12)

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

7.3. Спектры амплитуд и фаз периодического сигнала

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

, (7.13)

то есть двумя в общем случае бесконечными наборами чисел.

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

, (7.14)

а второй – спектром фаз,

. (7.15)

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

, (7.16)

147

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

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

Рис. 7.2 тельности прямоугольных

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

, (7.17)

а амплитуды синфазной и квадратурной составляющих -

, (7.18)

. (7.19)

Для амплитуды и начальной фазы -й гармоники получим

, (7.20)

(7.21)

148

Графики спектров амплитуд и фаз при условии ,мс,мс показаны на рис. 7.3а и рис. 7.3б соответственно. Каждая гармоника отображается вертикальной линией, длина которой равна величине амплитуды или фазы.

Рис. 7.3

Переменная является номером гармоники. Ее можно рассматривать как нормированную частоту гармоники,

,

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

Спектры имеют дискретный (линейчатый) характер, интервал частот междусоседнимигармоникамиодинаков и

149

равен .

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

Рис. 7.4 амплитуд гармоник. Формулу

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

, (7.22)

где непрерывная переменная.

В примере (7.20) получим

, (7.23)

график показан на рис. 7.3а пунктирной линией. Характерной особенностью огибающей спектра амплитуд сигнала рис. 7.2 является наличие точек с нулевым значением (нулей огибающей), определяемых из уравнения

, (7.24)

решение которого имеет вид

150

, (7.25)

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

7.4. Синтез сигнала по его спектру

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

. (7.26)

Соответствующие кривые при ,ипоказаны на рис. 7.5а, рис. 7.5б, и рис. 7.5в соответственно.

Как видно, с увеличением форма синтезированного сигнала приближается к исходной (рис. 7.2).

7.5. Ряд Фурье в комплексной форме

Гармоники сигнала могут быть представлены своими комплексными амплитудамив виде

, (7.27)

тогда исходный сигнал можно представить в виде ряда Фурье,

. (7.28)

151

Рис. 7.5

Амплитуда -й гармоникиравна модулю комплексной амплитуды,

, (7.29)

а ее начальная фаза – аргументус противоположным знаком,

152

(7.30)

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

Рассмотрим сигнал, показанный на рис. 7.6, тогда

(7.31)

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

7.6. Влияние формы сигнала на спектры амплитуд и фаз

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

На рис. 7.6 показана последовательность прямоугольных импульсов длительностьюи с амплитудой 1 на интервале периодав нормированных координатах времени. Для этого сигнала характерныкрутые (с нулевой продолжи-

153

тельностью) фронт и срез импульса. Величину

(7.32)

называют скважностьюимпульсов. На рис. 7.7 приведены спектры амплитуд (рис. 7.7а) и фаз (рис. 7.7б)

Рис. 7.6 при .

Рис. 7.7

На рис. 7.8 приведены аналогичные зависимости при .

Рис. 7.8.

154

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

Рассмотрим трапециидальный импульс, программа исследование которого в среде MathCADпоказана на рис. 7.9. Спектральный анализ проводится с помощью стандартной процедуры спектрального анализаfft(s). Она построена на основе алгоритма быстрого преобразования Фурье (БПФ) и позволяет получить комплексные коэффициенты, с помощью которых комплексная амплитуда-й гармоники определяется выражением

. (7.33)

Период сигнала выбран равным 1,- число отсчетов сигнала на периоде. Результаты расчета спектров амплитуд и фаз приведены в листинге программы на рис. 7.9 (повторите расчеты самостоятельно для различных параметров сигнала).

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

155

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

, (7.34)

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

Рис. 7.9.

График сигнала и его спектр амплитуд показаны в листинге программы на рис. 7.10. Как видно, спектр «гладкого» сигнала сосредоточен в области нижних частот, высшие гармоники практически отсутствуют.

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

156

Рис. 7.10.

7.7. Свойства спектров сигналов

Свойства спектров сигналов часто формулируются в виде теорем.

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

157

Взяв модули левой и правой частей (7.30), получим

, (7.36)

то есть спектр амплитуд не изменяется при задержке сиг-

нала во времени.

Вычислим аргументы обеих частей выражения (7.30),

, (7.37)

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

Для доказательства теоремы смещения запишем

. (7.38)

Проведем замену переменных , тогда получим

. (7.39)

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

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

, (7.40)

158

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

, (7.41)

а начальная фаза равна 0 или в зависимости от знака.

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

, (7.42)

комплексная амплитуда -й гармоникиявляетсямнимым числом,

, (7.43)

а начальная фаза равна 0 или в зависимости от знака.

Эти свойства иллюстрирует пример четного сигнала на рис. 7.2, для которого имеет место равенство (7.19). Его фазовый спектр со значениями 0 или показан на рис. 7.3б.

Рассмотрим комплексные спектры двух сигналов (рис. 7.11а) и(рис. 7.11б), и их сумму(рис. 7.11в).

Рис. 7.11.

Сигнал получен изсдвигом во времени на

159

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

Комплексная амплитуда -й гармоникиопределена ранее (7.31) и равна

(7.44)

По теореме смещенияможно найти комплексную амплитуду-й гармоники сигналав виде

. (7.45)

Тогда согласно свойству линейностикомплексная амплитуда-й гармоники сигналаравна

(7.46)

160

С другой стороны, при прямом вычислении (проведите расчеты самостоятельно) комплексная амплитуда-й гармоники сигналаравна

, (7.47)

что полностью совпадает с (7.46).

7.8. Мощность периодического сигнала

Пусть имеется сигнал (ток или напряжение) в сопротивленииОм, тогда средняя мощность сигнала равна

. (7.48)

Эту же величину можно выразить через гармоники сигнала с помощью равенства (теоремы) Парсеваляв виде

. (7.49)

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

. (7.50)

7.9. Ширина спектра

Как видно по графикам спектров амплитуд рассмотренных сигналов, в целом наблюдается тенденция уменьшения

161

амплитуд гармоник с ростом их номера (частоты). Графики на рис. 7.5 показывают, что форма сигнала определяется сравнительно небольшим числом гармоник. Все это свидетельствует

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

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

. (7.51)

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

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

выражения

. (7.52)

В результате можно определить ширину спектра в виде

. (7.53)

162

В качестве примера рассмотрим последовательность прямоугольных импульсов, показанную на рис. 7.2 со спектром амплитуд, показанном на рис. 7.3а. (скважность импульсов ) Зависимость нормированной мощностиот числа учитываемых гармоник показана на рис. 7.12. Как видно, функцияявляетсянеубывающейи достигает уровняпри(,), тогда ширина спектра сигнала определяется выражением (7.44).

Рис. 7.12

Этот же график в области значений от 0,9 до 1 показан на рис. 7.13. С ростом кривая очень медленно приближается к 1 и достигает значения 0,99 уже при.

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

163

пример, рис. 7.2) ширина спектра определяется выражением

(рад/с) или(Гц) (7.54)

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

Рис. 7.13

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

Оценки ширины спектра можно выразить через число гармоник,

, (7.55)

где требуемое число гармоник, равное

. (7.56)

164

На практике чаще всего используются соотношения с единичным множителем вида

В рассмотренном примере сигнала на рис. 7.2 скважность и для обеспечения 90% мощности необходимо учитыватьгармоник (рис. 7.12), по оценке (7.58) требуется учитывать 10 гармоник.

7.10. Задания для самостоятельного решения

Задание 7.1. Определите и постройте графики спектров амплитуд и фаз сигналов вида:

,

,

,

,

.

Задание 7.2. Определите спектры амплитуд и фаз сигналов, показанных на рис. 7.14, постройте их графики. Проведите расчет ширины спектра при,и, сравните полученные результаты.

165

Рис. 7.14

Задание 7.3. С помощьютеоремы смещенияпроведите расчет спектров амплитуд и фаз сигнала, показанного на рис. 7.14а, воспользовавшись результатами, полученными для сигнала на рис.7.2.

Задание 7.4. Определите спектры амплитуд и фаз сигнала, показанного на рис.7.15, постройте их графики. Сравните спектр амплитуд со спектром гармонического сигнала, проанализируйте результаты.

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

Рис. 7.15

Задание 7.5. Определите спектры амплитуд и фаз сигналjd, показанного на рис.7.16, постройте их графики.

166

Рис. 7.16

Задание 7.6. Определите спектры амплитуд и фаз сигнала, показанного на рис.7.17, постройте их графики.

Рис. 7.17

Проведите тот же расчет, представив сигнал на рис. 7.17 в виде суммы двух импульсных последовательностей, показанных на рис. 7.18, и используя свойство линейности.

Рис. 7.18

167