线性Hamilton系统边值问题的保辛数值方法

以Hamilton系统的正则变换和生成函数为基础研究线性时变Hamilton系统边值问题的保辛数值求解算法.根据第二类生成函数系数矩阵与状态传递矩阵的关系,构造了生成函数系数矩阵的区段合并递推算法,并进一步将递推算法推广到线性非齐次边值问题中;然后利用生成函数的性质将边值问题转化为初值问题,最后采用初值问题的保辛算法求解以达到整个Hamilton系统保辛的目的.数值算例表明该方法能够有效地求解线性齐次与非齐次问题,并能很好地保持Hamilton系统的固有特性.     

广义系统的Hamilton矩阵与H\-2代数Riccati方程的稳定化解

研究线性连续广义系统的Hamilton矩阵及H\-2代数Riccati方程.　提出一个标准的广义H\-2代数Riccati方程及对应的Hamilton矩阵，给出该Hamilton矩阵的几个重要性质. 在此基础上，得到该广义H\-2代数Riccati方程的稳定化解存在的一个充分条件并给出求解方法.此条件具有一般性, 主要定理是正常系统相应结果的推广.     

不等式约束优化一个可行序列线性方程组算法

提出了求解非线性不等式约束优化问题的一个可行序列线性方程组算法. 在每次迭代中, 可行下降方向通过求解两个线性方程组产生, 系数矩阵具有较好的稀疏性. 在较为温和的条件下, 算法具有全局收敛性和强收敛性, 数值试验表明算法是有效的.     

人工蜂群算法在LRIP机会约束模型中的应用

针对现有研究中未考虑配送阶段客户随机需求的问题,本文采用在一定置信区间上满足客户需求的方法,描述这种客户需求不确定的约束,在此基础上,建立了选址-路径-库存问题(Location-Routing-Inventory Problem,LRIP)的机会约束模型。提出人工蜂群算法(Artificial Bee Colony algorithm,ABC)对该问题模型进行优化求解。结合问题特征和邻域知识,提出了一种基于矩阵的编码方法,构造了启发式初始化方法,设计了2种基于矩阵编码的交换策略,在此基础上构造了5种蜂群搜索算子。通过仿真实验,分析比较了初始化方法和5种搜索策略;同时将人工蜂群算法与两阶段法进行了比较,优化结果证明人工蜂群算法是求解LRIP问题的有效方法。     

求解Hamilton约束的新算法

利用吴方法对多项式类型带约束的Hamilton系统作了研究.给出了判断系统是否正则的一个新算法.对于正则系统,可以得到Hamilton函数和运动方程,而对退化的系统给出了两个求解约束的新算法,得到带约束的Hamilton函数和运动方程.利用符号计算软件,这几个算法都可以在计算机上实现.     

非线性极大极小问题一个新的QP-free算法

本文研究非线性无约束极大极小优化问题. QP-free算法是求解光滑约束优化问题的有效方法之一,但用于求解极大极小优化问题的成果甚少.基于原问题的稳定点条件,既不需含参数的指数型光滑化函数,也不要等价光滑化,提出了求解非线性极大极小问题一个新的QP-free算法.新算法在每一次迭代中,通过求解两个相同系数矩阵的线性方程组获得搜索方向.在合适的假设条件下,该算法具有全局收敛性.最后,初步的数值试验验证了算法的有效性.     

无约束优化问题的对角稀疏拟牛顿法

对无约束优化问题提出了对角稀疏拟牛顿法,该算法采用了Armijo非精确线性搜索,并在每次迭代中利用对角矩阵近似拟牛顿法中的校正矩阵,使计算搜索方向的存贮量和工作量明显减少,为大型无约束优化问题的求解提供了新的思路．在通常的假设条件下,证明了算法的全局收敛性,线性收敛速度并分析了超线性收敛特征。数值实验表明算法比共轭梯度法有效,适于求解大型无约束优化问题．     

多重纳什均衡解的粒子群优化算法

提出了一种求解双矩阵对策多重纳什均衡解的粒子群优化算法。该算法通过随机初始点以及迭代粒子的归一化,保证粒子群始终保持在对策的可行策略空间内,避免了在随机搜索中产生无效的粒子,提高了粒子群优化算法求解纳什均衡解的计算性能。最后给出了几个数值例子,说明了粒子群优化算法的高效性。     

求解大规模Hamilton矩阵特征问题的辛Lanczos算法的误差分析

对求解大规模稀疏Hamilton矩阵特征问题的辛Lanczos算法给出了舍入误差分析.分析表明辛Lanczos算法在无中断时,保Hamilton结构的限制没有破坏非对称Lanczos算法的本质特性.本文还讨论了辛Lanczos算法计算出的辛Lanczos向量的J一正交性的损失与Ritz值收敛的关系.结论正如所料,当某些Ritz值开始收敛时.计算出的辛Lanczos向量的J-正交性损失是必然的.以上结果对辛Lanczos算法的改进具有理论指导意义.     

用奇异值分解方法计算具有重特征值矩阵的特征矢量

若当(Jordan)形是矩阵在相似条件下的一个标准形,在代数理论及其工程应用中都具有十分重要的意义．针对具有重特征值的矩阵,提出了一种运用奇异值分解方法计算它的特征矢量及若当形的算法．大量数值例子的计算结果表明,该算法在求解具有重特征值的矩阵的特征矢量及若当形上效果良好,优于商用软件MATLAB和MATHEMATICA．     

Numerically Stable Generation of Correlation Matrices and Their Factors

Correlation matrices—symmetric positive semidefinite matrices with unit diagonal—are important in statistics and in numerical linear algebra. For simulation and testing it is desirable to be able to generate random correlation matrices with specified eigenvalues (which must be nonnegative and sum to the dimension of the matrix). A popular algorithm of Bendel and Mickey takes a matrix having the specified eigenvalues and uses a finite sequence of Givens rotations to introduce 1s on the diagonal. We give improved formulae for computing the rotations and prove that the resulting algorithm is numerically stable. We show by example that the formulae originally proposed, which are used in certain existing Fortran implementations, can lead to serious instability. We also show how to modify the algorithm to generate a rectangular matrix with columns of unit 2-norm. Such a matrix represents a correlation matrix in factored form, which can be preferable to representing the matrix itself, for example when the correlation matrix is nearly singular to working precision.     

A FEAST algorithm with oblique projection for generalized eigenvalue problems

The contour integral‐based eigensolvers are the recent efforts for computing the eigenvalues inside a given region in the complex plane. The best‐known members are the Sakurai–Sugiura method, its stable version CIRR, and the FEAST algorithm. An attractive computational advantage of these methods is that they are easily parallelizable. The FEAST algorithm was developed for the generalized Hermitian eigenvalue problems. It is stable and accurate. However, it may fail when applied to non‐Hermitian problems. Recently, a dual subspace FEAST algorithm was proposed to extend the FEAST algorithm to non‐Hermitian problems. In this paper, we instead use the oblique projection technique to extend FEAST to the non‐Hermitian problems. Our approach can be summarized as follows: (a) construct a particular contour integral to form a search subspace containing the desired eigenspace and (b) use the oblique projection technique to extract desired eigenpairs with appropriately chosen test subspace. The related mathematical framework is established. Comparing to the dual subspace FEAST algorithm, we can save the computational cost roughly by a half if only the eigenvalues or the eigenvalues together with their right eigenvectors are needed. We also address some implementation issues such as how to choose a suitable starting matrix and design‐efficient stopping criteria. Numerical experiments are provided to illustrate that our method is stable and efficient.     

Target re-identification based on adaptive incremental KISS measure learning

Target re-identification from across cameras is a difficult problem in multi-camera surveillance, which needs to be urgently solved. Traditional solutions, in addition to relying on the statistical characteristics of targets’ appearance, are more often using excellent measurement algorithms. Among many such algorithms, the Keep It Simple and Stupid Measure Learning (KISSME) algorithm based on statistical probability is an outstanding one. But it has a problem that the eigenvalue is not stable, and the actual matching rate is relatively low. So, in this paper, we optimize the measurement algorithms based on large scale Keep It Simple and Stupid (KISS) measure learning. From elements, such as inadequate sample, size and smaller or larger eigenvalues, we introduce eigenvalue stabilization technique, and finally form our algorithm which can be called Adaptive Incremental Keep It Simple and Stupid Measure Learning (AIKISSME). Finally, through many experiments based on Viewpoint Invariant Pedestrian Recognition (VIPeR) and by comparing with other algorithms, this work concludes that AIKISSME achieves the best overall performance.     

Superharmonic instability of stokes waves

A stability of nearly limiting Stokes waves to superharmonic perturbations is considered numerically in approximation of an infinite depth. Investigation of the stability properties can give one an insight into the evolution of the Stokes wave. The new, previously inaccessible branches of superharmonic instability were investigated. Our numerical simulations suggest that eigenvalues of linearized dynamical equations, corresponding to the unstable modes, appear as a result of a collision of a pair of purely imaginary eigenvalues at the origin, and a subsequent appearance of a pair of purely real eigenvalues: a positive and a negative one that are symmetric with respect to zero. Complex conjugate pairs of purely imaginary eigenvalues correspond to stable modes, and as the steepness of the underlying Stokes wave grows, the pairs move toward the origin along the imaginary axis. Moreover, when studying the eigenvalues of linearized dynamical equations we find that as the steepness of the Stokes wave grows, the real eigenvalues follow a universal scaling law, that can be approximated by a power law. The asymptotic power law behavior of this dependence for instability of Stokes waves close to the limiting one is proposed. Surface elevation profiles for several unstable eigenmodes are made available through  http://stokeswave.org website.     

A modified bisection algorithm for the determination of the eigenvalues of a symmetric tridiagonal matrix

Summary A modification to the well known bisection algorithm [1] when used to determine the eigenvalues of a real symmetric matrix is presented. In the new strategy the terms in the Sturm sequence are computed only as long as relevant information on the required eigenvalues is obtained. The resulting algorithm usingincomplete Sturm sequences can be shown to minimise the computational work required especially when only a few eigenvalues are required.The technique is also applicable to other computational methods which use the bisection process.     

The {\cal {D Q R}} algorithm, basic theory, convergence, and conditional stability

Summary. We describe a fast matrix eigenvalue algorithm that uses a matrix factorization and reverse order multiply technique involving three factors and that is based on the symmetric matrix factorization as well as on –orthogonal reduction techniques where is computed from the given matrix . It operates on a similarity reduction of a real matrix to general tridiagonal form and computes all of 's eigenvalues in operations, where the part of the operations is possibly performed over , instead of the 7–8 real flops required by the eigenvalue algorithm. Potential breakdo wn of the algorithm can occur in the reduction to tridiagonal form and in the –orthogonal reductions. Both, however, can be monitored during the computations. The former occurs rather rarely for dimensions and can essentially be bypassed, while the latter is extremely rare and can be bypassed as well in our conditionally stable implementation of the steps. We prove an implicit theorem which allows implicit shifts, give a convergence proof for the algorithm and show that is conditionally stable for general balanced tridiagonal matrices . Received April 25, 1995 / Revised version received February 9, 1996     

TheAB-algorithm and its properties

One presents a new algorithm, called the -algorithm, for solving the generalized eigenvalue problem Ax=λBx, where det (A—λB) ≠ 0, relative to A. The algorithm is iterative, it is based on the application of plane rotations and allows us to pass from the initial problem to the solving of a similar problem having simpler matrices, whose eigenvalues can be easily computed and coincide with the eigenvalues of the initial problem. Thus, if all the eigenvalues of the initial problem are distinct, then the application of the -algorithm leads to the computation of the eigenvalues of a pencil with triangular matrices. In the case of an arbitrary initial pencil A—λB, the problem reduces to solving the eigenvalue problem for a pencil of quasitriangular form. One proves the convergence of the algorithm. One establishes its properties which in many respects are similar with the properties of the known algorithms QR and QZ, the first of which solves the usual eigenvalue problem while the second one solves the generalized problem of the above-mentioned form.     

Numerically stable deflation of hessenberg and symmetric tridiagonal matrices

While numerically stable techniques have been available for deflating a fulln byn matrix, no satisfactory finite technique has been known which preserves Hessenberg form. We describe a new algorithm which explicitly deflates a Hessenberg matrix in floating point arithmetic by means of a sequence of plane rotations. When applied to a symmetric tridiagonal matrix, the deflated matrix is again symmetric tridiagonal. Repeated deflation can be used to find an orthogonal set of eigenvectors associated with any selection of eigenvalues of a symmetric tridiagonal matrix.     

Numerical Determination of a Canonical Form of a Symplectic Matrix

We propose an algorithm that transforms a real symplectic matrix with a stable structure to a block diagonal form composed of three main blocks. The two extreme blocks of the same size are associated respectively with the eigenvalues outside and inside the unit circle. Moreover, these eigenvalues are symmetric with respect to the unit circle. The central block is in turn composed of several diagonal blocks whose eigenvalues are on the unit circle and satisfy a modification of the Krein-Gelfand-Lidskii criterion. The proposed algorithm also gives a qualitative criterion for structural stability.     

Eigenvalues of a symmetric tridiagonal matrix: A divide-and-conquer approach

Summary The Symmetric Tridiagonal Eigenproblem has been the topic of some recent work. Many methods have been advanced for the computation of the eigenvalues of such a matrix. In this paper, we present a divide-and-conquer approach to the computation of the eigenvalues of a symmetric tridiagonal matrix via the evaluation of the characteristic polynomial. The problem of evaluation of the characteristic polynomial is partitioned into smaller parts which are solved and these solutions are then combined to form the solution to the original problem. We give the update equations for the characteristic polynomial and certain auxiliary polynomials used in the computation. Furthermore, this set of recursions can be implemented on a regulartree structure. If the concurrency exhibited by this algorithm is exploited, it can be shown that thetime for computation of all the eigenvalues becomesO(nlogn) instead ofO(n ²) as is the case for the approach where the order is increased by only one at every step. We address the numerical problems associated with the use of the characteristic polynomial and present a numerically stable technique for the eigenvalue computation.