.m file to solve harmonic load flow for power systems (Matlab 2009).
Line Matrix is line impedance matrix.
load_power Matrix is P+jQ matrix of the loads.
source matrix is the complex value of the harmonic sources and they considered as a harmonic current source for each harmonic frequency.
the code contains 3 example including IEEE 34 bus test system.
Please note that this algorithm first solves load flow for fundamental frequency and then solves harmonic load flow using the results of previous step, individually for each harmonic order.
**I'm not the main writer of this code and also I'm not professional programmer, I've just modified this code for my own purpose and just share it here for your help.
Masoud Farhoodnea, Azah Mohamed, Hussain Shareef. Harmonic Source Localization in Power Systems: Single and Multiple Harmonic Sources. 2011, LAMBERT Academic Publishing (LAP). ISBN: 978-3-8443-2639-0