We have performed a numerical investigation of the influence of disorder on the dynamical non-equilibrium evolution of a 3D site-diluted Ising model from a low-temperature initial state with magnetization m0 = 1. It is shown that two-time dependences of the autocorrelation and integrated response functions for systems with spin concentrations p = 1.0, 0.95, 0.8, 0.6 and 0.5 demonstrate ageing properties with anomalous slowing-down relaxation and violation of the fluctuation-dissipation ratio. It was revealed that during non-equilibrium critical dynamics in the long-time regime the autocorrelation functions for diluted systems are extremely slow due to the pinning of domain walls on impurity sites. We have found that the fluctuation-dissipation ratio for diluted systems with spin concentration p < 1 while the pure system is characterized by . The autocorrelation function power-law delay becomes the same as for the time dependence of the magnetization in the critical point and is characterized by exponent . Also, for diluted systems we reveal memory effects for critical evolution in the ageing regime with realization of cyclic temperature change and quenching at .