Applications of frequency-domain analysis in pipelines and pipe networks include resonance analysis, time-domain simulation and fault detection. Current frequency-domain analysis methods are restricted to series pipelines, single-branching pipelines and single-loop networks and are not suited to complex networks. This paper presents a number of formulations for the frequency-domain solution in pipe networks of arbitrary topology and size. The formulations focus on the topology of arbitrary networks and do not consider any complex network devices or boundary conditions, other than head and flow boundaries. The frequency-domain equations are presented for node elements and pipe elements, which correspond to the continuity of flow at a node and the unsteady flow in a pipe, respectively. Additionally, a pipe-node-pipe and reservoir-pipe pair set of equations are derived. A matrix-based approach is used to display the solution to entire networks in a systematic and powerful way. Three different formulations are derived based on the unknown variables of interest that are to be solved for, being the head-formulation, flow-formulation and the head-flow-formulation. These hold significant analogies to different steady-state network solutions. The frequencydomain models are tested against the method of characteristics (a commonly used timedomain model), with good result. The computational efficiency of each formulation is discussed with the most efficient formulation being the head-formulation.