Pythonda raqamli signalni qayta ishlash dsp haqida o'ylang

Haqiqiy signallarning 7.9 DFT
Machine Translated by Google

---

Dominant chastota 500 Gts da, garmoniklar 1/ f ga o'xshab tushadi.
Bu nuqtadan keyin hech qanday qo'shimcha ma'lumot yo'q, biz DFT ning faqat
birinchi yarmini baholab, vaqtni tejashimiz mumkin va aynan shu narsa np.fft.rfft
Haqiqiy signalning DFT katlama chastotasi atrofida simmetrikdir. beri
0 dan 10 kHz gacha. Shaklning chap yarmi biz kutishimiz kerak bo'lgan narsadir: the
,
ko'p ilovalar uchun, lekin tezroq algoritm mavjud, Tez Furye Trans shakli (FFT), bu
N log N ga proportsional vaqtni oladi.
7.3-rasmda har bir chastota komponenti uchun signalning amplitudasi ko'rsatilgan
6000 Gts chastotasi 4000 Gts chastotasi bilan bir xil va hokazo.
5500 Gts da DFT, biz 4500 Gts bilan bir xil qiymatga ega bo'lamiz. Xuddi shunday, qiymat
ylabel='amplituda')
bu erda N - to'lqinlar qatorining uzunligi. Bu yetarlicha tez
4500 Gts chastotadagi komponentdan farq qilmaydi. Biz baholaganimizda
Thinkplot.config(xlabel='chastota (Hz)',
7.2 -mashq Ushbu bobda men DFT va matritsalarni ko'paytirish sifatida DFT oyatida
qanday ifodalash mumkinligini ko'rsatdim. Bu operatsiyalar N2 ga mutanosib vaqt
oladi
Mashq 7.1 Ushbu bob uchun daftar, chap07.ipynb, qo'shimcha misollar va
tushuntirishlarni o'z ichiga oladi. Uni o'qing va kodni ishga tushiring.
Thinkplot.plot(fs, amper)
chastotasi 5000 Hz. 2.3-bo'limda ko'rganimizdek, 5500 Gts chastotada komponent
Ushbu lemma DFT uchun rekursiv algoritmni taklif qiladi:
Javob: taxallus. Shuni esda tutingki, kvadrat tezligi 10000 Gts bilan katlama
88
Ushbu mashqlarning yechimlari chap07soln.ipynb da.
Bu yerda DFT(y)[n] y ning DFT ning n-elementi; e - y ning juft elementlarini o'z
ichiga olgan to'lqinli massiv, o esa y ning toq elementlarini o'z ichiga oladi.