When solving hyperbolic Boltzmann model equation with discrete velocity models (DVM), the strong discontinuity of the velocity distribution function can be captured well by utilizing the non-oscillatory and non-free parameter dissipation (NND) finite difference scheme. However, most NND scheme solvers march in time explicitly, which compromise the computation efficiency due to the limitation of stability condition, especially when solving unsteady problems. In order to improve the efficiency, an implicit scheme based on NND is presented in this paper. Linearization factors are introduced to construct the implicit scheme and to reduce the stencil size. With the help of dual time-stepping method, the convergence rate of unsteady rarefied flow simulation can be massively improved. Numerical tests of steady and unsteady supersonic flow around cylinders are computed in different flow regimes. Results are shown to prove the validity and efficiency of the  implicit scheme.