tf2ca

Передаточная функция к связанному allpass

Синтаксис

[d1,d2] = tf2ca(b,a)
[d1,d2] = tf2ca(b,a)
[d1,d2,beta] = tf2ca(b,a)

Описание

[d1,d2] = tf2ca(b,a) где b - действительный, симметричный вектор коэффициентов числителя и a - вектор действительных чисел коэффициентов знаменателя, соответствующий стабильному цифровому фильтру, возвращает векторы действительных чисел d1 и d2 содержащие коэффициенты знаменателя фильтров allpass H1(z) и H2(z) такие, что

H(z)=B(z)A(z)=(12)[H1(z)+H2(z)]

представление связанного разложения альпаса.

[d1,d2] = tf2ca(b,a) где b является действительным, антисимметричным вектором коэффициентов и a числителя - вектор действительных чисел коэффициентов знаменателя, соответствующий стабильному цифровому фильтру, возвращает векторы действительных чисел d1 и d2 содержащие коэффициенты знаменателя фильтров allpass H1(z) и H2(z) такие, что

H(z)=B(z)A(z)=(12)[H1(z)H2(z)]

В некоторых случаях разложение невозможно при действительных H1(z) и H2(z). В этих случаях может быть возможно обобщенное связанное разложение альпаса, как описано в следующем синтаксисе.

[d1,d2,beta] = tf2ca(b,a) возвращает комплексные векторы d1 и d2 содержащий коэффициенты знаменателя фильтров H1(z) и H2(z) альпаса и комплексный скаляр beta, удовлетворяющие |<reservedrangesplaceholder0>| = 1, таким что

H(z)=B(z)A(z)=(12)[β¯H1(z)+βH2(z)]

представление обобщенного разложения альпаса.

В вышеприведенных уравнениях H1(z) и H2(z) являются действительными или комплексными фильтрами allpass БИХ, заданными как

H1(z)=fliplr((D1(z))¯)D1(z),H2(1)(z)=fliplr((D2(1)(z)¯))D2(1)(z)

где D1(z) и D2(z) являются полиномами, коэффициенты которых заданы d1 и d2.

Примечание

Связанное разложение альпаса не всегда возможно. Тем не менее, фильтры Butterworth, Chebyshev, и Elliptic БИХ, среди прочих, могут быть факторизованы таким образом. Для получения дополнительной информации см. Signal Processing Toolbox™ руководство пользователя.

Примеры

[b,a]=cheby1(9,.5,.4);
[d1,d2]=tf2ca(b,a); 	% TF2CA returns denominators of the allpass.
num = 0.5*conv(fliplr(d1),d2)+0.5*conv(fliplr(d2),d1);
den = conv(d1,d2); 	% Reconstruct numerator and denonimator.
MaxDiff=max([max(b-num),max(a-den)]); % Compare original and reconstructed
							% numerator and denominators.

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

.

См. также

| | | |

Введенный в R2011a