Search

Finding the optimal parameters and functions of iterative methods is among the main problems of the Numerical Analysis. For this aim, a technique of the stochastic arithmetic (SA) is used to control of accuracy on Taylor-collocation method for solving first kind weakly regular integral equations (IEs). Thus, the CESTAC (Controle et Estimation Stochastique des Arrondis de Calculs) method is applied and instead of usual mathematical softwares the CADNA (Control of Accuracy and Debugging for Numerical Applications) library is used. Also, the convergence theorem of presented method is illustrated. In order to apply the CESTAC method we will prove a theorem that it will be our licence to use the new termination criterion instead of traditional absolute error. By using this theorem we can show that number of common significant digits (NCSDs) between two successive approximations are almost equal to NCSDs between exact and numerical solution. Finally, some examples are solved by using the Taylor-collocation method based on the CESTAC method. Several tables of numerical solutions based on the both arithmetics are presented. Comparison between number of iterations are demonstrated by using the floating point arithmetic (FPA) for different values of $\varepsilon$., Comment: 24 pages

An inverse problem in spectroscopy is considered. The objective is to restore the discrete spectrum from observed spectrum data, taking into account the spectrometer's line spread function. The problem is reduced to solution of a system of linear-nonlinear equations (SLNE) with respect to intensities and frequencies of the discrete spectral lines. The SLNE is linear with respect to lines' intensities and nonlinear with respect to the lines' frequencies. The integral approximation algorithm is proposed for the solution of this SLNE. The algorithm combines solution of linear integral equations with solution of a system of linear algebraic equations and avoids nonlinear equations. Numerical examples of the application of the technique, both to synthetic and experimental spectra, demonstrate the efficacy of the proposed approach in enabling an effective enhancement of the spectrometer's resolution., submitted to Applied Spectroscopy Journal

The linear PDE ${\mathbf B} {\mathbf L} (\frac{\partial}{\partial x}) u ={\mathbf L}_1(\frac{\partial}{\partial x})u +f(x)$ with nonclassic conditions on boundary $\partial \Omega$ is considered. Here ${\mathbf B}$ is linear noninvertible bounded operator acting from linear space $E$ into $E,$ $x=(t,x_1,\dots, x_m) \in \Omega, $ $\Omega \subset {\mathbb R}^{m+1}.$ It is assumed that ${\mathbf B}$ enjoys the skeleton decomposition ${\mathbf B}={\mathbf A}_1 {\mathbf A}_2,$ ${\mathbf A}_2 \in {\mathcal L}(E\rightarrow E_1),$ ${\mathbf A}_1 \in {\mathcal L}(E_1\rightarrow E)$ where $E_1$ is linear normed space. Differential operators ${\mathbf L}, \, {\mathbf L}_1$ are partial differential operators. In the concrete cases the domains of definition of operators ${\mathbf L}, {\mathbf L}_1$ consist of linear manifolds $E_{\partial}$ of sufficiently smooth abstract functions $u(x)$ with domain in $\Omega$ and their ranges in $E,$ which satisfy certain system of homogeneous boundary conditions. The abstract function $f: \Omega \subset {\mathbb R}^{m+1} \rightarrow E $ is assumed to be given. It is requested to find the solution $u: \Omega \subset {\mathbb R}^{m+1} \rightarrow E_{\partial},$ which satisfy certain condition on boundary $\partial \Omega.$ The concept of a skeleton chains is introduced as sequence of linear operators ${\mathbf B}_i \in {\mathcal L}(E_i \rightarrow E_i), \, i=1,2,\dots, p,$ where $E_i$ are linear spaces corresponding to the skeleton decomposition of operator ${\mathbf B}.$ It is assumed that irreversible operator ${\mathbf B}$ generates skeleton chain of the finite length $p.$ The problem is reduced to a regular split system with respect to higher order derivative terms with certain initial and boundary conditions.

Development of reliable methods for optimised energy storage and generation is one of the most imminent challenges in moder power systems. In this paper an adaptive approach to load leveling problem using novel dynamic models based on the Volterra integral equations of the first kind with piecewise continuous kernels. These integral equations efficiently solve such inverse problem taking into account both the time dependent efficiencies and the availability of generation/storage of each energy storage technology. In this analysis a direct numerical method is employed to find the least-cost dispatch of available storages. The proposed collocation type numerical method has second order accuracy and enjoys self-regularization properties, which is associated with confidence levels of system demand. This adaptive approach is suitable for energy storage optimisation in real time. The efficiency of the proposed methodology is demonstrated on the Single Electricity Market of Republic of Ireland and Sakhalin island in the Russian Far East., arXiv admin note: text overlap with arXiv:1507.06484

We propose the generalized quadrature methods for numerical solution of singular integral equation of Abel type. We overcome the singularity using the analytical calculation of the singular integral expression. The problem of solution of singular integral equation is reduced to nonsingular system of linear algebraic equations without shift meshes techniques employment. We also propose generalized quadrature method for solution of Abel equation using the singular integral. Relaxed errors bounds are derived. In order to improve the accuracy we use Tikhonov regularization method. We demonstrate the efficiency of proposed techniques on infrared tomography problem. Numerical experiments show that it make sense to apply regularization in case of highly noisy sources only. That is due to the fact that singular integral equations enjoy selfregularization property.

We propose the numerical methods for solution of the weakly regular linear and nonlinear evolutionary (Volterra) integral equation of the first kind. The kernels of such equations have jump discontinuities along the continuous curves (endogenous delays) which starts at the origin. In order to linearize these equations we use the modified Newton-Kantorovich iterative process. Then for linear equations we propose two direct quadrature methods based on the piecewise constant and piecewise linear approximation of the exact solution. The accuracy of proposed numerical methods is $\mathcal{O}(1/N)$ and $\mathcal{O}(1/N^2)$ respectively. We also suggest a certain iterative numerical scheme enjoying the regularization properties. Furthermore, we adduce generalized numerical method for nonlinear equations. We employ the midpoint quadrature rule in all the cases. In conclusion we include several numerical examples in order to demonstrate the efficiency of proposed numerical methods, Submitted to APNUM Journal

The numerical method for solution of the weakly regular scalar Volterra integral equation of the 1st kind is proposed. The kernels of such equations have jump discontinuities on the continuous curves which starts at the origin. The mid-rectangular quadrature rule is employed. The accuracy of proposed numerical method is $\mathcal{O}(1/N).$

The sufficient conditions are obtained for existence of the main solution of the nonlinear Volterra integral equation of the second kind on the semi-axis and on a finite interval. The method for computation of this boundary interval is designed. Beyond such integral the solution has the blow-up. The efficiency of proposed technique is demonstrated on concrete examples.

In this text matrix Volterra integral equation of the first kind is addressed. It is assumed that kernels of the equation have jump discontinuities on non-intersecting curves. Such equations appear in the theory of evolving dynamic systems. Differentiation of such equations with jump discontinue kernels yields the new class of the Volterra integral equations with functionally perturbed argument. The algorithm for construction of the logarithmic power asymptotics of the desired continuous solutions is proposed. The theorem of existance of the parametric families of solutions is proved. Finally the sufficient conditions for existence and uniqueness of continuous solution are derived., Comment: submitted to Russian Mathematics Journal

Sufficient conditions for existence and uniqueness of the solution of the Volterra integral equations of the first kind with piecewise continuous kernels are derived in framework of Sobolev-Schwartz distribution theory. The asymptotic approximation of the parametric family of generalized solutions is constructed. The method for the solution's regular part refinement is proposed using the successive approximations method.

The sufficient conditions for existence and uniqueness of continuous solutions of the Volterra operator equations of the first kind with piecewise continuous kernel are derived. The asymptotic approximation of the parametric family of solutions are constructed in case of non-unique solution. The algorithm for the solution's improvement is proposed using the successive approximations method., Comment: 17 pages

The Volterra integral equations of the first kind with piecewise smooth kernel are considered. Such equations appear in the theory of optimal control of the evolving systems. The existence theorems are proved. The method for constructing approximations of parametric families of solutions of such equations is suggested. The parametric family of solutions is constructed in terms of a logarithmic-power asymptotics., Comment: 11 pages, bibl 14, Submitted to Studia Informatica Universalis