Walsh function

From HandWiki
Natural ordered Hadamard matrix (middle matrix) of order 16 that is sequency ordered to output a Walsh matrix (right matrix).
Both contain the 16 Walsh functions of order 16 as rows (and columns).
In the right matrix, the number of sign changes per row is consecutive.

In mathematics, more specifically in harmonic analysis, Walsh functions form a complete orthogonal set of functions that can be used to represent any discrete function—just like trigonometric functions can be used to represent any continuous function in Fourier analysis.[1] They can thus be viewed as a discrete, digital counterpart of the continuous, analog system of trigonometric functions on the unit interval. But unlike the sine and cosine functions, which are continuous, Walsh functions are piecewise constant. They take the values −1 and +1 only, on sub-intervals defined by dyadic fractions.

The system of Walsh functions is known as the Walsh system. It is an extension of the Rademacher system of orthogonal functions.[2]

Walsh functions, the Walsh system, the Walsh series,[3] and the fast Walsh–Hadamard transform are all named after the American mathematician Joseph L. Walsh. They find various applications in physics and engineering when analyzing digital signals.

Historically, various numerations of Walsh functions have been used; none of them is particularly superior to another. This articles uses the Walsh–Paley numeration.

Definition

We define the sequence of Walsh functions [math]\displaystyle{ W_k : [0,1] \rightarrow \{-1,1\} }[/math], [math]\displaystyle{ k\in\mathbb{N} }[/math] as follows.

For any natural number k, and real number [math]\displaystyle{ x \in [0,1] }[/math], let

[math]\displaystyle{ k_j }[/math] be the jth bit in the binary representation of k, starting with [math]\displaystyle{ k_0 }[/math] as the least significant bit, and
[math]\displaystyle{ x_j }[/math] be the jth bit in the fractional binary representation of [math]\displaystyle{ x }[/math], starting with [math]\displaystyle{ x_1 }[/math] as the most significant fractional bit.

Then, by definition

[math]\displaystyle{ W_k(x) = (-1)^{\sum_{j=0}^\infty k_jx_{j+1}} }[/math]

In particular, [math]\displaystyle{ W_0(x) = 1 }[/math] everywhere on the interval, since all bits of k are zero.

Notice that [math]\displaystyle{ W_{2^m} }[/math] is precisely the Rademacher function rm. Thus, the Rademacher system is a subsystem of the Walsh system. Moreover, every Walsh function is a product of Rademacher functions:

[math]\displaystyle{ W_k(x) = \prod_{j=0}^\infty r_j(x)^{k_j} }[/math]

Comparison between Walsh functions and trigonometric functions

Walsh functions and trigonometric functions are both systems that form a complete, orthonormal set of functions, an orthonormal basis in the Hilbert space [math]\displaystyle{ L^2[0,1] }[/math] of the square-integrable functions on the unit interval. Both are systems of bounded functions, unlike, say, the Haar system or the Franklin system.

Both trigonometric and Walsh systems admit natural extension by periodicity from the unit interval to the real line. Furthermore, both Fourier analysis on the unit interval (Fourier series) and on the real line (Fourier transform) have their digital counterparts defined via Walsh system, the Walsh series analogous to the Fourier series, and the Hadamard transform analogous to the Fourier transform.

Properties

The Walsh system [math]\displaystyle{ \{W_k\}, k \in \mathbb{N}_0 }[/math] is an abelian multiplicative discrete group isomorphic to [math]\displaystyle{ \coprod_{n=0}^\infty \mathbb{Z}/2\mathbb{Z} }[/math], the Pontryagin dual of the Cantor group [math]\displaystyle{ \prod_{n=0}^\infty \mathbb{Z}/2\mathbb{Z} }[/math]. Its identity is [math]\displaystyle{ W_0 }[/math], and every element is of order two (that is, self-inverse).

The Walsh system is an orthonormal basis of the Hilbert space [math]\displaystyle{ L^2[0,1] }[/math]. Orthonormality means

[math]\displaystyle{ \int_0^1 W_k(x)W_l(x) dx = \delta_{kl} }[/math],

and being a basis means that if, for every [math]\displaystyle{ f \in L^2[0,1] }[/math], we set [math]\displaystyle{ f_k = \int_0^1 f(x)W_k(x)dx }[/math] then

[math]\displaystyle{ \int_0^1 ( f(x) - \sum_{k=0}^N f_k W_k(x) )^2 dx \;\ \xrightarrow[N\rightarrow\infty]{}\;\ 0 }[/math]

It turns out that for every [math]\displaystyle{ f \in L^2[0,1] }[/math], the series [math]\displaystyle{ \sum_{k=0}^\infty f_k W_k(x) }[/math] converges to [math]\displaystyle{ f(x) }[/math] for almost every [math]\displaystyle{ x \in [0,1] }[/math].

The Walsh system (in Walsh-Paley numeration) forms a Schauder basis in [math]\displaystyle{ L^p[0,1] }[/math],   [math]\displaystyle{ 1 \lt p \lt \infty }[/math]. Note that, unlike the Haar system, and like the trigonometric system, this basis is not unconditional, nor is the system a Schauder basis in [math]\displaystyle{ L^1[0,1] }[/math].

Generalizations

Walsh-Ferleger systems

Let [math]\displaystyle{ \mathbb{D} = \prod_{n=1}^\infty \mathbb{Z}/2\mathbb{Z} }[/math] be the compact Cantor group endowed with Haar measure and let [math]\displaystyle{ \hat {\mathbb D} = \coprod_{n=1}^\infty \mathbb Z / 2\mathbb Z }[/math] be its discrete group of characters. Elements of [math]\displaystyle{ \hat {\mathbb D} }[/math] are readily identified with Walsh functions. Of course, the characters are defined on [math]\displaystyle{ \mathbb D }[/math] while Walsh functions are defined on the unit interval, but since there exists a modulo zero isomorphism between these measure spaces, measurable functions on them are identified via isometry.

Then basic representation theory suggests the following broad generalization of the concept of Walsh system.

For an arbitrary Banach space [math]\displaystyle{ (X, ||\cdot||) }[/math] let [math]\displaystyle{ \{ R_t \}_{t \in \mathbb D} \subset \operatorname{Aut}X }[/math] be a strongly continuous, uniformly bounded faithful action of [math]\displaystyle{ \mathbb D }[/math] on X. For every [math]\displaystyle{ \gamma \in \hat{\mathbb D} }[/math], consider its eigenspace [math]\displaystyle{ X_\gamma = \{x \in X : R_t x = \gamma(t)x \} }[/math]. Then X is the closed linear span of the eigenspaces: [math]\displaystyle{ X = \overline{\operatorname{Span}}(X_\gamma, \gamma \in \hat {\mathbb D}) }[/math]. Assume that every eigenspace is one-dimensional and pick an element [math]\displaystyle{ w_\gamma \in X_\gamma }[/math] such that [math]\displaystyle{ \|w_\gamma\| = 1 }[/math]. Then the system [math]\displaystyle{ \{w_\gamma\}_{\gamma \in \hat {\mathbb D}} }[/math], or the same system in the Walsh-Paley numeration of the characters [math]\displaystyle{ \{w_k\}_{k \in {\mathbb N}_0} }[/math] is called generalized Walsh system associated with action [math]\displaystyle{ \{ R_t \}_{t \in \mathbb D} }[/math]. Classical Walsh system becomes a special case, namely, for

[math]\displaystyle{ R_t: x = \sum_{j=1}^\infty x_j2^{-j} \mapsto \sum_{j=1}^\infty (x_j \oplus t_j)2^{-j} }[/math]

where [math]\displaystyle{ \oplus }[/math] is addition modulo 2.

In the early 1990s, Serge Ferleger and Fyodor Sukochev showed that in a broad class of Banach spaces (so called UMD spaces[4]) generalized Walsh systems have many properties similar to the classical one: they form a Schauder basis[5] and a uniform finite-dimensional decomposition[6] in the space, have property of random unconditional convergence.[7] One important example of generalized Walsh system is Fermion Walsh system in non-commutative Lp spaces associated with hyperfinite type II factor.

Fermion Walsh system

The Fermion Walsh system is a non-commutative, or "quantum" analog of the classical Walsh system. Unlike the latter, it consists of operators, not functions. Nevertheless, both systems share many important properties, e.g., both form an orthonormal basis in corresponding Hilbert space, or Schauder basis in corresponding symmetric spaces. Elements of the Fermion Walsh system are called Walsh operators.

The term Fermion in the name of the system is explained by the fact that the enveloping operator space, the so-called hyperfinite type II factor [math]\displaystyle{ \mathcal R }[/math], may be viewed as the space of observables of the system of countably infinite number of distinct spin [math]\displaystyle{ 1/2 }[/math] fermions. Each Rademacher operator acts on one particular fermion coordinate only, and there it is a Pauli matrix. It may be identified with the observable measuring spin component of that fermion along one of the axes [math]\displaystyle{ \{x,y,z\} }[/math] in spin space. Thus, a Walsh operator measures the spin of a subset of fermions, each along its own axis.

Vilenkin system

Fix a sequence [math]\displaystyle{ \alpha = (\alpha_1,\alpha_2,...) }[/math] of integers with [math]\displaystyle{ \alpha_k \geq 2, k=1,2,\dots }[/math] and let [math]\displaystyle{ \mathbb{G} = \mathbb{G}_\alpha = \prod_{n=1}^\infty \mathbb{Z}/\alpha_k\mathbb{Z} }[/math] endowed with the product topology and the normalized Haar measure. Define [math]\displaystyle{ A_0 = 1 }[/math] and [math]\displaystyle{ A_k = \alpha_1 \alpha_2 \dots \alpha_{k-1} }[/math]. Each [math]\displaystyle{ x \in \mathbb G }[/math] can be associated with the real number

[math]\displaystyle{ \left|x\right| = \sum_{k=1}^{\infty} \frac{x_k}{A_{k}} \in \left[0,1\right]. }[/math]

This correspondence is a module zero isomorphism between [math]\displaystyle{ \mathbb G }[/math] and the unit interval. It also defines a norm which generates the topology of [math]\displaystyle{ \mathbb G }[/math]. For [math]\displaystyle{ k=1,2,\dots }[/math], let [math]\displaystyle{ \rho_k: \mathbb{G}\to\mathbb{C} }[/math] where

[math]\displaystyle{ \rho_k(x) = \exp\left(i\frac{2 \pi x_k}{\alpha_k}\right) = \cos\left(\frac{2 \pi x_k}{\alpha_k}\right) + i \sin\left(\frac{2 \pi x_k}{\alpha_k}\right). }[/math]

The set [math]\displaystyle{ \{\rho_k\} }[/math] is called generalized Rademacher system. The Vilenkin system is the group [math]\displaystyle{ \hat{\mathbb G} = \coprod_{n=1}^\infty \mathbb{Z}/\alpha_k\mathbb{Z} }[/math] of (complex-valued) characters of [math]\displaystyle{ \mathbb G }[/math], which are all finite products of [math]\displaystyle{ \{\rho_k\} }[/math]. For each non-negative integer [math]\displaystyle{ n }[/math] there is a unique sequence [math]\displaystyle{ n_0, n_1, \dots }[/math] such that [math]\displaystyle{ 0 \leq n_k \lt \alpha_{k+1}, k=0,1,2,\dots }[/math] and

[math]\displaystyle{ n = \sum_{k=0}^{\infty} n_k A_k. }[/math]

Then [math]\displaystyle{ \hat{\mathbb G} = {\chi_n | n=0,1,\dots} }[/math] where

[math]\displaystyle{ \chi_n = \sum_{k=0}^{\infty} \rho_{k+1}^{n_k}. }[/math]

In particular, if [math]\displaystyle{ \alpha_k = 2, k=1,2... }[/math], then [math]\displaystyle{ \mathbb G }[/math] is the Cantor group and [math]\displaystyle{ \hat{\mathbb G} = \left\{\chi_n | n=0,1,\dots\right\} }[/math] is the (real-valued) Walsh-Paley system.

The Vilenkin system is a complete orthonormal system on [math]\displaystyle{ \mathbb G }[/math] and forms a Schauder basis in [math]\displaystyle{ L^p(\mathbb{G}, \mathbb{C}) }[/math],  [math]\displaystyle{ 1 \lt p \lt \infty }[/math].[8]

Binary Surfaces

Romanuke showed that Walsh functions can be generalized to binary surfaces in a particular case of function of two variables.[9] There also exist eight Walsh-like bases of orthonormal binary functions,[10] whose structure is nonregular (unlike the structure of Walsh functions). These eight bases are generalized to surfaces (in the case of the function of two variables) also. It was proved that piecewise-constant functions can be represented within each of nine bases (including the Walsh functions basis) as finite sums of binary functions, when weighted with proper coefficients.[11]

Nonlinear Phase Extensions

Nonlinear phase extensions of discrete Walsh-Hadamard transform were developed. It was shown that the nonlinear phase basis functions with improved cross-correlation properties significantly outperform the traditional Walsh codes in code division multiple access (CDMA) communications.[12]

Applications

Applications of the Walsh functions can be found wherever digit representations are used, including speech recognition, medical and biological image processing, and digital holography.

For example, the fast Walsh–Hadamard transform (FWHT) may be used in the analysis of digital quasi-Monte Carlo methods. In radio astronomy, Walsh functions can help reduce the effects of electrical crosstalk between antenna signals. They are also used in passive LCD panels as X and Y binary driving waveforms where the autocorrelation between X and Y can be made minimal for pixels that are off.

See also

Notes

  1. ↑ Walsh 1923.
  2. ↑ Fine 1949.
  3. ↑ Schipp, Wade & Simon 1990.
  4. ↑ Pisier 2011.
  5. ↑ Sukochev & Ferleger 1995.
  6. ↑ Ferleger & Sukochev 1996.
  7. ↑ Ferleger 1998.
  8. ↑ Young 1976
  9. ↑ Romanuke 2007.
  10. ↑ Romanuke 2010b.
  11. ↑ Romanuke 2010c.
  12. ↑ A.N. Akansu and R. Poluri, "Walsh-Like Nonlinear Phase Orthogonal Codes for Direct Sequence CDMA Communications," IEEE Trans. Signal Process., vol. 55, no. 7, pp. 3800–3806, July 2007.

References

  • Ferleger, Sergei V.; Sukochev, Fyodor A. (March 1996). "On the contractibility to a point of the linear groups of reflexive non-commutative Lp-spaces". Mathematical Proceedings of the Cambridge Philosophical Society 119 (3): 545–560. doi:10.1017/s0305004100074405. Bibcode: 1996MPCPS.119..545F. 
  • Romanuke, V. V. (2007). "On the Point of Generalizing the Walsh Functions to Surfaces". Herald of Khmelnytskyi National University 6 (1): 187–193. 
  • Romanuke, V. V. (2010b). "Generalization of the Eight Known Orthonormal Bases of Binary Functions to Surfaces". [full citation needed]
  • Schipp, Ferenc; Wade, W.R.; Simon, P. (1990). Walsh series. An introduction to dyadic harmonic analysis. Akadémiai Kiadó. 
  • Sukochev, Fyodor A.; Ferleger, Sergei V. (December 1995). "Harmonic analysis in (UMD)-spaces: Applications to the theory of bases". Mathematical Notes 58 (6): 1315–1326. doi:10.1007/bf02304891. 

External links



Retrieved from "https://handwiki.org/wiki/index.php?title=Walsh_function&oldid=3364918"