A new method based on semi-tensor product of matrices for solving Generalized Lyapunov Equation on the quaternion skew-filed