Abstract

Jacket-Haar transform has been recently generalized from Haar transform and Jacket transform, but, unfortunately, it is not available in a case where the length is not a power of 2. In this paper, we have proposed an arbitrary-length Jacket-Haar transform which can be conveniently constructed from the 2-point generalized Haar transforms with the fast algorithm, and thus it can be constructed with any sizes. Moreover, it can be further extended with elegant structures, which result in the fast algorithms for decomposing. We show that this approach can be practically applied for the electrocardiogram (ECG) signal processing. Simulation results show that it is more efficient than the conventional fast Fourier transform (FFT) in signal processing.

1. Introduction

Hadamard transform, Haar transform, discrete Fourier transform, and their derivatives are discrete orthogonal transforms with extensive applications in signal and image processing [14]. Jacket transform, motivated by the center weighted Hadamard transform [5], is a special transform with its inverse transform matrix being determined by the element or block-wise inverse of an original matrix [6]. Jacket transform has been extensively applied in many fields, such as signal and data processing [1], digital wireless communications [7], cryptography [8], and encoding designs [9]. Meanwhile, several interesting matrices, such as Hadamard matrices and DFT matrices, belong to the Jacket transform matrix family. Besides all, lots of widely used matrices, such as unitary matrices and Hermitian matrices, have tightly relations with Jacket matrices.

Recently, literatures relevant to discrete orthogonal matrices and their transforms mainly have involved explorations of other categories of orthogonal transform matrices and exploitations of their applications. In the later aspect, block Jacket transforms have been tentatively applied to quantum signal processing [10], Big-Data processing [11], emerging new-generation mobile communication, and so on. In addition, new orthogonal transforms, such as complex Hadamard transform [1214], fractional Hadamard or Jacket transforms [15], parametric transforms [1618], hybrid transforms [7], and the generalized orthogonal transforms [14], have been gradually proposed while enriching the orthogonal transform family. Particularly, with advancement of digital systems and widespread availability in the recent few decades, there exist some urgent demands for seeking a scheme to achieve compromise between the generalization efficiencies of Hadamard or Haar transforms and their the extended transforms to adaptively meet the practical implementation requirements. Subsequent research can be seen as an attempt.

Haar matrix, which is useful for localized signal analysis [19], edge detection [20], OFDM, and filter design and electrocardiogram (ECG) analysis, has been generalized for Jacket-Haar matrix [21], whose entries are and compared with entries of the original Haar matrix being , , and . Although the -point Jacket-Haar matrices are successfully proposed in [21], there is still a problem on how to construct the arbitrary-length Jacket-Haar transform as the arbitrary-length Walsh-Jacket transform has already done with high efficiency [22]. Unfortunately, until this paper, the method to solve this problem is preliminary and not comprehensive except that some original results are shown [23]. Consequently, in this paper, we focus on the systematic construction of the generalized Jacket-Haar transforms of any lengths. Compared to the existing literatures, the contributions are listed as follows. Firstly, the multilevel structure of the arbitrary-length Jacket-Haar transform matrix can be derived with the recursive algorithm. Secondly, a novel approach for designing the fast generation of the Jacket-Haar transform matrices of the arbitrary length is derived in a successive fashion. Finally, the arbitrary length Jacket-Haar transforms can be generated and decomposed for extensions with the fast algorithms with several specified applications.

The remaining of this paper is organized as follows. In Section 2, the elegant structure of Jacket-Haar transform is investigated, followed by the generalized form of Jacket-Haar transform of size . In Section 3, the -point Jacket-Haar transform is extensively derived from the -point Jacket-Haar transform. Together with the proposed algorithm for generating the -point Jacket-Haar transform, an arbitrary Jacket-Haar transform can be derived in a similar fashion. After that, the fast algorithm and its applications in signal processing will be shown in Section 4. In Section 5, the fast Jacket-Haar transforms are designed on the basis of the Kronecker product of the traditional Jacket and Jacket-Haar transforms. Finally, conclusion is drawn in Section 6.

2. Jacket-Haar Transform

Walsh transform can be extended for the arbitrary-point Walsh-Jacket transform of an arbitrary length [22]. Similarly, it is theoretically feasible to extend the Haar transform for the design of the Haar transform of any size. In order to verify the correctness of this judgement, we show an elegant structure of the Jacket-Haar transform.

Definition 1. Jacket-Haar transform is a discrete transform and it must satisfy three constraints as follows. (1)If and denote the forward and inverse -point Jacket-Haar transform matrices, respectively, then all entries of both matrices are equal to or , where is an arbitrary integer.(2)The number of zero crossing point of any row of is the same as the conventional Haar matrix. Namely, zero crossing point number of any row of is one except the first row, where zero crossing point does not exist.(3)Fast algorithm of the conventional Haar transform is preserved.

According to the definition of Jacket-Haar transform, all the conventional Haar transform matrices can be seen as a special case of Jacket-Haar transform matrices since the above three constraints are all satisfied. Elements of the Haar transform matrices can only be chosen from , , or , and order of any Haar matrix must be a power of 2. However, Jacket-Haar transform relaxes these limits with elements equal to or . For this reason, Jacket-Haar transform can still be efficiently implemented by bit-shifting without multiplication just like the conventional Haar transform. In addition, the Jacket-Haar transform is an extension of the Jacket transform since the constraint that where is a constant is not required any longer.

Before generating a Jacket-Haar transform of any size, the Jacket-Haar matrices of size are extensively derived, based on which the Jacket-Haar matrix of any size can be efficiently achieved.

Theorem 2. Any -point Jacket-Haar transform matrices should conform to the following forms:where , , , and are equal to or a power of with .

Proof. It is straightforward that the third constraint is always satisfied for the -point Jacket-Haar transform. Besides, the -point Jacket-Haar matrix can be generally denoted as follows:and its inverse isAccording to the first constraint, there existsMeanwhile owing to the second constraint, we can obtainSo there exist three cases as follows. Firstly, if , then it is easy to get that holds. So the -point Jacket-Haar matrix possesses the form of (1). Secondly, if , then similar to case , the -point Jacket-Haar matrix has the form of (2). Thirdly, if , according to (6), we can obtain . In the first case, if and with , then . Since (6) holds, must be equal to . Namely, . In the second case, we can get the same result. So, the -point Jacket-Haar transform conforms to (3). This completes the proof of this theorem.

3. Universal Jacket-Haar Transform Matrices

General forms of the 2-point Jacket-Haar transform matrices have been derived in the previous section. In this section, two algorithms are proposed to construct Jacket-Haar transform matrices of any large sizes.

Theorem 3. If the size of the matrix is an odd number, then one can use generation algorithm presented below to derive the -point Jacket-Haar transform matrix from both an -point Jacket-Haar transform matrix and an -point identity matrix. Let be an -point Jacket-Haar transform matrix with inverse . The , column of and the row of are described by and , respectively. Namely, one hasSuppose are the 2-point Jacket-Haar transform matrices with the respective inverse matrices . Let and denote the first and second row of , respectively. Meanwhile, and are the first and second column of , respectively. Then one haswhere and , . Let denote the column of the -point identity matrix; that is,where . Then the -point Jacket-Haar transform matrix can be derived aswith the inverse -point Jacket-Haar transform matrix

Proof. The product of and can be calculated asAccording to (20), we have , , , and that can be, respectively, calculated asHence, we have . Furthermore, it is easy to find that the first constraint and the second constraint are both satisfied. In fact, the third constraint can be easily derived. This completes the proof of this theorem.

According to Theorem 3, the -point Jacket-Haar transform matrix with matrix size being a power function of 2 or not can be conveniently generated.

Based on the above-mentioned construction, any odd integer order Jacket-Haar transform matrix can be elegantly derived, while at the same time, another construction method will be similarly proposed to construct Jacket-Haar transform matrix with the size being even numbers.

Theorem 4. If the size of the matrix is even, then one can use the generating algorithm presented below to derive the -point Jacket-Haar transform matrix from both the -point Jacket-Haar matrix and the -point identity matrix. The -point Jacket-Haar matrix can be calculated asand its inverse can be computed aswhere is the -point Jacket-Haar matrix and is its inverse matrix. The , , column of and the row of are denoted by and , respectively. Then one has and such that . Similarly, let and denote the first and second row of , respectively. Let and denote the first and second column of , respectively. Then, one haswhere and , . Moreover, the notation denotes the column of the -point identity matrix with the following constraints:where .

Proof. Since the proof is similar to that of Theorem 3, it is omitted here.

For the further illustration and clarity, we give an example in what follows. Suppose that a 5-point Jacket-Haar matrix isMeanwhile, let , , , , and , respectively. Then, we obtain a 10-point Jacket-Haar matrixwith the inverse matrix given by

In a word, based on Theorems 3 and 4, we can derive the Jacket-Haar transform matrix of any sizes in the recursive and alternate way. When is odd, we can choose Theorem 3 to derive the -point Jacket-Haar transform matrix from the -point Jacket-Haar transform. While is even, we can choose Theorem 4 to get the -point Jacket-Haar transform matrix from the -point Jacket-Haar transform. For example, when , the value of is 4 and then Theorem 3 is chosen. It means the 9-point Jacket-Haar matrix can be derived from the 5-point Jacket-Haar matrix. Next, the 5-point Jacket-Haar matrix can be gotten from the 3-point Jacket-Haar matrix with Theorem 3. In the end, the 3-point Jacket-Haar matrix can be constructed from the 2-point Jacket-Haar matrix by using Theorem 3. Therefore, the 9-point Jacket-Haar matrix can be ultimately obtained from the 2-point Jacket-Haar matrices. In Table 1 we generate the -point Jacket-Haar matrices with size . In fact, Jacket-Haar transform matrix of any sizes ultimately originates from the 2-point Jacket-Haar transform matrices.

Based on the -point Jacket-Haar transform matrices, the Jacket-Haar transform matrices of any sizes can be conveniently and elegantly constructed alternately and recursively by applying Theorems 3 and 4. Consequently, certain available sparse Jacket-Haar transform matrices can be theoretically derived.

Corollary 5. There exists an arbitrary-length sparse Jacket-Haar transform matrix with a maximum of two nonzero elements in each of the rows. In addition, there exists an arbitrary-length Jacket-Haar transform matrix, and the transpose of its inverse matrix is also a Jacket-Haar transform matrix with the following constraints:

Proof. The inductive method will be used for proving this corollary. At first, it is obvious that all the 2-point Jacket-Haar transform matrices meet the prescribed request. At the next step, there exist two cases based on which theorem has adopted generating a relatively higher order adjacent Jacket-Haar transform matrix. For the first case, suppose that is a sparse -point Jacket-Haar transform matrix. For the two-point Jacket-Haar transform matrices in the form described in (1) or (2), the -point Jacket-Haar transform matrix derived from (11) is a sparse matrix since each of the rows possesses a maximum of two nonzero elements. For the other case, if there exists a sparse -point Jacket-Haar transform matrix with the given matrices , the -point Jacket-Haar transform matrix from (18) is also a sparse matrix. This completes the proof of one part of corollary. This kind of sparse Jacket-Haar transform matrices is interesting and available in practical applications that are provided with the priority of acceptable performances promised. For another part of the proof of this corollary, the inductive method will be similarly used. At the first step, it is easy to check that the transpose of two-point inverse matrix in the form of (5) is also a Jacket-Haar transform matrix, and that (25) is satisfied. At the next step, based on the selected theorems for the construction of a relatively adjacent Jacket-Haar matrix, there exist two cases. For the first case with an -point Jacket-Haar transform matrix , the transpose of whose inverse matrix is also a Jacket-Haar transform matrix that satisfies the constraint in (25). For the 2-point Jacket-Haar transform matrices , it is obvious that derived from (12) is a -point Jacket-Haar transform matrix and the transpose of whose inverse matrix is also a Jacket-Haar transform matrix. Similarly, if is an -point Jacket-Haar transform matrix, the transpose of whose inverse matrix is also a Jacket-Haar transform matrix. For the 2-point Jacket-Haar transforms in the form described in (3), it is straightforward that derived from (18) is a -point Jacket-Haar transform matrix and the transpose of whose inverse matrix is also a Jacket-Haar transform matrix. Therefore, the whole corollary holds and the proof is complete.

It is interesting to note the Jacket-Haar transform referred to in the afore-derived corollary is similar to that of the Jacket transform. As we all know, if is an -point Jacket transform matrix, then is also an -point Jacket transform matrix with the following constraints:

4. Fast Jacket-Haar Transform

4.1. Fast Implementation Algorithms

Just like the conventional Haar transform, the Jacket-Haar transform of large sizes can be performed with fast algorithms. In order to give a more clear demonstration with simple comparison, fast algorithms of both conventional Haar transform and Jacket-Haar transform are presented simultaneously.

In Figure 1, we show the implementation structures of the conventional Haar transforms with sizes 2, 4, and , respectively. In Figure 2, the implementation structures of the 2-point, the 3-point and the 4-point Jacket-Haar transform are presented, from which the 3-point and the 4-point Jacket-Haar transforms can be decomposed into 2 and 3 butterflies, respectively. For , we can decompose the -point Jacket-Haar transform into combination of the 2-point Jacket-Haar transforms by using the methods described as in Figure 3. From Figures 1 to 3, it is obvious that the structures of the fast algorithms of the original Haar transform and the proposed Jacket-Haar transform are similar. For example, the whole implementation structure of a 9-point Jacket-Haar transform matrix shown as has been illustrated in Figure 4. It is shown that there exist layers, with , , , and butterflies being included in the adjacent layers from the left to the right. It needs arithmetic addition operations without any multiplication.

Compared to the conventional Haar transform, if both categories transform matrices possess the same matrix size which is a power of , they have the same implementation layers, numbers of butterflies, and arithmetic addition operations without any multiplication operations needed, but with different numbers of bit-shift operations. When the matrix size does not equal a power of , there do not exist the corresponding Haar transform matrices. For the proposed Jacket-Haar transform matrices, there are implementation layers and arithmetic addition operations with butterflies included, but without any multiplication operations. The details are shown in Table 2.

4.2. Generalized Jacket-Haar Transform

Definition 6. For a given square matrix , if its inverse matrix can be simply obtained from its element-wise inverse, that is, , for , where is a nonzero constant, one calls matrix a Jacket matrix; that is,and its inverse matrix iswhere denotes the operation of matrix transpose.

Definition 7. Suppose is a Jacket matrix of order and is a Jacket-Haar matrix of order . If , one calls a generalized Jacket-Haar matrix.
According to the above-mentioned definition, it is obvious that both Jacket matrix and Jacket-Haar matrix belong to generalized Jacket-Haar matrix family.

Theorem 8. Since both Jacket transform and Jacket-Haar transform can be designed with the fast algorithms, the generalized Jacket-Haar transform which is the Kronecker product of Jacket transform and Jacket-Haar transform can be composed or decomposed with the fast algorithms.

Proof. Suppose can be decomposed until and ; thenwhere is the identity matrix. The detailed procedure can be seen in [6]. Then, by using the property of Kronecker product, can be decomposed aswithDue to the fact that can be decomposed into 2-point Jacket-Haar transforms, can be decomposed or constructed with the fast algorithms. This completes the proof of this theorem.

For example, if , we can construct the 3-point Jacket-Haar transform and 4-point Jacket transform as follows:Then we can construct a 12-point generalized Jacket-Haar transformThe fast algorithm is shown in Figure 5.

5. Applications in Signal Processing

The proposed arbitrary-point Jacket-Haar transform can also be theoretically applied in areas where the conventional Haar transforms have been gained, such as signal analysis, image processing, OFDM, and filter design. In this section, two representative Jacket-Haar transforms are applied to electrocardiogram (ECG) analysis compared to the discrete Fourier transform (DFT).

There are approaches for comprehensive analysis with even and odd point transforms, respectively. The length of the ECG signal in Figure 6 is 202, which is even and not a power of 2. It is not convenient to analyze this case by directly using the original Haar transform, but it can be efficiently analyzed by the proposed arbitrary-point Jacket-Haar transform. In Figure 7, we show the normalized mean square error (NMSE) of the reconstructed signal when using part of the coefficients of the 202-point generalized Haar transformwhere is the original signal; and are the forward and inverse Jacket-Haar transform matrices. The vector preserves coefficients of and others are set to zero. In Figure 7, we also show the NMSE of the reconstructed signal when using the 202-point fast Fourier transform (FFT). The results in Figure 7 show that, with the Jacket-Haar transform, we can achieve less approximation error when using only terms () to expand the ECG signal.

As shown in Figures 8 and 9, another simulation experiment is arranged. The length of the ECG in Figure 8 is 321, which is odd and not a power of 2. We show the NMSE of the reconstructed signal in Figure 9 when using part of the coefficients of the 321-point Jacket-Haar transform compared with the 321-point FFT. It is obvious that the Jacket-Haar transform is more efficient for analyzing the ECG signal than that of FFT. In the future, the proposed Jacket-Haar transform will be applied to the more practical applications with the extensive comparisons.

6. Conclusion

We have investigated the fast construction of the Jacket-Haar transform of the arbitrary length, which overcomes of problem of the traditional constraint that the number of the points is the power of 2. With the proposed fast generation algorithms, the arbitrary length Jacket-Haar transform can be derived in a successive fashion. Moreover, we show the possible implementations of the fast algorithms and applications in signal processing. On the basis of the structures of the traditional Jacket transform and the Jacket-Haar transform with any size, the fast algorithms of Jacket-Haar transform are derived for the arbitrary length. Compared to the traditional FFT and its extensions, the proposed Jacket-Haar transform is more efficient in signal reconstruction. However, more properties of the Jacket-Haar transform and its practical applications may be further investigated in the future work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (61379153, 61572529) and partly by MEST 2015R1A2A1A05000977, NRF, Korea.