Multivariate time series forecasting method based on nonlinear spiking neural P systems and non-subsampled shearlet transform. (August 2022)