Pythonda raqamli signalni qayta ishlash dsp haqida o'ylang

92
fitna.
8-bob. Filtrlash va konvolyutsiya
Endi men birinchi 11 elementning o'rtacha qiymatini ko'paytirish orqali hisoblay olaman
oyna /= yig'indisi (oyna)
to'lqin - signalning 1 soniyali qismi; segment men foydalanadigan qisqaroq bo'lakdir
segment = to'lqin.segment(davomiylik=0,01)
N
a
rxi
($
)
kundalik yopilish
O'rtacha 30 kun
Machine Translated by Google

---

mahsulot = to'ldirilgan * ys
def silliq (ys, oyna):
o'rtacha.
qaysi birlarni o'ngga siljitadi va oxiridan nollardan birini o'radi
N = len(ys)
Harakatlanuvchi o'rtachaning keyingi elementini hisoblash uchun biz oynani aylantiramiz ,
8.2-rasm: 400 Gts chastotali kvadrat signal (kulrang) va 11 elementli harakatlanuvchi
massiv.
Bu elementlarning barchasi -1 bo'lgani uchun ularning o'rtacha ko'rsatkichlari -1 ga teng.
Qolgan elementlarni ham xuddi shunday hisoblashimiz mumkin. Quyidagi funktsiya biz
ko'rgan kodni tsiklga o'rab oladi va natijalarni saqlaydi
Natija yana -1 ga teng.
elementwise mahsuloti - massivning dastlabki 11 elementining o'rtacha qiymati.
prod oyna va to'lqin massivining mahsulotidir. yig'indisi
summa (mahsulot)
mahsulot = o'ralgan * ys
segment.ys bilan bir xil uzunlik. Bu kabi nollarni qo'shish to'ldirish deb ataladi .
93
padded - oxiriga nollar qo'shilgan oynaning versiyasi, shuning uchun u
8.1. Silliqlash
o'ralgan = np.roll(o'ralgan, 1)
summa (mahsulot)
Oyna va to'lqin massivini ko'paytirganda biz ikkinchisidan boshlab to'lqin massivining
keyingi 11 elementining o'rtacha qiymatini olamiz.
boshigacha.
1.0
0,5
1.0
Vaqt(lar)
0,006
0,5
0,010
0,004
0,0
0,008
0,002
0.000
Machine Translated by Google

---

silliqlangan = np.zero(N)
to'ldirilgan = thinkdsp.zero_pad(oyna, N) o'ralgan
= to'ldirilgan
diapazondagi i
uchun(N): tekislangan[i] =
yig‘indisi(o‘ralgan * ys) o‘ralgan = np.roll(o‘ralgan, 1)
Loop ichida biz ys ni rulonga ko'paytiramiz va 11 ta elementni tanlaymiz va ularni
qo'shamiz.
konvolved = np.convolve(ys, window, mode='valid') silliq2 =
thinkdsp.Wave(convolved, framerate=wave.framerate)
94
Konvolyutsiya shu qadar keng tarqalgan operatsiyaki, NumPy mening versiyamdan
ko'ra soddaroq va tezroq bo'lgan dasturni taqdim etadi:
silliqlangan - natijalarni o'z ichiga olgan massiv; padded - uzunligi N bo'lishi uchun
oyna va etarli nollarni o'z ichiga olgan massiv; va rulonli - padded nusxasi bo'lib, har
safar halqa bo'ylab bir element bilan o'ngga siljiydi.
Joriy rejim bayrog'i faqat oyna va to'lqinlar massivi to'liq bir-biriga mos tushganda
qiymatlarni hisoblashi kerakligini ko'rsatadi, shuning uchun oynaning o'ng qirrasi
to'lqin qatorining oxiriga yetganda to'xtaydi. Bundan tashqari, natija 8.2-rasmdagi
bilan bir xil.