Pythonda raqamli signalni qayta ishlash dsp haqida o'ylang

7-bob. Diskret Furye transformatsiyasi
harmonikalar yana o'sishni boshlaydi va 9500 Gts ga etadi. Nima gaplar?
qiladi.
Ammo raqamning o'ng yarmi ajablanib. 5000 Gts dan o'tgan, amplitudasi
DFT(y)[n] = DFT(e)[n] + exp(ÿ2pin/N)DFT(o)[n]
FFT kaliti Danielson-Lanzos lemmasi:
7.10 Mashqlar
Machine Translated by Google

---

FFT haqida ko'proq ma'lumotni https://en.wikipedia.org/wiki/ Fast_Fourier_transform sahifasida
o'qishingiz mumkin.
Bunday holda, DFT(y) = y. Yoki y uzunligi etarlicha kichik bo'lsa, siz uning DFT ni matritsani
ko'paytirish orqali, ehtimol oldindan hisoblangan matritsadan foydalanib hisoblashingiz mumkin.
Ushbu rekursiyaning asosiy holati uchun y uzunligi 1 bo'lguncha kutishingiz mumkin.
elementlar, o.
7.10. Mashqlar
Maslahat: Men sizga ushbu algoritmni rekursiv bo'lmagan versiyadan boshlab bosqichma-bosqich
amalga oshirishni taklif qilaman. 2-bosqichda rekursiv qo'ng'iroq qilish o'rniga, 7.7-bo'limda
ta'riflanganidek dft yoki np.fft.fft dan foydalaning. 3-bosqichni ishga tushiring va natijalar boshqa
ilovalar bilan mos kelishini tasdiqlang.
2. Rekursiv chaqiruvlarni amalga oshirish orqali e va o larning DFT larini hisoblang.
89
Keyin asosiy qutini qo'shing va uning ishlashini tasdiqlang. Nihoyat, 2-bosqichni rekursiv qo'ng'iroqlar
bilan almashtiring.
3. Danielson-Lanzos yordamida n ning har bir qiymati uchun DFT(y) ni hisoblang
1. To‘lqin massivi berilgan y, uni juft elementlarga, e va toq elementlarga ajrating.
Yana bir maslahat: DFT davriy ekanligini unutmang; np.tile ni foydali deb topishingiz mumkin.
lemma.
Machine Translated by Google

---

90
7-bob. Diskret Furye transformatsiyasi
Machine Translated by Google

---