Яка найгостріша частотна характеристика для позапричинного низькочастотного фільтра, ступінчаста відповідь якого не перевищує?


13

Фільтри низького пропускання Баттерворта, Бесселя, Чебичева та sinc використовуються в різних випадках, коли є різні компроміси між тим, що мають рівномірно зменшувану частотну характеристику, рівномірний фазовий відбиток, крутий обрив або реакцію "цегляна стіна". Я вважаю, що у всіх таких фільтрах можуть в деяких випадках спостерігатися перехід на ступінчасту реакцію, тобто імпульсна реакція в деяких місцях негативна.

Яка була б оптимальна частотна характеристика або які типи частотних відповідей були б доступні у фільтрі, єдиним обмеженням якого було те, що імпульсна характеристика не може бути негативною ні в якому разі? Звичайно, можливо, щоб фільтр низьких частот відповідав такому обмеженню, оскільки базовий RC-фільтр це зробить (хоча реакція такого фільтра трохи крихка). Оптимальна реакція імпульсу - це нормальна крива розподілу чи щось інше?


1
@supercat, якщо включити цифрову фільтрацію, дивно, як відповідь цегляної стіни ви можете отримати без промахування.
Кортук

2
@Kortuk: Дійсно? Я б подумав, що важко уникнути прострілу, оскільки квадратна хвиля, відфільтрована цегляною стіною, має невеликі шипи, ширина яких наближається до нуля, коли частота обрізання зростає, але амплітуда якої - ні. Яка б була хороша довідка?
supercat

2
Ви говорите "непричинно" у питанні, але всі ваші приклади є причинними. Що ти маєш на увазі? Безпричини вимагають записати всю форму хвилі, а потім застосувати фільтр до запису. (Або, можливо, використання флюсових конденсаторів та великих джерел живлення.)
ендоліт

1
@endolith: Що було б оптимальним фільтром, якщо припустити, що це не обов'язково, щоб бути причинним.
supercat

1
@Kortuk: Відсікання сигналу в нулі повністю знищить будь-які переваги фільтрації. І поки я публікую повідомлення в DSP, мені також цікаво такі речі, як магнітофони (з інтелектуального, а не практичного аспекту), де можна робити будь-яку негативну функцію імпульсу, за умови обмеження ширини.
supercat

Відповіді:


6

Я збираюся перелічити купу "фільтрів, які не перекривають". Я сподіваюся, що ви знайдете цю часткову відповідь краще, ніж взагалі ніякої відповіді. Сподіваємось, люди, які шукають "фільтр, який не перекриває", вважають цей список таких фільтрів корисним. Можливо, один із цих фільтрів буде працювати належним чином у вашій програмі, навіть якщо ми ще не знайшли математично оптимального фільтра.

Причинно-наслідкові фільтри першого та другого порядку

Етапна відповідь фільтра першого порядку ("RC фільтр") ніколи не перевищує.

Етапна відповідь фільтра другого порядку ("біквад") може бути спроектована таким чином, що він ніколи не перестарається. Існує кілька рівнозначних способів опису цього класу фільтра другого порядку, який не перевершує на крок введення:

  • вона критично затухає або перекривається.
  • це не завалюється.
  • коефіцієнт демпфування (зета) становить 1 або більше
  • коефіцієнт якості (Q) становить 1/2 або менше
  • параметр швидкості занепаду (альфа) є щонайменше неврівноваженою природною кутовою частотою (омега_0) або більше

Зокрема, критично демпфірується коефіцієнт посилення єдності фільтра Саллена – Клю з рівними конденсаторами та рівними резисторами: Q = 1/2, а отже, не перевищує ступінчастий вхід.

Фільтр Бесселя другого порядку трохи недоотриманий: Q = 1 / sqrt (3), тому він має невеликий переріз.

Фільтр Баттерворта другого порядку більш занепокоєний: Q = 1 / sqrt (2), тому він має більше промахів.

З усіх можливих фільтрів LTI першого порядку та другого порядку, які є причинними та не перевищують, той, що має "найкращу" (найбільш круту) частотну характеристику, є "критично затухаючими" фільтрами другого порядку.

Причинно-наслідкові фільтри вищого порядку

Найбільш часто використовуваний причинний фільтр вищого порядку, який має імпульсну характеристику, яка ніколи не є негативною (і, отже, ніколи не переступає на ступінчастий вхід) - це "фільтр середнього бігу", який також називають "фільтром боксерського автомобіля" або " фільтром ковзних середніх" ".

Деякі люди люблять запускати дані через один фільтр коробки, а вихід з цього фільтра в інший фільтр боксерського вагона. Після кількох таких фільтрів, результат є гарним наближенням фільтра Гаусса. (Чим більше фільтрів ви каскадуєте, тим ближче кінцевий вихід наближається до Гаусса, незалежно від того, з якого фільтра ви починаєте - боксер, трикутник, RC першого порядку або будь-який інший - через центральну граничну теорему).

Практично всі віконні функції мають імпульсну характеристику, яка ніколи не є негативною, і тому в принципі може використовуватися як фільтри FIR, які ніколи не запускають на крок вхід. Зокрема, я чую хороші речі про вікно Ланцоса , що є центральною (позитивною) долею функції sinc () (і нуль поза цією часточкою). Кілька фільтрів, що формують імпульс, мають імпульсну характеристику, яка ніколи не є негативною, і тому їх можна використовувати як фільтри, які ніколи не переступають на ступінчастий вхід.

Я не знаю, який із цих фільтрів найкращий для вашої програми, і я підозрюю, що математично оптимальний фільтр може бути дещо кращим, ніж будь-який із них.

нелінійні причинно-наслідкові фільтри

Медіанний фільтр є популярним нелінійним фільтром , який ніколи не проскакує на вході ступінчастої функції.

EDIT: безаварійні фільтри LTI

Функція sech (t) = 2 / (e ^ (- t) + e ^ t) - це його власне перетворення Фур'є, і я вважаю, що він може бути використаний як своєрідний некоузальний низькочастотний фільтр LTI, який ніколи не переступає на a крок введення.

Непричинний фільтр LTI, який має імпульсну відповідь (sinc (t / k)) ^ 2, має частотну характеристику "abs (k) * трикутник (k * w)". Якщо дано ступінчастий вхід, у нього багато пульсацій часової області, але він ніколи не перекриває остаточну точку відстоювання. Вище високочастотного кута цього трикутника він дає ідеальне відхилення смуги зупинки (нескінченне ослаблення). Таким чином, в області стоп-діапазону він має кращу частотну характеристику, ніж фільтр Гаусса.

Тому я сумніваюся, що фільтр Гаусса дає "оптимальну частотну характеристику".

У наборі всіх можливих «фільтрів, які не перекривають», я підозрюю, що немає жодного єдиного «оптимальної частотної характеристики» - деякі мають кращу відмову від зупинки, а інші мають більш вузькі смуги переходу тощо.


Спасибі за Вашу відповідь. Я нехтував обмеженням свого питання лише лінійними фільтрами, хоча, звичайно, характеристика частотної характеристики нелінійного фільтра може бути мутною пропозицією. Як зазначалося, каскадний фільтр боксерів змушує його наближатися до гаусса. Мені було цікаво, чи має фільтр Гаусса оптимальну частотну характеристику, яку можна отримати без прокрутки. Пишучи питання, я думав про різні аналогові процеси, які виконують щось на кшталт визначеного імпульсного фільтра, наприклад, розмивання камери або пікселів дисплея, щоб мінімізувати згладжування.
supercat

Можна створити камеру так, щоб кожен піксель набирав різну кількість світла з різних точок по центру. В ідеалі, камера відфільтрувала б усе вище Nyquist, не розмиваючи нічого внизу, але на практиці це, мабуть, не відбудеться.
supercat

1
LTI? Ви ніколи цього не визначаєте. Додавання, що це означає "лінійна інваріантність у часі", можливо, буде корисним.
Коннор Вольф

1
Тож Q = 0,5 критично затухає? Чи є для даного замовлення кілька систем, які критично затухають? Біквад з Q = 0,5 називається фільтром LR2 Linkwitz-Riley . Схоже, версії LR-фільтра вищого порядку мають дзвінок у відповіді на ступінь.
ендоліт

"каскадні фільтри, що критично демпфірують, призведуть до появи іншого фільтру з критичним демпфіруванням". Так що просто тримайте демпінгові полюси на рівні -1, і це завжди буде критично затухати? (І підходити до реакції фільтра Гаусса, коли число збільшується?)
ендоліт

2

Більшість фільтрів, що використовуються в цифровому світі, є лише вибірковою версією аналогового аналога. Великою причиною цього є те, що перед аналогічним фільтруванням було зроблено багато роботи, перш ніж з'явилася цифра, тож замість того, щоб винаходити колесо, більшість щойно використовуваних попередніх конструкцій. Перевага цифрового в тому, що фільтр вищого порядку можна досягти набагато простіше, ніж в аналоговому світі. Уявіть собі складну схему, отримуючи щоразу, коли ви додаєте чергове замовлення до дизайну.

Якщо ви збираєтесь фільтрувати тип цегляної стіни, крива Гаусса - непогане місце для початку. Якщо ви знаєте про часовий домен <-> частотний домен; Гаусс перетворюється на гаусса в іншій області. Оскільки в одній стає вітреше, в іншому стає вужчою. Отже, щоб отримати ідеальний сплеск частотної області, вам знадобиться нескінченна кількість зразків.

Якщо у вас є Matlab доступний для використання, вам слід ознайомитися з деякими вбудованими інструментами дизайну фільтрів. Ось посилання, що говорить про Баттерворта та Бесселя . Засоби проектування дозволяють задати певні аспекти фільтра. Ці аспекти змінюються для кожного типу фільтра, але деякі приклади - пропускна здатність, зупинка, пульсація тощо. Якщо ви надаєте дизайнеру потрібні обмеження, він або видасть помилку (тобто він не може зробити цей фільтр із таким типом фільтра ) або він дасть вам фільтр з мінімальним замовленням, необхідним для виконання специфікацій.


Gaussian відповідає вимогам для неотрицательного імпульсної характеристики, але це не велика частина цегляної стіни. З іншого боку, Баттерворт, Бессель і Чебишев мають імпульсні реакції, подібні до синхронних коливань, що призводять до перегину. З них фільтр Бесселя має найменший пробіг, оскільки він має майже рівну групову затримку (лінійна фаза) в смузі пропускання.
Ерік Вс

Окрім гауссових, ці фільтри є причинними. Для офлайн-обробки лінійна фаза NNFIR (негативний FIR) буде добре працювати, або ви можете скасувати фазове спотворення причинно-рекурсивного фільтра, використовуючи техніку filtfilt .... Звичайно, вам все ще потрібен спосіб проектування LPF NNIR для уникнути перекриття / підкреслення. Пропозиції кому? Список літератури?
Ерик ВС

@eryksun, як бічна примітка, якщо він буде коливатися в 1,05 рази більше значення максимуму, просто зволожте, щоб зупинитися на 1,00, і ваша крокова реакція буде трохи меншою, як-от .96, коли стабільна. Проблема вирішена.
Кортук

@Kortuk: Проблема вирішена у часовій області, але виконання цього відсікання є не тільки нелінійним, а й відкриває певні частини частотної області для передачі сигналу, який раніше не був. Він хоче максимально жорсткий фільтр пропуску / без пропуску за частотою без прокрутки в часовій області. Немає перевищення часового домену - це те саме, що сказати, що імпульсна відповідь ніколи не є негативною.
Олін Латроп

1
@Kortuk: У деяких областях нуль знаходиться поблизу середньої точки між максимумом і мінімумом, і масштабування сигналу до середини не дозволить уникнути проблем із простроченням. В інших областях, таких як візуалізація, нуль - це мінімум, а динамічний діапазон є найважливішим там; як правило, краще мати «нечіткий» фільтр, який не проганяє, ніж гостріший.
supercat
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.