Birim adım dizisi nin DTFT'sinin bu türetilmesindeki kusur nerede ?


11

Bu soru, birim adım dizisi in ayrık zamanlı Fourier dönüşümünün (DTFT) türevlerini istediğim diğer bir soru ile ilgilidir . Türev aramam sırasında inanılmaz derecede basit bir tane buldum. İlk olarak BA Shenoi tarafından bu kitabın 138. sayfasında gördüm . Ben de matematikle karşılaştım.Bu cevapta SE .u[n]

Argüman kısa ve basit olduğundan kolaylık sağlamak için burada tekrarlayacağım.

Birim adım sırası, olarak Açıkçası, DTFT'yi her iki tarafa da uygulama verir burada arasında DTFT olan . Kaynaktan elde ederiz kaynaktan ve biz DTFT almak için

(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)(1)u[n]U(ω)=F(ω)+πδ(ω)=1
(5)F(ω)=11ejω
(5)(1)u[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
burada , .- π ω < πDTFT{1}=2πδ(ω)πω<π

Denk. nin DTFT'si için kuşkusuz doğrudur. Ancak, türetme kusurludur.u [ n ](6)u[n]

Soru şudur: yukarıdaki türevdeki kusuru bulmak ve açıklamak.

Lütfen yanıtınızı spoiler etiketi ile değiştirin >!.


1
beni rahatsız eden şey, nin, sınırlı bir enerji sinyali değil, sınırlı bir güç sinyali olduğudur , bu da bu iki sonsuz enerji sinyalini bir araya getirdiğimizde elde ederiz. f[n]
robert bristow-johnson

ayrıca, mi?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
robert bristow-johnson

Yanıtlarınız için teşekkürler çocuklar! Hepsini iptal ettim ve her biri, DTFT'nin garip sinyallerin (yani veya olmayanlar) çok iyi bilinmeyen yönleri hakkında güzel bir tartışma ile sonuçlanıyor . Yalnızca bir tanesini kabul edebilirim ve yeni yanıtlar veya mevcut yanıtlardaki değişiklikler için biraz daha bekleyeceğim. Daha sonra kendi cevabımı da ekleyeceğim. 212
Matt L.

1
Mat, kesinlikle sınırlı enerji değildir . kare olacak sonsuz sayıda örnek sonlu sayı olarak . 1f[n]14
robert bristow-johnson

1
@ robertbristow-johnson: Bununla ilgili rahatsız edici ne buluyorsun? Eğer sinyaller sınırlı sayıda nokta dışında her yerde birbirini iptal ederse, o zaman elde ederiz.
Matt L.

Yanıtlar:


7

Aşağıdaki eşitliği sağlayan sonsuz sayıda sinyal vardır: Önemli olan tek şey ve sonra katsayılarının geri kalanı Denk. durumlar (yani için ardışık örneklerin çıkarılması olmalıdır ). Başka bir deyişle, Denk. herhangi bir sinyal elde edilecektir öyle ki Başka bir şekilde görmek bu, temelde olan ve bir ofset (sabit bir değer eklenmiş) olan herhangi bir işlevin

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). Bu, robert bristow-johnson tarafından cevabında yapılan açıklamayı açıklıyor : farklılaştırıcılar bu bilgileri yok ediyor (sürekli bir zamanda bir türev almak, orijinal fonksiyondaki sabit bir değerin kanıtını yok ediyor).

Özetle, ispatın kusurlu olduğuna inanıyorum çünkü takip edilen prosedür formunun ile herhangi bir fonksiyonunu kullanabilir ve bu aynı Fourier dönüşümüne sahip birçok fonksiyona yol açacaktır Fourier dönüşümü bir bijection olduğu için gerçekten yanlıştır. Belki de yazar, 'nın nin DTFT olduğunu göstermek için (en popüler kanıtı DTFT'den elde edilen) birikim özelliğine ihtiyaç duyacağının bilincinde olarak DC değerleriyle ilgili herhangi bir şeyi göz ardı etmeye karar verdi . birim adım - ergo, oldukça dairesel bir kanıt). Kanıt kesinlikle yanlış değildir , belirttiği her şey ( ve formülleriu[n]+CCRF(ω)f[n]F ( ω ) U ( ω ) F ( ω )F(ω)U(ω), birim adımın ayrışması, fark denklemi) doğrudur, ancak birikim özelliğinin neden nın herhangi bir Dirac deltasına sahip olmadığını göstermesi gerekir.F(ω)


Tamamen doğru yoldasınız! Bu kusurun nasıl çözülebileceğine dair bir fikriniz var mı, yani nasıl doğru yapılacağını?
Matt L.

@MattL. için bir başlangıç ​​koşulu ayarlamak hile yapar ve sinyali tek taraflı olarak belirler. Bu başlangıç ​​koşulu , DTFT'de bir Dirac impulsunu (birikme özelliğine göre) çarpan bir sabit olarak görünen sinyalinin DC değerini belirleyecektir . Verilen kanıtta bunun işe yaradığını düşünüyorum, çünkü sinyalinin civarında simetrik olduğu için DC değeri yoktur ve bu nedenle DTFT sadece bu durumda doğrudur. Ancak sinyalin DC'si olmadığı gerçeği, temel olarak, inanıyorum. y [ n ] f [ n ] 0y[n]y[n]f[n]0
Tendero

Çok iyi cevaplar var ve hangisinin kabul edileceğini seçmek zor. Ama bu en çok topluluk tarafından takdir edildi ve ayrıca türetmedeki hatayı en açık şekilde işaret ettiğini düşünüyorum. Hepinize teşekkür ederim!
Matt L.

4

Aldığım yanıtların sayısı beni çok etkiledi (şimdiye kadar 10 cevap!). Tabii ki, hepsi benim oyumu aldı. Bu eğlenceliydi, düşünceleriniz, yorumlarınız vb. İnsanlar bir şeyleri farklı şekilde ifade ederler ve her zaman yanlış anlamalar için yer vardır, bu yüzden bu türevdeki en önemli kusur olduğunu düşündüğüm şeyi açıkça formüle etmeye çalışacağım. Herkesin aynı fikirde olamayacağının farkındayım ve bu iyi. Bu tür ezoterik DSP konularını hepiniz gibi keskin zihinlerle tartışabildiğim için mutluyum! İşte başlıyoruz.

İlk iddiam, sorumdaki her bir denklemin doğru olduğu. Bununla birlikte, bazılarının türetilmesi ve motivasyonu tamamen yanlış ve yanıltıcıdır ve bu "türetme" ancak yazarın sonucun neye benzemesi gerektiğini bildiği için var olabilir.

Denk. (3) soruda ( ) verilen ( sorudaki Denk. dizisi için doğrudur , ancak aynı zamanda açıktır bazı keyfi sabit ile formundaki tüm diziler için doğrudur . Dolayısıyla, derivasyona göre, elde edilen DTFT , sabitinin değerine bakılmaksızın , formun tüm sekanslarının DTFT'si olmalıdır.f[n]f[n1]=δ[n]( 2 ) f [ n ] = u [ n ] + c c F ( ω ) ( 1 ) c F ( ω ) ( 5 ) ( 3 ) u [ n ] u [ n ] ( 1 )f[n](2)

(1)f[n]=u[n]+c
cF(ω)(1)c. Bu elbette anlamsız çünkü DTFT benzersizdir. Özellikle, bu çok "kanıtı" kullanarak, yı Eq. sorum (veya aşağıdaki Denk. ) aslında aradığımız nin DTFT'sidir . Peki neden yukarı bölme rahatsız denklemde verilen. sorusu?F(ω)(5)(3)u[n]u[n](1)

Bununla birlikte, tüm sekansların DTFT'lerinin Denklemi karşıladığı doğrudur . Söz konusu (kolaylık olması açısından burada tekrarlanan): Ama şimdi gerçek matematiksel kusur gelir: Gönderen Denk. sonuçlandırmak için yanlış Denk. arasında sonsuz sayıda muhtemel çözümler sadece bir tanesidir , ve uygun bir doğru son sonucu elde etmek için, yazar tarafından ihtiyaç duyulan bir olması umulur. Denk. , ile deki nin( 4 ) F ( ω ) ( 1 - e - j ω ) = 1 ( 2 ) F ( ω ) = 1(1)(4)

(2)F(ω)(1ejω)=1
(2) (3)(2)(3)f[n](1)c=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1)c=12, ancak verilen türevden bunu bilmenin bir yolu yoktur.

Peki, bu matematiksel hatayı önleyebilir ve herhangi bir sabit ile dizilerin DTFT'lerini türetmek için nasıl kullanabiliriz ? Doğru bir sonuç olup bazı henüz belirlenmemiş sabiti . Takma sol taraftaki içine verir Yani tarafından verilen tüm işlevler gerektiği gibi tatmin edin .a 1 l ( 1 ) c ( 2 ) F ( ω ) = 1(2)all(1)c(2)α(4)(2)1+α(1-e-jω)δ(ω)=1+α(1-e-jω)| ω=0δ(ω)=1+0δ

(4)F(ω)=11ejω+αδ(ω)
α(4)(2)F ( ω ) [ 4 ) ( 2 )
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2)

Sabit in değerinden belirlenebilir ile : Bu gösterilebilir ve ayrıca Wolfram alpha kabul de integrali Cauchy temel değer olduğu, olduğu kaynaktan ve elde ederiz Yani için( 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)
α=0f[n]c=0f[n]=u[n]α=πu[n] U ( ω ) = 1c=12α=0( ispatın yazarı tarafından kullanılan orijinal sırasına karşılık gelir ) ve (yani, ) nihayetinde bize nin istenen DTFT'sini verir :f[n]c=0f[n]=u[n]α=πu[n]
(9)U(ω)=11ejω+πδ(ω)


" (4) tarafından verilen tüm işlevler tatmin eder (2)", fakat "tüm işlevler tatmin eder (2) 'nin (4) formuna sahip olduğunu kanıtlamalı mıyız ?" F ( ω )F(ω)F(ω)
AlexTP

@AlexTP: Yani formun fonksiyonları sadece tatmin edici fonksiyonların bir alt kümesi olabilir ? Bu geçerli bir nokta. Ama bu sadece çünkü o, başka işlevler var olamayacağı çok açık olduğunu düşünüyorum en biz ek bir katkısı fonksiyonları ihtiyaç yüzden, sorunu nereden ile çarpılması zaman kaybolur . Bu fonksiyonlar (aslında dağılımlar) Dirac delta impulse ve türevleridir. Bununla birlikte, türevler çarpıldığında yok olmaz , bu yüzden kalan sadece Dirac delta impulsu( 2 ) ω = 0 ω = 0 ( 1 - e - j ω ) ( 1 - e - j ω )(4)(2)ω=0ω=0(1ejω)(1ejω)
Matt L.

Gerçekten bu özelliğe sahip Dirac delta dürtü (ve türevleri) dışında herhangi bir işlevi olamaz emin değilim. Ama sorun değil, cevabınız iyi yazılmış. Ben oy kullanıyorum. Teşekkürler.
AlexTP

2

Eğer Dirac Delta İşlevi olması gerekiyorsa, kusur "Açıkçası" sözcüğünü izler.

Asla sormadığım diğer sorunuz için bir cevap taslağı:

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

Kanıtın mümkün olduğunu düşünmüyorum. Bu, istenen özelliklere sahip bir "fonksiyonel tanım" durumu olabilir.

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=1
U=limN[1ejωN1ejω]
ω=0πω0
U=11ejωlimN[ejωN1ejω]
Son sınır değerine bakıyor. İçin bunun bir Dirac Delta gibi davranır olduğu açıktır. Katsayının neden olması gerektiğini bilmiyorum. Birim çemberin alanı ile ilgili olabilir. Ne zamanω=0πω0payda limitin dışına çıkarılabilir ve pay sadece birim daire boyunca atlar ve asla bir sınıra ulaşmaz. Sıfıra ayarlamak tanımlayıcı bir eylemdir.

Tanımın istenen şekilde çalıştığını kanıtlamak farklı bir konudur.

Sayfa 138 kanıtı yanlış (en azından) çünkü:

δ(t)=lima012a[u(t+a)u(ta)]=dudt
Bu, tanımlandığı şekliyle ile hiçbir şekilde benzer değildir . δ(n)=u2(n)u2(n1)

İlginç bir durum, umarım bu yardımcı olur. Söyleyeceklerini dört gözle bekliyorum.

ced


Cevabınız için teşekkürler! Ancak burada Dirac delta impulsundan değil, ayrık zamanlı dizilerden bahsettiğimizi unutmayın. Ayrık zamanda her yerde sıfıra eşittir , olduğu hariç . Yani sürekli zamanda olduğu gibi çirkin şeyler yok. Yani Denk. (3) benim sorum aslında doğru (ayrık zamanda!). δ[n]n=01
Matt L.

Cedron, bence bu denklem: ilginizi çekebilecek bu diğer söz konusu Belki bunu biraz daha düşünce veren ve isterseniz bunu orada gönderme düşünmelisiniz .. @MattL arasında.
U=11ejωlimN[ejωN1ejω]
Tendero

2

Eğer sıfıra bölmeme izin verirseniz, size olduğunu kanıtlayabilirim . dediğinizde , bir şeyi sıfıra çarpma konusunda bir sorun vardır ( ) ve ürünün eşit bir ürün olması bekleniyor.1=2

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

3
hikayenin ahlakı: farklılaştırıcılar bilgiyi yok eder. bir farklılaştırıcı ve arasındaki farkı bilmiyor . kayıp bilgileri elde etmek için sıfıra bölmeyi denemeyin. u[n]u[n]12
robert bristow-johnson

yazar kusur düzeltmek için belirtmelidir? w2πk
Yağ32

nin DTFT'si olduğunda bile herhangi bir için 1'e eşittir . δ[n] ωω=2kπ
robert bristow-johnson

neden ? ... nin kısıtlamasını yapan ' ın etki alanıdır ? δ[n]1ejwF(w)
Fat32

çünkü Matt ve elde etmek için her iki tarafı da birleştirdiğini söylediğinde (Matt'in Denk. sorta "1 elde etmek için aynı miktarda numara ile çarpın ve bölün . " ama bazen sıfıra çarpar ve böler. sıfıra çarpmak bilgiyi yok eder. bunu sıfıra bölmek bilgiyi geri getirmeyecektir.
f[n]f[n1]=δ[n]
u[n]
robert bristow-johnson

2

Denklem (4) için , değil (5). Kaçınılmadan kanıtın nasıl düzeltileceğini bilmiyorum (3).

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

Birisi iddia edebilir ne zaman var tamsayı için . F(ω)ω2kπk
robert bristow-johnson

@ robertbristow-johnson Düzenlediğiniz için teşekkürler. Yani yoktur çünkü değil, her için yakınsama yapmaz . İçin , hatta doğru formülü için bir Cauchy asıl değer olarak yorumlanmalıdırF(ω)ω=2kπn=1sin(ωn)ωω=2kπU(ω)
AlexTP

@AlexTP Neden olmadığını söylemek için bu toplama nedeninin sapması ? Birçok sinyal, "toplama DTFT'leri" birbirlerinden uzaklaşır, ancak yine de (Dirac deltaları açısından) tanımlanmış bir DTFT'ye sahiptir. Akıl yürütmenin yanlış olduğunu söylemiyorum, sadece anlamaya çalışıyorum (dağıtım teorisine çok fazla aşina değilim). F(ω)
Tendero

Toplamların geleneksel anlamda birleşmediği konusunda haklısınız, ancak yine de DTFT'ler dağıtımsal anlamda var. Sonuçta, aynı argüman bu anlamda da var olan nin DTFT'si için geçerlidir . Yani, alıştığımız ıraksama probleminden daha az matematiksel olarak karmaşık olan başka bir kusur var. u[n]
Matt L.

1
@MattL. Açık olup olmadığını bilmiyorum ama ve arasındaki farkı göremiyorum çünkü ! u[n]f[n]u[n]u[n1]=δ[n]
AlexTP

2

Sanırım bu kanıttaki kusuru ifade etmenin en iyi yolunu buldum. Bu yüzden ona bir bıçak daha vereceğim.

(1) 'de seçimi isteğe bağlıdır. Bunu ile değiştirelim . Kanıtı takip edin ve aşağıdakilerle bitirin: 12x

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

Kanıtta, olarak sınırlayan hiçbir şey yoktur, herhangi bir sonlu değer alabilir ve kanıt aynı şekilde çalışır. x12

Ayrıca, son cevabımda attığım adımı atarsanız ve bul (4)

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

Ardından (5) ve (6) 'yı ekleyerek şunları elde edersiniz:

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

Daha önce de belirttiğim gibi, oraya ulaşmak tanımıyla tutarsız.

Bu kanıt i gösterememektedir ve herhangi bir tanımlanmış x için tutarsız sonuçların geleceğini göstermektedir. Bu yüzden ilk cevabımdaki ifadeye geri dönüp katsayısı için değerinin matematiksel bir gerçekçilik değil tanımlayıcı bir eylem olduğunu görüyorum. x=12πδ(ω)

Belki de i doğru değer yapan başka bir durum daha vardır , ancak bu kanıt bunu sağlamaz. x=12

ced


1

Bu benim ilk cevabımdaki yorumlara yanıt olarak. Spoiler gizleme yüzünden ayrı bir cevap olarak gönderiyorum.

Diğer cevabımı diğer soruya da gönderecektim, ancak bu alandaki deneyim eksikliğimden dolayı olmadı. Dün yayınladım, sildim, sonra silindim, sonra spoiler etiketlerini nasıl kullanacağımı anladım.

Açıkça sorunda tanımlanan işlevi Dirac Delta işlevi değildir. Wikipedia'da DTFT'yi aradım ve Dirac Delta fonksiyonu için DTFT de bir tane. Arayacağım problem . δδδp

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

Sol ve sağ parçaların DTFT'sinin alınması. Gösterimlerin doğru olduğundan emin değilim, ama matematik açık olmalı. Kanıtlanan tanımı kullanarak.

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

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

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

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

Bu nedenle, (4) 'ün RHS değeri hariç . [Düzenle: Doh, Dirac Delta, bu yüzden bu ifade yanlış. Sanırım "undefined" dışında doğru olmalıdır . Gerçek Analiz benim en sevdiğim matematikti. Bunu şimdi yalnız bırakıyorum.] ω=2kπω=2kπ

ced

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

Takip et:

O bir DTFT açıktır bir DTFT tanımına takıldığında 1 olmalıdır. Bu nedenle, kanıtlanacak tanımı kullanırken farklı bir cevap aldığım için kanıtlanacak tanımın doğru olmadığı anlamına gelir (matematiksel anlamda). Ayrıca, düzeltmeyi kanıtın sonuna kadar götürürseniz farklı bir tanıma ulaşırsınız. İddianın doğru olduğunu varsaymak, yanlış olduğunu kanıtlamak için kullanılır.δp


Aslında terimi sıfıra eşit olduğu için olduğunu gösterdiniz . sürekli olan herhangi bir fonksiyonu için ve (burada durum budur), tüm terim kaybolur. Fp(ω)=1(1ejω)δ(ω)f(ω)ω=0f(ω)δ(ω)=f(0)δ(ω)f(0)=0
Matt L.

1

Bana göre, (3) ve (4) arasında bir ilk kusur ortaya çıkıyor: bu, klasik dikkatsiz integral / sonsuz toplam bölünmesinin bir örneğidir. izin vermek için koşullar gereklidir: Standart veya koşulları yeterince keskin olmayabilir. Bu durum, bir şekilde buraya, ilgili olabilir için, Fubini türetme teoremi veya: sonsuz toplamı ve kesikli türevi alışverişinde bulunabildiği zaman? Etraftaki yollar monotonluk veya Cesaro benzeri meblağlar etrafında dönebilir, ancak bunu daha uzun düşüneceğim.1 2 f [ n ] - f [ n - 1 ]

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

1

yani Matt,

Neden güç sinyallerini enerji sinyallerine göre karşılaştırmanın sorunlu olduğunu düşünmüyorsunuz, ancak tanımını biraz değiştirdiğimizi varsayalım :f[n]

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

bazı .α>0

şimdi sonlu enerji sinyallerimiz var ve DTFT'lerin karşılaştırılabilir olması gerekiyor.

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

DTFT'lerin ne olduğunu merak ediyorum? ve sonra bıraktığımızda ne olur ? Ben hala bir sorun olduğunu bilgileri (ve frekans alanında 0 ile çarparak bilginin karşılık gelen imha) farklılaştırıcılar sorunu olduğunu düşünüyorum. ama belki aynı Hilbert alanını paylaşmayan sinyal sınıflarını karşılaştırma problemini kaybedebiliriz.α0

ama, ne yazık ki, neredeyse 2 am ve şimdi onunla uğraşmayacağım.


Bu iyi bir şeydir ve bu olmayan sinyallerin DTFT'sini sınırını kullanarak hesaplamak için bir seçenektir . Deneyin ve başarılı olacağınızdan eminim, ama acı verici. Aynı sonucu almanın daha kolay yolları vardır. Verilen kanıt IMHO çalışacak şekilde değiştirilebilir (cevabımı görün). α 0αα0
Matt L.
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.