A computer algorithm for the harmonic analysis on the steady state of singl
e phase power systems containing a transformer saturation model is presente
d in this paper. The model proposed has been developed on the finite differ
ence numerical method. The algorithm illustrates the application of the mod
el. The numerical results have been validated with PSpice simulation progra
m. The method proposed exploits the double-diagonal structure of the matrix
and a fast resolution algorithm is obtained. Different forms of representi
ng transformer saturation have been analysed for improvement of the converg
ence. (C) 1999 Elsevier Science Ltd. All rights reserved.