Ayrık zamanlı Fourier dönüşümü - Discrete-time Fourier transform

Fourier dönüşümleri
Sürekli Fourier dönüşümü
Fourier serisi
Ayrık zamanlı Fourier dönüşümü
Ayrık Fourier dönüşümü
Bir halka üzerinde ayrık Fourier dönüşümü
Sonlu gruplar üzerinde Fourier dönüşümü
Fourier analizi
İlgili dönüşümler

İçinde matematik, ayrık zamanlı Fourier dönüşümü (DTFT) bir biçimdir Fourier analizi bu bir dizi değer için geçerlidir.

DTFT genellikle sürekli bir fonksiyonun örneklerini analiz etmek için kullanılır. Dönem ayrık zaman Dönüşümün, genellikle aralığı zaman birimlerine sahip olan örnekler olan ayrı veriler üzerinde çalıştığı gerçeğini ifade eder. Düzgün aralıklı örneklerden bir frekans işlevi üretir. periyodik toplama of sürekli Fourier dönüşümü orijinal sürekli işlevin. Belirli teorik koşullar altında, örnekleme teoremi orijinal sürekli fonksiyon, DTFT'den ve dolayısıyla orijinal ayrık örneklerden mükemmel şekilde kurtarılabilir. DTFT'nin kendisi sürekli bir frekans fonksiyonudur, ancak bunun ayrık örnekleri kolayca hesaplanabilir. ayrık Fourier dönüşümü (DFT) (bkz. § DTFT'yi Örnekleme ), modern Fourier analizinin açık ara en yaygın yöntemidir.

Her iki dönüşüm de tersine çevrilebilir. Ters DTFT, orijinal örneklenmiş veri dizisidir. Ters DFT, orijinal dizinin periyodik bir toplamıdır. hızlı Fourier dönüşümü (FFT), DFT'nin bir döngüsünü hesaplamak için bir algoritmadır ve tersi, ters DFT'nin bir döngüsünü üretir.

Tanım

Ayrık gerçek veya karmaşık sayılar kümesinin ayrık zamanlı Fourier dönüşümü x[n], hepsi için tamsayılar n, bir Fourier serisi, bir frekans değişkeninin periyodik bir fonksiyonunu üreten. Frekans değişkeni ω, normalleştirilmiş birimler nın-nin radyan / örnek, periyodiklik ve Fourier serisi:[1]:s. 147

 

 

 

 

(Denklem.1)

Bu frekans etki alanı işlevinin faydası, Poisson toplama formülü. İzin Vermek X(f) herhangi bir fonksiyonun Fourier dönüşümü olabilir, x(t), belirli aralıklarla örnekleri T (saniye) eşittir (veya orantılıdır) x[n] dizi, yani Tx(nT) = x[n]. O halde, Fourier serisi tarafından temsil edilen periyodik fonksiyon, periyodik bir toplamıdır. X(f) sıklık açısından f içinde hertz (döngü / saniye):[a]

 

 

 

 

(Denklem.2)

Şekil 1. Bir Fourier dönüşümünün (sol üst) ve sol alt köşede periyodik toplamının (DTFT) tasviri. Sağ alt köşe, ayrık bir Fourier dönüşümü (DFT) ile hesaplanan DTFT örneklerini gösterir.

Tamsayı k birimleri var döngü / örnek, ve 1/T örnekleme oranı, fs (örnek / saniye). Yani X1/T(f) tam kopyalarından oluşur X(f) katları tarafından kaydırılan fs hertz ve ekleme ile birleştirilir. Yeterince büyük fs k = 0 bölgede dönem gözlemlenebilir [−fs/ 2, fs/2] az veya hiç distorsiyon olmadan (takma ad ) diğer terimlerden. Şekil 1'de, sol üst köşedeki dağılımın uçları, periyodik toplamda (sol alt) örtüşme ile maskelenmiştir.

Ayrıca şunu da not ediyoruz ei2πfTn Fourier dönüşümüdür δ(tnT). Bu nedenle, DTFT'nin alternatif bir tanımı şöyledir:[A]

 

 

 

 

(Denklem 3)

Modüle edilmiş Dirac tarağı işlev, bazen şu şekilde ifade edilen matematiksel bir soyutlamadır dürtü örneklemesi.[2]

Ters dönüşümü

DTFT işlevinden ayrık veri dizisini kurtaran bir işleme, ters DTFT. Örneğin, her iki tarafın ters sürekli Fourier dönüşümü Denklem 3 diziyi modüle edilmiş bir Dirac tarak işlevi biçiminde üretir:

Ancak bunu not ederek X1/T(f) periyodiktir, gerekli tüm bilgiler herhangi bir uzunluk aralığında bulunur 1/T. Hem de Denklem.1 ve Denklem.2, n üzerindeki toplamlar a Fourier serisi katsayılarla x[n]. Fourier katsayılarının standart formülleri aynı zamanda ters dönüşümlerdir:

 

 

 

 

(Denklem.4)

Periyodik veriler

Giriş veri dizisi x[n] dır-dir N-periyodik, Denklem.2 hesaplamalı olarak ayrık bir Fourier dönüşümüne (DFT) indirgenebilir, çünkü:

  • Mevcut tüm bilgiler içinde bulunur N örnekler.
  • X1/T(f) tam sayı katları dışında her yerde sıfıra yakınsar 1/(NT), olarak bilinir harmonik frekanslar. Bu frekanslarda, DTFT farklı frekansa bağlı oranlarda sapar. Ve bu oranlar, bir döngünün DFT'si tarafından verilmektedir. x[n] sıra.
  • DTFT periyodiktir, bu nedenle maksimum benzersiz harmonik genlik sayısı (1/T) / (1/(NT)) = N

DFT katsayıları şu şekilde verilir::

ve DTFT:
      [b]

Bu ifadeyi ters dönüşüm formülüne dönüştürmek,:

(tüm tam sayılar )

beklenildiği gibi. Yukarıdaki satırdaki ters DFT, bazen bir Ayrık Fourier serileri (DFS).[1]:s 542

DTFT'yi örnekleme

DTFT sürekli olduğunda, yaygın bir uygulama, rastgele sayıda örneği hesaplamaktır (N) periyodik fonksiyonun bir döngüsü X1/T: [1]:s. 557–559 ve 703

nerede bir periyodik toplama:

(görmek Ayrık Fourier serileri )

dizi ters DFT'dir. Böylece, DTFT örneklememiz ters dönüşümün periyodik olmasına neden olur. Dizisi |Xk|2 değerler olarak bilinir periodogramve parametre N aynı isimli Matlab fonksiyonunda NFFT olarak adlandırılır.[3]

Bir döngüyü değerlendirmek için sayısal olarak, sonlu uzunlukta bir x[n] sıra. Örneğin, uzun bir dizi, bir pencere işlevi uzunluk L özel olarak anılmaya değer üç davayla sonuçlandı. Notasyonel basitlik için, x[n] pencere işlevi tarafından değiştirilen değerleri temsil etmek için aşağıdaki değerler.

Durum: Frekans katsayısı. L = Nben, bir tam sayı için ben (tipik olarak 6 veya 8)

Bir döngü bir toplamına indirgenir ben uzunluk segmentleri N. DFT daha sonra çeşitli isimlerle anılır, örneğin:

  • çok fazlı DFT[8][9]
  • çok fazlı filtre bankası[11]
  • çoklu blok pencereleme ve zaman örtüşme.[12]

Örneklenen verilerin tek bir alanda (zaman veya sıklık) kesilmesinin örtüşme oluşturduğunu hatırlayın (bazen takma ad ) diğerinde ve tam tersi. İle karşılaştırıldığında L-uzunluk DFT, toplama / üst üste binme, sıklıkta azalmaya neden olur,[1]:s. 558 sadece en az etkilenen DTFT örneklerini bırakarak spektral sızıntı. Bu genellikle bir FFT uygularken bir önceliktir filtre bankası (kanal oluşturucu). Geleneksel bir pencere uzunluğu fonksiyonu ile L, taraklanma kaybı kabul edilemez. Böylece çok bloklu pencereler kullanılarak oluşturulur FIR filtresi tasarım araçları.[13][14] Frekans profilleri en yüksek noktada düzdür ve kalan DTFT örnekleri arasındaki orta noktada hızla düşer. Parametrenin değeri ne kadar büyükse benpotansiyel performans o kadar iyi olur.

Durum: L = N+1.

Simetrik olduğunda, L-uzunluk pencere işlevi () 1 katsayı ile kesilir, denir periyodik veya DFT-çift. Kesme, DTFT'yi etkiler. Kesilmiş dizinin bir DFT'si, DTFT'yi aşağıdaki frekans aralıklarında örnekler. 1/N. Örneklemek için aynı frekanslarda, karşılaştırma için, DFT periyodik toplamanın bir döngüsü için hesaplanır, [D]

Şekil 2. DFT ei2πn / 8 için L = 64 ve N = 256
Şekil 3. DFT ei2πn / 8 için L = 64 ve N = 64

Durum: Frekans enterpolasyonu. LN

Bu durumda, DFT daha tanıdık bir biçime sadeleştirilir:

DFT'yi hesaplamak için hızlı bir Fourier dönüşüm algoritmasından yararlanmak için, toplama genellikle tüm N şartlar olsa bile NL bunlardan sıfır. Bu nedenle durum L < N genellikle şu şekilde anılır sıfır dolgu.

Spektral sızıntı, L azalırsa, birden çok frekans bileşeninin çözünürlüğü ve her DTFT örneği tarafından ölçülen gürültü miktarı gibi bazı önemli performans ölçütleri için zararlıdır. Ancak bu şeyler her zaman önemli değildir, örneğin x[n] dizi, bir pencere işlevi tarafından şekillendirilen gürültüsüz bir sinüzoiddir (veya sabittir). O zaman kullanmak yaygın bir uygulamadır sıfır dolgu pencere işlevlerinin ayrıntılı sızıntı modellerini grafiksel olarak görüntülemek ve karşılaştırmak. Dikdörtgen bir pencere için bunu göstermek için şu sırayı göz önünde bulundurun:

ve

Şekil 2 ve 3 etiketlerinde belirtildiği gibi, iki farklı boyutlu DFT'nin büyüklüğünün çizimidir. Her iki durumda da, baskın bileşen sinyal frekansındadır: f = 1/8 = 0.125. Ayrıca şurada da görülebilir İncir. 2 spektral sızıntı modelidir L = 64 dikdörtgen pencere. İçinde illüzyon Şek. 3 DTFT'yi sadece sıfır geçişlerinde örneklemenin bir sonucudur. Sonlu uzunlukta bir dizinin DTFT'sinden ziyade, sonsuz uzunlukta bir sinüzoidal dizi izlenimi verir. İllüzyona katkıda bulunan faktörler, dikdörtgen bir pencerenin kullanılması ve 64 örnek başına tam olarak 8 (bir tam sayı) döngü ile bir frekans (1/8 = 8/64) seçimidir. Bir Hann penceresi tepe noktasının 3 örneğe genişletilmesi dışında benzer bir sonuç üretecektir (bkz. DFT-hatta Hann penceresi ).

Evrişim

evrişim teoremi diziler için:

[16]:s. 297[c]

Önemli bir özel durum, dairesel evrişim dizilerin x ve y tarafından tanımlandı nerede periyodik bir toplamdır. Ayrık frekans doğası sürekli işlevi olan ürünün ayrıca ayrıktır, bu da ters dönüşümün önemli ölçüde basitleştirilmesine neden olur:

[17][1]:s. 548

İçin x ve y sıfır olmayan süresi şuna eşit veya daha az olan diziler Nson bir basitleştirme ise:

Bu sonucun önemi şu adreste açıklanmıştır: Dairesel evrişim ve Hızlı evrişim algoritmaları.

Simetri özellikleri

Karmaşık bir işlevin gerçek ve hayali kısımları, çift ​​ve tek parçalar Aşağıda RE, RO, IE ve IO alt simgeleriyle gösterilen dört bileşen vardır. Ve karmaşık bir zaman fonksiyonunun dört bileşeni ile karmaşık frekans dönüşümünün dört bileşeni arasında bire bir eşleştirme vardır.:[16]:s. 291

Bundan çeşitli ilişkiler, örneğin:

  • Gerçek değerli bir fonksiyonun dönüşümü (xYENİDEN+ xRO) hatta simetrik işlevi XYENİDEN+ i XIO. Tersine, eşit simetrik bir dönüşüm, gerçek değerli bir zaman alanını ifade eder.
  • Hayali değerli bir fonksiyonun dönüşümü (ben xIE+ i xIO) garip simetrik işlevi XRO+ i XIEve sohbet doğrudur.
  • Eşit simetrik bir fonksiyonun dönüşümü (xYENİDEN+ i xIO) gerçek değerli fonksiyondur XYENİDEN+ XROve sohbet doğrudur.
  • Garip simetrik bir fonksiyonun dönüşümü (xRO+ i xIE) hayali değerli bir fonksiyondur i XIE+ i XIOve sohbet doğrudur.

Z-dönüşümü ile ilişki

bir Fourier serisi bu aynı zamanda ikili olarak da ifade edilebilir Z-dönüşümü. Yani:

nerede gösterim, Z-dönüşümünü Fourier dönüşümünden ayırır. Bu nedenle, Z-dönüşümünün bir kısmını Fourier dönüşümü cinsinden de ifade edebiliriz:

Unutmayın ki parametre T değişiklikleri, şartları sabit bir ayrılık olarak kalmak ayrı ve genişlikleri yukarı veya aşağı ölçeklenir. Şartları X1/T(f) sabit bir genişlik ve ayrılıkları kalır 1/T yukarı veya aşağı ölçeklenir.

Ayrık zamanlı Fourier dönüşümleri tablosu

Bazı yaygın dönüşüm çiftleri aşağıdaki tabloda gösterilmektedir. Aşağıdaki gösterim geçerlidir:

  • sürekli açısal frekansı temsil eden gerçek bir sayıdır (örnek başına radyan cinsinden). ( döngü / saniye cinsindendir ve saniye / örnek cinsindendir.) Tablodaki tüm durumlarda, DTFT 2π-periyodiktir ( ).
  • tanımlanmış bir işlevi belirtir .
  • tanımlanmış bir işlevi belirtir ve başka yerde sıfır. Sonra:
  • ... Dirac delta işlevi
  • normalleştirilmiş mi sinc işlevi
  • ... dikdörtgen işlevi
  • ... üçgen işlevi
  • n ayrık zaman alanını temsil eden bir tamsayıdır (örneklerde)
  • ayrık zamandır birim adım işlevi
  • ... Kronecker deltası
Zaman alanı
x[n]
Frekans alanı
X(ω)
UyarılarReferans
[16]:s. 305
tamsayı

garip M
hatta M

tamsayı

terim olarak yorumlanmalıdır dağıtım anlamında Cauchy ana değeri çevresinde kutuplar -de .
[16]:s. 305

gerçek Numara

gerçek Numara ile
gerçek Numara ile
tamsayılar ve
gerçek sayılar ile
gerçek Numara ,
olarak çalışır farklılaştırıcı filtre
gerçek sayılar ile
Hilbert dönüşümü
Trapezoid signal.svggerçek sayılar
karmaşık

Özellikleri

Bu tablo, zaman alanındaki bazı matematiksel işlemleri ve frekans alanındaki karşılık gelen etkileri gösterir.

EmlakZaman alanı
x[n]
Frekans alanı
UyarılarReferans
DoğrusallıkKarışık sayılar [16]:s. 294
Zamanın tersine çevrilmesi / Frekansın tersine çevrilmesi[16]:s. 297
Zaman konjugasyonu[16]:s. 291
Zamanın tersine çevrilmesi ve konjugasyon[16]:s. 291
Zamanın gerçek kısmı[16]:s. 291
Zamanın hayali bölümü[16]:s. 291
Frekansta gerçek kısım[16]:s. 291
Frekansta hayali kısım[16]:s. 291
Zamanda kayma / frekansta modülasyontamsayı k[16]:s. 296
Frekansta kayma / Zaman içinde modülasyongerçek Numara [16]:s. 300
Decimation  [E]tamsayı
Zaman Genişletmetamsayı [1]:s. 172
Frekansta türev[16]:s. 303
Frekans entegrasyonu
Zaman farklılığı
Zaman içinde toplama
Zamanda evrişim / frekansta çarpma[16]:s. 297
Zamanda çarpma / frekansta evrişimPeriyodik evrişim[16]:s. 302
Çapraz korelasyon
Parseval teoremi[16]:s. 302

Ayrıca bakınız

Notlar

  1. ^ Aslında Denklem.2 genellikle şu şekilde gerekçelendirilir:[1]:s sayfa 143
  2. ^ WOLA ile karıştırılmamalıdır Örtüşme ekleme yöntemi parçalı evrişim.
  3. ^ WOLA örneği: Dosya: WOLA channelizer example.png
  4. ^ Bir örnek şekil DTFT'yi örnekleme. Gerçek değerli DFT örnekleri aşağıdakilerin bir sonucudur: DFT-eşit simetri[15]:s sayfa 52
  5. ^ Bu ifade şu şekilde türetilmiştir:[1]:s sayfa 168

Sayfa alıntıları

  1. ^ Oppenheim ve Schafer, p 147 (4.20), p 694 (10.1) ve Prandoni ve Vetterli, p 255, (9.33), burada:   ve
  2. ^ Oppenheim ve Schafer, p 551 (8.35) ve Prandoni ve Vetterli, s 82, (4.43), nerede:      ve
  3. ^ Oppenheim ve Schafer, s 60, (2.169) ve Prandoni ve Vetterli, s 122, (5.21)

Referanslar

  1. ^ a b c d e f g h Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). "4.2, 8.4". Ayrık zamanlı sinyal işleme (2. baskı). Upper Saddle River, NJ: Prentice Hall. ISBN  0-13-754920-2. Bir periyodik olmayan sekans x [n] 'nin Fourier dönüşümünün örnekleri, x [n]' nin periyodik kopyalarının toplanmasıyla elde edilen periyodik bir sekansın DFS katsayıları olarak düşünülebilir. url =https://d1.amobbs.com/bbs_upload782111/files_24/ourdev_523225.pdf
  2. ^ Rao, R. (2008). Sinyaller ve Sistemler. Prentice-Hall Of India Pvt. Sınırlı. ISBN  9788120338593.
  3. ^ "Periodogram güç spektral yoğunluk tahmini - MATLAB periodogram".
  4. ^ Gumas, Charles Constantine (Temmuz 1997). "Window-presum FFT, yüksek dinamik aralık, çözünürlük sağlar". Kişisel Mühendislik ve Enstrümantasyon Haberleri: 58–64. 2001-02-10 tarihinde orjinalinden arşivlendi.CS1 bakimi: BOT: orijinal url durumu bilinmiyor (bağlantı)
  5. ^ Crochiere, R.E .; Rabiner, L.R. (1983). "7.2". Çok Oranlı Dijital Sinyal İşleme. Englewood Kayalıkları, NJ: Prentice-Hall. sayfa 313–326. ISBN  0136051626.
  6. ^ Wang, Hong; Lu, Youxin; Wang, Xuegang (16 Ekim 2006). "WOLA Filterbank ile Kanallaştırılmış Alıcı". 2006 CIE Uluslararası Radar Konferansı. Şangay, Çin: IEEE: 1–3. doi:10.1109 / ICR.2006.343463. ISBN  0-7803-9582-4.
  7. ^ Lyons, Richard G. (Haziran 2008). "DSP Püf Noktaları: Pratik bir spektrum analizörü oluşturma". EE Times. Alındı 2020-02-20. Bununla birlikte, etiketli bir bağlantı içerdiğini unutmayın. ağırlıklı örtüşme-ekleme yapısı yanlış giden Örtüşme ekleme yöntemi.
  8. ^ a b Lillington, John (Mart 2003). "Geniş Bant Kanallaştırma Mimarilerinin Karşılaştırması" (PDF). Dallas: Uluslararası Sinyal İşleme Konferansı. s. 4 (şekil 7). Alındı 2020-09-06. "Weight Overlap and Add" veya WOLA veya onun alt grubu "Polyphase DFT", daha yerleşik hale geliyor ve büyük, yüksek kaliteli filtre kümelerinin gerekli olduğu yerlerde kesinlikle çok verimli.
  9. ^ a b Lillington, John. "Filtre Bankası Tekniklerinin Gözden Geçirilmesi - RF ve Dijital" (PDF). armms.org. Wight Adası, İngiltere: Libra Design Associates Ltd. s. 11. Alındı 2020-09-06. Neyse ki, aşağıdaki Şekil 20'de gösterildiği gibi, Çok Fazlı veya WOLA (Ağırlık, Örtüşme ve Ekleme) FFT olarak bilinen çok daha zarif bir çözüm var.
  10. ^ Hochgürtel, Stefan (2013). "Yüksek çözünürlüklü geniş bantlı FFT-spektrometrelerin verimli uygulamaları ve bunların APEX Galaktik Merkez hat araştırmalarına uygulanması" (PDF). hss.ulb.uni-bonn.de. Bonn: Rhenish Friedrich Wilhelms Bonn Üniversitesi. s. 26–27. Alındı 2020-09-06. N noktalı DFT için M-katlama WOLA gerçekleştirmek için, M · N gerçek giriş örnekleri aj ilk olarak bir pencere işlevi ile çarpılır wj aynı büyüklükte
  11. ^ Chennamangalam, Jayanth (2016-10-18). "Çok Fazlı Filtre Bankası Tekniği". CASPER Grubu. Alındı 2016-10-30.
  12. ^ Dahl, Jason F. (2003-02-06). Spektrum Tahminin Zaman Takma Yöntemleri (Doktora). Brigham Young Üniversitesi. Alındı 2016-10-31.
  13. ^ Lin, Yuan-Pei; Vaidyanathan, P.P. (Haziran 1998). "Kosinüs Modüle Edilmiş Filtre Bankalarının Prototip Filtrelerinin Tasarımı için Kaiser Pencere Yaklaşımı" (PDF). IEEE Sinyal İşleme Mektupları. 5 (6): 132–134. Bibcode:1998ISPL .... 5..132L. doi:10.1109/97.681427. Alındı 2017-03-16.
  14. ^ Harris, Frederic J. (2004-05-24). "9". İletişim Sistemleri için Çok Oranlı Sinyal İşleme. Upper Saddle River, NJ: Prentice Hall PTR. s. 226–253. ISBN  0131465112.
  15. ^ Harris, Fredric J. (Ocak 1978). "Ayrık Fourier Dönüşümü ile Harmonik Analiz için Windows'un Kullanımı hakkında" (PDF). IEEE'nin tutanakları. 66 (1): 51–83. Bibcode:1978IEEEP.66 ... 51H. CiteSeerX  10.1.1.649.9880. doi:10.1109 / PROC.1978.10837.
  16. ^ a b c d e f g h ben j k l m n Ö p q r Proakis, John G .; Manolakis, Dimitri G. (1996). Sayısal Sinyal İşleme: İlkeler, Algoritmalar ve Uygulamalar (3 ed.). New Jersey: Prentice-Hall Uluslararası. Bibcode:1996dspp.book ..... P. ISBN  9780133942897. sAcfAQAAIAAJ.
  17. ^ Rabiner, Lawrence R.; Altın, Bernard (1975). Dijital sinyal işleme teorisi ve uygulaması. Englewood Cliffs, NJ: Prentice-Hall, Inc. s. 59 (2.163). ISBN  978-0139141010.
  1. Prandoni, Paolo; Vetterli, Martin (2008). İletişim için Sinyal İşleme (PDF) (1 ed.). Boca Raton, FL: CRC Press. s. 72, 76. ISBN  978-1-4200-7046-0. Alındı 4 Ekim 2020. periyodik sinyal için DFS katsayıları, DTFT'si için ayrı bir değer kümesidir.

daha fazla okuma

  • Porat, Boaz (1996). Dijital Sinyal İşleme Kursu. John Wiley and Sons. sayfa 27–29 ve 104–105. ISBN  0-471-14961-6.
  • Siebert, William M. (1986). Devreler, Sinyaller ve Sistemler. MIT Elektrik Mühendisliği ve Bilgisayar Bilimleri Serisi. Cambridge, MA: MIT Press. ISBN  0262690950.
  • Lyons Richard G. (2010). Dijital Sinyal İşlemeyi Anlamak (3. baskı). Prentice Hall. ISBN  978-0137027415.