Cycle-to-cycle variability of switching parameters inherent to memristive devices introduces significant problems in the design of neuromorphic systems and non-volatile memory. This study investigates the dynamics of a second-order memristive system incorporating capacitive effects that model parasitic charge within individual memristors, addressing both the technical need for accurate analysis of complex regimes and the demand for exploratory environments. Simulations were performed using CUDAynamics, an interactive software suite developed by the authors, which utilizes parallel computing, primarily via NVIDIA Compute Unified Device Architecture (CUDA). It integrates multiple analysis tools for dynamical systems, including bifurcation diagrams, the largest Lyapunov exponent and periodicity mapping, and interactive navigation in multidimensional parameter spaces. The memristive system was discretized applying multiple integration methods with a fixed time step and various waveforms of the input signal. Analysis tools revealed well-defined regions of chaotic dynamics in the memristor resistance parameter space as functions of input signal properties. Sinusoidal and triangular waveforms produced topologically similar distributions of dynamical regimes, whereas the square waveform, mimicking digital inputs, generated distinct dynamical patterns while still preserving chaotic trajectories under specific conditions. Interactive visualization capabilities of CUDAynamics effectively demonstrate attractor evolution and hysteresis deformation, providing immediate visual feedback that significantly enhances conceptual comprehension of nonlinear feedback mechanisms. Beyond its practical implications for the design of analog and digital memristive devices, CUDAynamics offers a scalable, open-source toolkit to aid researchers and engineers in exploring complex dynamical phenomena.
Khanov et al. (2026) studied this question.