blogy logo
login PRIHLÁS SA
BLOG deadawp
ČLÁNKY
DISKUSIE
2
SLEDOVAŤ BLOG
Programátor
deadawp



Frekvenčné spektrum - MATLAB
pridal deadawp 10.1. 2020 o 2:15 (naposledy upravené 2.3. 2022 o 10:24)



Frekvenčné spektrum môže reprezentovať jednotlivé spektrá signálu a analyzovať rôzne stavy, napríklad či dochádza k aliasingu, t.j. prekrývaniu spektier signálu. Využitie tejto metódy nájde uplatnenie pri návrhu číslicových filtrov FIR, IIR pri skúmaní ich vlastností. V tejto implementácii použijeme interpolátor (respektíve operáciu interpolácie) na pôvodné spektrum signálu. Interpolácia dokáže stlačiť pôvodné spektrum signálu. Na voľné miesto je možné vložiť kópiu spektra (nazývame to aj nultou hodnotou). Interpolátor teda predchádza aliasingu, nakoľko v ňom nedochádza k prekrývaniu spektier. V prípade vyšších interpolácii môžeme spektrum stlačiť ešte viac a vložiť vždy  L-1 (L mínus 1) kópii do spektra namiesto pôvodnej vzorky. L označuje stupeň interpolácie.

Script v prostredí MATLAB pracuje s menovitou frekvenciou signálu 1Hz. Nakoľko budeme pracovať s uhlovou (kruhovou) frekvenciou, potrebujeme si frekvenciu vyjadriť  v rad/s (možno vyjadriť aj ako s-1). Následne si určíme vzorkovaciu frekvenciu s hodnotou 2*uhlová frekvencia. Následne násobky vzorkovacej frekvencie do + i do - *2, *3. Tieto hodnoty budú následne slúžiť na reprezentáciu na x osi grafu. MATLAB umožňuje upraviť funkčné hodnoty osi x a y použitím xticks funkcie a následne funkciou xticklabels priradiť aj textové hodnoty - labely k jednotlivým hodnotám. Následne môžeme do čiarového grafu vykresliť spektrum pôvodného signálu. 

Vidíme, že v spektre nie je aliasing - prekrývanie spektier signálu. Následne vykonáme operáciu s interpolátorom, ktorý umožňuje stlačiť pôvodné spektrum o faktor L-1. Interpolátor má faktor L, ktorý sa na spektrum aplikuje. Všeobecne platí, že interpolátor do spektra vloží L-1 kópii, respektíve nultú zložku. Teda aby interpolátor zmenil spektrum, musí mať hodnotu minimálne, pri L=1 nedôjde ku zmene - stlačeniu pôvodného spektra signálu. V tejto implementácii sme použili faktor L=2 a L=3. Výsledkom sú spektrá, do ktorých sú vložené kópie a pôvodné spektrum je stlačené. Vzorkovacia frekvencia je nadvzorkovaná - vynásobená * faktor L.

Interpolátor má využitie predovšetkým v elektotechnike pri spracovaní a prenose signálov, napríklad pri odosielaní nadvzorkuje signál, čím predchádza prípadnej chybe, nakoľko rovnaký signál je nadvzorkovaný, t.j. môže ho reprodukovať vernejšie. 

Výstup celého grafu s priebehom hlavného spektra a interpolácii L=2, L=3:

Programová impementácia dostupná na Githube: 
https://github.com/martinius96/MATLAB-scripty/blob/main/Cislicove_spracovanie_signalov/frekv_spektrum_interpolacia.m



Prístupov 14292
Kvalita článku
hlasov 0

PRÍSPEVKY
SLEDUJETE
Prosím prihláste sa pre možnosť pridania komentáru.
Prihláste sa, alebo použite facebook login facebook login
ĎALŠIE ČLÁNKY V BLOGU
Update grafov, ukážkových kódov - Hladin...
[ 16.11.2024] (príspevkov 0)
XIAO ESP32-C6 od Seeed Studio
[ 10.11.2024] (príspevkov 0)
ESP-IDF v4.4.2 prechod na v5.2
[ 5.11.2024] (príspevkov 0)
Krabička pre RFID DOMINATOR 2.0
[ 18.10.2024] (príspevkov 0)
Známe neduhy - Ecotec 1.6E Opel/GM Chevr...
[ 8.10.2024] (príspevkov 0)
Cold-start BQ25570 vs BQ25504 od Texas I...
[ 2.8.2024] (príspevkov 0)
EG21-G - MQTT pripojenie na Thingsboard
[ 5.6.2024] (príspevkov 0)
Quectel EG21-G - HTTP request
[ 17.5.2024] (príspevkov 0)
Tip na darček k jubileu 60
[ 29.4.2024] (príspevkov 0)