Zeitdiskrete Fourier-Transformation der Einheitsschrittfolge

10

Aus Lehrbüchern wissen wir, dass die DTFT von gegeben ist durchu[n]

(1)U(ω)=πδ(ω)+11ejω,πω<π

Ich habe jedoch kein DSP-Lehrbuch gesehen, das zumindest vorgibt, eine mehr oder weniger solide Ableitung von .(1)

Proakis [1] leitet die rechte Hälfte der rechten Seite von indem er in der -Transformation von setzt und sagt, dass es gültig ist außer (was natürlich richtig ist). Er gibt dann an, dass wir am Pol der -Transformation einen Delta-Impuls mit einer Fläche von hinzufügen müssen , aber das erscheint mir eher wie ein Rezept als alles andere.(1)z=ejωZu[n]ω=2πkZπ

Oppenheim und Schafer [2] erwähnen in diesem Zusammenhang

Obwohl es nicht ganz einfach zu zeigen ist, kann diese Sequenz durch die folgende Fourier-Transformation dargestellt werden:

Darauf folgt eine Formel, die . Leider haben sie sich nicht die Mühe gemacht, uns diesen "nicht ganz einfachen" Beweis zu zeigen.(1)

Ein Buch, das ich eigentlich nicht kannte, das ich aber bei der Suche nach einem Beweis für ist die Einführung in die digitale Signalverarbeitung und das Filterdesign von BA Shenoi. Auf Seite 138 gibt es eine "Ableitung" von , aber leider ist es falsch. Ich habe eine "DSP-Puzzle" -Frage gestellt , damit die Leute zeigen, was mit diesem Beweis nicht stimmt.](1)(1)

Also meine Frage ist:

Kann jemand einen Beweis / eine Ableitung von liefern?, der solide oder sogar streng ist, während er für mathematisch veranlagte Ingenieure zugänglich ist? Es spielt keine Rolle, ob es nur aus einem Buch kopiert wurde. Ich denke, es wäre trotzdem gut, es auf dieser Seite zu haben.(1)

Beachten Sie, dass selbst in math.SE fast nichts Relevantes zu finden ist: Diese Frage hat keine Antworten, und eine hat zwei Antworten, von denen eine falsch ist (identisch mit Shenois Argument), und die andere verwendet die "Akkumulationseigenschaft". , mit dem ich mich freuen würde, aber dann muss man diese Eigenschaft beweisen, was Sie wieder an den Anfang bringt (weil beide Beweise im Grunde dasselbe beweisen).

Als letzte Anmerkung habe ich mir so etwas wie einen Beweis ausgedacht (nun, ich bin Ingenieur), und ich werde ihn in einigen Tagen auch als Antwort veröffentlichen, aber ich würde gerne andere veröffentlichte oder unveröffentlichte Beweise sammeln Das ist einfach und elegant und vor allem für DSP-Ingenieure zugänglich.

PS: Ich bezweifle nicht die Gültigkeit von (1) , ich möchte nur einen oder mehrere relativ einfache Beweise sehen.


[1] Proakis, JG und DG Manolakis, Digitale Signalverarbeitung: Prinzipien, Algorithmen und Anwendungen , 3. Auflage, Abschnitt 4.2.8

[2] Oppenheim, AV und RW Schafer, Zeitdiskrete Signalverarbeitung , 2. Auflage, S. 54.



Inspiriert von einem Kommentar von Marcus Müller möchte ich zeigen, dass U(ω) nach Gl. (1) erfüllt die Anforderung

u[n]=u2[n]U(ω)=12π(UU)(ω)

Wenn U(ω) die DTFT von u[n] , dann

V(ω)=11ejω

muss die DTFT von sein

v[n]=12sign[n]

(wo wir ), weilsign[0]=1

V(ω)=U(ω)πδ(ω)u[n]12=12sign[n]

Also haben wir

12π(VV)(ω)(12sign[n])2=14

woraus folgt das

12π(VV)(ω)=DTFT{14}=π2δ(ω)

Damit bekommen wir

12π(UU)(ω)=12π[(πδ(ω)+V(ω))(πδ(ω)+V(ω))]=12π[π2δ(ω)+2πV(ω)+(VV)(ω)]=π2δ(ω)+V(ω)+π2δ(ω)=U(ω)q.e.d.
Matt L.
quelle
waaah. Mach meine Welt nicht kaputt. Der Zweifel an dieser Formel führt zu einem Bereich des Chaos. Zum Beispiel ist und daher (mit einem Cont. FT-Definitionsvorfaktor in Abhängigkeit von der Konstante c ) DTFT ( u 2 ) ( ω )u2(t)=u(t)c
DTFT(u2)(ω)=cU(ω)U(ω)=cπU(ω)+c1+ejωU(ω)=cπ(πδ(ω)+11ejω)+cπ1+ejω+c11+ejω11+ejω=cπ2δ(ω)+2cπ1+ejω+c11+ejω11+ejω=magic?U
Marcus Müller
@ MarcusMüller: Kein Zweifel an dieser Formel, sie ist richtig. Die Frage ist nur, wie man es auf eine Weise zeigt, die ein einfacher Ingenieur verstehen kann. Und funktioniert für die gegebene DTFT, kein Problem. u2[n]=u[n]
Matt L.
Ich halte mich für sehr einfältig, und das bedeutet, dass ich mir Sorgen mache, wenn sich die Dinge nicht "sicher" anfühlen, wenn ich nicht sehen kann, wie sie abgeleitet werden.
Marcus Müller
1
Ich sehe, dass es nicht darum geht, zu beweisen, ob die Gleichung korrekt ist oder nicht, sondern U ( w ) rigoros und direkt aus den ersten Prinzipien und der Definition von DTFT abzuleiten . Wenn man dann einen strengen Beweis mit Impulsen erbringen will, sollte man sich besser auf die zitierten Bücher aus der verallgemeinerten Funktionstheorie beziehen: Lighthill-1958 wird in Opp & Schafer für eine Diskussion der Impulsfunktion und ihrer Verwendung in Fourier-Transformationen zitiert. Alle anderen Beweise stützen sich unweigerlich auf die Beweise, die auf diesen Referenzen gemacht wurden, und reichen nicht aus, um einen strengen Beweis zu ersetzen. U(w)
Fat32
1
@ Fat32: Das ist ein gültiger Standpunkt. Ich denke jedoch, dass eine einigermaßen solide Ableitung möglich ist, wenn wir grundlegende Transformationen wie akzeptieren und wenn wir uns damit zufrieden geben, Integrale durch ihren Cauchy-Hauptwert zu definieren. DTFT{1}=2πδ(ω)
Matt L.

Antworten:

3

Cedron Dawg hat in dieser Antwort einen interessanten Ausgangspunkt gepostet . Es beginnt mit diesen Schritten:

U(ω)=n=0+ejωn=limNn=0N1ejωn=limN[1ejωN1ejω]=11ejωlimN[ejωN1ejω]

Es stellt sich heraus, dass der Begriff innerhalb des Grenzwerts wie folgt erweitert werden kann :

ejωN1ejω=1sin2(ω)+(1cos(ω))2[cos(ω)cos(Nω)+cos(Nω)sin(ω)sin(Nω)+j(sin(ω)cos(Nω)+cos(ω)sin(ω)sin(Nω))]

Der gemeinsame Faktor außerhalb der Klammern kann ausgedrückt werden als :

1sin2(ω)+(1cos(ω))2=14sin2(ω/2)

Der Realteil in den Klammern ist auch gleich :

cos(ω)cos(Nω)+cos(Nω)sin(ω)sin(Nω)=2sin(ω/2)sin[ω(N+1/2)]

Andererseits kann der Imaginärteil wie folgt umgeschrieben werden :

sin(ω)cos(Nω)+cos(ω)sin(ω)sin(Nω)=2sin(ω/2)cos[ω(N+1/2)]

Wenn wir den ursprünglichen Begriff umschreiben, erhalten wir Folgendes:

ejωN1ejω=2sin(ω2)4sin2(ω2)(sin[ω(N+1/2)]jcos[ω(N+1/2)])=sin[ω(M+1/2)]2sin(ω2)jcos[ω(M+1/2)]2sin(ω2)

M=N1M

Gemäß der 7. Definition auf dieser Site :

limM12sin(ω/2)sin[ω(M+1/2)]=πδ(ω)

Bisher haben wir das:

limMejω(M+1)1ejω=πδ(ω)jlimMcos[ω(M+1/2)]2sin(ω/2)

0

U(ω)=11ejωlimN[ejωN1ejω]=11ejω+πδ(ω)+jlimMcos[ω(M+1/2)]2sin(ω/2)=11ejω+πδ(ω)
Tendero
quelle
Das ist sehr nett! Ich habe es überprüft und alles scheint korrekt zu sein, so dass der Imaginärteil in gewissem Sinne gegen Null tendieren muss. Ich werde ein bisschen darüber nachdenken.
Matt L.
@ MattL. Lassen Sie mich wissen, ob Sie Fortschritte machen können!
Tendero
@ MattL. Der Beweis ist endlich vollständig!
Tendero
ω=0ω01/(1ejω)ω=0
2

Ich werde zwei relativ einfache Beweise liefern, die keine Kenntnisse der Verteilungstheorie erfordern. Einen Beweis, der die DTFT durch einen Grenzwertprozess unter Verwendung von Ergebnissen aus der Verteilungstheorie berechnet, finden Sie in dieser Antwort von Tendero .

Ich werde hier nur den ersten Beweis erwähnen (und nicht näher darauf eingehen), da ich ihn als Antwort auf diese Frage veröffentlicht habe , um zu zeigen, dass ein bestimmter veröffentlichter Beweis fehlerhaft ist.

u[n]

(1)ue[n]=12(u[n]+u[n])=12+12δ[n]

(1)

(2)DTFT{ue[n]}=πδ(ω)+12

u[n]

(3)UR(ω)=Re{U(ω)}=πδ(ω)+12

u[n]U(ω)UR(ω)U(ω)h[n]=h[n]u[n]h[n]H(ω)=12π(HU)(ω)u[n]U(ω)

UI(ω)=Im{U(ω)}u[n]

(4)uo[n]=12(u[n]u[n])=u[n1]12+12δ[n]

(4)

jUI(ω)=ejωU(ω)πδ(ω)+12=ejω(UR(ω)+jUI(ω))πδ(ω)+12=ejω(πδ(ω)+12)+ejωjUI(ω)πδ(ω)+12(5)=12(1+ejω)+ejωjUI(ω)

(3)(5)

(6)jUI(ω)(1ejω)=12(1+ejω)

(6)

(7)jUI(ω)=121+ejω1ejω+cδ(ω)

UI(ω)ωu[n]c=0(3)(7)

U(ω)=UR(ω)+jUI(ω)=πδ(ω)+12+121+ejω1ejω=πδ(ω)+12(1+1+ejω1ejω)(8)=πδ(ω)+11ejω
Matt L.
quelle