  Bosonization approach to the classical supersymmetric systems is presented.
By introducing the multi-fermionic parameters in the expansions of the
superfields, the $\mathcal {N}=1$ supersymmetric KdV (sKdV) equations are
transformed to a system of coupled bosonic equations. The method can be applied
to any fermionic systems. By solving the coupled bosonic equations, some novel
types of exact solutions can be explicitly obtained. Especially, the richness
of the localized excitations of the supersymmetric integrable system are
discovered. The rich multi-soliton solutions obtained here have not yet been
obtained by using other methods. Unfortunately, the traditional known
multi-soliton solutions can also not be obtained by the bosonization approach
of this paper. Some open problems on the bosonization of the supersymmetric
integrable models are proposed in the both classical and quantum levels.
