Presents physics-based models of bipolar power semiconductor devices and their implementation in MATLAB and Simulink. The devices are subdivided into different regions, and the operation in each region, along with the interactions at the interfaces which are analysed using basic semiconductor physics equations that govern their behaviour.