# 5.2 Greedy algorithms

 Page 1 / 3
In this module we provide an overview of some of the most common greedy algorithms and their application to the problem of sparse recovery.

## Setup

As opposed to solving a (possibly computationally expensive) convex optimization program, an alternate flavor to sparse recovery is to apply methods of sparse approximation . Recall that the goal of sparse recovery is to recover the sparsest vector $x$ which explains the linear measurements $y$ . In other words, we aim to solve the (nonconvex) problem:

$\underset{\mathcal{I}}{min}\left\{|\mathcal{I}|:y,=,\sum _{i\in \mathcal{I}},{\phi }_{i},{x}_{i}\right\},$

where $\mathcal{I}$ denotes a particular subset of the indices $i=1,...,N$ , and ${\phi }_{i}$ denotes the ${i}^{\mathrm{th}}$ column of $\Phi$ . It is well known that searching over the power set formed by the columns of $\Phi$ for the optimal subset ${\mathcal{I}}^{*}$ with smallest cardinality is NP-hard. Instead, classical sparse approximation methods tackle this problem by greedily selecting columns of $\Phi$ and forming successively better approximations to $y$ .

## Matching pursuit

Matching Pursuit (MP), named and introduced to the signal processing community by Mallat and Zhang  [link] , [link] , is an iterative greedy algorithm that decomposes a signal into a linear combination of elements from a dictionary. In sparse recovery, this dictionary is merely the sampling matrix $\Phi \in {\mathbb{R}}^{M×N}$ ; we seek a sparse representation ( $x$ ) of our “signal” $y$ .

MP is conceptually very simple. A key quantity in MP is the residual $r\in {\mathbb{R}}^{M}$ ; the residual represents the as-yet “unexplained” portion of the measurements. At each iteration of the algorithm, we select a vector from the dictionary that is maximally correlated with the residual $r$ :

${\lambda }_{k}=arg\underset{\lambda }{max}\frac{⟨{r}_{k},{\phi }_{\lambda }⟩{\phi }_{\lambda }}{\parallel {\phi }_{\lambda }{\parallel }^{2}}.$

Once this column is selected, we possess a “better” representation of the signal, since a new coefficient indexed by ${\lambda }_{k}$ has been added to our signal approximation. Thus, we update both the residual and the approximation as follows:

$\begin{array}{cc}\hfill {r}_{k}& ={r}_{k-1}-\frac{⟨{r}_{k-1},{\phi }_{{\lambda }_{k}}⟩{\phi }_{{\lambda }_{k}}}{\parallel {\phi }_{{\lambda }_{k}}{\parallel }^{2}},\hfill \\ \hfill {\stackrel{^}{x}}_{{\lambda }_{k}}& ={\stackrel{^}{x}}_{{\lambda }_{k}}+⟨{r}_{k-1},{\phi }_{{\lambda }_{k}}⟩.\hfill \end{array}$

and repeat the iteration. A suitable stopping criterion is when the norm of $r$ becomes smaller than some quantity. MP is described in pseudocode form below.

Inputs: Measurement matrix $\Phi$ , signal measurements $y$ Outputs: Sparse signal $\stackrel{^}{x}$ initialize: ${\stackrel{^}{x}}_{0}=0$ , $r=y$ , $i=0$ . while ħalting criterion false do 1. $i←i+1$ 2. $b←{\Phi }^{T}r$ {form residual signal estimate} 3. ${\stackrel{^}{x}}_{i}←{\stackrel{^}{x}}_{i-1}+\mathbf{T}\left(1\right)$ {update largest magnitude coefficient} 4. $r←r-\Phi {\stackrel{^}{x}}_{i}$ {update measurement residual} end while return $\stackrel{^}{x}←{\stackrel{^}{x}}_{i}$ 

Although MP is intuitive and can find an accurate approximation of the signal, it possesses two major drawbacks: (i) it offers no guarantees in terms of recovery error; indeed, it does not exploit the special structure present in the dictionary $\Phi$ ; (ii) the required number of iterations required can be quite large. The complexity of MP is $O\left(MNT\right)$   [link] , where $T$ is the number of MP iterations

## Orthogonal matching pursuit (omp)

Matching Pursuit (MP) can prove to be computationally infeasible for many problems, since the complexity of MP grows linearly in the number of iterations $T$ . By employing a simple modification of MP, the maximum number of MP iterations can be upper bounded as follows. At any iteration $k$ , Instead of subtracting the contribution of the dictionary element with which the residual $r$ is maximally correlated, we compute the projection of $r$ onto the orthogonal subspace to the linear span of the currently selected dictionary elements. This quantity thus better represents the “unexplained” portion of the residual, and is subtracted from $r$ to form a new residual, and the process is repeated. If ${\Phi }_{\Omega }$ is the submatrix formed by the columns of $\Phi$ selected at time step $t$ , the following operations are performed:

what is math number
x-2y+3z=-3 2x-y+z=7 -x+3y-z=6
Need help solving this problem (2/7)^-2
x+2y-z=7
Sidiki
what is the coefficient of -4×
-1
Shedrak
the operation * is x * y =x + y/ 1+(x × y) show if the operation is commutative if x × y is not equal to -1
An investment account was opened with an initial deposit of $9,600 and earns 7.4% interest, compounded continuously. How much will the account be worth after 15 years? Kala Reply lim x to infinity e^1-e^-1/log(1+x) given eccentricity and a point find the equiation Moses Reply 12, 17, 22.... 25th term Alexandra Reply 12, 17, 22.... 25th term Akash College algebra is really hard? Shirleen Reply Absolutely, for me. My problems with math started in First grade...involving a nun Sister Anastasia, bad vision, talking & getting expelled from Catholic school. When it comes to math I just can't focus and all I can hear is our family silverware banging and clanging on the pink Formica table. Carole I'm 13 and I understand it great AJ I am 1 year old but I can do it! 1+1=2 proof very hard for me though. Atone hi Adu Not really they are just easy concepts which can be understood if you have great basics. I am 14 I understood them easily. Vedant hi vedant can u help me with some assignments Solomon find the 15th term of the geometric sequince whose first is 18 and last term of 387 Jerwin Reply I know this work salma The given of f(x=x-2. then what is the value of this f(3) 5f(x+1) virgelyn Reply hmm well what is the answer Abhi If f(x) = x-2 then, f(3) when 5f(x+1) 5((3-2)+1) 5(1+1) 5(2) 10 Augustine how do they get the third part x = (32)5/4 kinnecy Reply make 5/4 into a mixed number, make that a decimal, and then multiply 32 by the decimal 5/4 turns out to be AJ how Sheref can someone help me with some logarithmic and exponential equations. Jeffrey Reply sure. what is your question? ninjadapaul 20/(×-6^2) Salomon okay, so you have 6 raised to the power of 2. what is that part of your answer ninjadapaul I don't understand what the A with approx sign and the boxed x mean ninjadapaul it think it's written 20/(X-6)^2 so it's 20 divided by X-6 squared Salomon I'm not sure why it wrote it the other way Salomon I got X =-6 Salomon ok. so take the square root of both sides, now you have plus or minus the square root of 20= x-6 ninjadapaul oops. ignore that. ninjadapaul so you not have an equal sign anywhere in the original equation? ninjadapaul hmm Abhi is it a question of log Abhi 🤔. Abhi I rally confuse this number And equations too I need exactly help salma But this is not salma it's Faiza live in lousvile Ky I garbage this so I am going collage with JCTC that the of the collage thank you my friends salma Commplementary angles Idrissa Reply hello Sherica im all ears I need to learn Sherica right! what he said ⤴⤴⤴ Tamia hii Uday hi salma hi Ayuba Hello opoku hi Ali greetings from Iran Ali salut. from Algeria Bach hi Nharnhar A soccer field is a rectangle 130 meters wide and 110 meters long. The coach asks players to run from one corner to the other corner diagonally across. What is that distance, to the nearest tenths place. Kimberly Reply Jeannette has$5 and \$10 bills in her wallet. The number of fives is three more than six times the number of tens. Let t represent the number of tens. Write an expression for the number of fives.
What is the expressiin for seven less than four times the number of nickels
How do i figure this problem out.
how do you translate this in Algebraic Expressions
why surface tension is zero at critical temperature
Shanjida
I think if critical temperature denote high temperature then a liquid stats boils that time the water stats to evaporate so some moles of h2o to up and due to high temp the bonding break they have low density so it can be a reason
s.
Need to simplify the expresin. 3/7 (x+y)-1/7 (x-1)=
. After 3 months on a diet, Lisa had lost 12% of her original weight. She lost 21 pounds. What was Lisa's original weight?
how did you get the value of 2000N.What calculations are needed to arrive at it
Privacy Information Security Software Version 1.1a
Good
Got questions? Join the online conversation and get instant answers!