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

2799-up_ch2

.pdf
Скачиваний:
26
Добавлен:
11.05.2015
Размер:
7.25 Mб
Скачать

31

На рисунках 3.2 и 3.3 представлены графики модуля и аргумента комплексной спектральной плотности сигнала.

Е 2

S

Тоc

4

 

2

 

2

 

4

 

 

 

 

 

 

 

Рисунок 3.2 – График модуля спектральной плотности

( )

 

 

2

 

 

 

 

 

2

4

 

 

 

 

4

 

2

 

-2

 

 

 

 

а)

( )

 

 

2

 

 

4

 

 

 

2

 

 

 

 

 

4

 

2

 

 

-2

 

 

 

 

 

б)

Рисунок 3.3 – Графики аргумента комплексной спектральной плотности: ненормированный (а) и нормированный (б) относительно 2 n

Анализируя рисунки 3.2 и 3.3, можно сказать, что сдвиг сигнала во временной области не изменяет значения модуля спектральной плотности, а приводит к изменению фазовой характеристики.

32

3.2Дискретизация аналогового сигнала по времени

Руководствуясь теоремой Котельникова, частоту дискретизации д выбираем так, чтобы д 2 в . Для определения верхней частоты в спектра

сигнала воспользуемся пороговым критерием: для частот выше «верхней» модуль спектральной плотности не превышает уровня 0,1 от максимального значения:

 

 

 

E E

 

sin

в

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

0,1

 

 

 

 

 

 

 

 

 

.

 

 

2

2

 

 

 

в

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

Так как sin

в

не превышает 1,

 

 

 

положим в 20 6 . Таким

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

образом частота дискретизации д и интервал дискретизации Тд будут равны

д 2 в 12 , Т д 2 д в Тос 12 .

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

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

Число степеней свободы N будет равно

N Tос 12.

Tд

На рисунке 3.4 приведены отсчеты прямоугольного импульса.

1,2

Е

0,8

sn 0,6 0,4

0,2

0,2 0

1

2

3

4

5

6

7

8

9

10

11

 

 

 

 

 

 

 

 

 

 

 

n

Рисунок 3.4 – Отсчеты исходного аналогового сигнала

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

33

Е{1,1,1,1,1,1,0,0,0,0,0,0}.

Запишем математическую модель дискретного сигнала sд t с помощью дискретных функций Хэвисайда:

sд t E n E n 6 .

3.3 Расчет спектральной плотности дискретизированного сигнала

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

 

 

 

 

 

 

 

 

 

Sд

 

 

 

E

Tос t t T

... t 5T

 

e j t dt .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя фильтрующее свойство -функции, получим

 

 

 

 

 

 

Sд ( )

 

E

(e j0 Tд

 

 

e j Tд e j2 Tд e j3 Tд e j4 Tд e j5 Tд )

 

 

 

 

 

 

 

 

 

N

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

j

5

T

 

 

 

 

j

3

T

 

j

1

T

 

j

1

T

 

 

j

3

T

 

 

 

j

5

T

j

5

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

2

 

д e

2

 

 

 

 

д e

 

2

 

д

e

2

 

д

e

 

2

д

 

e

2

д e

2

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

5

 

 

 

 

 

3

 

 

 

 

 

 

 

1

 

 

 

j

5

Tд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

Tд

cos

 

Tд

cos

 

 

Tд e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

2

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sд

 

 

2E

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

j

5

Tд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

T

2cos 2 T

1 e

 

2

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

12

 

 

 

2

 

д

 

 

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проанализируем

полученное

выражение.

 

Для

 

этого

 

 

на

рисунке

3.5

изобразим

множители выражения

 

S

д

 

,

 

 

 

отвечающие

за

модуль

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

спектральной плотности.

34

3

2 1

0

2cos 2 Tд 1

 

 

 

 

 

1

T

 

cos

 

 

2

д

Рисунок 3.5 Представление множителей, входящих в выражение Sд

N

На рисунке 3.6 изображен модуль спектральной плотности дискретного непериодического сигнала Sд .

Sд ( )

0.6

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

N

0.4

 

 

 

 

S ()

0.3

 

 

 

 

0.2

 

 

 

 

Т ос

 

 

 

 

0.1

 

 

 

 

- д

2

 

2

 

д

 

 

 

 

 

 

Рисунок 3.6 – Модули спектральных плотностей исходного аналогового (пунктирная линия) и дискретизированного (сплошная линия) сигналов

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

3.4Расчет коэффициентов ДПФ

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

35

отсчетами его спектральной плотности с помощью прямого дискретного преобразования Фурье (ПДПФ):

 

 

N 1

j

2 nk

 

 

 

 

 

 

1

s nTд e

N

 

 

 

C k

 

 

 

.

 

N n 0

 

 

 

Расчет комплексных коэффициентов ряда Фурье Ck с использованием

ЦВМ не составляет труда. Однако необходимо представлять себе, каким образом происходит векторное суммирование компонент.

На рисунках 3.7 – 3.9 изображены векторные диаграммы, отражающие расчет комплексных коэффициентов ряда Фурье. Рисунку 3.7 соответствует аналитическое выражение

 

E

 

 

j

2

j

4

j

6

j

8

j

10

 

 

 

С

 

e

12

12

12

12

12

0 0 0 0 0 0

 

 

 

1

 

 

e

 

e

 

e

 

e

 

 

1

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е 0,083 j0,311 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90

 

 

 

 

 

 

 

 

 

 

 

 

 

 

120

 

 

 

60

 

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

180

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0.05

 

0.1

 

 

210 330

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

240

 

 

 

 

 

300

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

270

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 3.7 – Векторная диаграмма определения коэффициента C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунку 3.8 соответствует аналитическое выражение

 

 

 

 

 

 

 

E

 

j

4

j

8

j

12

j

16

 

j

20

 

 

 

 

 

 

 

 

 

12

12

12

12

 

 

12

 

 

 

 

 

С2

 

12

1

e

e

e

 

e

 

 

e

 

 

0 0

0 0

0 0

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

36

Рисунок 3.8 – Векторная диаграмма определения коэффициента C2

 

 

E

 

j

6

j

12

j

18

j

24

j

30

 

 

 

 

 

 

12

12

12

12

12

 

 

 

С

 

 

1

e

 

e

 

e

 

e

 

e

 

0 0 0 0 0 0

 

 

3

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е 0,083 j0,083 .

Рисунок 3.9 – Векторная диаграмма определения коэффициента C3

В таблице 3.1 приведены все коэффициенты ряда Фурье.

37

Таблица 3.1 – Комплексные коэффициенты ряда Фурье Ck

 

 

 

 

 

C0

0,5Е

C6

0

C1

Е(0,083-j0,311)

C7

Е(0,083+j0,022)

 

 

 

 

C2

0

C8

0

C3

Е(0,083-j0,083)

C9

Е(0,083+j0,083)

C4

0

C10

0

 

 

 

 

C5

Е(0,083-j0,022)

C11

Е(0,083+j0,311)

 

 

 

 

На рисунке 3.10 представлены графически модули комплексных коэффициентов Фурье

. 0,6 Сk 0,5

Е0,4

0,3

0,2

0,1

k ·2π

Tоc

Рисунок 3.10 – Отсчеты спектральной плотности, полученные по ДПФ

Дискретный периодический сигнал обладает дискретным периодическим спектром. Отсчеты во временной и частотной областях связаны парой ДПФ. При этом число степеней свободы во временной и частотной областях имеет одно и то же значение N:

N

в

2 T

.

 

оc

3.5Восстановление исходного сигнала по ДПФ

Если на основании совокупности отсчетов s0, s1, s2, ..., sN-1 некоторого сигнала найдены коэффициенты ДПФ С012,...,СN-1, то по ним всегда можно восстановить исходный сигнал s(t) с ограниченным спектром, который был подвергнут дискретизации. Ряд Фурье такого сигнала (для четного числа N) принимает вид конечной суммы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s t C0

 

 

 

 

2

 

 

 

4

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

C1

 

cos

T

t 1

2

C2

 

cos

T

t 2

...

C N

 

cos

T

t

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оc

 

 

 

оc

 

2

 

оc

 

2

38

В данном случае i argCi - фазовый угол коэффициента ДПФ.

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

На рисунке 3.11 представлен результат восстановления сигнала по отсчетам его спектральной плотности.

s(t)

1,2

Е

0,8

0,6

0,4

0,2

0,2

t

Рисунок 3.11 – Аналоговый периодический сигнал, восстановленный по коэффициентам ДПФ

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

3.6Z-преобразование дискретной последовательности

Прямое Z-преобразование последовательности s nTд определяется формулой

S z s nTд z n .

n 0

Функцию S z называют Z-образом последовательности s nTд .

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

Найдем Z-образ дискретной функции Хэвисайда (n):

Z (n) (z) ,

 

 

 

 

 

 

 

 

 

39

 

 

 

 

 

 

 

 

 

 

 

(z) 1 z 1

z 2 z 3 ... z n...

 

 

 

 

1 (z) z 1 z 2 z 3 z 4 ... z n 1... .

 

 

 

 

 

z

 

 

 

 

 

 

(1 z 1) (z) 1 z (n 1)

 

Если

 

z

 

1, то

lim z (n 1)

0 , откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

(z)

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 z 1

Сдвиг дискретной последовательности приводит к умножению Z-образа

на оператор сдвига z 6 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z (n 6) (z)z6 .

Таким образом, Z-образ дискретизированного прямоугольного

импульса равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S(z) E (z) (z)z 6 E

1 z 6

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 z 1

Преобразуем Z-образ к виду, удобному для получения спектральной

плотности дискретного сигнала:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S z Е

 

z3

z 3

z 5 2 .

 

 

 

 

 

 

z

1 2

 

1 2

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

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

ze j Tд .

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

Sд Е sin 3 Tд e j5 Tд 2 . sin Tд 2

40

На рисунке 3.12 изображен модуль спектральной плотности дискретизированного сигнала, полученной на основе Z-преобразования (с учетом нормирования относительно числа степеней свободы N =12).

Таким образом, пара Z-преобразований позволяет связать частотный и временной образы дискретного сигнала. Причем выборке отсчетов сигнала во временной области соответствует периодическая спектральная плотность в частотной области с периодом повторения, равным д.

Sд ( )

0,6

 

 

 

0,5

 

 

 

 

 

 

 

N

0,4

 

 

 

 

0 ,3

 

 

 

 

0,2

 

 

 

 

0,1

 

 

 

- д

2 0

2

д

 

Рисунок 3.12 – Спектральная плотность последовательности отсчетов

сучетом нормирования относительно N

3.7Восстановление аналогового сигнала использованием ряда Котельникова

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

 

 

sin Tд (t nTд )

 

s(t)

s(nTд )

.

Tд (t nTд )

 

 

n

 

Другими словами, восстановленный сигнал представляет собой сумму функций Котельникова с весами, равными отсчетам сигнала. На рисунке 3.13 приведен результат восстановления.

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