In this paper, we present a novel extension of the well-known split-step Fourier transform (SSFT) approach for solving the one-dimensional nonlinear Schrödinger equation (NLSE), which incorporates the fiber loss term. While this essential equation governs the pulse propagation in a lossy optical fiber, it is not supported by an exact analytical solution. In this regard, extended versions of the Fourier pseudospectral method (FPSM) and Hopscotch method (HSM) are effectively established as well to cope with the fiber losses effects associated with the pulses’ propagation through the fiber optics, and thus, numerous comparisons are exhaustively conducted among these three compelling numerical approaches to validate their reliability, stability, and accuracy. Based on this, the MATLAB numerical findings bolster that the extended version of the SSFT approach demonstrates superior performance over the other suggested schemes in simulating the solitons propagation in a lossy optical fiber.