<< Chapter < Page Chapter >> Page >
Linear versus nonlinear approximation in the wavelet domain. (a) Linear approximation of the Cameraman test image obtained by keeping the K = 4096 lowest-frequency wavelet coefficients from the five-level wavelet decomposition. The MSE with respect to the original image is 353 . (b) Nonlinear approximation of the Cameraman test image obtained by keeping the K = 4096 largest wavelet coefficients from the five-level wavelet decomposition. The MSE with respect to the original image is 72 . Compared with linear approximation, more high frequency coefficients are included, which allows better representation of features such as edges.

Nonlinear approximation

A related question often arises in settings involving signal dictionaries. Rather than finding the best approximation to asignal f using a fixed collection of K elements from the dictionary Ψ , one may often seek the best K -term representation to f among all possible expansions that use K terms from the dictionary. Compared to linear approximation, this type of nonlinearapproximation  [link] , [link] utilizes the ability of the dictionary to adapt: different elements may be important forrepresenting different signals.

The K -term nonlinear approximation problem corresponds to the optimization

s K , p * : = arg min s Σ K s - f p .
(For the sake of generality, we consider general L p and p norms in this section.) Due to the nonlinearity of the set Σ K for a given dictionary, solving this problem can be difficult. Supposing Ψ is an orthonormal basis and p = 2 , the solution to [link] is easily obtained by thresholding: one simply computes the coefficients α and keeps the K largest (setting the remaining coefficients to zero). The approximation error is then given simplyby the coefficients that are discarded:
s K , 2 * - f 2 = k > K α ˜ k 2 1 / 2 .
When Ψ is a redundant dictionary, however, the situation is much more complicated. We mention more on this below (see also [link] (b)).

Measuring approximation quality

One common measure for the quality of a dictionary Ψ in approximating a signal class is the fidelity of its K -term representations. Often one examines theasymptotic rate of decay of the K -term approximation error as K grows large. Defining

σ K ( f ) p : = s K , p * - f p ,
for a given signal f we may consider the asymptotic decay of σ K ( f ) p as K . (We recall the dependence of [link] and hence [link] on the dictionary Ψ .) In many cases, the function σ K ( f ) p will decay as K - r for some r , and when Ψ represents a harmonic dictionary, faster decay rates tend to correspond to smoother functions. Indeed, onecan show that when Ψ is an orthonormal basis, then σ K ( f ) 2 will decay as K - r if and only if α ˜ k decays as k - r + 1 / 2 [link] .

Nonlinear approximation of piecewise smooth functions

Let f C H be a 1-D function. Supposing the wavelet dictionary has more than H vanishing moments, then f can be well approximated using its K largest coefficients (most of which are at coarse scales). As K grows large, the nonlinear approximation error will decay We use the notation f ( α ) g ( α ) , or f ( α ) = O ( g ( α ) ) , if there exists a constant C , possibly large but not dependent on the argument α , such that f ( α ) C g ( α ) . as σ K ( f ) 2 K - H .

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Concise signal models. OpenStax CNX. Sep 14, 2009 Download for free at http://cnx.org/content/col10635/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Concise signal models' conversation and receive update notifications?

Ask