Преобразуйте lowpass аналоговые фильтры к bandstop
[bt,at] = lp2bs(b,a,Wo,Bw)
[At,Bt,Ct,Dt] = lp2bs(A,B,C,D,Wo,Bw)
lp2bs преобразовывает аналог lowpass прототипы фильтра с сокращением угловая частота 1 рад/с в заграждающие фильтры с желаемой пропускной способностью и центральная частота. Преобразование является одним шагом в процессе создания цифровых фильтров для butter, cheby1, cheby2 и функций ellip.
lp2bs может выполнить преобразование на двух различных представлениях линейной системы: форма передаточной функции и форма пространства состояний. В обоих случаях входная система должна быть аналоговым прототипом фильтра.
[bt,at] = lp2bs(b,a,Wo,Bw) преобразовывает аналог lowpass прототип фильтра, данный полиномиальными коэффициентами в заграждающий фильтр с центральной частотой Wo и пропускная способность Bw. Векторы - строки b и a задают коэффициенты числителя и знаменатель прототипа в убывающих степенях s.
Скаляры Wo и Bw задают центральную частоту и пропускную способность в модулях рад/с. Для фильтра с ребром нижней полосы w1 и ребро верхней полосы w2, используйте Wo = sqrt(w1*w2) и Bw = w2-w1.
lp2bs возвращается, преобразованная частота просачиваются векторы - строки bt и at.
[At,Bt,Ct,Dt] = lp2bs(A,B,C,D,Wo,Bw) преобразовывает непрерывно-разовое пространство состояний lowpass прототип фильтра в матрицах A, B, C, D, показанный ниже
в заграждающий фильтр с центральной частотой Wo и пропускная способность Bw. Для фильтра с ребром нижней полосы w1 и ребро верхней полосы w2, используйте Wo = sqrt(w1*w2) и Bw = w2-w1.
Заграждающий фильтр возвращен в матрицах At, Bt, Ct, Dt.
lp2bs является очень точной формулировкой пространства состояний классического аналогового преобразования частоты фильтра. Если заграждающий фильтр должен иметь центральную частоту ω0 и пропускная способность B w, стандартный s - доменное преобразование
где Q = ω0/Bw и p = s/ω0. Версия пространства состояний этого преобразования
Q = Wo/Bw; At = [Wo/Q*inv(A) Wo*eye(ma);-Wo*eye(ma) zeros(ma)]; Bt = -[Wo/Q*(A\B); zeros(ma,n)]; Ct = [C/A zeros(mc,ma)]; Dt = D - C/A*B;
Смотрите lp2bp для деривации полосовой версии этого преобразования.
