Modern production processes frequently require steady-state analysis of continuous dynamic systems. Traditional numerical approaches, however, fall short in efficiency when tasked with addressing large-scale or dynamic problems. To tackle the inverse problem inherent in stability analysis, this study presents an innovative approach by integrating a combined excitation function into the foundational zeroing neural network (ZNN) model. This integration constrains the ZNN model, evolving it into an enhanced EZNN model specifically designed for solving the inverse of dynamic complex matrices. Additionally, this paper conducts a rigorous theoretical analysis of the robust performance of the EZNN model when excited by the combined function, both in the presence and absence of noise interference. The model solution process is promoted by using a class of high-dimensional continuous dynamic systems as an example, and numerical simulation experiments are used for validation. Considering the dynamic system satisfying

{ A(t)=(4+sin(2t) 4-cos(2t) 5+sin(2t))C(t)=(cos(t)sin(t)-cos(t)-sin(t)cos(t)sin(t))b(t)=4+cos(4t),d(t)=(cos(2t),cos(2t))T
\left\{ \matrix{ A(t) = \left( {4 + \sin (2t)\quad 4 - \cos (2t)\quad 5 + \sin (2t)} \right) \hfill \cr C(t) = \left( {\matrix{{\cos (t)} & {\sin (t)} & { - \cos (t)} \cr { - \sin (t)} & {\cos (t)} & {\sin (t)} \cr } } \right) \hfill \cr b(t) = 4 + \cos (4t),d(t) = {(\cos (2t),\cos (2t))^T} \hfill \cr} \right.

, the error E1(x(t),t) obtained by the EZNN model with combinatorial function excitation always remains negative or tends rapidly to 0. The x (t) obtained by the model converges rapidly to an exact solution of the system. Through the discussion of parametric conditions, it is also found that increasing the value of parameter γ increases the rate of convergence of the ZNN model.