In this work, we consider a homotopic principle for solving large-scale and dense l1underdetermined problems and its applications in image processing and classification. We solve the face recognition problem where the...In this work, we consider a homotopic principle for solving large-scale and dense l1underdetermined problems and its applications in image processing and classification. We solve the face recognition problem where the input image contains corrupted and/or lost pixels. The approach involves two steps: first, the incomplete or corrupted image is subject to an inpainting process, and secondly, the restored image is used to carry out the classification or recognition task. Addressing these two steps involves solving large scale l1minimization problems. To that end, we propose to solve a sequence of linear equality constrained multiquadric problems that depends on a regularization parameter that converges to zero. The procedure generates a central path that converges to a point on the solution set of the l1underdetermined problem. In order to solve each subproblem, a conjugate gradient algorithm is formulated. When noise is present in the model, inexact directions are taken so that an approximate solution is computed faster. This prevents the ill conditioning produced when the conjugate gradient is required to iterate until a zero residual is attained.展开更多
A joint two-dimensional(2D)direction-of-arrival(DOA)and radial Doppler frequency estimation method for the L-shaped array is proposed in this paper based on the compressive sensing(CS)framework.Revised from the conven...A joint two-dimensional(2D)direction-of-arrival(DOA)and radial Doppler frequency estimation method for the L-shaped array is proposed in this paper based on the compressive sensing(CS)framework.Revised from the conventional CS-based methods where the joint spatial-temporal parameters are characterized in one large scale matrix,three smaller scale matrices with independent azimuth,elevation and Doppler frequency are introduced adopting a separable observation model.Afterwards,the estimation is achieved by L1-norm minimization and the Bayesian CS algorithm.In addition,under the L-shaped array topology,the azimuth and elevation are separated yet coupled to the same radial Doppler frequency.Hence,the pair matching problem is solved with the aid of the radial Doppler frequency.Finally,numerical simulations corroborate the feasibility and validity of the proposed algorithm.展开更多
This paper presents a new method for image separation through employing a combined dictionary consisting of wavelets and complex shearlets. Because the combined dictionary sparsely represents points and curvilinear si...This paper presents a new method for image separation through employing a combined dictionary consisting of wavelets and complex shearlets. Because the combined dictionary sparsely represents points and curvilinear singularities respectively, the image can be decomposed into pointlike and curvelike parts as accurate as possible. The proposed method based on the geo- metric separation theory introduced by Donoho in 2005 shows that accurate geometric separation of the morphologically distinct fea- tures of points and curves can be achieved by l1 minimization. The experimental results show that the proposed method can not only be effective but also greatly reduce the computing time.展开更多
Radial functions have become a useful tool in numerical mathematics. On the sphere they have to be identified with the zonal functions. We investigate zonal polynomials with mass concentration at the pole, in the sens...Radial functions have become a useful tool in numerical mathematics. On the sphere they have to be identified with the zonal functions. We investigate zonal polynomials with mass concentration at the pole, in the sense of their L1-norm is attaining the minimum value. Such polynomials satisfy a complicated system of nonlinear e-quations (algebraic if the space dimension is odd, only) and also a singular differential equation of third order. The exact order of decay of the minimum value with respect to the polynomial degree is determined. By our results we can prove that some nodal systems on the sphere, which are defined by a minimum-property, are providing fundamental matrices which are diagonal-dominant or bounded with respect to the ∞-norm, at least, as the polynomial degree tends to infinity.展开更多
We consider efficient methods for the recovery of block sparse signals from underdetermined system of linear equations. We show that if the measurement matrix satisfies the block RIP with δ2s 〈 0.4931, then every bl...We consider efficient methods for the recovery of block sparse signals from underdetermined system of linear equations. We show that if the measurement matrix satisfies the block RIP with δ2s 〈 0.4931, then every block s-sparse signal can be recovered through the proposed mixed l2/ll-minimization approach in the noiseless case and is stably recovered in the presence of noise and mismodeling error. This improves the result of Eldar and Mishali (in IEEE Trans. Inform. Theory 55: 5302-5316, 2009). We also give another sufficient condition on block RIP for such recovery method: 58 〈 0.307.展开更多
This paper gives new bounds for restricted isometry constant(RIC)in compressed sensing.LetΦbe an m×n real matrix and k be a positive integer with k≤n.The main results of this paper show that if the restricted i...This paper gives new bounds for restricted isometry constant(RIC)in compressed sensing.LetΦbe an m×n real matrix and k be a positive integer with k≤n.The main results of this paper show that if the restricted isometry constant ofΦsat-isfiesδ8ak<1 andδk+ak<3/2−1+√(4a+3)^(2)−8/8aforα>3/8,then k-sparse solution can be recovered exactly via l1 minimization in the noiseless case.In particular,whenα=1,1.5,2 and3,we haveδ2k<0.5746 andδ8k<1,orδ2.5k<0.7046 andδ12k<1,orδ3k<0.7731 andδ16k<1 orδ4k<0.8445 andδ24k<1.展开更多
An algorithmic framework, based on the difference of convex functions algorithm (D- CA), is proposed for minimizing a class of concave sparse metrics for compressed sensing problems. The resulting algorithm iterates...An algorithmic framework, based on the difference of convex functions algorithm (D- CA), is proposed for minimizing a class of concave sparse metrics for compressed sensing problems. The resulting algorithm iterates a sequence ofl1 minimization problems. An exact sparse recovery theory is established to show that the proposed framework always improves on the basis pursuit (l1 minimization) and inherits robustness from it. Numerical examples on success rates of sparse solution recovery illustrate further that, unlike most existing non-convex compressed sensing solvers in the literature, our method always out- performs basis pursuit, no matter how ill-conditioned the measurement matrix is. Moreover, the iterative l1 (ILl) algorithm lead by a wide margin the state-of-the-art algorithms on l1/2 and logarithimic minimizations in the strongly coherent (highly ill-conditioned) regime, despite the same objective functions. Last but not least, in the application of magnetic resonance imaging (MRI), IL1 algorithm easily recovers the phantom image with just 7 line projections.展开更多
In this paper,we study large m asymptotics of the l 1 minimal m-partition problem for the Dirichlet eigenvalue.For any smooth domainΩ⊂R^(n)such that|Ω|=1,we prove that the limit lim_(m→∞)l^(1)_(m)(Ω)=c 0 exists,a...In this paper,we study large m asymptotics of the l 1 minimal m-partition problem for the Dirichlet eigenvalue.For any smooth domainΩ⊂R^(n)such that|Ω|=1,we prove that the limit lim_(m→∞)l^(1)_(m)(Ω)=c 0 exists,and the constant c 0 is independent of the shape ofΩ.Here,l^(1)_(m)(Ω)denotes the minimal value of the normalized sum of the first Laplacian eigenvalues for any m-partition ofΩ.展开更多
In this paper,we investigate truncated l2\l1-2 minimization and its associated alternating direction method of multipliers(ADMM)algorithm for recovering the block sparse signals.Based on the block restricted isometry ...In this paper,we investigate truncated l2\l1-2 minimization and its associated alternating direction method of multipliers(ADMM)algorithm for recovering the block sparse signals.Based on the block restricted isometry property(Block-RIP),a theoretical analysis is presen ted to guarantee the validity of proposed method.Our theore tical resul ts not only show a less error upper bound,but also promote the former recovery condition of truncated l1-2 method for sparse signal recovery.Besides,the algorithm has been compared with some state-of-the-art algorithms and numerical experiments have shown excellent performances on recovering the block sparse signals.展开更多
Based on the range space property (RSP), the equivalent conditions between nonnegative solutions to the partial sparse and the corresponding weighted l1-norm minimization problem are studied in this paper. Different...Based on the range space property (RSP), the equivalent conditions between nonnegative solutions to the partial sparse and the corresponding weighted l1-norm minimization problem are studied in this paper. Different from other conditions based on the spark property, the mutual coherence, the null space property (NSP) and the restricted isometry property (RIP), the RSP- based conditions are easier to be verified. Moreover, the proposed conditions guarantee not only the strong equivalence, but also the equivalence between the two problems. First, according to the foundation of the strict complemenrarity theorem of linear programming, a sufficient and necessary condition, satisfying the RSP of the sensing matrix and the full column rank property of the corresponding sub-matrix, is presented for the unique nonnegative solution to the weighted l1-norm minimization problem. Then, based on this condition, the equivalence conditions between the two problems are proposed. Finally, this paper shows that the matrix with the RSP of order k can guarantee the strong equivalence of the two problems.展开更多
文摘In this work, we consider a homotopic principle for solving large-scale and dense l1underdetermined problems and its applications in image processing and classification. We solve the face recognition problem where the input image contains corrupted and/or lost pixels. The approach involves two steps: first, the incomplete or corrupted image is subject to an inpainting process, and secondly, the restored image is used to carry out the classification or recognition task. Addressing these two steps involves solving large scale l1minimization problems. To that end, we propose to solve a sequence of linear equality constrained multiquadric problems that depends on a regularization parameter that converges to zero. The procedure generates a central path that converges to a point on the solution set of the l1underdetermined problem. In order to solve each subproblem, a conjugate gradient algorithm is formulated. When noise is present in the model, inexact directions are taken so that an approximate solution is computed faster. This prevents the ill conditioning produced when the conjugate gradient is required to iterate until a zero residual is attained.
文摘A joint two-dimensional(2D)direction-of-arrival(DOA)and radial Doppler frequency estimation method for the L-shaped array is proposed in this paper based on the compressive sensing(CS)framework.Revised from the conventional CS-based methods where the joint spatial-temporal parameters are characterized in one large scale matrix,three smaller scale matrices with independent azimuth,elevation and Doppler frequency are introduced adopting a separable observation model.Afterwards,the estimation is achieved by L1-norm minimization and the Bayesian CS algorithm.In addition,under the L-shaped array topology,the azimuth and elevation are separated yet coupled to the same radial Doppler frequency.Hence,the pair matching problem is solved with the aid of the radial Doppler frequency.Finally,numerical simulations corroborate the feasibility and validity of the proposed algorithm.
基金supported by the Aviation Science Foundation(201120M5007)the Natural Science Foundation of Beijing(4102050)
文摘This paper presents a new method for image separation through employing a combined dictionary consisting of wavelets and complex shearlets. Because the combined dictionary sparsely represents points and curvilinear singularities respectively, the image can be decomposed into pointlike and curvelike parts as accurate as possible. The proposed method based on the geo- metric separation theory introduced by Donoho in 2005 shows that accurate geometric separation of the morphologically distinct fea- tures of points and curves can be achieved by l1 minimization. The experimental results show that the proposed method can not only be effective but also greatly reduce the computing time.
文摘Radial functions have become a useful tool in numerical mathematics. On the sphere they have to be identified with the zonal functions. We investigate zonal polynomials with mass concentration at the pole, in the sense of their L1-norm is attaining the minimum value. Such polynomials satisfy a complicated system of nonlinear e-quations (algebraic if the space dimension is odd, only) and also a singular differential equation of third order. The exact order of decay of the minimum value with respect to the polynomial degree is determined. By our results we can prove that some nodal systems on the sphere, which are defined by a minimum-property, are providing fundamental matrices which are diagonal-dominant or bounded with respect to the ∞-norm, at least, as the polynomial degree tends to infinity.
基金Supported by National Natural Science Foundation of China (Grant Nos. 11171299 and 91130009)Natural Science Foundation of Zhejiang Province of China (Grant No. Y6090091)
文摘We consider efficient methods for the recovery of block sparse signals from underdetermined system of linear equations. We show that if the measurement matrix satisfies the block RIP with δ2s 〈 0.4931, then every block s-sparse signal can be recovered through the proposed mixed l2/ll-minimization approach in the noiseless case and is stably recovered in the presence of noise and mismodeling error. This improves the result of Eldar and Mishali (in IEEE Trans. Inform. Theory 55: 5302-5316, 2009). We also give another sufficient condition on block RIP for such recovery method: 58 〈 0.307.
基金This work was partially supported by the National Basic Research Program of China(No.2010CB732501)the National Natural Science Foundation of China(No.11171018)+1 种基金d the Fundamental Research Funds for the Central Universities(No.2013JBM095)We thank the two anonymous referees for their very useful comments.
文摘This paper gives new bounds for restricted isometry constant(RIC)in compressed sensing.LetΦbe an m×n real matrix and k be a positive integer with k≤n.The main results of this paper show that if the restricted isometry constant ofΦsat-isfiesδ8ak<1 andδk+ak<3/2−1+√(4a+3)^(2)−8/8aforα>3/8,then k-sparse solution can be recovered exactly via l1 minimization in the noiseless case.In particular,whenα=1,1.5,2 and3,we haveδ2k<0.5746 andδ8k<1,orδ2.5k<0.7046 andδ12k<1,orδ3k<0.7731 andδ16k<1 orδ4k<0.8445 andδ24k<1.
文摘An algorithmic framework, based on the difference of convex functions algorithm (D- CA), is proposed for minimizing a class of concave sparse metrics for compressed sensing problems. The resulting algorithm iterates a sequence ofl1 minimization problems. An exact sparse recovery theory is established to show that the proposed framework always improves on the basis pursuit (l1 minimization) and inherits robustness from it. Numerical examples on success rates of sparse solution recovery illustrate further that, unlike most existing non-convex compressed sensing solvers in the literature, our method always out- performs basis pursuit, no matter how ill-conditioned the measurement matrix is. Moreover, the iterative l1 (ILl) algorithm lead by a wide margin the state-of-the-art algorithms on l1/2 and logarithimic minimizations in the strongly coherent (highly ill-conditioned) regime, despite the same objective functions. Last but not least, in the application of magnetic resonance imaging (MRI), IL1 algorithm easily recovers the phantom image with just 7 line projections.
基金supported by National Science Foundation of USA(Grant Nos.DMS1501000 and DMS-1955249)。
文摘In this paper,we study large m asymptotics of the l 1 minimal m-partition problem for the Dirichlet eigenvalue.For any smooth domainΩ⊂R^(n)such that|Ω|=1,we prove that the limit lim_(m→∞)l^(1)_(m)(Ω)=c 0 exists,and the constant c 0 is independent of the shape ofΩ.Here,l^(1)_(m)(Ω)denotes the minimal value of the normalized sum of the first Laplacian eigenvalues for any m-partition ofΩ.
基金The authors would like to thank reviewers for valuable comments.This work was supported by Natural Science Foundation of China(Grant Nos.61673015,61273020)Fundamental Research Funds for the Central Universities(Grant Nos.XDJK2015A007,XDJK 2018C076,SWU1809002).
文摘In this paper,we investigate truncated l2\l1-2 minimization and its associated alternating direction method of multipliers(ADMM)algorithm for recovering the block sparse signals.Based on the block restricted isometry property(Block-RIP),a theoretical analysis is presen ted to guarantee the validity of proposed method.Our theore tical resul ts not only show a less error upper bound,but also promote the former recovery condition of truncated l1-2 method for sparse signal recovery.Besides,the algorithm has been compared with some state-of-the-art algorithms and numerical experiments have shown excellent performances on recovering the block sparse signals.
基金Research supported by the National Natural Science Foundation of China under Grant 61672005
文摘Based on the range space property (RSP), the equivalent conditions between nonnegative solutions to the partial sparse and the corresponding weighted l1-norm minimization problem are studied in this paper. Different from other conditions based on the spark property, the mutual coherence, the null space property (NSP) and the restricted isometry property (RIP), the RSP- based conditions are easier to be verified. Moreover, the proposed conditions guarantee not only the strong equivalence, but also the equivalence between the two problems. First, according to the foundation of the strict complemenrarity theorem of linear programming, a sufficient and necessary condition, satisfying the RSP of the sensing matrix and the full column rank property of the corresponding sub-matrix, is presented for the unique nonnegative solution to the weighted l1-norm minimization problem. Then, based on this condition, the equivalence conditions between the two problems are proposed. Finally, this paper shows that the matrix with the RSP of order k can guarantee the strong equivalence of the two problems.