In [1]:
using PyPlot
using WaveletScattering
In [3]:
ws = WaveletScattering
spec = ws.Spec1D(log2_size = 9, n_filters_per_octave = 4)
@time bank = ws.Bank1D(spec);
In [4]:
N = 1 << spec.log2_size
halfN = N >> 1
ωs = (-halfN):(halfN-1)
for ψ in bank.ψs
semilogx(fftshift(abs2(ψ[ωs])))
end
semilogx(fftshift(abs2(bank.ϕ[ωs])))
Out[4]:
In [6]:
lp = ws.renormalize!(bank.ϕ, bank.ψs, bank.spec)
semilogx(lp)
Out[6]:
In [ ]: