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