Импульсная характеристика дискретной цепи

Обычно нас интересуют линейные дискретные цепи. Особенностью этих дискретных цепей является то, что на выходе не образуется новых частотных составляющих сигнала. Если дискретная цепь является фильтром, то соотношение частотных компонент спектра входного сигнала на выходе изменяется. Это удается сделать при помощи дискретной свертки сигнала. Формула дискретной свертки входной последовательности отсчетов с импульсной характеристикой дискретной цепи (фильтра) записывается следующим образом:

(1)

где h(k) — импульсная характеристика дискретной цепи. Импульсную характеристику можно определить как отклик дискретной цепи на воздействие единичного импульса (δ-функция). Формулу свертки (1) можно записать в сокращенном виде, используя символ операции свертки '*':

y(k) = x(k)*h(k) (2).

Связь переходной характеристики h(t) с операторной передаточной функцией Н(р).

По определению, операторная передаточная функция

(2)

где - изображение воздействия цепи ( или );

- изображение реакции цепи ( или ).

Из (2) следует, что

. (3)

Если в качестве воздействия на цепь подать единичную ступенчатую функцию , то ее изображение будет . При этом реакцией цепи будет переходная характеристика цепи . Изображение реакции цепи можно найти, воспользовавшись выражением (3):

(4)

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

Быстрым преобразованием Фурье (БПФ) называют набор алгоритмов, реализация которых приводит к существенному уменьшению вычислительной сложности ДПФ. Основная идея БПФ состоит в том, чтобы разбить исходный N-отсчетный сигнал x(n) на два более коротких сигнала, ДПФ которых могут быть скомбинированы таким образом, чтобы получить ДПФ исходного N-отсчетного сигнала.

Так, если исходный N-отсчетный сигнал разбить на два N/2-отсчетных сигнала, то для вычисления ДПФ каждого из них потребуется около (N/2)2 комплексных умножений. Тогда для вычисления искомого N-отсчетного ДПФ потребуется порядка 2(N/2)2=N2/2 комплексных умножений , т.е. вдвое меньше по сравнению с прямым вычислением. Операцию разбиения можно повторить, вычисляя вместо (N/2)-отсчетного ДПФ два (N/4)-отсчетных ДПФ и сокращая тем сасым объем вычислений еще в два раза. Выигрыш в два раза является приблизительным, поскольку не учитывается, каким образом из ДПФ меньшего размера образуется искомое N-отсчетное ДПФ.

Определение 1 .

Дана конечная последовательность x₀ , x₁ , x₂ ,...,x_N-1 (в общем случае комплексных). Дискретное преобразование Фурье (ДПФ) заключается в поиске другой последовательности X₀ , X₁ , X₂ ,...,X_N-1 элементы которой вычисляются по формуле:

(1).

Определение 2 .

Дана конечная последовательность X₀ , X₁ , X₂ ,...,X_N-1 (в общем случае комплексных). Обратное дискретное преобразование Фурье (ДПФ) заключается в поиске другой последовательности x₀ , x₁ , x₂ ,...,x_N-1 элементы которой вычисляются по формуле:

(2).

Основным свойством этих преобразований (которое доказывается в соответствующих разделах математики) является тот факт, что из последовательности {x} получается (при прямом преобразовании) последовательность {X}, а если потом применить к {X} обратное преобразование, то снова получится исходная последовательность {x}.

Определение 3 .

Величина

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

Рассмотрим ряд свойств поворачивающих множителей, которые нам понадобятся в дальнейшем.

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

Прямое преобразование Фурье можно выразить через поворачивающие множители. В результате формула (1) примет вид:

(3).

Эти коэффициенты действительно оправдывают свое название. Нарисуем на комплексной плоскости любое комплексное число, в виде вектора, исходящего из начала координат. Представим это комплексное число в показательной форме: re^jφ , где r - модуль числа, а φ - аргумент. Модуль соответствует длине вектора, а аргумент - углу поворота:

Теперь возьмем какой-нибудь поворачивающий множитель

. Его модуль равен единице, а фаза - 2π/N. Как известно, при умножении комплексных чисел, представленных в показательной форме, их модули перемножаются, а аргументы суммируются. Тогда умножение исходного числа на поворачивающий множитель не изменит длину вектора, но изменит его угол. То есть, произойдет поворот вектора на угол 2π/N (см. предыдущий рисунок).

Если теперь посмотреть на формулу (3), то станет ясен геометрический смысл преобразования Фурье: он состоит в том, чтобы представить N комплексных чисел-векторов из набора {x}, каждое в виде суммы векторов из набора {X}, повернутых на углы, кратные 2π/N.

Теорема 0 .

Если комплексное число представлено в виде e ^j2πN , где N - целое, то это число e ^j2πN = 1.

Доказательство :

По формуле Эйлера, и ввиду периодичности синуса и косинуса:

e ^{j2 π N} = cos(2πN) + j sin(2πN) = cos 0 + j sin 0 = 1 + j0 = 1

Цифровые фильтры

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

Если определена функция передачи F(z) цифрового фильтра, свой­ства которого повторяют свойства аналогового фильтра-прототипа с функцией передачи F(p), то по ней просто строится каноническая схема, которую можно рассматривать как алгоритм обработки сигнала. Цифровую фильтрацию сигналов осуществляют на основе выполнения операций только трех типов: задержки, сложения и умножения. По­этому алгоритм цифровой обработки может быть реализован двумя способами: универсальной ЭВМ, выполняющей цифровую обработку по специальной программе, или специализированным вычислительным устройством, выполняющим только три указанные выше операции. Первый способ реализации алгоритма цифровой фильтрации называют программным, второй — аппаратурным.

Программный способ реализации эффективен при модели­ровании различных систем цифровой обработки сигналов, так как поз­воляет легко изменять алгоритм фильтрации. Цифровые фильтры, предназначенные для работы в системах автоматики, телемеханики и связи, должны обрабатывать сигналы в реальном масштабе времени, т. е. за время, не большее периода дискретизации входных сигналов, что является их особенностью. Использование универсальных ЭВМ для этих целей практически невозможно из-за их сравнительно низ­кого быстродействия и значительной стоимости. Аппаратур­ная реализация цифровых фильтров основана на использовании циф­ровых интегральных схем, представляющих собой регистры сдвига, сумматоры, умножители и т. п. Она стала возможной в связи с появле­нием в последние годы больших интегральных схем, имеющих большие функциональные возможности и высокое быстродействие.

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