Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
КСД-12 ОСИ / ЛАБ 6.doc
Скачиваний:
44
Добавлен:
03.03.2016
Размер:
69.63 Кб
Скачать

Freqz (b, a, n)

(точки n равномерно распределяются в диапазоне 0… π).

Для анализа фильтров с комплексными коэффициентами, частотные характеристики которые несимметричны, интервал частот равномерно распределяют от 0 до с помощью четвертого входного параметра – строки whole:

freqz (b, a, n, whole)

1.7. отображение нулей и полюсов фильтра.

Для отображения нулей и полюсов функции передачи фильтра на комплексной плоскости используется функция zplane:

zplane (z, p)

zplane (b, a)

т.е. входными параметрами могут быть векторы – столбцы нулей и полюсов (z, p), либо векторы – строки коэффициентов полиномов числителя и знаменателя функции передачи (b, a) фильтра.

>> [b, a] = (n, Rs, Rw)

параметры фильтра зависят от типа фильтра в части уровня пульсаций в полосе пропускания и задерживания.

Нули функций на графике отображаются кружочками, полюсы – крестиками.

1.8. Преобразование способов описания дискретных фильтров.

Преобразование описаний дискретных фильтров между тремя вариантами: (b, a); (z, p, k); (A, B, C, D) выполняется теми же шестью функциями, что и для аналоговых фильтров см. лабораторную работу №5. единственная особенность их использования - векторы b и a должны иметь одинаковую длину. Для удобства преобразования описаний дискретных фильтров в пакете есть функция eqtflength, дополняющая более короткий вектор нулями:

[b1, a1] = eqtflength (b, a)

Разложение на простые дроби в дискретной системе отличается от аналогового случая, поэтому используется специальная функция:

Reziduez

Для дискретных систем большое практическое значение имеет последовательная или каскадная реализация, когда система представляется в виде последовательного соединения секций второго порядка (second oder sections - sos).

Для преобразования описания дискретной системы в эту форму и из неё служат функции:

sos2ss, sos2tf, sos2zp, zp2sos, ss2sos, tf2sos

1.8.1. Разложение на простые дроби.

Функция residuez может выполнять преобразование в обе стороны:

[z, p, k] = residuez (b, a)

[b, a] = residuez (r, p, k)

Здесь b и a – векторы коэффициентов полиномов числителя и знаменателя функции передачи; r - вектор вычета; p – вектор полюсов; k – вектор коэффициентов целой части функции передачи.

Разложение функции передачи на простые дроби соответствует параллельной схеме реализации фильтра. Каждой простой дроби соответствует рекурсивный фильтр первого порядка. Выходные сигналы этих фильтров суммируются.

1.8.2. Функции, работающие с секциями второго порядка.

Секции второго порядка используют для преобразования двух последовательно соединенных комплексных фильтров первого порядка в вещественный фильтр второго порядка, что уменьшает количество операций при расчете выходного сигнала.

В пакете имеется шесть функций, осуществляющих преобразование описания дискретного фильтра между набором секций второго порядка и другими формами описания.

Преобразование секций второго порядка в коэффициенты полиномов числителя и знаменателя функции передачи:

[b, a] = sos2tf (sos, g)

Преобразование в параметры состояния:

[A, B, C, D] = sos2ss (sos, g)

Преобразование в нули, полюсы и коэффициенты усиления:

[z, p, k] = , sos2zp (sos, g)

И обратные преобразования:

[sos, g] = tf2sos (b, a, ‘order’, ‘scale’)

[sos, g] = zp2sos (z, p, k, ‘order’, ‘scale’)

[sos, g] = ss2sos (A, B, C, D, ‘order’, ‘scale’)

Соседние файлы в папке КСД-12 ОСИ