Con questo programma verrà dimensionato un banco di filtri come combinazione di un filtro polifase e FFT.

Il procedimento seguito consiste nello scomporre lo spettro di un generico segnale in un numero di canali da noi fissato a priori, utilizzando un opportuno filtro passabasso e sottocampionando sia la risposta all’impulso del segnale che del filtro per il numero di canali considerati.

La risposta all’impulso del filtro è introdotta utilizzando una workspace ossia un file di dati di MATLAB  ( hanno estensione ‘mat’ ).

 

banco.m

 

% PROGETTO DI UN BANCO DI FILTRI ( decomposizione multicanale del segnale )

 

% La risposta all'impulso del filtro passabasso verrà introdotta dall'esterno

 

% Lo spettro del segnale introdotto quindi sarà suddiviso in un numero di blocchi pari

% a quelli richiesti dalle nostre specifiche per cui non ci sarà alcuna elaborazione

% del segnale in queste bande

 

clear all

 

err=0;

while err==0 

   disp('inserisci il nome del file (.mat) contenente i campioni ');

   H=input('della risposta all''impulso del filtro passabasso : ','s');

   H=[H,'.mat'];

   if exist(H)==0

      disp('LA WORKSPACE INTRODOTTA NON ESISTE');

   else

      err=1;

   end

end

load(H);

N=length(h);

strN=int2str(N);

disp(['la risposta all''impulso del filtro è composta da ',strN,' campioni']);

err=0;

while err==0

   M=input('introduci il numero di canali : ');

   if M>N

      disp('Il numero di canali richiesto è troppo elevato! ');

   else

      err=1;

   end

end

if floor(N/M)~=ceil(N/M)

   h=[h',zeros(1,M-mod(N,M))];

   disp('ATTENZIONE : il rapporto tra la lunghezza della risposta all''impulso ');

   disp('             del filtro e il numero di canali non è un numero intero! ');

   N=length(h);

end

K=N/M;

err=0;

while err==0

   F=input('inserisci il nome della function contenente il segnale : ','s');

   if exist(F)==0

      disp('LA FUNZIONE INTRODOTTA NON ESISTE');

   else

      err=1;

   end

end

f=feval(F);

L=length(f);

[H,w]=freqz(h,1,L);

 

figure

plot(w/(2*pi),abs(H))

xlabel('f/fc')

set(gcf,'color',[1,1,1]);

title('spettro del filtro utilizzato')

axis([min(w/(2*pi)) max(w/(2*pi)) 0 max(abs(H))])

for m=1:M

   for k=0:K-1

      hh(k+1)=h(k*M+m);

   end

   pass1=[zeros(1,m-1) f];

   for k=0:floor((L-1)/M)

      ff(k+1)=pass1(k*M+1);

   end

   pass2=conv(ff,hh);

   lung=length(ff);

   convol(:,m)=pass2(1:lung)';

   clear ff hh

end

for i=1:lung

   dft(i,:)=fft(convol(i,:));

end

for i=1:M

   spettro(:,i)=abs(fft(dft(:,i)))*M;

end

 

mas=max(max(spettro));

 

freq=(-1/(2*M):1/(M*lung):(1/(2*M))*(1-1/(M*lung)))+1/M;

for i=1:M

   figure

   freq=freq-1/M;

   plot(freq,[spettro(ceil(lung/2)+1:lung,i)' spettro(1:floor(length(spettro)/2)+1,i)'])

   xlabel('f/fc')

   set(gcf,'color',[1,1,1]);

   stri=int2str(i);

   title(['spettro del segnale nella banda ',stri])

   axis([min(freq) max(freq) 0 mas])

   fine=[spettro(ceil(lung/2)+1:lung,i)' spettro(1:floor(length(spettro)/2)+1,i)' fine];

end

figure

freq=((1/(2*M))-1:1/(M*lung):(1/(2*M))*(1-1/(M*lung)));

plot(freq,fine)

xlabel('f/fc')

set(gcf,'color',[1,1,1]);

title('spettro del segnale rimettendo insieme tutte le bande')