Jaka jest najostrzejsza odpowiedź częstotliwościowa dla niepochodzącego z przyczyn filtra dolnoprzepustowego, którego odpowiedź skokowa nie jest zbyt wysoka?


13

Filtry dolnoprzepustowe Butterwortha, Bessela, Czebycheva i cynku stosuje się w różnych przypadkach, w których występują różne kompromisy między posiadaniem jednorodnie malejącej odpowiedzi częstotliwościowej, jednolitej odpowiedzi fazowej, stromego odcięcia lub odpowiedzi „ściany z cegieł”. Wierzę, że wszystkie takie filtry mogą, w niektórych przypadkach, przekroczyć swoją odpowiedź krokową, co oznacza, że ​​ich odpowiedź impulsowa jest w niektórych miejscach negatywna.

Jaka byłaby optymalna odpowiedź częstotliwościowa lub jakie typy odpowiedzi częstotliwości byłyby dostępne w filtrze, którego jedynym ograniczeniem było to, że odpowiedź impulsowa nie może być nigdzie ujemna? Z pewnością możliwe jest, aby filtr dolnoprzepustowy spełniał takie ograniczenie, ponieważ zrobi to podstawowy filtr RC (choć odpowiedź takiego filtra jest nieco kiepska). Czy optymalna odpowiedź impulsowa byłaby krzywą rozkładu normalnego, czy coś innego?


1
@ superupat, jeśli włączysz filtrowanie cyfrowe, to zadziwiające, jak ceglana ściana może uzyskać odpowiedź bez przekroczenia.
Kortuk

2
@Kortuk: Naprawdę? Wydaje mi się, że trudno byłoby uniknąć przekroczenia, ponieważ kwadratowa fala o przefiltrowanej ścianie z cegieł ma małe skoki, których szerokość zbliża się do zera wraz ze wzrostem częstotliwości odcięcia, ale której amplituda nie. Jakie byłoby dobre odniesienie?
supercat

2
Mówisz „nie przyczynowo” w pytaniu, ale wszystkie twoje przykłady są przyczynowe. Co masz na myśli W przypadku braku związku przyczynowego należy zarejestrować cały przebieg, a następnie zastosować filtr do nagrania. (Lub, być może, zastosowanie kondensatorów topnikowych i dużych źródeł zasilania.)
endolith

1
@endolith: Jaki byłby optymalny filtr, zakładając, że nie musi być przyczynowy.
supercat,

1
@Kortuk: Obcięcie sygnału na zero całkowicie zniweczy wszelkie korzyści z filtrowania. Podczas gdy publikuję w DSP, ciekawi mnie także takie rzeczy, jak nagrywarki filmów audio (co prawda z intelektualnego, a nie praktycznego punktu widzenia), w których można zastosować dowolną nieujemną funkcję impulsu, którą chcemy, z zastrzeżeniem ograniczenia szerokości.
supercat

Odpowiedzi:


6

Mam zamiar wyświetlić listę „filtrów, które nie przekraczają”. Mam nadzieję, że ta częściowa odpowiedź będzie dla ciebie lepsza niż brak odpowiedzi. Mam nadzieję, że osobom szukającym „filtra, który nie przereguluje”, ta lista takich filtrów będzie pomocna. Być może jeden z tych filtrów będzie działał odpowiednio w Twojej aplikacji, nawet jeśli nie znaleźliśmy jeszcze matematycznie optymalnego filtra.

filtry przyczynowe pierwszego i drugiego rzędu LTI

Odpowiedź krokowa filtra pierwszego rzędu („filtr RC”) nigdy nie jest zbyt wysoka.

Reakcja krokowa filtra drugiego rzędu („biquad”) może być zaprojektowana w taki sposób, aby nigdy nie przekraczała. Istnieje kilka równoważnych sposobów opisania tej klasy filtra drugiego rzędu, który nie jest przekroczony po wprowadzeniu kroku:

  • jest krytycznie tłumiony lub nadmiernie tłumiony.
  • nie jest zbyt wilgotne.
  • współczynnik tłumienia (zeta) wynosi 1 lub więcej
  • współczynnik jakości (Q) wynosi 1/2 lub mniej
  • parametr szybkości rozpadu (alfa) jest co najmniej nie tłumioną naturalną częstotliwością kątową (omega_0) lub więcej

W szczególności topologia filtra Sallen-Key z równymi kondensatorami i jednakowymi rezystorami jest krytycznie tłumiona: Q = 1/2, a zatem nie przekracza wartości wejściowej kroku.

Filtr Bessela drugiego rzędu jest nieco niedostatecznie tłumiony: Q = 1 / sqrt (3), więc ma niewielkie przekroczenie.

Filtr Butterwortha drugiego rzędu jest bardziej niedokładny: Q = 1 / sqrt (2), więc ma większe przekroczenie.

Ze wszystkich możliwych filtrów LTI pierwszego i drugiego rzędu, które są przyczynowe i nie przekraczają, te z „najlepszą” (stromą) odpowiedzią częstotliwościową są „krytycznie tłumione” filtry drugiego rzędu.

filtry przyczynowe wyższego rzędu LTI

Najczęściej stosowanym filtrem przyczynowym wyższego rzędu, który ma odpowiedź impulsową, która nigdy nie jest ujemna (a zatem nigdy nie przekracza parametrów wejściowych kroku), jest „filtr średniej bieżącej”, zwany również „filtrem Boxcara” lub „ filtrem średniej ruchomej ” „.

Niektórzy ludzie lubią przepuszczać dane przez jeden filtr Boxcar, a dane wyjściowe z tego filtra do innego filtra Boxcar. Po kilku takich filtrach wynikiem jest dobre przybliżenie filtra Gaussa. (Im więcej filtrów kaskadujesz, tym bliższe końcowe wyjście jest w przybliżeniu gaussowskie, bez względu na filtr, z którego zaczynasz - boxcar, trójkąt, RC pierwszego rzędu lub inne - ze względu na centralne twierdzenie graniczne).

Praktycznie wszystkie funkcje okna mają odpowiedź impulsową, która nigdy nie jest ujemna, a zatem w zasadzie mogą być używane jako filtry FIR, które nigdy nie przekraczają wartości wejściowej kroku. W szczególności dobrze słyszę o oknie Lanczosa , które jest środkowym (dodatnim) płatem funkcji sinc () (i zerowym poza tym płatem). Kilka filtrów kształtujących impulsy ma odpowiedź impulsową, która nigdy nie jest ujemna, a zatem może być stosowana jako filtry, które nigdy nie przekraczają wartości wejściowej kroku.

Nie wiem, który z tych filtrów jest najlepszy dla twojej aplikacji i podejrzewam, że optymalny matematycznie filtr może być nieco lepszy niż którykolwiek z nich.

nieliniowe filtry przyczynowe

Filtr mediana to popularny filtr nieliniowy, że nigdy przeregulowania na wejście krok funkcyjnego.

EDYCJA: filtry przyczynowe LTI

Funkcja sech (t) = 2 / (e ^ (- t) + e ^ t) jest własną transformacją Fouriera, i przypuszczam, że mogłaby zostać użyta jako rodzaj bez przyczynowego dolnoprzepustowego filtru LTI, który nigdy nie jest zbyt wysoki wejście krokowe.

Bez przyczynowy filtr LTI, który ma odpowiedź impulsową (sinc (t / k)) ^ 2, ma odpowiedź częstotliwościową „abs (k) * trójkąt (k * w)”. Po wprowadzeniu danych krokowych ma wiele tętnień w dziedzinie czasu, ale nigdy nie przekracza ostatecznego punktu osiadania. Powyżej narożnika tego trójkąta o wysokiej częstotliwości zapewnia idealne odrzucenie pasma zatrzymania (nieskończone tłumienie). Zatem w obszarze pasma zatrzymania ma lepszą odpowiedź częstotliwościową niż filtr Gaussa.

Dlatego wątpię, aby filtr Gaussa zapewniał „optymalną odpowiedź częstotliwościową”.

Podejrzewam, że w zestawie wszystkich możliwych „filtrów, które nie przekraczają”, nie ma jednego „optymalnego pasma przenoszenia” - niektóre mają lepsze odrzucanie pasma zatrzymania, podczas gdy inne mają węższe pasma przejściowe itp.


Dzięki za odpowiedź. Zlekceważyłem ograniczanie mojego pytania do filtrów liniowych, chociaż oczywiście charakteryzowanie odpowiedzi częstotliwościowej filtra nieliniowego może być mętną propozycją. Jak wspomniano, kaskadowanie filtra boxcar powoduje, że zbliża się on do gaussowskiego. Zastanawiałem się, czy filtr Gaussa ma optymalną charakterystykę częstotliwościową, którą można uzyskać bez przekroczenia. Pisząc pytanie, myślałem o różnych procesach analogowych, które wykonują coś w rodzaju określonego filtra odpowiedzi impulsowej, np. Rozmycie kamery lub wyświetlenie pikseli w celu zminimalizowania aliasingu.
supercat

Możliwe jest zbudowanie kamery, aby każdy piksel wychwytywał różne ilości światła z różnych punktów wokół środka. Najlepiej byłoby, gdyby kamera odfiltrowała wszystko powyżej Nyquist bez rozmycia niczego poniżej, ale w praktyce tak się nie stanie.
supercat

1
LTI? Nigdy tego nie definiujesz. Dodanie, że oznacza to „liniowy niezmiennik czasowy”, prawdopodobnie byłoby pomocne.
Connor Wolf,

1
Więc Q = 0,5 jest krytycznie tłumione? Czy dla danego zamówienia istnieje wiele systemów, które są krytycznie tłumione? Biquad o Q = 0,5 nazywa się filtrem LR2 Linkwitz-Riley . Wygląda jednak na to, że wersje filtra LR wyższego rzędu dzwonią w odpowiedzi krokowej.
endolith

„Kaskadowanie filtrów z tłumieniem krytycznym spowoduje powstanie kolejnego filtra z tłumieniem krytycznym” Więc po prostu utrzymuj bieguny zrzutu na poziomie -1 i zawsze będzie ono tłumione krytycznie? (I zbliża się odpowiedź filtra Gaussa, gdy liczba rośnie?)
endolith

2

Większość filtrów używanych w świecie cyfrowym to tylko próbkowana wersja analogowego odpowiednika. Głównym powodem tego jest to, że przed wprowadzeniem technologii cyfrowej wykonano wiele pracy w zakresie filtrowania analogowego, więc zamiast wynaleźć koło, większość z nich była używana wcześniej. Zaletą technologii cyfrowej jest to, że filtr wyższego rzędu można osiągnąć znacznie łatwiej niż w świecie analogowym. Wyobraź sobie, jak skomplikowany obwód dostajesz za każdym razem, gdy dodasz kolejne zamówienie do projektu.

Jeśli wybierasz filtr typu ceglanego muru, krzywa Gaussa jest całkiem dobrym miejscem do rozpoczęcia. Jeśli wiesz o domenie czasu <-> domena częstotliwości; Gaussian przekształca się w Gaussa w drugiej domenie. Gdy robi się coraz bardziej nawijający w jednym, w drugim staje się węższy. Aby uzyskać doskonały skok w dziedzinie częstotliwości, potrzebujesz nieskończonej ilości próbek.

Jeśli akurat masz dostęp do Matlaba, powinieneś sprawdzić niektóre z wbudowanych narzędzi do projektowania filtrów. Oto link mówiący o Butterworth i Bessel . Narzędzia do projektowania pozwalają określić pewne aspekty filtra. Te aspekty zmieniają się dla każdego typu filtra, ale niektóre przykłady to Pasmo przepustowe, pasmo, tętnienie itp. Jeśli podasz projektantowi wymagane ograniczenia, albo spowoduje błąd (co oznacza, że ​​nie może utworzyć tego filtru z tym typem filtra) ) lub da Ci filtr o minimalnej kolejności wymaganej do spełnienia specyfikacji.


Gaussa pasuje wymóg nieujemną odpowiedzi impulsowej, ale to nie jest dużo murem. Z drugiej strony Butterworth, Bessel i Czebyshev mają drgania podobne do cynkowych w odpowiedzi impulsowej, co prowadzi do przeregulowania. Spośród nich filtr Bessela ma najmniejsze przekroczenie, ponieważ ma prawie płaskie opóźnienie grupowe (faza liniowa) w paśmie przepustowym.
Eryk Sun

Oprócz filtrów gaussowskich filtry te są przyczynowe. Do przetwarzania w trybie offline dobrze działałaby liniowa faza NNFIR (nieujemna FIR), lub można anulować zniekształcenie fazowe przyczynowego filtra rekurencyjnego za pomocą techniki filtrowania ... Oczywiście nadal potrzebujesz sposobu zaprojektowania NNIR LPF do unikać przeregulowania / niedopracowania. Sugestie ktoś? Bibliografia?
Eryk Sun

@eryksun, na marginesie, jeśli ma oscylować przy 1,05 razy wartość max, po prostu wytłum go, aby zatrzymać się na 1,00, a twoja odpowiedź krokowa będzie nieco mniejsza, jak .96, gdy stabilny. Problem rozwiązany.
Kortuk

@Kortuk: Problem rozwiązany w dziedzinie czasu, ale robienie tego obcinania jest nie tylko nieliniowe, ale także otwiera pewne części dziedziny częstotliwości, aby przekazać sygnał, który wcześniej nie był. Chce jak najściślejszego filtru pass / no pass w częstotliwości bez przekroczenia w dziedzinie czasu. Żadne przekroczenie dziedziny czasu nie jest tym samym, co stwierdzenie, że odpowiedź impulsowa nigdy nie jest negatywna.
Olin Lathrop,

1
@Kortuk: W niektórych domenach zero znajduje się w pobliżu punktu środkowego między maksimum a minimum, a skalowanie sygnału w kierunku punktu środkowego pozwoli uniknąć problemów z przeregulowaniem. W innych dziedzinach, takich jak obrazowanie, zero jest minimum, a zakres dynamiczny jest tam najważniejszy; ogólnie lepiej byłoby mieć filtr „bardziej zamglony”, który nie jest wyregulowany, niż filtr ostrzejszy.
supercat
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.