Хм,это работает только на основной частоте. По сути это полный аналог исходного варианта,только с двумя параллельными ветвями. Но вопрос то не в этом,нужно снизить в два раза тактовую частоту и за один такт обрабатывать сразу два входных отсчёта. С КИХ всё понятно,разделили на чётный и нечётный и вперёд. А здесь пока не выходит.

Четные отсчеты и нечетные фильтруются параллельно фильтрами с H(z) = 1 / (1 - a^2 * z^-1); Получаем y1(n), y2(n)

Четный выход формируется как y1(n)+a*y2(n-1)
Нечетные - y2(n)+a*y1(n)

Сообщение отредактировал thermit - Feb 17 2016, 14:25

Спасибо за ответ! Пока взял а = 1. Да,для z^-1 по вашим формулам всё работает замечательно,но для z^-2 всё валится.
На картинке сверху реализация для z^-1
На картинке снизу реализация для z^-2
Может я чего не так сделал? Кстати,для нечётного канала никакая задержка на такт не нужна

Как-то так должно быть.



При чем тут z^-2? Вы сами понижаете частоту дискретизации в 2 раза и z^-2 становится z^-1.

Спасибо! Беру тайм-аут для переваривания.



Очень боюсь что мы тут друг друга недопонимаем.А чем тогда становится z^-1? Ведь тот узел что я сделал прекрасно выполняет свою функцию.
Вот этот узел у Вас является эталонным для сравнения? Там стоит z^-1,а у меня в эталоне стоит z^-2.

Мне нужна полифазная реализация БИХ фильтра H(z) = 1 / (1 - a*z^-2) для двух фаз. Начать решил с z^-1.
И я уже сказал что для z^-1 по Вашим формулам всё прекрасно работает. А вот перейти к z^-2 не удаётся.

В ваших формулах показано, как организовать вычисления для рекурсивного фильтра 1-го порядка на половинной частоте дискретизации при помощи параллельной обработки рекурсивными фильтрами 1-го порядка пары отсчетов. Зачем в эталоне лепить z^-2 совершенно не понятно. Эталон ведь первого порядка.

Ещё раз по полочкам.
Конечным результатом мне нужна полифазная реализация БИХ фильтра H(z) = 1 / (1 - a*z^-2) для двух фаз.
На первом этапе для начала я решил упростить задачу и реализовать БИХ фильтр H(z) = 1 / (1 - a*z^-1) и с помощью ваших рекомендаций мне это удалось.Для этой реализации в качестве эталона я брал БИХ фильтр H(z) = 1 / (1 - a*z^-1). Этот этап был промежуточным.
На втором этапе нужно сделать то же, но уже для фильтра H(z) = 1 / (1 - a*z^-2) и для этой реализации конечно нужен эталон с H(z) = 1 / (1 - a*z^-2).

Понятно. Проще некуда.
Четные и нечетные отсчеты фильтруются параллельно фильтрами H(z) = 1 / (1 - a*z^-1). Все.

Или распараллелить на 4 фазы

Странно, но я же уже делал параллельно фильтры H(z) = 1 / (1 - a*z^-1), о чём сообщал выше и выход с этих фильтров сравнивался с эталонным фильтром работающем на основной частоте. Эталонный фильтр имел H(z) = 1 / (1 - a*z^-1) и все данные совпадали. Не может же эта самая схема выдавать отклик совпадающий и с эталонным фильтром H(z) = 1 / (1 - a*z^-2)???

Значит не совсем совпадали.

o'k! Завтра вышлю модельку.

Вот моделька двухфазной реализации фильтра H(z) = 1 / (1 - a*z^-1) при а = 1. Похоже что совпадает со всем.
Огромное Вам спасибо за модель parallel_vs_serial_4ph! Буду анализировать. Не сможете ли поделиться своими выкладками(формулами) по синтезу этой схемы?
А для двух фаз вариант с z^-2 не реализуем?

Всё у меня получилось! Огромное Вам спасибо за потраченное время!
Вот вариант двухфазного фильтра для z^-2.

Тогда, резюме:

1.
Параллельное вычисление 2-х фаз





2.
Параллельное вычисление 2-х фаз


Параллельное вычисление 4-х фаз

Спасибо! Добавим в копилку.