1,636 research outputs found
Centrosymmetric Matrices in the Sinc Collocation Method for Sturm-Liouville Problems
Recently, we used the Sinc collocation method with the double exponential
transformation to compute eigenvalues for singular Sturm-Liouville problems. In
this work, we show that the computation complexity of the eigenvalues of such a
differential eigenvalue problem can be considerably reduced when its operator
commutes with the parity operator. In this case, the matrices resulting from
the Sinc collocation method are centrosymmetric. Utilizing well known
properties of centrosymmetric matrices, we transform the problem of solving one
large eigensystem into solving two smaller eigensystems. We show that only
1/(N+1) of all components need to be computed and stored in order to obtain all
eigenvalues, where (2N+1) corresponds to the dimension of the eigensystem. We
applied our result to the Schr\"odinger equation with the anharmonic potential
and the numerical results section clearly illustrates the substantial gain in
efficiency and accuracy when using the proposed algorithm.Comment: 11 pages, 4 figure
Two-parameter Sturm-Liouville problems
This paper deals with the computation of the eigenvalues of two-parameter
Sturm- Liouville (SL) problems using the Regularized Sampling Method, a method
which has been effective in computing the eigenvalues of broad classes of SL
problems (Singular, Non-Self-Adjoint, Non-Local, Impulsive,...). We have shown,
in this work that it can tackle two-parameter SL problems with equal ease. An
example was provided to illustrate the effectiveness of the method.Comment: 9 page
Information Transmission using the Nonlinear Fourier Transform, Part II: Numerical Methods
In this paper, numerical methods are suggested to compute the discrete and
the continuous spectrum of a signal with respect to the Zakharov-Shabat system,
a Lax operator underlying numerous integrable communication channels including
the nonlinear Schr\"odinger channel, modeling pulse propagation in optical
fibers. These methods are subsequently tested and their ability to estimate the
spectrum are compared against each other. These methods are used to compute the
spectrum of various signals commonly used in the optical fiber communications.
It is found that the layer-peeling and the spectral methods are suitable
schemes to estimate the nonlinear spectra with good accuracy. To illustrate the
structure of the spectrum, the locus of the eigenvalues is determined under
amplitude and phase modulation in a number of examples. It is observed that in
some cases, as signal parameters vary, eigenvalues collide and change their
course of motion. The real axis is typically the place from which new
eigenvalues originate or are absorbed into after traveling a trajectory in the
complex plane.Comment: Minor updates to IEEE Transactions on Information Theory, vol. 60,
no. 7, pp. 4329--4345, July 201
A system of ODEs for a Perturbation of a Minimal Mass Soliton
We study soliton solutions to a nonlinear Schrodinger equation with a
saturated nonlinearity. Such nonlinearities are known to possess minimal mass
soliton solutions. We consider a small perturbation of a minimal mass soliton,
and identify a system of ODEs similar to those from Comech and Pelinovsky
(2003), which model the behavior of the perturbation for short times. We then
provide numerical evidence that under this system of ODEs there are two
possible dynamical outcomes, which is in accord with the conclusions of
Pelinovsky, Afanasjev, and Kivshar (1996). For initial data which supports a
soliton structure, a generic initial perturbation oscillates around the stable
family of solitons. For initial data which is expected to disperse, the finite
dimensional dynamics follow the unstable portion of the soliton curve.Comment: Minor edit
- …