typ='db2';
[h_,g_,h,g]=wfilters(typ);
Hom=fft(h);
for i=1:poceth/2
    suma22(i)=(abs(Hom(i)))^2+(abs(Hom(i+poceth/2))^2);
end
suma22'
