13.2. Алгоритм дискретной свертки (во временной области)
Схема устройства дискретной обработки континуального сигнала (без цифрового кодирования) представлена на рис. 13.2.
Рис. 13.2. Схема дискретной обработки сигнала
Процедуру дискретизации входного сигнала (взятие выборок) удобно рассматривать как умножение функции s(t) на периодическую последовательность yT(t) тактовых импульсов. В качестве таких импульсов обычно рассматривают дельта-функции, так что функцию yT(t) можно определить выражением
График этой функции показан на рис. 13.3, б (цифра единица обозначает площадь тактового импульса).
Рис. 13.3. Дискретизация сигнала: а - исходный континуальный сигнал; б - последовательность тактовых импульсов; в - дискретные отсчеты сигнала
Тогда дискретизированный сигнал
где
При таком обозначении sk имеет смысл бесконечно короткого импульса, расположенного на оси времени в точке t = kT и обладающего площадью, численно равной выборке s(kT). Таким образом, выборки s(kT) из сигнала s(t) (рис. 13.3, а), являются весовыми коэффициентами дельта-функции (рис. 13.3, в)*.
* (Следует иметь в виду, что по своей размерности sT(t) отличается от s(t), поскольку функция δ(t) имеет размерность 1/t. Для сохранения размерности выражение (13.3) следовало бы записать в форме sk = Ts(kT)×δ(t - kT), что равносильно приравнению площади импульса sk площади прямоугольного импульса с основанием Т и амплитудой s(kТ) [4]. )
Дискретизированный сигнал ST(t) подается на дискретный фильтр, влияние которого на ST(t) должно быть эквивалентно влиянию некоторого аналогового фильтра на s(t). Пусть импульсная характеристика последнего g(t) задана. Тогда сигнал на выходе аналогового фильтра определяется сверткой (см. выражение (6.12)]:
Дискретный эквивалент свертки при надлежащем выборе шага Т можно записать в форме
где sвых(nТ) - n-я выборка из sвых(t), а соответствующий этой выборке импульс согласно (13.3)
В выражении (13.5) g(kT) имеет смысл k-й выборки из импульсной характеристики аналогового фильтра, показанной на рис. 13.4 штриховой линией.
Рис. 13.4. Дискретизация импульсной характеристики аналогового фильтра
Выражение (13.5) можно рассматривать как алгоритм дискретной обработки, эквивалентной пропусканию сигнала через аналоговый фильтр с заданной импульсной характеристикой g(t).
При отсчете времени от начала входного сигнала выражение (13.5) принимает вид
Верхний предел суммирования следует из условия, что при k > n s[(n - k)T] = 0, а нижний - из условия, что при t < 0 g(t) = 0.
Схема устройства, обладающего требуемой дискретной импульсной характеристикой gT(t), представлена на рис. 13.5, На этой схеме буквами а0, а1, а2, ..., аН обозначены не зависящие от частоты коэффициенты усиления, а T - идеальная линия задержки; задержка совпадает с темпом поступления выборок сигнала s(t). Подбором коэффициентов а0 = g(0), а1 = g(T), а2 = g(2T), ... можно, в принципе, осуществить весьма сложные импульсные характеристики gT(t).
Очевидно, что импульсная характеристика рассматриваемого устройства должна записываться в форме
Таким образом, ak = g(kT) является лишь весовым коэффициентом при дельта-функции δ(t - kT). При подаче на вход схемы, обведенной на рис. 13.5 штриховой линией, дискретного сигнала sT(t) на выходе сумматора ∑ получается дискретная последовательность импульсов sT вых(t), а на выходе синтезирующего фильтра - профильтрованный непрерывный сигнал sвых(t).
Рис. 13.5. Дискретный фильтр
Осуществить дискретный фильтр непосредственно по схеме на рис. 13.5 можно лишь для относительно простых сигналов, обладающих небольшой базой N = 2fmTc. (Тс = NT - длительность обрабатываемой реализации сигнала.) Однако алгоритм (13.5) широко используется при моделировании фильтров с помощью ЭВМ.
Основная трудность реализации дискретного фильтра заключается также в осуществлении элемента памяти Т. Эта трудность отпадает при переходе к цифровой фильтрации, когда запоминание сигнала на любое необходимое время осуществляется с помощью двоичных элементов (триггеров). Важным параметром дискретного фильтра является память Т. Этот параметр оказывает решающее влияние на основные характеристики фильтра - передаточную функцию и импульсную характеристику.