% DTFT priklad
% ukazku DTFT spravime naraz pre 4 diskretne signaly v case - vysledkom su styri "spojite" funkcie vo frekvencii
% 4 konecne diskretne signaly na vstupe ziskame pomocou funkcie "wfilters" ...
% (ako sa da uz z nazvu funkcie vytusit, jedna sa o impulzove charakteristiky "waveletovych" filtrov)
% dalsim postupom ziskame ich DTFT
% vidime, ze sme dostali magitudove frekvencne charakteristiky (vid dodatok skripta o Z transformacii)

SIGSIZE=512

[LO_D,HI_D,LO_R,HI_R] = wfilters('Bior1.3')
   n=0:1:SIGSIZE-1;
   F1=zeros(1,SIGSIZE);
   for a=1:length(LO_D)
      F1=F1+LO_D(a)*exp(-i*2*3.1415926*(a-1)*n/SIGSIZE);
   end
   plot(n,abs(F1));
   hold on
 
   n=0:1:SIGSIZE-1;
   F2=zeros(1,SIGSIZE);
   for a=1:length(LO_R)
      F2=F2+LO_R(a)*exp(-i*2*3.1415926*(a-1)*n/SIGSIZE);
   end
   plot(n,abs(F2), 'red');
   hold on
   
   n=0:1:SIGSIZE-1;
   F3=zeros(1,SIGSIZE);
   for a=1:length(HI_D)
      F3=F3+HI_D(a)*exp(-i*2*3.1415926*(a-1)*n/SIGSIZE);
   end
   plot(n,abs(F3),'green');
   hold on
   
   n=0:1:SIGSIZE-1;
   F4=zeros(1,SIGSIZE);
   for a=1:length(HI_R)
      F4=F4+HI_R(a)*exp(-i*2*3.1415926*(a-1)*n/SIGSIZE);
   end
   plot(n,abs(F4),'yellow');
   hold on
      
%   Ukazeme niektore zakladne vlastnosti waveletovych filtrov
%   F_LOW=abs(F1.*F2);
%   F_HI =abs(F3.*F4);
%   plot(n,F_LOW,'red');
%   plot(n,F_HI,'green');
%   plot(n,F_LOW+F_HI,'yellow');   
%   hold on
 

 return