exponenta event banner

lp2bp

Преобразование низкочастотных аналоговых фильтров в полосовые

Описание

пример

[bt,at] = lp2bp(b,a,Wo,Bw) преобразует прототип аналогового фильтра нижних частот, заданный полиномиальными коэффициентами (заданными векторами строк b и a) в полосовой фильтр с центральной частотой Wo и полоса пропускания Bw. Система ввода должна быть прототипом аналогового фильтра.

[At,Bt,Ct,Dt] = lp2bp(A,B,C,D,Wo,Bw) преобразует прототип фильтра нижних частот состояния и пространства непрерывного времени (задается матрицами A, B, C, и D) к полосовому фильтру с центральной частотой Wo и полоса пропускания Bw. Система ввода должна быть прототипом аналогового фильтра.

Примеры

свернуть все

Разработка прототипа аналогового фильтра Butterworth 14-го порядка.

n = 14;
[z,p,k] = buttap(n);

Преобразование прототипа в форму функции переноса. Отображение его величины и частотных характеристик.

[b,a] = zp2tf(z,p,k);
freqs(b,a)

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

Преобразование прототипа в полосовой фильтр с полосой пропускания от 30 Гц до 100 Гц. Укажите центральную частоту и полосу пропускания в радиусе/с.

fl = 30;
fh = 100;

Wo = 2*pi*sqrt(fl*fh); % center frequency
Bw = 2*pi*(fh-fl); % bandwidth

[bt,at] = lp2bp(b,a,Wo,Bw);

Отображение амплитуды и частотных характеристик преобразованного фильтра.

freqs(bt,at)

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

Входные аргументы

свернуть все

Коэффициенты числителя и знаменателя прототипа, указанные как векторы строк. b и a задать коэффициенты числителя и знаменателя прототипа в нисходящих степенях s:

B (s) A (s) = b (1) sn+⋯+b (n) s + b (n + 1) a (1) sm+⋯+a (m) s + a (m + 1)

Типы данных: single | double

Представление состояния-пространства прототипа, указанное как матрицы. Матрицы состояния-пространства связывают вектор состояния x, вход u и выход y через

x˙=Ax+Buy=Cx+Du

Типы данных: single | double

Центральная частота, заданная как скаляр. Для фильтра с нижним краем полосы w1 и верхний край полосы w2, использовать Wo = sqrt (w1*w2). Экспресс Wo в единицах рад/с.

Типы данных: single | double

Полоса пропускания, заданная как скаляр. Для фильтра с нижним краем полосы w1 и верхний край полосы w2, использовать Bw = w2w1. Экспресс Bw в единицах рад/с.

Типы данных: single | double

Выходные аргументы

свернуть все

Преобразованные числительные и знаменательные коэффициенты, возвращаемые в виде векторов строк.

Преобразованное представление состояния-пространства, возвращаемое в виде матриц.

Алгоритмы

lp2bp преобразует аналоговые прототипы фильтров нижних частот с угловой частотой отсечки 1 рад/с в полосовые фильтры с требуемой полосой пропускания и центральной частотой. Преобразование является одним из этапов процесса проектирования цифрового фильтра для butter, cheby1, cheby2, и ellip функции.

lp2bp - высокоточная формулировка состояния-пространства классического преобразования частоты аналогового фильтра. Рассмотрим государственно-космическую систему

x˙=Ax+Buy=Cx+Du

где u - входной сигнал, x - вектор состояния, а y - выходной сигнал. Преобразование Лапласа первого уравнения (предполагающее нулевое начальное условие) равно

sX (s) = AX (s) + BU (s)

Теперь, если полосовой фильтр имеет центральную частоту λ 0 и полосу Bw, стандартное преобразование s-области

s = Q (p2 + 1 )/p

где Q = ω0/Bw и p = s/λ 0. Подставляя это для s в уравнении состояния-пространства Лапласа и рассматривая оператор p как d/dt приводит к

Qx¨+Qx=A˙x+Bu˙

или

Qx¨−A˙x−Bu˙=−Qx

Теперь определите

Qω˙=−Qx

который при замене приводит к

Qx˙=Ax+Qω+Bu

Последние два уравнения дают уравнения состояния. Запишите их в стандартном виде и умножьте дифференциальные уравнения на λ 0, чтобы восстановить временной или частотный масштаб, представленный p, и найти матрицы состояний для полосового фильтра:

Q = start0Bw

A Q глаз (ma, m); eye (ma, m) нули (ma, m)]

Bt = λ 0 [BQ; нули (ma, n)]

Ct = [C нули (mc, ma)]

Dt = D

где [m, m] = размер (А).

lp2bp может выполнять преобразование для двух различных представлений линейной системы: формы передаточной функции и формы состояния-пространства. Если вход в lp2bp находится в форме передаточной функции, функция преобразует его в форму state-space перед применением этого алгоритма.

Расширенные возможности

.

См. также

| | | |

Представлен до R2006a