Large-scale sparse recovery (SR) by solving -norm relaxations over Big Dictionary is a very challenging task. Plenty of greedy methods have therefore been proposed to address big SR problems, but most of them require restricted conditions for the convergence. Moreover, it is non-trivial for them to incorporate the -norm regularization that is required for robust signal recovery. We address these issues in this paper by proposing a Matching Pursuit LASSO (MPL) algorithm, based on a novel quadratically constrained linear program (QCLP) formulation, which has several advantages over existing methods. Firstly, it is guaranteed to converge to a global solution. Secondly, it greatly reduces the computation cost of the -norm methods over Big Dictionaries. Lastly, the exact sparse recovery condition of MPL is also investigated.Index Terms-Sparse recovery, compressive sensing, LASSO, matching pursuit, big dictionary, convex programming.(4) where denotes the regularization parameter.Many algorithms have been proposed to address the nonsmooth LASSO problem. An interior-point method is proposed in [33], and a gradient projection (GPSR) method was subsequently proposed by transforming LASSO into a quadratic programming problem [27]. A proximal gradient (PG) method is proposed to solve (4) based on a shrinkage-threshold operator [34]. To speed up the PG, a fast iterative shrinkage-threshold algorithm (FISTA) (also known as the accelerated proximal gradient method) is proposed in [29]. Recently, to tackle large-scale 1053-587X