Numerical analysis of Bragg fibers using a compact 1D finite-difference frequency-domain method