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

dsd11-12 / dsd-11=ТКС / full

.pdf
Скачиваний:
128
Добавлен:
05.06.2015
Размер:
3.1 Mб
Скачать

ω1 = 2Tπ - основнаячастота.

Отыскание ряда Фурье для сложного гармонического сигнала, называется гармоническим или спектральным анализом. Составляющие этого ряда - циклические частоты, которые равны ω,2ω,3ω и т.д., называются соответственно первой (основной), второй, третьей и т.д. гармониками сложногосигнала.

Каждая гармоника ряда Фурье характеризуется амплитудой Am и начальной фазой ϕm, и для коэффициентов разложения (2.3а) можно запи-

сать: am=Amcosϕm, bm=Amsinϕm.

Откуда

A

=

 

a

2 + b2 , и

bm

= tgϕ

m

.

 

 

 

 

 

 

m

 

 

m

m

am

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда единственное решение (2.3)

представляется в более удобном

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S(t) =

a

0

+

 

A

 

cos(m ω t ϕ

 

)

 

( 2.4 )

 

 

 

m

 

2

 

 

 

 

 

m =

1

m

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример1.1 Найти спектральное преобразование последовательности прямо-

угольных импульсов спараметрами(см. рис.2.1а) четнымиотносительно t = 0 , T - период сигнала, τ - длительность импульса, A - амплитуда им-

пульса.

 

 

 

 

 

 

 

 

 

 

 

Отношение

T

 

= q

называетсяскважностьимпульсов.

τ

 

 

 

 

 

 

 

 

 

 

Используя(2.3а), имеем:

 

 

 

 

 

 

 

 

 

 

2 T

2 τ / 2

2

 

2A

;

 

a

=

 

S(t)dt =

 

Adt =

 

τ A =

 

 

 

 

 

 

 

0

T T

T τ./ 2

T

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

 

 

a

 

=

 

2 T

S(t)cos(mω t)dt =

 

2 τ / 2

Acos(mω t)dt

=

2A

sin mω t

1

 

 

 

 

τ / 2

=

 

ω

=

2π

 

=

A

m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

T

 

T

 

 

T

mω

 

 

Т

mπ

T

m

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

τ / 2

 

 

 

 

 

 

 

 

 

1

 

τ / 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

mπτ

 

mπτ

 

 

 

2A

 

 

mπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

sin

 

 

sin

 

 

 

 

 

=

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

mπ

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mπ

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

=

2

T

S(t)sin(mω t)dt

2

 

τ / 2

Asin(mω t)dt = −

2A

cos mω t

 

1

 

τ / 2

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

T

 

 

 

T

 

mω

 

 

 

 

 

 

 

 

 

 

 

 

m

 

T

 

 

1

 

τ / 2

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

τ / 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя выражение (2.3) для представления спектрального разложениясигнала, получим:

S(t)

=qA

=qА

 

2 A

 

 

2 A

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mπ

 

 

 

+ 0]

 

 

 

 

 

 

 

 

 

=

 

 

+

 

 

 

 

 

 

cos m

 

 

t

=

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

2

 

m =1 mπ

 

q

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mπ

 

 

 

 

 

 

 

 

 

mπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

A

 

 

q

 

 

 

 

sin q

 

 

 

 

 

sin

 

 

 

1

+ 2

 

 

 

 

 

cos m

 

 

t

=

 

 

1

+ 2

 

 

 

 

 

 

 

cos(m ωt)

=

 

 

 

mπ

 

T

q

 

 

mπ

 

 

 

m =1

 

 

 

 

 

 

 

m =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+ 2An cos(m ωt)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mπ

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

q

 

 

где

An =

 

 

. График такого спектрального разложения, для двух

mπ

 

 

 

 

 

q

крайнихслучаев, когда q >> 1(T >>τ) и q 1(T τ) , показаннарис.2.1б.

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

U

0

= 1 ;

U

= exp(jω1t);

 

 

T

 

1

T

 

 

 

 

 

U

2

= exp(2 jω1t);

U

3

= exp(3 jω1t);

 

T

 

T

 

 

 

 

 

 

 

 

19

 

 

 

T

U

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-τ/2

τ/2

 

τ

t

S(t)

S(t)

 

A/q

A/q

 

 

 

а)

2π / T

2π /τ

ω

q ≈1,

T ≈τ

ω

q >> 1,

T >> τ

 

 

 

 

короткиеимпульсы, идущиеоченьредко

 

 

длинныеимпульсы, идущиечасто

б)

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

20

 

 

 

 

 

 

 

U

k

= exp(kjω1t)

,

 

 

 

 

 

гдеk = 0, ±1, ±2...

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ТогдарядФурьепринимаетвид:

 

 

 

 

 

 

 

 

 

 

 

 

 

S(t) = 1

 

C

 

e jnω1t

,

(2.5)

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

T n = −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cn

= 1

T / 2

 

 

 

 

 

 

 

 

 

 

где

S(t) ejnω1tdt , Cn - коэффициентырядаФурье.

 

 

 

 

T T / 2

 

 

 

 

 

 

 

 

 

 

 

Разложение(2.5) иногдапредставляютвболееудобномвиде

 

 

 

 

 

 

 

 

 

S(t) =

C

 

e jnω1t ,

 

(2.5а)

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

n = −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

C

 

=

1

 

T / 2

S(t) ejnω1tdt

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T T / 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражения (2.5) и (2.5а) являются разложением ряда Фурье в комплекснойформе.

Так как n принимает значения от - до +, то из (2.5а) видно, что спектр сигнала содержит компоненты как с положительными, так и с отрицательными частотами. Обычно их объединяют в пары, т.е. по формулеЭйлера

 

 

 

 

 

 

 

 

j ω t + ϕ

 

 

 

 

j ω t + ϕ

 

 

jωt

jωt

 

 

 

e jω1t + C

 

ejϖ1t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ e

 

 

C

 

 

 

 

e

1

 

 

n +

 

C

 

e

1

n =

cosωt =

e

 

 

n

n

=

C

n

 

n

 

=

 

2

 

 

 

 

 

 

 

 

 

 

2

 

Cn

 

cos(nω1t + ϕn )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

напряжения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отрицательной

 

частоте

 

 

соответствует

вектор

(см.рис.2.2) из разложения (2.5), вращающийся по часовой стрелки, а положительнойчастотевекторвращающийсяпротивчасовойстрелки.

21

Im

Cn=|Cn|

+ω

Re

-ω

C-n=|Cn|

Рис.2.2 КомплексноепредставлениегармоникрядаФурье

2.2. Спектральнаяплотностьсигнала

Пусть S(t) - одиночный импульсный сигнал. Чтобы найти его спек-

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

 

 

S (t) =

C

 

e jnω1t

 

 

 

 

 

 

 

 

 

 

n = −∞

 

n

 

 

 

 

 

скоэффициентамиразложения:

 

 

 

 

 

 

 

 

 

 

 

1

T / 2

 

 

 

 

jn

ω

1t dt

C

n

=

S (t) e

 

 

T

 

 

 

 

 

T / 2

 

 

 

 

 

 

 

Чтобы вернуться к одиночному импульсу надо период Т устремить

в бесконечность Т → ∞ , тогда частоты nω1

и (n +1)ω1 будут сколь угод-

ноблизкидругкдругуиихможнозаменитьнатекущуючастоту ω

22

 

 

 

 

 

 

 

 

1

 

 

 

ejωt dt

 

( 2.6)

 

 

 

 

C

n

=

S(t)

 

 

 

 

 

T

 

 

 

 

 

 

 

 

− ∞

 

 

 

 

 

 

 

 

Тогда коэффициенты ряда Фурье образуют комплексно-

сопряженные пары C

n

= A e jϕn

и C

n

= A ejϕn , и каждой паре соответ-

 

 

 

 

n

 

 

 

 

 

 

 

n

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j nω t + ϕ

n

 

A e

j nω t + ϕ

 

 

 

 

 

 

A e

 

1

 

+

 

1

n =

 

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

 

 

или осуществляя предельный переход, т.е. когда nω1 мало отличает-

сяот (n +1)ω1 , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ϕ

 

 

 

 

 

 

 

 

 

 

 

 

j ωt + ϕ

 

j ωt

 

 

 

 

 

cos(ωt

+ ϕ

 

)

 

 

A e

n

+ A e

 

 

 

 

n = 2A

 

n

 

 

n

 

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

Рассмотрим бесконечно малый интервал Δω. В рамках этого интер-

вала будет содержаться N пар спектральных составляющих N = ω

= ωT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(так как ω1=2π/T) с коэффициентами разложения Сn (2.6 ), представляющими собой комплексную амплитуду соответствующего единичного гармоническогоколебания, т.е.

 

 

1

 

A

=

S(t) e jωtdt

(2.7)

T

ω

 

− ∞

 

 

 

 

 

Так как в рамках интервала Δω таких колебаний будет N, то полную комплексную эквивалентную амплитуду получим, домножив (2.7) на 2 и N, т.к. спектрсигналасодержиткомплексносопряженныепары.

 

 

2N

2

ω T

A

=

S(t) e jωtdt =

S(t) e jωtdt

 

 

 

ω

T

− ∞

2

π T

− ∞

 

 

23

Aω

 

1

 

 

=

S(t) e jωtdt ;

S(ω)=

S(t) e jωtdt

(2. 8 )

 

 

ω

π − ∞

 

−∞

 

Выражение S(ω) (2.8) называется спектральной плотностью сигна-

ла S(t) . Сучетомвышесказанного(2.8) принимаетвид

Aω

=

1

S(ω)

( 2.9)

 

π

ω

 

 

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

го интервала частот

ω и отвечающей ему амплитудой гармонического

сигналасчастотой ω0

лежащейвнутри ω.

Пример 2.2. Найти спектральную плотность прямоугольного импульсаизображенногонарис.2.3а.

Используясоотношение(2.8) имеем:

 

 

 

 

 

τ/ 2

 

 

τ/ 2

+ j sin ωt)dt

 

S(ω) =

 

U e jωtdt

 

U e jωtdt

= U (cosωt

=

 

 

−∞

 

− τ/ 2

 

 

− τ/ 2

 

 

 

 

 

 

 

τ/ 2

 

τ/ 2

 

 

 

 

τ/ 2

 

2U

 

 

τ

 

= U

cos(ωt)dt +

j sin(ωt)dt = 2U

cos(ωt)dt =

 

 

sin

ω

 

 

 

ω

2

 

− τ/ 2

 

− τ/ 2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть ξ =

ωτ

 

, тогда

S(ξ) =

U τ

sin ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

ξ

 

 

 

 

 

 

 

График нормированной спектральной плотности

S(ξ)

 

представлен

нарис.2.3б.

24

S(ξ)/S(0)

S

U

t

-τ/2

a)

τ/2

-π

π

ξ

-2π

 

2π

 

б)

 

Рис.2.3 Cпектральнаяплотностьпрямоугольногосигнала: а- форма ипараметрыпрямоугольногосигнала; б– графикнормированнойспектральнойплотности

25

3.Случайныесигналыишумы

3.1Математическийаппаратслучайныхсигналов

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

Трезистора), топодразумеваютегостатистическуюхарактеристику.

Втеориивероятностейвведеныпонятия:

-функции распределения случайной величины F(x) , т.е. вероятность

того, что случайная величина x из множества X примет значение, равноеилименьшеечем x0 . F(x) = P(xi x0 )

Дляфункциираспределениясправедливыдвапредельныхравенства:

F ( −∞ ) = 0 F (+ ∞)=1 (3.1)

- производная от функции распределения - есть плотность вероятности:

 

P(x)=

dF(x)

.

 

(3.2)

 

 

 

 

 

dx

 

 

Плотностьвероятностиудовлетворяетусловиюнормировки,

Pi =1

т.е.

P(x)dx =1

(3.3)

i

 

 

 

− ∞

 

Дляслучайныхвеличин x измножества X введеныпонятия:

Математическое ожидание ( m ) - среднее значение, которое может приниматьслучайнаявеличина x :

26

_

 

(3.4)

m = x

=

xP(x)dx

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

δ2x = ∫ (x m)2 P(x)dx (3.5)

− ∞

Пример1 Равномерноераспределение.

Пусть случайная величина x может принимать значения из интервала xi x x2 , причем вероятности попадания в любые внутренние интервалы равны.

 

 

 

 

 

 

 

 

0, x

< x1

 

 

 

 

 

 

 

P(x)=

1

 

, x x x

 

 

 

 

2

 

 

1

 

x2 x1

 

 

 

 

 

 

 

 

 

0, x

> x

 

 

 

 

2

 

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

 

 

 

 

 

 

 

 

 

0, x < x1

xТЕК

 

x

x1

 

F(x) =

P(x)dx =

, x1 x x2

 

 

− ∞

x2

x1

 

 

 

 

x > x2

 

 

1,

 

 

 

 

 

 

График плотности вероятности и функции распределения показан на рис.3.1.

27

Соседние файлы в папке dsd-11=ТКС