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

Все что касается фильтрации
Аватара пользователя
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
Сообщения: 71
Зарегистрирован: 06 апр 2015, 10:09
Откуда: г, Нижний Новгород (rif27@bk.ru)

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

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

Задача по доплеровской обработке локационных сигналов от
быстродвижущихся объектов определяет весьма жёсткие требования
к целочисленному цифровому обрабатывающему фильтру - своеобразный
гибрид гауссова фильтра и ФВЧ или узкополосный гауссовый БИХ-фильтр на
"пьедестале" ФВЧ. Вот его ЦНП-синтез и результаты измерения
в частотной области.
Первые два графика - решение по синтезу, вторые два - измерение АЧХ
этого доплеровского фильтра на реальном сигнале.
Вложения
gain_all.gif
АЧХ доплеровского фильтра по синтезу
gain_all.gif (7.77 КБ) 1034 просмотра
gain_all_dB_2.gif
ЦНП-решение, дБ
gain_all_dB_2.gif (7.64 КБ) 1034 просмотра
панорама_АЧХ_нормир_1.gif
Измерение АЧХ
панорама_АЧХ_нормир_1.gif (9.8 КБ) 1034 просмотра
панорама_АЧХ_нормир_дБ_Fн.gif
Измерения АЧХ в Дб

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

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

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

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

Ответить

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

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

Создано на основе phpBB® Forum Software © phpBB Limited Русская поддержка phpBB