После задания функции fvtool (filter visualization tool) нужно нажать клавишу . Аргументы функции – коэффициенты числителя и знаменателя передаточной функции, расположенные в порядке возрастания отрицательных степеней z (см. выражение (3)). Для просмотра характеристик нескольких фильтров одновременно нужно указать несколько пар векторов в списке входных параметров функции fvtool:

>> fvtool (b1, a1, b2, a2, b3, a3)

- 20 -

После вызова функции fvtool откроется графическое окно. На панели инструментов в верхней части этого окна имеются значки характеристик фильтра. Подводя курсор к одному из значков и нажимая левую кнопку мыши, можно вывести на экран график желаемой характеристики. Можно просмотреть

• амплитудно-частотную характеристику (АЧХ, magnitude response),

• фазочастотную характеристику (ФЧХ, phase response),

• АЧХ и ФЧХ совместно,

• xарактеристику группового времени запаздывания (group delay),

• импульсную характеристику (impulse response),

• переходную характеристику (step response),

• диаграмму полюсов и нулей (pole/zero plot),

• коэффициенты фильтра (filter coefficients).

Для детального просмотра участков графика используется кнопка zoom in панели инструментов. Щёлкнув по этой кнопке, нужно затем подвести курсор мыши к той точке графика, которая должна оказаться в центре увеличенного изображения и нажать левую кнопку мыши. Можно щёлкать левой кнопкой мыши многократно, что будет приводить ко всё большему увеличению. Возврат к исходному состоянию осуществляется щелчком по правой кнопке мыши, что соответствует нажатию кнопки zoom out панели инструментов.

8. Синтез цифрового фильтра с использованием программы fdatool

Расчёт коэффициентов передаточной функции цифрового фильтра при заданных требованиях к частотной характеристике может быть произведён с использованием тех же самых функций, которые применяются для расчёта аналоговых фильтров (см. раздел 4). Отличие заключается в том, что строковый параметр ‘s’ вводить не нужно. Кроме того, все частоты (w0, wn, wp, ws) указываются в долях от частоты Найквиста (Fs/2), то есть лежат в интервале от 0 до 1. Есть и специальный пакет программ, где собраны многие из функций, рассмотренных нами по отдельности. Он носит название fdatool (filter design & analysis tool). Этот пакет содержит удобный пользовательский интерфейс и позволяет производить расчёт передаточной функции рекурсивных и нерекурсивных фильтров разнообразными методами синтеза, просматривать характеристики фильтра, анализировать изменение характеристик при квантовании коэффициентов фильтра, отсчётов входного сигнала и результатов промежуточных вычислений. Возможна работа с различными структурами фильтра. Вызов пакета осуществляется путём ввода его имени в командном окне MatLab:

>> fdatool

• 21 -

После ввода имени с клавиатуры следует нажать клавишу .

На экране монитора появится окно программы fdatool. Кроме основного меню и панели инструментов здесь содержится текущая информация о структуре фильтра, его порядке, устойчивости (Current Filter Information); график допусков для АЧХ (в дБ) (Filter Specifications), а также вкладка Design Filter для задания типа фильтра, его класса, метода синтеза, порядка фильтра, частоты дискретизации, граничных частот полос пропускания и задерживания и допустимых затуханий в этих полосах. Кроме того, имеется вкладка для исследования эффектов квантования (Set Quantization Parameters).

8.1. Задание требований к АЧХ и расчёт фильтра

На вкладке Design Filter в нижней части окна установите переключатель Filter Type в одно из следующих положений: Lowpass (ФНЧ), Highpass (ФВЧ), Bandpass (ППФ) или Bandstop (ПЗФ). Затем используйте переключатель Design Method. Если выбрать рекурсивный фильтр, иначе БИХ-фильтр (IIR – Infinite Impulse Response), то далее в раскрывающемся списке нужно указать класс фильтра (Batterworth (Баттерворта), Chebyshev Type I (Чебышёва), Chebyshev Type II (инверсный Чебышёва), Elliptic (эллиптический)). При синтезе этих фильтров используется метод билинейного z-преобразования. В случае синтеза нерекурсивного фильтра (КИХ-фильтра, FIR – Finite Impulse Response) возможны методы: Equiripple (метод Ремеза, обеспечивающий равномерные пульсации АЧХ), Least-Squares (обеспечение минимума среднего квадратического отклонения АЧХ от заданной), Window (использование окон в качестве весовых функций при синтезе фильтра) и др.

В разделе Filter Order укажите требуемый порядок фильтра или установите переключатель в положение Minimum order (наименьший возможный порядок).

Далее перейдите к разделам Filter Specifications и Magnitude Specifications. Последовательно подводите курсор мыши к полям ввода параметров и вводите желаемые значения с клавиатуры. Смысл параметров можно понять из расположенного в верхней части окна графика допусков (Filter Specifications). Следует ввести частоту дискретизации Fs, граничные частоты полосы пропускания и полосы задерживания (Fpass и Fstop), допустимые затухания в полосе пропускания и в полосе задерживания (Apass и Astop).

После задания всех параметров щёлкните по кнопке Design Filter, расположенной в самом низу. Будет произведён расчёт, после чего можно просмотреть характеристики синтезированного фильтра.

- 22 -

8.2 Просмотр характеристик фильтра

Вывести на экран частотные и временные характеристики фильтра, диаграмму полюсов и нулей, коэффициенты фильтра можно точно так же, как это делается в программе fvtool (см. раздел 7). Для перехода к просмотру характеристик удобно вывести на экран специальное окно просмотра путём нажатия кнопки Full View Analysis панели инструментов.

8.3 Экспорт и импорт описания фильтра

Выберите команду Export меню File (+E). В появившемся окне экспорта укажите область, куда будут переданы коэффициенты фильтра (Workspace – рабочая область MatLab, Text-file – текстовый файл, MAT-file – МАТ-файл, который затем можно загрузить в MatLab командой load). Задайте имена переменных для записи векторов коэффициентов числителя (Numerator) и знаменателя (Denominator) передаточной функции (по умолчанию Num и Den). Если фильтр имеет каскадную структуру (second-order-sections), то экспортируется матрица SOS и коэффициент усиления G. Изменить структуру фильтра можно щёлкнув по кнопке Convert Structure и далее произведя выбор из списка (в версии MatLab 6.1). В версии MatLab 6.5 изменение структуры осуществляется путём выбора соответствующей команды в меню Edit. Если в рабочей области MatLab уже есть переменные с указанными именами, то установите флажок Overwrite existing variables, чтобы их значения были заменены новыми. Щёлкните по кнопке Apply. При экспорте в файл будет запрошено имя файла.

Если программу fdatool нужно использовать не для расчёта, а только для анализа характеристик фильтра, в том числе для анализа их изменения при квантовании коэффициентов и переменных фильтра, то применяется процедура импорта описания фильтра. Выберите команду Import Filter (+ I) меню Filter (в версии MatLab 6.1) или меню File (в версии MatLab 6.5). Вкладка Design Filter окна программы fdatool будет заменена на вкладку Import Filter. Укажите тип структуры (Filter Structure), выбрав нужную строку в раскрывающемся списке. В полях ввода укажите коэффициенты числителя и знаменателя передаточной функции фильтра. Список заключается в квадратные скобки, значения коэффициентов разделяются пробелами. Например, если передаточная функция имеет вид

0.2 – 0.42z ^–1 +0.05z ^–2

K (z) =

1+ 0.18z ^–1 – 0.24z ^–2 + 0.081z ^–3 ,

то в поле Numerator следует ввести [0.2 –0.42 0.05], а в поле Denominator

[1 0.18 –0.24 0.081].

- 23 -

В поле Sampling Frequency введите частоту дискретизации Fs, указав единицу измерения в поле Units.

Если векторы коэффициентов фильтра b и a, а также частота дискретизации Fs уже существуют в рабочей области MatLab, то вместо ввода численных значений в соответствующих полях ввода нужно просто указать имена переменных.

Щёлкните по кнопке Import Filter. Далее можно просматривать характеристики фильтра.

Для возврата в режим расчёта фильтров используется команда Design Filter меню Filter (в версии MatLab 6.1) или кнопка с тем же названием, расположенная в вертикальном ряду кнопок у левой границы окна (в версии

MatLab 6.5).

9. Исследование влияния квантования сигналов и коэффициентов фильтра на характеристики фильтра

9.1. Влияние квантования коэффициентов на АЧХ и другие характеристики, а также на устойчивость фильтра

Округление коэффициентов b и a передаточной функции цифрового фильтра, т.е. представление их числами с ограниченным количеством двоичных разрядов, приводит к искажению частотных и временных характеристик фильтров. Искажения могут быть как незначительными, так и более существенными. При проектировании цифровых фильтров желательно установить, какова минимально допустимая разрядность представления коэффициентов фильтра, при которой характеристики не выходят за пределы заданных полей допусков, а дисперсия шума квантования при обработке какого-либо цифрового сигнала данным фильтром не превышает допустимой величины. Если коэффициенты и переменные фильтра представлены числами с меньшим количеством двоичных разрядов, то достигается экономия памяти и повышается быстродействие фильтра, поэтому нахождение наименьшей разрядности – весьма важная задача.

Для исследования эффектов квантования можно использовать программу fdatool, речь о которой уже шла в разделе 8.

Вызовите программу fdatool:

>>fdatool

Произведите расчёт фильтра или импортируйте его описание (см. разд. 8). Установите флажок Turn quantization on в поле Quantization, расположенном в левой верхней части окна (в версии MatLab 6.1) или нажмите соответствующую кнопку на панели инструментов (в версии MatLab 6.5). Перейдите на вкладку Set quantization parameters (в нижней части окна). Задайте нужную структуру, щёлкнув по кнопке Convert structure,

- 24 -

расположенной в поле Current Filter Information, и далее произведя выбор из раскрывающегося списка (в версии Matlab 6.1), или задайте команды Convert Structure, Convert to Second-Order Sections или Convert to Single Section меню Edit (в версии MatLab 6.5). Возможно исследование структур: Direct form I (прямая), Direct form II (каноническая), Direct form I transposed (транспонированная каноническая), Direct form II transposed (транспонированная прямая), Second-order sections (каскадная). В правой части раздела Set Quantization Parameters приведён набор параметров квантователей, которые используются программой. Изменять нужно только значение в поле Format строки Coefficient (Convert coefficient to) (первая строка, последний столбец). Здесь задаётся формат двоичного представления коэффициентов фильтра. Например, формат [16 15] означает, что всего разрядов 16, а 15 из них отводится под представление дробной части (нужно учесть, что старший разряд – знаковый). Такой формат применим, если коэффициенты представлены в нормализованной форме (их модули не превышают единицы). Просмотрите список коэффициентов, нажав кнопку

[b a] панели инструментов, расположенной в верхней части окна, и оцените, сколько двоичных разрядов нужно отводить на представление целой части (с учётом знакового), чтобы правильно представить целую часть коэффициента с наибольшим модулем. Последовательно изменяйте формат, уменьшая количество разрядов (разность между первым и вторым числами формата должна при этом сохраняться неизменной). Для изменения формата подведите курсор мыши к полю ввода (курсор преобразуется в вертикальную черту) и нажмите левую кнопку мыши. Затем, используя клавишу  и цифровые клавиши, введите новые значения чисел формата. После ввода щёлкните левой кнопкой мыши вне таблицы параметров квантователей, чтобы активизировать кнопку Apply. Щёлкните по кнопке Apply. Просмотрите характеристики фильтра, нажимая соответствующие кнопки панели инструментов. Выводятся сразу две характеристики: исходная (Reference) и полученная в результате квантования (Quantized). Просмотрев критические участки АЧХ в укрупнённом масштабе (используйте кнопки изменения масштаба Zoom In и Zoom Out), определите, не выходит ли АЧХ за пределы заданного поля допусков. Если АЧХ ещё удовлетворяет требованиям, продолжите изменение формата представления коэффициентов. В результате анализа найдите, при каком наименьшем количестве двоичных разрядов АЧХ ещё удовлетворяет предъявляемым к ней требованиям. Посмотрите, как изменяется диаграмма полюсов и нулей, какими стали коэффициенты фильтра. Интересно пронаблюдать и изменения других характеристик (ФЧХ, импульсной, переходной, группового времени запаздывания).