2. Запустить Matlab, создать m-файл в соответствии с приложением 2. Ввести значения коэффициентов числителя и знаменателя передаточной функции D(z) в соответствии с вариантом задания (табл. 4.1), начиная с коэффициентов при нулевой степени z. Таким образом, первый элемент вектора a должен быть равен 1.
3. Ввести значение числа двоичных разрядов N, использующихся для хранения дробной части коэффициентов разностных уравнений (при анализе схемных реализаций цифровых фильтров) либо промежуточных результатов вычислений (при анализе предельных циклов) в соответствии с табл. 4.1. При этом цена младшего разряда будет равна 2–N.
Табл. 4.1. Варианты заданий
Вариант
D(z)
N
4. Нажатием клавиши [F5] запустить программу на выполнение. В первом графическом окне пронаблюдать графики:
- АЧХ фильтра с заданными коэффициентами разностных уравнений (красная кривая);
- АЧХ фильтра с округленными коэффициентами, реализованного в прямой форме (синяя кривая);
- АЧХ фильтра с округленными коэффициентами, реализованного в последовательной форме (зеленая кривая).
Убедиться, что АЧХ квантованного фильтра, представленного в последовательной форме, практически полностью совпадает с исходной АЧХ.
5. Во втором графическом окне пронаблюдать графики импульсных характеристик:
- исходного фильтра (красная кривая);
- фильтра с округлением результатов промежуточных вычислений (синяя кривая).
Определить амплитуду предельного цикла. Сопоставить ее значение с ценой единицы младшего разряда вычислителя.
Определить период предельного цикла. Сопоставить его значение с порядком фильтра.