Wo liegt der Fehler bei dieser Ableitung der DTFT der Einheitsschrittfolge


11

Diese Frage bezieht sich auf meine andere Frage, bei der ich nach Ableitungen der zeitdiskreten Fourier-Transformation (DTFT) der Einheitsschrittfolge frage u[n]. Bei meiner Suche nach Ableitungen habe ich eine gefunden, die erstaunlich einfach ist. Ich habe es zum ersten Mal auf Seite 138 dieses Buches von BA Shenoi gesehen. Ich bin in dieser Antwort auch auf Mathematik gestoßen .

Da das Argument kurz und einfach ist, werde ich es hier der Einfachheit halber wiederholen.

Die Einheitsschrittfolge kann geschrieben werden als

(1)u[n]=f[n]+12
mit
(2)f[n]={12,n012,n<0
Offensichtlich ergibt
(3)f[n]f[n1]=δ[n]
Anwenden der DTFT auf beiden Seiten von(3)ergibt
(4)F(ω)(1ejω)=1
wobeiF(ω)ist die DTFT vonf[n]. Von(4)bekommen wir
(5)F(ω)=11ejω
Aus(5)und(1)wir für die DTFT vonu[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
wo ichDTFT{1}=2πδ(ω),πω<π.

Gl. (6) für die DTFT von u[n] ist zweifellos richtig. Die Ableitung ist jedoch fehlerhaft.

Die Frage ist: Finden und erklären Sie den Fehler in der obigen Ableitung.

Bitte stellen Sie Ihrer Antwort das Spoiler-Tag voran >!.


1
Was mich stört, ist, dass ein endliches Leistungssignal ist , kein endliches Energiesignal , was wir erhalten, wenn wir diese beiden unendlichen Energiesignale addieren. f[n]
Robert Bristow-Johnson

ist nicht auch ?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
Robert Bristow-Johnson

Danke Jungs für eure Antworten! Ich habe sie alle positiv bewertet, und jeder führt zu einer netten Diskussion über nicht so bekannte Aspekte der DTFT von seltsamen Signalen (dh diejenigen, die nicht in oder 2 enthalten sind ). Ich kann nur eine akzeptieren und warte etwas länger auf neue Antworten oder Änderungen an vorhandenen Antworten. Ich werde später auch meine eigene Antwort hinzufügen. 12
Matt L.

1
Matt, ist definitiv keine endliche Energie. eine unendliche Anzahl von Proben, deren Quadrat 1 istf[n] Addiere nicht, um eine endliche Zahl zu sein. 14
Robert Bristow-Johnson

1
@ robertbristow-johnson: Was stört dich daran? Wenn sich die Signale bis auf eine endliche Anzahl von Punkten überall gegenseitig aufheben, dann bekommen wir das.
Matt L.

Antworten:


6

Es gibt unendlich viele Signale, die die folgende Gleichheit gelten lassen:

y[n]y[n1]=δ[n](1)
Das einzige, was zählt, ist, dassy[0]y[1]=1 , und dann kann der Rest der Koeffizienten vony unter der Einschränkung bestimmt werden, dass Gl. (1) Zustände (dh die Substraktion von aufeinanderfolgenden Proben müssen0 fürn0 ). Mit anderen Worten, Gl. (1) wird durch jedes Signaly[n] so dass
y[0]=y[1]+1y[n]=y[n1] n0
Eine andere Möglichkeit, dies zu sehen, besteht darin, dass jede Funktion, die im Grunde genommenu[n] mit einem Versatz (einer konstanten Wertschöpfung) ist, erfüllt(1) . Dies erklärt die Aussage vonRobert Bristow-Johnson in seiner Antwort: Unterscheidungsmerkmale zerstören diese Informationen (z. B. wenn eine Ableitung in kontinuierlicher Zeit genommen wird, werden Hinweise auf einen konstanten Wert in der ursprünglichen Funktion zerstört).

Zusammenfassend glaube ich, dass der Beweis fehlerhaft ist, weil das angewandte Verfahren jede Funktion der Form u[n]+C mit C R verwenden könnte , und dies würde dazu führen, dass viele Funktionen dieselbe Fourier-Transformation haben, was in der Tat falsch ist da die Fourier-Transformation eine Bijektion ist. Vielleicht hat der Autor absichtlich beschlossen, alles zu ignorieren, was mit DC-Werten zusammenhängt, um zu zeigen, dass F ( ω ) die DTFT von f [ n ] ist. er würde die Akkumulationseigenschaft benötigen (deren populärster Beweis aus der DTFT des Einheitsschrittes abgeleitet wird - ergo, ein ziemlich kreisförmiger Beweis). Der Beweis ist nicht streng falsch , wie alles, was er sagt (die Formeln für F ( ω)CRF(ω)f[n]F(ω) undU(ω) , die Zerlegung des Einheitsschritts, die Differenzgleichung), wahr ist, aber es würde die Akkumulationseigenschaft erfordern, um zu zeigen, warumF(ω) hat keine Dirac-Deltas.


Sie sind auf dem richtigen Weg! Haben Sie eine Idee, wie dieser Fehler behoben werden könnte, dh wie Sie es richtig machen?
Matt L.

@ MattL. Das Einstellen einer Anfangsbedingung für würde den Trick machen und das Signal eindeutig bestimmen. Diese Anfangsbedingung würde den Gleichstromwert des Signals y [ n ] bestimmen , der in der DTFT als Konstante erscheint, die einen Dirac-Impuls multipliziert (gemäß der Akkumulationseigenschaft). Ich denke, dass dies im gegebenen Beweis funktioniert, weil das Signal f [ n ] keinen Gleichstromwert hat, da es um 0 symmetrisch ist , und daher ist die DTFT nur in diesem Fall korrekt. Aber die Tatsache, dass das Signal keinen Gleichstrom hat, sollte angegeben werden, da es meiner Meinung nach grundlegend ist.y[n]y[n]f[n]0
Tendero

Es gibt viele gute Antworten und es ist schwer zu entscheiden, welche man akzeptiert. Aber dieser wurde von der Community am meisten geschätzt, und ich denke auch, dass er am deutlichsten auf den Fehler bei der Ableitung hinweist. Danke euch allen!
Matt L.

4

Ich war überwältigt von der Anzahl der Antworten (10 Antworten bisher!). Natürlich haben alle meine Gegenstimme bekommen. Das hat Spaß gemacht, danke Jungs für eure Gedanken, Kommentare usw. Ich weiß, dass die meisten von euch inzwischen wissen, was der Fehler ist, zumindest der, den ich gemeint habe. Die Leute drücken die Dinge anders aus und es gibt immer Raum für Missverständnisse. Deshalb werde ich versuchen, klar zu formulieren, was meiner Meinung nach der wichtigste Fehler in dieser Ableitung ist. Mir ist bewusst, dass nicht jeder zustimmen wird und das ist in Ordnung. Ich bin froh, diese Art von esoterischen DSP-Themen mit so scharfem Verstand wie Sie diskutieren zu können! Auf geht's.

Meine erste Behauptung ist, dass jede Gleichung in meiner Frage richtig ist. Die Ableitung und Motivation einiger von ihnen ist jedoch völlig falsch und irreführend, und diese "Ableitung" kann nur existieren, weil der Autor wusste, wie das Ergebnis aussehen sollte.

Gl. (3) in der Frage ( ) ist für die gegebene Folge f [ n ] korrekt (Gleichung ( 2 ) in der Frage), aber es ist eindeutig auch korrekt für alle Folgen der Form f [ n ] = u [ n ] + c mit einer beliebigen Konstante c . Entsprechend der Ableitung ergibt sich also die resultierende DTFT F (f[n]f[n1]=δ[n]f[n](2)

(1)f[n]=u[n]+c
c sollte die DTFT aller Sequenzen der Form ( 1 ) sein n ] wie in Gl. ( 1 ) der Frage?F(ω)(1)unabhängig vom Wert der Konstanten . Das ist natürlich unsinnig, weil die DTFT einzigartig ist. Insbesondere könnte ich unter Verwendung dieses "Beweises" "zeigen", dass F ( ω ) wie in Gl. ( 5 ) meiner Frage (oder Gleichung ( 3 ) unten) ist tatsächlich die DTFT von u [ n ] , nach der wir suchen. Warum also die Mühe machen, dich aufzuteilen ?cF(ω)(5)(3)u[n]u[n](1)

Es ist jedoch wahr, dass die DTFTs aller Sequenzen Gl. ( 4 ) in der Frage (hier der Einfachheit halber wiederholt): F ( ω ) ( 1 - e - j ω ) = 1 Nun kommt aber der eigentliche mathematische Fehler: Aus Gl. ( 2 ) Es ist falsch zu schließen(1)(4)

(2)F(ω)(1ejω)=1
(2) Gl. (3)ist nur eine von unendlich vielen möglichen Lösungen von(2), und es ist zweckmäßigerweise diejenige, die der Autor benötigt, um zum richtigen Endergebnis zu gelangen. Gl. (3)ist die DTFT vonf[n]in(1)mitc=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1) , aber aus der gegebenen Ableitung gibt es keine Möglichkeit, das zu wissen.c=12

Wie können wir also diesen mathematischen Fehler vermeiden und , um die DTFTs von a l l Sequenzen ( 1 ) mit einer beliebigen Konstante c abzuleiten ? Die korrekte Schlussfolgerung aus ( 2 ) ist F ( ω ) = 1(2)all(1)c(2)mit einer noch unbestimmten Konstanteα. Das Einstecken von(4)in die linke Seite von(2)ergibt1+α(1-e-jω)δ(ω)=1+α(1-e-jω)| ω=0

(4)F(ω)=11ejω+αδ(ω)
α(4)(2) also alle Funktionen F ( ω ) gegebendurch ( 4 ) erfüllen ( 2 ) ,nach Bedarf.
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2)

Die Konstante in ( 4 ) kann aus dem Wert von f [ n ] bei n = 0 bestimmt werden : f [ 0 ] = 1 + c = 1α(4)f[n]n=0 Es kann gezeigt werdenund auchWolfram stimmt, dass der Cauchy Hauptwert des Integrals in(6)istPV π - π dω

(6)f[0]=1+c=12πππF(ω)dω=12πππdω1ejω+α2π
(6)Aus(6)und(7) erhaltenwirα=π(1+2c)Also fürc=-1
(7)PVππdω1ejω=π
(6)(7)
(8)α=π(1+2c)
wirα=0(was der ursprünglichen Sequenzf[n] entspricht,wie sie vom Autor des Beweises verwendet wird), und fürc=0(dh fürf[n]=u[n]) haben wirα=π, was uns schließlich die gewünschte DTFT vonu[n] gibt: U ( ω ) = 1c=12α=0f[n]c=0f[n]=u[n]α=πu[n]
(9)U(ω)=11ejω+πδ(ω)


F(ω)F(ω)

(4) could only be a subset of functions satisfying (2)? That's a valid point. But I think it's pretty clear that there can't be any other functions, because it's only at ω=0 where we get the problem, so we need functions that have an additional contribution at ω=0 which disappears when multiplied by (1ejω). Such functions (actually distributions) are the Dirac delta impulse and its derivatives. However, the derivatives do not disappear when multiplied by (1ejω), so it's only the Dirac delta impulse that's left
Matt L.

I am not really sure that there can't be any function other than Dirac delta impulse (and its derivatives) that has this property. But it is ok, your answer is well-written. I upvote. Thanks.
AlexTP

2

The flaw follows the word "Obviously", if that is supposed to be the Dirac Delta Function.

Here is the draft of an answer for your other question that I never posted:

---------------------------------------------------------------

I don't think a proof is possible. This may be a case of a "functional definition" having desired properties.

X2π(ω)=n=+x[n]ejωn
U=n=0+ejωn
U=limNn=0N1ejωn
U=limN[1ejωN1ejω]
U=11ejωlimN[ejωN1ejω]
Looking at the last limit value. For ω=0 it is clear that it acts like a Dirac Delta. Why the coefficient should be π, I don't know. It may have to do with the area of the unit circle. When ω0, the denominator can be factored out of the limit and the numerator just jumps along the unit circle and never reaches a limit. Setting it to zero is a definitional act.

Proving the definition works in a desirable way is a different matter.

The page 138 proof is wrong (at least) because:

δ(t)=lima012a[u(t+a)u(ta)]=dudt
Which is not similar in any way to δ(n)=u2(n)u2(n1) as they define it.

Interesting situation, I hope this helps. I am looking forward to what you have to say.

Ced


Thanks for your answer! But note that we're talking about discrete-time sequences here, not about the Dirac delta impulse. In discrete-time, δ[n] equals zero everywhere, except for n=0 where it is 1. So no ugly stuff as in continuous-time. So Eq. (3) in my question is actually correct (in discrete time!).
Matt L.

Cedron, I think this equation:
U=11ejωlimN[ejωN1ejω]
may be of interest in this other question of @MattL.. Maybe you should consider giving it some more thought and posting it there if you want to.
Tendero

2

if you allow me to divide by zero, i can prove to you that 1=2. when you say

F(ω)(1ejω)=1
there is a problem about multiplying something by zero (when ω=2kπ for kZ) and expecting the product to equal one.

3
moral of the story: differentiators destroy information. a differentiator doesn't know the difference between u[n] and u[n]12. don't try dividing by zero to get that lost information.
robert bristow-johnson

so should the author state that w2πk to correct the flaw?
Fat32

well, the DTFT of δ[n] is equal to 1 for any ω, even when ω=2kπ.
robert bristow-johnson

why δ[n] ?... it's the domain of 1ejw which restricts that of F(w) ?
Fat32

because when Matt says
f[n]f[n1]=δ[n]
and says to integrate both sides to get u[n] (which is what Matt says essentially in Eq.5), then we have a sorta "multiply and divide by the same quantity trick to get 1". but sometimes he is multiplying and dividing by zero. multiplying by zero destroys information. dividing that by zero will not get the information back.
robert bristow-johnson

2

The Equation (4) should be written as

limNn=Nn=Nf[n]ejωn(1ejωn)+(ejωNf[N]+ejωNf[N])ejω=1
For f[n]=u[n],
limNn=Nn=Nf[n]ejωn(1ejωn)+ejωNejω=1
that is not (5). I don't know how to fix the proof without avoiding (3).

someone might argue that F(ω) exists when ω2kπ for integer k.
robert bristow-johnson

@robertbristow-johnson Thanks for editting. I mean that F(ω) does not exist not because ω=2kπ but because n=1sin(ωn) does not converge for every ω. For ω=2kπ, it must be interpreted as a Cauchy principal value even for the correct formula of U(ω)
AlexTP

@AlexTP Why is divergence of that summation reason to say that F(ω) doesn't exist? Many signals are such that their "summation DTFTs" diverge, but they have a DTFT defined nevertheless (in terms of Dirac deltas). I'm not saying the reasoning its wrong, I'm just trying to understand (I'm not too familiar with distribution theory and that).
Tendero

You're right that the sums do not converge in the conventional sense, but nevertheless, the DTFTs exist in a distributional sense. After all, the same argument would hold for the DTFT of u[n], which also exists in that sense. So there's another flaw which is less mathematically complex than the divergence problem, which we've kind of gotten used to.
Matt L.

1
@MattL. I dont know if it is obvious but I dont see the difference between u[n] and f[n] because we also have u[n]u[n1]=δ[n] !
AlexTP

2

I think I have figured out the best way to express the flaw in this proof. So I am going to give it another stab.

The choice of 12 in (1) is arbitrary. Let's replace it with x. Follow the proof through, and end up with:

U(ω)=11ejω+2πxδ(ω)

There is nothing in the proof which constrains x to be 12, it can take any finite value and the proof works the same.

Furthermore, if you take the step I did in my last answer and find (4) is expressed as

F(ω)(1ejω)=1+2πx(1ejω)δ(ω)

Followed by including it in (5) and (6) you get:

U(ω)=11ejω+4πxδ(ω)

Which, as I pointed out earlier, is inconsistent with the definition to get there.

This proof fails to show that x=12, and it seems to indicate that for any defined x, inconsistent results will follow. Therefore I go back to the statement in my first answer that the value of π for the coefficient of δ(ω) is a definitional act, not an mathematical truism.

Perhaps there is some other situation which makes x=12 the correct value, but this proof doesn't provide it.

Ced


1

This is in response to the comments in my first answer. Because of the spoiler cloaking I am posting it as a separate answer.

I was going to post my other answer to the other question, but I didn't due to my lack of experience in this area. I posted it yesterday, deleted it, then undeleted, then figured out how to employ spoiler tags.

Clearly the δ function defined in the problem is not the Dirac Delta function. I looked up DTFT in Wikipedia and the DTFT for the Dirac Delta function is one. I will call the δ of the problem δp.

δp[n]=f[n]f[n1]=u[n]u[n1]

Taking the DTFT of the left and right parts. I'm not sure I have the notations right, but the math should be clear. Using the definition that is being proved.

Fp(ω)=Fu(ω)Fu(ω)ejω

Fp(ω)=[11ejω+πδ(ω)][ejω1ejω+(πejω)δ(ω)]

Fp(ω)=1ejω1ejω+π(1ejω)δ(ω)

Fp(ω)=1+π(1ejω)δ(ω)1

Thus the RHS of (4) is incorrect except when ω=2kπ. [Edit: Doh, it's the Dirac Delta, so this statement is wrong. I guess it should be correct except "undefined" at ω=2kπ. Real Analysis was my least favorite math. I am leaving this alone now.]

Ced

==============================

Followup:

It is clear that the DTFT of δp should be 1 when it is plugged into the definition of a DTFT. Therefore, since I got a different answer when using the definition to be proved means that the definition to be proved is not correct (in a mathematical sense). Furthermore, if you carry the correction through to the end of the proof you arrive at a different definition. Assuming the assertion to be true is used to prove that it is false.


You've actually shown that Fp(ω)=1 because the term (1ejω)δ(ω) equals zero. That's the case because for any function f(ω) that is continuous at ω=0 you have f(ω)δ(ω)=f(0)δ(ω), and if f(0)=0 (which is the case here), the whole term disappears.
Matt L.

1

To me, a first flaw appears between (3) and (4): this is an instance of the classical careless integral/infinite sum splitting. Conditions are required to allow the equation:

(a[n]b[n])cω[n]=a[n]cω[n]b[n]cω[n]
Standard 1 or 2 conditions might not be sharp enough. This could be related here, due to the form f[n]f[n1], to the Fubini derivation theorem, or: When can we exchange infinite sum and discrete derivative? Ways around could revolve around monotony or Cesaro-like sums, but I shall think about this longer.

1

so Matt,

i dunno why you don't think it's problematic comparing power signals to energy signals, but suppose we modify the definition of f[n] slightly:

f[n]{ 12eαnn012eαnn<0

for some α>0.

now we have finite energy signals and the DTFTs should all be comparable.

f[n]f[n1]={12(eαneα(n1))n>012(1+eα)n=012(eαneα(n1))n<0={12(1eα)eαnn>012(1+eα)n=012(eα1)eαnn<0

i wonder what the DTFTs are? and then what happens when we let α0 ? i think there is still the problem of differentiators destroying information (and the corresponding destruction of information by multiplying by 0 in the frequency domain) that is a problem. but maybe we can lose the problem of comparing signal classes that don't share the same Hilbert space.

but, alas, it's nearly 2 a.m. and i ain't gonna deal with it now.


That α thing is good, and it's one option to compute the DTFT of those non-decaying signals, by taking the limit α0. Try it, and I'm sure you'll succeed, but it's painful. There are easier ways to get the same result. The given proof can actually be modified such that it works IMHO (see my answer).
Matt L.
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.