In power systems there are complex transformer structures, whose accurate analysis is not possible using the techniques available today. This paper presents a systematic data driven analysis method for coupled inductors of arbitrary complexity. The method first establishes a winding matrix N mapping the windings to the limbs of the transformer. A permeance matrix P is created from the reluctance network of the magnetic core. A generalized inductance matrix L mapping currents in the transformer windings to the induced voltages is generated based on the winding (N) and permeance (P) matrices. The inductance matrix representation of a coupled inductor is then transformed to an admittance matrix, which can be integrated to the nodal analysis of the electrical circuit surrounding the coupled inductor. The method presented is validated by simulations with real transformer structures using electromagnetic transient program (EMTP/ATP).