In this paper, a class of (3+1)-dimensional generalized Korteweg-de Vries-Zakharov-Kuznetsov (KdV-ZK) equation is studied by utilizing the bifurcation theory of the planar dynamical systems and the Fan sub-function method. This model can be used to explain the effects of magnetic fields on weakly nonlinear ion-acoustic waves investigated in plasma fields composed of cold and hot electrons. Under the different parameter conditions, the phase portraits and bifurcations are derived, and new exact solutions including soliton, periodic, kink and breaking wave solutions for the model are constructed. Moreover, some exact solutions, which contain soliton, kink, trigonometric function, hyperbolic function, Jacobi elliptic function solutions, are derived via the improved Fan sub-function method. The types of solutions obtained completely correspond to the types of the orbits acquired
above, which verifies the validity of the method. Finally, the physical structures of some exact solutions are analyzed in graphical forms.