This article is devoted to simulation of classical phase-locked loop (PLL). Based on new analytical method for computation of phase detector characteristics (PD), an realization in Simulink for simulation of classical PLL in phase space for general types of signal waveforms is done. This enables to avoid a number of numerical problems in the simulation of PLL in signal.