Abstract-A novel method based on the hybrid volume-surface integral equation (VSIE) and the impedance matrix interpolation technique is presented for the fast analysis of microstrip antennas in frequency sweeps. A novel impedance matrix interpolation scheme is extended to the impedance matrix associated with VSIE, thus providing high accuracy, high efficiency, and large interpolation bandwidth for metal-dielectric composite problems. To demonstrate the effectiveness and accuracy of the proposed technique, numerical results for typical rectangular patch antennas and a broadband U-slot rectangular patch antenna are presented. Good agreement among the interpolation results, the exact method of moments (MoM) solutions, the finite element method (FEM) solutions, and measured data is observed over the bandwidth. The interpolation bandwidth is further investigated through a scattering problem. Numerical results show that high accuracy is obtainable within 10 : 1 bandwidth.