Где недостаток в этом выводе DTFT последовательности единичных шагов ?


11

Этот вопрос связан с другим моим вопросом, где я спрашиваю о выводах дискретного преобразования Фурье (DTFT) последовательности единичных шагов . Во время поиска производных я нашел удивительно простую. Впервые я увидел это на странице 138 этой книги Б. А. Шеноя. Я тоже наткнулся на это на mathematics.SE в этом ответе .u[n]

Поскольку аргумент короткий и простой, я повторю его здесь для удобства.

Последовательность единичных шагов может быть записана в виде с помощью Очевидно, Применение DTFT с обеих сторон от дает где - DTFT для . Из мы получаем Из и получаем DTFT для

(1)u[n]=f[n]+12
(2)f[n]={12,n012,n<0
(3)f[n]f[n1]=δ[n]
(3)
(4)F(ω)(1ejω)=1
F(ω)f[n](4)
(5)F(ω)=11ejω
(5)(1)u[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
где я использовал DTFT{1}=2πδ(ω) , πω<π .

Eq. (6) для DTFT из u[n] , без сомнения, правильно. Тем не менее, вывод является ошибочным.

Вопрос в том, чтобы найти и объяснить недостатки в приведенном выше выводе.

Пожалуйста, добавьте ваш ответ с тегом спойлера >!.


1
меня беспокоит то, что - это сигнал конечной мощности , а не сигнал конечной энергии , что мы и получаем, когда складываем эти два сигнала бесконечной энергии вместе. f[n]
Роберт Бристоу-Джонсон

кроме того, не ?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
Роберт Бристоу-Джонсон

Спасибо, ребята, за ваши ответы! Я проголосовал за всех из них, и каждый из них приводит к хорошему обсуждению не очень известных аспектов DTFT странных сигналов (то есть тех, которые не в или ). Я могу принять только один, и я буду ждать немного дольше новых ответов или изменений в существующих ответах. Я также добавлю свой собственный ответ позже. 212
Мэтт Л.

1
Мэтт, определенно не конечная энергия. бесконечное число выборок, которые должны быть квадратными, чтобы быть , не добавляют к конечному числу. 1f[n]14
Роберт Бристоу-Джонсон

1
@ robertbristow-johnson: Что вас беспокоит? Если сигналы отменяют друг друга везде, кроме конечного числа точек, то это то, что мы получаем.
Мэтт Л.

Ответы:


7

Существует бесконечно много сигналов, обеспечивающих выполнение следующего равенства: Единственное, что имеет значение, это то, что , и тогда остальные коэффициенты могут быть определены при условии, что уравнение состояния (т. Е. Вычитание последовательных отсчетов должно быть для ). Другими словами, уравнение будет достигнуто любым сигналом таким, что Еще один способ увидеть это то, что любая функция, которая в основном является со смещением (добавленное постоянное значение), будет удовлетворять

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), Это объясняет утверждение, сделанное Робертом Бристоу-Джонсоном в его ответе : дифференциаторы уничтожают эту информацию (например, получение производной в непрерывном времени уничтожает свидетельство любой постоянной величины в исходной функции).

Подводя итог, я полагаю, что доказательство является ошибочным, потому что использованная процедура могла использовать любую функцию вида с , и это привело бы ко многим функциям, имеющим одно и то же преобразование Фурье , что действительно неправильно, так как преобразование Фурье является биекцией. Возможно, автор намеренно решил игнорировать все, что связано со значениями DC, сознавая, что для того, чтобы показать, что является DTFT для ему потребуется свойство накопления (наиболее популярное доказательство которого получено из DTFT шаг единицы - эрго, довольно круговое доказательство). Доказательство не является строго неверным , так как все, что оно утверждает (формулы для иu[n]+CCRF(ω)f[n]F(ω)U(ω)Разложение единичного шага (разностное уравнение) верно, но для накопления свойства нет дельта Дирака, поэтому требуется свойство накопления .F(ω)


Вы на правильном пути! У вас есть идея, как можно устранить этот недостаток, то есть как сделать это правильно?
Мэтт Л.

@MattL. Установка начального условия для сделает свое дело и определит сигнал однозначно. Это начальное условие будет определять значение постоянного тока сигнала , который появляется в DTFT как постоянная, умножающая импульс Дирака (согласно свойству накопления). Я думаю, что в данном доказательстве это работает, потому что сигнал не имеет значения постоянного тока, так как он симметричен относительно , и поэтому DTFT является правильным только в этом случае. Но тот факт, что сигнал не имеет постоянного тока, должен быть заявлен, так как я считаю, что это фундаментально. y [ n ] f [ n ] 0y[n]y[n]f[n]0
Тендеро

Есть много хороших ответов, и трудно выбрать, какой из них принять. Но этот вопрос был оценен сообществом больше всего, и я думаю, что он наиболее четко указывает на ошибку в выводе. Спасибо вам всем!
Мэтт Л.

4

Я был ошеломлен количеством полученных ответов (пока 10 ответов!). Конечно, все они получили мое одобрение. Это было весело, спасибо, ребята, за ваши мысли, комментарии и т. Д. Я знаю, что к настоящему времени большинство из вас знает, что это за недостаток, по крайней мере тот, который я имел в виду. Люди выражают вещи по-разному, и всегда есть место для недопонимания, поэтому я постараюсь четко сформулировать то, что я считаю наиболее важным недостатком в этом выводе. Я осознаю тот факт, что не все согласятся, и это нормально. Я рад возможности обсуждать такие эзотерические темы DSP с такими острыми умами, как вы! Вот так.

Мое первое утверждение состоит в том, что каждое уравнение в моем вопросе является правильным. Однако вывод и мотивация некоторых из них совершенно неверны и вводят в заблуждение, и этот «вывод» может существовать только потому, что автор знал, как должен выглядеть результат.

Eq. (3) в вопросе ( ) верно для данной последовательности f [ n ] (уравнение ( 2 ) в вопросе), но это также очевидно корректно для всех последовательностей вида f [ n ] = u [ n ] + c с некоторой произвольной константой c . Итак, согласно выводу, результирующий DTFT F (f[n]f[n1]=δ[n]f[n](2)

(1)f[n]=u[n]+c
c должно быть DTFT всех последовательностей вида ( 1 ) независимо от значения константы c . Это, конечно, бессмысленно, потому что DTFT уникален. В частности, используя это самое «доказательство», я мог бы «показать», что F ( ω ), как указано в формуле. ( 5 ) моего вопроса (или уравнения ( 3 ) ниже) на самом деле является DTFT от u [ n ], который мы ищем. Так зачем ломать голову u [ n ], как в формуле. ( 1 ) вопроса?F(ω)(1)cF(ω)(5)(3)u[n]u[n](1)

Однако верно, что DTFT всех последовательностей удовлетворяют уравнению. ( 4 ) в вопросе (повторяется здесь для удобства): F ( ω ) ( 1 - e - j ω ) = 1 Но теперь возникает фактический математический недостаток: Из уравнения. ( 2 ) неверно заключить, что F ( ω ) = 1(1)(4)

(2)F(ω)(1ejω)=1
(2) уравнение (3)является лишь одним из бесконечно многих возможных решений(2), и это удобно для того, чтобы автор пришел к правильному конечному результату. Eq. (3)является DTFT дляf[n]в(1)сc=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1) , но из данного вывода нет способа узнать это.c=12

Так как мы можем избежать , что математическая ошибка и использование вывести DTFTs из через л л последовательностей ( 1 ) с любой постоянной с ? Правильный вывод из ( 2 ) : F ( ω ) = 1(2)all(1)c(2)с некоторой пока неопределенной константойα. Подстановка(4)в левую часть(2)дает1+α(1-e-jω)δ(ω)=1+α(1-e-jω)| ω=0

(4)F(ω)=11ejω+αδ(ω)
α(4)(2) Таким образом, все функции F ( ω ), заданные ( 4 ), удовлетворяют ( 2 ) , как требуется.
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2)

Константа в ( 4 ) может быть определена из значения f [ n ] при n = 0 : f [ 0 ] = 1 + c = 1α(4)f[n]n=0 Можно показать, иВольфрам Альфа согласен с тем, что главным значением интеграла Коши в(6)являетсяPV π - π dω

(6)f[0]=1+c=12πππF(ω)dω=12πππdω1ejω+α2π
(6)Из(6)и(7)получаемα=π(1+2c).Поэтому дляc=-1
(7)PVππdω1ejω=π
(6)(7)
(8)α=π(1+2c)
мы получаемα=0(что соответствует исходной последовательностиf[n],использованной автором доказательства), а дляc=0(т. Е. Дляf[n]=u[n]) имеемα=π, что в итоге дает нам желаемый DTFT дляu[n]: U ( ω ) = 1c=12α=0f[n]c=0f[n]=u[n]α=πu[n]
(9)U(ω)=11ejω+πδ(ω)


«все функции заданные (4), удовлетворяют (2)», но должны ли мы доказать, что «все функции F ( ω ) удовлетворяют (2), имеют вид (4)»? F(ω)F(ω)
AlexTP

@AlexTP: То есть вы имеете в виду, что функции вида могут быть только подмножеством функций, удовлетворяющих ( 2 ) ? Это верный момент. Но я думаю, что довольно ясно, что не может быть никаких других функций, потому что только при ω = 0 мы получаем проблему, поэтому нам нужны функции, которые имеют дополнительный вклад при ω = 0, который исчезает при умножении на ( 1 - e - j ω ) . Такими функциями (собственно распределениями) являются дельта-импульс Дирака и его производные. Однако производные не исчезают при умножении на ((4)(2)ω=0ω=0(1ejω) , так что остался только дельта-импульс Дирака(1ejω)
Мэтт Л.

Я не совсем уверен, что не может быть никакой функции, кроме дельта-импульса Дирака (и его производных), обладающего этим свойством. Но это нормально, ваш ответ хорошо написан. Я приветствую. Спасибо.
AlexTP

2

Недостаток следует за словом «Очевидно», если предполагается, что это дельта-функция Дирака.

Вот черновик ответа на ваш другой вопрос, который я никогда не публиковал:

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

Я не думаю, что доказательства возможны. Это может быть случай «функционального определения», имеющего желаемые свойства.

U = + Е п = 0 е - J ω п U = Нт N N - 1 Е п = 0 е - J ω п U = Нт N [

X2π(ω)=n=+x[n]ejωn
U=n=0+ejωn
U=limNn=0N1ejωn
U=1
U=limN[1ejωN1ejω]
Глядя на последнее предельное значение. Дляω=0ясно, что он действует как дельта Дирака. Почему коэффициент должен бытьπ, я не знаю. Это может быть связано с площадью круга юнитов. Когдаω0, знаменатель может быть разложен вне предела, и числитель просто перепрыгивает вдоль единичного круга и никогда не достигает предела. Установка его в ноль является определяющим актом.
U=11ejωlimN[ejωN1ejω]
ω=0πω0

Доказательство определения работает желательным образом, это другой вопрос.

Страница 138 доказательство неверно (по крайней мере), потому что:

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

Интересная ситуация, надеюсь, это поможет. Я с нетерпением жду того, что вы скажете.

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

Я думаю, что я нашел лучший способ выразить недостаток в этом доказательстве. Так что я собираюсь дать ему еще один удар.

12x

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

x12

Кроме того, если вы сделаете шаг, который я сделал в моем последнем ответе, и обнаружите, что (4) выражается как

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

Затем, включив его в (5) и (6), вы получите:

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

Что, как я указывал ранее, несовместимо с определением, как туда попасть.

x=12πδ(ω)

x=12

CED


1

Это в ответ на комментарии в моем первом ответе. Из-за маскировки спойлера я выкладываю это как отдельный ответ.

Я собирался опубликовать другой ответ на другой вопрос, но не сделал этого из-за недостатка опыта в этой области. Я отправил это вчера, удалил это, затем удалил, затем выяснил, как использовать признаки спойлера.

δδδp

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

Взяв DTFT из левой и правой частей. Я не уверен, что у меня есть правильные обозначения, но математика должна быть ясной. Используя определение, которое доказывается.

Fп(ω)знак равноFU(ω)-FU(ω)е-Jω

Fп(ω)знак равно[11-е-Jω+πδ(ω)]-[е-Jω1-е-Jω+(πе-Jω)δ(ω)]

Fп(ω)знак равно1-е-Jω1-е-Jω+π(1-е-Jω)δ(ω)

Fп(ω)знак равно1+π(1-е-Jω)δ(ω)1

ω=2kπω=2kπ

CED

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

Следовать за:

δp


Fп(ω)знак равно1(1-е-Jω)δ(ω)е(ω)ωзнак равно0е(ω)δ(ω)знак равное(0)δ(ω)е(0)знак равно0

1

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

1

так Мэтт,

е[N]

е[N]{ 12е-αNN0-12еαNN<0

α>0

Теперь у нас есть сигналы конечной энергии, и DTFT должны быть сопоставимы.

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

α0

но, увы, уже почти 2 часа ночи, и я не собираюсь сейчас с этим разбираться.


αα0
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.