Filtru IIR de tip Butterworth trece banda
Cerinte
Parte teoretica
Filtre Butterworth (au caracteristica plata in banda de trecere, banda de tranzitie nu foarte ingusta, iar atenuarea in banda de oprire este mare).
Functia folosita este butter( ),
[b, a] = butter(ord, wn, optiuni), unde:
ord reprezinta ordinul filtrului, avand valori uzuale in gama 10.15,
wn reprezinta frecventa de taiere normalizata, f / (fe/2), cu valori intre 0 si 1; pentru filtrele de tip trece jos si trece sus, wn este un scalar, iar pentru cele de tip trece banda si opreste banda wn este un vector linie cu 2 elemente, reprezentand banda de trecere, respectiv cea de taiere; toate frecventele de taiere sunt date pentru o atenuare de 1/ 3dB, daca caracteristica este considerata in coordonate semilogaritmice), in caracteristica magnitudine-frecventa.
optiuni:
-daca acest camp nu contine nici un caracter, filtrul este trece jos sau trece banda, functie de forma argumentului wn;
-daca contine caracterele 'high', filtrul este trece sus,
-daca contine caracterele 'stop', filtrul este opreste banda.
Se precizeaza ca pentru filtrele trece banda si opreste banda, filtrul obtinut este de ordinul 2ord, adica vectorii b si a au 2ord +1 elemente.
Pentru obtinerea ordinului minim necesar ca un filtru sa aibe caracteristici impuse referitoare la banda (benzile) de tranzitie se foloseste functia buttord(), astfel:
[ord, wn] = buttord(wp, ws, Rp, Rs), unde
Rp reprezinta atenuarea, in banda de trecere, la wp, in dB,
Rs reprezinta atenuarea, in banda de oprire, la ws, in dB,
iar ord si wn, rezultatele furnizate de aceasta functie, au aceeasi semnificatie ca si la functia butter.
Exemple:
1. Pentru un filtru trece jos, este necesar ca wp<ws.
2. Pentru un filtru trece sus, este necesar ca wp>ws.
3. Pentru un filtru trece banda este necesar ca elementele vectorului wp sa fie intercalate intre cele ale vectorului ws, de exemplu wp si ws=[0.4 0.8], ca in fig.5. Atenuarile pot fi de exemplu, Rp=1 si Rs=20 (valorile sunt considerate in dB, avand de fapt, cum se vede si in figura, semnul minus).
Pentru un filtru opreste banda este necesar ca elementele vectorului ws sa fie intercalate intre cele ale vectorului wp, de exemplu wp si ws Apoi, pentru toate cazurile, avand parametrii ord si wn, se determina coeficientii filtrului cu functia % definitia vectorilor de frecventa M=[-30 -30 -2 -2 -30 -30]; definitia vectorilor de magnitudine subplot(212); plot(F*fe/2,M);%trasarea caracteristicii magnitudine-frecventa figure(2); subplot(211) plot(t,x);%semnal suma de intrare subplot(212); plot(t,y);%semnalul de iesire Nota: Wp -frecventa limita normalizata pentru banda de trecere; Ws -frrecventa limita normalizata pentru banda de oprire In urma rularii programului s-au obtinut urmatoarele rezultate:
se poate observa ca frecventele de trece, respective de oprire, a benzii se incadreaza in cele specificate in cerintele proiectului (500,700,1200,1400); totodata se poate observa faptul ca, caracteristica magnitudine-frecventa este asemanatoare cu cea impusa. in figura de mai sus putem observa semnalul de la intrare format din suma a doua semnale cu frecvente diferite, unul avand frecventa de 700Hz iar celalat de 2000Hz in semnalul de iesire filtrat este prezent doar semnalul care are frecventa de 700Hz. am considerat doua semnale care au fecventele de 1800Hz, respective, 2000Hz; dupa trecerea semnalului respectiv prin filtru ,semalul de iesire, dupa filtrare, trebuie sa fie 0(zero). dupa simulare acesta este, intradevar,0, adica la iesirea filtrului nu avem nici un semnal. Document InfoAccesari: 5699 Apreciat: Comenteaza documentul:Nu esti inregistratTrebuie sa fii utilizator inregistrat pentru a putea comenta Creaza cont nou A fost util?Daca documentul a fost util si crezi ca meritasa adaugi un link catre el la tine in site in pagina web a site-ului tau.
Copyright © Contact (SCRIGROUP Int. 2024 ) |