analysis and optimization of these quantum processors. Here we introduce a new method of modelling that can be applied to superconducting structures involving multiple Josephson junctions, high-Q superconducting cavities, external ports, and voltage sources. Our technique, an extension of our previous work on single-port structures [1], permits the derivation of system Hamiltonians that are capable of representing every feature of the physical system over a wide frequency band and the computation of T1 times for qubits. We begin with a black box model of the linear and passive part of the system. Its response is given by its multiport impedance function Zsim(w), which can be obtained using a finite-element electormagnetics simulator. The ports of this black box are defined by the terminal pairs of Josephson junctions, voltage sources, and 50 Ohm connectors to high-frequency lines. We fit Zsim(w) to a positive-real (PR) multiport impedance matrix Z(s), a function of the complex Laplace variable s. We then use state-space techniques to synthesize a finite electric circuit admitting exactly the same impedance Z(s) across its ports; the PR property ensures the existence of this finite physical circuit. We compare the performance of state-space algorithms to classical frequency domain methods, justifying their superiority in numerical stability. The Hamiltonian of the multiport model circuit is obtained by using existing lumped element circuit quantization formalisms [2, 3]. Due to the presence of ideal transformers in the model circuit, these quantization methods must be extended, requiring the introduction of an extension of the Kirchhoff voltage and current laws.