Compresión de imágenes usando la transformada de wavelet y el algoritmo de Huffman

La implementación en matlab estupenda, pero me da un fallo, si alguién pudiera ayudarme....

??? Undefined function or method 'reconstruir' for input arguments of type 'double'.

Error in ==> waveletCarmen at 68
recuperada=reconstruir(fpbtpbd,fpbtpad,fpatpbd,fpatpad);

Ya que la función recontruir no esta definida, podría alguien darme la función definida???.
O si el fallo esta en otra parte, decidmelo por favor.
Muchas gracias!!!!!

la función reconstruir es la siguiente

function img = reconstruir(pbpb,pbpa,papb,papa)
pbpb=deshacacerDiezmado(pbpb,2);
pbpa=deshacacerDiezmado(pbpa,2);
papb=deshacacerDiezmado(papb,2);
papa=deshacacerDiezmado(papa,2);

fpb=[1/sqrt(2),1/sqrt(2)];
fpa=[-1/sqrt(2),1/sqrt(2)];

pb1=filter2(fpb',pbpb);
pb2=filter2(fpa',pbpa);

pa1=filter2(fpb',papb);
pa2=filter2(fpa',papa);

pb=pb1+pb2;
pa=pa1+pa2;

pb=deshacacerDiezmado(pb,1);
pa=deshacacerDiezmado(pa,1);

pb=filter2(fpb,pb);
pa=filter2(fpa,pa);

img=pb+pa;

necesitaran también esta funcion llamada deshacerDiezmado

function img=deshacacerDiezmado(imagen,sentido)%sentido, 1 para hacerlo verticalmente y 2 para hacerlo horizontal
[h,w]=size(imagen);
if sentido == 1
img=zeros(h,w*2);
else
img=zeros(h*2,w);
end

if sentido == 1
for i=1:2:w*2
img(:,i)=imagen(:,(i+1)/2);
end
else
for i=1:2:h*2
img(i,:)=imagen((i+1)/2,:);
end
end
end

si les hace falta algo dejen sus comentarios.

Hola podrias, subir ne nuevo el archivo. Es que no lo he logrado bajar.