function tu=xbf(B,l,n)
%把B矩陣做小波反變換，l為B的長度，n為小波變換時所用的層數，其值只有1，2，3
if n<1&n>3
   disp('n為小波變換時所用的層數，其值只有1，2，3');
else
   ca1=B(1:l/2,1:l/2);
   ch1=B(1:l/2,l/2+1:l);
   cv1=B(l/2+1:l,1:l/2);
   cd1=B(l/2+1:l,l/2+1:l);
   if n>1
      l=l/2;
      ca2=ca1(1:l/2,1:l/2);
      ch2=ca1(1:l/2,l/2+1:l);
      cv2=ca1(l/2+1:l,1:l/2);
      cd2=ca1(l/2+1:l,l/2+1:l);
   end
   if n==3
      l=l/2;
      ca3=ca2(1:l/2,1:l/2);
      ch3=ca2(1:l/2,l/2+1:l);
      cv3=ca2(l/2+1:l,1:l/2);
      cd3=ca2(l/2+1:l,l/2+1:l);
   end
   if n==3
      ca2=idwt2(ca3,ch3,cv3,cd3,'db1');
   end
   if n>1
      ca1=idwt2(ca2,ch2,cv2,cd2,'db1');
   end
   tu=idwt2(ca1,ch1,cv1,cd1,'db1');
end