Преобразуйте связанную решетку allpass в форму передаточной функции
[b,a] = cl2tf(k1,k2)
[b,a] = cl2tf(k1,k2,beta)
[b,a,bp] = cl2tf(k1,k2)
[b,a,bp] = cl2tf(k1,k2,beta)
[b,a] = cl2tf(k1,k2) возвращает числитель и векторы знаменателя коэффициентов b и a соответствие передаточной функции
где H1(z) и H2(z) являются передаточными функциями фильтров allpass, определенных k1 и k2, и k1 и k2 вектора действительных чисел отражательных коэффициентов, соответствующих allpass структуры решетки.
[b,a] = cl2tf(k1,k2,beta) где k1, k2 и beta являются комплексными, возвращает числитель и векторы знаменателя коэффициентов b и соответствия передаточной функции
[b,a,bp] = cl2tf(k1,k2) где k1 и k2 действительны, возвращает векторный bp из действительных коэффициентов, соответствующих числителю степени дополнительный фильтр G(z)
[b,a,bp] = cl2tf(k1,k2,beta) где k1, k2 и beta являются комплексными, возвращает вектор коэффициентов bp из возможно объединяют коэффициенты, соответствующие числителю степени дополнительный фильтр G(z)
tf2cl возвращает отражательный coefs
[b,a]=cheby1(10,.5,.4); [k1,k2,beta]=tf2cl(b,a);
Создайте оригинал и степень дополнительный фильтр.
[num,den,numpc]=cl2tf(k1,k2,beta); [h,w]=freqz(num,den); hpc = freqz(numpc,den);
Покажите частотную характеристику
subplot(211); plot(w./pi,20*log10(abs(h)),'k'); hold on; grid on; plot(w./pi,20*log10(abs(hpc)),'b'); xlabel('Normalized Frequency (x \pi radians/sample)'); ylabel('dB'); legend('Original Filter','Power Complementary Filter',... 'Location','best'); subplot(212); plot(w./pi,unwrap(angle(h)),'k'); hold on; grid on; xlabel('Normalized Frequency (x \pi radians/sample)'); ylabel('Phase (radians)'); plot(w./pi,unwrap(angle(hpc)),'b');
