In this paper, a numerical method for analyzing coupling between high-altitude electromagnetic pulse (HEMP) as external field and a shielded twisted pair (STP) cable is proposed, which is based on an expanded chain matrix. Load responses of electromagnetic (EM) field excitation in uniform transmission line (TL) are solved by Baum-Liu-Tesche (BLT) equations in frequency domain, however, it is difficult to apply BLT equations to solve load responses of STP cable because the iteratively changing configuration of each twisted pairs are involved in cable. To avoid this problem and decrease memory and CPU time, we proposed the expanded chain matrix modeling method that is calculated using ABCD parameters, and applied multi-conductor transmission line (MTL) theory to consider the EMP coupling effectiveness of each twisted pairs. The results implemented by the proposed method are presented and compared with those obtained by the finite-difference time domain (FDTD) method as a kind of 3D full wave analysis.