Целочисленное проектирование фильтров

Все что касается фильтрации
Vlad27
Сообщения: 76
Зарегистрирован: 06 апр 2015, 10:09
Откуда: г, Нижний Новгород (rif27@bk.ru)

Re: Целочисленное проектирование фильтров

Сообщение Vlad27 » 18 сен 2018, 16:04

Бахурин Сергей писал(а):
18 сен 2018, 13:15
Посмотрел. Ваш фильтр выродился. Числитель его H(z) практически стал равен знаменателю и он практически никак не изменяет сигнал. Реально амплитудные колебания сотые доли дБ и фаза поворачивается на доли градуса. Импульсная характеристика очень близка к дельта импульсу.
Ну почему-же выродился? Ведь это идеальный всепропускающий (ALL-Pass) БИХ-фильтр, который синтезировался методом ЦНП как минимальнофазовый фильтр. При этом высокий порядок фильтра N=8 позволяет обеспечить постоянство модуля коэфф.передачи и линейность его фазы на всём интервале Найквиста, а выполнение условия минимальной фазы - обеспечить минимальное временем обработки сигнала, минимальную инерционность всепропускающей системы. Но коли вы уже рассчитали этот фильтр в MATLAB, то наверное стоит выложить его характеристики на форуме для осмысления и изучения. Ведь как минимальфазовые можно проектировать и все селективные фильтры типа ФНЧ, ФВЧ и ППФ, обеспечивая их высокое быстродействие при работе в реальном времени.
Вложения
Faza_ALL_4.gif
Faza_ALL_4.gif (6.56 КБ) 167 просмотров
Последний раз редактировалось Vlad27 21 сен 2018, 10:33, всего редактировалось 1 раз.

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 772
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: Целочисленное проектирование фильтров

Сообщение Бахурин Сергей » 20 сен 2018, 22:42

Почему на вашем графике минимально-фазовый allpass имеет линейную ФЧХ? Ведь минимально фазовый allpass это
H(z) = A(z) / A(z) = 1, откуда фчх должна быть равна 0 для всех частот.

Vlad27
Сообщения: 76
Зарегистрирован: 06 апр 2015, 10:09
Откуда: г, Нижний Новгород (rif27@bk.ru)

Re: Целочисленное проектирование фильтров

Сообщение Vlad27 » 21 сен 2018, 10:41

Бахурин Сергей писал(а):
20 сен 2018, 22:42
Почему на вашем графике минимально-фазовый allpass имеет линейную ФЧХ? Ведь минимально фазовый allpass это
H(z) = A(z) / A(z) = 1, откуда фчх должна быть равна 0 для всех частот.
При пректировании фильтров высокого порядка численными методами ЦНП вы можете обеспечить любые требования по ФЧХ - линейность, выпуклость или вогнутость фазы для фазовой коррекции сигнальных трактов (см. рис). Видны пределы фазовой коррекции в All-Pass при гарантированной единичной передаче Кп=1.
Кстати, адрес почты мой изменился pbox27@yandex.ru - характеристик фильтра я от вас не получил, к сожалению
Вложения
Forum_6zv_int.jpg

Vlad27
Сообщения: 76
Зарегистрирован: 06 апр 2015, 10:09
Откуда: г, Нижний Новгород (rif27@bk.ru)

Re: Целочисленное проектирование фильтров

Сообщение Vlad27 » 21 сен 2018, 10:43

Забыл вложить график ALL-Pass корректора фазы. Это простого, не минимальфазового IIR.
На минимальнофазовом ALL-Pass корректоры фазы не делают.
Вложения
Faza_ALL_ZV_2.gif
Faza_ALL_ZV_2.gif (9.81 КБ) 165 просмотров
Последний раз редактировалось Vlad27 25 сен 2018, 11:49, всего редактировалось 2 раза.

Vlad27
Сообщения: 76
Зарегистрирован: 06 апр 2015, 10:09
Откуда: г, Нижний Новгород (rif27@bk.ru)

Re: Целочисленное проектирование фильтров

Сообщение Vlad27 » 21 сен 2018, 12:36

Бахурин Сергей писал(а):
20 сен 2018, 22:42
Почему на вашем графике минимально-фазовый allpass имеет линейную ФЧХ? Ведь минимально фазовый allpass это
H(z) = A(z) / A(z) = 1, откуда фчх должна быть равна 0 для всех частот.
Да, теоретически и фаза и ГВЗ равны нулю, это модельная абстракция. Я же привожу уже измерения, работу фильтра в реальном времени. А там нулевого фазового сдвига быть, конечно, не может. Но этот сдвиг самый минимальный, как видите. Одно дело - теория, а другое - реализация на цифровой платформе. В области минимальных фазовых сдвигов железо всегда добавляет малую задержку.
Последний раз редактировалось Vlad27 06 ноя 2018, 11:31, всего редактировалось 1 раз.

Аватара пользователя
хуссейн макассис
Сообщения: 14
Зарегистрирован: 28 ноя 2017, 15:59

Re: Целочисленное проектирование фильтров

Сообщение хуссейн макассис » 10 окт 2018, 10:45

What is a minimum phase filter? Is it possible to make a low-bit FIR filter minimal phase? :?:

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 772
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: Целочисленное проектирование фильтров

Сообщение Бахурин Сергей » 10 окт 2018, 11:10

Follow paper describes minimum phase filters

https://www.dsprelated.com/freebooks/fi ... lters.html

You can try to design minimum phase low bit filter too.

Аватара пользователя
хуссейн макассис
Сообщения: 14
Зарегистрирован: 28 ноя 2017, 15:59

Re: Целочисленное проектирование фильтров

Сообщение хуссейн макассис » 28 окт 2018, 14:24

Есть рассчитанные коэффициенты КИХ-фильтра. Как по ним определить и построить его графики АЧХ и ФЧХ? Можно ссылку если это где-то разжевано.

ЗЫ матлаб навскидку скачать не удалось, есть GNU Octave (с интерфейсом командной строки), GNU Plot.

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 772
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: Целочисленное проектирование фильтров

Сообщение Бахурин Сергей » 28 окт 2018, 15:37

Передаточная функция КИХ фильтра порядка



Комплексный коэффициент передачи:



Тогда АЧХ это модуль от , а ФЧХ - аргумент. В octave это выглядить так:

Код: Выделить всё

% b - FIR filter coeff
w = -pi:0.001:pi;
H = freqz(b,1,w);
MAG = abs(H); % Magnitude
PHI = angle(H) % Phase response



Ответить

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость