Struktura sprzężenia zwrotnego
Filtr FIR (skończona odpowiedź impulsowa) oblicza każdą próbkę wyjścia jako ważoną sumę obecnych i przeszłych wejść tylko. Bez sprzężenia zwrotnego. Odpowiedź impulsowa ma skończony czas trwania.
Filtr IIR (nieskończona odpowiedź impulsowa) wprowadza poprzednie wyjścia z powrotem do obliczeń:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
Wyraz rekurencyjny −Σ a_k · y_{n−k} tworzy sprzężenie zwrotne. Pojedynczy impuls na wejściu będzie odbijać się wokół pętli sprzężenia zwrotnego w nieskończoność (zanikając geometrycznie, jeśli jest stabilny).
Dlaczego używać sprzężenia zwrotnego?
Filtr IIR może osiągnąć ostrą selektywność częstotliwościową z dużo mniejszą liczbą współczynników niż filtr FIR. IIR z 2 biegunami może przybliżyć to, co wymaga 50-współczynnikowego FIR dla tego samego tłumienia pasma zaporowego.
Cena: potencjalna niestabilność. Bieguny H(z) określają stabilność. Wszystkie bieguny muszą leżeć ściśle wewnątrz koła jednostkowego.
Historia prysznika Hamminga
Hamming użył żywej osobistej opowieści, aby zilustrować niestabilność sprzężenia zwrotnego.
Zatrzymywał się w tym samym pokoju hotelowym wielokrotnie, ponieważ znajomość pomagała mu się zorientować, gdy był zmęczony. Hydraulik zainstalował rury gorącej wody o dużej średnicy w prysznicu. Stworzyło to znaczące opóźnienie między regulowaniem zaworu a odczuciem zmiany temperatury wody.
Każdego ranka Hamming podążał tym samym wzorem: woda za zimna → podnieś ciepłą → wciąż zimna → podnieś bardziej → nagle wrzący → wyskoczyć → obniż → powtórz.
Opóźnienie w ścieżce sprzężenia zwrotnego oznaczało, że jego poprawki zawsze się przesadniały. Nie mógł przystosować się do opóźnienia, nawet po wielu powtórzeniach.
Lekcja inżynierii: niestabilność wynika albo z nadmiernego wzmocnienia w ścieżce sprzężenia zwrotnego ALBO z nadmiernego opóźnienia w ścieżce sprzężenia zwrotnego. Oba przejawiają się tym samym szukającym zachowaniem. W kategoriach filtrów: bieguny na lub poza kołem jednostkowym wytwarzają dokładnie tę oscylacyjną lub rozbiegającą się odpowiedź.
Charakteryzowanie niestabilności
Hamming zauważył, że tę samą niestabilność prysznica można analizować na dwa sposoby:
1. Jego odpowiedź była za silna (nadmierne wzmocnienie w działaniu korekcyjnym).
2. Jego wykrywanie było zbyt opóźnione (zbyt pośpieszne wejście do wanny, zanim system się ustalił).
Oba opisy dają ten sam rezultat matematyczny: biegun pętli sprzężenia zwrotnego przesunął się poza koło jednostkowe.
Cztery klasyczne rodziny
Teoria filtrów analogowych rozwinęła się wokół czterech klasycznych rodzin projektowych, każda reprezentująca inny kompromis. Te rodziny transformują się na czas dyskretny poprzez transformację biliniową lub niezmienność impulsową.
Butterworth (maksymalnie płaskie)
Odpowiedź pasma przenoszenia: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Monotonicznie malejące. Bez tętnień w paśmie przenoszenia lub paśmie zaporowym. Bieguny leżą na okręgu o promieniu ω_c na płaszczyźnie s (lub przekształconym okręgu na płaszczyźnie z). Najpłaściejsze możliwe pasmo przenoszenia dla danego rzędu N.
Chebyshev typ I
Równe tętnienia w paśmie przenoszenia, monotonne w paśmie zaporowym. Dla danego rzędu N i poziomu tętnienia, osiąga ostrzejsze odcięcie niż Butterworth. Bieguny leżą na elipsie (na płaszczyźnie s).
Chebyshev typ II
Równe tętnienia w paśmie zaporowym, monotonne w paśmie przenoszenia. Lustrzane odbicie typu I w dziedzinie częstotliwości.
Eliptyczne (Cauer)
Równe tętnienia w OBU paśmie przenoszenia i paśmie zaporowym. Dla danego rzędu N i poziomów tętnień, osiąga najostrzejsze możliwe przejście z pasma przenoszenia do pasma zaporowego. Używa funkcji eliptycznych do optymalnego umieszczenia biegunów i zer. Hamming: nazwa pochodzi z faktu, że funkcje eliptyczne są używane w wyprowadzeniu.
Fundamentalny kompromis
Wszystkie cztery rodziny osiągają ten sam podstawowy kompromis w różny sposób: wyższy rząd N daje ostrzejsze przejście. Pozwolenie na tętnienia (Chebyshev, eliptyczne) osiąga ostrzejsze przejście dla tego samego N. Eliptyczne osiąga absolutnie najostrzejsze przejście dla dowolnego danego N i specyfikacji tętnień.
Wybór między rodzinami filtrów
Wybór między rodzinami zależy od tego, co aplikacja toleruje.
Kwestionowanie roszczenia eksperta
Hamming przypomniał sobie, że pewni eksperci twierdzili, że wszystkie filtry IIR (rekurencyjne) posiadały określoną właściwość. Zapytał siebie, czy to naprawdę jest prawdą — i znalazł kontrprzykład.
Jego punkt: eksperci często noszą twierdzenia wchłonięte w szkole bez kiedykolwiek ich ponownego zbadania w obecnym kontekście. Jeśli zapytasz siebie, czy to, co się Ci mówi, naprawdę jest prawdą, zaskakuje, jak wiele okazuje się fałszywe, nawet w dobrze ugruntowanej dziedzinie.
Kontrprzykład nie był rodzajem filtra, który normalnie projektowałbyś, ale udowodnił, że roszczenie jest powierzchowne. Jeden kontrprzykład wystarczy do obalenia uniwersalnego roszczenia.
Projektowanie IIR w praktyce
Hamming zauważył, że niezależnie opracował wiele teorii filtrów IIR podczas rozwiązywania innego problemu: wyprowadzenia stabilnych formuł korektora dla numerycznych równań różniczkowych zwyczajnych.
Forma formuły korektora: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
Sprzężenie zwrotne pojawia się zarówno w warunkach y (sprzężenie zwrotne liniowe) jak i warunkach f(y) (sprzężenie zwrotne nieliniowe poprzez równanie różniczkowe). Stabilność dla filtrów IIR jest szczególnym przypadkiem bardziej ogólnego problemu stabilności dla numerycznych integratorów ODE.
Łączenie sprzężenia zwrotnego w różnych domenach
Ta sama struktura matematyczna — sprzężenie zwrotne, bieguny, granica stabilności — pojawia się w filtrach cyfrowych, numerycznych solverach ODE, systemach kontroli, biologicznych rytmach i modelach ekonomicznych.
W każdej domenie: pętla sprzężenia zwrotnego oblicza nowy stan na podstawie poprzednich stanów. Stabilność wymaga, aby sprzężenie zwrotne nie amplifikowało zaburzeń w nieskończoność.
Granica stabilności koła jednostkowego na płaszczyźnie Z odpowiada: osi urojonej na płaszczyźnie Laplace'a s (czas ciągły), warunkowi promienia spektralnego ρ(A) < 1 dla iteracji liniowych i warunkowi wykładnika Lyapunova λ < 0 dla systemów nieliniowych.