This paper is devoted to the study of simulating a large class of self-similar processes. Since most current simulation approaches are limited to case-by-case studies, every existing approach has its constraints and flaws; hence a general and efficient simulation approach is in demand. Our study sheds some light in this direction. The paper's contributions are bi-fold. First, reviews and improvements are made to some existing methods for simulating specific self-similar processes. Second, we propose a novel method to simulate a general self-similar process, where we use a modified inverse Lamperti transformation to transform self-similarity to stationarity. Successful applications are made to simulate fractional Brownian motion and sub-fractional Brownian motion.