This paper presents the implementation of a novel nonlinear PMSM model for signal injection-based sensorless polarity detection algorithm development. The model extension incorporates a quadratic flux-current function to describe the polarity-dependent saliency of the machine. The model was implemented in the Mathworks MATLAB-Simulink environment, which is widely used for numerical simulations of electric machines and drives. The presented simulation results include non-modulated square-wave injection and modulated sinusoidal injection which are commonly used in sensorless algorithms. Furthermore, these simulation results were compared to corresponding measurements, and the comparison between the simulated and measured data is also provided. The models are open source and available on GitHub.