Введение в теорию сигналов. Цифровая обработка сигналов. Учебное пособие
.pdfM2 s; s(t) s(t kT ), t ( , ), k 0,1, 2, .. –
множество периодических сигналов.
Для анализа и обработки информации, содержащейся в сигнале, необходимо из глобального множества M выделить сигналы с определенными параметрами, чтобы, например, сравнивать их между собой или производить какие-либо другие операции. С этой целью в M выделяюттакназываемыеметрическиепространствасопределеннымисвойствамииединицамиизмерений. Например, прианализесигналовнаиболее удобно использовать пространство L2 . Это пространство кусочно-
непрерывных, интегрируемых со своим квадратом функций, т. е. таких, что
[s(t)]2 dt .
Метрическое пространство определяется своей метрикой (расстоянием (s1, s2 ) между двумя элементами этого пространства), обладающей следующими свойствами:
(s1, s2 ) 0, (s1, s2 ) 0, если s1 s2 ;(s1, s2 ) (s2 , s1);
(s1, s2 ) (s1, s2 ) (s2 , s3 ) .
Метрическое пространство M будем называть линейным L , если в нем определены некоторые операции, например сложение или умножение, в результате которых образуется новый элемент, принадлежащий этому жепространству. Линейноепространство L должно удовлетворять следующим требованиям.
1. Существование нулевого элемента такого, что s 0 s .
2. Для любых s L существует их сумма
s3 s1 s2 ,
11
которая тоже принадлежит L , причем операция суммирования: а) коммутативна – s1 s2 s2 s1 ,
б) ассоциативна – s1 (s2 s3 ) (s1 s2 ) s3 , в) однозначна – s1 ( s1) 0 .
3. Существование множества скалярных элементов a , для которых s2 a s1 L , причем операция умножения:
а) ассоциативна – a(b s) ab s;
б) дистрибутивна – a(s1 s2 ) a s1 a s2 , (a b) s a s b s; в) пропорциональна – 1s s, 0 s 0 .
В L всегда можно выделить некоторое множество элементов N , для которых выполняется условие
N
ai i 0, i 0
причем если все ai 0, i [0, N ] , то система функций i , i [0, N ] , будет линейно независимой. Если же хотя бы один ai 0, i [0, N ] , то система функций i , i [0, N ] , – линейно зависима.
Система линейно независимых элементов представляет собой базис пространства L . Посредством базиса в L можно единственной комбинацией выразить элемент s этого пространства в виде
s cr r .
r 0
Совокупность коэффициентов cr , r [0, ] , в теории сигналов
называют спектром элемента s в базисе { } .
Линейные пространства, как правило, имеют не единственный базис. Выбор базиса определяется простотой и удобством его использования.
Основнымиметрическимипараметрамилинейногопространстваявляются норма, метрика и скалярное произведение элементов этого пространства.
Определим указанные параметры в пространстве сигналов L2 .
12
Норма s(t) , t [0, T ] (sk , k [0, N 1]) |
определяется в виде |
|||||||||||||||
|
|
|
|
|
|
T |
|
|
|
|
N 1 |
|||||
|
s(t) |
|
|
s2 (t) dt , |
sk |
|
|
sk2 . |
||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
k 0 |
|
Метрика (s1(t), s2 (t)), |
|
или расстояние между s1(t), s2 (t), t [0, T ] , |
||||||||||||||
(s1k , s2k , k [0, N 1]), вычисляется как норма разности элементов: |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
(s1(t), s2 (t)) |
|
|
|
s1 (t) s2 (t) |
|
|
|
|
[s1(t) s2 (t)]2 dt , |
|||||||
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
N 1 |
||
(s1k , s2k ) |
|
s1k s2k |
|
|
|
|
|
[s1k s2k ]2 . |
||||||||
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
N k 0 |
|
||
Скалярное произведение в общем случае отражает степень связи (сходства) по форме и положению элементов в пространстве и обозначается как (s1(t), s2 (t)) , (s1k , s2k ) причем
T
(s1(t), s2 (t)) s1(t) s2 (t) dt
0
N 1
, (s1k , s2k ) s1k s2k . k 0
1.5. РЯДЫ ФУРЬЕ
Система функций 0 (t), 1(t), 2 (t), ... , n (t) L2 на [a, b] называется ортогональной, если
|
|
b |
|
|
||||
( i (t), j (t)) i (t) j (t) dt 0, i j , |
i, j 0,1, 2, ... , n . |
|||||||
|
|
a |
|
|
||||
При выполнении дополнительного условия |
|
|||||||
|
|
|
|
|
|
b |
|
|
|
i (t) |
|
|
|
|
[ i (t)]2 dt 1, |
i 0,1, 2, ..., n |
|
|
|
|
||||||
a
система i (t), i [0, n], называется ортонормированной.
13
Произвольная кусочно-непрерывная функция s(t) L2 на [a, b] мо-
жет быть аппроксимирована обобщенным рядом Фурье по ортогональному (ортонормированному) базису r (t), r [0, ) :
s(t) Cr r (t).
r 0
Практическое использование рядов Фурье предполагает их конечность
l
s(t) Cr r (t),
r0
врезультате чего возникает ошибка аппроксимации из-за отбрасывания остаточного члена
Rl (t) Cr r (t) .
r l 1
Минимизация ошибки в среднеквадратичном смысле позволила дать аналитическое выражение коэффициентов Фурье, которое для ортонормированного базиса имеет вид
b
Cr s(t) r (t) dt .
a
Ортогональная система функций r (t), r [0, ), называется полной, если для нее выполняется равенство Парсеваля
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 , |
||
Cr2 |
|
|
|
r (t) |
|
|
|
|
|
|
s(t) |
|
|
|
||
|
|
|
|
|
|
|
|
|||||||||
r 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
которое служит условием среднеквадратичной сходимости ряда Фурье к аппроксимируемой функции. В равенстве Парсеваля
b

s(t)
2 s2 (t) dt E
a
14
и если a 0 , а b T , то
|
|
|
|
l |
|
|
|
||
|
|
E Cr2 |
|
|
|
||||
|
|
|
|
r 0 |
|
|
|
||
представляет собой энергию сигнала s(t) , а |
|
|
|||||||
|
E |
|
1 |
l |
|
|
|
|
2 |
|
|
|
|
|
|||||
Pср |
|
|
|
Cr2 |
|
|
r |
(t) |
|
|
|
|
|||||||
|
T |
|
T r 0 |
|
|
|
|
|
|
среднюю на T мощность сигнала.
В задачах аппроксимации сигналов рядами Фурье можно выделить два направления.
1.Приближение сложного сигнала минимально возможным числом членов ряда при заданной погрешности аппроксимации.
2.Представление сложного сигнала в виде простейших функций. Использование в гармоническом анализе математических моделей
сигналов в виде рядов Фурье преследует вторую цель, а в качестве простейших функций выбираются ортонормированные системы тригонометрических функций sin t и cos t . Перечислим причины такого вы-
бора.
1. Гармонические колебания sin t и cos t являются простейшими функциями, определенными при всех t ( , ), и не поддаются дальнейшему разложению на составляющие.
2.Гармонические колебания – единственные функции времени, сохраняющие свою форму при прохождении через любую линейную стационарную систему. Изменяться могут только амплитуда и фаза.
3.Аппроксимация сложного сигнала гармоническим рядом Фурье используетхорошоразработанныйматематическийаппаратдляанализа передачи гармонических колебаний через линейные системы.
Следует отметить и факт представления сигналов их математической моделью в виде гармонического ряда Фурье с конечным числом слагаемых. Из теории приближения функций известно, что если на конечном интервале аппроксимируемая функция имеет точки разрыва первого рода (например, случай импульсного сигнала), то гармонический ряд Фурье в точке разрыва сходится к среднему арифметическому
15
значению функции слева и справа от точки разрыва. Причем в указаннойточкевозникаютколебания, затухающиепомереудаленияотточки разрыва влево и вправо. Это явлениев теории аппроксимации называют явлением Гиббса, оно увеличивает погрешность аппроксимации, и поэтому его пытаются устранить, что удается сделать лишь частично путем умножения коэффициентов Фурье на специальные сглаживающие множители (окна) Кайзера, Хэмминга, Бартлета, Парзена и др.
1.6. ГАРМОНИЧЕСКИЙ АНАЛИЗ ПЕРИОДИЧЕСКИХ СИГНАЛОВ
Запишем математическую модель сигнала ского ряда Фурье:
s(t) a0 |
|
|
|
|
|
|
|
2 t |
|
|
|
2 t |
A0 |
|
|
||||||
ar cos r |
br sin r |
|
|
||||||||||||||||||
|
|||||||||||||||||||||
|
2 |
r 1 |
|
|
|
T |
|
|
|
T |
|
|
|
2 |
|
|
|
||||
|
|
2 |
|
Т/2 |
s(t)cos (r t) dt , br |
|
2 |
Т/2 |
|||||||||||||
|
ar |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
T |
|
|
|
T |
|
|
|||||||||||||
|
|
|
|
Т/2 |
|
|
|
|
|
|
|
Т/2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
A |
a2 b2 |
, |
|
r |
arctg br , a |
r |
A cos |
r |
, |
||||||||||||
r |
r |
r |
|
|
|
|
|
ar |
|
r |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
s(t) в виде гармониче-
|
|
2 |
|
, |
|
||||
Ar cos r |
T |
t r |
||
r 1 |
|
|
|
|
s(t)sin (r t) dt ,
br Ar sin r , 2T .
Используя формулу Эйлера
|
cos exp ( j ) exp ( j ) |
, |
|||
|
|
|
|
2 |
|
полагая |
r t и учитывая |
обозначение |
ˆ |
||
Ar Ar exp( j r ), ряд |
|||||
Фурье можно записать в комплексной форме, а именно: |
|||||
|
s(t) |
1 |
|
ˆ |
|
|
|
|
Ar exp( j r t). |
|
|
|
|
2 r |
|
|
|
ˆ |
называют комплексной амплитудой r -й гармоники, содержа- |
||||
Здесь Ar |
|||||
щей информацию об амплитуде и начальной фазе этой гармоники.
16
При r 0 |
имеет место |
|
ˆ |
|
|
|
ˆ |
|||||
|
A r Ar exp ( j r ) Ar – комплексно-со- |
|||||||||||
пряженная с |
ˆ |
|
|
|
|
|
|
|
|
|
|
|
Ar величина. |
|
|
|
|
|
|
||||||
Комплексную амплитуду можно определить непосредственно по |
||||||||||||
сигналу s(t) |
в виде |
|
|
|
|
|
|
|
|
|||
ˆ |
|
|
|
|
|
Ar (cos r |
j sin r ) ar j br . |
|||||
Ar Ar exp( j r ) |
||||||||||||
Подставляя в последнюю формулу выражения для ar , br и объеди- |
||||||||||||
няя интегралы, получаем |
|
|
|
|
|
|
|
|||||
|
|
|
|
ˆ |
|
2 |
|
Т/2 |
s(t)exp( j r t) dt. |
|||
|
|
|
|
Ar |
|
|
|
|||||
|
|
|
T |
|
||||||||
|
|
|
|
|
|
Т/2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||
Если |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A0 |
|
|
|
|
|
2 |
|
|
A0 |
ˆ |
s(t) |
|
|
|
Ar |
cos r |
|
t r |
|
|
Re Ar exp( j r t), |
||
|
2 |
|
T |
2 |
||||||||
|
|
|
r 1 |
|
|
|
|
r 1 |
||||
то полученнаяпара выраженийпредставляет собойсоответственно прямое и обратное непрерывные преобразования Фурье, осуществляющие связь между временным и частотным представлением сигналов.
1.7. ГАРМОНИЧЕСКИЙ АНАЛИЗ НЕПЕРИОДИЧЕСКИХ СИГНАЛОВ
Пусть s(t) – непериодическийсигнал(рис. 1.2). Образуемвиртуаль-
ный периодический сигнал u(t) |
(рис. 1.3) с периодом T : |
||||||
u(t) |
|
|
|
1 |
|
ˆ |
|
s(t kT ) |
|
|
Ar exp( j r t). |
||||
|
k |
|
|
2r |
|
||
Для периодического сигнала u(t) |
справедливо |
||||||
ˆ |
ˆ |
2 |
Т/2 |
|
|
|
|
Ar A(r ) |
|
|
u(t)exp( j r t) dt |
||||
T |
|||||||
|
|
Т/2 |
|
|
|
||
|
|
|
|
|
|
||
17
Рис. 1.2
и, если T , то
lim u(t)
T
|
Рис. 1.3 |
|
|
|
s(t k T ) s(t) |
|
k |
для конечных значений t, и тогда ряд Фурье
|
|
|
|
s(t) |
1 |
ˆ |
|
||
|
|
|
|
|
Ar exp( j r t) . |
|
|||
|
|
|
|
|
|
|
2 r |
|
|
Подставив в последнюю формулу выражение для |
ˆ |
||||||||
Ar , получим |
|||||||||
|
|
1 |
|
2 |
Т/2 |
|
|
|
|
s(t) |
|
|
|
|
|
s(t)exp( j r t) dt exp( j r t) |
|||
|
|
|
|||||||
|
|
2 r T |
Т/2 |
|
|
||||
|
|
|
|
|
|
|
|
Т/2 |
|
|
|
|
exp( j r t) |
s(t)exp( j r t) dt. |
|||||
|
|
||||||||
|
2 r |
|
|
Т/2 |
|
||||
При T Ar |
0 , |
т. е. получаем бесконечно малые амплитуды |
|||||||
гармонических составляющих, сумма которых дает исходную непериодическую функцию s(t) , заданную на [ T / 2,T / 2] . С другой стороны,
при T 0 , а – это расстояние (шаг) между спектральными линиями Ar . В результате спектр вместо линейчатого (дискретного)
становится сплошным (непрерывным), а это позволяет считать
d , r ,
18
и тогда
s(t)
Обозначим
|
1 |
exp( j t) |
Т/2 |
s(t) exp( j t) dt d . |
||
2 |
|
|||||
|
|
|
|
|||
|
|
|
|
Т/2 |
|
|
S( j ) |
Т/2 |
s(t) exp( j t) dt |
|
||
|
Т/2 |
|
и, если об этом T ничего не оговаривается (в общем случае), то
S( j ) s(t) exp( j t) dt .
S( j ) будем называть спектральной плотностью или комплексной спектральной характеристикой сигнала s(t) . Тогда
s(t) 21 S( j )exp( j t) d .
Последние два соотношения образуют пару прямого и обратного преобразования Фурье. При этом можно установить связь между ком-
плексной амплитудой r -й гармоники ˆ и комплексной спектральной
Ar
характеристикой S( j ) , которая выражается как
ˆ |
2 |
S( j ), |
|
Ar |
|
||
T |
|||
|
|
а это значит, что для определения комплексной амплитуды любой гар-
моники ˆ периодического сигнала достаточно вычислить спектраль-
Ar
ную характеристику S( j ) на одном периоде, взять ее значение на частоте r искомой r -й гармоники и умножить на 2/Т.
19
1.8. НЕПРЕРЫВНЫЕ ПРЕОБРАЗОВАНИЯ ФУРЬЕ
Полученные в предыдущем разделе выражения пары прямого и обратного преобразования Фурье в общем случае будем записывать в следующем виде:
|
|
|
|
||
S( j ) Fr{s(t)} |
s(t)exp( j t) dt, |
||||
|
|
|
|||
s(t) Fr 1{S( j )} |
1 |
|
|
S( j ) exp( j t) d . |
|
|
|
||||
2 |
|||||
|
|
|
|||
Однако при решении практических задач приходится оперировать с конечнымиреализациямисигналов. Поэтомуеслипреобразуемыесигналы удовлетворяют условиям финитности, то, задавая уровни усечения s и S , определим эффективные длительности Ts – сигнала во
времени |
и S |
– |
частотной |
характеристики. |
|
Если sm |
sup |
|
s(t) |
|
|
||||||||||||||||||
|
|
||||||||||||||||||||||||||||
и |
|
s(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t [0, ) |
|
|
|
|
||
|
|
s sm |
при |
t Ts , то |
Ts – |
эффективная |
|
длительность |
|
сиг- |
|||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||
нала s(t). Аналогично, если Sm |
|
sup |
|
S( ) |
|
и |
|
S( j ) |
|
|
|
S Sm |
|
при |
|||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||
S , |
то S – |
|
|
|
|
( , ) |
|
|
|
|
|
|
|
|
|
S( j ) |
|
|
|
|
|
|
|||||||
эффективная |
|
длительность |
|
спектра |
|
|
|
|
сиг- |
||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||
нала s(t) |
(рис. 1.4, 1.5). В результате приходим к приближенным ана- |
||||||||||||||||||||||||||||
логам пары преобразований Фурье: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Ts |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
S( j ) Fr{s(t)} |
s(t)exp( j t) dt , |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
s(t) Fr 1{S( j )} |
|
|
|
S S( j ) exp( j t) d . |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Непрерывное преобразование Фурье обладает определенными свойствами, сходными с преобразованиями Лапласа, причем эти свойства рассматривались ранее в курсах «Математического анализа»
20
