Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Введение в теорию сигналов. Цифровая обработка сигналов. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

M2 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

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