Finite-difference-based mode-matching method for 3-D waveguide structures under semivectorial approximation