Bi-Directional Algorithm for Computing Discrete Spectral Amplitudes in the NFT

The nonlinear Fourier transform represents a signal in terms of its continuous spectrum, discrete eigenvalues, and the corresponding discrete spectral amplitudes. This paper presents a new bi-directional algorithm for computing the discrete spectral amplitudes, which addresses the significant problem of rounding errors inherent in previously known techniques. We use the proposed method to obtain accurate spectral-domain noise statistics for 2-soliton signals using numerical simulation.