Computing quaternion matrix pseudoinverse with zeroing neural networks