Applications of Mathematics最新文献

We investigate the Cauchy problem of the one dimensional Maxwell-Schrödinger (MS) system under the Lorenz gauge condition. Different from the classical case, we consider the electromagnetic and electrostatic potentials which are growing at space infinity. More precisely, the electrostatic potential is allowed to grow linearly, while for the electromagnetic potential the growth is sublinear. Based on the energy estimates and the gauge transformation, we prove the global existence and the uniqueness of the weak solutions to this system.

The aim of this article is to investigate the well-posedness, stability of solutions to the time-dependent Maxwell’s equations for electric field in conductive media in continuous and discrete settings, and study convergence analysis of the employed numerical scheme. The situation we consider would represent a physical problem where a subdomain is emerged in a homogeneous medium, characterized by constant dielectric permittivity and conductivity functions. It is well known that in these homogeneous regions the solution to the Maxwell’s equations also solves the wave equation, which makes computations very efficient. In this way our problem can be considered as a coupling problem, for which we derive stability and convergence analysis. A number of numerical examples validate theoretical convergence rates of the proposed stabilized explicit finite element scheme.

We consider the inverse nodal problem for Sturm-Liouville (S-L) equation with frozen argument. Asymptotic behaviours of eigenfunctions, nodal parameters are represented in two cases and numerical algorithms are produced to solve the given problems. Subsequently, solution of inverse nodal problem is calculated by the second Chebyshev wavelet method (SCW), accuracy and effectiveness of the method are shown in some numerical examples.

The Sturm-Liouville eigenvalue problem is symmetric if the coefficients are even functions and the boundary conditions are symmetric. The eigenfunction is expressed in terms of orthonormal bases, which are constructed in a linear space of trial functions by using the Gram-Schmidt orthonormalization technique. Then an n-dimensional matrix eigenvalue problem is derived with a special matrix A:= [a_ij], that is, a_ij = 0 if i + j is odd.

Based on the product formula, an integration method with a fictitious time, namely the fictitious time integration method (FTIM), is developed to obtain the higher-index eigenvalues. Also, we recover the symmetric potential function q(x) in the Sturm-Liouville operator by specifying a few lower-index eigenvalues, based on the product formula and the Newton iterative method.

We multiply both sides of the complex symmetric linear system Ax = b by 1 − iω to obtain a new equivalent linear system, then a dual-parameter double-step splitting (DDSS) method is established for solving the new linear system. In addition, we present an upper bound for the spectral radius of iteration matrix of the DDSS method and obtain its quasi-optimal parameter. Theoretical analyses demonstrate that the new method is convergent when some conditions are satisfied. Some tested examples are given to illustrate the effectiveness of the proposed method.

We study an n-dimensional system of ordinary differential equations with a constant matrix, a relay-type nonlinearity, and an external disturbance in the right-hand side. We consider a nonideal relay characteristic. The external disturbance is described by the product of an exponential function and a sine function with an initial phase as a parameter. We assume the matrix of the linear part and the vector at the relay characteristic such that, by a nonsingular transformation, the system is reduced to the form with the diagonal matrix and the vector being opposite to the unit vector. We establish a necessary and sufficient condition for the existence of two-point oscillatory solutions, i.e., the solutions with two fixed points on the hyperplanes of the relay switching in phase space. Also, we give the sufficient conditions under which such solutions do not exist. We provide a supporting example, which demonstrates how to apply the obtained results.

Initial value problem for three dimensional (3D) elastodynamic system in two dimensional (2D) inhomogeneous quasicrystals is considered. An analytical method is studied for the solution of this problem. The system is written in terms of Fourier images of displacements with respect to lateral variables. The resulting problem is reduced to integral equations of the Volterra type. Finally, using Paley Wiener theorem it is shown that the solution of the initial value problem can be found by the inverse Fourier transform. A numerical example is considered for the comparison of the exact solution with the computed solution obtained by using the method.

To have accuracy in the extracted information is the goal of the reliability theory investigation. In information theory, varentropy has recently been proposed to describe and measure the degree of information dispersion around entropy. Theoretical investigation on varentropy of past life has been initiated, however details on its stochastic properties are yet to be discovered. In this paper, we propose a novel stochastic order and introduce new classes of life distributions based on past varentropy. Further, we illustrate some of its applications in reliability modeling and in the diversity measure of Boltzmann distribution.

This paper deals with the problem of estimating a covariance matrix from the data in two classes: (1) good data with the covariance matrix of interest and (2) contamination coming from a Gaussian distribution with a different covariance matrix. The ridge penalty is introduced to address the problem of high-dimensional challenges in estimating the covariance matrix from the two-class data model. A ridge estimator of the covariance matrix has a uniform expression and keeps positive-definite, whether the data size is larger or smaller than the data dimension. Furthermore, the ridge parameter is tuned through a cross-validation procedure. Lastly, the proposed ridge estimator is verified with better performance than the existing estimator from the data in two classes and the traditional ridge estimator only from the good data.