**ABSTRACT**

Wavelets are regarded by many as primarily a new subject in pure and applied mathematics. Perhaps one of the most common applications of wavelet is in signal processing. In this study we obtain an algorithm to analyze and synthesize a signal or two-dimensional function s = f(t) by using two-dimensional wavelet method. We consider a sample point (t

_{i,j},s

_{i,j}) includes a value s

_{i,j}= f(t

_{i,j}) at height s

_{i,j}and abscissa (time or location) t

_{i,j}. We propose, obtaining an algorithm of two-dimensional wavelet decomposition to estimate a function by using MATLAB software for computing wavelet coefficients. Some relationships between wavelet coefficients are investigated.

PDF Abstract XML References Citation

**Received:**July 19, 2010;

**Accepted:**December 08, 2010;

**Published:**February 22, 2011

####
**How to cite this article**

*Asian Journal of Applied Sciences, 4: 414-422.*

**DOI:**10.3923/ajaps.2011.414.422

**URL:**https://scialert.net/abstract/?doi=ajaps.2011.414.422

**INTRODUCTION**

Wavelets are regarded by many as primarily a new subject in pure and applied math-ematics. Indeed, many papers published on wavelets contain esoteric-looking theorems with complicated proofs. Wavelet analysis was led by Daubechies (1988). Many colleagues contributed in different ways: Meyer (1990), Walter (1993), Vidakovic (1999), Cohen *et al*. (1993), Doosti *et al*. 2008, Afshari (2008), Antoniadis *et al*. (1994) and Clyed *et al*. (1998). Perhaps one of the most common applications of wavelets is in signal processing. A signal is a sequence of numerical measurements, typically obtained electronically.

To analyze and synthesize a signal- which can be any array of data- in terms of simple wavelets, we employ shifts and dilation of mathematical function, but do not involve either calculus or linear algebra.

The first step in applying wavelets to any signal consists in representing the signal under consideration by a mathematical function f. For example, a sound, the values s = f(t) measure the sound at each time t at a fixed location.

The first step in the analysis of a one-dimensional signal with wavelets consists in approximating its function by means of sample alone. One of the simplest methods of approximation uses a horizontal stair step extended through each sample point. The resulting steps form a new function denotes by and called a step function, which approximates the sampled function s = f(t). The analysis of the approximating function f in terms of wavelets requires a precise labeling of each step. By means of shifts and dilations of the basic unit step function, denoted by φ(0.1).

If a sample point (t_{i},s_{j}) includes a value s_{j} = f(t_{j}) at height s_{j} and abscissa (time or location) t_{j}, then the sample point corresponds to the step function- which approximates f at height s_{j} on the interval [t_{j},t_{(j+1)}] , where denotes the indicator function of set [t_{j},t_{(j+1)}].

Adding all step functions approximating corresponding to all the points in the sample yields the simple step function below:

(1) |

To analyze a signal or function in term of wavelets, the fast Haar wavelet transform begins with initialization of an array with 2^{n} entries, and then proceeds with n iterations of the basic transform explained in Eq. 1.

For each index jε{1,2,....,n}, before iteration number j, the array will consist of 2^{n–(j–1)} coefficients of 2^{n–(j–1)} step function φ_{k,n–(j–1)} defined below. After iteration number j, the array will consist of half as many 2^{n–j} coefficient of 2^{n–j} step function φ_{k,n–j} and 2^{n–j} coefficient ψ_{k,n–j}, such as:

(2) |

(3) |

**MAIN RESULTS**

**Two dimensional wavelets algorithm: **For any function fεL^{2}(R) we can write a formal expansion (Daubechies, 1992):

(4) |

Here φ(x) and ψ(x) are the scale function and the orthogonal wavelet, respectively as the following:

(5) |

Constitute an (inhomogeneous) orthonormal basis of L^{2}(R). It is clear that for Haar wavelet:

(6) |

So, we can write:

(7) |

Not that φεV_{0} therefore φεV_{1} , because V_{0}⊂V_{1}.

Since {φ_{l,k},(x), kεZ is an orthonorml basis for V_{1}, there exists a sequence b_{k} such that:

(8) |

Know we consider a two-dimensional function f(x,y)εL^{2}(R^{2}). For approximating this function we using two dimensional of wavelets and obtaining an algorithm for estimation of function f with computation of wavelet coefficients. For two dimensional wavelets, encodings can consist of matrices, indexed by rows and columns. To expand this argument, we need some definition as the following:

• | Definition 1: Let f(x,y)εL^{2}(R^{2}) denote the two dimensional signal. We said ψ_{j1},_{j2},_{k1},_{k2} = ψ_{j1},_{k1}(x).ψ_{j2},_{k2} (x) is orthogonal wavelet basis |

• | Definition 2: For each pair functions f and g, the tensor product of f and g is the function denoted by f⊗ g as follows: |

(9) |

• | Definition 3: Suppose that function space, V_{j},εz to be {V_{j} = f,gεL^{2}(R^{2}):f,g are piecewise constant on [k2^{–j},(k+1)2^{–j}],kεz} = V_{j}⊗ V_{j} |

If this sequence of subspaces possesses has the following properties:

• | –...⊂V_{–2} ⊂⊂V_{–1} ⊂V_{–0} ⊂V_{–1} ⊂V... |

• |

• | –f(x,y)εV_{j} ⇔f(2x,2y)εV_{j+1} |

• | f(x,y)εV_{0}⇒f(x–k_{1},y–k_{2})ε V_{0}∀k_{1},k_{2}εZ |

• | There exists a function φ(x,y)=φ(x).φ(y)εV_{0} such that the set {φ_{j,k1,k2}(x,y)=2^{j}φ(2^{j}x–k_{1},2^{j}y–k_{2}),j,k_{1},k_{2}εz} constitutes an orthonormal basis for V_{0}, then we said (V_{j})_{jεZ} to be form a multiresolution analysis (MRA) of L^{2}(R^{2}). |

**Remark 1:** Suppose that function space, w_{j},jεz to be then we can write as the following:

• | |

• | V_{j+1} = V_{j+1} ⊗ V_{j+1} = (V_{j}⊕ W_{j}) ⊗ (V_{j} ⊕ W_{j})=(V _{j} ⊗ V_{j}) ⊕ (V_{j} ⊗ W_{j})⊕ (W_{j} ⊗ V_{j})⊕ (W_{j} ⊗ W_{j}) = (V_{j} ⊕ W_{j}) |

**Theorem 1:** In Definition 1, The wavelet space {W_{j}, jεZ} and scale space {(V_{j}, jεZ} are mutually orthogonal.

**Proof:** First we prove that the scaling function and mother wavelet are orthogonal.

The last step follows since the summand for k is the opposite of the summand for k–1–2*l*, so each term is negated.

It can be seen similarly that each integer translation of the mother wavelet ψ is also orthogonal to φ as the following:

The last step follows, since the summand for k is the opposite of the summand k–1–2*l*, so each term is negated, because of the square sumability of the sequence b_{k}.

Straightforward extension of this argument will show that ψ_{j,k}⊥φ_{j,k} for all k,*l*εZ and complete the proof.

**Remark 2:** Each basic square-step function has value l in a selected square and 0 everywhere as the following:

(10) |

Similarly we can write:

(11) |

(12) |

**Lemmal:**

**Proof:** It is easy to see that ,

(13) |

By substituting φ_{(0,1)} and ψ_{(0,1)} instead of f and g in Eq. 9 and changing the place of them, we can write the tensor product of as follows:

(14) |

(15) |

(16) |

**Remark 3:** The subscript h and v and d indicates the correspondence of such wavelets with horizontal and vertical, and diagonal changes in the data, since the detail spaces tend to emphasize coefficient-cliques describing horizontal and vertical and diagonal futures the image. The spaces are spanned by translation of:

(17) |

(18) |

For nonnegative integers j which is denoting the frequency and k,*l* which are denoting the location.

**Algorithm:** We consider a function f with sample values approximated by a square-function and denote by a matrix:

We can write as follows:

(19) |

Algorithm begins with one dimensional wavelet transform as described in algorithm, for each row as follows:

We can repeat this algorithm and estimate function f.

**EXAMPLE**

**Example 1:** let for approximating f we chose sample values:

The square-step approximation is:

(20) |

Consider again the above data. By the algorithm, we can write:

The one-dimensional Fast wavelet transform extends to two- dimensional fast wavelet transform with tensor products, through alternating applications of the one-dimensional transform to each row and then to each new column.

So by the tensor product wavelet,

(21) |

**Example 2:** supposed that a function f sampled at 4x4 matrix values

On a square grid, and we consider approximating f, by square-step function as following:

The intermediate result just obtained corresponds to one basic two dimensional wavelet transform on each of the four adjacent 2x2 square matrixes.

**RESULTS**

• | All of the sums from the upper left hand corners of the four 2x2 blocks go into the upper left hand corner of the larger grade, corresponding to wavelets of the form φ_{m,n} ⊗ φ_{k,l}. |

• | All of the difference from the upper right-hand corners of the four 2x2 blocks go into the upper right-hand corner of the larger grade, corresponding to wavelets of the form φ_{m,n} ⊗ ψ_{k,l}. |

• | All of the difference from the lower right-hand corners of the four 2x2 blocks go into the lower right-hand corner of the larger grade, corresponding to wavelets of the form φ_{m,n} ⊗ φ_{k,l}. |

• | All of the difference from the lower left-hand corners of the four 2x2 blocks go into the lower left-hand corner of the larger grade, corresponding to wavelets of the form φ_{m,n} ⊗ ψ_{k,l}. |

**Proposition:** We can complete the two dimensional transform by the following method:

Finally, perform a two dimensional wavelet transform only the four entries iv the upper left-hand corner as follows:

Thus the matrix:

**ACKNOWLEDGMENT**

The support of Research Committee of Persian Gulf University is greatly acknowledged.

####
**REFERENCES**

- Doosti, H., M. Afshari and H.A. Niroumand, 2008. Wavelets for nonparametric stochastic regression with mixing stochastic process. Commun. Stat. Theory Methods, 37: 373-385.

Direct Link - Antoniadis, A., G. Gregoire and I. McKeague, 1994. Wavelet methods for curve estimation. J. Am. Stat. Assoc., 89: 1340-1353.

Direct Link - Clyed, M.A., G. Parmigiana and B. Vidakovic, 1998. Multiple shrinkage and subset selection in wavelets. Biometrika, 85: 391-401.

Direct Link - Cohen, A., I. Daubechies and P. Vial, 1993. Wavelets on the interval and fast wavelet transform. Applied Comput. Harmon. Anal., 1: 54-81.

Direct Link