New Predictor-Corrector Methods Based on Exponential Time Differencing forSystems of Nonlinear Differential Equations

We present the new predictor-corrector methods for systems of nonlinear differential equations, based on the method of exponential time differencing. We compare the present schemes with the explicit multistep exponential time differencing and Adams–Bashforth–Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.

Three-Step Predictor-Corrector of Exponential Fitting Method for Nonlinear Schroedinger Equations

We develop the three-step explicit and implicit schemes of exponential fitting methods. We use the three- step explicit exponential fitting scheme to predict an approximation, then use the three-step implicit exponential fitting scheme to correct this prediction. This combination is called the three-step predictor-corrector of exponential fitting method. The three-step predictor-corrector of exponential fitting method is applied to numerically compute the coupled nonlinear Schroedinger equation and the nonlinear Schroedinger equation with varying coefficients. The numerical results show that the scheme is highly accurate.

ON CONSTRUCTION OF HIGH ORDER EXPONENTIALLY FITTED METHODS BASED ON PARAMETERIZED RATIONAL APPROXIMATIONS TO exp

By the discussion of the formula and properties of (4,4) parametric form rational approximation to function exp(q), the fourth order derivative one_step exponentially fitted method and the third order derivative hybrid one_step exponentially fitted method are presented, their order p satisfying 6≤p≤8. The necessary and sufficient conditions for the two methods to be A_ stable are given. Finally, for the fourth order derivative method, the error bound and the necessary and sufficient conditions for it to be median are discussed.

PLANE WAVES NUMERICAL STABILITY OF SOME EXPLICIT EXPONENTIAL METHODS FOR CUBIC SCHRODINGER EQUATION

Numerical stability when integrating plane waves of cubic SchrSdinger equation is thor- oughly analysed for some explicit exponential methods. We center on the following second- order methods： Strang splitting and Lawson method based on a one-parameter family of 2-stage 2nd-order explicit Runge-Kutta methods. Regions of stability are plotted and numerical results are shown which corroborate the theoretical results. Besides, a tech- nique is suggested to avoid the possible numerical instabilities which do not correspond to continuous ones.

Improved precise integration method for differential Riccati equation

An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the

Improvement of scattering correction for in situ coastal and inland water absorption measurement using exponential fitting approach

The absorption coefficient of water is an important bio-optical parameter for water optics and water color remote sensing. However, scattering correction is essential to obtain accurate absorption coefficient values in situ using the nine-wavelength absorption and attenuation meter AC9. Establishing the correction always fails in Case 2 water when the correction assumes zero absorption in the near-infrared(NIR) region and underestimates the absorption coefficient in the red region, which affect processes such as semi-analytical remote sensing inversion. In this study, the scattering contribution was evaluated by an exponential fitting approach using AC9 measurements at seven wavelengths(412, 440, 488, 510, 532, 555, and 715 nm) and by applying scattering correction. The correction was applied to representative in situ data of moderately turbid coastal water, highly turbid coastal water, eutrophic inland water, and turbid inland water. The results suggest that the absorption levels in the red and NIR regions are significantly higher than those obtained using standard scattering error correction procedures. Knowledge of the deviation between this method and the commonly used scattering correction methods will facilitate the evaluation of the effect on satellite remote sensing of water constituents and general optical research using different scatteringcorrection methods.

2020—2025年广东省医疗机构床位需求预测

目的预测2020—2025年广东省医疗机构的床位需求总量。方法基于卫生服务需求法与Holt双参数指数平滑模型,结合年龄别人口数据预测床位需求。结果2025年,广东省住院人数为2425.11万人,床位需求数为70.04万张,每千常住人口床位需求数为5.63张、每千常住人口拥有床位数5.55张,供需比例为98.58%。预测模型的平均百分误差为1.63%(标准差=±1.90%,均方根误差=20.63)。结论结合人口的年龄结构进行预测结果更稳定、误差更小。2020年,广东省的床位配置量基本能满足床位需求,供需较为平衡。但2024年床位需求将超过床位配置总量。未来,广东省应加大床位资源的投入力度,提高基层卫生机构的床位利用率,全面落实分级诊疗制度。

EXPONENTIAL TIKHONOV REGULARIZATION METHOD FOR SOLVING AN INVERSE SOURCE PROBLEM OF TIME FRACTIONAL DIFFUSION EQUATION

In this paper,we mainly study an inverse source problem of time fractional diffusion equation in a bounded domain with an over-specified terminal condition at a fixed time.A novel regularization method,which we call the exponential Tikhonov regularization method with a parameter γ,is proposed to solve the inverse source problem,and the corresponding convergence analysis is given under a-priori and a-posteriori regularization parameter choice rules.Whenγis less than or equal to zero,the optimal convergence rate can be achieved and it is independent of the value of γ.However,when γ is greater than zero,the optimal convergence rate depends on the value of γ which is related to the regularity of the unknown source.Finally,numerical experiments are conducted for showing the effectiveness of the proposed exponential regularization method.

EXPONENTIAL FOURIER COLLOCATION METHODS FOR SOLVING FIRST-ORDER DIFFERENTIAL EQUATIONS

In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.

Exponential Runge-Kutta Methods for the Multispecies Boltzmann Equation

This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared to the single species Boltzmann equation that the method was originally applied on,this set of equation presents a new difficulty that comes from the lack of local conservation laws due to the interaction between different species.Hence extra stiff nonlinear source terms need to be treated properly to maintain the accuracy and the AP property.The method we propose does not contain any nonlinear nonlocal implicit solver,and can capture the hydrodynamic limit with time step and mesh size independent of the Knudsen number.We prove the positivity and strong AP properties of the scheme,which are verified by two numerical examples.

Predicting LTE Throughput Using Traffic Time Series

Throughput prediction is essential for congestion control and LTE network management. In this paper, the autoregressive integrated moving average （ARIMA） model and exponential smoothing model are used to predict the throughput in a single cell and whole region in an LTE network. The experimental results show that these two models perform differently in both scenarios. The ARIMA model is better than the exponential smoothing model for predicting throughput on weekdays in a whole region. The exponential smoothing model is better than the ARIMA model for predicting throughput on weekends in a whole region. The exponential smoothing model is better than the ARIMA model for predicting throughput in a single cell. In these two LTE network scenarios, throughput prediction based on traffic time series leads to more efficient resource management and better QoS.

AN EXPONENTIAL WAVE INTEGRATOR PSEUDOSPECTRAL METHOD FOR THE SYMMETRIC REGULARIZED-LONG-WAVE EQUATION

An efficient and accurate exponential wave integrator Fourier pseudospectral （EWI-FP） method is proposed and analyzed for solving the symmetric regularized-long-wave （SRLW） equation, which is used for modeling the weakly nonlinear ion acoustic and space-charge waves. The numerical method here is based on a Gautschi-type exponential wave integrator for temporal approximation and the Fourier pseudospectral method for spatial discretization. The scheme is fully explicit and efficient due to the fast Fourier transform. Numerical analysis of the proposed EWI-FP method is carried out and rigorous error estimates are established without CFL-type condition by means of the mathematical induction. The error bound shows that EWI-FP has second order accuracy in time and spectral accuracy in space. Numerical results are reported to confirm the theoretical studies and indicate that the error bound here is optimal.

Heat and mass flux estimation of modern seafloor hydrothermal activity

Research on heat and mass flux yielded by modern seafloor hydrothermal activity is very important, because it is involved not only in the base of ocean environment research, but also in the historical evolution of seawater properties. Currently, estimating heat flux is based on the observation data of hydrothermal smokers, low-temperature diffusive flow and mid-ocean ridge mainly. But there are some faults, for example, there is lack of a concurrent conductive item in estimating the heat flux by smokers and the error between the half-space cooling model and the observation data is too large. So, three kinds of methods are applied to re-estimating the heat flux of hydrothermal activity resepectively, corresponding estimation is 97. 359 GW by hydrothermal smoker and diffusive flow, 84.895 GW by hydrothermal plume, and 4. 11 TW by exponential attenuation method put forward by this paper. Research on mass flux estimation is relatively rare, the main reason for this is insufficient field observation data. Mass fluxes of different elements are calculated using hydrothermal vent fluid data from the TAG hydrothermal area on the Mid-Atlantic Ridge for the first time. Difference of estimations by different methods reflects the researching extent of hydrothermal activity, and systematically in - situ observation will help to estimate the contribution of hydrothermal activity to ocean chemical environment, ocean circulation and global climate precisely.

CONVERGENCE AND MEAN-SQUARE STABILITY OF EXPONENTIAL EULER METHOD FOR SEMI-LINEAR STOCHASTIC DELAY INTEGRO-DIFFERENTIAL EQUATIONS

In this paper,the numerical methods for semi-linear stochastic delay integro-difFerential equations are studied.The uniqueness,existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained.Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied.It is proved that the exponential Euler method is convergent with strong order 1/2 and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size.In addition,numerical experiments are presented to confirm the theoretical results.

EXPONENTIAL FITTED METHODS FOR THE NUMEMCAL SOLUTION OF THE SCHRDINGER EQUATION

A new sixth-order Runge-Kutta type method is developed for the numericalintegration of the radial Schrodinger equation and of the coupled differential equa-tions of the Schrodinger type. The formula developed contains certain free pa-rameters which allows it to be fitted automatically to exponential functions. Wegive a comparative error analysis with other sixth order exponentially fitted meth-ods. The theoretical and numerical results indicate that the new method is moreaccurate than the other exponentially fitted methods.

Efficient and Stable Exponential Runge-Kutta Methods for Parabolic Equations

In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations.By incorporating the linear splitting technique into the explicit exponential Runge-Kutta schemes,we are able to greatly improve the numerical stability.The proposed numerical methods could be fast implemented through use of decompositions of compact spatial difference operators on a regular mesh together with discrete fast Fourier transform techniques.The exponential Runge-Kutta schemes are easy to be adopted in adaptive temporal approximations with variable time step sizes,as well as applied to stiff nonlinearity and boundary conditions of different types.Linear stabilities of the proposed schemes and their comparison with other schemes are presented.We also numerically demonstrate accuracy,stability and robustness of the proposed method through some typical model problems.

An Exact l_(1) Exponential Penalty Function Method for Multiobjective Optimization Problems with Exponential-Type Invexity

The purpose of this paper is to devise exact l_(1) exponential penalty function method to solve multiobjective optimization problems with exponentialtype invexity.The conditions governing the equivalence of the(weak)efficient solutions to the vector optimization problem and the(weak)efficient solutions to associated unconstrained exponential penalized multiobjective optimization problem are studied.Examples are given to illustrate the obtained results.

An exponential for the solutions model of collocation method of the HIV infection CD4＋T cells

In this paper, an exponential method is presented for the approximate solutions of the HIV infection model of CD4＋T. The method is based on exponential polynomi- als and collocation points. This model problem corresponds to a system of nonlinear ordinary differential equations. Matrix relations are constructed for the exponential functions. By aid of these matrix relations and the collocation points, the proposed technique transforms the model problem into a system of nonlinear algebraic equations. By solving the system of the algebraic equations, the unknown coefficients are com- puted and thus the approximate solutions are obtained. The applications of the method for the considered problem are given and the comparisons are made with the other methods.

Improved error estimates for a modified exponential Euler method for the semilinear stochastic heat equation with rough initial data

A class of stochastic Besov spaces BpL^(2)(Ω;˙H^(α)(O)),1≤p≤∞andα∈[−2,2],is introduced to characterize the regularity of the noise in the semilinear stochastic heat equation du−Δudt=f(u)dt+dW(t),under the following conditions for someα∈(0,1]:||∫_(0)^(t)e−(t−s)^(A)dW(s)||L^(2)(Ω;L^(2)(O))≤C^(t^(α/2))and||∫_(0)^(t)e−(t−s)^(A)dW(s)||_B^(∞)L^(2)(Ω:H^(α)(O))≤C..The conditions above are shown to be satisfied by both trace-class noises(withα=1)and one-dimensional space-time white noises(withα=1/2).The latter would fail to satisfy the conditions withα=1/2 if the stochastic Besov norm||·||B∞L^(2)(Ω;˙H^(α)(O))is replaced by the classical Sobolev norm||·||L^(2)(Ω;˙Hα(O)),and this often causes reduction of the convergence order in the numerical analysis of the semilinear stochastic heat equation.In this paper,the convergence of a modified exponential Euler method,with a spectral method for spatial discretization,is proved to have orderαin both the time and space for possibly nonsmooth initial data in L^(4)(Ω;˙H^(β)(O))withβ>−1,by utilizing the real interpolation properties of the stochastic Besov spaces and a class of locally refined stepsizes to resolve the singularity of the solution at t=0.

A Fitted Numerov Method for Singularly Perturbed Parabolic Partial Differential Equation with a Small Negative Shift Arising in Control Theory

In this paper,a fitted Numerov method is constructed for a class of singularly perturbed one-dimensional parabolic partial differential equations with a small negative shift in the temporal variable.Similar boundary value problems are associated with a furnace used to process a metal sheet in control theory.Here,the study focuses on the effect of shift on the boundary layer behavior of the solution via finite difference approach.When the shift parameter is smaller than the perturbation parameter,the shifted term is expanded in Taylor series and an exponentially fitted tridiagonal finite difference scheme is developed.The proposed finite difference scheme is unconditionally stable.When the shift parameter is larger than the perturbation parameter,a special type of mesh is used for the temporal variable so that the shift lies on the nodal points and an exponentially fitted scheme is developed.This scheme is also unconditionally stable.The applicability of the proposed methods is demonstrated by means of two examples.