Pythonda raqamli signalni qayta ishlash dsp haqida o'ylang

Avtokorrelyatsiya
5-bob
5.1 Korrelyatsiya
Nuxsy
Machine Translated by Google

---

5.1-rasm: 1 radian faza siljishi bilan farq qiluvchi ikkita sinus to'lqin; ularning korrelyatsiya
koeffitsienti 0,54 ga teng.
Men faqat ikkita o'zgaruvchiga misol keltiraman. Birinchidan, men sinus to'lqinlarni turli xil faza
ofsetlari bilan tuzadigan funktsiyani aniqlayman: def make_sine(ofset): signal = thinkdsp.SinSignal
(freq=440, ofset=offset) wave = signal.make_wave(davomiylik=0,5, kadr tezligi=10000)
Men "ehtimol zaif" deyman, chunki korrelyatsiya koeffitsienti tomonidan ushlanmagan chiziqli
bo'lmagan munosabatlar ham bo'lishi mumkin. Quloq bo'lmagan aloqalar ko'pincha statistikada
muhim, ammo signalni qayta ishlash uchun kamroq ahamiyatga ega, shuning uchun men bu
erda ular haqida ko'proq gapirmayman.
52
Python korrelyatsiyalarni hisoblashning bir necha usullarini taqdim etadi. np.corrcoef har
qanday miqdordagi o'zgaruvchilarni oladi va har bir o'zgaruvchi juftligi orasidagi korrelyatsiyalarni
o'z ichiga olgan korrelyatsiya matritsasi hisoblaydi.
r ning kattaligi korrelyatsiyaning kuchini ko'rsatadi. Agar r 1 yoki -1 bo'lsa, o'zgaruvchilar
mukammal korrelyatsiya qilinadi, ya'ni agar siz birini bilsangiz, ikkinchisi haqida mukammal
bashorat qilishingiz mumkin. Agar r nolga yaqin bo'lsa, korrelyatsiya zaif bo'lishi mumkin,
shuning uchun agar siz bittasini bilsangiz, u sizga boshqalar haqida ko'p ma'lumot bermaydi,
to'lqin1 = make_sine (offset = 0) to'lqin2
= make_sine (offset = 1)
Keyin men turli xil ofsetlarga ega ikkita to'lqinni yarataman:
5-bob. Avtokorrelyatsiya
o'zgaruvchan yuqori, ikkinchisi esa yuqori bo'lishga intiladi. Agar r manfiy bo'lsa, korrelyatsiya
salbiy bo'ladi, shuning uchun bir o'zgaruvchi yuqori bo'lsa, ikkinchisi past bo'ladi.
qaytish to'lqini
0,5
0.000
Vaqt (lar)
1.0
0,002
0,008
0,5
0,004
0,0
0,010
1.0
0,006
to'lqin1
to'lqin2
Machine Translated by Google

---

5.2-rasm: Ikki sinus to'lqinning fazaga bog'liqligi
korrelyatsiya 1 ga teng.
N - 1 standartidan foydalanish o'rniga yuqoridagi tenglama .
5.2-rasmda a uchun korrelyatsiya va faza siljishi o'rtasidagi bog'liqlik ko'rsatilgan
ddof=0 varianti korrcoefni N ga bo'lish kerakligini ko'rsatadi, xuddi shunday
ofset 360 daraja farq qiladi. O'sha paytda biz to'liq doira va keldik
[0,54 1. ]]
180 daraja fazadan tashqarida, bu -1 korrelyatsiyasini beradi. Keyin u qadar ko'payadi
Faza ofsetining ortishi bilan bu korrelyatsiya to'lqinlar paydo bo'lguncha kamayadi
[[ 1. 0,54]
>>> corr_matrix = np.corrcoef(wave1.ys, wave2.ys, ddof=0)
korrelyatsiya o'rtacha.
Diagonaldan tashqari elementlar bizni qiziqtirgan qiymatni, to'lqin1 va to'lqin2 ning
korrelyatsiyasini o'z ichiga oladi. 0,54 qiymati ning kuchini ko'rsatadi
bir to'lqin baland, ikkinchisi odatda yuqori, shuning uchun biz ularning korrelyatsiya qilinishini
kutamiz.