This paper aims to study the model of Klein–Gordon–Schrödinger equations with a memory term and locally distributed damping : \( {\left\{ \begin{array}{ll} i\psi '+ \Delta \psi + i\alpha (|\psi |^{2} +1)\psi =- \phi \psi ~\text{ in }~\Omega \times (0, \infty ), \\ \phi '' - \Delta \phi +\displaystyle \int _{0}^{t} g(t-\tau ) div[a(x) \nabla \phi (\tau )] d \tau +b(x)\phi '= |\psi |^{2\theta }\chi _{\omega } ~\text{ in }~\Omega \times (0, \infty ),\\ \end{array}\right. } \) where \(\Omega \) is a bounded domain of \(\mathbb {R}^n\) with \(n\le 3\) and smooth boundary \(\partial \Omega =\Gamma \) . Here, \(\alpha \) and \(\theta \) are positive constants. In this work, \(\chi _{\omega }\) represents the cutoff function of \(\omega \) . Assuming that a and b are non-negative functions such that \(a(x) + b(x) \ge \delta > 0\) in \(\Omega \) , the exponential decay rate is demonstrated for each regular solution of the above system.