Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with rand...Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with random errors.However,in many geodetic applications,some elements are error-free and some random observations appear repeatedly in different positions in the augmented coefficient matrix.It is called the linear structured EIV(LSEIV)model.Two kinds of methods are proposed for the LSEIV model from functional and stochastic modifications.On the one hand,the functional part of the LSEIV model is modified into the errors-in-observations(EIO)model.On the other hand,the stochastic model is modified by applying the Moore-Penrose inverse of the cofactor matrix.The algorithms are derived through the Lagrange multipliers method and linear approximation.The estimation principles and iterative formula of the parameters are proven to be consistent.The first-order approximate variance-covariance matrix(VCM)of the parameters is also derived.A numerical example is given to compare the performances of our proposed three algorithms with the STLS approach.Afterwards,the least squares(LS),total least squares(TLS)and linear structured weighted total least squares(LSWTLS)solutions are compared and the accuracy evaluation formula is proven to be feasible and effective.Finally,the LSWTLS is applied to the field of deformation analysis,which yields a better result than the traditional LS and TLS estimations.展开更多
To improve the self-cleaning ability of aquaculture tank and the efficiency of circulating water,physical and numerical experiments were conducted on the influence of inlet structure on sewage discharge in a rounded s...To improve the self-cleaning ability of aquaculture tank and the efficiency of circulating water,physical and numerical experiments were conducted on the influence of inlet structure on sewage discharge in a rounded square aquaculture tank with a single inlet.Based on the physical model of the tank,analysis of how inlet structure adjustment affects sewage discharge efficiency and flow field characteristics was conducted to provide suitable flow field conditions for sinkable solid particle discharge.In addition,an internal flow field simulation was conducted using the RNG k-εturbulence model in hydraulic drive mode.Then a solid-fluid multiphase model was created to investigate how the inlet structure affects sewage collection in the rounded square aquaculture tank with single inlet and outlet.The finding revealed that the impact of inlet structure is considerably affecting sewage collection.The conditions of C/B=0.07-0.11(the ratio of horizontal distance between the center of the inlet pipe and the tank wall(C)to length of the tank(B))andα=25°(αis the angle between the direction of the jet and the tangential direction of the arc angle)resulted in optimal sewage collection,which is similar to the flow field experiment in the rounded square aquaculture tank with single inlet and outlet.An excellent correlation was revealed between sewage collection and fluid circulation stability in the aquaculture tank.The present study provided a reference for design and optimization of circulating aquaculture tanks in aquaculture industry.展开更多
Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted...Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.展开更多
One-class classification problem has become a popular problem in many fields, with a wide range of applications in anomaly detection, fault diagnosis, and face recognition. We investigate the one-class classification ...One-class classification problem has become a popular problem in many fields, with a wide range of applications in anomaly detection, fault diagnosis, and face recognition. We investigate the one-class classification problem for second-order tensor data. Traditional vector-based one-class classification methods such as one-class support vector machine (OCSVM) and least squares one-class support vector machine (LSOCSVM) have limitations when tensor is used as input data, so we propose a new tensor one-class classification method, LSOCSTM, which directly uses tensor as input data. On one hand, using tensor as input data not only enables to classify tensor data, but also for vector data, classifying it after high dimensionalizing it into tensor still improves the classification accuracy and overcomes the over-fitting problem. On the other hand, different from one-class support tensor machine (OCSTM), we use squared loss instead of the original loss function so that we solve a series of linear equations instead of quadratic programming problems. Therefore, we use the distance to the hyperplane as a metric for classification, and the proposed method is more accurate and faster compared to existing methods. The experimental results show the high efficiency of the proposed method compared with several state-of-the-art methods.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
This article explores the comparison between the probability method and the least squares method in the design of linear predictive models. It points out that these two approaches have distinct theoretical foundations...This article explores the comparison between the probability method and the least squares method in the design of linear predictive models. It points out that these two approaches have distinct theoretical foundations and can lead to varied or similar results in terms of precision and performance under certain assumptions. The article underlines the importance of comparing these two approaches to choose the one best suited to the context, available data and modeling objectives.展开更多
The primary aim of the power system grounding is to safeguard the person and satisfying the performance of the power systemtomaintain reliable operation.With equal conductor spacing grounding grid design,the distribut...The primary aim of the power system grounding is to safeguard the person and satisfying the performance of the power systemtomaintain reliable operation.With equal conductor spacing grounding grid design,the distribution of the current in the grid is not uniform.Hence,unequal grid conductor span in which grid conductors are concentrated more at the periphery is safer to practice than equal spacing.This paper presents the comparative analysis of two novel techniques that create unequal spacing among the grid conductors:the least-square curve fitting technique and the compression ratio techniquewith equal grid configuration for both square and rectangular grids.Particle Swarm Optimization(PSO)is adopted for finding out one optimal feasible solution among many feasible solutions of equal grid configuration for both square and rectangular grids.Comparative analysis is also carried out between square and rectangular grids using the least square curve fitting technique as it results in only one unequal grid configuration.Simulation results are obtained by theMATLAB software developed.Percentage of improvement in ground potential rise,step voltage,touch voltage,and grid resistancewith variation in compression ratios are plotted.展开更多
Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, ani...Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, anisotropic rectangular plate subjected to arbitrary boundary conditions on one edge and simply supported on other edge. The rectangular laminate consists of anisotropic and homogeneous laminae of arbitrary thicknesses. This study presents the elastic analysis of laminated composite plates subjected to sinusoidal mechanical loading under arbitrary boundary conditions. Least square finite element solutions for displacements and stresses are investigated using a mathematical model, called a state-space model, which allows us to simultaneously solve for these field variables in the composite structure’s domain and ensure that continuity conditions are satisfied at layer interfaces. The governing equations are derived from this model using a numerical technique called the least-squares finite element method (LSFEM). These LSFEMs seek to minimize the squares of the governing equations and the associated side conditions residuals over the computational domain. The model is comprised of layerwise variables such as displacements, out-of-plane stresses, and in- plane strains, treated as independent variables. Numerical results are presented to demonstrate the response of the laminated composite plates under various arbitrary boundary conditions using LSFEM and compared with the 3D elasticity solution available in the literature.展开更多
An S-N curve fitting approach is proposed based on the weighted least square method, and the weights are inversely proportional to the length of mean confidence intervals of experimental data sets. The assumption coin...An S-N curve fitting approach is proposed based on the weighted least square method, and the weights are inversely proportional to the length of mean confidence intervals of experimental data sets. The assumption coincides with the physical characteristics of the fatigue life scatter. Two examples demonstrate the method. It is shown that the method has better accuracy and reasonableness compared with the usual least square method.展开更多
Least squares projection twin support vector machine(LSPTSVM)has faster computing speed than classical least squares support vector machine(LSSVM).However,LSPTSVM is sensitive to outliers and its solution lacks sparsi...Least squares projection twin support vector machine(LSPTSVM)has faster computing speed than classical least squares support vector machine(LSSVM).However,LSPTSVM is sensitive to outliers and its solution lacks sparsity.Therefore,it is difficult for LSPTSVM to process large-scale datasets with outliers.In this paper,we propose a robust LSPTSVM model(called R-LSPTSVM)by applying truncated least squares loss function.The robustness of R-LSPTSVM is proved from a weighted perspective.Furthermore,we obtain the sparse solution of R-LSPTSVM by using the pivoting Cholesky factorization method in primal space.Finally,the sparse R-LSPTSVM algorithm(SR-LSPTSVM)is proposed.Experimental results show that SR-LSPTSVM is insensitive to outliers and can deal with large-scale datasets fastly.展开更多
Improving the focusing capability of pre-stack time migration allows the imaged section to reflect structural characteristics, depth, and interface shape and it is a key step for the preparation of the initial depth m...Improving the focusing capability of pre-stack time migration allows the imaged section to reflect structural characteristics, depth, and interface shape and it is a key step for the preparation of the initial depth migration velocity model. The traditional symmetrical travel time equation is derived based on the assumption of a layered model. It is difficult to achieve the desired effect of focusing in media with strong lateral variation. The nonsymmetrical travel time equation based on Lie algebra and a pseudo-differential operator contains a lateral velocity derivative which can improve the focusing capability even in strongly lateral variable media and also the computation precision of the weight coefficients for relative amplitude preservation. Compared with the symmetrical methods, the nonsymmetrical method is more effective. In this paper, we describe several key steps of nonsymmetric pre-stack travel time calculation and present some test results using synthetic and real data.展开更多
In regression, despite being both aimed at estimating the Mean Squared Prediction Error (MSPE), Akaike’s Final Prediction Error (FPE) and the Generalized Cross Validation (GCV) selection criteria are usually derived ...In regression, despite being both aimed at estimating the Mean Squared Prediction Error (MSPE), Akaike’s Final Prediction Error (FPE) and the Generalized Cross Validation (GCV) selection criteria are usually derived from two quite different perspectives. Here, settling on the most commonly accepted definition of the MSPE as the expectation of the squared prediction error loss, we provide theoretical expressions for it, valid for any linear model (LM) fitter, be it under random or non random designs. Specializing these MSPE expressions for each of them, we are able to derive closed formulas of the MSPE for some of the most popular LM fitters: Ordinary Least Squares (OLS), with or without a full column rank design matrix;Ordinary and Generalized Ridge regression, the latter embedding smoothing splines fitting. For each of these LM fitters, we then deduce a computable estimate of the MSPE which turns out to coincide with Akaike’s FPE. Using a slight variation, we similarly get a class of MSPE estimates coinciding with the classical GCV formula for those same LM fitters.展开更多
Scour around a submerged square pile was realized experimentally in a steady flow to study the effects of flow depth on local scour.Flow depth to pile height ratios ranging from 1.5 to 5 in uniform sand and 2 to 5 in ...Scour around a submerged square pile was realized experimentally in a steady flow to study the effects of flow depth on local scour.Flow depth to pile height ratios ranging from 1.5 to 5 in uniform sand and 2 to 5 in non-uniform sand were tested in the approaching flow velocity to critical velocity(larger than which the sediment particle is motivated)ratios of 0.56 and 1.03,respectively.The influences of flow depth were investigated on the basis of analysis of the three-dimensional topography,temporal maximum scour depth,bed profile development,and equilibrium scour depth.Results showed that the maximum scour depth was at the upstream corners of the pile other than at the stagnation point.The evolutions of the maximum scour depth data in non-uniform sand were well fitted with a recent exponential function,which characterized the initial,developing,and equilibrium stages of scour depth.The scour hole slopes upstream of the pile were found to be parallel to each other in the process of each test and were mainly governed by the sediment repose underwater.The equilibrium scour depth varied slightly with flow depth when the submergence ratio was larger than 1 in uniform sand while it was 2 in non-uniform sand.The armoring effects of coarse sediment particles markedly reduced the sediment transport in non-uniform sand despite the 0.34 increment in non-uniformity.展开更多
For the accurate extraction of cavity decay time, a selection of data points is supplemented to the weighted least square method. We derive the expected precision, accuracy and computation cost of this improved method...For the accurate extraction of cavity decay time, a selection of data points is supplemented to the weighted least square method. We derive the expected precision, accuracy and computation cost of this improved method, and examine these performances by simulation. By comparing this method with the nonlinear least square fitting (NLSF) method and the linear regression of the sum (LRS) method in derivations and simulations, we find that this method can achieve the same or even better precision, comparable accuracy, and lower computation cost. We test this method by experimental decay signals. The results are in agreement with the ones obtained from the nonlinear least square fitting method.展开更多
In factor analysis, a factor loading matrix is often rotated to a simple target matrix for its simplicity. For the purpose, Procrustes rotation minimizes the discrepancy between the target and rotated loadings using t...In factor analysis, a factor loading matrix is often rotated to a simple target matrix for its simplicity. For the purpose, Procrustes rotation minimizes the discrepancy between the target and rotated loadings using two types of approximation: 1) approximate the zeros in the target by the non-zeros in the loadings, and 2) approximate the non-zeros in the target by the non-zeros in the loadings. The central issue of Procrustes rotation considered in the article is that it equally treats the two types of approximation, while the former is more important for simplifying the loading matrix. Furthermore, a well-known issue of Simplimax is the computational inefficiency in estimating the sparse target matrix, which yields a considerable number of local minima. The research proposes a new rotation procedure that consists of the following two stages. The first stage estimates sparse target matrix with lesser computational cost by regularization technique. In the second stage, a loading matrix is rotated to the target, emphasizing on the approximation of non-zeros to zeros in the target by least squares criterion with generalized weighing that is newly proposed by the study. The simulation study and real data examples revealed that the proposed method surely simplifies loading matrices.展开更多
A two-dimensional numerical study of laminar natural convection in a square enclosure filled with air with a wall partially heated on the bottom is presented.The heat source is located on the lower wall with different...A two-dimensional numerical study of laminar natural convection in a square enclosure filled with air with a wall partially heated on the bottom is presented.The heat source is located on the lower wall with different heated widths varied from 20 to 80%(ε=0.2–0.8)of the total width of the lower wall and different heights h=H/4 and H/2 of the partition.The effect of the partition height on the main system dynamics is investigated through solution of the two-dimensional Navier-Stokes equations and the energy equation by means of a finite volume method based on the SIMPLE algorithm.The influence of the Rayleigh number(Ra=10^(3) to 10^(6))and the hot wall length is also examined.It is shown that the average Nusselt number grows whenεincreases and when h decreases.For a given value ofεand h,the average Nusselt number increases as Ra increases.It is concluded that the partition height causes a decrease in the average Nusselt number.展开更多
This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations in...This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.展开更多
Let be an injective function. For a vertex labeling f, the induced edge labeling is defined by, or;then, the edge labels are distinct and are from . Then f is called a root square mean labeling of G. In this paper, we...Let be an injective function. For a vertex labeling f, the induced edge labeling is defined by, or;then, the edge labels are distinct and are from . Then f is called a root square mean labeling of G. In this paper, we prove root square mean labeling of some degree splitting graphs.展开更多
基金the financial support of the National Natural Science Foundation of China(Grant No.42074016,42104025,42274057and 41704007)Hunan Provincial Natural Science Foundation of China(Grant No.2021JJ30244)Scientific Research Fund of Hunan Provincial Education Department(Grant No.22B0496)。
文摘Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with random errors.However,in many geodetic applications,some elements are error-free and some random observations appear repeatedly in different positions in the augmented coefficient matrix.It is called the linear structured EIV(LSEIV)model.Two kinds of methods are proposed for the LSEIV model from functional and stochastic modifications.On the one hand,the functional part of the LSEIV model is modified into the errors-in-observations(EIO)model.On the other hand,the stochastic model is modified by applying the Moore-Penrose inverse of the cofactor matrix.The algorithms are derived through the Lagrange multipliers method and linear approximation.The estimation principles and iterative formula of the parameters are proven to be consistent.The first-order approximate variance-covariance matrix(VCM)of the parameters is also derived.A numerical example is given to compare the performances of our proposed three algorithms with the STLS approach.Afterwards,the least squares(LS),total least squares(TLS)and linear structured weighted total least squares(LSWTLS)solutions are compared and the accuracy evaluation formula is proven to be feasible and effective.Finally,the LSWTLS is applied to the field of deformation analysis,which yields a better result than the traditional LS and TLS estimations.
基金Supported by the 2023 Central Government Finance Subsidy Project for Liaoning Fisheries,the Key Research Project of Liaoning Provincial Department of Education in 2022(No.LJKZZ20220091)the National Natural Science Foundation of China(No.31872609)+1 种基金the Innovation Support Program for High-level Talents of Dalian City(No.2019RD12)the earmarked fund for CARS-49。
文摘To improve the self-cleaning ability of aquaculture tank and the efficiency of circulating water,physical and numerical experiments were conducted on the influence of inlet structure on sewage discharge in a rounded square aquaculture tank with a single inlet.Based on the physical model of the tank,analysis of how inlet structure adjustment affects sewage discharge efficiency and flow field characteristics was conducted to provide suitable flow field conditions for sinkable solid particle discharge.In addition,an internal flow field simulation was conducted using the RNG k-εturbulence model in hydraulic drive mode.Then a solid-fluid multiphase model was created to investigate how the inlet structure affects sewage collection in the rounded square aquaculture tank with single inlet and outlet.The finding revealed that the impact of inlet structure is considerably affecting sewage collection.The conditions of C/B=0.07-0.11(the ratio of horizontal distance between the center of the inlet pipe and the tank wall(C)to length of the tank(B))andα=25°(αis the angle between the direction of the jet and the tangential direction of the arc angle)resulted in optimal sewage collection,which is similar to the flow field experiment in the rounded square aquaculture tank with single inlet and outlet.An excellent correlation was revealed between sewage collection and fluid circulation stability in the aquaculture tank.The present study provided a reference for design and optimization of circulating aquaculture tanks in aquaculture industry.
基金supported by the Key Project of Gansu Provincial National Science Foundation(23JRRA1022)the National Natural Science Foundation of China(12071431)+1 种基金the Fundamental Research Funds for the Central Universities(lzujbky-2021-ey18)the Innovative Groups of Basic Research in Gansu Province(22JR5RA391).
文摘Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.
文摘One-class classification problem has become a popular problem in many fields, with a wide range of applications in anomaly detection, fault diagnosis, and face recognition. We investigate the one-class classification problem for second-order tensor data. Traditional vector-based one-class classification methods such as one-class support vector machine (OCSVM) and least squares one-class support vector machine (LSOCSVM) have limitations when tensor is used as input data, so we propose a new tensor one-class classification method, LSOCSTM, which directly uses tensor as input data. On one hand, using tensor as input data not only enables to classify tensor data, but also for vector data, classifying it after high dimensionalizing it into tensor still improves the classification accuracy and overcomes the over-fitting problem. On the other hand, different from one-class support tensor machine (OCSTM), we use squared loss instead of the original loss function so that we solve a series of linear equations instead of quadratic programming problems. Therefore, we use the distance to the hyperplane as a metric for classification, and the proposed method is more accurate and faster compared to existing methods. The experimental results show the high efficiency of the proposed method compared with several state-of-the-art methods.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
文摘This article explores the comparison between the probability method and the least squares method in the design of linear predictive models. It points out that these two approaches have distinct theoretical foundations and can lead to varied or similar results in terms of precision and performance under certain assumptions. The article underlines the importance of comparing these two approaches to choose the one best suited to the context, available data and modeling objectives.
文摘The primary aim of the power system grounding is to safeguard the person and satisfying the performance of the power systemtomaintain reliable operation.With equal conductor spacing grounding grid design,the distribution of the current in the grid is not uniform.Hence,unequal grid conductor span in which grid conductors are concentrated more at the periphery is safer to practice than equal spacing.This paper presents the comparative analysis of two novel techniques that create unequal spacing among the grid conductors:the least-square curve fitting technique and the compression ratio techniquewith equal grid configuration for both square and rectangular grids.Particle Swarm Optimization(PSO)is adopted for finding out one optimal feasible solution among many feasible solutions of equal grid configuration for both square and rectangular grids.Comparative analysis is also carried out between square and rectangular grids using the least square curve fitting technique as it results in only one unequal grid configuration.Simulation results are obtained by theMATLAB software developed.Percentage of improvement in ground potential rise,step voltage,touch voltage,and grid resistancewith variation in compression ratios are plotted.
文摘Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, anisotropic rectangular plate subjected to arbitrary boundary conditions on one edge and simply supported on other edge. The rectangular laminate consists of anisotropic and homogeneous laminae of arbitrary thicknesses. This study presents the elastic analysis of laminated composite plates subjected to sinusoidal mechanical loading under arbitrary boundary conditions. Least square finite element solutions for displacements and stresses are investigated using a mathematical model, called a state-space model, which allows us to simultaneously solve for these field variables in the composite structure’s domain and ensure that continuity conditions are satisfied at layer interfaces. The governing equations are derived from this model using a numerical technique called the least-squares finite element method (LSFEM). These LSFEMs seek to minimize the squares of the governing equations and the associated side conditions residuals over the computational domain. The model is comprised of layerwise variables such as displacements, out-of-plane stresses, and in- plane strains, treated as independent variables. Numerical results are presented to demonstrate the response of the laminated composite plates under various arbitrary boundary conditions using LSFEM and compared with the 3D elasticity solution available in the literature.
文摘An S-N curve fitting approach is proposed based on the weighted least square method, and the weights are inversely proportional to the length of mean confidence intervals of experimental data sets. The assumption coincides with the physical characteristics of the fatigue life scatter. Two examples demonstrate the method. It is shown that the method has better accuracy and reasonableness compared with the usual least square method.
基金supported by the National Natural Science Foundation of China(6177202062202433+4 种基金621723716227242262036010)the Natural Science Foundation of Henan Province(22100002)the Postdoctoral Research Grant in Henan Province(202103111)。
文摘Least squares projection twin support vector machine(LSPTSVM)has faster computing speed than classical least squares support vector machine(LSSVM).However,LSPTSVM is sensitive to outliers and its solution lacks sparsity.Therefore,it is difficult for LSPTSVM to process large-scale datasets with outliers.In this paper,we propose a robust LSPTSVM model(called R-LSPTSVM)by applying truncated least squares loss function.The robustness of R-LSPTSVM is proved from a weighted perspective.Furthermore,we obtain the sparse solution of R-LSPTSVM by using the pivoting Cholesky factorization method in primal space.Finally,the sparse R-LSPTSVM algorithm(SR-LSPTSVM)is proposed.Experimental results show that SR-LSPTSVM is insensitive to outliers and can deal with large-scale datasets fastly.
基金This research was supported by the National Basic Research Program of China (Grant No. 2007CB209603), Key Project of the National Natural Science Foundation (Grant No. 40830424), State Key Laboratory of Geological Processes and Mineral Resources Geo-detection Laboratory of the Ministry of Education for their sponsorship (GPMR 200633, GDL0801).
文摘Improving the focusing capability of pre-stack time migration allows the imaged section to reflect structural characteristics, depth, and interface shape and it is a key step for the preparation of the initial depth migration velocity model. The traditional symmetrical travel time equation is derived based on the assumption of a layered model. It is difficult to achieve the desired effect of focusing in media with strong lateral variation. The nonsymmetrical travel time equation based on Lie algebra and a pseudo-differential operator contains a lateral velocity derivative which can improve the focusing capability even in strongly lateral variable media and also the computation precision of the weight coefficients for relative amplitude preservation. Compared with the symmetrical methods, the nonsymmetrical method is more effective. In this paper, we describe several key steps of nonsymmetric pre-stack travel time calculation and present some test results using synthetic and real data.
文摘In regression, despite being both aimed at estimating the Mean Squared Prediction Error (MSPE), Akaike’s Final Prediction Error (FPE) and the Generalized Cross Validation (GCV) selection criteria are usually derived from two quite different perspectives. Here, settling on the most commonly accepted definition of the MSPE as the expectation of the squared prediction error loss, we provide theoretical expressions for it, valid for any linear model (LM) fitter, be it under random or non random designs. Specializing these MSPE expressions for each of them, we are able to derive closed formulas of the MSPE for some of the most popular LM fitters: Ordinary Least Squares (OLS), with or without a full column rank design matrix;Ordinary and Generalized Ridge regression, the latter embedding smoothing splines fitting. For each of these LM fitters, we then deduce a computable estimate of the MSPE which turns out to coincide with Akaike’s FPE. Using a slight variation, we similarly get a class of MSPE estimates coinciding with the classical GCV formula for those same LM fitters.
基金the support of the National Natural Science Foundation of China(Nos.51679223 and 51739010)the 111 Project(No.B14028),the Shangdong Provincial Key Laboratory of Ocean Engineering(No.kl oe202009)+1 种基金the Ningbo Natural Science Foundation(No.2021J096)a grant from the 7th Generation Ultra-Deepwater Drilling Rig Innovation Project。
文摘Scour around a submerged square pile was realized experimentally in a steady flow to study the effects of flow depth on local scour.Flow depth to pile height ratios ranging from 1.5 to 5 in uniform sand and 2 to 5 in non-uniform sand were tested in the approaching flow velocity to critical velocity(larger than which the sediment particle is motivated)ratios of 0.56 and 1.03,respectively.The influences of flow depth were investigated on the basis of analysis of the three-dimensional topography,temporal maximum scour depth,bed profile development,and equilibrium scour depth.Results showed that the maximum scour depth was at the upstream corners of the pile other than at the stagnation point.The evolutions of the maximum scour depth data in non-uniform sand were well fitted with a recent exponential function,which characterized the initial,developing,and equilibrium stages of scour depth.The scour hole slopes upstream of the pile were found to be parallel to each other in the process of each test and were mainly governed by the sediment repose underwater.The equilibrium scour depth varied slightly with flow depth when the submergence ratio was larger than 1 in uniform sand while it was 2 in non-uniform sand.The armoring effects of coarse sediment particles markedly reduced the sediment transport in non-uniform sand despite the 0.34 increment in non-uniformity.
基金supported by the Preeminent Youth Fund of Sichuan Province,China(Grant No.2012JQ0012)the National Natural Science Foundation of China(Grant Nos.11173008,10974202,and 60978049)the National Key Scientific and Research Equipment Development Project of China(Grant No.ZDYZ2013-2)
文摘For the accurate extraction of cavity decay time, a selection of data points is supplemented to the weighted least square method. We derive the expected precision, accuracy and computation cost of this improved method, and examine these performances by simulation. By comparing this method with the nonlinear least square fitting (NLSF) method and the linear regression of the sum (LRS) method in derivations and simulations, we find that this method can achieve the same or even better precision, comparable accuracy, and lower computation cost. We test this method by experimental decay signals. The results are in agreement with the ones obtained from the nonlinear least square fitting method.
文摘In factor analysis, a factor loading matrix is often rotated to a simple target matrix for its simplicity. For the purpose, Procrustes rotation minimizes the discrepancy between the target and rotated loadings using two types of approximation: 1) approximate the zeros in the target by the non-zeros in the loadings, and 2) approximate the non-zeros in the target by the non-zeros in the loadings. The central issue of Procrustes rotation considered in the article is that it equally treats the two types of approximation, while the former is more important for simplifying the loading matrix. Furthermore, a well-known issue of Simplimax is the computational inefficiency in estimating the sparse target matrix, which yields a considerable number of local minima. The research proposes a new rotation procedure that consists of the following two stages. The first stage estimates sparse target matrix with lesser computational cost by regularization technique. In the second stage, a loading matrix is rotated to the target, emphasizing on the approximation of non-zeros to zeros in the target by least squares criterion with generalized weighing that is newly proposed by the study. The simulation study and real data examples revealed that the proposed method surely simplifies loading matrices.
文摘A two-dimensional numerical study of laminar natural convection in a square enclosure filled with air with a wall partially heated on the bottom is presented.The heat source is located on the lower wall with different heated widths varied from 20 to 80%(ε=0.2–0.8)of the total width of the lower wall and different heights h=H/4 and H/2 of the partition.The effect of the partition height on the main system dynamics is investigated through solution of the two-dimensional Navier-Stokes equations and the energy equation by means of a finite volume method based on the SIMPLE algorithm.The influence of the Rayleigh number(Ra=10^(3) to 10^(6))and the hot wall length is also examined.It is shown that the average Nusselt number grows whenεincreases and when h decreases.For a given value ofεand h,the average Nusselt number increases as Ra increases.It is concluded that the partition height causes a decrease in the average Nusselt number.
基金Supported by the NSF of Hubei Province(2022CFD042)。
文摘This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.
文摘Let be an injective function. For a vertex labeling f, the induced edge labeling is defined by, or;then, the edge labels are distinct and are from . Then f is called a root square mean labeling of G. In this paper, we prove root square mean labeling of some degree splitting graphs.