Application of neural network model coupling with the partial least-squares method for forecasting watre yield of mine

Scientific forecasting water yield of mine is of great significance to the safety production of mine and the colligated using of water resources. The paper established the forecasting model for water yield of mine, combining neural network with the partial least square method. Dealt with independent variables by the partial least square method, it can not only solve the relationship between independent variables but also reduce the input dimensions in neural network model, and then use the neural network which can solve the non-linear problem better. The result of an example shows that the prediction has higher precision in forecasting and fitting.

Adaptive multi-step piecewise interpolation reproducing kernel method for solving the nonlinear time-fractional partial differential equation arising from financial economics

This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.

LEAST-SQUARES METHOD-BASED FEATURE FITTING AND EXTRACTION IN REVERSE ENGINEERING

The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.

ON THE BREAKDOWNS OF THE GALERKIN AND LEAST-SQUARES METHODS

The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.

NEGATIVE NORM LEAST-SQUARES METHODS FOR THE INCOMPRESSIBLE MAGNETOHYDRODYNAMIC EQUATIONS

The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.

A partial least-squares regression approach to land use studies in the Suzhou-Wuxi-Changzhou region

In several LUCC studies, statistical methods are being used to analyze land use data. A problem using conventional statistical methods in land use analysis is that these methods assume the data to be statistically independent. But in fact, they have the tendency to be dependent, a phenomenon known as multicollinearity, especially in the cases of few observations. In this paper, a Partial Least-Squares （PLS） regression approach is developed to study relationships between land use and its influencing factors through a case study of the Suzhou-Wuxi-Changzhou region in China. Multicollinearity exists in the dataset and the number of variables is high compared to the number of observations. Four PLS factors are selected through a preliminary analysis. The correlation analyses between land use and influencing factors demonstrate the land use character of rural industrialization and urbanization in the Suzhou-Wuxi-Changzhou region, meanwhile illustrate that the first PLS factor has enough ability to best describe land use patterns quantitatively, and most of the statistical relations derived from it accord with the fact. By the decreasing capacity of the PLS factors, the reliability of model outcome decreases correspondingly.

Optimization of Random Feature Method in the High-Precision Regime

Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.

Monte Carlo method for evaluation of surface emission rate measurement uncertainty

The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.

Reduced-order method for nuclear reactor primary circuit calculation

Accurate real-time simulations of nuclear reactor circuit systems are particularly important for system safety analysis and design.To effectively improve computational efficiency without reducing accuracy,this study establishes a thermal-hydraulics reduced-order model(ROM)for nuclear reactor circuit systems.The full-order circuit system calculation model is first established and verified and then used to calculate the thermal-hydraulic properties of the circuit system under different states as snapshots.The proper orthogonal decomposition method is used to extract the basis functions from snapshots,and the ROM is constructed using the least-squares method,effectively reducing the difficulty in constructing the ROM.A comparison between the full-order simulation and ROM prediction results of the AP1000 circuit system shows that the proposed ROM can improve computational efficiency by 1500 times while achieving a maximum relative error of 0.223%.This research develops a new direction and perspective for the digital twin modeling of nuclear reactor system circuits.

THE INTERIOR TRANSMISSION EIGENVALUE PROBLEM FOR AN ANISOTROPIC MEDIUM BY A PARTIALLY COATED BOUNDARY

We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.

Factorized Smith Method for A Class of High-Ranked Large-Scale T-Stein Equations

We introduce a factorized Smith method(FSM)for solving large-scale highranked T-Stein equations within the banded-plus-low-rank structure framework.To effectively reduce both computational complexity and storage requirements,we develop techniques including deflation and shift,partial truncation and compression,as well as redesign the residual computation and termination condition.Numerical examples demonstrate that the FSM outperforms the Smith method implemented with a hierarchical HODLR structured toolkit in terms of CPU time.

An Adaptive Least-Squares Mixed Finite Element Method for Fourth Order Parabolic Problems

A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approximate. The a posteriori error estimator which is needed in the adaptive refinement algorithm is proposed. The local evaluation of the least-squares functional serves as a posteriori error estimator. The posteriori errors are effectively estimated. The convergence of the adaptive least-squares mixed finite element method is proved.

LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR A CLASS OF STOKES EQUATION

A least-squares mixed finite element method was formulated for a class of Stokes equations in two dimensional domains. The steady state and the time-dependent Stokes' equations were considered. For the stationary equation, optimal H-t and L-2-error estimates are derived under the standard regularity assumption on the finite element partition ( the LBB-condition is not required). Far the evolutionary equation, optimal L-2 estimates are derived under the conventional Raviart-Thomas spaces.

PARTIAL LEAST-SQUARES(PLS)REGRESSION AND SPECTROPHOTOMETRY AS APPLIED TO THE ANALYSIS OF MULTICOMPONENT MIXTURES

The UV absorption spectra of o-naphthol,α-naphthylamine,2,7-dihydroxy naphthalene,2,4-dimethoxy ben- zaldehyde and methyl salicylate,overlap severely;therefore it is impossible to determine them in mixtures by traditional spectrophotometric methods.In this paper,the partial least-squares(PLS)regression is applied to the simultaneous determination of these compounds in mixtures by UV spectrophtometry without any pretreatment of the samples.Ten synthetic mixture samples are analyzed by the proposed method.The mean recoveries are 99.4%,996%,100.2%,99.3% and 99.1%,and the relative standard deviations(RSD) are 1.87%,1.98%,1.94%,0.960% and 0.672%,respectively.

Least-Squares Finite Element Method for the Steady Upper-Convected Maxwell Fluid

In this paper, a least-squares finite element method for the upper-convected Maxell (UCM) fluid is proposed. We first linearize the constitutive and momentum equations and then apply a least-squares method to the linearized version of the viscoelastic UCM model. The L2 least-squares functional involves the residuals of each equation multiplied by proper weights. The corresponding homogeneous functional is equivalent to a natural norm. The error estimates of the finite element solution are analyzed when the conforming piecewise polynomial elements are used for the unknowns.

Complex derivatives valuation: applying the Least-Squares Monte Carlo Simulation Method with several polynomial basis

Background:This article investigates the Least-Squares Monte Carlo Method by using different polynomial basis in American Asian Options pricing.The standard approach in the option pricing literature is to choose the basis arbitrarily.By comparing four different polynomial basis we show that the choice of basis interferes in the option's price.Methods:We assess Least-Squares Method performance in pricing four different American Asian Options by using four polynomial basis:Power,Laguerre,Legendre and Hermite A.To every American Asian Option priced,three sets of parameters are used in order to evaluate it properly.Results:We show that the choice of the basis interferes in the option's price by showing that one of them converges to the option's value faster than any other by using fewer simulated paths.In the case of an Amerasian call option,for example,we find that the preferable polynomial basis is Hermite A.For an Amerasian put option,the Power polynomial basis is recommended.Such empirical outcome is theoretically unpredictable,since in principle all basis can be indistinctly used when pricing the derivative.Conclusion:In this article The Least-Squares Monte Carlo Method performance is assessed in pricing four different types of American Asian Options by using four different polynomial basis through three different sets of parameters.Our results suggest that one polynomial basis is best suited to perform the method when pricing an American Asian option.Theoretically all basis can be indistinctly used when pricing the derivative.However,our results does not confirm these.We find that when pricing an American Asian put option,Power A is better than the other basis we have studied here whereas when pricing an American Asian call,Hermite A is better.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

Initial costates derived by near-optimal reference sequence and least-squares method

In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training a Multilayer Perceptron(MLP)for generating reference sequence and Time of Flight(TOF)to the target,and deriving a system of linear algebraic equations for obtaining the initial costates.To overcome the challenge of generating training samples for the MLP,the backward generation method is proposed to obtain five different training databases.The training database and sample form are determined by analyzing the input and output correlation using the Pearson correlation coefficient.The best-performing MLP is obtained by analyzing the training results with various hyper-parameter combinations.A reference sequence starting from the initial states is obtained by integrating forward with the near-optimal control vector from the output of MLP.Finally,a system of linear algebraic equations for estimating the initial costates is derived using the reference sequence and the necessary conditions for optimality.Simulation results demonstrate that the proposed initial costates solver improves the convergence ratio and reduce the function calls of the shooting function.Furthermore,Monte-Carlo simulation illustrates that the initial costates solver is applicable to different initial velocities,demonstrating excellent generalization ability.