With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel al...With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel algorithms is a difficult issue.We introduce the idea of group search to the single-point search-based non-linear inversion algorithm, taking the quantum Monte Carlo method as an example for two-dimensional seismic wave velocity inversion and practical impedance inversion and test the calculation efficiency of using different node numbers.The results show the parallel algorithm in theoretical and practical data inversion is feasible and effective.The parallel algorithm has good versatility. The algorithm efficiency increases with increasing node numbers but the algorithm efficiency rate of increase gradually decreases as the node numbers increase.展开更多
A new limited memory symmetric rank one algorithm is proposed. It combines a modified self-scaled symmetric rank one (SSR1) update with the limited memory and nonmonotone line search technique. In this algorithm, th...A new limited memory symmetric rank one algorithm is proposed. It combines a modified self-scaled symmetric rank one (SSR1) update with the limited memory and nonmonotone line search technique. In this algorithm, the descent search direction is generated by inverse limited memory SSR1 update, thus simplifying the computation. Numerical comparison of the algorithm and the famous limited memory BFGS algorithm is given. Comparison results indicate that the new algorithm can process a kind of large-scale unconstrained optimization problems.展开更多
In order to study the sliding characteristics when the cable structures are connected with other rods, a string of sliding cable dements (SCE) consisting of one active threenode SCE passing through the sliding point...In order to study the sliding characteristics when the cable structures are connected with other rods, a string of sliding cable dements (SCE) consisting of one active threenode SCE passing through the sliding point and multiple inactive two-node SCEs is put forward. Based on the updated Lagrangian formulation, the geometric nonlinear stiffness matrix of the three-node straight sliding cable dement is deduced. The examples about two-span and three-span continuous cable structures are studied to verify the effectiveness of the derived SCE. Comparing the cable tension of SCE with the existing research results, the calculating results show that the error is less than 1%. The sliding characteristics should be considered in practical engineering because of the obvious difference between the cable tension of the SCE and that of the cable element without considering sliding characteristics.展开更多
A broadband microstrip patch antenna was analyzed and designed.Full wave analysis method(FWAM) was employed to show that a stacked microstrip dual patch antenna(SMDPA) might have a much wider bandwidth than that of ...A broadband microstrip patch antenna was analyzed and designed.Full wave analysis method(FWAM) was employed to show that a stacked microstrip dual patch antenna(SMDPA) might have a much wider bandwidth than that of the ordinanry uni patch one.By means of discrete complex image theory(DCIT),the Sommerfeld integrals (SI) involved were accurately calculated at a speed several hundred times faster than numerical integration method(NIM).The feeding structure of the SMDPA was then improved and the bandwidth was extended to about 22% or more for voltage standing wave ratio (VSWR)s≤2 Finally,a matching network was constructed to obtain a bandwidth of about 25% for s≤1.5.展开更多
Based on the analysis to the random sear ch algorithm of LUUS, a modified random directed integer search algorithm (MRDI SA) is given for first time. And a practical example is given to show that the adva ntage of th...Based on the analysis to the random sear ch algorithm of LUUS, a modified random directed integer search algorithm (MRDI SA) is given for first time. And a practical example is given to show that the adva ntage of this kind of algorithm is the reliability can’t be infuenced by the ini tial value X (0) and the start search domain R (0) . Besides, i t can be applied to solve the higher dimensional constrained nonlinear integer p rogramming problem.展开更多
The maximal entropy ordered weighted averaging (ME-OWA) operator is used to aggregate metasearch engine results, and its newly analytical solution is also applied. Within the current context of the OWA operator, the...The maximal entropy ordered weighted averaging (ME-OWA) operator is used to aggregate metasearch engine results, and its newly analytical solution is also applied. Within the current context of the OWA operator, the methods for aggregating metasearch engine results are divided into two kinds. One has a unique solution, and the other has multiple solutions. The proposed method not only has crisp weights, but also provides multiple aggregation results for decision makers to choose from. In order to prove the application of the ME-OWA operator method, under the context of aggregating metasearch engine results, an example is given, which shows the results obtained by the ME-OWA operator method and the minimax linear programming ( minimax-LP ) method. Comparison between these two methods are also made. The results show that the ME-OWA operator has nearly the same aggregation results as those of the minimax-LP method.展开更多
The observed images of the asteroid and the asteroid reference images are used to obtain the probe-to-asteroid direction and the location of the limb features of the asteroid in the inertial coordinate. These informa-...The observed images of the asteroid and the asteroid reference images are used to obtain the probe-to-asteroid direction and the location of the limb features of the asteroid in the inertial coordinate. These informa-tion in combination with the shape model of the asteroid and attitude information of the probe are utilized to ob-tain the position of the probe. The position information is then input to the UKF which determines the real-timeorbit of the probe. Finally, the autonomous orbit determination algorithm is validated using digital simulation.The determination of orbit using UKF is compared with that using extended Kalman filter (EKF), and the resultshows that UKF is superior to EKF.展开更多
In this paper,we present a smoothing Newton-like method for solving nonlinear systems of equalities and inequalities.By using the so-called max function,we transfer the inequalities into a system of semismooth equalit...In this paper,we present a smoothing Newton-like method for solving nonlinear systems of equalities and inequalities.By using the so-called max function,we transfer the inequalities into a system of semismooth equalities.Then a smoothing Newton-like method is proposed for solving the reformulated system,which only needs to solve one system of linear equations and to perform one line search at each iteration. The global and local quadratic convergence are studied under appropriate assumptions. Numerical examples show that the new approach is effective.展开更多
Migratory birds are known to use a number of different orientation cues to determine and maintain the direction of their movements. They are able to use at least three different sources of compass information, includi...Migratory birds are known to use a number of different orientation cues to determine and maintain the direction of their movements. They are able to use at least three different sources of compass information, including solar, stellar and geomagnetic cues. However, little is known about how these cues are calibrated into uniform reference direction information, while the hierarchy of these cues remains controversial. In recent studies, researchers suggest that avian migrants calibrate their geomagnetic compass on sunset cues, whereas others fail to find such patterns and insist on the prevalence of the magnetic compass. We carefully reviewed the existing literature and suggest that the conflicting results reported by different authors are due to genuine variation among species and propose hypotheses to explain this variation.展开更多
Double cost function linear quadratic regulator (DLQR) is developed from LQR theory to solve an optimal control problem with a general nonlinear cost function. In addition to the traditional LQ cost function, anothe...Double cost function linear quadratic regulator (DLQR) is developed from LQR theory to solve an optimal control problem with a general nonlinear cost function. In addition to the traditional LQ cost function, another free form cost function was introduced to express the physical need plainly and optimize weights of LQ cost function using the search algorithms. As an instance, DLQR was applied in determining the control input in the front steering angle compensation control (FSAC) model for heavy duty vehicles. The brief simulations show that DLQR is powerful enough to specify the engineering requirements correctly and balance many factors effectively. The concept and applicable field of LQR are expanded by DLQR to optimize the system with a free form cost function.展开更多
In order to study the sliding characteristics when the cable is connected with the other rods in the transmission line structures,a linear sliding cable element based on updated Lagrangian formulation and a sliding ca...In order to study the sliding characteristics when the cable is connected with the other rods in the transmission line structures,a linear sliding cable element based on updated Lagrangian formulation and a sliding catenary element considering the out-of-plane stiffness coefficient are put forward.A two-span and a three-span cable structures are taken as examples to verify the sliding cable elements.By comparing the tensions of the two proposed cable elements with the existing research results,the error is less than 1%,which proves the correctness of the proposed elements.The sliding characteristics should be considered in the practical engineering because of the significant difference between the tensions of sliding cable elements and those of cable element without considering sliding.The out-of-plane stiffness coefficient and friction characteristics do not obviously affect the cable tensions.展开更多
High dimensional data clustering,with the inherent sparsity of data and the existence of noise,is a serious challenge for clustering algorithms.A new linear manifold clustering method was proposed to address this prob...High dimensional data clustering,with the inherent sparsity of data and the existence of noise,is a serious challenge for clustering algorithms.A new linear manifold clustering method was proposed to address this problem.The basic idea was to search the line manifold clusters hidden in datasets,and then fuse some of the line manifold clusters to construct higher dimensional manifold clusters.The orthogonal distance and the tangent distance were considered together as the linear manifold distance metrics. Spatial neighbor information was fully utilized to construct the original line manifold and optimize line manifolds during the line manifold cluster searching procedure.The results obtained from experiments over real and synthetic data sets demonstrate the superiority of the proposed method over some competing clustering methods in terms of accuracy and computation time.The proposed method is able to obtain high clustering accuracy for various data sets with different sizes,manifold dimensions and noise ratios,which confirms the anti-noise capability and high clustering accuracy of the proposed method for high dimensional data.展开更多
In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wol...In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wolfe line search conditions.Numerical results show that the new method is efficient and stationary by comparing with PRP+ method,so it can be widely used in scientific computation.展开更多
In this paper,quadratic 0-1 programming problem (I) is considered, in terms of its features quadratic 0-1 programming problem is solved by linear approxity heurstic algrothm and a developed tabu search ahgrothm .
In this paper, we present a new form of successive approximation Broyden-like algorithm for nonlinear complementarity problem based on its equivalent nonsmooth equations. Under suitable conditions, we get the global c...In this paper, we present a new form of successive approximation Broyden-like algorithm for nonlinear complementarity problem based on its equivalent nonsmooth equations. Under suitable conditions, we get the global convergence on the algorithms. Some numerical results are also reported.展开更多
Based on the generalized Fischer-Burmeister function, Chen et al in 2008 put forward a regularization semismooth Newton method for solving the nonlinear complementarity problem with a P0-function. In this paper, we in...Based on the generalized Fischer-Burmeister function, Chen et al in 2008 put forward a regularization semismooth Newton method for solving the nonlinear complementarity problem with a P0-function. In this paper, we investigate the above algorithm with the monotone line search replaced by a non-monotone line search. It is shown that the non-monotone algorithm is well-defined, and is globally and locally superlinearly convergent under standard assumptions.展开更多
For physical ozone absorption without reaction,two parametric estimation methods,i.e.the common linear least square fitting and non-linear Simplex search methods,were applied,respectively,to determine the ozone mass t...For physical ozone absorption without reaction,two parametric estimation methods,i.e.the common linear least square fitting and non-linear Simplex search methods,were applied,respectively,to determine the ozone mass transfer coefficient during absorption and both methods give almost the same mass transfer coefficient.While for chemical absorption with ozone decomposition reaction,the common linear least square fitting method is not applicable for the evaluation of ozone mass transfer coefficient due to the difficulty of model linearization for describing ozone concentration dissolved in water.The nonlinear Simplex method obtains the mass transfer coefficient by minimizing the sum of the differences between the simulated and experimental ozone concentration during the whole absorption process,without the limitation of linear relationship between the dissolved ozone concentration and absorption time during the initial stage of absorption.Comparison of the ozone concentration profiles between the simulation and experimental data demonstrates that Simplex method may determine ozone mass transfer coefficient during absorption in an accurate and high efficiency way with wide applicability.展开更多
In this papert a recursive quadratic programming algorithm is proposed andstudied.The line search functions used are Han's nondifferentiable penalty functionswith a second order penalty term. In order to avoid mar...In this papert a recursive quadratic programming algorithm is proposed andstudied.The line search functions used are Han's nondifferentiable penalty functionswith a second order penalty term. In order to avoid maratos effect,Fukushima's mixeddirection is used as the direction of line search.Finallyt we prove the global convergenceand the local second order convergence of the algorithm.展开更多
Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameters in the search directions. In this note, by combining the nice numerical performance of PR and HS met...Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameters in the search directions. In this note, by combining the nice numerical performance of PR and HS methods with the global convergence property of the class of conjugate gradient methods presented by HU and STOREY(1991), a class of new restarting conjugate gradient methods is presented. Global convergences of the new method with two kinds of common line searches, are proved. Firstly, it is shown that, using reverse modulus of continuity function and forcing function, the new method for solving unconstrained optimization can work for a continously dif ferentiable function with Curry-Altman's step size rule and a bounded level set. Secondly, by using comparing technique, some general convergence properties of the new method with other kind of step size rule are established. Numerical experiments show that the new method is efficient by comparing with FR conjugate gradient method.展开更多
基金supported by National Key S&T Special Projects of Marine Carbonate(No.2008ZX05000-004)CNPC Projects(No.2008E-0610-10)
文摘With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel algorithms is a difficult issue.We introduce the idea of group search to the single-point search-based non-linear inversion algorithm, taking the quantum Monte Carlo method as an example for two-dimensional seismic wave velocity inversion and practical impedance inversion and test the calculation efficiency of using different node numbers.The results show the parallel algorithm in theoretical and practical data inversion is feasible and effective.The parallel algorithm has good versatility. The algorithm efficiency increases with increasing node numbers but the algorithm efficiency rate of increase gradually decreases as the node numbers increase.
基金the National Natural Science Foundation of China(10471062)the Natural Science Foundation of Jiangsu Province(BK2006184)~~
文摘A new limited memory symmetric rank one algorithm is proposed. It combines a modified self-scaled symmetric rank one (SSR1) update with the limited memory and nonmonotone line search technique. In this algorithm, the descent search direction is generated by inverse limited memory SSR1 update, thus simplifying the computation. Numerical comparison of the algorithm and the famous limited memory BFGS algorithm is given. Comparison results indicate that the new algorithm can process a kind of large-scale unconstrained optimization problems.
基金The National Natural Science Foundation of China (No.51308193)China Postdoctoral Science Foundation (No.20110491342)+1 种基金Jiangsu Planned Projects for Postdoctoral Research Funds(No.1101018C)the Science and Technology Project of State Grid Corporation of China(No.SGKJ[2007]116)
文摘In order to study the sliding characteristics when the cable structures are connected with other rods, a string of sliding cable dements (SCE) consisting of one active threenode SCE passing through the sliding point and multiple inactive two-node SCEs is put forward. Based on the updated Lagrangian formulation, the geometric nonlinear stiffness matrix of the three-node straight sliding cable dement is deduced. The examples about two-span and three-span continuous cable structures are studied to verify the effectiveness of the derived SCE. Comparing the cable tension of SCE with the existing research results, the calculating results show that the error is less than 1%. The sliding characteristics should be considered in practical engineering because of the obvious difference between the cable tension of the SCE and that of the cable element without considering sliding characteristics.
文摘A broadband microstrip patch antenna was analyzed and designed.Full wave analysis method(FWAM) was employed to show that a stacked microstrip dual patch antenna(SMDPA) might have a much wider bandwidth than that of the ordinanry uni patch one.By means of discrete complex image theory(DCIT),the Sommerfeld integrals (SI) involved were accurately calculated at a speed several hundred times faster than numerical integration method(NIM).The feeding structure of the SMDPA was then improved and the bandwidth was extended to about 22% or more for voltage standing wave ratio (VSWR)s≤2 Finally,a matching network was constructed to obtain a bandwidth of about 25% for s≤1.5.
文摘Based on the analysis to the random sear ch algorithm of LUUS, a modified random directed integer search algorithm (MRDI SA) is given for first time. And a practical example is given to show that the adva ntage of this kind of algorithm is the reliability can’t be infuenced by the ini tial value X (0) and the start search domain R (0) . Besides, i t can be applied to solve the higher dimensional constrained nonlinear integer p rogramming problem.
基金The National Natural Science Foundation of China(No.71171048)
文摘The maximal entropy ordered weighted averaging (ME-OWA) operator is used to aggregate metasearch engine results, and its newly analytical solution is also applied. Within the current context of the OWA operator, the methods for aggregating metasearch engine results are divided into two kinds. One has a unique solution, and the other has multiple solutions. The proposed method not only has crisp weights, but also provides multiple aggregation results for decision makers to choose from. In order to prove the application of the ME-OWA operator method, under the context of aggregating metasearch engine results, an example is given, which shows the results obtained by the ME-OWA operator method and the minimax linear programming ( minimax-LP ) method. Comparison between these two methods are also made. The results show that the ME-OWA operator has nearly the same aggregation results as those of the minimax-LP method.
文摘The observed images of the asteroid and the asteroid reference images are used to obtain the probe-to-asteroid direction and the location of the limb features of the asteroid in the inertial coordinate. These informa-tion in combination with the shape model of the asteroid and attitude information of the probe are utilized to ob-tain the position of the probe. The position information is then input to the UKF which determines the real-timeorbit of the probe. Finally, the autonomous orbit determination algorithm is validated using digital simulation.The determination of orbit using UKF is compared with that using extended Kalman filter (EKF), and the resultshows that UKF is superior to EKF.
基金supported by Guangdong Provincial Zhujiang Scholar Award Project,National Science Foundation of China(10671163,10871031)the National Basic Research Program under the Grant 2005CB321703Scientific Research Fund of Hunan Provincial Education Department(06A069,06C824)
文摘In this paper,we present a smoothing Newton-like method for solving nonlinear systems of equalities and inequalities.By using the so-called max function,we transfer the inequalities into a system of semismooth equalities.Then a smoothing Newton-like method is proposed for solving the reformulated system,which only needs to solve one system of linear equations and to perform one line search at each iteration. The global and local quadratic convergence are studied under appropriate assumptions. Numerical examples show that the new approach is effective.
基金supported by grant 12-04-00296-a from the Russian Foundation for Basic ResearchX.L. was supported by grant 60905060 from the National Natural Science Foundation of China+3 种基金grant BS2010DX025 from the Scientific Research Foundation for Excellent Middle-Aged and Youth Scientists of Shandong Province of Chinagrant 2011B11114 from the Fundamental Research Funds for Central UniversitiesThe writing of this review was inspired by the Bilateral Sino-Russian symposium ‘Animal Magnetic Navigation’ supported by grant 11-04-91216-NSFC-z from the Russian Foundation for Basic Researchgrant 61010164 from the National Natural Science Foundation of China
文摘Migratory birds are known to use a number of different orientation cues to determine and maintain the direction of their movements. They are able to use at least three different sources of compass information, including solar, stellar and geomagnetic cues. However, little is known about how these cues are calibrated into uniform reference direction information, while the hierarchy of these cues remains controversial. In recent studies, researchers suggest that avian migrants calibrate their geomagnetic compass on sunset cues, whereas others fail to find such patterns and insist on the prevalence of the magnetic compass. We carefully reviewed the existing literature and suggest that the conflicting results reported by different authors are due to genuine variation among species and propose hypotheses to explain this variation.
文摘Double cost function linear quadratic regulator (DLQR) is developed from LQR theory to solve an optimal control problem with a general nonlinear cost function. In addition to the traditional LQ cost function, another free form cost function was introduced to express the physical need plainly and optimize weights of LQ cost function using the search algorithms. As an instance, DLQR was applied in determining the control input in the front steering angle compensation control (FSAC) model for heavy duty vehicles. The brief simulations show that DLQR is powerful enough to specify the engineering requirements correctly and balance many factors effectively. The concept and applicable field of LQR are expanded by DLQR to optimize the system with a free form cost function.
基金Project(51308193)supported by the National Natural Science Foundation of ChinaProject(SGKJ[2007]116)supported by the Science and Technology Program of State Grid Corporation of China
文摘In order to study the sliding characteristics when the cable is connected with the other rods in the transmission line structures,a linear sliding cable element based on updated Lagrangian formulation and a sliding catenary element considering the out-of-plane stiffness coefficient are put forward.A two-span and a three-span cable structures are taken as examples to verify the sliding cable elements.By comparing the tensions of the two proposed cable elements with the existing research results,the error is less than 1%,which proves the correctness of the proposed elements.The sliding characteristics should be considered in the practical engineering because of the significant difference between the tensions of sliding cable elements and those of cable element without considering sliding.The out-of-plane stiffness coefficient and friction characteristics do not obviously affect the cable tensions.
基金Project(60835005) supported by the National Nature Science Foundation of China
文摘High dimensional data clustering,with the inherent sparsity of data and the existence of noise,is a serious challenge for clustering algorithms.A new linear manifold clustering method was proposed to address this problem.The basic idea was to search the line manifold clusters hidden in datasets,and then fuse some of the line manifold clusters to construct higher dimensional manifold clusters.The orthogonal distance and the tangent distance were considered together as the linear manifold distance metrics. Spatial neighbor information was fully utilized to construct the original line manifold and optimize line manifolds during the line manifold cluster searching procedure.The results obtained from experiments over real and synthetic data sets demonstrate the superiority of the proposed method over some competing clustering methods in terms of accuracy and computation time.The proposed method is able to obtain high clustering accuracy for various data sets with different sizes,manifold dimensions and noise ratios,which confirms the anti-noise capability and high clustering accuracy of the proposed method for high dimensional data.
基金Supported by the Fund of Chongqing Education Committee(KJ091104)
文摘In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wolfe line search conditions.Numerical results show that the new method is efficient and stationary by comparing with PRP+ method,so it can be widely used in scientific computation.
文摘In this paper,quadratic 0-1 programming problem (I) is considered, in terms of its features quadratic 0-1 programming problem is solved by linear approxity heurstic algrothm and a developed tabu search ahgrothm .
文摘In this paper, we present a new form of successive approximation Broyden-like algorithm for nonlinear complementarity problem based on its equivalent nonsmooth equations. Under suitable conditions, we get the global convergence on the algorithms. Some numerical results are also reported.
基金Supported by the Science Technology Development Plan of Tianjin (No.06YFGZGX05600)
文摘Based on the generalized Fischer-Burmeister function, Chen et al in 2008 put forward a regularization semismooth Newton method for solving the nonlinear complementarity problem with a P0-function. In this paper, we investigate the above algorithm with the monotone line search replaced by a non-monotone line search. It is shown that the non-monotone algorithm is well-defined, and is globally and locally superlinearly convergent under standard assumptions.
基金Project(2011467001)supported by the Ministry of Environment Protection of ChinaProject(2010DFB94130)supported by the Ministry of Science and Technology of China
文摘For physical ozone absorption without reaction,two parametric estimation methods,i.e.the common linear least square fitting and non-linear Simplex search methods,were applied,respectively,to determine the ozone mass transfer coefficient during absorption and both methods give almost the same mass transfer coefficient.While for chemical absorption with ozone decomposition reaction,the common linear least square fitting method is not applicable for the evaluation of ozone mass transfer coefficient due to the difficulty of model linearization for describing ozone concentration dissolved in water.The nonlinear Simplex method obtains the mass transfer coefficient by minimizing the sum of the differences between the simulated and experimental ozone concentration during the whole absorption process,without the limitation of linear relationship between the dissolved ozone concentration and absorption time during the initial stage of absorption.Comparison of the ozone concentration profiles between the simulation and experimental data demonstrates that Simplex method may determine ozone mass transfer coefficient during absorption in an accurate and high efficiency way with wide applicability.
文摘In this papert a recursive quadratic programming algorithm is proposed andstudied.The line search functions used are Han's nondifferentiable penalty functionswith a second order penalty term. In order to avoid maratos effect,Fukushima's mixeddirection is used as the direction of line search.Finallyt we prove the global convergenceand the local second order convergence of the algorithm.
文摘Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameters in the search directions. In this note, by combining the nice numerical performance of PR and HS methods with the global convergence property of the class of conjugate gradient methods presented by HU and STOREY(1991), a class of new restarting conjugate gradient methods is presented. Global convergences of the new method with two kinds of common line searches, are proved. Firstly, it is shown that, using reverse modulus of continuity function and forcing function, the new method for solving unconstrained optimization can work for a continously dif ferentiable function with Curry-Altman's step size rule and a bounded level set. Secondly, by using comparing technique, some general convergence properties of the new method with other kind of step size rule are established. Numerical experiments show that the new method is efficient by comparing with FR conjugate gradient method.