tf2cl

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

Синтаксис

[k1,k2] = tf2cl(b,a)
[k1,k2] = tf2cl(b,a)
[k1,k2,beta] = tf2cl(b,a)

Описание

[k1,k2] = tf2cl(b,a) где b - действительный, симметричный вектор коэффициентов числителя и a - вектор действительных чисел коэффициентов знаменателя, соответствующий стабильному цифровому фильтру, будет выполнять связанное allpass разложение

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

стабильного БИХ H(z) и преобразует передаточные функции allpass H1(z) и H2(z) в связанную структуру allpass решетки с коэффициентами, заданными в векторах k1 и k2.

[k1,k2] = tf2cl(b,a) где b является действительным, антисимметричным вектором коэффициентов и a числителя - вектор действительных чисел коэффициентов знаменателя, соответствующий стабильному цифровому фильтру, выполняет связанное allpass разложение

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

стабильного БИХ H(z) и преобразует передаточные функции allpass H1(z) и H2(z) в связанную структуру allpass решетки с коэффициентами, заданными в векторах k1 и k2.

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

[k1,k2,beta] = tf2cl(b,a) выполняет обобщенное allpass разложение стабильного БИХ фильтра H(z) и преобразует комплексные allpass передаточных функций H1(z) и H2(z) в соответствующие решетчатые allpass фильтры

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

где beta является комплексным скаляром величины, равной 1.

Примечание

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

Примеры

[b,a]=cheby1(9,.5,.4);
[k1,k2]=tf2cl(b,a); % Get the reflection coeffs. for the lattices.
[num1,den1]=latc2tf(k1,'allpass'); % Convert each allpass lattice
[num2,den2]=latc2tf(k2,'allpass'); % back to transfer function.
num = 0.5*conv(num1,den2)+0.5*conv(num2,den1);
den = conv(den1,den2); % Reconstruct numerator and denonimator.
MaxDiff=max([max(b-num),max(a-den)]); % Compare original and reconstructed
							% numerator and denominators.

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

.

См. также

| | | | |

Введенный в R2011a