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

Метода ЦОС

.pdf
Скачиваний:
160
Добавлен:
24.02.2016
Размер:
2.7 Mб
Скачать

.

Замечание. Если выполняется программная реализация КИХ фильтра, входные отсчеты должны быть предварительно записаны в ЗУ. Далее необходимо выбирать блоки из 53 отсчетов из массива входных данных , выполнять вычисления и запоминать результирующую выходную последовательность в ЗУ.

6. Ди етное п еоб зов ние

ье

6.1. Конечные дискретные комплексные экспоненциальные функции

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

(6.1)

Данная система функций составляет счетное бесконечное множество и определена на бесконечном интервале частоты

Экспоненциальная последовательность может быть задана на конечном

интервале

времени

, где

целое положительное

число. Тогда

величина

определяет

основной

период дискретной

комплексной

экспоненциальной последовательности. В этом случае значение

 

– основная линейная частота последовательности. Основная круговая частота

определяет период дискретизации по частоте. Абсолютное значение непрерывной частоты с учетом дискретизации переменной частоты равна

Далее преобразуем (6.1) следующим образом:

 

 

 

 

 

 

.

(6.2)

 

 

 

 

 

 

 

 

 

 

 

 

В дискретном

преобразовании Фурье используется

система

комплексных дискретных экспоненциальные функции (ДЭФ), определяемых выражением

Введем обозначение . Тогда

 

 

 

 

 

.

 

Переменная

называется поворачивающим множителем. Переменные

и

принимают

целочисленные

значения

Так

как

показатель степени комплексного числа

со знаком „плюс“, то функция

 

 

 

описывает точку,

которая

движется по окружности

в

 

 

 

направлении часовой стрелки. В выражении (6.2) переменные времени и частоты изменяются дискретно, в отличие от (6.1), где время изменяется дискретно, а частота непрерывно. Заметим, что огибающая дискретных

значений функции

 

соответствует функции

Рисунок 6.1

 

иллюстрирует графическое представление этой функции

 

.

 

Рисунок 6.1– Графическая иллюстрация функции

(формулы Эйлера)

Если переменная последовательно принимает значения то через шагов, комплексный вектор проходит радиан или совершает

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

равен

 

а его аргумент

 

 

 

 

 

 

 

 

 

 

 

Пример 6.1. Пусть

Значения фазы

вектора ДЭФ для

 

 

соответственно

равны

 

 

Следовательно, при увеличении

фаза ДЭФ нарастает по линейному закону.

Пример 6.2. Пусть

 

Значения фазы вектора ДЭФ для

 

 

соответственно равны

 

 

Через 8 шагов комплексный вектор проходит

радиан или совершает два

оборота на комплексной плоскости за тоже время , что и в примере 6.1, т.е. за

где – интервал дискретизации.

 

 

 

 

Скорость нарастания фазы

вектора ДЭФ определяет номер

Можно

сказать, что фаза ДЭФ нарастает со скоростью

 

радиан. В примере 6.2 со

 

скоростью

н в е енной от чет

 

Величина полной фазы за дискретное время определяется как

где – скорость изменения фазы ДЭФ или частота этой функции. Таким

образом, частота ДЭФ – это число оборотов, совершаемых вектором ДЭФ на интервале ее определения

Пример 6.3. Вычисление значений ДЭФ.

1.

Решение.

2.

Решение.

3.

Решение.

4.

Решение.

На рисунке 6.2 показаны положения вектора ДЭФ на комплексной плоскости примера 6.3.

Рисунок 6.2 – Положение вектора ДЭФ

 

Систему ДЭФ записывают в виде матрицы

строки которой

нумеруются переменной столбцы переменной . В пересечении k-й строки и -го столбца записывается величина :

 

.

 

Например, для

матрица имеет следующий вид:

 

 

.

(6.3)

Если подставить в эту матрицу числовые значения степенного ряда

, то

получим

 

 

 

.

(6.4)

На рисунке 6.3 показаны положения вектора ДЭФ и ее значения на комплексной плоскости, соответствующие матрице (6.4).

Рисунок 6.3 – Функции ДЭФ для

6.1.1. Свойства дискретных экспоненциальных функций

1. Функции

ортогональны, т.е.

 

е ли

 

е ли

Так как

, то

 

е ли

 

е ли

Следствием свойства ортогональности является:

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

(6.5)

, одна из

– скалярное произведение одинаковых двух строк матрицы

, одна из

которых должна быть комплексно сопряженной, равно .

 

Действительно,

 

 

 

Сумма единиц даст число

 

 

 

Матричная запись свойства ортогональности имеет вид

 

 

 

,

 

 

(6.6)

где

единичная матрица.

 

 

 

 

 

 

2. Периодичность:

 

 

 

 

 

если

то

 

 

 

 

 

 

 

.

(6.7)

Поскольку ДЭФ являются периодическими функциями, матрицу (6.3) можно

переписать с минимальными фазами

, образующимися после

вычитания из значения целого числа периодов

т.е.

 

.

 

Для

матрица ДЭФ (6.3) с минимальными фазами

3. Симметричность.

ДЭФ является функцией двух переменных и Выводы относительно одной из переменных справедливы и для другой. Тогда

4. Обратная матрица ДЭФ.

 

Из свойства ортогональности

. Умножим обе части этого

равенства слева на

 

.

(6.8)

5. Мультипликативность:

по строкам

по столбцам

 

 

 

 

 

Действительно,

 

 

 

 

 

. При умножении любых двух

строк (столбцов) матрицы

получается строка (столбец) той же матрицы.

Номер полученной строки (столбца) равен сумме номеров сомножителей.

6.2. Определение дискретного преобразования Фурье

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

области (экспоненциальная форма)

(6.9)

где – индекс ДПФ в частотной области, – индекс временной входной последовательности отсчетов сигнала.

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

Обратное ДПФ (ОДПФ) имеет следующий вид:

(6.10)

Взаимная обратимость выражений (6.9) и (6.10) доказывается подстановкой

вт.е.

(6.11)

Так как

не зависит от , изменяем порядок суммирования в (6.11),

(6.12)

Всилу ортогональности ДЭФ внутренняя сумма отлична от нуля только при

Вэтом случае правая часть выражения (4.24) равна

Тригонометрическая форма ДПФ:

прямое ДПФ

обратное ДПФ

Замечание. Принципиальное различие между дискретизированным по времени преобразованием Фурье и ДПФ обусловлено характером системы

функций

и {

}, а именно:

 

 

– огибающая дискретных значений функции

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

функции

 

 

 

 

– конечный интервал времени

задания функции

;

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

6.3. Свойства дискретного преобразования Фурье

1. Периодичность. Свойство периодичности ДЭФ приводит к выражениям

и

Действительно, Обычно ограничиваются рассмотрением одного периода длиной во

временной и в частотной области. Это позволяет определить матричную форму ДПФ:

– прямое ДПФ

где

 

 

 

 

и

 

– векторы

отсчетов

последовательности

спектральных

коэффициентов

и сигнала

соответственно;

 

 

 

 

 

 

 

– обратное ДПФ

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

Используя формулу (4.20), получаем

 

 

 

 

 

 

 

 

.

 

 

 

 

2. Линейность. Класс линейных систем определяется линейными

операциями или принципом

 

суперпозиции.

Если

и

входные

последовательности, а

 

и

соответственно их ДПФ, то при подаче

на

вход

последовательности

 

 

систему

называют

линейной тогда и только тогда, когда выполняется

 

 

где

и

произвольные

постоянные параметры

(константы). Спектр

последовательности

равен

 

 

 

3. Инвариантность ДПФ относительно сдвига по времени и частоте:

1. Инвариантность относительно циклического сдвига по времени. Если последовательность имеет ДПФ , то ДПФ последовательности

равно

Рассмотрим две последовательности

и

. Формы

последовательностей показаны на рисунке 6.4а,б.

 

а)

x(n)

n

x(n-h)

б)

 

h

 

n

 

 

 

 

 

 

 

 

Рисунок 6.4

 

 

 

 

ДПФ последовательности

 

равно

 

 

 

 

.

 

Заменяя индекс суммирования

и

введя новую переменную

,

получим

 

 

 

 

где

 

 

 

 

 

 

 

). Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, при сдвиге дискретного сигнала по времени изменениям подвергаются только фазы дискретных функций (фазовый спектр), амплитудный спектр не изменяется.

2. Инвариантность относительно сдвига по частоте. Если спектральной

последовательности

соответствует последовательность

то при

сдвиге последовательности

исходная последовательность

получит фазовый сдвиг, т.е.

 

 

Пусть

Обратное ДПФ последовательности

равно

 

.

 

Заменяя индекс

суммирования , и введя новую переменную

,

получим

 

 

где

 

 

 

 

 

).

 

 

 

4. Теорема о свертке. Если исходные последовательности отсчетов сигналов и имеют конечные периоды , их циклическая свертка определяется формулой

, = 0, 1,…, –1.

Вычислим ДПФ последовательности

(6.13)

Так как

не зависит от , изменяем порядок суммирования в (6.13).

. (6.14)

Используя свойство инвариантности относительно циклического сдвига по времени, можно записать составляющую выражения (6.14) как

Тогда

(6.15)

Таким образом, спектр свертки равен произведению спектров сворачиваемых последовательностей. Коэффициенты свертки вычисляются на основе ОДПФ по формуле

Теорема (6.15) позволяет вычислить коэффициенты свертки при помощи ДПФ по формуле

Д

Д

Д

.

При больших величинах на практике применяют эффективные алгоритмы вычисления свертки с использованием быстрых преобразований Фурье.

5. Теорема о корреляции. По определению (2.13)

корреляционная

функция двух конечных последовательностей

и

равна

 

, для

= 0, 1,…,

–1.

Вычислим ДПФ последовательности

 

 

 

 

 

 

 

(6.16)

Так как

не зависит от , изменяем порядок суммирования в (6.16).

 

 

 

.

(6.17)

Используя свойство инвариантности относительно циклического сдвига по времени, можно записать составляющую выражения (6.17) как

Тогда

(6.18)

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

Коэффициенты корреляционной функции вычисляются на основе ОДПФ по формуле

Теорема (6.18) позволяет вычислить коэффициенты корреляционной функции при помощи ДПФ по формуле

Д Д Д .

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

6. Теорема Парсеваля. Пусть последовательности

и

будут

идентичными. В этом случае теорема о корреляции записывается как

 

 

.

 

Коэффициенты корреляционной функции, вычисляются на основе выражения ОДПФ, т.е.