Gibt es eine Möglichkeit, die Impulsantwort eines diskreten Systems zu erhalten, indem man nur weiß, wie es auf die diskrete Einheitsschrittfunktion reagiert?


10

In ununterbrochener Zeit war es möglich;

u(t)systemy(t)δ(t)=du(t)dtsystemdy(t)dt=h(t)

Gilt das auch für das diskrete Zeitsystem, dh

δ[t]=du[t]dtwhere:{δ[t]is the discrete time deltau[t]is the discrete time unit step function

Gibt es eine Möglichkeit, die Impulsantwort eines diskreten Systems zu erhalten, indem nur die Antwort des diskreten Einheitsschritts bekannt ist?


1
Tolle Frage! Willkommen bei DSP.SE. Bleib dabei und trage dazu bei!
Phonon

Antworten:


7

Eine einfachere Version von Phonons Antwort lautet wie folgt.

Angenommen, bezeichnet die Reaktion des Systems auf die Einheitsschrittfunktion. Dann ist , wie in dieser Antwort diskutiert , im Allgemeinen die Summe von skalierten und zeitverzögerten Kopien der Impulsantwort, und in diesem speziellen Fall ist keine Skalierung erforderlich; nur Zeitverzögerungen. Somit wo jeweils Die rechte Spalte enthält eine (nicht skalierte und) zeitverzögerte Impulsantwort. Somit erhalten wir leicht, dass y y [ 0 ]y yh [ 0 ]

y[0]=h[0]y[1]=h[1]+h[0]y[2]=h[2]+h[1]+h[0]y[3]=h[3]+h[2]+h[1]+h[0] = 
h[0]=y[0]h[1]=y[1]y[0]h[2]=y[2]y[1] = h[n] =y[n]y[n1] = 
Ohne Erwähnung von Filtern, Inversen, Faltungen, Integration, Operatoren und dergleichen, nur einfache Konsequenzen der Definition eines linearen zeitinvarianten Systems.

Sie haben dies eindeutig länger getan als ich =)
Phonon

6

Ja, dies gilt auch für diskrete Systeme. Die Differenzierungsoperation wird in diesem Fall durch eine Differenz erster Ordnung ersetzt. Es glaubt nicht, dass es ein universelles Symbol hat, aber nennen wir es . Diese Operation entspricht dem Filtern Ihres Signals mit . Nennen wir diesen Filter . Ich werde Faltung als das Symbol bezeichnen.D()y[n]=x[n]x[n1]d[n]

Wenden wir nun das, was wir über Faltung wissen, auf diesen Operator an. Wir wissen, dass wir mit einer laufenden Summe (diskreter Integrator) auf . Tatsächlich erweist sich das durch selbst dargestellte System als dieser diskrete Integrator. Beachten Sie auch, dass diese beiden Operatoren invers zueinander sind und dass speziell .u[n]δ[n]u[n]u[n]d[n]=δ[n]

Jetzt wissen wir, dass die Faltung kommutativ ist, dh

a[n]b[n]=b[n]a[n]

und assoziativ, das heißt

(a[n]b[n])c[n]=a[n](b[n]c[n])

Also ist

x[n]=δ[n]x[n]=u[n]d[n]x[n]=d[n]u[n]x[n]=d[n](u[n]x[n])

Sie können also sehen, dass Sie von wiederherstellen können, indem Sie die Differenz erster Ordnung anwenden, genau wie im kontinuierlichen Fall.x[n](u[n]x[n])


2

Annahmen:

  • Kontinuierlicher Zeitbereich : Sei - Impulsantwort und - Schrittantworth(t)s(t)
  • Diskreter Zeitbereich : Es sei - Einheitsimpulsantwort und - Einheitsschrittantworth[n]s[n]

Intuitiv entspricht die Integration in den kontinuierlichen Zeitbereich der Summierung im diskreten Zeitbereich. In ähnlicher Weise entspricht die Ableitung im kontinuierlichen Zeitbereich der endlichen Differenz im diskreten Bereich.

Betrachten Sie mit diesem intuitiven Verständnis die Beziehung zwischen und (linke Seite der zweiten Gleichung in Ihrem Beitrag):uδ

  • Kontinuierlicher Zeitbereich: u(t)=δ(t)
  • Diskreter Zeitbereich: u[n]=k=0δ[nk]

Betrachten Sie in ähnlicher Weise die Beziehung zwischen und (rechte Seite der zweiten Gleichung in Ihrem Beitrag):sh

  • Kontinuierlicher Zeitbereich: s(t)=h(t)
  • Diskreter Zeitbereich: s[n]=k=0h[nk]

Wenn Sie sich nun die letzte Gleichung genau ansehen:

s[n]=k=0h[nk]

Nun kann aus dieser Gleichung unter Verwendung der endlichen Differenz von mit einer verzögerten Version von sich selbst, dh mit :h[n]s[n]s[n1]

h[n]=s[n]s[n1]
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.