08.10.2016 Views

Foundations of Data Science

2dLYwbK

2dLYwbK

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Suppose ˜x 0 were another minimum. Then ∇ is also a subgradient at ˜x 0 as it is at x 0 .<br />

To see this, for ∆x such that A∆x = 0 ,<br />

‖˜x 0 + ∆x‖ 1<br />

=<br />

∥ x 0 + ˜x 0 − x 0 + ∆x<br />

} {{ }<br />

≥ ‖x 0 ‖ 1<br />

+ ∇ T (˜x 0 − x 0 + ∆x) .<br />

∥<br />

1<br />

α<br />

The above equation follows from the definition <strong>of</strong> ∇ being a subgradient for the one norm<br />

function, ‖‖ 1<br />

, at x 0 . Thus,<br />

But<br />

‖˜x 0 + ∆x‖ 1<br />

≥ ‖x 0 ‖ 1<br />

+ ∇ T (˜x 0 − x 0 ) + ∇ T ∆x.<br />

∇ T (˜x 0 − x 0 ) = w T A (˜x 0 − x 0 ) = w T (b − b) = 0.<br />

Hence, since ˜x 0 being a minimum means || ˜x 0 || 1 = ||x 0 || 1 ,<br />

‖˜x 0 + ∆x‖ 1<br />

≥ ‖x 0 ‖ 1<br />

+ ∇ T ∆x = || ˜x 0 || 1 + ∇ T ∆x.<br />

This implies that ∇ is a sub gradient at ˜x 0 .<br />

Now, ∇ is a subgradient at both x 0 and ˜x 0 . By Proposition 10.2, we must have that<br />

(∇) i = sgn((x 0 ) i ) = sgn((˜x 0 ) i ), whenever either is nonzero and |(∇) i | < 1, whenever either<br />

is 0. It follows that x 0 and ˜x 0 have the same sparseness pattern. Since Ax 0 = b and<br />

A˜x 0 = b and x 0 and ˜x 0 are both nonzero on the same coordinates, and by the assumption<br />

that the columns <strong>of</strong> A corresponding to the nonzeros <strong>of</strong> x 0 and ˜x 0 are independent, it<br />

must be that x 0 = ˜x 0 .<br />

10.3.3 Restricted Isometry Property<br />

Next we introduce the restricted isometry property that plays a key role in exact<br />

reconstruction <strong>of</strong> sparse vectors. A matrix A satisfies the restricted isometry property,<br />

RIP, if for any s-sparse x there exists a δ s such that<br />

(1 − δ s ) |x| 2 ≤ |Ax| 2 ≤ (1 + δ s ) |x| 2 . (10.1)<br />

Isometry is a mathematical concept; it refers to linear transformations that exactly preserve<br />

length such as rotations. If A is an n × n isometry, all its eigenvalues are ±1 and<br />

it represents a coordinate system. Since a pair <strong>of</strong> orthogonal vectors are orthogonal in all<br />

coordinate system, for an isometry A and two orthogonal vectors x and y, x T A T Ay = 0.<br />

We will prove approximate versions <strong>of</strong> these properties for matrices A satisfying the restricted<br />

isometry property. The approximate versions will be used in the sequel.<br />

A piece <strong>of</strong> notation will be useful. For a subset S <strong>of</strong> columns <strong>of</strong> A, let A S denote the<br />

submatrix <strong>of</strong> A consisting <strong>of</strong> the columns <strong>of</strong> S.<br />

Lemma 10.4 If A satisfies the restricted isometry property, then<br />

337

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!