In this paper, we introduce a fourth-order iterative scheme tailored for a parametric matrix family to calculate the inverse square root of matrices. Initially, we will design a fourth-order algorithm for the matrix equation, followed by the development of a coupled scheme to address stability concerns. Moreover, we will conduct a comprehensive analysis of convergence and stability, substantiated by numerical examples.