<p>We propose a general fitting formula of energy spectrum of suprathermal particles, J=AE<sup>-&#946;1</sup>[1+(E/E<sub>0</sub>)<sup>&#945;</sup>]<sup>(&#946;1-&#946;2)/&#945;</sup>, where J is the particle flux (or intensity), E is the particle energy, A is the amplitude coefficient, E<sub>0</sub> represents the spectral break energy, &#945; (>0) describes the sharpness of energy spectral break around E<sub>0</sub>, and the power-law index &#946;<sub>1</sub> (&#946;<sub>2</sub>) gives the spectral shape before (after) the break. &#160;When &#945; tends to infinity (zero), this spectral formula becomes a classical double-power-law (logarithmic-parabola) spectrum. When both &#946;<sub>2</sub> and E<sub>0</sub> tend to infinity, this formula can be simplified to an Ellison-Ramaty-like equation. Under some other specific parameter conditions, this formula can be transformed to a Kappa or Maxwellian function. Considering &#160;the uncertainties both in particle intensity and energy, we fit this general formula well to the representative energy spectra of various suprathermal particle phenomena including solar energetic particles (electrons, protons, &#160;<sup>3</sup>He and heavier ions), shocked particles, anomalous cosmic rays, hard X-rays, solar wind suprathermal particles, etc. Therefore, this general spectrum fitting formula would help us to comparatively examine the energy spectrum of different suprathermal particle phenomena and understand their origin, acceleration and transportation.</p>