The neural network partial least square (NNPLS) method was used to establish a robust reaction model for a multi-component catalyst of methane oxidative coupling. The details, including the learning algorithm, the num...The neural network partial least square (NNPLS) method was used to establish a robust reaction model for a multi-component catalyst of methane oxidative coupling. The details, including the learning algorithm, the number of hidden units of the inner network, activation function, initialization of the network weights and the principal components, are discussed. The results show that the structural organizations of inner neural network are 1-10-5-1, 1-8-4-1, 1-8-5-1, 1-7-4-1, 1-8-4-1, 1-8-6-1, respectively. The Levenberg-Marquardt method was used in the learning algorithm, and the central sigmoidal function is the activation function. Calculation results show that four principal components are convenient in the use of the multi-component catalyst modeling of methane oxidative coupling. Therefore a robust reaction model expressed by NNPLS succeeds in correlating the relations between elements in catalyst and catalytic reaction results. Compared with the direct network modeling, NNPLS model can be adjusted by experimental data and the calculation of the model is simpler and faster than that of the direct network model.展开更多
Based upon a generalized variational principle, which relaxed the inter element continuity requirements, a novel refined hybrid Mindlin plate element is developed, its non linear element stiffness matrices are decom...Based upon a generalized variational principle, which relaxed the inter element continuity requirements, a novel refined hybrid Mindlin plate element is developed, its non linear element stiffness matrices are decomposed into a series of matrices with respect to the assumed strain modes. The formulation presented in this paper is different from any other non linear mixed/hybrid element formulation all successful experience of linear hybrid formulation is absorbed into the formulation(adding non conforming modes and realizing orthogonalization) Numerical results show that the present approach is more effective than any other non linear hybrid element formulation over the accuracy and computational efficiency. In addition, non conforming modes can also overcome the shear locking effect.展开更多
Aim To determine the global optimal solution for a mine ventilation network under given network topology and airway characteristics. Methods\ The genetic algorithm was used to find the global optimal solution of the ...Aim To determine the global optimal solution for a mine ventilation network under given network topology and airway characteristics. Methods\ The genetic algorithm was used to find the global optimal solution of the network. Results\ A modified genetic algorithm is presented with its characteristics and principle. Instead of working on the conventional bit by bit operation, both the crossover and mutation operators are handled in real values by the proposed algorithms. To prevent the system from turning into a premature problem, the elitists from two groups of possible solutions are selected to reproduce the new populations. Conclusion\ The simulation results show that the method outperforms the conventional nonlinear programming approach whether from the viewpoint of the number of iterations required to find the optimum solutions or from the final solutions obtained.展开更多
Starting from a 3 × 3 matrix spectral problem, we derive a hierarchy of nonlinear equations. It is shown that the hierarchy possesses bi-Hamiltonian structure. Under the symmetry constraints between the potential...Starting from a 3 × 3 matrix spectral problem, we derive a hierarchy of nonlinear equations. It is shown that the hierarchy possesses bi-Hamiltonian structure. Under the symmetry constraints between the potentials and the eigenfunctions, Lax pair and adjoint Lax pairs including partial part and temporal part are nonlinearied into two finitedimensional Hamiltonian systems (FDHS) in Liouville sense. Moreover, an explicit N-fold Darboux transformation for CDNS equation is constructed with the help of a gauge transformation of the spectral problem.展开更多
In a recent article [Physics Letters A 372 (2008) 417], Wang et al. proposed a method, which is called the (G′/G)-expansion method, to look for travelling wave solutions of nonlinear evolution equations. The trav...In a recent article [Physics Letters A 372 (2008) 417], Wang et al. proposed a method, which is called the (G′/G)-expansion method, to look for travelling wave solutions of nonlinear evolution equations. The travelling wave solutions involving parameters of the KdV equation, the mKdV equation, the variant Boussinesq equations, and the Hirota-Satsuma equations are obtained by using this method. They think the (G′/G)-expansion method is a new method and more travelling wave solutions of many nonlinear evolution equations can be obtained. In this paper, we will show that the (G′/G)-expansion method is equivalent to the extended tanh function method.展开更多
In this paper, an improved nonlinear process fault detection method is proposed based on modified kernel partial least squares(KPLS). By integrating the statistical local approach(SLA) into the KPLS framework, two new...In this paper, an improved nonlinear process fault detection method is proposed based on modified kernel partial least squares(KPLS). By integrating the statistical local approach(SLA) into the KPLS framework, two new statistics are established to monitor changes in the underlying model. The new modeling strategy can avoid the Gaussian distribution assumption of KPLS. Besides, advantage of the proposed method is that the kernel latent variables can be obtained directly through the eigen value decomposition instead of the iterative calculation, which can improve the computing speed. The new method is applied to fault detection in the simulation benchmark of the Tennessee Eastman process. The simulation results show superiority on detection sensitivity and accuracy in comparison to KPLS monitoring.展开更多
Data reconciliation is an effective technique for providing accurate and consistent value for chemical process. However, the presence of gross errors can severely bias the reconciled results. Robust estimators can sig...Data reconciliation is an effective technique for providing accurate and consistent value for chemical process. However, the presence of gross errors can severely bias the reconciled results. Robust estimators can significantly reduce the effect of gross errors and yield less-biased results. In this article, a new method is proposed to solve the robust data reconciliation problem of nonlinear chemical process. By using several technologies including linearization method, penalty function, virtual observation equation, and equivalent weights method, the robust data reconciliation problem can be transformed into least squares estimator problem which leads to the convenience in computation. Simulation results in a nonlinear chemical process demonstrate the efficiency of the proposed method.展开更多
In this work,a variable structure control(VSC)technique is proposed to achieve satisfactory robustness for unstable processes.Optimal values of unknown parameters of VSC are obtained using Whale optimization algorithm...In this work,a variable structure control(VSC)technique is proposed to achieve satisfactory robustness for unstable processes.Optimal values of unknown parameters of VSC are obtained using Whale optimization algorithm which was recently reported in literature.Stability analysis has been done to verify the suitability of the proposed structure for industrial processes.The proposed control strategy is applied to three different types of unstable processes including non-minimum phase and nonlinear systems.A comparative study ensures that the proposed scheme gives superior performance over the recently reported VSC system.Furthermore,the proposed method gives satisfactory results for a cart inverted pendulum system in the presence of external disturbance and noise.展开更多
The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swa...The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swarms such as birds and ants when searching for food. In this article, first the particle swarm optimization algorithm was described in detail, and ant colony algorithm improved. Then the methods were applied to three different kinds of geophysical inversion problems: (1) a linear problem which is sensitive to noise, (2) a synchronous inversion of linear and nonlinear problems, and (3) a nonlinear problem. The results validate their feasibility and efficiency. Compared with the conventional genetic algorithm and simulated annealing, they have the advantages of higher convergence speed and accuracy. Compared with the quasi-Newton method and Levenberg-Marquardt method, they work better with the ability to overcome the locally optimal solutions.展开更多
Instead of the usual Hirota ansatz,i.e.,the functions in bilinear equations being chosen as exponentialtypes,a generalized Hirota ansatz is proposed for a (3+1)-dimensional nonlinear evolution equation.Based on theres...Instead of the usual Hirota ansatz,i.e.,the functions in bilinear equations being chosen as exponentialtypes,a generalized Hirota ansatz is proposed for a (3+1)-dimensional nonlinear evolution equation.Based on theresulting generalized Hirota ansatz,a family of new explicit solutions for the equation are derived.展开更多
Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolu...Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolution equation.展开更多
In this paper, a new 7×7 matrix spectral problem, which is associated with the super AKNS equation isconstructed.With the use of the binary nonlinearization method, a new integrable decomposition of the super AKN...In this paper, a new 7×7 matrix spectral problem, which is associated with the super AKNS equation isconstructed.With the use of the binary nonlinearization method, a new integrable decomposition of the super AKNSequation is presented.展开更多
The purpose of the present paper is twofold. First, the projective Riccati equations (PREs for short) are resolved by means of a linearized theorem, which was known in the literature. Based on the signs and values o...The purpose of the present paper is twofold. First, the projective Riccati equations (PREs for short) are resolved by means of a linearized theorem, which was known in the literature. Based on the signs and values of coeffcients of PREs, the solutions with two arbitrary parameters of PREs can be expressed by the hyperbolic functions, the trigonometric functions, and the rational functions respectively, at the same time the relation between the components of each solution to PREs is also implemented. Second, more new travelling wave solutions for some nonlinear PDEs, such as the Burgers equation, the mKdV equation, the NLS^+ equation, new Hamilton amplitude equation, and so on, are obtained by using Sub-ODE method, in which PREs are taken as the Sub-ODEs. The key idea of this method is that the travelling wave solutions of nonlinear PDE can be expressed by a polynomial in two variables, which are the components of each solution to PREs, provided that the homogeneous balance between the higher order derivatives and nonlinear terms in the equation is considered.展开更多
In this paper, an improved hybrid differential evolution-estimation of distribution algorithm (IHDE-EDA) is proposed for nonlinear programming (NLP) and mixed integer nonlinear programming (MINLP) models in engineerin...In this paper, an improved hybrid differential evolution-estimation of distribution algorithm (IHDE-EDA) is proposed for nonlinear programming (NLP) and mixed integer nonlinear programming (MINLP) models in engineering optimization fields. In order to improve the global searching ability and convergence speed, IHDE-EDA takes full advantage of differential information and global statistical information extracted respectively from differential evolution algorithm and annealing mechanism-embedded estimation of distribution algorithm. Moreover, the feasibility rules are used to handle constraints, which do not require additional parameters and can guide the population to the feasible region quickly. The effectiveness of hybridization mechanism of IHDE-EDA is first discussed, and then simulation and comparison based on three benchmark problems demonstrate the efficiency, accuracy and robustness of IHDE-EDA. Finally, optimization on an industrial-size scheduling of two-pipeline crude oil blending problem shows the practical applicability of IHDE-EDA.展开更多
In this work, by means of a new more general ansatz and the symbolic computation system Maple, we extend the Riccati equation rational expansion method [Chaos, Solitons & Fractals 25 (2005) 1019] to uniformly const...In this work, by means of a new more general ansatz and the symbolic computation system Maple, we extend the Riccati equation rational expansion method [Chaos, Solitons & Fractals 25 (2005) 1019] to uniformly construct a series of stochastic nontravelling wave solutions for nonlinear stochastic evolution equation. To illustrate the effectiveness of our method, we take the stochastic mKdV equation as an example, and successfully construct some new and more general solutions including a series of rational formal nontraveling wave and coefficient functions' soliton-like solution.s and trigonometric-like function solutions. The method can also be applied to solve other nonlinear stochastic evolution equation or equations.展开更多
In this paper, a global optimization algorithm is proposed for nonlinear sum of ratios problem (P). The algorithm works by globally solving problem (P1) that is equivalent to problem (P), by utilizing linearizat...In this paper, a global optimization algorithm is proposed for nonlinear sum of ratios problem (P). The algorithm works by globally solving problem (P1) that is equivalent to problem (P), by utilizing linearization technique a linear relaxation programming of the (P1) is then obtained. The proposed algorithm is convergent to the global minimum of (P1) through the successive refinement of linear relaxation of the feasible region of objective function and solutions of a series of linear relaxation programming. Numerical results indicate that the proposed algorithm is feasible and can be used to globally solve nonlinear sum of ratios problems (P).展开更多
Differential evolution (DE) is a global optimizer for continuous design variables. To enhance DE, it is necessary to handle discrete design variables. In this paper, a discrete differential evolution (DDE) algorit...Differential evolution (DE) is a global optimizer for continuous design variables. To enhance DE, it is necessary to handle discrete design variables. In this paper, a discrete differential evolution (DDE) algorithm is proposed to handle discrete design variables The proposed DDE is based on the DE/l/rand/bin method. In the proposed DDE, the mutation ratio is regarded as the exchange probability, and thus, no modifications of DE/l/rand/bin are required. In addition, in order to maintain diversity through the search process, we initialize all search points. By introducing the initialization of all search points, global or quasi-optimum solution can be found. We validate the proposed DDE by applying it to several benchmark problems.展开更多
Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameter in the search directions. In this note, conditions are given on the parameter in the conjugate gradie...Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameter in the search directions. In this note, conditions are given on the parameter in the conjugate gradient directions to ensure the descent property of the search directions. Global convergence of such a class of methods is discussed. It is shown that, using reverse modulus of continuity function and forcing function, the new method for solving unconstrained optimization can work for a continuously differentiable function with a modification of the Curry-Altman's step-size rule and a bounded level set. Combining PR method with our new method, PR method is modified to have global convergence property.Numerical experiments show that the new methods are efficient by comparing with FR conjugate gradient method.展开更多
文摘The neural network partial least square (NNPLS) method was used to establish a robust reaction model for a multi-component catalyst of methane oxidative coupling. The details, including the learning algorithm, the number of hidden units of the inner network, activation function, initialization of the network weights and the principal components, are discussed. The results show that the structural organizations of inner neural network are 1-10-5-1, 1-8-4-1, 1-8-5-1, 1-7-4-1, 1-8-4-1, 1-8-6-1, respectively. The Levenberg-Marquardt method was used in the learning algorithm, and the central sigmoidal function is the activation function. Calculation results show that four principal components are convenient in the use of the multi-component catalyst modeling of methane oxidative coupling. Therefore a robust reaction model expressed by NNPLS succeeds in correlating the relations between elements in catalyst and catalytic reaction results. Compared with the direct network modeling, NNPLS model can be adjusted by experimental data and the calculation of the model is simpler and faster than that of the direct network model.
文摘Based upon a generalized variational principle, which relaxed the inter element continuity requirements, a novel refined hybrid Mindlin plate element is developed, its non linear element stiffness matrices are decomposed into a series of matrices with respect to the assumed strain modes. The formulation presented in this paper is different from any other non linear mixed/hybrid element formulation all successful experience of linear hybrid formulation is absorbed into the formulation(adding non conforming modes and realizing orthogonalization) Numerical results show that the present approach is more effective than any other non linear hybrid element formulation over the accuracy and computational efficiency. In addition, non conforming modes can also overcome the shear locking effect.
文摘Aim To determine the global optimal solution for a mine ventilation network under given network topology and airway characteristics. Methods\ The genetic algorithm was used to find the global optimal solution of the network. Results\ A modified genetic algorithm is presented with its characteristics and principle. Instead of working on the conventional bit by bit operation, both the crossover and mutation operators are handled in real values by the proposed algorithms. To prevent the system from turning into a premature problem, the elitists from two groups of possible solutions are selected to reproduce the new populations. Conclusion\ The simulation results show that the method outperforms the conventional nonlinear programming approach whether from the viewpoint of the number of iterations required to find the optimum solutions or from the final solutions obtained.
基金The project supported by National Natural Science Foundation of China under Grant No. 10371023 and Shanghai Shuguang Project of China under Grant No. 02SG02
文摘Starting from a 3 × 3 matrix spectral problem, we derive a hierarchy of nonlinear equations. It is shown that the hierarchy possesses bi-Hamiltonian structure. Under the symmetry constraints between the potentials and the eigenfunctions, Lax pair and adjoint Lax pairs including partial part and temporal part are nonlinearied into two finitedimensional Hamiltonian systems (FDHS) in Liouville sense. Moreover, an explicit N-fold Darboux transformation for CDNS equation is constructed with the help of a gauge transformation of the spectral problem.
基金Supported by National Natural Science Foundation of China under Grant No. 10671172
文摘In a recent article [Physics Letters A 372 (2008) 417], Wang et al. proposed a method, which is called the (G′/G)-expansion method, to look for travelling wave solutions of nonlinear evolution equations. The travelling wave solutions involving parameters of the KdV equation, the mKdV equation, the variant Boussinesq equations, and the Hirota-Satsuma equations are obtained by using this method. They think the (G′/G)-expansion method is a new method and more travelling wave solutions of many nonlinear evolution equations can be obtained. In this paper, we will show that the (G′/G)-expansion method is equivalent to the extended tanh function method.
基金Supported by the Special Scientific Research of Selection and Cultivation of Excellent Young Teachers in Shanghai Universities(YYY11076)
文摘In this paper, an improved nonlinear process fault detection method is proposed based on modified kernel partial least squares(KPLS). By integrating the statistical local approach(SLA) into the KPLS framework, two new statistics are established to monitor changes in the underlying model. The new modeling strategy can avoid the Gaussian distribution assumption of KPLS. Besides, advantage of the proposed method is that the kernel latent variables can be obtained directly through the eigen value decomposition instead of the iterative calculation, which can improve the computing speed. The new method is applied to fault detection in the simulation benchmark of the Tennessee Eastman process. The simulation results show superiority on detection sensitivity and accuracy in comparison to KPLS monitoring.
基金Supported by the Funds for 0utstanding Young Researchers from the National Natural Science Foundation of China (No.60025308) and the Key Technologies R&D Program in the National "10th 5-year Plan" (No.2001BA204B07).
文摘Data reconciliation is an effective technique for providing accurate and consistent value for chemical process. However, the presence of gross errors can severely bias the reconciled results. Robust estimators can significantly reduce the effect of gross errors and yield less-biased results. In this article, a new method is proposed to solve the robust data reconciliation problem of nonlinear chemical process. By using several technologies including linearization method, penalty function, virtual observation equation, and equivalent weights method, the robust data reconciliation problem can be transformed into least squares estimator problem which leads to the convenience in computation. Simulation results in a nonlinear chemical process demonstrate the efficiency of the proposed method.
文摘In this work,a variable structure control(VSC)technique is proposed to achieve satisfactory robustness for unstable processes.Optimal values of unknown parameters of VSC are obtained using Whale optimization algorithm which was recently reported in literature.Stability analysis has been done to verify the suitability of the proposed structure for industrial processes.The proposed control strategy is applied to three different types of unstable processes including non-minimum phase and nonlinear systems.A comparative study ensures that the proposed scheme gives superior performance over the recently reported VSC system.Furthermore,the proposed method gives satisfactory results for a cart inverted pendulum system in the presence of external disturbance and noise.
基金supported by the 973 Program(Grant No 2007CB209600)Open Fund(No.GDL0706) of the Key Laboratory of Geo-detection(China University of Geosciences,Beijing),Ministry of Education
文摘The inversions of complex geophysical data always solve multi-parameter, nonlinear, and multimodal optimization problems. Searching for the optimal inversion solutions is similar to the social behavior observed in swarms such as birds and ants when searching for food. In this article, first the particle swarm optimization algorithm was described in detail, and ant colony algorithm improved. Then the methods were applied to three different kinds of geophysical inversion problems: (1) a linear problem which is sensitive to noise, (2) a synchronous inversion of linear and nonlinear problems, and (3) a nonlinear problem. The results validate their feasibility and efficiency. Compared with the conventional genetic algorithm and simulated annealing, they have the advantages of higher convergence speed and accuracy. Compared with the quasi-Newton method and Levenberg-Marquardt method, they work better with the ability to overcome the locally optimal solutions.
文摘Instead of the usual Hirota ansatz,i.e.,the functions in bilinear equations being chosen as exponentialtypes,a generalized Hirota ansatz is proposed for a (3+1)-dimensional nonlinear evolution equation.Based on theresulting generalized Hirota ansatz,a family of new explicit solutions for the equation are derived.
文摘Using a polynomial expansion method, the general exact solitary wave solution and singular one areconstructed for the non-linear KS equation. This approach is obviously applicable to a large variety of nonlinear evolution equation.
基金Supported by the National Natural Science Foundation of China under Grant No.10926036the Education Department of Zhejiang Province under Grant No.Y200906909the Zhejiang Provincial Natural Science Foundation of China under Grant No.Y6090172
文摘In this paper, a new 7×7 matrix spectral problem, which is associated with the super AKNS equation isconstructed.With the use of the binary nonlinearization method, a new integrable decomposition of the super AKNSequation is presented.
基金The project supported in part by the Natural Science Foundation of Education Department of Henan Province of China under Grant No. 2006110002 and the Science Foundations of Henan University of Science and Technology under Grant Nos. 2004ZD002 and 2006ZY001
文摘The purpose of the present paper is twofold. First, the projective Riccati equations (PREs for short) are resolved by means of a linearized theorem, which was known in the literature. Based on the signs and values of coeffcients of PREs, the solutions with two arbitrary parameters of PREs can be expressed by the hyperbolic functions, the trigonometric functions, and the rational functions respectively, at the same time the relation between the components of each solution to PREs is also implemented. Second, more new travelling wave solutions for some nonlinear PDEs, such as the Burgers equation, the mKdV equation, the NLS^+ equation, new Hamilton amplitude equation, and so on, are obtained by using Sub-ODE method, in which PREs are taken as the Sub-ODEs. The key idea of this method is that the travelling wave solutions of nonlinear PDE can be expressed by a polynomial in two variables, which are the components of each solution to PREs, provided that the homogeneous balance between the higher order derivatives and nonlinear terms in the equation is considered.
基金Supported by the National Basic Research Program of China (2012CB720500)the National Natural Science Foundation of China (60974008)
文摘In this paper, an improved hybrid differential evolution-estimation of distribution algorithm (IHDE-EDA) is proposed for nonlinear programming (NLP) and mixed integer nonlinear programming (MINLP) models in engineering optimization fields. In order to improve the global searching ability and convergence speed, IHDE-EDA takes full advantage of differential information and global statistical information extracted respectively from differential evolution algorithm and annealing mechanism-embedded estimation of distribution algorithm. Moreover, the feasibility rules are used to handle constraints, which do not require additional parameters and can guide the population to the feasible region quickly. The effectiveness of hybridization mechanism of IHDE-EDA is first discussed, and then simulation and comparison based on three benchmark problems demonstrate the efficiency, accuracy and robustness of IHDE-EDA. Finally, optimization on an industrial-size scheduling of two-pipeline crude oil blending problem shows the practical applicability of IHDE-EDA.
基金The author would like to thank the referees very much for their careful reading of the manuscript and many valuable suggestions.
文摘In this work, by means of a new more general ansatz and the symbolic computation system Maple, we extend the Riccati equation rational expansion method [Chaos, Solitons & Fractals 25 (2005) 1019] to uniformly construct a series of stochastic nontravelling wave solutions for nonlinear stochastic evolution equation. To illustrate the effectiveness of our method, we take the stochastic mKdV equation as an example, and successfully construct some new and more general solutions including a series of rational formal nontraveling wave and coefficient functions' soliton-like solution.s and trigonometric-like function solutions. The method can also be applied to solve other nonlinear stochastic evolution equation or equations.
基金Foundation item: Supported by the National Natural Science Foundation of China(10671057) Supported by the Natural Science Foundation of Henan Institute of Science and Technology(06054)
文摘In this paper, a global optimization algorithm is proposed for nonlinear sum of ratios problem (P). The algorithm works by globally solving problem (P1) that is equivalent to problem (P), by utilizing linearization technique a linear relaxation programming of the (P1) is then obtained. The proposed algorithm is convergent to the global minimum of (P1) through the successive refinement of linear relaxation of the feasible region of objective function and solutions of a series of linear relaxation programming. Numerical results indicate that the proposed algorithm is feasible and can be used to globally solve nonlinear sum of ratios problems (P).
文摘Differential evolution (DE) is a global optimizer for continuous design variables. To enhance DE, it is necessary to handle discrete design variables. In this paper, a discrete differential evolution (DDE) algorithm is proposed to handle discrete design variables The proposed DDE is based on the DE/l/rand/bin method. In the proposed DDE, the mutation ratio is regarded as the exchange probability, and thus, no modifications of DE/l/rand/bin are required. In addition, in order to maintain diversity through the search process, we initialize all search points. By introducing the initialization of all search points, global or quasi-optimum solution can be found. We validate the proposed DDE by applying it to several benchmark problems.
文摘Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameter in the search directions. In this note, conditions are given on the parameter in the conjugate gradient directions to ensure the descent property of the search directions. Global convergence of such a class of methods is discussed. It is shown that, using reverse modulus of continuity function and forcing function, the new method for solving unconstrained optimization can work for a continuously differentiable function with a modification of the Curry-Altman's step-size rule and a bounded level set. Combining PR method with our new method, PR method is modified to have global convergence property.Numerical experiments show that the new methods are efficient by comparing with FR conjugate gradient method.