Linearized ridge-path method for function minimization

A modification based on a linearization of a ridge-path optimization method is presented. The linearized ridge-path method is a nongradient, conjugate direction method which converges quadratically in half the number of search directions required for Powell's method of conjugate directions. The ridge-path method and its modification are compared with some basic algorithms, namely, univariate method, steepest descent method, Powell's conjugate direction method, conjugate gradient method, and variable-metric method. The assessment indicates that the ridge-path method, with modifications, could present a promising technique for optimization.This work was in partial fulfillment of the requirements for the MS degree of the first author at Cairo University, Cairo, Egypt. The authors would like to acknowledge the helpful and constructive suggestions of the reviewer.     

双层规划问题的一种模式搜索Filter方法

双层规划在工程设计和经济管理中应用广泛,结合模式搜索方法和Filter方法提出了一种解决双层规划问题的算法—模式搜索Filter方法.算法以Filter法思想构造接受准则,以模式搜索提供迭代方向和步长,能够有效的解决一类双层规划问题.     

A new hybrid method for solving global optimization problem

In this paper we present a new hybrid method, called the SASP method. The purpose of this method is the hybridization of the simulated annealing (SA) with the descent method, where we estimate the gradient using simultaneous perturbation. Firstly, the new hybrid method finds a local minimum using the descent method, then SA is executed in order to escape from the currently discovered local minimum to a better one, from which the descent method restarts a new local search, and so on until convergence.The new hybrid method can be widely applied to a class of global optimization problems for continuous functions with constraints. Experiments on 30 benchmark functions, including high dimensional functions, show that the new method is able to find near optimal solutions efficiently. In addition, its performance as a viable optimization method is demonstrated by comparing it with other existing algorithms. Numerical results improve the robustness and efficiency of the method presented.     

NGLM:一类全局收敛的Newton-GMRES方法

本文提出了一类具有全局收敛性质的Newton-GMRES方法—NGLM方法．该方法是对经典Newton—GMRES方法的推广．NGLM方法的全局策略是当在非精确Newton方向上后退不能成功时，转而在一个子空间上运用信赖域方法确定迭代步长．理论分析与数值实验均表明，NGLM方法改善了Newton—GMRES方法的强健性．     

Finding roots of a real polynomial simultaneously by means of Bairstow's method

Aberth's method for finding the roots of a polynomial was shown to be robust. However, complex arithmetic is needed in this method even if the polynomial is real, because it starts with complex initial approximations. A novel method is proposed for real polynomials that does not require any complex arithmetic within iterations. It is based on the observation that Aberth's method is a systematic use of Newton's method. The analogous technique is then applied to Bairstow's procedure in the proposed method. As a result, the method needs half the computations per iteration than Aberth's method. Numerical experiments showed that the new method exhibited a competitive overall performance for the test polynomials.     

用不确定德尔菲法预测GDP

德尔菲法是一种建立在专家意见基础上的预测评估方法.不确定统计是利用不确定理论收集和整理分析专家数据的一种统计方法,其中关键的一点是如何构造不确定变量的不确定分布.把德尔菲法和不确定统计相结合,就得到了一种估计不确定分布的新方法——不确定德尔菲法.对该方法的估计误差进行了改进,得到了一种预测GDP的新方法,并利用其预测邯郸市的生产总值(GDP).     

一种求解非线性互补问题的方法及其收敛性

本文将Newton方法和外梯度方法相结合,提出了一种求解非线性互补问题的方法,证明了此方法的全局收敛性和超线性收敛性,在适当的条件下给出了一个有限终止结果。数值实验表明,此方法是有效的。     

一种基于矩阵初等变换的Schmidt正交化方法

对向量组的Schmidt正交化法和合同变换法的关系进行了分析,指出Schmidt正交化法就是合同变换法中利用规范化初等变换后的一种特殊情况,由此给出一种基于矩阵初等变换的Schmidt正交化方法——Schmidt初等变换正交化法,以及这一方法在软件Matlab上实现的程序.     

拟Raabe判别法与拟对数判别法的强弱关系

拟Raabe判别法是新近提出的关于正项级数收敛性的一种比较细致的判别法.对通项递减的正项级数来说,此判别法强于传统的Raabe判别法与Gauss判别法.通过对拟Raabe判别法与另一个细致的判别法——拟对数判别法强弱关系的探讨,得出了后一判别法强于前者的结论.     

An Improved Hybrid Method for Inverse Obstacle Scattering Problems

An improved hybrid method is introduced in this paper as a numerical method to reconstruct the scatterer by far-field pattern for just one incident direction with unknown physical properties of the scatterer. The improved hybrid method inherits the idea of the hybrid method by Kress and Serranho which is a combination of Newton and decomposition method, and it improves the hybrid method by introducing a general boundary condition. The numerical experiments show the feasibility of this method.     

Comparative analysis of the cutting angle and simulated annealing methods in global optimization

This article presents a comparative analysis of two methods of global optimization: the simulated annealing method and a method based on a combination of the cutting angle method and a local search. This analysis is carried out using results of numerical experiments. These results demonstrate that the combined method is more effective than the simulated annealing method.     

Iterative Methods of Richardson-Lucy-Type for Image Deblurring

Image deconvolution problems with a symmetric point-spread function arise in many areas of science and engineering. These problems often are solved by the Richardson-Lucy method, a nonlinear iterative method. We first show a convergence result for the Richardson-Lucy method. The proof sheds light on why the method may converge slowly. Subsequently, we describe an iterative active set method that imposes the same constraints on the computed solution as the Richardson-Lucy method. Computed examples show the latter method to yield better restorations than the Richardson-Lucy method and typically require less computational effort.     

求线性规划问题可行基的一种方法

文章给出了一般情形下从线性规划问题的标准型求可行基的一种方法，并通过与大Ｍ法、两阶段法及文［１］方法进行对比分析，说明这是一种有效可行且有可能较简便的方法     

Efficient numerical methods for pricing American options under stochastic volatility

Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two‐dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M‐matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one‐dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid‐based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower. ©John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

一个解凸二次规划的预测-校正光滑化方法

本文为凸二次规划问题提出一个光滑型方法,它是Engelke和Kanzow提出的解线性规划的光滑化算法的推广。其主要思想是将二次规划的最优性K-T条件写成一个非线性非光滑方程组,并利用Newton型方法来解其光滑近似。本文的方法是预测-校正方法。在较弱的条件下,证明了算法的全局收敛性和超线性收敛性。     

基于非线性尺度不变性的CCG法

１．引言 ＣＧ法对于变量个数很多的问题,是很有用的．１９７０年后它有了许多改进和发展,ＣＣＧ法以正定圆锥函数为基础［１］,它的一般方法是：设圆锥函数为 ２］其中： Ｖ＝ Ｖ（ｘ）＝１＋ ａＴｘ ≠ ０;, ｒ ∈Ｒ１为常量; ａ,ｇ ∈ Ｒｎ为常向量;ｘ ∈ Ｒｎ为变向量;Ａ∈Ｒｎ×ｎ为对称正定矩阵．算法［１］：预先给出初始近似点ｘ０∈ Ｒｎ及初始搜索方向 ｐ０;满足：其中“Ｉ”是单位矩阵, Ｖ０＝ Ｖ（ｘ０）＝ １＋ ａｔｘ０及记号“”是函数的梯度．迭代格式为： ｘｋ＋１＝ ｘｋ ＋λｋｐｋ,ｋ＝ ０,１,２,…（３…     

SOLVING TRUST REGION PROBLEM IN LARGE SCALEOPTIMIZATION

1.IntroductionLetHbeagivennxnsymmetricpositivesemidefinitematrixandcER".Inthispaperweconsiderthefollowingquadraticprogrammingwithasimplequadraticalconstraint1r--XTH. CTX=min2(1)s.tIlxlls5a,wheretheparameteraisprescribed.Thisproblemoccursfrequentlyintrustrationmethodforunconstrainedoptimization[11.Anumberofapproachesforsolving(1)havebeenproposedintheliterature[2,3,6,12--16].OnetechniqueistoapproximateaLagrangemultiplierAbyNewton'smethod.Theapproximationofthisparametermaybequitedelicate,ho…     

Derivative-Free Filter Simulated Annealing Method for Constrained Continuous Global Optimization

In this paper, a simulated-annealing-based method called Filter Simulated Annealing (FSA) method is proposed to deal with the constrained global optimization problem. The considered problem is reformulated so as to take the form of optimizing two functions, the objective function and the constraint violation function. Then, the FSA method is applied to solve the reformulated problem. The FSA method invokes a multi-start diversification scheme in order to achieve an efficient exploration process. To deal with the considered problem, a filter-set-based procedure is built in the FSA structure. Finally, an intensification scheme is applied as a final stage of the proposed method in order to overcome the slow convergence of SA-based methods. The computational results obtained by the FSA method are promising and show a superior performance of the proposed method, which is a point-to-point method, against population-based methods.     

求解一类非线性方程的改进的平均法

本文研究确定非线性方程一致有效近似解的平均法,得到B方法(Krylov-Bogoliubov method)与KBM方法(Krylov-Bogoliubov-Mitropolski method)的改进形式,通过两个例子与多尺度方法的比较,说明改进的平均法的有效性,从而拓广了平均法的应用范围。     

A domain decomposition method based on the iterative operator splitting method

In this article a new approach is proposed for constructing a domain decomposition method based on the iterative operator splitting method. The convergence properties of such a method are studied. The main feature of the proposed idea is the decoupling of space and time. We present a multi-iterative operator splitting method that combines iteratively the space and time splitting. We confirm with numerical applications the effectiveness of the proposed iterative operator splitting method in comparison with the classical Schwarz waveform relaxation method as a standard method for domain decomposition. We provide improved results and convergence rates.