met_r_02 (557391), страница 4
Текст из файла (страница 4)
Если выбрать рекурсивный фильтр, иначеБИХ-фильтр (IIR – Infinite Impulse Response), то далее в раскрывающемсясписке нужно указать класс фильтра (Batterworth (Баттерворта), ChebyshevType I (Чебышёва), Chebyshev Type II (инверсный Чебышёва), Elliptic(эллиптический)). При синтезе этих фильтров используется методбилинейного z-преобразования. В случае синтеза нерекурсивного фильтра(КИХ-фильтра, FIR – Finite Impulse Response) возможны методы: Equiripple(метод Ремеза, обеспечивающий равномерные пульсации АЧХ), LeastSquares (обеспечение минимума среднего квадратического отклонения АЧХот заданной), Window (использование окон в качестве весовых функций присинтезе фильтра) и др.В разделе Filter Order укажите требуемый порядок фильтра илиустановите переключатель в положение Minimum order (наименьшийвозможный порядок).Далее перейдите к разделам Filter Specifications и MagnitudeSpecifications.
Последовательно подводите курсор мыши к полям вводапараметров и вводите желаемые значения с клавиатуры. Смысл параметровможно понять из расположенного в верхней части окна графика допусков(Filter Specifications). Следует ввести частоту дискретизации Fs, граничныечастоты полосы пропускания и полосы задерживания (Fpass и Fstop),допустимые затухания в полосе пропускания и в полосе задерживания (Apassи Astop).После задания всех параметров щёлкните по кнопке Design Filter,расположенной в самом низу.
Будет произведён расчёт, после чего можнопросмотреть характеристики синтезированного фильтра.- 22 8.2 Просмотр характеристик фильтраВывести на экран частотные и временные характеристики фильтра,диаграмму полюсов и нулей, коэффициенты фильтра можно точно так же,как это делается в программе fvtool (см. раздел 7). Для перехода к просмотрухарактеристик удобно вывести на экран специальное окно просмотра путёмнажатия кнопки Full View Analysis панели инструментов.8.3 Экспорт и импорт описания фильтраВыберите команду Export меню File (<Ctrl>+E). В появившемся окнеэкспорта укажите область, куда будут переданы коэффициенты фильтра(Workspace – рабочая область MatLab, Text-file – текстовый файл, MAT-file –МАТ-файл, который затем можно загрузить в MatLab командой load).Задайте имена переменных для записи векторов коэффициентов числителя(Numerator) и знаменателя (Denominator) передаточной функции (поумолчанию Num и Den).
Если фильтр имеет каскадную структуру (secondorder-sections), то экспортируется матрица SOS и коэффициент усиления G.Изменить структуру фильтра можно щёлкнув по кнопке Convert Structure идалее произведя выбор из списка (в версии MatLab 6.1). В версии MatLab 6.5изменение структуры осуществляется путём выбора соответствующейкоманды в меню Edit.
Если в рабочей области MatLab уже есть переменные суказанными именами, то установите флажок Overwrite existing variables,чтобы их значения были заменены новыми. Щёлкните по кнопке Apply. Приэкспорте в файл будет запрошено имя файла.Если программу fdatool нужно использовать не для расчёта, а только дляанализа характеристик фильтра, в том числе для анализа их изменения приквантовании коэффициентов и переменных фильтра, то применяетсяпроцедура импорта описания фильтра.
Выберите команду Import Filter(<Ctrl>+ I) меню Filter (в версии MatLab 6.1) или меню File (в версии MatLab6.5). Вкладка Design Filter окна программы fdatool будет заменена на вкладкуImport Filter. Укажите тип структуры (Filter Structure), выбрав нужную строкув раскрывающемся списке. В полях ввода укажите коэффициенты числителяи знаменателя передаточной функции фильтра.
Список заключается вквадратные скобки, значения коэффициентов разделяются пробелами.Например, если передаточная функция имеет вид0.2 – 0.42z –1 +0.05z –2K(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.
Далее можно просматриватьхарактеристики фильтра.Для возврата в режим расчёта фильтров используется команда DesignFilter меню 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), или задайте команды ConvertStructure, Convert to Second-Order Sections или Convert to Single Section менюEdit (в версии MatLab 6.5).
Возможно исследование структур: Direct form I(прямая), Direct form II (каноническая), Direct form I transposed(транспонированнаяканоническая),DirectformIItransposed(транспонированная прямая), 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), определите, не выходит ли АЧХза пределы заданного поля допусков. Если АЧХ ещё удовлетворяеттребованиям,продолжитеизменениеформатапредставлениякоэффициентов.