title: Lineaire fase
text: Lineaire fase is een eigenschap die bij FIR-filters kan worden ingebouwd. Deze eigenschap houdt in dat de faserespons van zo'n filter een lineaire functie van de frequentie is. Concreet wil dit zeggen dat een faseverschil tussen twee sinuscomponenten met verschillende frequentie hetzelfde blijft nadat het signaal door het filter gepasseerd is. Deze eigenschap, die kan worden vertaald in voorwaarden op de filtercoëfficiënten, is bij een aantal toepassingen gewenst. Bij IIR-filters is lineaire fase niet mogelijk.
  
  Voorwaarden op de filtercoëfficiënten
  
  Inleiding
  Een digitaal filter is een voorbeeld van een LTD-systeem. Bij een FIR-filter hoort een differentievergelijking die enkel rekening houdt met het invoersignaal (de excitatie) en niet met reeds vroeger berekende samples van het uitvoersignaal (de respons) zoals het geval is bij IIR-filters. De differentievergelijking is dus van de gedaante:
  
  waarbij  de filtercoëfficiënten zijn. Deze bepalen de wijze waarop het filter een invoersignaal omzet in een uitvoersignaal. Doorgaans is het filter bandselectief. Dit betekent dat bepaalde frequentiebanden worden doorgelaten en andere worden onderdrukt. Zo kan het laagdoorlaat, hoogdoorlaat, banddoorlaat of bandstop zijn. De bovengrens  van de sommatie noemt men de filterorde.
  
  Symmetrieën
  De filtercoëfficiënten kunnen eventueel aan bepaalde symmetrieën voldoen, namelijk:
   Even symmetrie:  voor 
   Oneven symmetrie:  voor 
  Daarnaast kan men opteren voor een waarde van  die zelf een even of een oneven getal is.
  Door combinatie van deze twee keuzes kan men dus vier gevallen bekomen. In elk van deze vier gevallen zal het filter een lineaire fase hebben.
  
  Geval 1:  is even,  is even
  Stel bijvoorbeeld dat , met als gevolg dat , ,  en  is willekeurig. De frequentierespons van een digitaal signaal is in het algemeen:
  
  waarbij  de genormaliseerde frequentie is, die gelijk is aan:
  
  met  de bemonsteringsperiode.
  In bovenstaand geval met even symmetrie en  oneven wordt dit:
  
  en dus:
  
  Door de laatste exponent voorop te zetten:
  
  Door de formules van Euler:
  
  De uitdrukking tussen haakjes is zuiver reëel en bevat dus enkel frequentie-intervallen waar de fase constant is, hetzij nul indien de uitdrukking positief is, hetzij  indien ze negatief is. Maar op de sporadische sprongen tussen deze twee waarden na is de fase constant. De complexe exponent vooraan het rechterlid bevat wel fase, maar dan een fase die duidelijk evenredig is met de genormaliseerde frequentie. Globaal is de fase dus, op sporadische sprongen na, een lineaire functie van de frequentie:
  
  In het algemeen geval wordt de uitdrukking voor de frequentierespons:
  
  waarbij:
  
  Andere gevallen
  In de andere gevallen kan een soortgelijke redenering gemaakt worden, die tot een soortgelijke uitdrukking leidt. Echter, in die gevallen treden beperkingen op waardoor die gevallen niet zo algemeen zijn als het eerste (bovenstaand uitgewerkte) geval.
  
   Type 2:  is even,  is oneven
  In dat geval vindt men:
  
  waarbij:
  
  Merk nu op dat, indien  deze uitdrukking steeds nul is, gezien het feit dat  enkel een halve waarden (0.5, 1.5, 2.5...) aanneemt. De cosinustermen worden dus, indien , steeds in  berekend zodat ze nul zijn. Dit type kan dus niet gebruikt worden voor hoogdoorlaat en bandstopdoorlaat filters. Deze filters hebben immers een niet-nulle frequentierespons aan de rand van het fundamenteel frequentie-interval.
  
   Type 3:  is oneven,  is oneven
  Dit type kent soortgelijke beperkingen als het vorige, en geeft steeds een nulle frequentierespons, zowel op  als op . Het kan dus enkel gebruikt worden voor banddoorlaat filters.
  
   Type 4:  is oneven,  is even
  Dit type kent eveneens soortgelijke beperkingen als de twee vorige en geeft steeds een nulle frequentierespons op . Het kan dus niet gebruikt worden voor laagdoorlaat filters en bandstopfilters.
  
  Samenvatting
  Van de vier gevallen biedt enkel het eerste een volledige vrijheid om het type van bandselectief filter te kiezen.
  
  Gevolgen voor de nullen
  De polen van een FIR-filter liggen steeds in . De nullen kunnen in principe overal liggen, maar bij een FIR-filter met een van bovenstaande symmetrieën voldoen de nullen aan een bijkomende eigenschap:
   indien  een nul is van de systeemfunctie, is  dit ook. 
  Deze eigenschap betreffende de posities van nullen kan als volgt aangetoond worden (voor de eenvoud voor een systeem van het eerste type, met  is even, maar de redenering kan zonder probleem worden veralgemeend).
  Stel bijvoorbeeld dat , zodat
  , 
  en stel dat Z een nul van de systeemfunctie is. Dit wil zeggen dat:
  
  Indien de systeemfunctie wordt geëvalueerd in  bekomt men:
  
  Dit toont aan dat, indien  een nul van het systeem is,  dit ook is.
  
  Deze redenering verklaart tevens waarom lineaire fase niet mogelijk is bij een IIR-systeem. Zo'n systeem heeft polen die niet in nul liggen, en opdat hun aanwezigheid de lineaire fase niet zou verstoren zouden ze ook aan een soortgelijke voorwaarde moeten voldoen: indien  een pool is moet  dit ook zijn. Echter, dit zou betekenen dat er polen buiten de eenheidscirkel liggen want indien een complex getal binnen deze cirkel ligt, ligt zijn invers complex getal er buiten. Polen buiten de eenheidscirkel zijn echter niet toegelaten omdat het systeem dat onstabiel is.
  
  Groepsvertraging
  De negatieve afgeleide van de frequentierespons noemt met de groepsvertraging of group delay:
  
  In het geval van een lineaire faserespons is de groepsvertraging dus constant.
  
  Voorbeeld
  
  Neem het filter van type 1, met , en met differentievergelijking:
  
  De nullen van het systeem zijn dan:
    en : dit zijn onderling gereflecteerde posities want 
   : deze nullen liggen op de eenheidscirkel
   : deze nullen liggen op de eenheidscirkel
   : dit is een koppel complex toegevoegde nullen
   : dit is een koppel complex toegevoegde nullen op de gereflecteerde posities van het vorig koppel; , z8=1/z10
  De amplituderespons en faserespons zijn te zien op bijgaande plots. Het systeem is een laagdoorlaat systeem. De twee nullen op de eenheidscirkel zijn terug te vinden als de nullen in de amplituderespons. De faserespons vertoont een lineair verloop, op enkele sprongen na. De sprongen van  in het faseverloop stemmen overeen met de nullen op de complexe eenheidscirkel. De sprongen van  zijn enkel een gevolg van het feit dat de fase tussen  en  geplot worden.
  
  Signaalanalyse