Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математические основы теории систем / Лекции / Павлова А.В. Математические основа теории систем. Лекции Часть 1.pdf
Скачиваний:
159
Добавлен:
02.06.2015
Размер:
10.66 Mб
Скачать

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

6.2. Спектральное представление сигналов. Разложение произвольного сигнала по заданной системе функций

 

 

Известно, что любая кусочно-непрерывная

функция f (t) , для

которой

ò

 

f (t)

 

2dt < ¥ (интеграл вычисляется по области

определения функцииf (t)) ,

 

 

 

 

 

 

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

¥

 

 

 

 

 

 

 

 

 

 

 

f (t) = C1j1(t) + C2j2 (t) + ... + Cnjn (t) + ... = åCk jk (t),

(6.1)

k =1

т. е. любой сигнал может рассматриваться как совокупность элементарных сиг-

налов jk (t) , умноженных на коэффициент Ck.

 

 

 

 

Множество

{jk (t)} составляет некоторую заранее заданную систему

функций, которые называются базисными.

 

 

 

 

Например:

 

 

 

(t) = tn-1, то f (t)

 

 

 

Если j (t) =1; j (t) = t,...,j

n

– полином по степеням t.

 

1

2

 

 

 

 

 

 

Если

= 1;

=

;

=

;…

=

;

=

, то f (t) – тригонометрический полином.

 

 

В качестве базисных выбирают такие функцииjk (t) , которые на заданном интервале [a, b] обладают свойством ортогональности, т. е.

b

òjmjk (t)dt = 0 при m ¹ k.

a

Выбор jk (t) из класса ортогональных даёт наилучшее приближение к ряду: среднеквадратичная ошибка (СКО) аппроксимации s на [a, b] минимальна при фиксированном числе N слагаемых ряда (6.1).

 

 

 

 

 

b

Если для каждой jk (t)

выполняется условие, что òjk2 (t)dt =1, то система

 

 

 

 

 

a

функций

jk (t) (k =

 

) ,

называется

ортонормированной. Величина

1, N

 

 

 

 

 

95

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

jk

 

 

 

= òjk2 (t) dt называется

нормой

функцииjk (t) . Квадрат нормы имеет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¥

 

смысл энергии сигнала E f =

 

 

 

 

f (t)

 

 

 

2 = ò

f 2 (t)dt.

 

 

 

 

 

 

 

 

 

Именно такая энергия выделяется в резисторе с сопротивлением1 Ом, если на его зажимах существует напряжение f (t).

Два сигнала u(t) и v(t) ортогональны, если их скалярное произведение, а, значит, и взаимная энергия равны нулю

¥

(u,v) = ò u(t)v(t)dt = 0.

Разнесённые во времени сигналы ортогональны.

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

1)точное разложение на простейшие ортогональные функции;

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

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

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

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

купностью коэффициентов Ck.

Умножим обе части выражения (6.1) на jn (t) и проинтегрируем его в пределах от a до b.

b

b ¥

 

 

 

 

 

ò f (t) × jn (t)dt = òåCk jk (t)jn (t)dt.

(6.2)

a

a k =1

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

Все слагаемые вида òCk jk (t)jn (t)dt при k n обращаются в нуль в силу

a

 

 

 

 

 

 

 

 

 

 

ортогональности функций jk

и jn , тогда в правой части остаётся одно слагаемое

 

b

 

 

 

2 .

 

 

òCnjn2 (t)dt = Cn

 

 

 

jn

 

 

 

(6.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

a

96

С учётом (6.3) выразим Cn из (6.2), при этом

 

1

 

 

 

 

b

 

Cn =

 

 

 

 

ò f (t)jn (t)dt.

(6.4)

 

 

jn

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

a

 

Ряд, в котором коэффициенты определяются в соответствии с выражением (6.4), называется обобщённым рядом Фурье.

6.3. Гармонический анализ периодических сигналов

Пусть задан периодический сигнал f (t) с периодом T. Выберем в качестве базисных систему ортогональных тригонометрических функций1; cos w t; sin w t; cos 2w t; sin 2w t;...; cos nw t;sin nw t. Интервал ортогональности совпадает

с периодом функции f (t) и равен T = 2wp .

При этом ряд (6.1) с коэффициентами (6.4) запишется следующим образом:

где

=

, jn = arctg

 

, а коэффициенты

и

определяются

 

выражениями:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если f (t) представляет

собой

чётную

относительноt

функцию, т. е.

f (t) = f (-t), то

в разложении в ряд

будут

присутствовать

только косинусои-

дальные члены,

так как все bn

обращаются в нуль. Для нечётной функции f (t)

ряд состоит только из синусоидальных членов, так как все

 

Таким образом функция f (t) представлена в виде суммы гармонических ко-

лебаний с разными частотами. Наименьшая частота w =

2p

это частота пер-

 

T

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

Спектры являются функциями, зависящими от номера гармоникиn, как независимой переменной. Графически спектры удобно изображать в виде -от резков длины и , проведённых перпендикулярно к оси w .

97

Число n может принимать как положительные, так и отрицательные значения, поэтому графики спектров и имеют смысл и при положительных и при отрицательных частотах nw . Следует отметить, что амплитудно-частотный спектр является чётно-симметричной, а фазочастотный спектр– нечётносимметричной функцией частоты nw . Поэтому при построении спектров достаточно отображать лишь их половину при nw > 0.

Пример 6.1. Найти и построить спектры прямоугольного сигнала, называемого меандром (рис. 6.8), который широко применяется в измерительной технике.

f (t) C

t

=

-C 0

 

 

 

 

 

 

Рис. 6.8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассматриваемая

функция

является

нечётной, поэтому коэффициенты

 

(как интеграл в симметричных пределах от нечётной функции).

Определим коэффициенты bn, учитывая, что w =

2p

.

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

T

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b =

2

2

 

f (t)sin wtdt =

4 2

 

C sin wtdt =

4C

(-cos nwt)

 

=

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò

 

T ò

 

 

n

 

T

 

 

 

 

 

 

 

 

Tnw

 

 

 

 

 

 

 

 

 

-

T

 

 

 

 

 

-

T

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

=

4C

(-cos nw

T

+1) =

4C

(1 - cos np) =

2C

(1 - cos np).

 

 

 

np

 

nwT

 

2

 

 

n2p

 

 

 

 

 

 

 

При чётных значениях n и n = 0 cosnπ = 1 и выражение в скобках обращается в нуль, так что нулевой составляющей и чётных гармоник данное выражение не содержит.

Вычислим коэффициенты bn:

при n =1

b =

2C

 

(1 - cos p) =

4C

;

 

 

 

 

 

 

 

 

 

 

 

1

 

 

p

 

 

p

 

 

 

 

 

 

 

 

при n = 3

b =

2C

(1 - cos 3p) =

4C

;

 

 

 

 

 

 

2

 

 

3p

 

 

 

 

3p

 

 

 

 

 

 

 

 

 

 

при n = 5

b =

2C

 

(1 - cos 5p) =

4C

 

и т. д.

 

 

 

 

3

 

 

5p

 

 

 

 

5p

 

 

 

 

 

 

 

 

 

 

98

Представим функцию f (t) в виде суммы гармонических составляющих

С увеличением числа суммируемых гармоник сумма ряда приближается к функции f (t).

Спектр коэффициентов тригонометрического ряда рассматриваемого сигнала определяется выражением и приведен на рис. 6.9, а. Фазоча-

стотный спектр определяется выражением

j (nw) = arctg

bn

= arctg¥

p

и пред-

 

 

 

n

an

2

 

 

 

 

ставлен на рис. 6.9, б. Для всех гармонических составляющих значения началь-

ных фаз .

а

б

Рис. 6.9

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

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

Вычислим среднюю мощность, используя разложение f (t) в ряд. Учитывая, что и , получим

99

Мощность первой гармоники и составляет 81% от полной мощ-

ности Pср = С2 . Мощность двух первых гармоник составит уже ≈90%. Сред-

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

Техническая ширина спектра – полоса частот, в которой сосредоточена основная часть энергии сигнала, как правило, это 90%. Для рассматриваемого сигнала ширина спектра составляет 3w.

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

f (t) = t, "t Î

é

T

,

T ù

f (t)

ê-

 

 

ú

2

2

 

ë

 

û

 

t

Рис. 6.10

Так как f (t) – нечётная функция, ряд Фурье содержит только синусоидальные составляющие с коэффициентами bn

T T

 

 

2

2

 

4

2

4 æ sin nwt

 

t cos nwt ö

 

bn

=

 

ò f (t) sin nwt dt =

 

 

ò C sin nwt dt =

 

ç

 

 

-

 

÷

=

T

T

T

2

2

nw

 

 

-

T

0

è

n

w

 

ø

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

2p

 

 

T

 

 

 

 

 

4

ç sin n

×

 

 

 

 

 

 

 

2

 

 

=

ç

 

 

 

 

 

T

-

 

 

 

 

 

2 4p 2

 

 

 

T ç

 

n

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 2

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

тогда

b =

T

;

b = -

T

;

p

2p

 

1

 

 

 

2

 

 

 

f (t) =

T cos n 2p × T 2 T 2

n2p

T

T

b3 = 3p ; b4 =

4C (sin wt + 1 p 3

ö

 

 

 

 

 

 

 

 

÷

4

æ

T 2

ö

 

T

 

÷ =

 

 

ç

 

- cos n

=

 

cos np,

 

 

 

 

÷

T

ç

4np

÷

 

np

 

è

ø

 

 

÷

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

T

4p ;...

sin 3wt + 1 sin 5wt + ...). 5

Амплитудно-частотный спектр An представлен на рис. 6.11, а, а фазоча-

стотный спектр Yn представлен на рис. 6.11, б.

100

An

 

 

 

Yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

б

Рис. 6.11

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

(рис. 6.12).

f (t)

, при

t

Рис. 6.12

Так как f (t) – чётная функция все коэффициенты bn = 0. Вычислим коэффициенты an, используя выражение

 

 

 

 

T

 

 

 

 

2

2

 

 

a

=

 

ò

 

f (t) cos nwt dt.

 

 

n

T

 

 

 

 

 

-

T

 

 

 

 

2

 

 

 

 

 

 

 

Тогда

 

 

 

 

T

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

4

2

 

 

4

 

T 2

 

 

 

 

 

 

 

 

2

a

=

 

ò

 

t dt =

 

×

 

 

 

 

 

 

 

 

 

 

 

0

T

 

 

T

 

2

 

0

 

 

 

-

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

=T . 2

101

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 2

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

æ cos nwt

 

 

t sin nwt ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

an

=

 

 

 

ò

 

f cos nwt dt

=

 

 

 

ç

 

 

 

 

-

 

÷

 

 

 

=

 

T

 

T

 

2 p2

nw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è n

 

 

ø

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

ç

cos np

 

 

 

 

 

1

 

 

÷

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

ç

 

-

 

 

 

 

 

÷ =

 

 

(cos np - 1).

 

 

 

 

 

 

 

 

 

 

 

 

2 4p2

 

 

 

2 4p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T ç

n

 

 

 

n

 

 

÷ n2 p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 2

 

 

T 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

Для различных n получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

=

 

T

(cos np - 1) =

2T

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

p2

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2 = a4 = a6... = 0 (для чётных n cosnπ = 1);

 

 

 

 

 

 

 

 

a3

=

T

 

(cos 3p - 1) =

 

 

2T

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9 p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9 p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представим f (t) в виде суммы гармонических составляющих

 

T

 

T

æ

cos 3w t

 

cos 25w t

ö

f (t) =

 

-

 

ç cos w t +

 

+

 

+ ...÷.

4

p2

9

25

 

 

è

 

ø

Амплитудно-частотный спектр представлен на рис. 6.13. Здесь учтено, что . Начальные фазы всех гармоник будут равны нулю, так как

5

Рис. 6.13

102

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

(рис. 6.14).

Cf (t)

τ

t

T

Рис. 6.14

Так как f (t) – чётная функция, все коэффициенты . Найдём коэффициенты :

Представим f (t) в виде ряда

Умножим и разделим выражение под знаком суммы на nwt , тогда

2

 

 

 

 

æ

¥

sin

nwt

ö

 

 

 

C t ç

 

÷

 

 

 

+ 2å

2

f (t) =

 

ç1

 

cos nwt ÷.

T

 

 

 

 

 

ç

n =1

 

nwt

 

÷

 

 

 

2

 

 

 

 

 

è

 

 

ø

Множитель типа

sin x

 

называется функцией отсчётов.

x

 

 

 

 

 

 

 

 

 

 

График этой функции представлен на рис. 6.15. sin x = 0 при x = kp, а

при x ® 0 sin x » x, тогда sin x = 1. Эта функция является огибающей спектра x

коэффициентов ряда.

103

Рис. 6.15

Амплитудно-частотный спектр представлен на рис. 6.16. Начальные

фазы всех гармонических составляющих равны нулю, так как arctg bn = 0 . an

 

Точки прохода

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

nwt

= kp, тогда

 

 

 

2kp

 

2kp

 

2

 

nw =

или W =

. На интервале

будет столько дискрет, во сколь-

 

 

t

 

t

 

 

 

 

 

ко раз

 

 

 

 

 

 

 

 

На рис. 6.16

 

 

= 6.

 

 

 

 

 

 

 

 

ω

w1 w3

Рис. 6.16

Величина N = Tt называется скважностью импульсной последовательно-

сти. Чем больше длительность импульса τ, тем уже спектр сигнала, и наоборот. чем меньше длительность сигнала, тем шире спектр.

104