The numerical mode matching method (NMM) with the perfectly matching layer (PML) has been presented in this paper to solve electromagnetic fields in stratified media with some possible lossless layers. When the 2D inhomogeneous media is lossy media, the field in this inhomogeneous media can be solved efficiently by using the NMM with the concepts of the general reflection matrix and general transmission matrix. Unfortunately, this method failed in the case of lossless media. In order to make the NMM available in the media with lossless layers, the PML has been introduced into the NMM in this paper to build the truncated absorbing boundary in lossless layers. Some numerical examples have also been given in this paper.