Least-squares reverse time migration in visco-acoustic media based on symplectic stereo-modeling method

The authors proposed a symplectic stereo-modeling method(SSM)in the Birkhoffian dynam-ics and apply it to the visco-acoustic least-squares reverse time migration(LSRTM).The SSM adopts ste-reo-modeling operator in space and symplectic Runge-Kutta scheme in time,resulting in great ability in suppressing numerical dispersion and long-time computing.These advantages are further confirmed by numerical dispersion analysis,long-time computation test and computational efficiency comparison.After these theoretical analyses and experiments,acoustic and visco-acoustic LSRTM are tested and compared between SSM method and the conventional symplectic method(CSM)using the fault and marmousi models.Meanwhile,dynamic source encoding and exponential decay moving average gradients method are adopted to reduce the computation cost and improve the convergence rate.The imaging results show that LSRTM based on visco-acoustic wave equations effectively takes into account the influence of viscosity can therefore compensate for the amplitude attenuation.Besides,SSM method not only has high numerical accuracy and computational efficiency,but also performs effectively in LSRTM.

Least-Squares及Galerkin谱元方法求解环形区域内的泊松方程

为研究基于Least-Squares变分及Galerkin变分两种形式的谱元方法的求解特性,推导了极坐标系中采用两种变分方法求解环形区域内Poisson方程时对应的弱解形式,采用Chebyshev多项式构造插值基函数进行空间离散,得到两种谱元方法对应的代数方程组,由此分析了系数矩阵结构的特点。数值计算结果显示:Least-Squares谱元方法为实现方程的降阶而引入新的求解变量,使得代数方程组形式更为复杂,但边界条件的处理比Galerkin谱元方法更为简单;两种谱元方法均能求解极坐标系中的Poisson方程且能获得高精度的数值解,二者绝对误差分布基本一致;固定单元内的插值阶数时,增加单元数可减小数值误差,且表现出代数精度的特点,误差降低速度较慢,而固定单元数时,在一定范围内数值误差随插值阶数的增加而减小的速度更快,表现出谱精度的特点;单元内插值阶数较高时,代数方程组系数矩阵的条件数急剧增多,方程组呈现病态,数值误差增大,这一特点限制了单元内插值阶数的取值。研究内容对深入了解两种谱元方法在极坐标系中求解Poisson方程时的特点、进一步采用相关分裂算法求解实际流动问题具有参考价值。

Regularized least-squares migration of simultaneous-source seismic data with adaptive singular spectrum analysis

Simultaneous-source acquisition has been recog- nized as an economic and efficient acquisition method, but the direct imaging of the simultaneous-source data produces migration artifacts because of the interference of adjacent sources. To overcome this problem, we propose the regularized least-squares reverse time migration method （RLSRTM） using the singular spectrum analysis technique that imposes sparseness constraints on the inverted model. Additionally, the difference spectrum theory of singular values is presented so that RLSRTM can be implemented adaptively to eliminate the migration artifacts. With numerical tests on a fiat layer model and a Marmousi model, we validate the superior imaging quality, efficiency and convergence of RLSRTM compared with LSRTM when dealing with simultaneoussource data, incomplete data and noisy data.

Characterizing and estimating rice brown spot disease severity using stepwise regression,principal component regression and partial least-square regression

Detecting plant health conditions plays a key role in farm pest management and crop protection. In this study, measurement of hyperspectral leaf reflectance in rice crop (Oryzasativa L.) was conducted on groups of healthy and infected leaves by the fungus Bipolaris oryzae (Helminthosporium oryzae Breda. de Hann) through the wavelength range from 350 to 2 500 nm. The percentage of leaf surface lesions was estimated and defined as the disease severity. Statistical methods like multiple stepwise regression, principal component analysis and partial least-square regression were utilized to calculate and estimate the disease severity of rice brown spot at the leaf level. Our results revealed that multiple stepwise linear regressions could efficiently estimate disease severity with three wavebands in seven steps. The root mean square errors (RMSEs) for training (n=210) and testing (n=53) dataset were 6.5% and 5.8%, respectively. Principal component analysis showed that the first principal component could explain approximately 80% of the variance of the original hyperspectral reflectance. The regression model with the first two principal components predicted a disease severity with RMSEs of 16.3% and 13.9% for the training and testing dataset, respec-tively. Partial least-square regression with seven extracted factors could most effectively predict disease severity compared with other statistical methods with RMSEs of 4.1% and 2.0% for the training and testing dataset, respectively. Our research demon-strates that it is feasible to estimate the disease severity of rice brown spot using hyperspectral reflectance data at the leaf level.

Single Point Positioning with Sequential Least-Squares Filter and Estimated Real-Time Stochastic Model

To obtain higher accurate position estimates, the stochastic model is estimated by using residual of observations, hence, the stochastic model describes the noise and bias in measurements more realistically. By using GPS data and broadcast ephemeris, the numerical results indicating the accurate position estimates at sub-meter level are obtainable.

AN ASSESSMENT OF THE MESHLESS WEIGHTED LEAST-SQUARE METHOD

The meshless weighted least-square (MWLS) method was developed based on the weighted least-square method. The method possesses several advantages, such as high accuracy, high stability and high e?ciency. Moreover, the coe?cient matrix obtained is symmetric and semi- positive de?nite. In this paper, the method is further examined critically. The e?ects of several parameters on the results of MWLS are investigated systematically by using a cantilever beam and an in?nite plate with a central circular hole. The numerical results are compared with those obtained by using the collocation-based meshless method (CBMM) and Galerkin-based meshless method (GBMM). The investigated parameters include the type of approximations, the type of weight functions, the number of neighbors of an evaluation point, as well as the manner in which the neighbors of an evaluation point are determined. This study shows that the displacement accuracy and convergence rate obtained by MWLS is comparable to that of the GBMM while the stress accuracy and convergence rate yielded by MWLS is even higher than that of GBMM. Furthermore, MWLS is much more e?cient than GBMM. This study also shows that the instability of CBMM is mainly due to the neglect of the equi- librium residuals at boundary nodes. In MWLS, the residuals of all the governing equations are minimized in a weighted least-square sense.

Mixed-Weights Least-Squares Stable Predictive Control Algorithm with Soft and Hard Constraints

Mixed-weight least-squares (MWLS) predictive control algorithm, compared with quadratic programming (QP) method, has the advantages of reducing the computer burden, quick calculation speed and dealing with the case in which the optimization is infeasible. But it can only deal with soft constraints. In order to deal with hard constraints and guarantee feasibility, an improved algorithm is proposed by recalculating the setpoint according to the hard constraints before calculating the manipulated variable and MWLS algorithm is used to satisfy the requirement of soft constraints for the system with the input constraints and output constraints. The algorithm can not only guarantee stability of the system and zero steady state error, but also satisfy the hard constraints of input and output variables. The simulation results show the improved algorithm is feasible and effective.

Multisource least-squares reverse-time migration with structure-oriented filtering

The technology of simultaneous-source acquisition of seismic data excited by several sources can significantly improve the data collection efficiency. However, direct imaging of simultaneous-source data or blended data may introduce crosstalk noise and affect the imaging quality. To address this problem, we introduce a structure-oriented filtering operator as preconditioner into the multisource least-squares reverse-time migration （LSRTM）. The structure-oriented filtering operator is a nonstationary filter along structural trends that suppresses crosstalk noise while maintaining structural information. The proposed method uses the conjugate-gradient method to minimize the mismatch between predicted and observed data, while effectively attenuating the interference noise caused by exciting several sources simultaneously. Numerical experiments using synthetic data suggest that the proposed method can suppress the crosstalk noise and produce highly accurate images.

Quantitative energy-dispersive X-ray fluorescence analysis for unknown samples using full-spectrum least-squares regression

The full-spectrum least-squares(FSLS) method is introduced to perform quantitative energy-dispersive X-ray fluorescence analysis for unknown solid samples.Based on the conventional least-squares principle, this spectrum evaluation method is able to obtain the background-corrected and interference-free net peaks, which is significant for quantization analyses. A variety of analytical parameters and functions to describe the features of the fluorescence spectra of pure elements are used and established, such as the mass absorption coefficient, the Gi factor, and fundamental fluorescence formulas. The FSLS iterative program was compiled in the C language. The content of each component should reach the convergence criterion at the end of the calculations. After a basic theory analysis and experimental preparation, 13 national standard soil samples were detected using a spectrometer to test the feasibility of using the algorithm. The results show that the calculated contents of Ti, Fe, Ni, Cu, and Zn have the same changing tendency as the corresponding standard content in the 13 reference samples. Accuracies of 0.35% and 14.03% are obtained, respectively, for Fe and Ti, whose standard concentrations are 8.82% and 0.578%, respectively. However, the calculated results of trace elements (only tens of lg/g) deviate from the standard values. This may be because of measurement accuracy and mutual effects between the elements.

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.

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.

Adaptive Control of a Two-Link Robot Using Batch Least-Square Identifier

We design a regulation-triggered adaptive controller for robot manipulators to efficiently estimate unknown parameters and to achieve asymptotic stability in the presence of coupled uncertainties.Robot manipulators are widely used in telemanipulation systems where they are subject to model and environmental uncertainties.Using conventional control algorithms on such systems can cause not only poor control performance,but also expensive computational costs and catastrophic instabilities.Therefore,system uncertainties need to be estimated through designing a computationally efficient adaptive control law.We focus on robot manipulators as an example of a highly nonlinear system.As a case study,a 2-DOF manipulator subject to four parametric uncertainties is investigated.First,the dynamic equations of the manipulator are derived,and the corresponding regressor matrix is constructed for the unknown parameters.For a general nonlinear system,a theorem is presented to guarantee the asymptotic stability of the system and the convergence of parameters'estimations.Finally,simulation results are discussed for a two-link manipulator,and the performance of the proposed scheme is thoroughly evaluated.

Practical constrained least-square algorithm for moving source location using TDOA and FDOA measurements

By utilizing the time difference of arrival （TDOA） and frequency difference of arrival （FDOA） measurements of signals received at a number of receivers, a constrained least-square （CLS） algorithm for estimating the position and velocity of a moving source is proposed. By utilizing the Lagrange multipliers technique, the known relation between the intermediate variables and the source location coordinates could be exploited to constrain the solution. And without requiring apriori knowledge of TDOA and FDOA measurement noises, the proposed algorithm can satisfy the demand of practical applications. Additionally, on basis of con- volute and polynomial rooting operations, the Lagrange multipliers can be obtained efficiently and robustly allowing real-time imple- mentation and global convergence. Simulation results show that the proposed estimator achieves remarkably better performance than the two-step weighted least square （WLS） approach especially for higher measurement noise level.

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 Quadratic Constraint Total Least-squares Algorithm for Hyperbolic Location

A novel algorithm for source location by utilizing the time difference of arrival (TDOA) measurements of a signal received at spatially separated sensors is proposed. The algorithm is based on quadratic constraint total least-squares (QC-TLS) method and gives an explicit solution. The total least-squares method is a generalized data fitting method that is appropriate for cases when the system model contains error or is not known exactly, and quadratic constraint, which could be realized via Lagrange multipliers technique, could constrain the solution to the location equations to improve location accuracy. Comparisons of performance with ordinary least-squares are made, and Monte Carlo simulations are performed. Simulation results indicate that the proposed algorithm has high location accuracy and achieves accuracy close to the Cramer-Rao lower bound (CRLB) near the small TDOA measurement error region.

Two new least-squares mixed finite element procedures for convection-dominated Sobolev equations

Two new convection-dominated are derived under the approximate solutions least-squares mixed finite element procedures are formulated for solving Sobolev equations. Optimal H（div;Ω）×H1（Ω） norms error estimates standard mixed finite spaces. Moreover, these two schemes provide the with first-order and second-order accuracy in time increment, respectively.

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.

Multiparameter least-squares reverse time migration for acoustic–elastic coupling media based on ocean bottom cable data

In marine seismic exploration,ocean bottom cable technology can record multicomponent seismic data for multiparameter inversion and imaging.This study proposes an elastic multiparameter lease-squares reverse time migration based on the ocean bottom cable technology.Herein,the wavefield continuation operators are mixed equations:the acoustic wave equations are used to calculate seismic wave propagation in the seawater medium,whereas in the solid media below the seabed,the wavefields are obtained by P-and S-wave separated vector elastic wave equations.At the seabed interface,acoustic–elastic coupling control equations are used to combine the two types of equations.P-and S-wave separated elastic migration operators,demigration operators,and gradient equations are derived to realize the elastic least-squares reverse time migration based on the P-and S-wave mode separation.The model tests verify that the proposed method can obtain high-quality images in both the P-and S-velocity components.In comparison with the traditional elastic least-squares reverse time migration method,the proposed method can readily suppress imaging crosstalk noise from multiparameter coupling.

Global interpolating meshless shape function based on generalized moving least-square for structural dynamic analysis

A global interpolating meshless shape function based on the generalized moving least-square （GMLS） is formulated by the transformation technique. Both the shape function and its derivatives meet the Kronecker delta function property. With the interpolating GMLS （IGMLS） shape function, an improved element-free Galerkin （EFG） method is proposed for the structural dynamic analysis. Compared with the conven- tional EFG method, the obvious advantage of the proposed method is that the essential boundary conditions including both displacements and derivatives can be imposed by the straightforward way. Meanwhile, it can greatly improve the ill-condition feature of the standard GMLS approximation, and provide good accuracy at low cost. The dynamic analyses of the Euler beam and Kirchhoff plate are performed to demonstrate the feasi- bility and effectiveness of the improved method. The comparison between the numerical results of the conventional method and the improved method shows that the proposed method has better stability, higher accuracy, and less time consumption.