harmonic balance (h b ) analysis is an established technique for the analysis of nonlinear circuits. in this thesis, we present a method to find the transient and steady-state behavior of oscillators based on a harmonic balance implementation which proves to be faster than traditional time domain simulations. it is derived from the warped multi-time partial differential equation (wampde) approach developed in recent years. this approach efficiently separates the frequency modulation (fm) and amplitude modulation (am) in oscillators. a review of multi-time partial differential equation (mpde) and wampde is presented along with the motivation and interest in warped multiple time axes. the first proposed method shows how to obtain initial boundary conditions for a wampde system consistently with arbitrary physical initia l conditions in the system of ordinary differential equations in transient analysis. the second proposed method improves steady-state analysis since it does not require a good initial guess of the oscillation frequency and exploits the frequency-domain latency of circuits by using a different number of harmonics in each state variable. the wampde approach is used to simultaneously bring the circuit state to the region of convergence of the hb analysis and determine the optimum number of harmonics required at each node in the circuit. in both transient and steady-state, an adaptive time step control technique is employed in one of the time axes and this considerably reduces the computational effort. simulation results of different applications are described to demonstrate the performance of the proposed method. finally, the proposed methods are validated with experimental results.