function f2=rekonstruktionper(f1,h,k0,g,k1,J); 

% periodische Wavelet-Transformation
% f1 enthŠlt am Schluss den Skalierungsanteil und die Waveletanteile in der Reihenfolge
% (a^J,c^J,c^(J-1),...,c^1)
% h Tiefpassfilter
% k0 Index des Anfangskoeffizienten von h
% g Hochpassfilter
% k1 Index des Anfangskoeffizienten von g
% J Anzahl der Level

[L1,L2]=size(h);
[M1,M2]= size(g);
[N1,N2] = size(f1);



a=zeros(1,N2);  % Hilfsvektor
a=f1;
a1=zeros(1,N2); % Hilfsvektor
b=zeros(1,N2); % Hilfsvektor
for j=J:(-1):1
   Nneu=2*N2/2^j;
   b(1,1:Nneu/2)=f1(1,Nneu/2+1:Nneu);
   for n=0:Nneu-1
     r1=ceil((-k0-L2+1+n)/2);
	 r2=floor((n-k0)/2);
     for r=r1:r2 a1(1,n+1)=a1(1,n+1)+h(n-2*r-k0+1)*a(1,mod(r,Nneu/2)+1);
	 end;
	 r1=ceil((-k1-M2+1+n)/2);
	 r2=floor((n-k1)/2);
     for r=r1:r2 a1(1,n+1)=a1(1,n+1)+g(n-2*r-k1+1)*b(1,mod(r,Nneu/2)+1);
	 end;
   end;
   for l=1:Nneu a(l)=a1(l); end;
   b=zeros(1,N2);
   a1=zeros(1,N2);
end;
f2=a;


%f2 enthŠlt den rekonstruierten Vektor der LŠnge N2



 
