In this work, we determine the non-perturbative fragmentation function (FF) of B-mesons (meson including a bottom quark flavor), phenomenologically, at next-to-leading order QCD perturbation theory using data from electron-positron annihilation. In this regard, we apply all current experimental data from OPAL, ALEPH, SLD2002 and DELPHI collaborations. We also apply the Hessian method to calculate the error in determination of FFs. As an application, using our extracted results we shall make our predictions for energy distribution of B-mesons produced through top quark decay. Study of the B-meson energy spectrum can be considered as a convenient channel to search for the properties of top quarks at CERN.