% DTFT waveletovych filtorv a podmienka uplnej rekonstrukcie

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));
   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));
   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));
   hold on
      
   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