This paper provides a necessary and sufficient condition for guaranteeing exponential stability of the linear difference equationwhere a > 0, b > 0 are constants and A, B are n × n square matrices, in terms of a linear matrix inequality (LMI) of size (k + 1) n × (k + 1) n where k ≥ 1 is some integer. Different from an existing condition where the coefficients (A, B) appear as highly nonlinear functions, the proposed LMI condition involves matrices that are linear functions of (A, B) . Such a property is further used to deal with the robust stability problem in case of norm bounded uncertainty and polytopic uncertainty, and the state feedback stabilization problem. Solutions to these two problems are expressed by LMIs. A time domain interpretation of the proposed LMI condition in terms of Lyapunov-Krasovskii functional is given, which helps to reveal the relationships among the existing methods. Numerical example demonstrates the effectiveness of the proposed method.