Based on the well-known analogy of propagation modes between metal-insulator-metal (MIM) waveguides and microwave transmission lines, an improved transmission line model (ITLM) have been proposed to describe reflection coefficients of a MIM waveguide stub. The Fresnel formula have been employed to predict the phase shift resulting from the reflection of surface plasmon polaritons (SPP) from the stub's end, which allows the ITLM to describe the influence of the stub-shaped resonators arranged in the side-coupled configuration on the transmission spectrum of main waveguide more precisely. In this paper, the improved transmission line models of MIM waveguide side coupled to single stub and double-side symmetric stubs have been proposed, and the predicted transmission spectrum by the ITLM is well agreed with simulation results of the finite element method. Specifically, transmission minimum predicted by ITLM is more close to that calculated by numerical methods than the classical transmission line model. Moreover, it is verified that ITLM possess much higher efficiency in calculating the transmission spectrum of the MIM waveguide connected with stub-shaped resonators than the finite element method employed in this work.