3. С использованием функции fvtool построить АЧХ, ФЧХ, импульсные и переходные характеристики цифровых фильтров, приведенных в табл. 2.1. Значения коэффициентов a и b выбирать по табл. 2.2 в соответствии с вариантом задания.
Аргументами функции fvtool являются векторы, составленные из коэффициентов передаточной функции, начиная с коэффициента при старшей степени z. Следует обратить внимание, что порядок фильтра (старшая степень z) определяется как максимальное значение порядков числителя и знаменателя передаточной функции, при этом недостающие слагаемые со старшими степенями z необходимо доопределить, считая соответствующие коэффициенты нулевыми.
Табл. 2.1. Цифровые фильтры
№
Тип фильтра
Передаточная функция
Интегрирующий
Изодромный
Аналог апериодического
звена 1 порядка
Дифференцирующий
Дифференцирующий с пропорциональной составляющей
Интегродифференцирующий
Интегрирующий по формуле трапеций
Интегрирующий по формуле Симпсона "1/3"
Полосовой
Табл. 2.2. Варианты заданий
Вариант
a
0,5
0,6
0,7
0,8
0,9
b
0,3
0,4
0,5
0,6
0,7
Например, фильтр с передаточной функцией имеет порядок 1 (в знаменателе присутствует z в первой степени), так что вектор коэффициентов числителя будем формировать как [0 1] (в числителе коэффициент при z в первой степени равен нулю, при z в нулевой степени – единице); соответственно вектор коэффициентов знаменателя будет иметь вид [1 –1]. Таким образом, в командной строке Matlab следует ввести
>> fvtool([0 1], [1 -1])
В то же время нулевые коэффициенты при младших степенях z вводить не обязательно, например, применительно к фильтру с передаточной функцией записи fvtool([1 0 1 0 1], [1 0 0 0 0]), fvtool([1 0 1 0 1], [1]) и fvtool([1 0 1 0 1]) будут эквивалентны.
Функция fvtool вызывает специальное окно, в котором предусмотрена возможность выбора интересующей характеристики фильтра (с помощью меню Analysis либо панели инструментов), формы ее представления, изменения масштаба графика и т.п. Для изменения вида отображаемого графика следует воспользоваться меню Analysis, выбрав пункт Analysis Parameters …, либо вызвать щелком правой кнопки мыши одно из контекстных меню. В частности:
- графики АЧХ (Magnitude Response) для всех фильтров, кроме интегрирующих, следует строить в линейном масштабе, для этого щелчком правой кнопки мыши по надписи Magnitude (dB), располагающейся вдоль вертикальной оси графика, вызвать контекстное меню и выбрать Magnitude; АЧХ интегрирующих фильтров рассчитывать в децибелах (Magnitude Response, dB);
- при построении импульсных и переходных характеристик желательно соединять полученные дискретные отсчеты прямыми линиями, для этого щелчком правой кнопки мыши по области построения соответствующего графика вызвать контекстное меню и выбрать Line;
- если число отображаемых отсчетов импульсной либо переходной характеристики недостаточно для выявления всех ее особенностей, щелчком правой кнопки мыши по области построения соответствующего графика вызвать контекстное меню, выбрать пункт Analysis Parameters, Specify length: Specified и задать количество отсчетов в строке Length.