Où est la faille dans cette dérivation de la DTFT de la séquence de pas unitaire


11

Cette question est liée à cette autre question à moi où je demande des dérivations de la transformée de Fourier à temps discret (DTFT) de la séquence de pas unitaire . Au cours de ma recherche de dérivations, j'en ai trouvé une qui est incroyablement simple. Je l'ai vu pour la première fois à la page 138 de ce livre de BA Shenoi. Je l'ai également rencontré sur les mathématiques.SE dans cette réponse .u[n]

Puisque l'argument est court et simple, je le répéterai ici pour plus de commodité.

La séquence de pas d'unité peut être écrite comme avec Évidemment, Application de la DTFT des deux côtés de donne où est la DTFT de . De nous obtenons De et nous obtenons pour la DTFT de

(1)u[n]=f[n]+12
f[n]-f[n-1]=δ[n](3)F(ω)(1-e-jω)=1F(ω)f[n](4)F(ω)=1
(2)f[n]={12,n012,n<0
(3)f[n]f[n1]=δ[n]
(3)
(4)F(ω)(1ejω)=1
F(ω)f[n](4) (5)(1)u[n]U(ω)=F(ω)+πδ(ω)=1
(5)F(ω)=11ejω
(5)(1)u[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
où j'ai utilisé , .- π ω < πDTFT{1}=2πδ(ω)πω<π

Eq. pour le DTFT de est sans aucun doute correct. Cependant, la dérivation est défectueuse.u [ n ](6)u[n]

La question est: trouver et expliquer la faille dans la dérivation ci-dessus.

Veuillez ajouter votre réponse avec la balise spoiler >!.


1
ce qui me dérange, c'est que est un signal de puissance finie , pas un signal d'énergie finie , ce que nous obtenons lorsque nous additionnons ces deux signaux d'énergie infinis. f[n]
robert bristow-johnson

également, ?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
robert bristow-johnson

Merci les gars pour vos réponses! Je les ai tous votés positivement, et chacun donne lieu à une discussion intéressante sur les aspects peu connus de la DTFT des signaux étranges (c'est-à-dire ceux qui ne sont pas dans ou ). Je ne peux en accepter qu'une et j'attendrai un peu plus longtemps pour de nouvelles réponses ou des changements dans les réponses existantes. J'ajouterai également ma propre réponse plus tard. 212
Matt L.

1
Matt, n'est décidément pas une énergie finie. un nombre infini d'échantillons qui carré pour être ne s'ajoute pas pour être un nombre fini. 1f[n]14
robert bristow-johnson

1
@ robertbristow-johnson: Qu'est-ce qui vous dérange à ce sujet? Si les signaux s'annulent partout sauf pour un nombre fini de points, c'est ce que nous obtenons.
Matt L.

Réponses:


7

Il existe une infinité de signaux qui maintiennent l'égalité suivante: La seule chose qui compte est que , puis le reste des coefficients de peut être déterminé sous la restriction que Eq. indique (c'est-à-dire que la soustraction d'échantillons consécutifs doit être pour ). En d'autres termes, l'Eq. sera atteint par tout signal tel que Une autre façon de voir c'est que toute fonction qui est fondamentalement avec un décalage (une valeur ajoutée constante) satisfera

y[n]y[n1]=δ[n](1)
y[0]y[1]=1y(1)0n0(1)y[n]
y[0]=y[1]+1y[n]=y[n1] n0
u[n](1). Cela explique la déclaration faite par Robert Bristow-Johnson dans sa réponse : les différenciateurs détruisent ces informations (comme prendre un dérivé en temps continu détruit la preuve de toute valeur constante dans la fonction d'origine).

Pour résumer, je pense que la preuve est erronée car la procédure suivie pourrait utiliser n'importe quelle fonction de la forme avec , et cela conduirait à de nombreuses fonctions ayant la même transformée de Fourier , ce qui est effectivement faux car la transformée de Fourier est une bijection. L'auteur a peut-être délibérément décidé d'ignorer tout ce qui était lié aux valeurs DC, conscient que pour montrer que est le DTFT de il aurait besoin de la propriété d'accumulation (dont la preuve la plus populaire est dérivée du DTFT de la marche unitaire - ergo, une jolie épreuve circulaire). La preuve n'est pas strictement fausse , car tout ce qu'elle énonce (les formules pour etu[n]+CCRF(ω)f[n]F ( ω ) U ( ω ) F ( ω )F(ω)U(ω), la décomposition du pas unitaire, l'équation de différence) est vraie, mais il faudrait la propriété d'accumulation pour montrer pourquoi n'a pas de deltas Dirac.F(ω)


Vous êtes totalement sur la bonne voie! Avez-vous une idée de la façon dont cette faille pourrait être résolue, c'est-à-dire, comment le faire correctement?
Matt L.

@MattL. Définir une condition initiale pour ferait l'affaire et déterminerait le signal de manière univoque. Cette condition initiale déterminerait la valeur DC du signal y [ n ] , qui apparaît dans le DTFT comme une constante multipliant une impulsion de Dirac (selon la propriété d'accumulation). Je pense que dans la preuve donnée, cela fonctionne parce que le signal f [ n ] n'a pas de valeur DC car il est symétrique autour de 0 , et donc le DTFT est correct juste dans ce cas. Mais le fait que le signal n'a pas de DC doit être déclaré, car c'est fondamental, je crois. y[n]y[n]f[n]0
Tendero

Il existe de nombreuses bonnes réponses et il est difficile de choisir laquelle accepter. Mais celui-ci a été le plus apprécié par la communauté, et je pense aussi qu'il souligne le plus clairement l'erreur dans la dérivation. Merci à tous!
Matt L.

4

J'ai été submergé par le nombre de réponses que j'ai reçues (10 réponses jusqu'à présent!). Bien sûr, tous ont obtenu mon vote positif. C'était amusant, merci les gars pour vos réflexions, commentaires, etc. Je sais que la plupart d'entre vous savent maintenant quel est le défaut, du moins celui que je voulais dire. Les gens expriment les choses différemment, et il y a toujours place à des malentendus, alors j'essaierai de formuler clairement ce que je pense être le défaut le plus important dans cette dérivation. Je sais que tout le monde ne sera pas d'accord et c'est très bien. Je suis heureux de pouvoir discuter de ce genre de sujets DSP ésotériques avec des esprits aussi vifs que vous tous! Et c'est parti.

Ma première affirmation est que chaque équation de ma question est correcte. Cependant, la dérivation et la motivation de certains d'entre eux sont totalement fausses et trompeuses, et cette "dérivation" ne peut exister que parce que l'auteur savait à quoi devait ressembler le résultat.

Eq. (3) dans la question ( ) est correct pour la séquence f [ n ] donnée (Eq. ( 2 ) dans la question), mais il est clairement aussi correct pour toutes les séquences de la forme f [ n ] = u [ n ] + c avec une constante arbitraire c . Ainsi, selon la dérivation, le DTFT F résultant (f[n]f[n1]=δ[n]f[n](2)

(1)f[n]=u[n]+c
c doit être la DTFT de toutes les séquences de la forme ( 1 ) , quelle que soit la valeur de la constante c . C'est évidemment insensé car le DTFT est unique. Plus précisément, en utilisant cette très "preuve", je pourrais "montrer" que F ( ω ) comme indiqué dans l'équation. ( 5 ) de ma question (ou l'équation ( 3 ) ci-dessous) est en fait la DTFT de u [ n ] que nous recherchons. Alors pourquoi s'embêter à diviser u [ n ] comme dans l'équation. ( 1 ) de la question?F(ω)(1)cF(ω)(5)(3)u[n]u[n](1)

Cependant, il est vrai que les DTFT de toutes les séquences satisfont à l'Eq. ( 4 ) dans la question (répétée ici pour plus de commodité): F ( ω ) ( 1 - e - j ω ) = 1 Mais vient maintenant le défaut mathématique réel: De l'équation. ( 2 ) il est incorrect de conclure F ( ω ) = 1(1)(4)

(2)F(ω)(1ejω)=1
(2) Éq. (3)n'est qu'une des infiniment nombreuses solutions possibles de(2), et il se trouve que c'est celle dont l'auteur a besoin pour arriver au résultat final correct. Eq. (3)est la DTFT def[n]dans(1)avecc=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1) , mais à partir de la dérivation donnée, il n'y a aucun moyen de le savoir.c=12

Alors, comment pouvons-nous éviter cette erreur mathématique et utiliser pour dériver les DTFT de a l l séquences ( 1 ) , avec une constante c ? La conclusion correcte de ( 2 ) est F ( ω ) = 1(2)all(1)c(2)avec une constanteαencore indéterminée. Le branchement(4)sur le côté gauche de(2)donne1+α(1-e-jω)δ(ω)=1+α(1-e-jω)| ω=0

(4)F(ω)=11ejω+αδ(ω)
α(4)(2) Donc toutes les fonctions F ( ω ) données par ( 4 ) satisfont ( 2 ) , comme requis.
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2)

La constante dans ( 4 ) peut être déterminée à partir de la valeur de f [ n ] à n = 0 : f [ 0 ] = 1 + c = 1α(4)f[n]n=0 On peut montrer, etWolframAlpha convientégalement, que la valeur principale de Cauchy de l'intégrale dans(6)estPV π - π dω

(6)f[0]=1+c=12πππF(ω)dω=12πππdω1ejω+α2π
(6)De(6)et(7)on obtientα=π(1+2c)Donc pourc=-1
(7)PVππdω1ejω=π
(6)(7)
(8)α=π(1+2c)
on obtientα=0(ce qui correspond à la séquence d'originef[n]telle qu'utilisée par l'auteur de la preuve), et pourc=0(c'est-à-dire pourf[n]=u[n]) on aα=π, ce qui nous donne finalement la DTFT souhaitée deu[n]: U ( ω ) = 1c=12α=0f[n]c=0f[n]=u[n]α=πu[n]
(9)U(ω)=11ejω+πδ(ω)


"toutes les fonctions données par (4) satisfont (2)", mais faut-il prouver que "toutes les fonctions F ( ω ) satisfont (2) ont la forme (4)"? F(ω)F(ω)
AlexTP

@AlexTP: Vous voulez donc dire que les fonctions de la forme ne peuvent être qu'un sous-ensemble de fonctions satisfaisant ( 2 ) ? Voilà un point valable. Mais je pense qu'il est assez clair qu'il ne peut pas y avoir d'autres fonctions, car ce n'est qu'à ω = 0 que nous obtenons le problème, nous avons donc besoin de fonctions qui ont une contribution supplémentaire à ω = 0 qui disparaît lorsqu'elle est multipliée par ( 1 - e - j ω ) . Ces fonctions (en fait les distributions) sont l'impulsion delta de Dirac et ses dérivés. Cependant, les dérivés ne disparaissent pas lorsqu'ils sont multipliés par ((4)(2)ω=0ω=0(1ejω) , il ne reste donc que l'impulsion du delta de Dirac(1ejω)
Matt L.

Je ne suis pas vraiment sûr qu'il ne puisse y avoir d'autre fonction que l'impulsion delta de Dirac (et ses dérivés) qui ait cette propriété. Mais ça va, votre réponse est bien écrite. Je vote positivement. Merci.
AlexTP

2

La faille suit le mot "Evidemment", si cela est censé être la fonction delta de Dirac.

Voici le brouillon d'une réponse à votre autre question que je n'ai jamais postée:

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

Je ne pense pas qu'une preuve soit possible. Cela peut être le cas d'une "définition fonctionnelle" ayant les propriétés souhaitées.

U = + n = 0 e - j ω n U = lim N N - 1 n = 0 e - j ω n U = lim N [

X2π(ω)=n=+x[n]ejωn
U=n=0+ejωn
U=limNn=0N1ejωn
U=limN[1ejωN1ejω]
Examen de la dernière valeur limite. Pourω=0,il est clair qu'il agit comme un delta de Dirac. Pourquoi le coefficient devrait êtreπ, je ne sais pas. Cela peut avoir à voir avec l'aire du cercle unitaire. Lorsqueω0, le dénominateur peut être factorisé hors de la limite et le numérateur saute simplement le long du cercle unitaire et n'atteint jamais une limite. Le mettre à zéro est un acte de définition.
U=11ejωlimN[ejωN1ejω]
ω=0πω0

Prouver la définition fonctionne d'une manière souhaitable est une autre affaire.

La preuve de la page 138 est fausse (au moins) parce que:

δ(t)=lima012a[u(t+a)u(ta)]=dudt
δ(n)=u2(n)u2(n1)

Situation intéressante, j'espère que cela vous aidera. J'attends avec impatience ce que vous avez à dire.

Ced


δ[n]n=01

U=11ejωlimN[ejωN1ejω]

2

1=2

F(ω)(1ejω)=1
ω=2kπ for kZ

3
u[n]u[n]12

w2πk

δ[n] ωω=2kπ

δ[n]1ejwF(w)

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

2

limNn=Nn=Nf[n]ejωn(1ejωn)+(ejωNf[N]+ejωNf[N])ejω=1
f[n]=u[n]
limNn=Nn=Nf[n]ejωn(1ejωn)+ejωNejω=1

F(ω)ω2kπk

F(ω)ω=2kπn=1sin(ωn)ωω=2kπU(ω)

F(ω)

u[n]

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

2

Je pense avoir trouvé la meilleure façon d'exprimer la faille dans cette preuve. Je vais donc lui donner un autre coup de couteau.

12x

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

x12

De plus, si vous prenez l'étape que j'ai faite dans ma dernière réponse et que la conclusion (4) est exprimée comme suit:

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

Suivi en l'incluant dans (5) et (6) vous obtenez:

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

Ce qui, comme je l'ai souligné plus tôt, n'est pas conforme à la définition pour y arriver.

x=12πδ(ω)

x=12

Ced


1

C'est en réponse aux commentaires de ma première réponse. En raison du camouflage des spoilers, je le poste comme une réponse distincte.

J'allais poster mon autre réponse à l'autre question, mais je ne l'ai pas fait en raison de mon manque d'expérience dans ce domaine. Je l'ai posté hier, l'ai supprimé, puis supprimé, puis compris comment utiliser les balises de spoiler.

δδδp

δp[n]=F[n]-F[n-1]=u[n]-u[n-1]

Prendre la DTFT des parties gauche et droite. Je ne suis pas sûr d'avoir les bonnes notations, mais les calculs devraient être clairs. Utilisation de la définition en cours de preuve.

Fp(ω)=Fu(ω)-Fu(ω)e-jω

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

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

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

ω=2kπω=2kπ

Ced

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

Suivre:

δp


Fp(ω)=1(1ejω)δ(ω)f(ω)ω=0f(ω)δ(ω)=f(0)δ(ω)f(0)=0

1

(a[n]b[n])cω[n]=a[n]cω[n]b[n]cω[n]
12f[n]f[n1]

1

alors Matt,

f[n]

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

α>0

nous avons maintenant des signaux d'énergie finie et les DTFT devraient tous être comparables.

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.
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.