In this paper, we consider the split quaternion matrix equation X−Af(X)B=C, f(X)∈{X,XH,XiH,XjHXkH}. The H representation method has the characteristics of transforming a matrix with a special structure into a column vector with independent elements. By using the real representation of split quaternion matrices, H representation method, the Kronecker product of matrices and the Moore-Penrose generalized inverse, we convert the split quaternion matrix equation into the real matrix equation, and derive the sufficient and necessary conditions and the general solution expressions for the (skew) bisymmetric solution of the original equation. Moreover, we provide numerical algorithms and illustrate the efficiency of our method by two numerical examples.