Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.

THE LINEAR SAMPLING METHOD FOR RECONSTRUCTING A PENETRABLE CAVITY WITH UNKNOWN EXTERNAL OBSTACLES

We consider the interior inverse scattering problem for recovering the shape of a penetrable partially coated cavity with external obstacles from the knowledge of measured scattered waves due to point sources.In the first part,we obtain the well-posedness of the direct scattering problem by the variational method.In the second part,we establish the mathematical basis of the linear sampling method to recover both the shape of the cavity,and the shape of the external obstacle,however the exterior transmission eigenvalue problem also plays a key role in the discussion of this paper.

Development of a Quantitative Prediction Support System Using the Linear Regression Method

The development of prediction supports is a critical step in information systems engineering in this era defined by the knowledge economy, the hub of which is big data. Currently, the lack of a predictive model, whether qualitative or quantitative, depending on a company’s areas of intervention can handicap or weaken its competitive capacities, endangering its survival. In terms of quantitative prediction, depending on the efficacy criteria, a variety of methods and/or tools are available. The multiple linear regression method is one of the methods used for this purpose. A linear regression model is a regression model of an explained variable on one or more explanatory variables in which the function that links the explanatory variables to the explained variable has linear parameters. The purpose of this work is to demonstrate how to use multiple linear regressions, which is one aspect of decisional mathematics. The use of multiple linear regressions on random data, which can be replaced by real data collected by or from organizations, provides decision makers with reliable data knowledge. As a result, machine learning methods can provide decision makers with relevant and trustworthy data. The main goal of this article is therefore to define the objective function on which the influencing factors for its optimization will be defined using the linear regression method.

Solving Multi-Objective Linear Programming Problem by Statistical Averaging Method with the Help of Fuzzy Programming Method

A multi-objective linear programming problem is made from fuzzy linear programming problem. It is due the fact that it is used fuzzy programming method during the solution. The Multi objective linear programming problem can be converted into the single objective function by various methods as Chandra Sen’s method, weighted sum method, ranking function method, statistical averaging method. In this paper, Chandra Sen’s method and statistical averaging method both are used here for making single objective function from multi-objective function. Two multi-objective programming problems are solved to verify the result. One is numerical example and the other is real life example. Then the problems are solved by ordinary simplex method and fuzzy programming method. It can be seen that fuzzy programming method gives better optimal values than the ordinary simplex method.

Stability analysis of a liquid crystal elastomer self-oscillator under a linear temperature field

Self-oscillating systems abound in the natural world and offer substantial potential for applications in controllers,micro-motors,medical equipments,and so on.Currently,numerical methods have been widely utilized for obtaining the characteristics of self-oscillation including amplitude and frequency.However,numerical methods are burdened by intricate computations and limited precision,hindering comprehensive investigations into self-oscillating systems.In this paper,the stability of a liquid crystal elastomer fiber self-oscillating system under a linear temperature field is studied,and analytical solutions for the amplitude and frequency are determined.Initially,we establish the governing equations of self-oscillation,elucidate two motion regimes,and reveal the underlying mechanism.Subsequently,we conduct a stability analysis and employ a multi-scale method to obtain the analytical solutions for the amplitude and frequency.The results show agreement between the multi-scale and numerical methods.This research contributes to the examination of diverse self-oscillating systems and advances the theoretical analysis of self-oscillating systems rooted in active materials.

Underdetermined direction of arrival estimation with nonuniform linear motion sampling based on a small unmanned aerial vehicle platform

Uniform linear array(ULA)radars are widely used in the collision-avoidance radar systems of small unmanned aerial vehicles(UAVs).In practice,a ULA's multi-target direction of arrival(DOA)estimation performance suffers from significant performance degradation owing to the limited number of physical elements.To improve the underdetermined DOA estimation performance of a ULA radar mounted on a small UAV platform,we propose a nonuniform linear motion sampling underdetermined DOA estimation method.Using the motion of the UAV platform,the echo signal is sampled at different positions.Then,according to the concept of difference co-array,a virtual ULA with multiple array elements and a large aperture is synthesized to increase the degrees of freedom(DOFs).Through position analysis of the original and motion arrays,we propose a nonuniform linear motion sampling method based on ULA for determining the optimal DOFs.Under the condition of no increase in the aperture of the physical array,the proposed method obtains a high DOF with fewer sampling runs and greatly improves the underdetermined DOA estimation performance of ULA.The results of numerical simulations conducted herein verify the superior performance of the proposed method.

A Study of EM Algorithm as an Imputation Method: A Model-Based Simulation Study with Application to a Synthetic Compositional Data

Compositional data, such as relative information, is a crucial aspect of machine learning and other related fields. It is typically recorded as closed data or sums to a constant, like 100%. The statistical linear model is the most used technique for identifying hidden relationships between underlying random variables of interest. However, data quality is a significant challenge in machine learning, especially when missing data is present. The linear regression model is a commonly used statistical modeling technique used in various applications to find relationships between variables of interest. When estimating linear regression parameters which are useful for things like future prediction and partial effects analysis of independent variables, maximum likelihood estimation (MLE) is the method of choice. However, many datasets contain missing observations, which can lead to costly and time-consuming data recovery. To address this issue, the expectation-maximization (EM) algorithm has been suggested as a solution for situations including missing data. The EM algorithm repeatedly finds the best estimates of parameters in statistical models that depend on variables or data that have not been observed. This is called maximum likelihood or maximum a posteriori (MAP). Using the present estimate as input, the expectation (E) step constructs a log-likelihood function. Finding the parameters that maximize the anticipated log-likelihood, as determined in the E step, is the job of the maximization (M) phase. This study looked at how well the EM algorithm worked on a made-up compositional dataset with missing observations. It used both the robust least square version and ordinary least square regression techniques. The efficacy of the EM algorithm was compared with two alternative imputation techniques, k-Nearest Neighbor (k-NN) and mean imputation (), in terms of Aitchison distances and covariance.

Solution of Nonlinear Integro Differential Equations by Two-Step Adomian Decomposition Method (TSAM)

The Adomian decomposition method (ADM) can be used to solve a wide range of problems and usually gets the solution in a series form. In this paper, we propose two-step Adomian Decomposition Method (TSAM) for nonlinear integro-differential equations that will facilitate the calculations. In this modification, compared to the standard Adomian decomposition method, the size of calculations was reduced. This modification also avoids computing Adomian polynomials. Numerical results are given to show the efficiency and performance of this method.

Image reconstruction based on total-variation minimization and alternating direction method in linear scan computed tomography

Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.

MAOR method for the generalized—order linear complementarity problems

The modified AOR method for solving linear complementarity problem(LCP(M,p))was proposed in literature,with some convergence results.In this paper,we considered the MAOR method for generalized-order linear complementarity problem(ELCP(M,N,p,q)),where M,N are nonsingular matrices of the following form:M=[D11H1K1D2],N=[D12H2K2D22],D11,D12,D21 and D22 are square nonsingular diagonal matrices.

A linearly-independent higher-order extended numerical manifold method and its application to multiple crack growth simulation

The numerical manifold method(NMM)can be viewed as an inherent continuous-discontinuous numerical method,which is based on two cover systems including mathematical and physical covers.Higher-order NMM that adopts higher-order polynomials as its local approximations generally shows higher precision than zero-order NMM whose local approximations are constants.Therefore,higherorder NMM will be an excellent choice for crack propagation problem which requires higher stress accuracy.In addition,it is crucial to improve the stress accuracy around the crack tip for determining the direction of crack growth according to the maximum circumferential stress criterion in fracture mechanics.Thus,some other enriched local approximations are introduced to model the stress singularity at the crack tip.Generally,higher-order NMM,especially first-order NMM wherein local approximations are first-order polynomials,has the linear dependence problems as other partition of unit(PUM)based numerical methods does.To overcome this problem,an extended NMM is developed based on a new local approximation derived from the triangular plate element in the finite element method(FEM),which has no linear dependence issue.Meanwhile,the stresses at the nodes of mathematical mesh(the nodal stresses in FEM)are continuous and the degrees of freedom defined on the physical patches are physically meaningful.Next,the extended NMM is employed to solve multiple crack propagation problems.It shows that the fracture mechanics requirement and mechanical equilibrium can be satisfied by the trial-and-error method and the adjustment of the load multiplier in the process of crack propagation.Four numerical examples are illustrated to verify the feasibility of the proposed extended NMM.The numerical examples indicate that the crack growths simulated by the extended NMM are in good accordance with the reference solutions.Thus the effectiveness and correctness of the developed NMM have been validated.

ON THE MONOTONE CONVERGENCE OF THE PROJECTED ITERATION METHODS FOR LINEAR COMPLEMENTARITY PROBLEMS

Under suitable conditions,the monotone convergence about the projected iteration method for solving linear complementarity problem is proved and the influence of the involved parameter matrix on the convergence rate of this method is investigated.

On-Machine Measurement of the Straightness and Tilt Errors of a Linear Slideway Using a New Four-Sensor Method

Although there are some multi-sensor methods for measuring the straightness and tilt errors of a linear slideway, they need to be further improved in some aspects, such as suppressing measurement noise and reducing precondition.In this paper, a new four-sensor method with an improved measurement system is proposed to on-machine separate the straightness and tilt errors of a linear slideway from the sensor outputs, considering the influences of the reference surface profile and the zero-adjustment values. The improved system is achieved by adjusting a single sensor to di erent positions. Based on the system, a system of linear equations is built by fusing the sensor outputs to cancel out the e ects of the straightness and tilt errors. Three constraints are then derived and supplemented into the linear system to make the coe cient matrix full rank. To restrain the sensitivity of the solution of the linear system to the measurement noise in the sensor outputs, the Tikhonov regularization method is utilized. After the surface profile is obtained from the solution, the straightness and tilt errors are identified from the sensor outputs. To analyze the e ects of the measurement noise and the positioning errors of the sensor and the linear slideway, a series of computer simulations are carried out. An experiment is conducted for validation, showing good consistency. The new four-sensor method with the improved measurement system provides a new way to measure the straightness and tilt errors of a linear slideway, which can guarantee favorable propagations of the residuals induced by the noise and the positioning errors.

Piecewise linear recursive convolution FDTD method for magnetized plasmas

The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method greatly improves accuracy over the original recursive convolution （RC） FDTD approach but retains its speed and efficiency advantages. A PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time is presented, enabled the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations the reflection and transmission coefficients through a magnetized plasma layer. The results show that the PLRC-FDTD method has significantly improved the accuracy over the original RC method.

Linearization Learning Method of BP Neural Networks

Feedforward multi layer neural networks have very strong mapping capability that is based on the non linearity of the activation function, however, the non linearity of the activation function can cause the multiple local minima on the learning error surfaces, which affect the learning rate and solving optimal weights. This paper proposes a learning method linearizing non linearity of the activation function and discusses its merits and demerits theoretically.

ASYNCHRONOUS RELAXED ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS OF EQUATIONS

In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.

THE STABILITY OF LINEAR MULTISTEP METHODS FOR SYSTEMS OF DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.

MHD flow of upper-convected Maxwell fluid over porous stretching sheet using successive Taylor series linearization method

This paper investigates the magnetohydrodynamic （MHD） boundary layer flow of an incompressible upper-convected Maxwell （UCM） fluid over a porous stretching surface. Similarity transformations are used to reduce the governing partial differential equations into a kind of nonlinear ordinary differential equations. The nonlinear prob- lem is solved by using the successive Taylor series linearization method （STSLM）. The computations for velocity components are carried out for the emerging parameters. The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem.