Accurate upper-lower bounds on homogenized matrix by FFT-based Galerkin method