Simulink models, MATLAB scripts and output plots for power electronics, electric drives and power systems: DC/induction/PMSM/BLDC drives, load flow, fault analysis, transient stability, load frequency control and relay characteristics.
Made by Harish M (@harishkernel), B.E. Electrical and Electronics Engineering, Easwari Engineering College, Chennai.
- MATLAB with Simulink
- Simscape Electrical (Specialized Power Systems), needed for the thyristor, IGBT, DC machine, induction machine and PMSM blocks
- Clone or download the repo.
- In MATLAB, set the current folder to the topic folder you want.
- Open the
.slx(or.mdl) file and press Run. - Double-click the scopes to see speed, torque, current and voltage. Each folder's
output.pdfshows the expected result.
Two Power Systems scripts call functions that are not part of the scripts themselves. They are included in Power_Systems/Helper_Functions/ and come from the MATLAB files that accompany Hadi Saadat's Power System Analysis. All credit for these functions goes to the author.
| Script | Helper functions it calls |
|---|---|
Power_Systems/Fault_analysis_LL_LG_LLG/Fault_analysis.m |
zbuild, symfault, lgfault, llfault, dlgfault |
Power_Systems/Newton_Raphson_Load_Flow_Analysis/Load_flow_analysis.m |
lfybus, lfnewton, busout, lineflow |
The other Power Systems scripts run on their own.
Add the helper folder to the MATLAB path before running those two scripts:
% 1. Set MATLAB's current folder to the repo root (EE_Coursework)
% 2. Add the helper functions (and subfolders) to the path
addpath(genpath('Power_Systems/Helper_Functions'));
% 3. Move into the topic folder and run the script
cd Power_Systems/Fault_analysis_LL_LG_LLG
Fault_analysisaddpath lasts only for the current MATLAB session. To keep it permanently, use Home > Set Path > Add with Subfolders, select Helper_Functions, then Save. If MATLAB reports Undefined function 'zbuild' (or lfybus etc.), the helper folder is not on the path.
Several scripts also ask for inputs at the prompt (input(...)); enter values when asked.
| Topic | Folder |
|---|---|
| Single-phase fully controlled rectifier circuit | Open |
| Topic | Folder |
|---|---|
| Y-bus calculation | Open |
| Z-bus calculation (building algorithm) | Open |
| Gauss-Seidel load flow analysis | Open |
| Newton-Raphson load flow analysis | Open |
| Fault analysis: LG, LL and LLG faults | Open |
| Transient stability of SMIB system: equal area criterion and critical clearing angle | Open |
| Two-area load frequency control (Simulink) | Open |
| Tie-line power and frequency deviation (Simulink) | Open |
| Inverse definite minimum time (IDMT) overcurrent relay characteristics | Open |
| Transmission line parameters: inductance and capacitance with ACSR conductors | Open |
| CT and PT testing (ratio and percent error) | Open |
| Topic | Folder |
|---|---|
| Single-phase fully controlled bridge converter fed separately excited DC motor | Open |
| Three-phase fully controlled converter fed DC motor | Open |
| Open-loop chopper fed (MOSFET) separately excited DC motor drive | Open |
| Closed-loop PI speed control of chopper fed DC motor | Open |
| Topic | Folder |
|---|---|
| Open-loop VSI fed 3-phase induction motor drive (SPWM) | Open |
| Vector control (FOC with SVPWM) of 3-phase induction motor drive | Open |
| Topic | Folder |
|---|---|
| Speed control of wound-field synchronous motor drive | Open |
| Open-loop speed control of PMSM drive | Open |
| BLDC / PMSM drive with Hall-sensor 6-step commutation | Open |
- Circuit: four thyristors on a single-phase AC source feeding a separately excited DC machine.
- Control: a ramp (repeating sequence) is compared with a control level to set the firing angle.
- Output: armature voltage, armature current, speed and electromagnetic torque.
- Circuit: six-thyristor bridge on a three-phase AC supply driving a DC machine with a load torque.
- Control: gate pulses T1 to T6 in sequence.
- Output: speed build-up, armature current, torque and output voltage ripple.
- Circuit: MOSFET chopper with pulse generator and freewheeling diode feeding the armature of a separately excited DC motor.
- Output: chopped voltage, armature current ripple, torque, speed and constant field current.
- Circuit: GTO chopper with series RLC branch and freewheeling diode feeding a DC machine.
- Control: speed reference (120 rad/s) minus measured speed goes to a PI controller (proportional gain 1.6, integral gain 16), then a saturation block and relay that switch the GTO.
- Output: speed and torque response.
- Circuit: three-phase IGBT voltage source inverter feeding a 460 V induction motor, with RMS voltage meters on the lines.
- Control: sinusoidal PWM reference generation.
- Output: rotor speed, electromagnetic torque and stator currents during start-up.
- Circuit: 50 HP, 460 V induction motor on a three-phase inverter from a DC link.
- Control: field-oriented (vector) control generating alpha-beta voltage references for space-vector PWM, with selectable speed and torque references.
- Output: phase currents, rotor speed and electromagnetic torque.
- Circuit: HV battery, converter and wound-field synchronous machine, with a two-quadrant chopper for field excitation.
- Control: speed request converted to torque command with limits.
- Output: rotor speed against reference, phase currents, torque and rotor electrical angle.
- Circuit: PMSM fed by a three-phase IGBT inverter.
- Control: sine-triangle comparison with complementary gate signals, plus a load torque step.
- Output: PWM pole voltages, phase currents, rotor speed and torque.
- Circuit: universal bridge on a DC source driving a permanent magnet synchronous machine.
- Control: rotor position (theta) converted to Hall signals, which are decoded into six switching signals.
- Y-bus / Z-bus: builds the bus admittance matrix from line impedance and line charging data, and the bus impedance matrix with the building algorithm.
- Load flow: Gauss-Seidel and Newton-Raphson solutions on small test systems, with bus voltages and line flows.
- Fault analysis: line-to-ground, line-to-line and double line-to-ground faults using sequence networks.
- Transient stability: critical clearing angle of a single-machine infinite-bus system using the equal area criterion.
- Load frequency control: two-area system and tie-line power flow in Simulink, with frequency deviation plots.
- IDMT relay: operating time against fault current for a standard inverse overcurrent relay.
- Transmission lines: inductance and capacitance per km for three-phase lines with bundled ACSR conductors.
- Instrument transformers: CT and PT ratio error calculations from test data.
Found a bug or want to suggest an improvement? Open an issue. You're welcome to reuse the models for learning; see the LICENSE file for terms.