Step 3

In this step we will build a slimpy script to solve the problem:

$\tilde{\textbf{x}} = \arg\min \|\textbf{x}\|_1\quad\mbox{s.t.}\quad \textbf{C} \textbf{x}+\epsilon=\textbf{y}$

where C is the cosine transform and x is the swell noise and $\epsilon$ is the signal. I will use the existing thresholded landweber solver with a linear threshold cooling scheme.

Objectives:

Use an exsiting iterative solver

Walkthrough

Using the code from the previous example I have chenged the code in the body of the function SwellSeparation.

>>> thresh = LinearCooling( lmax, lmin, nouter )
    
>>> solver = GenThreshLandweber( nouter, ninner, thresh  )
    
>>> x = solver.solve( C.H, y )

>>> res = C.H * x