Progetto filtro Matlab

[richqueen19]

Ciao a tutti, dovrei realizzare un filtro passabanda con bande di transizione strette il più possibile con Matlab, ho pensato al metodo delle finestre, sapreste indicarmi quale finestra è più adatta per ottenere questa specifica??Grazie!

[dimaios]

La frase "strette il piu' possibile" non ha molto senso se non inquadrata in un discorso di specifiche del problema.

Fissa i seguenti parametri :

[1] Frequenza di campionamento
[2] Banda passante

A questo punto scegli la struttura del filtro digitale che intendi progettare ( FIR o IIR ).
Si dovrebbe anche fissare la struttura del filtro ( Forma diretta, lattice a variabili di stato .... ecc. ).

Puoi utilizzare il metodo delle finestre anche se magari sarebbe piu' immediato mettere in serie un filtro passa alto con un filtro passa basso.


Se la banda che vuoi ottenere e' molto stretta puoi fare quanto segue.
Una tecnica valida e' quella di sintetizzare un notch filter e poi trasformarlo in un passa banda tramite la trasformazione ( questo nel dominio della s-Trasformata ma analogamente per la z-Trasformata ).

A riguardo ti invito a leggere i seguenti documenti per scegliere la tecnica che preferischi :

http://djj.ee.ntu.edu.tw/Notch_Filter.pdf
http://euler.fd.cvut.cz/publikace/files/vlcek/ieeenoc.pdf
http://euler.fd.cvut.cz/publikace/files/vlcek/ISCASnew.pdf
http://ntur.lib.ntu.edu.tw/bitstream/246246/142473/1/63.pdf


Esistono molti metodi per la sintesi dei filtri digitali ma senza ulteriori informazioni/specifiche e' difficile consigliarti.

[richqueen19]

Grazie mille..comunque hai ragione, le specifiche sono:
Fs=62500Hz;
Fc1=350Hz
Fc2=20000Hz;
avevo pensato ai filtri notch però siccome sono alle prime armi mi risulta un po' complicato implementarli in matlab, tuttavia ho realizzato un programmino con la finestra di kaiser, però per avere una banda di transizione stretta mi serve un numero di coefficienti molto elevato..quindi non so se può essere sbagliato...

[dimaios]

Se non ti serve fase lineare hai pensato di utilizzare un IIR ?

Potresti per esempio utilizzare una serie di due filtri ellittici ( la transizione e' molto pronunciata ).

Guarda questo documento per il confronto tra i vari filtri ( sezione "Comparison with other linear filters" ).
http://en.wikipedia.org/wiki/Chebyshev_filter

Se pero' in banda passante vuoi una risposta piatta devi utilizzare il filtro di Butterworth oppure Chebyshev del secondo tipo.