The paper develops a phase-coordinates model for analyzing unbalanced electrical networks, enabling the specification of single-phase or three-phase harmonic injections of any sequence. The model is applied to a four-wire delta system, demonstrating its ability to represent structural asymmetries and identify resonances at different frequencies.