sqrt_spec = sqrt(spec / (nfft))  # sqrt(spec(wn)*dw )