Abstract:A numerical method to solve the direct scattering problem for the Zakharov-Shabat system associated to the initial value problem for the nonlinear Schroedinger equation is proposed. The method involves the numerical solution of Volterra integral systems with structured kernels and the identification of coefficients and parameters appearing in monomial-exponential sums. Numerical experiments confirm the effectiveness of the proposed technique.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to propose a numerical method to calculate the scattering data of the Zakharov - Shabat (ZS) system related to the initial - value problem of the nonlinear Schrödinger (NLS) equation. Specifically, the author aims to solve the following problems:
1. **Numerical solution of the direct scattering problem**: Given the initial potential \(u_0(x)\), obtain the initial scattering data by solving the ZS system. These scattering data include the transmission coefficient \(T(\lambda)\), the left and right reflection coefficients \(L(\lambda)\) and \(R(\lambda)\), the discrete spectrum (i.e., the bound states \(\{\lambda_j\}\) and their multiplicities \(\{m_j\}\)), and the normalization constants \(\{(\Gamma^\ell_{j,s}), (\Gamma^r_{j,s})\}\).
2. **Numerical solution of the Volterra integral equation**: In order to calculate the above - mentioned scattering data, it is necessary to solve the Volterra integral equation with a structured kernel and identify the coefficients and parameters appearing in the polynomial - exponential sum.
3. **Calculation of the Marchenko kernel**: Based on the auxiliary function, solve two Volterra integral equations to obtain the left Marchenko kernel \(\Omega_\ell(\alpha)\) and the right Marchenko kernel \(\Omega_r(\alpha)\).
4. **Calculation of the Fourier transform**: By solving the Fredholm integral equation, calculate the inverse Fourier transform \(\rho(\alpha)\) of the reflection coefficients \(R(\lambda)\) and \(L(\lambda)\) and the Fourier transform \(\ell(\alpha)\).
5. **Algorithm stability and efficiency**: Ensure the effectiveness and stability of the proposed numerical method when dealing with smooth initial values, and at the same time consider how to deal with the situation of non - smooth or jump - discontinuous initial values.
The main contribution of the paper lies in providing a novel and effective numerical method to calculate the complete scattering data, which also has independent application value in some engineering fields. In addition, this method provides the necessary input data for subsequent steps (such as time evolution and solving the Marchenko equation).