- •Оглавление
- •1. Общие вопросы цифровой обработки сигналгов
- •1.1. Основные расчетные алгоритмы для цифровых фильтров
- •2. Дискретные линейные системы
- •2.1. Модель дискретной линейной системы
- •2.2. Линейное разностное уравнение первого порядка
- •2.3. Частотная характеристика цепи первого порядка
- •2.4. Геометрическая интерпретация частотной характеристики
- •2.6. Обратное z-преобразование
- •2.7. Теорема о свертке
- •2.8. Теорема о комплексной свертке
- •2.9. Решение разностных уравнений первого порядка с помощью z-преобразования
- •2.10. Решение разностных уравнений второго порядка с помощью z-преобразования
- •2.11. Двустороннее z-преобразование
- •2.12. Цепи для разностного уравнения второго порядка
- •3. Расчет цифровых фильтров в частотной области
- •3.1. Синтез цифровых фильтров
- •3.2. Различные методы расчета цифровых фильтров
- •3.3. Применение принципа инвариантности импульсной характеристики
- •3.4. Коэффициент передачи цифровых резонаторов
- •3.5. Расчет цифровых фильтров на основе непрерывных фильтров с нулями на бесконечности
- •3.6. Определение цифрового фильтра с помощью квадрата модуля передаточной функции
- •3.7. Расчет цифровых фильтров путем билинейного преобразования функции непрерывного фильтра
- •3.8. Фильтры на основе частотной выборки
- •3.9 Метод частотной выборки
- •4. Эффекты квантования в цифровых фильтрах
- •4.1. Постановка задачи
- •4.2. Ошибки, вызываемые неточными значениями постоянных параметров
- •4.3. Ошибки, вызываемые аналого-цифровым преобразованием
- •4.4. Ошибки, вызываемые квантованием произведений
- •4.5. Эффект мертвой зоны
- •4.6. Формулы для шума округления при различных реализациях цифровых цепей
- •4.7. Пример. Различные структуры цепи с двумя полюсами и одним нулем
- •5. Дискретные преобразования фурье
- •5.1. Дискретное преобразование Фурье
- •5.2. Алгоритм Герцеля
- •5.3. Быстрое преобразование Фурье
- •Прореживание по времени
- •Прореживание по частоте
- •5.4. Соотношение между прореживанием по времени и прореживанием по частоте
- •Заключение
- •394026 Воронеж, Московский просп., 14
2.2. Линейное разностное уравнение первого порядка
В качестве примера простого линейного разностного уравнения первого порядка рассмотрим
(2.1)
На рис 2.1 представляет собой дискретный входной сигнал. Исключив квантизатор, мы сделаем равным , так что может не рассматриваться. Цифровая машина на рис. 2.1 принимает последовательность ; каждый раз поступает новое входное значение, производится вычисление и образуется выходной сигнал , который подается на декодер или запоминается для дальнейшей обработки в вычислительной машине. Можно легко представить себе структуру машины, которая выполняет операции, требуемые в (2.1). Необходимо иметь три регистра памяти: один для , другой для и третий для К. В машине образуется произведение , содержимое регистра для суммируется с этим произведением, а результат запоминается в регистре для . Теперь цифровая машина снова готова к приему следующего входного отсчета.
Сразу после того, как поступит новый входной отсчет {и начинается реализация алгоритма), состояние системы полностью определится числом в регистре. Без потери общности можно назвать этим числом в предположении, что на цифровую машину должно поступить . Из (2.1) видно, что , и по индукции можно определить общий результат:
(2.2)
Таким образом, решение уравнения (2.1) полностью определяется выражением (2.2).
Если входным сигналом на рис. 2.1 является единичная ступенчатая функция, то последовательность представляет собой единицу для .В этом случае (2.2) принимает вид
(2.3)
Для больших значений n решение стремится к установившемуся значению . Однако это решение справедливо только для ; в противоположном случае, как это можно видеть из (2.2), безгранично возрастает.
Рис. 2,2 показывает в зависимости от n для различных значений К при . Мы видим сходство между этими кривыми и откликом RC-цепочки на ступенчатую функцию, показанным на рис. 2.3. Сделать равносильно введению отрицательного сопротивления в аналоговую RC-цепочку. В дальнейшем мы убедимся в том, имеется большое сходство между аналоговыми и цифровыми цепями. Однако оно означает не эквивалентность, а только подобие.
Рис. 2.2 Отклик фильтра первого порядка на ступенчатый входной сигнал
Удобная иллюстрация вычислений но формуле (2.1) показана на рис. 2.4. Алгоритм вычислений можно объяснить следующий образом: когда поступает новый отсчет данных , предыдущее вычисленное значение выходного сигнала для новой итерации становится . Выполняется умножение и сложение, и результат, вводимый в выходной регистр, становится равным .
Назовем цепь на рис. 2.4 цифровой, или дискретной цепью, или, если не будет недоразумений, просто цепью. В этой книге будет обсуждаться и анализироваться много таких цепей, поэтому полезно сделать несколько замечаний об обозначениях.
Рис. 2.3 Отклик RC-цепи на ступенчатый входной сигнал
Рис. 2.4 Цифровая цепь первого порядка
Прямоугольник с внутри будет обозначать блок задержки на время, равное интервалу дискретизации Т. Круг со знаком плюс означает сложение всех входных сигналов (стрелки, направленные в круг) с одним выходным сигналом. Умножение на постоянный коэффициент представляется коэффициентом, записанным около линии, изображающей путь сигнала.
Интересно выяснить, какова же будет действительная схемная реализация (2.1). Возьмем несколько регистров: K — для хранения коэффициента K, Y—для хранения и X—длz хранения . Затем, при поступлении нового отсчета в X, осуществляется операция , и на этом завершается одна итерация.
Таким образом, мы видим, что необходимо иметь три регистра, а также осуществлять операции сложения и умножения. Ясно, что реализация (2.1) является более дорогой по сравнению с пассивной аналоговой цепью, за исключением, может быть, случая работы на очень низких частотах, когда величины емкости становятся большими. Однако если необходимо реализовать несколько, а не одну-единственную цепь, то в каждую из них необходимо ставить два регистра, эквивалентные K и Y, а остальные элементы могут быть общими. Так, если должно быть реализовано 500 цепей и если быстродействие арифметического блока достаточно высокое, то все цепи могут использовать один и тот же сумматор и умножитель. Например, пусть частота дискретизации 100 Гц и 500 разностных уравнений должны бить обработаны за 10 мсек. Тогда упомянутый выше ряд операций должен быть произведен менее чем за 20 мксек каждый для того, чтобы получить возможность работать в реальном масштабе времени.