Тема появилась благодаря пользователю Vlad27, предоставившему ссылки на ряд интересных статей:
Попытаемся оценить, насколько хороши предложенные фильтры.
Пока опираться буду на статью "ctat_3.pdf", где подробно описаны характеристики фильтра.
Для начала я планирую посмотреть, что собой представляют "звенья" этого фильтра.
Для каждого звена посмотрим H(z).
Потом "тупо" умножим полиномы и посмотрим, что из себя представляет полученный фильтр 16 порядка...
Может он ничем от Кауэра 10 порядка ( с округлёнными коэффициентами) и не отличается?
(Это моё такое скептическое предположение ) Надеюсь, что я не прав!
Целочисленное проектирование фильтров
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Ну вот, начало положено:
О, как интересно - 1 звено - "расходящийся фильтр" - Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Вот характеристики фильтра по звеньям:
Интересно, что "фильтром" на частоту близкую к 0.5 являются только звенья с номерами 1,3,8.
Значит остальные 5 звеньев "заняты" линеаризацией ФЧХ ?
Надо будет еще ФЧХ построить....
Результат перемножения H1(z)...H8(z):Интересно, что "фильтром" на частоту близкую к 0.5 являются только звенья с номерами 1,3,8.
Значит остальные 5 звеньев "заняты" линеаризацией ФЧХ ?
Надо будет еще ФЧХ построить....
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
И вот сейчас появились вопросы к Vlad27:
1. Хотелось бы увидеть фильтр не на 0.5
Дело в том, что частота 0.5 - имеет некий "сакральный" смысл... И все нечётные коэффициенты "классических" полосовых фильтров (Баттерворта, Чебышева, Кауэра) на этой частоте равны 0. Т.е. число "реальных" коэффициентов почти в 2 раза меньше.
2. Хотелось бы увидеть фильтр с полосой, хотя бы 0.2 (и на центральную частоту не равную 0.5)
1. Хотелось бы увидеть фильтр не на 0.5
Дело в том, что частота 0.5 - имеет некий "сакральный" смысл... И все нечётные коэффициенты "классических" полосовых фильтров (Баттерворта, Чебышева, Кауэра) на этой частоте равны 0. Т.е. число "реальных" коэффициентов почти в 2 раза меньше.
2. Хотелось бы увидеть фильтр с полосой, хотя бы 0.2 (и на центральную частоту не равную 0.5)
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Конечно интересно сравнить ЦНП-фильтр с другими известными типами фильтров, ну например:
Фильтр 1:
b0=1. b1=0. b2=0.939
a0=6.98E+5 a1=0. a2=6.98E+5
Фильтр 1:
b0=1. b1=0. b2=0.939
a0=6.98E+5 a1=0. a2=6.98E+5
-
- Сообщения: 90
- Зарегистрирован: 06 апр 2015, 10:09
- Откуда: г, Нижний Новгород (rif27@bk.ru)
Re: Целочисленное проектирование фильтров
Конечно, Александр, гауссов фильтр можно синтезировать на любую
частоту, а не только 0,5. При 0,5 можно перейти к фильтру Бесселя
без боковых лепестков вообще - синтезировать его гораздо проще.
Когда я решал эту залдачу (довольно давно) серединка была выбрана
случайно, по умолчанию что-ли. Но сужать полосу фильтра уже не удаётся,
даже более высоким порядком фильтра. Сейчас получено решение
двухполосного гауссова фильтра, позволяющего осуществлять гауссову
фильтрацию сразу двух сигналов.
частоту, а не только 0,5. При 0,5 можно перейти к фильтру Бесселя
без боковых лепестков вообще - синтезировать его гораздо проще.
Когда я решал эту залдачу (довольно давно) серединка была выбрана
случайно, по умолчанию что-ли. Но сужать полосу фильтра уже не удаётся,
даже более высоким порядком фильтра. Сейчас получено решение
двухполосного гауссова фильтра, позволяющего осуществлять гауссову
фильтрацию сразу двух сигналов.
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Это понятно. Вопрос в том , если с теми же параметрами фильтр создать не на частоту 0.5, а на 0.3 у него порядок 16 сохраниться? Или порядок получится более высокий?Vlad27 писал(а):Конечно, Александр, гауссов фильтр можно синтезировать на любую
частоту, а не только 0,5....
-
- Сообщения: 90
- Зарегистрирован: 06 апр 2015, 10:09
- Откуда: г, Нижний Новгород (rif27@bk.ru)
Re: Целочисленное проектирование фильтров
В ЦЦФ, как и в ЦФ ведь всё наоборот по сравнению с аналоговыми фильтрами.
В аналог.фильтрах - чем выше частота, тем сложнее получить хорошую селективность
(прямоугольность АЧХ) . А в ЦЦФ - чем меньше угол (цифровая частота), тем сложнее.
То есть при той же ширине полосы пропускания гауссова фильтра, его легче реализовать
на более высоких цифровых частотах, чем на малых углах (или с меньшим порядком).
А сдвинуть можно куда угодно, главное подавить боковые лепестки АЧХ.
В аналог.фильтрах - чем выше частота, тем сложнее получить хорошую селективность
(прямоугольность АЧХ) . А в ЦЦФ - чем меньше угол (цифровая частота), тем сложнее.
То есть при той же ширине полосы пропускания гауссова фильтра, его легче реализовать
на более высоких цифровых частотах, чем на малых углах (или с меньшим порядком).
А сдвинуть можно куда угодно, главное подавить боковые лепестки АЧХ.
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Да я понял теперь. Всё равно Ваша H(z) - это произведение Hn(z) - сдвинутых относительно 0.5, причем ни одна на частоту 0.5 не попадает (все а1, b1 не нулевые)!
Было бы интересно для фильтра 0.5 "насильно" задать в одном из звеньев а1, b1 =0 и оптимизировать все остальные коэффициенты.
Было бы интересно для фильтра 0.5 "насильно" задать в одном из звеньев а1, b1 =0 и оптимизировать все остальные коэффициенты.
- Santik
- Сообщения: 609
- Зарегистрирован: 28 дек 2010, 08:04
- Откуда: Мирный (Якутия)
- Контактная информация:
Re: Целочисленное проектирование фильтров
Продолжаем "рассматривать" фильтр: