非线性半定规划一个全局收敛的无罚无滤子SSDP算法

提出了一个求解非线性半定规划的无罚函数无滤子序列二次半定规划(SSDP)算法. 算法每次迭代只需求解一个二次半定规划子问题确定搜索方向; 非单调线搜索保证目标函数或约束违反度函数的充分下降, 从而产生新的迭代点. 在适当的假设条件下, 证明了算法的全局收敛性. 最后给出了初步的数值实验结果.     

不等式约束优化一个基于滤子思想的广义梯度投影算法

讨论非线性不等式约束优化问题, 借鉴于滤子算法思想,提出了一个新型广义梯度投影算法.该方法既不使用罚函数又无真正意义下的滤子.每次迭代通过一个简单的显式广义投影法产生搜索方向,步长由目标函数值或者约束违反度函数值充分下降的Armijo型线搜索产生.算法的主要特点是: 不需要迭代序列的有界性假设;不需要传统滤子算法所必需的可行恢复阶段;使用了ε积极约束集减小计算量.在合适的假设条件下算法具有全局收敛性, 最后对算法进行了初步的数值实验.     

无罚函数和滤子的QP-free非可行域方法

提出了求解光滑不等式约束最优化问题的无罚函数和无滤子QP-free非可行域方法. 通过乘子和非线性互补函数, 构造一个等价于原约束问题一阶KKT条件的非光滑方程组. 在此基础上, 通过牛顿-拟牛顿迭代得到满足KKT最优性条件的解, 在迭代中采用了无罚函数和无滤子线搜索方法, 并证明该算法是可实现,具有全局收敛性. 另外, 在较弱条件下可以证明该方法具有超线性收敛性.     

一种无恢复过程的SQP-滤子法

本文研究非线性不等式约束优化问题,构造一个新的SQP-滤子法.该方法将滤子技术有机融合到简金宝提出的可行SQP方法中,利用转轴运算的思想,产生一个近似积极约束集,当QP子问题不相容时,利用广义投影技术获得可行搜索方向.该算法既能避免罚函数的选择,又能避免常规滤子算法中的恢复算法,一定程度上简化了计算.最后,在合理的条件下,证明了算法的全局收敛性.     

一般约束优化基于识别函数的模松弛算法

借助于半罚函数和产生工作集的识别函数以及模松弛SQP算法思想, 本文建立了求解带等式及不等式约束优化的一个新算法. 每次迭代中, 算法的搜索方向由一个简化的二次规划子问题及一个简化的线性方程组产生. 算法在不包含严格互补性的温和条件下具有全局收敛性和超线性收敛性. 最后给出了算法初步的数值试验报告.     

无罚函数和滤子的一个新的QP-free方法

通过构造一个等价于原约束问题一阶KKT条件的非光滑方程组, 提出一类新的QP-free方法. 在迭代中采用了无罚函数和无滤子线搜索方法, 在此基础上, 通过牛顿-拟牛顿迭代得到满足KKT最优条件的解, 并证明该算法是可实现、具有全局收敛性. 另外, 在较弱条件下可以证明该方法具有超线性收敛性.     

一类求解非线性约束优化问题的线搜索渐缩滤子算法（英文）

针对非线性等式约束优化问题,本文给出一种新的线搜索滤子算法.算法中将非线性等式约束优化问题的最优性条件作为滤子,并在接受准则中加入渐缩函数,使得当线搜索试探步长减小时时滤子包络的越来越薄,从而使得试探步被接受程度更有弹性,不会被当前的迭代点拒绝.在适当的假设下,证明算法的全局收敛性,并给出算法初步的数值实验结果.     

一个解界约束非线性方程组的无导数回溯线搜索仿射内点信赖域方法（英文）

文章给出了一个求解界约束非线性方程组的无导数回溯线搜索仿射内点信赖域方法.该方法利用非线性方程组的特点,对方程组中每一个函数建立插值模型.通过利用信赖域模型和回溯先搜索技术的结合,利用插值信赖域子问题子问题求解搜索方向,并利用回溯先搜索技术保证可行性.在合理的假设条件下,证明了算法的全局和快速局部收敛性.并且,通过数值实验表明该种无导数算法对求解界约束非线性方程组问题是有效的.     

不等式约束优化全局收敛的滤子线性方程组算法

设计了求解不等式约束非线性规划问题的一种新的滤子序列线性方程组算法,该算法每步迭代由减小约束违反度和目标函数值两部分构成.利用约束函数在某个中介点线性化的方法产生搜索方向.每步迭代仅需求解两个线性方程组,计算量较小.在一般条件下,证明了算法产生的无穷迭代点列所有聚点都是可行点并且所有聚点都是所求解问题的KKT点.     

带等式约束二次规划子问题的滤子SQP算法

一类求解非线性规划问题的滤子序列二次规划(SQP)方法被提出.为了提高收敛速度,给目标函数和约束违反度函数都设置了斜边界.二次规划子问题(QP)设置为两项:不等式约束QP和等式约束QP.两个子问题产生的搜索方向进行线性迭加后为算法的搜索方向.这样的设置可以改善收敛性,并调节算法运行中的一些不良效果.在较温和的条件下,可得到全局收敛性.     

Determining a Rotation of a Tetrahedron from a Projection

The following problem, arising from medical imaging, is addressed: Suppose that T is a known tetrahedron in ?³ with centroid at the origin. Also known is the orthogonal projection U of the vertices of the image ?T of T under an unknown rotation ? about the origin. Under what circumstances can ? be determined from T and U?     

On a mapping of a projective space into a sphere

We obtain an exact estimate for the minimum multiplicity of a continuous finite-to-one mapping of a projective space into a sphere for all dimensions. For finite-to-one mappings of a projective space into a Euclidean space, we obtain an exact estimate for this multiplicity for n = 2, 3. For n ≥ 4, we prove that this estimate does not exceed 4. Several open questions are formulated.     

Calculating a function of a matrix with a real spectrum

Let T be a square matrix with a real spectrum, and let f be an analytic function. The problem of the approximate calculation of f(T) is discussed. Applying the Schur triangular decomposition and the reordering, one can assume that T is triangular and its diagonal entries t_ii are arranged in increasing order. To avoid calculations using the differences t_ii ? t_jj with close (including equal) t_ii and t_jj, it is proposed to represent T in a block form and calculate the two main block diagonals using interpolating polynomials. The rest of the f(T) entries can be calculated using the Parlett recurrence algorithm. It is also proposed to perform some scalar operations (such as the building of interpolating polynomials) with an enlarged number of significant decimal digits.

     

Evolution to a steady state for a rarefied gas flowing from a tank into a vacuum through a plane channel

A kinetic equation (S-model) is used to solve the nonstationary problem of a monatomic rarefied gas flowing from a tank of infinite capacity into a vacuum through a long plane channel. Initially, the gas is at rest and is separated from the vacuum by a barrier. The temperature of the channel walls is kept constant. The flow is found to evolve to a steady state. The time required for reaching a steady state is examined depending on the channel length and the degree of gas rarefaction. The kinetic equation is solved numerically by applying a conservative explicit finite-difference scheme that is firstorder accurate in time and second-order accurate in space. An approximate law is proposed for the asymptotic behavior of the solution at long times when the evolution to a steady state becomes a diffusion process.     

Torsional Oscillations of a Punch on a Layer Bonded to a Half-Space Containing a Cylindrical Cavity

A previously published algorithm [9] is implemented in application to the Reissner–Sagoci problem for a layer bonded to a half-space containing a cylindrical cavity. The influence of the mechanical and geometrical parameters of the layer and the half-space on the amplitude-frequency response curves of the punch oscillations is analyzed. Practical applications of the results are proposed for ensuring the seismic isolation of buildings on the investigated foundation in the presence of dynamic torsional excitations.     

Symmetry-based solution of a model for a combination of a risky investment and a riskless investment

Benth and Karlsen [F.E. Benth, K.H. Karlsen, A note on Merton's portfolio selection problem for the Schwartz mean-reversion model, Stoch. Anal. Appl. 23 (2005) 687-704] treated a problem of the optimisation of the selection of a portfolio based upon the Schwartz mean-reversion model. The resulting Hamilton-Jacobi-Bellman equation in 1+2 dimensions is quite nonlinear. The solution obtained by Benth and Karlsen was very ingenious. We provide a solution of the problem based on the application of the Lie theory of continuous groups to the partial differential equation and its associated boundary and terminal conditions.