A real m × n matrix A and a vector y ∈ R m determine the discrete l 1 -regularization (DLR) problemwhere | · | 1 denotes the l 1 -norm of a vector and ρ ≥ 0 is a nonnegative parameter. In this paper, we provide a detailed analysis of this problem which include a characterization of all solutions to (0.1), remarks about the geometry of the solution set and an effective iterative algorithm for numerical solution of (0.1).We are specially interested in the behavior of the solution of (0.1) as a function 442 C.A. Micchelli, T. Zhao of ρ and in this regard, we prove in general the existence of critical values of ρ between which the l 1 -norm of any solution remains constant. These general remarks are significantly refined when A is a strictly totally positive (STP) matrix. The importance of STP matrices is well-established [5,14]. Under this setting, the relationship between the number of nonzero coordinates of a distinguished solution of the DLR problem is precisely explained as a function of the regularization parameter for a certain class of vectors in R m . Throughout our analysis of the DLR problem, we emphasize the importance of the dual maximum problem by demonstrating that any solution of it leads to a solution of the DLR problem, and vice versa.