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

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

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

Сообщение Santik » 18 фев 2016, 18:47

Коэффициенты 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 » 19 фев 2016, 09:04

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

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

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

Сообщение trofim2 » 19 фев 2016, 09:35

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

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

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

Сообщение Santik » 19 фев 2016, 10:10

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 » 19 фев 2016, 11:05

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

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

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

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

Сообщение Santik » 19 фев 2016, 11:22

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

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

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

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

Сообщение Vlad27 » 18 апр 2016, 11:29

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

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

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

Сообщение Vlad27 » 13 мар 2017, 12:02

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

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

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

Сообщение Vlad27 » 18 сен 2018, 12:42

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

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

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

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

Сообщение Бахурин Сергей » 18 сен 2018, 13:15

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

Ответить

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

Сейчас этот форум просматривают: Алексей Климов и 1 гость