The fast solution of linear equations has always been one of the hot spots in scientific computing.A kind of the diagonal matrix splitting iteration methods are provided,which is different from the classical matrix sp...The fast solution of linear equations has always been one of the hot spots in scientific computing.A kind of the diagonal matrix splitting iteration methods are provided,which is different from the classical matrix splitting methods.Taking the decomposition of the diagonal elements for coefficient matrix as the key point,some new preconditioners are constructed.Taking the tri-diagonal coefficient matrix as an example,the convergence domains and optimal relaxation factor of the new method are analyzed theoretically.The presented new iteration methods are applied to solve linear algebraic equations,even 2D and 3D diffusion problems with the fully implicit discretization.The results of numerical experiments are matched with the theoretical analysis,and show that the iteration numbers are reduced greatly.The superiorities of presented iteration methods exceed some classical iteration methods dramatically.展开更多
An intuitive method for circle fitting is proposed. Assuming an approximate circle(CA,n) for the fitting of some scattered points, it can be imagined that every point would apply a force to CA,n, which all together fo...An intuitive method for circle fitting is proposed. Assuming an approximate circle(CA,n) for the fitting of some scattered points, it can be imagined that every point would apply a force to CA,n, which all together form an overall effect that "draws" CA,n towards best fitting to the group of points. The basic element of the force is called circular attracting factor(CAF) which is defined as a real scalar in a radial direction of CA,n. An iterative algorithm based on this idea is proposed, and the convergence and accuracy are analyzed. The algorithm converges uniformly which is proved by the analysis of Lyapunov function, and the accuracy of the algorithm is in accord with that of geometric least squares of circle fitting. The algorithm is adopted to circle detection in grayscale images, in which the transferring to binary images is not required, and thus the algorithm is less sensitive to lightening and background noise. The main point for the adaption is the calculation of CAF which is extended in radial directions of CA,n for the whole image. All pixels would apply forces to CA,n, and the overall effect of forces would be equivalent to a force from the centroid of pixels to CA,n. The forces from would-be edge pixels would overweigh that from noisy pixels, so the following approximate circle would be of better fitting. To reduce the amount of calculation, pixels are only used in an annular area including the boundary of CA,n just in between for the calculation of CAF. Examples are given, showing the process of circle fitting of scattered points around a circle from an initial assuming circle, comparing the fitting results for scattered points from some related literature, applying the method proposed for circular edge detection in grayscale images with noise, and/or with only partial arc of a circle, and for circle detection in BGA inspection.展开更多
Dualcolor systems were used to reduce the collinearity of multicomponent spectra, which is described by the angles between spectra vectors. Combined with iterative target transformation factor analysis, single rare ea...Dualcolor systems were used to reduce the collinearity of multicomponent spectra, which is described by the angles between spectra vectors. Combined with iterative target transformation factor analysis, single rare earth element was determined in its mixture. The calculated results show that the average angle between rare earth spectra in one color system(trichloroarsenazorare earths, pH 34) is 45, and that in two color systems(trichloroarsenazorare earths, pH 34, 14) is 215. This technique makes it easy to select the real number of the components in mixtures, and the determination results show dualcolor system method is an effective technique in rare earth mixture analysis.展开更多
For the electro-hydraulic servo vibrating system（ESVS） with the characteristics of non-linearity and repeating motion, a novel method, PI-type iterative learning control（ILC）, is proposed on the basis of tradition...For the electro-hydraulic servo vibrating system（ESVS） with the characteristics of non-linearity and repeating motion, a novel method, PI-type iterative learning control（ILC）, is proposed on the basis of traditional PID control. By using memory ability of computer, the method keeps last time＇s tracking error of the system and then applies the error information to the next time＇s control process. At the same time, a forgetting factor and a D-type learning law of feedforward fuzzy-inferring referenced displacement error under the optimal objective are employed to enhance the systemic robustness and tracking accuracy. The results of simulation and test reveal that the algorithm has a trait of high repeating precision, and could restrain the influence of nonlinear factors like leaking, external disturbance, aerated oil, etc. Compared with traditional PID control, it could better meet the requirement of nonlinear electro -hydraulic servo vibrating system.展开更多
The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of...The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of boundary integral equation is proposed in this paper, together with the solution procedures for multiple crack problems in plane elasticity. With the proposed approach, the multiple crack problems can be solved with the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix as that in the numerical Green’s function （NGF） approach but without the trouble to determine the complementary solutions since the standard boundary element discretization on the crack surface is no longer required with the proposed approach. Some numerical examples computing the stress intensity factors are presented and compared with those in literature to show the accuracy and the effectiveness of the proposed approach.展开更多
In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared wi...In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.展开更多
In the digital image correlation research of fatigue crack growth rate,the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor,thereby affecting the life predic...In the digital image correlation research of fatigue crack growth rate,the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor,thereby affecting the life prediction.This paper proposes a Gauss-Newton iteration method for solving the crack tip position.The conventional linear fitting method provides an iterative initial solution for this method,and the preconditioned conjugate gradient method is used to solve the ill-conditioned matrix.A noise-added artificial displacement field is used to verify the feasibility of the method,which shows that all parameters can be solved with satisfactory results.The actual stress intensity factor solution case shows that the stress intensity factor value obtained by the method in this paper is very close to the finite element result,and the relative error between the two is only−0.621%;The Williams coefficient obtained by this method can also better define the contour of the plastic zone at the crack tip,and the maximum relative error with the test plastic zone area is−11.29%.The relative error between the contour of the plastic zone defined by the conventional method and the area of the experimental plastic zone reached a maximum of 26.05%.The crack tip coordinates,stress intensity factors,and plastic zone contour changes in the loading and unloading phases are explored.The results show that the crack tip change during the loading process is faster than the change during the unloading process;the stress intensity factor during the unloading process under the same load condition is larger than that during the loading process;under the same load,the theoretical plastic zone during the unloading process is higher than that during the loading process.展开更多
A SINS/GNSS location method based on factor diagram is proposed to meet the requirement of accurate location of substation construction personnel. In this paper, the inertial autonomous positioning, carrier motion inf...A SINS/GNSS location method based on factor diagram is proposed to meet the requirement of accurate location of substation construction personnel. In this paper, the inertial autonomous positioning, carrier motion information acquisition and satellite positioning technologies are integrated. The factor graph method is adopted to abstract the measurement information received by inertial navigation and satellite into factor nodes, and the state information into variable nodes, so as to construct the SINS/GNSS construction personnel positioning fusion factor graph model. The Gauss-Newton iterative method is used to implement the recursive updating of variable nodes, and the optimal estimate of the location information of the construction personnel is calculated, which realized the high precision location of the construction personnel. The factor graph method is verified by pedestrian navigation data. The results show that the factor graph method can continuously and stably output high-precision positioning results, and realize non-equidistant fusion of SINS and GNSS. The positioning accuracy is better than Kalman filter algorithm, and the horizontal positioning accuracy is less than 1 m. Therefore, the factor graph method proposed can provide accurate location information for substation construction personnel.展开更多
基金The National Natural Science Foundations of China (12202219)the Natural Science Foundations of Ningxia (2024AAC02009, 2023AAC05001)the Ningxia Youth Top Talents Training Project。
文摘The fast solution of linear equations has always been one of the hot spots in scientific computing.A kind of the diagonal matrix splitting iteration methods are provided,which is different from the classical matrix splitting methods.Taking the decomposition of the diagonal elements for coefficient matrix as the key point,some new preconditioners are constructed.Taking the tri-diagonal coefficient matrix as an example,the convergence domains and optimal relaxation factor of the new method are analyzed theoretically.The presented new iteration methods are applied to solve linear algebraic equations,even 2D and 3D diffusion problems with the fully implicit discretization.The results of numerical experiments are matched with the theoretical analysis,and show that the iteration numbers are reduced greatly.The superiorities of presented iteration methods exceed some classical iteration methods dramatically.
基金This paper is partially supported by the Bulgarian Mutistry of Education SciencesTechnologies Project MM-515/97+1 种基金partially supported by the National Natural Science Foundation of Chinaby Natu-ral Science Foundation of Zhejiang Province.
文摘In this paper we consider some parallel iterations for splitting quadratic factors of polynomials and their convergence.
基金Project(2013CB035504) supported by the National Basic Research Program of ChinaProject(2012zzts078) supported by the Fundamental Research Funds for the Central Universities of Central South University,ChinaProject(2009ZX02038) supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China
文摘An intuitive method for circle fitting is proposed. Assuming an approximate circle(CA,n) for the fitting of some scattered points, it can be imagined that every point would apply a force to CA,n, which all together form an overall effect that "draws" CA,n towards best fitting to the group of points. The basic element of the force is called circular attracting factor(CAF) which is defined as a real scalar in a radial direction of CA,n. An iterative algorithm based on this idea is proposed, and the convergence and accuracy are analyzed. The algorithm converges uniformly which is proved by the analysis of Lyapunov function, and the accuracy of the algorithm is in accord with that of geometric least squares of circle fitting. The algorithm is adopted to circle detection in grayscale images, in which the transferring to binary images is not required, and thus the algorithm is less sensitive to lightening and background noise. The main point for the adaption is the calculation of CAF which is extended in radial directions of CA,n for the whole image. All pixels would apply forces to CA,n, and the overall effect of forces would be equivalent to a force from the centroid of pixels to CA,n. The forces from would-be edge pixels would overweigh that from noisy pixels, so the following approximate circle would be of better fitting. To reduce the amount of calculation, pixels are only used in an annular area including the boundary of CA,n just in between for the calculation of CAF. Examples are given, showing the process of circle fitting of scattered points around a circle from an initial assuming circle, comparing the fitting results for scattered points from some related literature, applying the method proposed for circular edge detection in grayscale images with noise, and/or with only partial arc of a circle, and for circle detection in BGA inspection.
文摘Dualcolor systems were used to reduce the collinearity of multicomponent spectra, which is described by the angles between spectra vectors. Combined with iterative target transformation factor analysis, single rare earth element was determined in its mixture. The calculated results show that the average angle between rare earth spectra in one color system(trichloroarsenazorare earths, pH 34) is 45, and that in two color systems(trichloroarsenazorare earths, pH 34, 14) is 215. This technique makes it easy to select the real number of the components in mixtures, and the determination results show dualcolor system method is an effective technique in rare earth mixture analysis.
文摘For the electro-hydraulic servo vibrating system（ESVS） with the characteristics of non-linearity and repeating motion, a novel method, PI-type iterative learning control（ILC）, is proposed on the basis of traditional PID control. By using memory ability of computer, the method keeps last time＇s tracking error of the system and then applies the error information to the next time＇s control process. At the same time, a forgetting factor and a D-type learning law of feedforward fuzzy-inferring referenced displacement error under the optimal objective are employed to enhance the systemic robustness and tracking accuracy. The results of simulation and test reveal that the algorithm has a trait of high repeating precision, and could restrain the influence of nonlinear factors like leaking, external disturbance, aerated oil, etc. Compared with traditional PID control, it could better meet the requirement of nonlinear electro -hydraulic servo vibrating system.
基金supported by the National Natural Science Foundation of China (Grant No.10972131)the Graduate Innovation Foundation of Shanghai University (Grant No.SHUCX102351)
文摘The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of boundary integral equation is proposed in this paper, together with the solution procedures for multiple crack problems in plane elasticity. With the proposed approach, the multiple crack problems can be solved with the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix as that in the numerical Green’s function （NGF） approach but without the trouble to determine the complementary solutions since the standard boundary element discretization on the crack surface is no longer required with the proposed approach. Some numerical examples computing the stress intensity factors are presented and compared with those in literature to show the accuracy and the effectiveness of the proposed approach.
基金Project supported by the National Natural Science Foundation of China （Grant No.10271074）, and the Special Funds for Major Specialities of Shanghai Education Commission （Grant No.J50101）
文摘In this paper, an improved gradient iterative （GI） algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.
基金Supported by National Natural Science Foundation of China(Grant No.51675446)Independent Research Project of State Key Laboratory of Traction Power(Grant No.2019TPL-T13).
文摘In the digital image correlation research of fatigue crack growth rate,the accuracy of the crack tip position determines the accuracy of the calculation of the stress intensity factor,thereby affecting the life prediction.This paper proposes a Gauss-Newton iteration method for solving the crack tip position.The conventional linear fitting method provides an iterative initial solution for this method,and the preconditioned conjugate gradient method is used to solve the ill-conditioned matrix.A noise-added artificial displacement field is used to verify the feasibility of the method,which shows that all parameters can be solved with satisfactory results.The actual stress intensity factor solution case shows that the stress intensity factor value obtained by the method in this paper is very close to the finite element result,and the relative error between the two is only−0.621%;The Williams coefficient obtained by this method can also better define the contour of the plastic zone at the crack tip,and the maximum relative error with the test plastic zone area is−11.29%.The relative error between the contour of the plastic zone defined by the conventional method and the area of the experimental plastic zone reached a maximum of 26.05%.The crack tip coordinates,stress intensity factors,and plastic zone contour changes in the loading and unloading phases are explored.The results show that the crack tip change during the loading process is faster than the change during the unloading process;the stress intensity factor during the unloading process under the same load condition is larger than that during the loading process;under the same load,the theoretical plastic zone during the unloading process is higher than that during the loading process.
文摘A SINS/GNSS location method based on factor diagram is proposed to meet the requirement of accurate location of substation construction personnel. In this paper, the inertial autonomous positioning, carrier motion information acquisition and satellite positioning technologies are integrated. The factor graph method is adopted to abstract the measurement information received by inertial navigation and satellite into factor nodes, and the state information into variable nodes, so as to construct the SINS/GNSS construction personnel positioning fusion factor graph model. The Gauss-Newton iterative method is used to implement the recursive updating of variable nodes, and the optimal estimate of the location information of the construction personnel is calculated, which realized the high precision location of the construction personnel. The factor graph method is verified by pedestrian navigation data. The results show that the factor graph method can continuously and stably output high-precision positioning results, and realize non-equidistant fusion of SINS and GNSS. The positioning accuracy is better than Kalman filter algorithm, and the horizontal positioning accuracy is less than 1 m. Therefore, the factor graph method proposed can provide accurate location information for substation construction personnel.
