LTD-systeem

Uit Wikipedia, de vrije encyclopedie
Ga naar: navigatie, zoeken

Een Lineair Tijdinvariant Discreet Systeem, kortweg LTD-systeem is de digitale of discrete tegenhanger van het analoge of continue LTC-systeem. Een LTD-systeem verwerkt één of meerdere digitale ingangssignalen (excitaties genoemd) tot één of meerdere uitgangssignalen (responsen of responsies genoemd), en dit op een lineaire manier. Dit wil zeggen dat een lineaire combinatie van excitaties wordt omgezet in dezelfde lineaire combinatie van de afzonderlijke responsen. Tijdinvariantie betekent dat, indien de excitaties in de tijd worden verschoven, de responsen ongewijzigd blijven, behalve dat ze over een gelijk tijdsinterval worden verschoven als de excitaties.Een digitaal signaal bestaat uit een reeks opeenvolgende getallen, die men sampels noemt. LTD-systemen worden toegepast in digitale signaal verwerking, onder andere bij digitale filters en digitale regelaars.

Differentievergelijking[bewerken]

Een digitaal signaal met één excitatie x[n] en één respons y[n] wordt beschreven door een lineaire differentievergelijking van de vorm:

y[n] \, = \, \sum_{k=0}^{N} a_k \, x[n-k] \, + \, \sum_{i=1}^{M} b_i \, y[n-i]

In het rechterlid staat:

  • steeds een lineaire combinatie van een aantal sampels van het excitatiesignaal. Deze sampels kunnen de actuele sampel x[n] zijn, en/of een aantal vorige sampels.
  • eventueel ook nog een lineaire combinatie van een aantal reeds berekende sampels van het responssignaal, dus sampels y[n-1] en/of vorige.

Het feit dat dit in beide gevallen lineaire combinaties zijn garandeert de lineariteit van het systeem. Het feit dat de coëfficiënten van de lineaire combinaties niet van de parameter n afhangen, en dus niet van de tijd, garandeert de tijdinvariantie.

Soorten systemen[bewerken]

Bovenstaand onderscheid suggereert een opdeling in:

  • Niet-recursief LTD-systeem : hierbij is de uitgangsampel y[n] enkel functie van de excitatie. De tweede sommatie in de bovenstaande differentievergelijking valt dus weg.
  • Recursief LTD-systeem : hierbij bevat de differentievergelijking wel beide sommaties

In de praktijk wordt een andere opdeling gebruikt die bijna perfect op hetzelfde neerkomt. Deze opdeling is gebaseerd op de aard van de respons op een digitale Diracpuls:

 \delta[n] = \left\{ \begin{matrix} 1 & als & n = 0 \\ 0 & als & n \neq 0 \end{matrix}\right.

De respons op dit excitatiesignaal, de impulsrespons kan al dan niet eindig zijn in lengte of oneindig. Men onderscheidt dus:

  • FIR-systeem : Finite Impulse Response-systeem. In dit geval is de impulsrespons eindig
  • IIR-systeem : Infinite Impulse Response-systeem. In dit geval is de impulsrespons oneindig

Doorgaans worden de twee opdelingen als perfect met elkaar overeenstemmend beschouwd waarbij FIR-systemen als niet-recursief worden beschouwd, en IIR-systemen als recursief. In dit artikel zal deze gewoonte worden aangehouden, hoewel het wiskundig mogelijk is IIR-systemen te maken die wel een eindige impulsrespons hebben.

De impulsrespons en de systeemfunctie[bewerken]

Indien de differentievergelijking van een LTD-systeem wordt geëxciteerd met een digitale Diracpuls wordt de bijhorende respons per definitie de impulsrespons h[n]. De systeemfunctie is dan per definitie de z-getransformeerde van de impulsrespons. Men kan deze twee bewerkingen ook omwisselen: Men z-transformeert eerst de differentievergelijking en vervangt vervolgens de z-getransformeerde X(z) van de excitatie vervangt door 1, de z-getransformeerde van een digitale Diracpuls). Hierdoor wordt de z-getransformeerde Y(z) van de excitatie automatisch de gevraagde systeemfunctie H(z).

Voor een FIR-systeem wordt dit:

H_{FIR}(z) \, = \sum_{k=0}^{N} \, a_k \, z^{-k}

of nog:

H_{FIR}(z) \, = \frac{\sum_{k=0}^{N} \, a_k \, z^{N-k}}{z^N}

Dit toont aan dat het FIR-systeem N nullen heeft, de oplossingen van de teller. Verder zijn er ook N polen, de oplossingen van de noemer, die allemaal in z = 0 liggen. Om stabiliteitsredenen moeten de polen binnen de eenheidscirkel in het complexe vlak liggen. Bij een FIR-systeem is dit steeds voldaan zodat een FIR-systeem steeds stabiel is.

Voor een IIR-systeem wordt de systeemfunctie:

H_{IIR}(z) \, = \frac{\sum_{k=0}^{N} \, a_k \, z^{-k}}{1-\sum_{i=1}^{M} \, b_i \, z^{-i}}

of nog:

H_{IIR}(z) \, = z^{M-N} \frac{\sum_{k=0}^{N} \, a_k \, z^{N-k}}{z^M-\sum_{i=1}^{M} \, b_i \, z^{M-i}}

In dit geval heeft het systeem opnieuw evenveel polen als nullen, maar liggen de polen niet noodzakelijk binnen de eenheidscirkel. De recursieve coëfficiënten b_k moeten dus zodanig bepaald worden dat dit wel het geval is.

De frequentierespons[bewerken]

De frequentierespons van een digitaal systeem ontstaat door de complexe variabele z van de systeemfunctie te beperken tot de eenheidscirkel. Dit is het digitale equivalent van wat men bij analoge systemen doet: daar wordt de complexe variabele s van de Laplacetransformatie beperkt tot de imaginaire as, om toegang te krijgen tot de reële frequenties. De frequentierespons van een FIR-systeem wordt dus:

H_{FIR}(\theta) \, = \sum_{k=0}^{N} \, a_k \, e^{-jk\theta}

waarbij de frequentievariabele \theta is. Dit is een genormaliseerde frequentievariabele die het interval [-\pi,\pi] doorloopt. Fysisch stemt dit interval overeen met een reëel frequentiebereik in het interval [-fs/2,fs/2], waar fs de bemonsteringsfrequentie van de signalen is. Nullen op de eenheidscirkel hebben dus tot gevolg dat een reële frequentie volledig wordt onderdrukt. Een pool op de eenheidscirkel zou betekenen dat een reële frequentie oneindig wordt versterkt. De absolute waarde van de frequentierespons heet de amplituderespons, de complexe hoek of argument de faserespons.

Lineaire fase[bewerken]

Bij een FIR-systeem kan een lineaire faserespons bekomen worden, indien de rij filtercoëffficënten [ a_0 ... a_N ] ongewijzigd blijft indien men ze achterstevoren leest. Dus moet:

a_j \, = \, a_{N-j} \quad \forall j = 0..N \!

Een ander gevolg en tevens voorwaarde voor lineaire fase is de onderlinge ligging van de nullen. Voor elke (complexe) nul op afstand R van oorsprong moet er een andere nul op de gereflecteerde positie liggen, dus op een afstand 1/R van de oorsprong, onder dezelfde complexe hoek. Een nul op de eenheidscirkel voldoet hier automatisch aan, en is zijn eigen gereflecteerd punt. Bij een IIR-systeem is een perfect lineaire fase niet mogelijk.

Voorbeelden[bewerken]

Een FIR systeem[bewerken]

Amplitude en faserespons van het FIR-voorbeeld. De frequentieas op beide figuren is de gereduceerde frequente theta. Onderaan staat de reële as in Hertz.

Een eenvoudig FIR-systeem is bijvoorbeeld:

y[n] \, = \, \frac{1}{6} (x[n] \, + \, 2x[n-1] \, + \, 2x[n-2] \, + \, x[n-3])

Dit is een voorbeeld van een zogenaamde averager, een systeem dat een gemiddelde berekent. In dit geval is de n-de sampel van de respons een gewogen gemiddelde van de n-de sampel van de input en zijn drie voorgangers. De impulsrespons is:

 h[n=0,1,...] \, = \, \frac{1}{6} \quad \frac{1}{3} \quad \frac{1}{3} \quad \frac{1}{6} \quad 0 \quad 0 \quad 0 \quad ... \!

Zoals verwacht is heeft de impulsrespons een eindige lengte, en zijn de sampels van de impulsrespons gewoon de filtercoëfficiënten. De systeemfunctie is:

H(z) \, = \, \frac{1+2z+2z^2+z^3}{6z^3}

De drie polen liggen in de oorsprong z = 0, en de drie nullen zijn -1, e^{j2\Pi/3} en e^{-j2\Pi/3}. Deze liggen alle drie op de eenheidscirkel, en stemmen overeen met de helft van de bemonsteringsfrequentie, en met plus en min één derde van de bemonsteringsfrequentie. De nevenstaande figuur toont de amplituderespons, die voor hogere frequentie klein wordt. Dit is steeds het geval bij een averager want die middelt snelle veranderingen uit. De faserespons is lineair, gezien de waarden van de filtercoëfficiënten [1/6,1/3,1/3,1/6] voldoen aan de eis voor lineaire fase.

Een IIR-systeem[bewerken]

De differentievergelijking:

y[n] \, = \, \frac{1}{2}x[n] \, + \, \frac{1}{2}y[n-1] \!

is een eerste orde IIR-systeem. De impulsrespons is:

 h[n=0,1,...] \, = \, \frac{1}{2} \quad \frac{1}{4} \quad \frac{1}{8} \quad ... \quad \frac{1}{2^{n+1}} \quad ... \!

De systeemfunctie is:

H(z) \, = \, \frac{1}{2-z^{-1}} \, = \, \frac{z}{2z-1}

Dit systeem heeft een nul in de oorsprong, en een pool in z = 1/2. Het is ook een milde averager. Een constant signaal zal met gelijke waarde worden doorgelaten, want de frequentierespons op frequentie gelijk aan nul is gelijk aan 1. Hoogfrequentie schommelingen worden gedeeltelijk onderdrukt. Ter hoogte van de helft van de bemonsteringsfrequentie is de doorlating van de amplituderespons nog 1/3.

Bronnen[bewerken]

Een inleiding tot digitale systemen is te vinden in vrijwel alle boeken en bronnen over digitale signaalverwerking (DSP, digital signal processing).

  • E.C. Ifeachor, B.W. Jervis "Ditigal Signal Processing, a practical approach" Addison-Wesley Publ.Company ISBN 0-201-54431-8
  • S.K. Mitra "Digital Signal Processing, a computer based approach" McGraw-Hill International Edition ISBN 007-124467-0
  • A.W.M. Van den Enden, N.A.M. Verhoeckx "Discrete Time Signal Processing" Prentice Hall, ISBN 0-13-216755-7

Zie ook[bewerken]