Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide ra...Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide range of potential applications in computer vision and pattern recognition. In this paper, we propose a new model to deal with the TRPCA problem by an alternation minimization algorithm along with two adaptive rankadjusting strategies. For the underlying low-rank tensor, we simultaneously perform low-rank matrix factorizations to its all-mode matricizations; while for the underlying sparse tensor,a soft-threshold shrinkage scheme is applied. Our method can be used to deal with the separation between either an exact or an approximate low-rank tensor and a sparse one. We established the subsequence convergence of our algorithm in the sense that any limit point of the iterates satisfies the KKT conditions. When the iteration stops, the output will be modified by applying a high-order SVD approach to achieve an exactly low-rank final result as the accurate rank has been calculated. The numerical experiments demonstrate that our method could achieve better results than the compared methods.展开更多
The present paper spreads the principal axis intrinsic method to the highdimensional case and discusses the solution of the tensor equation AX --XA = C
Tensor complementarity problem (TCP) is a special kind of nonlinear complementarity problem (NCP). In this paper, we introduce a new class of structure tensor and give some examples. By transforming the TCP to the sys...Tensor complementarity problem (TCP) is a special kind of nonlinear complementarity problem (NCP). In this paper, we introduce a new class of structure tensor and give some examples. By transforming the TCP to the system of nonsmooth equations, we develop a semismooth Newton method for the tensor complementarity problem. We prove the monotone convergence theorem for the proposed method under proper conditions.展开更多
We propose a generalized Lanczos method to generate the many-body basis states of quantum lattice models using tensor-network states （TNS）. The ground-state wave function is represented as a linear superposition com...We propose a generalized Lanczos method to generate the many-body basis states of quantum lattice models using tensor-network states （TNS）. The ground-state wave function is represented as a linear superposition composed from a set of TNS generated by Lanczos iteration. This method improves significantly the accuracy of the tensor-network algorithm and provides an effective way to enlarge the maximal bond dimension of TNS. The ground state such obtained contains significantly more entanglement than each individual TNS, reproducing correctly the logarithmic size dependence of the entanglement entropy in a critical system. The method can be generalized to non-Hamiltonian systems and to the calculation of low-lying excited states, dynamical correlation functions, and other physical properties of strongly correlated systems.展开更多
This paper studies the problem of tensor principal component analysis (PCA). Usually the tensor PCA is viewed as a low-rank matrix completion problem via matrix factorization technique, and nuclear norm is used as a c...This paper studies the problem of tensor principal component analysis (PCA). Usually the tensor PCA is viewed as a low-rank matrix completion problem via matrix factorization technique, and nuclear norm is used as a convex approximation of the rank operator under mild condition. However, most nuclear norm minimization approaches are based on SVD operations. Given a matrix , the time complexity of SVD operation is O(mn2), which brings prohibitive computational complexity in large-scale problems. In this paper, an efficient and scalable algorithm for tensor principal component analysis is proposed which is called Linearized Alternating Direction Method with Vectorized technique for Tensor Principal Component Analysis (LADMVTPCA). Different from traditional matrix factorization methods, LADMVTPCA utilizes the vectorized technique to formulate the tensor as an outer product of vectors, which greatly improves the computational efficacy compared to matrix factorization method. In the experiment part, synthetic tensor data with different orders are used to empirically evaluate the proposed algorithm LADMVTPCA. Results have shown that LADMVTPCA outperforms matrix factorization based method.展开更多
Based on the equivalence between the Sylvester tensor equation and the linear equation obtained by discretization of partial differential equations(PDEs),an overlapping Schwarz alternative method based on the tensor f...Based on the equivalence between the Sylvester tensor equation and the linear equation obtained by discretization of partial differential equations(PDEs),an overlapping Schwarz alternative method based on the tensor format and an overlapping parallel Schwarz method based on the tensor format for solving high-dimensional PDEs are proposed.The complexity of the new algorithms is discussed.Finally,the feasibility and effectiveness of the new methods are verified by some numerical examples.展开更多
Recovering the low-rank structure of data matrix from sparse errors arises in the principal component pursuit （PCP）. This paper exploits the higher-order generalization of matrix recovery, named higher-order princip...Recovering the low-rank structure of data matrix from sparse errors arises in the principal component pursuit （PCP）. This paper exploits the higher-order generalization of matrix recovery, named higher-order principal component pursuit （HOPCP）, since it is critical in multi-way data analysis. Unlike the convexification （nuclear norm） for matrix rank function, the tensorial nuclear norm is stil an open problem. While existing preliminary works on the tensor completion field provide a viable way to indicate the low complexity estimate of tensor, therefore, the paper focuses on the low multi-linear rank tensor and adopt its convex relaxation to formulate the convex optimization model of HOPCP. The paper further propose two algorithms for HOPCP based on alternative minimization scheme： the augmented Lagrangian alternating direction method （ALADM） and its truncated higher-order singular value decomposition （ALADM-THOSVD） version. The former can obtain a high accuracy solution while the latter is more efficient to handle the computationally intractable problems. Experimental results on both synthetic data and real magnetic resonance imaging data show the applicability of our algorithms in high-dimensional tensor data processing.展开更多
The modified matrix method of construction of wavefield on the free surface of an anisotropic medium is proposed. The earthquake source represented by a randomly oriented force or a seismic moment tensor is placed on ...The modified matrix method of construction of wavefield on the free surface of an anisotropic medium is proposed. The earthquake source represented by a randomly oriented force or a seismic moment tensor is placed on an arbitrary boundary of a layered anisotropic medium. The theory of the matrix propagator in a homogeneous anisotropic medium by introducing a ＂wave propagator＂ is presented. It is shown that the matrix propagator can be represented by a ＂wave propagator＂ in each layer for anisotropic layered medium. The matrix propagator P（z, z0=0） acts on the free surface of the layered medium and generates stress-displacement vector at depth z. The displacement field on the free surface of an anisotropic medium is obtained from the received system of equations considering the radiation condition and that the free surface is stressless. The new method determining source time function in anisotropic medium for three different types of seismic source is validated.展开更多
This paper studies the problem of recovering low-rank tensors, and the tensors are corrupted by both impulse and Gaussian noise. The problem is well accomplished by integrating the tensor nuclear norm and the l1-norm ...This paper studies the problem of recovering low-rank tensors, and the tensors are corrupted by both impulse and Gaussian noise. The problem is well accomplished by integrating the tensor nuclear norm and the l1-norm in a unified convex relaxation framework. The nuclear norm is adopted to explore the low-rank components and the l1-norm is used to exploit the impulse noise. Then, this optimization problem is solved by some augmented-Lagrangian-based algorithms. Some preliminary numerical experiments verify that the proposed method can well recover the corrupted low-rank tensors.展开更多
基金Supported by the National Natural Science Foundation of China(Grant Nos.6157209961320106008+2 种基金91230103)National Science and Technology Major Project(Grant Nos.2013ZX040050212014ZX04001011)
文摘Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide range of potential applications in computer vision and pattern recognition. In this paper, we propose a new model to deal with the TRPCA problem by an alternation minimization algorithm along with two adaptive rankadjusting strategies. For the underlying low-rank tensor, we simultaneously perform low-rank matrix factorizations to its all-mode matricizations; while for the underlying sparse tensor,a soft-threshold shrinkage scheme is applied. Our method can be used to deal with the separation between either an exact or an approximate low-rank tensor and a sparse one. We established the subsequence convergence of our algorithm in the sense that any limit point of the iterates satisfies the KKT conditions. When the iteration stops, the output will be modified by applying a high-order SVD approach to achieve an exactly low-rank final result as the accurate rank has been calculated. The numerical experiments demonstrate that our method could achieve better results than the compared methods.
文摘The present paper spreads the principal axis intrinsic method to the highdimensional case and discusses the solution of the tensor equation AX --XA = C
文摘Tensor complementarity problem (TCP) is a special kind of nonlinear complementarity problem (NCP). In this paper, we introduce a new class of structure tensor and give some examples. By transforming the TCP to the system of nonsmooth equations, we develop a semismooth Newton method for the tensor complementarity problem. We prove the monotone convergence theorem for the proposed method under proper conditions.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11190024 and 11474331)
文摘We propose a generalized Lanczos method to generate the many-body basis states of quantum lattice models using tensor-network states （TNS）. The ground-state wave function is represented as a linear superposition composed from a set of TNS generated by Lanczos iteration. This method improves significantly the accuracy of the tensor-network algorithm and provides an effective way to enlarge the maximal bond dimension of TNS. The ground state such obtained contains significantly more entanglement than each individual TNS, reproducing correctly the logarithmic size dependence of the entanglement entropy in a critical system. The method can be generalized to non-Hamiltonian systems and to the calculation of low-lying excited states, dynamical correlation functions, and other physical properties of strongly correlated systems.
文摘This paper studies the problem of tensor principal component analysis (PCA). Usually the tensor PCA is viewed as a low-rank matrix completion problem via matrix factorization technique, and nuclear norm is used as a convex approximation of the rank operator under mild condition. However, most nuclear norm minimization approaches are based on SVD operations. Given a matrix , the time complexity of SVD operation is O(mn2), which brings prohibitive computational complexity in large-scale problems. In this paper, an efficient and scalable algorithm for tensor principal component analysis is proposed which is called Linearized Alternating Direction Method with Vectorized technique for Tensor Principal Component Analysis (LADMVTPCA). Different from traditional matrix factorization methods, LADMVTPCA utilizes the vectorized technique to formulate the tensor as an outer product of vectors, which greatly improves the computational efficacy compared to matrix factorization method. In the experiment part, synthetic tensor data with different orders are used to empirically evaluate the proposed algorithm LADMVTPCA. Results have shown that LADMVTPCA outperforms matrix factorization based method.
基金supported by the National Natural Science Foundation of China(12161027)the Guangxi Natural Science Foundation of China(2020GXNSFAA159143)partially supported by the Science and Technology Project of Guangxi of China(Guike AD23023002).
文摘Based on the equivalence between the Sylvester tensor equation and the linear equation obtained by discretization of partial differential equations(PDEs),an overlapping Schwarz alternative method based on the tensor format and an overlapping parallel Schwarz method based on the tensor format for solving high-dimensional PDEs are proposed.The complexity of the new algorithms is discussed.Finally,the feasibility and effectiveness of the new methods are verified by some numerical examples.
基金supported by the National Natural Science Foundationof China(51275348)
文摘Recovering the low-rank structure of data matrix from sparse errors arises in the principal component pursuit （PCP）. This paper exploits the higher-order generalization of matrix recovery, named higher-order principal component pursuit （HOPCP）, since it is critical in multi-way data analysis. Unlike the convexification （nuclear norm） for matrix rank function, the tensorial nuclear norm is stil an open problem. While existing preliminary works on the tensor completion field provide a viable way to indicate the low complexity estimate of tensor, therefore, the paper focuses on the low multi-linear rank tensor and adopt its convex relaxation to formulate the convex optimization model of HOPCP. The paper further propose two algorithms for HOPCP based on alternative minimization scheme： the augmented Lagrangian alternating direction method （ALADM） and its truncated higher-order singular value decomposition （ALADM-THOSVD） version. The former can obtain a high accuracy solution while the latter is more efficient to handle the computationally intractable problems. Experimental results on both synthetic data and real magnetic resonance imaging data show the applicability of our algorithms in high-dimensional tensor data processing.
文摘The modified matrix method of construction of wavefield on the free surface of an anisotropic medium is proposed. The earthquake source represented by a randomly oriented force or a seismic moment tensor is placed on an arbitrary boundary of a layered anisotropic medium. The theory of the matrix propagator in a homogeneous anisotropic medium by introducing a ＂wave propagator＂ is presented. It is shown that the matrix propagator can be represented by a ＂wave propagator＂ in each layer for anisotropic layered medium. The matrix propagator P（z, z0=0） acts on the free surface of the layered medium and generates stress-displacement vector at depth z. The displacement field on the free surface of an anisotropic medium is obtained from the received system of equations considering the radiation condition and that the free surface is stressless. The new method determining source time function in anisotropic medium for three different types of seismic source is validated.
文摘This paper studies the problem of recovering low-rank tensors, and the tensors are corrupted by both impulse and Gaussian noise. The problem is well accomplished by integrating the tensor nuclear norm and the l1-norm in a unified convex relaxation framework. The nuclear norm is adopted to explore the low-rank components and the l1-norm is used to exploit the impulse noise. Then, this optimization problem is solved by some augmented-Lagrangian-based algorithms. Some preliminary numerical experiments verify that the proposed method can well recover the corrupted low-rank tensors.
