MATLAB.Exponenta
MATLAB и Simulink на русском
Технологии разработки и отладки
		сложных технических систем

Обработка сигналов и изображений\Signal Processing Toolbox

Список функций Signal Processing Toolbox: Функции синтеза фильтров

В оглавление \ К следующему разделу \ К предыдущему разделу

CHEBY2
Синтез фильтра Чебышева второго рода (пульсации в полосе задерживания)

Синтаксис

[b,a] = cheby2(n,Rs,Wn)

[b,a] = cheby2(n,Rs,Wn,'ftype')

[b,a] = cheby2(n,Rs,Wn,'s')

[b,a] = cheby2(n,Rs,Wn,'ftype','s')

[z,p,k] = cheby2(...)

[A,B,C,D] = cheby2(...)

Описание

Функция cheby2 позволяет синтезировать аналоговые и дискретные фильтры Чебышева второго рода нижних частот, верхних частот, полосовые и режекторные. Фильтры Чебышева второго рода имеют АЧХ, монотонно спадающую в полосе пропускания и равномерно пульсирующую в полосе задерживания. По сравнению с фильтрами Чебышева первого рода фильтры второго рода обеспечивают менее крутой спад АЧХ в переходной зоне между полосами пропускания и задерживания, зато их АЧХ не имеет пульсаций в полосе пропускания.

Синтез дискретных фильтров

[b,a] = cheby2(n,Rs,Wn)

Производит синтез дискретного фильтра Чебышева второго рода n-го порядка, имеющего АЧХ фильтра нижних частот, нормированную частоту среза Wn, и обеспечивающего подавление сигнала в полосе задерживания, равное Rs дБ. Функция возвращает описание фильтра в виде векторов-строк b и a, имеющих длину n+1 и содержащих коэффициенты полиномов числителя и знаменателя функции передачи в порядке убывания степеней переменной z:

Частотой среза фильтра Чебышева второго рода называется начало полосы задерживания, то есть частота, на которой модуль коэффициента передачи равен -Rs дБ. Функция cheby2 использует значение частоты среза Wn, нормированное к частоте Найквиста (p радиан на отсчет). Таким образом, нормированное значение частоты среза Wn должно лежать в пределах от 0 до 1, при этом 1 соответствует частоте Найквиста (то есть половине частоты дискретизации). При увеличении подавления сигнала в полосе задерживания Rs происходит расширение переходных зон между полосами пропускания и задерживания (то есть уменьшается крутизна спада АЧХ).

Если параметр Wn является двухэлементным вектором, Wn = [w1 w2], функция cheby2 производит синтез дискретного полосового фильтра порядка 2*n, полоса пропускания которого лежит в пределах от w1 до w2.

[b,a] = cheby2(n,Rs,Wn,'ftype')

Производит синтез фильтра верхних частот или режекторного фильтра, в зависимости от значения строкового параметра 'ftype'. Возможны следующие варианты:

  • 'high' - для синтеза дискретного фильтра верхних частот с нормированной частотой среза Wn;
  • 'stop' - для синтеза дискретного режекторного фильтра порядка 2*n. В этом случае входной параметр Wn должен быть двухэлементным вектором, Wn = [w>1 w2]. Полоса задерживания синтезированного фильтра лежит в пределах от w1 до w2.

Функция cheby2 может представлять синтезированный дискретный фильтр различными способами, в зависимости от числа указанных при вызове функции выходных параметров. Чтобы получить описание синтезированного фильтра в виде коэффициента усиления и наборов нулей и полюсов (zero-pole-gain form), используйте три выходных параметра, как показано ниже:

[z,p,k] = cheby2(n,Rs,Wn)

или

[z,p,k] = cheby2(n,Rs,Wn,'ftype')

Возвращает нули и полюсы функции передачи фильтра в виде векторов-столбцов z и p длины n (или 2*n, в зависимости от типа фильтра), а коэффициент усиления - в виде скаляра k.

Чтобы получить описание синтезированного фильтра в виде параметров пространства состояния (state-space form), используйте четыре выходных параметра, как показано ниже:

[A,B,C,D] = cheby2(n,Rs,Wn)

или

[A,B,C,D] = cheby2(n,Rs,Wn,'ftype')

Возвращает параметры пространства состояния A, B, C и D, описывающие фильтр следующим образом:

Здесь u - входной сигнал, x - вектор состояния, y - выходной сигнал.

Синтез аналоговых фильтров

[b,a] = cheby2(n,Rs,Wn,'s')

Производит синтез аналогового фильтра Чебышева второго рода n-го порядка, имеющего АЧХ фильтра нижних частот, круговую частоту среза, равную Wn рад/с, и обеспечивающего подавление сигнала в полосе задерживания, равное Rs дБ. Функция возвращает описание фильтра в виде векторов-строк b и a, имеющих длину n+1 и содержащих коэффициенты полиномов числителя и знаменателя функции передачи в порядке убывания степеней переменной s:

Частотой среза фильтра Чебышева второго рода называется частота, на которой модуль коэффициента передачи равен -Rs дБ. Значение круговой частоты среза Wn должно быть положительным.

Если параметр Wn является двухэлементным вектором, Wn = [w1 w2], вызов функции cheby2(n,Rs,Wn,'s') производит синтез аналогового полосового фильтра порядка 2*n, полоса пропускания которого лежит в пределах от w1 до w2.

[b,a] = cheby2(n,Rs,Wn,'ftype','s')

Производит синтез аналогового фильтра верхних частот или режекторного фильтра, в зависимости от значения строкового параметра 'ftype'. Возможные значения этого параметра такие же, как и при расчете дискретных фильтров (см. выше).

Функция cheby2 может представлять синтезированный аналоговый фильтр различными способами, в зависимости от числа указанных при вызове функции выходных параметров. Чтобы получить описание синтезированного фильтра в виде коэффициента усиления и наборов нулей и полюсов (zero-pole-gain form), используйте три выходных параметра, как показано ниже:

[z,p,k] = cheby2(n,Rs,Wn,'s')

или

[z,p,k] = cheby2(n,Rs,Wn,'ftype','s')

Возвращает нули и полюсы функции передачи фильтра в виде векторов-столбцов z и p длины n (или 2*n, в зависимости от типа фильтра), а коэффициент усиления - в виде скаляра k.

Чтобы получить описание синтезированного фильтра в виде параметров пространства состояния (state-space form), используйте четыре выходных параметра, как показано ниже:

[A,B,C,D] = cheby2(n,Rs,Wn,'s')

или

[A,B,C,D] = cheby2(n,Rs,Wn,'ftype','s')

Возвращает параметры пространства состояния A, B, C и D, описывающие фильтр следующим образом:

Здесь u - входной сигнал, x - вектор состояния, y - выходной сигнал.

Примеры Пример 1: фильтр нижних частот

Пусть частота дискретизации равна 1000 Гц. Синтезируем дискретный фильтр нижних частот Чебышева второго рода девятого порядка с частотой среза, равной 300 Гц (что соответствует нормированному значению 0,6), и подавлением сигнала в полосе задерживания, равным 20 дБ:

[b,a] = cheby2(9,20,300/500);

Построим графики АЧХ и ФЧХ синтезированного фильтра:

freqz(b,a,512,1000)

title('n=9 Chebyshev Type II Lowpass Filter')

Пример 2: полосовой фильтр

Синтезируем полосовой дискретный фильтр Чебышева второго рода 10-го порядка с полосой пропускания, простирающейся от 100 до 200 Гц, и построим график его импульсной характеристики:

n = 5;

r = 20;

Wn = [100 200]/500;

[b,a] = cheby2(n,r,Wn);

[y,t] = impz(b,a,101);

stem(t,y)

title('Impulse Response of n=5 Chebyshev Type II Filter')

Ограничения

Для фильтров высокого порядка наибольшая численная точность достигается при использовании представления в пространстве состояний. Следующей по степени точности идет форма представления в виде нулей и полюсов. Представление фильтра в виде коэффициентов полиномов функции передачи является наименее точным; проблемы могут возникать уже для фильтров 15-го порядка.

Алгоритм

Функция cheby2 использует алгоритм, состоящий из пяти шагов:

  1. Производится расчет нулей, полюсов и коэффициента усиления фильтра-прототипа с помощью функции cheb2ap.
  2. Полученное описание фильтра трансформируется в параметры пространства состояний.
  3. В пространстве состояний производится трансформация фильтра нижних частот в фильтр верхних частот, полосовой или режекторный фильтр с требуемыми частотами среза (согласно параметрам, указанным при вызове функции).
  4. При синтезе дискретных фильтров функция cheby2 использует билинейное преобразование (функция bilinear) для дискретизации аналогового фильтра. При этом используется предыскажение частот среза, гарантирующее совпадение коэффициентов передачи аналогового и дискретного фильтров на этих частотах.
  5. Если при вызове функции указано два или три выходных параметра, производится преобразование описания фильтра из параметров пространства состояний в коэффициенты полиномов функции передачи или в наборы нулей и полюсов.

Сопутствующие функции: besself, butter, cheb2ap, cheb1ord, cheby1, ellip

В оглавление \ К следующему разделу \ К предыдущему разделу


Поиск по сайту:

Система Orphus

Яндекс.Метрика