On the low rank solution of the Q‐weighted nearest correlation matrix problem

The low rank solution of the Q‐weighted nearest correlation matrix problem is studied in this paper. Based on the property of Q‐weighted norm and the Gramian representation, we first reformulate the Q‐weighted nearest correlation matrix problem as a minimization problem of the trace function with quadratic constraints and then prove that the solution of the minimization problem the trace function is the stationary point of the original problem if it is rank‐deficient. Finally, the nonmonotone spectral projected gradient method is constructed to solve them. Numerical examples illustrate that the new method is feasible and effective. Copyright © 2015 John Wiley & Sons, Ltd.     

公路工程评标中的模糊一致矩阵方法

公路工程评标是一项多目标的复杂决策过程,在公路工程多方案分析评价的层次分析决策中应用模糊一致矩阵方法,有效的避开了模糊综合评判中隶属度确定问题.为公路工程评价方案决策这类定性与定量因素并存、正相关与负相关混杂的多因素、多层次评价提供了另一类有效的评判方法,通过公路工程评标方案的分析得到公路工程方案的排序结果,实例分析证明模糊一致矩阵方法是可行的.     

随机P矩阵和随机P0矩阵线性互补问题

定义了随机P矩阵和随机P0矩阵,给出了矩阵为随机P矩阵或随机P0矩阵的充要条件.研究了随机线性互补问题(SLCP)的矩阵为随机P矩阵时,期望残差方法(ERM)解集的有界性.得到了期望矩阵为P矩阵时,(ERM)解集非空有界.并且研究离散情形(ERM)与期望值方法(EV)解的关系,给出了(ERM)解唯一的条件.     

双对称线性矩阵方程的最佳逼近解

本文讨论了wang和Chang的双线件矩阵方程(ATXA,BTXB):(C,D)对称解的一致性条件.利用Hilbert空间的投影定理、商奇异值分解及其通解表达式和典型相关分解(CCD)的有效工具,获得了关于这个矩形方阵对的最小二乘问题的明确的解析表达式反对称(或最小Frobenius范数反对称解作为特例)最佳逼近解.     

Anderson acceleration of the alternating projections method for computing the nearest correlation matrix

In a wide range of applications it is required to compute the nearest correlation matrix in the Frobenius norm to a given symmetric but indefinite matrix. Of the available methods with guaranteed convergence to the unique solution of this problem the easiest to implement, and perhaps the most widely used, is the alternating projections method. However, the rate of convergence of this method is at best linear, and it can require a large number of iterations to converge to within a given tolerance. We show that Anderson acceleration, a technique for accelerating the convergence of fixed-point iterations, can be applied to the alternating projections method and that in practice it brings a significant reduction in both the number of iterations and the computation time. We also show that Anderson acceleration remains effective, and indeed can provide even greater improvements, when it is applied to the variants of the nearest correlation matrix problem in which specified elements are fixed or a lower bound is imposed on the smallest eigenvalue. Alternating projections is a general method for finding a point in the intersection of several sets and ours appears to be the first demonstration that this class of methods can benefit from Anderson acceleration.     

有资格限制的指派问题的求解方法

在实际的指派工作中，常会遇到某个人有没有资格去承担某项工作的问题，因此，本建立了有资格限制的指派问题的数学模型。在此数学模型中，将效益矩阵转化为判定矩阵，由此给出了判定此种指派问题是否有解的方法；在有解的情况下，进一步将效益矩阵转化为求解矩阵，从而将有资格限制的指派问题化为传统的指派问题来求解。最后给出了一个数值例子来说明这样的处理方法是有效的。     

On the reduction of Tikhonov minimization problems and the construction of regularization matrices

Tikhonov regularization replaces a linear discrete ill-posed problem by a penalized least-squares problem, whose solution is less sensitive to errors in the data and round-off errors introduced during the solution process. The penalty term is defined by a regularization matrix and a regularization parameter. The latter generally has to be determined during the solution process. This requires repeated solution of the penalized least-squares problem. It is therefore attractive to transform the least-squares problem to simpler form before solution. The present paper describes a transformation of the penalized least-squares problem to simpler form that is faster to compute than available transformations in the situation when the regularization matrix has linearly dependent columns and no exploitable structure. Properties of this kind of regularization matrices are discussed and their performance is illustrated.     

Structure methods for solving the nearest correlation matrix problem

The nearest correlation matrix problem is to find a positive semidefinite matrix with unit diagonal, that is, nearest in the Frobenius norm to a given symmetric matrix A. This problem arises in the finance industry, where the correlations are between stocks. In this paper, we formulate this problem as a smooth unconstrained minimization problem, for which rapid convergence can be obtained. Other methods are also studied. Comparative numerical results are reported.     

Minimizing the condition number of a positive definite matrix by completion

Summary. We consider the problem of minimizing the spectral condition number of a positive definite matrix by completion: \noindent where is an Hermitian positive definite matrix, a matrix and is a free Hermitian matrix. We reduce this problem to an optimization problem for a convex function in one variable. Using the minimal solution of this problem we characterize the complete set of matrices that give the minimum condition number. Received October 15, 1993     

POD reduced-order modeling for evolution equations utilizing arbitrary finite element discretizations

The main focus of the present work is the inclusion of spatial adaptivity for the snapshot computation in the offline phase of model order reduction utilizing proper orthogonal decomposition (POD-MOR) for nonlinear parabolic evolution problems. We consider snapshots which live in different finite element spaces, which means in a fully discrete setting that the snapshots are vectors of different length. From a numerical point of view, this leads to the problem that the usual POD procedure which utilizes a singular value decomposition of the snapshot matrix, cannot be carried out. In order to overcome this problem, we here construct the POD model/basis using the eigensystem of the correlation matrix (snapshot Gramian), which is motivated from a continuous perspective and is set up explicitly, e.g., without the necessity of interpolating snapshots into a common finite element space. It is an advantage of this approach that the assembly of the matrix only requires the evaluation of inner products of snapshots in a common Hilbert space. This allows a great flexibility concerning the spatial discretization of the snapshots. The analysis for the error between the resulting POD solution and the true solution reveals that the accuracy of the reduced-order solution can be estimated by the spatial and temporal discretization error as well as the POD error. Finally, to illustrate the feasibility of our approach, we present a test case of the Cahn–Hilliard system utilizing h-adapted hierarchical meshes and two settings of a linear heat equation using nested and non-nested grids.     

偏最小二乘(PartialLeast Square)方法的拟合指标及其在满意度研究中的应用

本文在对顾客满意度模型及PLS方法进行简单介绍的基础上,对PLS的拟合指标,包括共同因子、多元相关平方和冗余,进行了讨论。     

Weighted least squares solutions to general coupled Sylvester matrix equations

This paper is concerned with weighted least squares solutions to general coupled Sylvester matrix equations. Gradient based iterative algorithms are proposed to solve this problem. This type of iterative algorithm includes a wide class of iterative algorithms, and two special cases of them are studied in detail in this paper. Necessary and sufficient conditions guaranteeing the convergence of the proposed algorithms are presented. Sufficient conditions that are easy to compute are also given. The optimal step sizes such that the convergence rates of the algorithms, which are properly defined in this paper, are maximized and established. Several special cases of the weighted least squares problem, such as a least squares solution to the coupled Sylvester matrix equations problem, solutions to the general coupled Sylvester matrix equations problem, and a weighted least squares solution to the linear matrix equation problem are simultaneously solved. Several numerical examples are given to illustrate the effectiveness of the proposed algorithms.     

The Inverse Problem Method on the Half-Line and a Nested System of Riccati Equations

A mixed problem for the nonlinear Bogoyavlenskii system on the half-line is studied by the inverse problem method. The solution of the mixed problem is reduced to the solution of the inverse spectral problem of recovering a forth-order differential operator on the half-line from the Weyl matrix. We derive evolution equations for the elements of the Weyl matrix and give an algorithm for the solution of the mixed problem. Evolution equations of the elements of the Weyl matrix are nonlinear. It is shown that they can be reduced to a nested system of three successively solvable matrix Riccati equations.     

Perturbation Identities for Regularized Tikhonov Inverses and Weighted Pseudoinverses

We consider the perturbation analysis of two important problems for solving ill-conditioned or rank-deficient linear least squares problems. The Tikhonov regularized problem is a linear least squares problem with a regularization term balancing the size of the residual against the size of the weighted solution. The weight matrix can be a non-square matrix (usually with fewer rows than columns). The minimum-norm problem is the minimization of the size of the weighted solutions given by the set of solutions to the, possibly rank-deficient, linear least squares problem.It is well known that the solution of the Tikhonov problem tends to the minimum-norm solution as the regularization parameter of the Tikhonov problem tends to zero. Using this fact and the generalized singular value decomposition enable us to make a perturbation analysis of the minimum-norm problem with perturbation results for the Tikhonov problem. From the analysis we attain perturbation identities for Tikhonov inverses and weighted pseudoinverses.     

Studying the stability of solutions to systems of linear inequalities and constructing separating hyperplanes

We consider the methods for matrix correction or correction of all parameters of systems of linear equations and inequalities. We show that the problem of matrix correction of an inconsistent system of linear inequalities with the nonnegativity condition is reduced to a linear programming problem. Some stability measure is defined for a given solution to a system of linear inequalities as the minimal possible variation of parameters under which this solution does not satisfy the system. The problem of finding the most stable solution to the system is considered. The results are applied to constructing an optimal separating hyperplane in the feature space that is the most stable to the changes of features of the objects.     

Equivalence of the Generalized Vertical Block Linear Complementarity Problems and Linear Complementarity Problems

In this paper, generalization of a vertical block linear complementarity problem associated with two different types of matrices, one of which is a square matrix and the other is a vertical block matrix, is proposed. The necessary and sufficient conditions for the existence of the solution of the generalized vertical block linear complementarity problem is derived and the relationship between the solution set of the generalized vertical block linear complementarity problem and the linear complementarity problem is established. It is proved that the generalized vertical block linear complementarity problem has the P-property if and only if the vertical block linear complementarity problem has the P-property.     

Mean–variance asset–liability management with asset correlation risk and insurance liabilities

Consider an insurer who invests in the financial market where correlations among risky asset returns are randomly changing over time. The insurer who faces the risk of paying stochastic insurance claims needs to manage her asset and liability by taking into account of the correlation risk. This paper investigates the impact of correlation risk to the optimal asset–liability management (ALM) of an insurer. We employ the Wishart process to model the stochastic covariance matrix of risky asset returns. The insurer aims to minimize the variance of the terminal wealth given an expected terminal wealth subject to the risk of paying out random liabilities of compound Poisson process. This ALM problem then becomes a linear–quadratic stochastic optimal control problem with stochastic volatilities, stochastic correlations and jumps. The recognition of an affine form in the solution process enables us to derive the explicit closed-form solution to the optimal ALM portfolio policy, obtain the efficient frontier, and identify the condition that the solution is well behaved.     

Factorization of moving-average spectral densities by state-space representations and stacking

To factorize a spectral density matrix of a vector moving average process, we propose a state space representation. Although this state space is not necessarily of minimal dimension, its associated system matrices are simple and most matrix multiplications involved are nothing but index shifting. This greatly reduces the complexity of computation. Moreover, in this article we stack every q consecutive observations of the original process MA(q) and generate a vector MA(1) process. We consider a similar state space representation for the stacked process. Consequently, the solution hinges on a surprisingly compact discrete algebraic Riccati equation (DARE), which involves only one Toeplitz and one Hankel block matrix composed of autocovariance functions. One solution to this equation is given by the so-called iterative projection algorithm. Each iteration of the stacked version is equivalent to q iterations of the unstacked one. We show that the convergence behavior of the iterative projection algorithm is characterized by the decreasing rate of the partial correlation coefficients for the stacked process. In fact, the calculation of the partial correlation coefficients via the Whittle algorithm, which takes a very simple form in this case, offers another solution to the problem. To achieve computational efficiency, we apply the general Newton procedure given by Lancaster and Rodman to the DARE and obtain an algorithm of quadratic convergence rate. One immediate application of the new algorithms is polynomial stabilization. We also discuss various issues such as check of positivity and numerical implementation.     

Regularization of the solution of the cauchy problem for the generalized Moisil-Theodoresco system

We consider the problem of analytic continuation of a solution of the generalized Moisil-Theodoresco system in a spatial domain on the basis of its values on part of the boundary of this domain, that is, the Cauchy problem. We construct an approximate solution of this problem on the basis of the Carleman matrix method.