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

Ris_DSP / TZOS_red4_ch1_for_Tablet

.pdf
Скачиваний:
142
Добавлен:
18.04.2015
Размер:
1.97 Mб
Скачать

 

 

Цифр.

 

 

4q

сигнал

 

 

3q

 

h(t)

а)

2q

 

 

 

q

 

x(t)

 

- 4q -2q

q 2q

Дискр.

 

 

- 2q

сигнал

 

 

- 3q

 

 

 

- 4q

 

 

 

e(n)

 

б)

q/2

 

n

 

-q/2

11

а)

б)

в)

Рис.1.6. Механизм образования шума квантования. а) характеристика преобразования идеального АЦП, б) шум

квантования

Рис.1.7. Шума квантования гармонического сигнала а) аналоговый сигнал, б) квантованный сигнал, в) шум квантования

Для примера на рис.1.7 показан шум квантования, образующийся после дискретизации гармонического сигнала.

Кодер ставит в соответствие каждому уровню входного сигнала двоичный В- разрядный код (рис.1.8).

Рис.1.8. Двоичное кодирование цифрового сигнала

1.3. КОДИРОВАНИЕ ЧИСЕЛ В УСТРОЙСТВАХ ЦИФРОВОЙ ОБРАБОТКИ СИГНАЛОВ

Для кодирования чисел в устройствах ЦОС используются три основных двоичных кода: прямой, обратный, дополнительный.

12

Код числа содержит В разрядов. Старший разряд служит для фиксации знака числа и называется знаковым разрядом. Остальные разряды - числовые.

Формат представления чисел может быть целочисленный или дробный. Различия этих форматов заключается в форме представления числовых разрядов (мантиссы).

В целочисленном формате значащие числовые биты выравниваются по правому краю, а «лишние» старшие разряды заполняются значением знакового разряда (происходит расширение знака). Расширение знака характерно для чисел двойной длины, представленных в дополнительном коде.

Пример представления целых 8-разрядных чисел со знаком в прямом коде.

а) положительное число

7

6

5

4

3

2

1

0

0

0

1

0

0

1

1

1

знак

26

25

24

23

22

21

20

десятичный эквивалент: + (25+22+21+20) = 39

а) отрицательное число

7

6

5

4

3

2

1

0

1

0

1

0

0

1

1

1

знак

26

25

24

23

22

21

20

десятичный эквивалент: - (25+22+21+20) = - 39

В дробном формате значащие числовые биты выравниваются по левому краю, а «лишние» младшие разряды обнуляются. После знакового бита фиксируется условная точка, отделяющая целую часть (равную нулю) от дробной.

Пример представления дробного 8-разрядного числа со знаком в прямом коде.

7

6

5

4

3

2

1

0

0

0

1

0

0

1

1

1

знак

2-1

2-2

2-3

2-4

2-5

2-6

2-7

десятичный эквивалент: + (2-2+2-5+2-6+2-7) = 0.3046875

13

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

Далее рассмотрим целочисленный формат. Прямой код

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

Aпр =

0.

a

B 2

... a a

;

A 0

 

 

1 0

;

A 0

 

1.

aB 2

... a1a0

Пример: 0.10100 = + 2010 1.1101 = - 1310

Обратный код Правило кодирования: для положительных чисел совпадает с прямым

кодом Аобр = Апр; для отрицательных чисел в знаковый разряд кода записывается «1», а числовые разряды инвертируются:

 

0.

a

B 2

... a a

0

;

A 0

Аобр =

 

 

 

1

 

 

aB 2 ...a1a0 ;

A 0

 

1.

Пример: 0.10100 = + 2010 1.1101 = - 210

Дополнительный код Правило кодирования: для положительных чисел совпадает с прямым

кодом Адоп = Апр; для отрицательных чисел в знаковый разряд записывается «1», числовые разряды инвертируются и к младшему разряду добавляется 1:

Адоп =

0.

aB 2

... a1a0

;

A 0

1.

 

 

...

 

 

 

20

;

A 0

a

B 2

a

a

 

 

 

1 0

 

 

 

Пример:

Представить число (-25) в дополнительном коде.

прямой код дополнительный код

 

 

14

1.11001

 

+ 1.00110

1

1.00111

- 2510 = 1.00111

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

 

Смещенный код

Дополнительный код

+Вся шкала

1.1111111

0.1111111

+Вся шкала - 1

1.1111110

0.1111110

 

 

 

0+1 мл. разряда

1.0000001

0.0000001

0

1.0000000

0.0000000

0- 1 мл. разряда

0.1111111

1.1111111

 

 

 

- Вся шкала +1

0.0000001

1.0000001

- Вся шкала

0.0000000

1.0000000

1.4. ОГРАНИЧЕНИЕ РАЗРЯДНОСТИ ОБРАБАТЫВАЕМЫХ ДАННЫХ

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

Округление

Для округления числа x(n) до B разрядов исходное М-разрядное число (В<M< ) заменяется на ближайшее В-разрядное xорк(n).

Т.е. округление соответствует выбору ближайшего уровня квантования. При округлении возникает ошибка, равная

eорк (n) = xорк(n) - x(n) , причем - q/2 eорк (n) q/2 .

Ошибка округления распределена равномерно, математическое ожидание (МО) ошибки округления равно нулю (рис.1.9).

15

Рис.1.9. Распределение ошибки округления

B A

Пример: x(n)=0.00111110 , после округления xокр(n)=0.0100. M

Усечение При усечении М-разрядного числа до В разрядов младшие (М-В)

разрядов исходного числа отбрасываются. При усечении также возникает ошибка:

eус(n) = xус(n) - x(n) , причем - q eус(n) 0 .

При усечении разрядности большого массива данных ошибка усечения накапливается, математическое ожидание смещается влево (рис.1.10).

Рис.1.10. Распределение ошибки усечения

B A

Пример: x(n)=0.00111110 , после усечения xус(n)=0.0011. M

2. АЛГОРИТМЫ ЦИФРОВОЙ ОБРАБОТКИ СИГНАЛОВ

Среди многообразия алгоритмов цифровой обработки сигналов основными являются:

цифровая фильтрация;

корреляционная обработка;

согласованная фильтрация;

спектральный анализ;

адаптивная обработка.

16

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

2.1. СВЕРТКА

Известно, что линейная стационарная система (рис.2.1) преобразует непрерывный входной сигнал x(t) в соответствии с интегралом Дюамеля таким образом, что на ее выходе возникает колебание y(t), равное

y(t) x( ) h(t )d (2.1)

0

Рис.2.1. Линейная система Цифровая же система (рис.2.2), представляющая собой дискретную

систему, преобразует последовательность дискретных входных отсчетов x(n) в последовательность дискретных выходных отсчетов y(n) в соответствии с выражением свертка:

y(n) x(m) h(n m) (2.2)

m 0

Рис.2.2. Дискретная система

или

y(n) x(n m) h(m), (2.3)

m 0

где h(t), h(m) - импульсные характеристики, соответственно аналоговой и дискретной систем.

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

Наиболее эффективные алгоритмы разработаны для вычисления циклической свертки. Она определяется соотношением:

N 1

 

y(n) x(n m) h(m),

(2.4)

m 0

где x(n), h(m) - конечные последовательности, n 0,N 1;m 0,N 1. Свертку можно вычислять прямыми методами и на основе

преобразования Фурье (дискретного преобразования Фурье - ДПФ).

17

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

Таким образом, свертка двух дискретных последовательностей x(n) и h(m) может быть вычислена в результате выполнения следующих операций:

1. Вычисление ДПФ последовательности x(n):

N 1

j

2

nk

 

 

X(k) x(n) e N

 

,n 0,N 1;k 0,N 1.

n0

2.Вычисление ДПФ последовательности h(m):

N 1

j

2

mk

 

 

H(k) h(m) e N

 

,m 0,N 1;k 0,N 1.

m0

3.Вычисление произведения результатов ДПФ:

Y(k) X(k) H(k)

4. Вычисление обратного ДПФ последовательности Y(k):

 

1

N 1

j

2

nk

 

 

 

 

 

 

 

 

 

 

 

y(n)

Y(k) e N

,n 0,N 1;k 0,N 1,

N

 

k 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если для вычисления свертки используется алгоритм быстрого преобразования Фурье (БПФ), повышающий скорость вычисления преобразования Фурье при больших массивах данных на несколько порядков, то этот метод вычисления свертки называют высокоскоростным. На рис.2.3 представлена структурная схема вычислителя высокоскоростной свертки.

Рис.2.3. Структурная схема вычисления высокоскоростной свертки

18

2.2.ЦИФРОВЫЕ ФИЛЬТРЫ

Сфизической точки зрения цифровая фильтрация - это выделение в определенном частотном диапазоне с помощью цифровых методов полезного сигнала на фоне мешающих помех (рис.2.4).

По своим частотным свойствам фильтры делятся на: фильтры нижних частот (ФНЧ) – Low pass – рис.2.5а; фильтры верхних частот (ФВЧ) – High pass – рис.2.5б; полосовые фильтры (ПФ) – Band pass – рис.2.5в; режекторные фильтры (РФ) – Band stop – рис.2.5г.

Идеальные частотные характеристики цифровых фильтров представлены на рис.2.5, соответственно.

На рис.2.5 приняты следующие обозначения:

ПП – полоса пропускания – частотная область, внутри которой сигналы проходят через фильтр практически без затухания;

ПЗ – полоса задерживания – выбирается разработчиком такой, чтобы обеспечить затухание сигнала не хуже заданного;

Переходная полоса – частотная область между ПП и ПЗ (характеризуется скоростью спада, обычно измеряется в дБ/декаду);

fп - частота среза полосы пропускания – определяется допустимым уровнем пульсаций ЧХ в полосе пропускания;

fз - частота среза полосы задерживания – определяется уровнем пульсаций ЧХ в ПЗ;

fнп , fвп – нижняя и верхняя частоты среза полосы пропускания; fнз , fвз – нижняя и верхняя частоты среза полосы задерживания.

Для ЦФ частоты среза являются условными границами (не на уровне - 3дБ!), располагаются между частотой среза полосы пропускания и частотой среза полосы задерживания и определяются границами нестабильности ЧХ в ПП и ПЗ.

АЧХ реальных фильтров (рис.2.6, на примере ФНЧ) имеют пульсации в полосе пропускания п и задерживания з (нестабильность ЧХ в ПП и ПЗ). Часто в литературе они имеют другое название:

Rз – минимальное подавление в полосе задерживания, дБ; Rп – максимальное подавление в полосе пропускания, дБ.

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

19

Рис.2.4. Фильтрация мешающих помех с помощью цифрового ПФ

Рис.2.5. Идеальные частотные характеристики фильтров

20

Рис.2.6 Реальная АЧХ цифрового фильтра (на примере ФНЧ)

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

M 1 N 1

y(n t) b( j) y(n t j t) a(i) x(n t i t), (2.5)

j 1 i 0

где x(n t) и y(n t) - n- тые отсчеты входного и выходного сигналов фильтра, взятые через интервал t ( t= td); a(i) и b(j) - постоянные коэффициенты цифрового фильтра.

Цифровые фильтры принято делить на два класса:

нерекурсивные фильтры;

рекурсивные фильтры.

Нерекурсивные фильтры называют еще фильтрами с конечной импульсной характеристикой (КИХ-фильтры), а рекурсивные фильтры - фильтрами с бесконечной импульсной характеристикой (БИХ-фильтры). В иностранной литературе их называют:

FIR (Finite Impulse Response) – фильтр с конечной импульсной характеристикой;

IIR (Infinite Impulse Response) – фильтр с бесконечной импульсной характеристикой.

Если в выражении (2.5) положить коэффициенты b(j)=0, то фильтр, реализующий этот алгоритм, называется нерекурсивным. Его работа описывается уравнением:

N 1

 

y(n t) a(i) x(n t i t)

(2.6)

i0

ипредставляет собой свертку двух последовательностей: коэффициентов

a(i) и дискретных отсчетов входного сигнала x(n t).

Соседние файлы в папке Ris_DSP