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

Все что касается фильтрации
Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

Коэффициенты b:
9.985882797899660000e-001
-1.997176559578946400e+000
9.985882797899660000e-001
Коэффициенты a:
1.000000000000000000e+000
-1.997176559578946400e+000
9.971765595799320100e-001
Но как из них сделать целые я не знаю :-) На 4096 умножить?

trofim2
Сообщения: 14
Зарегистрирован: 24 апр 2015, 13:44
Откуда: Гродно By

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

Сообщение trofim2 »

Ув. Santik.
А Вы не перепутали к-ты А и B?
Обычно В0 = 1.

trofim2
Сообщения: 14
Зарегистрирован: 24 апр 2015, 13:44
Откуда: Гродно By

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

Сообщение trofim2 »

У меня такой фильтр возбуждается.

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

trofim2 писал(а): А Вы не перепутали к-ты А и B?
В обозначениях, принятых на dsplib

Соответственно:
Y[n]=1/a_0{b_0X[n]+b_1X[n-1+b_2X[n-2]- a_1Y[n-1]-a_2Y[n-2]}

trofim2
Сообщения: 14
Зарегистрирован: 24 апр 2015, 13:44
Откуда: Гродно By

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

Сообщение trofim2 »

Я примерно так и понял.
При частоте среза свыше 50 Гц (при частоте квантования 3500 Гц) фильтр работает нормально, а при понижении начинаются непонятки. Коэффициенты В я принимаю равными 1,2 и 1. К-ты А умножаю на 256 и округляю.
Наверное не хватает точности.
Остановлюсь на предыдущем варианте.

В последней формуле при к-те а_2 Y[n-2]?

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

trofim2 писал(а): В последней формуле при к-те а_2 Y[n-2]?
Да. Исправил.

Да, на частотах среза близких к 0 и частоте Найквиста - представление коэффициентов - float, желательно...

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

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

Сообщение Vlad27 »

Задача по доплеровской обработке локационных сигналов от
быстродвижущихся объектов определяет весьма жёсткие требования
к целочисленному цифровому обрабатывающему фильтру - своеобразный
гибрид гауссова фильтра и ФВЧ или узкополосный гауссовый БИХ-фильтр на
"пьедестале" ФВЧ. Вот его ЦНП-синтез и результаты измерения
в частотной области.
График - решение по синтезу
Вложения
АЧХ доплеровского фильтра по синтезу
АЧХ доплеровского фильтра по синтезу
gain_all.gif (7.77 КБ) 7667 просмотров
Последний раз редактировалось Vlad27 06 ноя 2018, 11:28, всего редактировалось 2 раза.

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

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

Сообщение Vlad27 »

С новыми возможностями проектирования цифровых фильтров со сложными селективными требованиями в
целочисленном пространстве параметров малой разрядности можно ознакомиться
в книге "В.Н.Бугров. Целочисленные цифровые фильтры" издательства LAMBERT (ФРГ) 2017 года, в открытом доступе on-line.
Целочисленные_ЦФ.rar
(6.5 МБ) 189 скачиваний

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

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

Сообщение Vlad27 »

Вопрос к Сергею Бахурину - как специалисту по МАТЛАБу

Методом ЦНП рассчитан AllPass 8-го порядка, рекурсивный - его коэффициенты на рисунке
Не могли бы вы, Сергей, просчитать его характеристики в МАЛЛАБе или другом альт.пакете,
просто просчитать, не синтезировать. У нас в ЦНП-пакете удивительные характеристики
получаются. Я их позже выложу, чтобы вас не смущать. Заранее благодарен
Вложения
Forum_AllPass.gif
Forum_AllPass.gif (3.16 КБ) 5543 просмотра

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

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

Сообщение Бахурин Сергей »

Посмотрел. Ваш фильтр выродился. Числитель его H(z) практически стал равен знаменателю и он практически никак не изменяет сигнал. Реально амплитудные колебания сотые доли дБ и фаза поворачивается на доли градуса. Импульсная характеристика очень близка к дельта импульсу.

Ответить