Use The Coiflets and Daubechies Wavelet Transform To Reduce Data Noise For a Simple Experiment

Taher, Mahmood M; Redha, Sabah Manfi

doi:10.33899/iqjoss.2022.176225

Journals List

Use The Coiflets and Daubechies Wavelet Transform To Reduce Data Noise For a Simple Experiment

IRAQI JOURNAL OF STATISTICAL SCIENCES

Article 9, Volume 19, Issue 2, December 2022, Pages 91-103 PDF (1.31 M)

Document Type: Research Paper

DOI: 10.33899/iqjoss.2022.176225

Authors

Mahmood M Taher^* ¹; Sabah Manfi Redha²

¹Department of Informatics & Statistic, College of Computer & Mathematical Science, University of Mosul, Mosul, Iraq

²Department of Statistics, College of Administration And Economics , Baghdad University, Iraq.

Abstract

In this research, a simple experiment in the field of agriculture was studied, in terms of the effect of out-of-control noise as a result of several reasons, including the effect of environmental conditions on the observations of agricultural experiments, through the use of Discrete Wavelet transformation, specifically (The Coiflets transform of wavelength 1 to 2 and the Daubechies transform of wavelength 2 To 3) based on two levels of transform (J-4) and (J-5), and applying the hard threshold rules, soft and non-negative, and comparing the wavelet transformation methods using real data for an experiment with a size of 2⁶ observations. The application was carried out through a program in the language of MATLAB. The researcher concluded that using the wavelet transform with the Suggested threshold reduced the noise of observations through the comparison criteria.

Highlights

1- When applying The Discrete wavelet transform at the level (j-4), it led to a decline in the value of MSe and CV for all cases, but at the level (j-5), there was an increase and decrease in the values of the Criteria MSe and CV, in this experiment.

2- Obtaining the best results when applying the hard threshold rule with the Universal and suggested threshold according to standards.

3-When processing the observations noise (wheat crop 2⁶) and the level (j-4), the suggested threshold gave better results than the Universal Threshold based on the criteria MSe, CV, MSe_(w), SNR. and The best results were obtained for the first-order Coiflet wavelet transformation filter with the hard threshold rule, the second-order Daubechies wavelet transformation filter with the non-negative threshold rule, and second-order Coiflet wavelet transformation filter with the non-negative threshold rule, in this experiment.

4-When processing the observations noise (wheat crop 2⁶) and the level (j-5), the suggested threshold gave better results than the Universal Threshold based on the criteria MSe, CV, MSe_(w), SNR. and The best results were obtained for the first-order Coiflet wavelet transformation filter with the hard threshold rule, the second-order Daubechies wavelet transformation filter with the hard threshold rule, and the second-order Coiflet wavelet transformation filter with the hard threshold rule, in this experiment.

Keywords

Data Noise; Wavelet; Coiflets; Daubechies; Random Complete Blocks Design

Full Text

1. Introduction

The Discrete Wavelet transformation is one of the very important topics used in several application fields, especially in data noise processing. In general, noise represents an unexplained variance within the data (Reid & Reading,2010); that is, it cannot be eliminated. But it can be reduced in several ways according to the study or its field of application. In agricultural experiments, the experiment is divided into blocks. It replicates according to the different experimental units, which directly or indirectly affect the observation value, resulting from the mathematical model according to the applied design. The method or procedure cannot be controlled in most experiments due to out-of-control conditions during the application of the experiment.

The wavelet transform has been discussed in addressing pollution and heterogeneity (Ali & Mawlood,2010) for the complete randomized design using wavelet filters and some types of threshold rules. In addition, a threshold was Suggested to reduce the observations noise for a factorial experiment by (Taher&Sabah,2022) compared with the Universal Threshold. In this research, the application of different levels of analysis, through the use of the Daubechies transform of wavelengths from 2 to 3 and the transformation of Coiflets of wavelengths from 1 to 2 with hard, soft, and non-negative threshold rules, and comparison of the results.

2. Discrete Wavelet transform (DWT)

The Discrete wave transform is one of the most Transfers used in the wavelet due to its multiple applications in various practical fields and its theoretical uses in various sciences. The researcher will give a comprehensive idea of this transformation and focus on its use in designing experiments Through the application of a simple experiment.

The work of the discrete wavelet transform depends on the Mallat pyramidal algorithm, which is an efficient algorithm proposed by the researcher Mallat (1989) (Nason, 2008) to calculate the wavelet coefficients for a set of data containing noise

The principle of The work of this algorithm is to create filters for smoothing and heterogeneity from the wavelet coefficients, and these filters are used frequently to obtain data for all scales, meaning that the wavelet transform splits the data into two components, the first component is called detail, which includes high frequencies and can be calculated from the mother wavelet by the following formula ( In & Kim, 2013)

Whereas

Represents a high-pass filter (In & Kim, 2013)

he Second component is called Approximate, which includes low frequencies (noise) or (anomalous) values according to the nature of the study and its application, and it can be calculated from the father wavelet by the following formula: (In, F., & Kim, S 2013)

Whereas Represents a low-pass filter, (In & Kim, 2013).and it is related to through

The following figure shows the division of data into two components

Low pass filter

Approximation coefficients

X[n]

High pass filter

Detail coefficients

Data filter coefficients

Figure (1): High-pass filter and low-pass filter for x vector

In general, the discrete wavelet transform is used with data that contain discrete variables and have discrete outputs.

3. The Transfer levels Multi Resolution Analysis

We will give a brief overview of the key concepts of multiscale analysis before attempting formal definitions of wavelets and the wavelet transform and How we extract multiscale "information" from the vector y . We identify the "detail" in the sequence at various scales and places as the essential information.

Transform levels are determined from the design observation, and through the application side, we will have a simple experiment containing sixteen treatments and four blocks represented by the following vector.

Through the observations vector, the levels of analysis for this experiment will be (Nason, 2008).

The next stage is how to extract information from vector , where the information extracted from vector It is called (detail), which can be obtained from different locations and levels, and in general, the word "detail" means "degree of difference" or "variance" in the observations of the vector. This information is calculated based on the following two equations (Nason,2008).

The next step is calculating the elaboration and measurement coefficients for the other levels (Nason,2008).

This process is called the multiscale transform algorithm, Through the simple experiment applied, we will have the following levels { (j) , (j-1) , (j-2) , (j-3) , (j-4) , (j-5)}.

4. Stages of Discrete Wavelet Transform Using Orthogonality

4.1. Representing the observations of a randomized complete block design (RCBD) with a vector That contains all the observations of the experiment, which the following mathematical model represents

Whereas

observation value (j) from treatment (i), The general arithmetic mean, The effect of the i-treatment for this observation, The amount of random error , number of treatments number of blocks

The vector x contains all the observations of the experiment. One of the critical conditions in the wavelet transformation is The size of the observations fulfills the following condition.

4.2.Applying Coiflets transform and the wavelet Daubechies on observation  the Experiment, we will have wavelet coefficients that can be represented Using an orthogonal matrix 
 
, and multiplying it by the observations vector, the Wavelet Coiflets and Daubechies coefficients vector of the following form is Obtained.

whereas

The wavelet coefficients vector

Orthogonal matrix

Vector observation

From the formula (9), we get a vector coefficient of a discrete wavelet which can be represented in the following form.

Approximation coefficients

Whereas, Represents the first component of the transform, which is the detail coefficients computed from the rate of the difference of the data at each measurement and is symbolized by CD. As for Represents the second component of the transform, which is the approximation coefficients and represents the rate of The measurement is symbolized by CA.

4.3. We note that the formula (6) can be obtained through which the values of the original Data Depending on the orthogonality condition of the Discrete wavelet coefficients.

To clarify what was mentioned above, we take the following figure, which shows the discrete wavelet transformation coefficients for the data for four levels .

DWT₁

Approximation

CAA

CAD

Level 1

Level 2

DWT₂

CAAA

CAAD

CAD

Level 3

DWT₃

CAAAA

CAAAD

CAAD

CAD

Level 4

DWT₄

Details

Figure(2): The wavelet transform levels for size N = 32

Several types of Wavelets exist through the above offer About Discrete Wavelet transformation and orthogonality. In this research, we will address: Daubechies Wavelet, Coiflets Wavelet

4.4. Threshold rules (hard, soft, non-negative) are applied with threshold limit Universal Threshold and suggested threshold.

Equation (7) represents the Universal Threshold (UT) value (Gençay et al.,2001) (Zhang et al., 2021), while equation (8) represents the Suggested Threshold (ST) (Taher & Ridha, 2022).

5. Daubechies Wavelet

The name came in relation to the researcher Ingrid Daubechies (Tammireddy & Tammu, 2014). It has made a boom in the wavelet theory, as it is generated from a group of intermittent wavelets. The most crucial characteristic of this family of wavelets is their smoothness. It is abbreviated as follows.

where's

Db: Abbreviation for the name of the researcher Daubechies

N: Wavelet rank

considered wavelet haar of one the members of this family of wavelets and is symbolized by the symbol db1 or The wavelet is called Daubechies the first order Because built from the function of the father (father wave), the function of the mother (mother wave), as follows (Tammireddy & Tammu, 2014).

father wavelet Daubechies

Where Daubechies scaling function coefficients

Mother wavelet Daubechies

Where Daubechies wavelet function coefficients

Figure Next represents a Daubechies Wavelet of several lengths

Figure(3): the wave function and The scale function of a wave Daubechies

6. Coiflets Wavelet

They are discrete wavelets designed by Ingrid Daubechies, at the request of Ronald Coifman (Tammireddy & Tammu, 2014).

where's

: Abbreviation for the name of the researcher Ronald Coifman

: Wavelet rank

This family of wavelets is characterized by the presence of a relationship that relates the length of the filter with its rank

It is also built from the father function (father wavelet) and mother function (mother wavelet), in the following formula (Tammireddy & Tammu, 2014)

father wavelet Coiflets

Where Coiflets scaling function coefficients

Mother wavelet Coiflets

Where Coiflets wavelet function coefficients

Figure(4): the wave function and The scale function of a wave Coiflets

This family of wavelets is considered orthogonal and is close to symmetry, as it connects the mother and father function through the high-pass filter and the low-pass filter by obtaining the vanishing torque, unlike each filter separately.

7. Threshold Rules

Three types of threshold rules will be applied in this research, namely:

7.1. Hard Threshold Rule

One of the types of threshold rules, it is applied in discrete wavelet transform and takes the following form (Dehda & Melkemi, 2017),( Zaeni et al.,2018).

Through the formula (9) where the coefficients whose values are greater and equal to the threshold do not change, and the coefficients whose values are less than the threshold are replaced by the value zero

7.2. Soft Threshold Rule

It is another type of threshold rule and can be written as (Bruce & Gao , 1995).

And It is considered one of the standard techniques for processing observational noise (Tang et al., 2013) (Han & Xu, 2016).

Through the above formulas (10) and (11), we note that if the coefficients of the wavelet are less than the threshold value, it goes to zero, but in the case of being greater than the threshold value, it preserves its value. Using a Shrinking wavelet based on a soft-threshold rule tends to bias because all large coefficients Shrink towards zero.

7.3 Non-Negative Garrote Threshold Rule

This threshold is characterized by the small samples. It is less sensitive to observation than the hard threshold, especially in small fluctuations, and it is less biased than the soft threshold and can be written as follows (Gao,1997).

It is also the base of wavelet Shrink introduced by Gao

8. Evaluation Criteria

For the purpose of comparing the results between the Coiflets Wavelet Transformation and the Daubechies Wavelet Transformation, several criteria were applied,

8.1.The mean square error for design is defined by the following formula(Montgomery, 2020):

Where total sum of squares , treatments Sum of squares , blocks Sum of squares

8.2. The coefficient of variation(cv) is defined by the following formula: (Montgomery, 2020):

8.3. The mean of squares for original and transform observation is defined by the following
Formula (He et al.,2014):

, observation value of experiment , Experiment data after transformation

8.4. The Signal-to-noise ratio (SNR) is defined by the following formula(He et al.,2014):

Where’s: ,

and The following diagram shows the steps of the wavelet transform, with the comparison criteria

start

Input Data

Coiflets Wavelet

Daubchies Wavelet

Calculation of the Universal Threshold value

Calculation of the Suggested Threshold value

Evaluation Criteria's

MSe ,CV. MSe_(w),SNR

Best Vector

Application Hard rule with Universal

Application Hard rule with Suggested

Application Non N rule with Universal

Application Non N rule with Suggested

Application Soft rule with Universal

Application Soft rule with Suggested

* The Scheme was prepared by the researcher

Figure (5): Stages of The Wavelet Transform in this research

9. The practical side

The field experiment on wheat cultivation was conducted in one of the stations of the National Program for the Development of Wheat Cultivation in Iraq, and sixteen wheat varieties were included.(al sds (12) A1, al sahel (1)A2 , al sds (1)A3 , Egypt (1)A4 .Egypt (2) A5, Sakha (93)A6, al Geza (11)A7, al Geza (168)A8, Apaa (99)A9, Italia(1)A10, Italia (2)A11, Caronia A12, gold kernels A13, aom al rabee A14, Smitto A15, Waha A16).

According to the complete random blocks design (CRBD) of four blocks, each block contained 16 experimental units. Then the characteristics of the field yield were taken, which are (number of branches, plant height cm, dry weight of g, number of spikes / m², number of seeds/spike, weight of 1000 grains of/ g, Grain yield / m²), where the trait was studied: grain yield / m².

By applying equations (10-19) on experiment observations, The results are shown in table (1), Which represents a summary of wavelet transformation when the level of analysis is (J-4)

Table (1): The best results of wavelet transformation when the

transformation level is (j-4) for The first experiment

wavelets	Data	MSe	CV	MSe(w)	SNR	Test Normal data	Test Homogeneity of variance
	X	16769	0.30			normal	Pv>0.05
DbN	XD2h3st	16439	0.30	4599.2	10.21	P.v>0.15		P.v=0.93
	XD2h4st	13560	0.27	7310.5	8.19	P.v=0.14		P.v=0.89
	XD2s2st	10373	0.27	7595.1	8.03	P.v=0.09		P.v=0.91
	XD2s3st	8634	0.27	15902	4.82	P.v>0.15		P.v=0.91
	XD2n2st	13908	0.28	2767.4	12.41	P.v>0.15		P.v=0.92
	XD2n3st	11832	0.27	7496.4	8.08	P.v>0.15		P.v=0.91
	XD3h3st	16313	0.30	5589.8	9.36	P.v=0.09		P.v=0.82
	XD3s2st	10953	0.28	11498.9	6.23	P.v=0.02		P.v=0.17
	XD3n2st	14023	0.28	4695.9	10.12	P.v=0.13		P.v=0.96
	XD3n3st	12544	0.28	10820.9	6.49	P.v=0.02		P.v=0.97
Coif N	XC1h2st	16207	0.29	2099.5	13.61	P.v>0.15		P.v=0.94
	XC1h3st	16762	0.30	5108.8	9.75	P.v>0.15		P.v=0.99
	XC1s2st	9374	0.26	10088.8	6.79	P.v>0.15		P.v=0.96
	XC1n2st	12623	0.27	3812.5	11.02	P.v>0.15		P.v=0.97
	XC1n3st	9621	0.24	10004.1	6.83	P.v=0.08		P.v=0.96
	XC2h3st	15856	0.29	3893.4	10.93	P.v>0.15		P.v=0.98
	XC2s2st	9109	0.26	9520.9	7.05	P.v=0.1		P.v=0.95
	XC2n2st	12508	0.26	3457.2	11.45	P.v>0.15		P.v=0.95
	XC2n3st	9154	0.23	9257.2	7.17	P.v=0.02		P.v=0.98

By applying equations (10-19) on experiment observations, The results are shown in table (2), Which represents a summary of wavelet transformation when the level of analysis is (J-5)

Table (2): The best results of wavelet transformation when the

transformation level is (j-5) for The first experiment

wavelets	Data	MSe	CV	MSe(w)	SNR	Test Normal data	Test Homogeneity of variance
	X	16769	0.30			normal	Pv>0.05
DbN	XD2h2th	17011	0.30	1587.1	14.83	P.v>0.15	P.v=0.89
	XD2h3th	17315	0.31	5925.5	9.11	P.v>0.15	P.v=0.98
	XD2s2th	10950	0.31	13444.8	5.55	P.v>0.15	P.v=0.96
	XD2n2th	14574	0.30	4127.1	10.68	P.v>0.15	P.v=0.95
	XD2n3th	13487	0.32	10955.3	6.44	P.v>0.15	P.v=0.95
	XD3h2th	18447	0.32	1448.4	15.22	P.v>0.15	P.v=0.92
	XD3h3th	14503	0.28	4554.9	10.25	P.v>0.15	P.v=0.98
	XD3s2th	12318	0.33	12287.2	5.94	P.v>0.15	P.v=0.95
	XD3n2th	15620	0.31	3806.6	11.03	P.v>0.15	P.v=0.96
	XD3n3th	15619	0.34	9804.4	6.92	P.v>0.15	P.v=0.97
Coif N	XC1h2th	16246	0.29	1522.9	15.01	P.v>0.15	P.v=0.85
	XC1s2th	13728	0.35	13167.1	5.64	P.v>0.15	P.v=0.88
	XC1n2th	16864	0.32	4029.3	10.78	P.v>0.15	P.v=0.92
	XC1n3th	16400	0.35	11277.3	6.31	P.v>0.15	P.v=0.91
	XC2h2th	15950	0.29	1486.1	15.11	P.v>0.15	P.v=0.86
	XC2h3th	18164	0.32	6210.8	8.90	P.v>0.15	P.v=0.95
	XC2s2th	12904	0.34	13800.4	5.43	P.v>0.15	P.v=0.86
	XC2n2th	16326	0.32	4230	10.57	P.v>0.15	P.v=0.89
	XC2n3th	15105	0.34	11486.1	6.23	P.v>0.15	P.v=0.91

9.1 Results

Table (1) represents a summary of the wavelet transform results for all cases when the transform level is (J-4), where led to decline in the MSe value of the used design and a significant improvement in the CV criterion. In addition to obtaining low values for the standard .with an increase in the SNR criterion value, Especially for ( XD2h3st , XD2n2st , XD3n2st , XC1h2st , XC2h3st , XC2n2st )

Table(2) represents a summary of the wavelet transform results for all cases when the transform level is (J-5), where led to decline in the MSe value of the used design and a significant improvement in the CV criterion. In addition to obtaining low values for the standard .with an increase in the SNR criterion value, Especially for( XD2n2th, XD3h3th , XC1h2th, XC2h2th).

Where XD2h3st : Represents the second-order Daubechies wavelet transformation filter with the hard rule and suggested threshold, And so on for the rest of the vectors

10. Conclusion

2- Obtaining the best results when applying the hard threshold rule with the Universal and suggested threshold according to standards.

References

[1] Ali, Taha, Husean & Mawlod , Kurdestan ,Ibrahem (2010)." Addressing the problem of contamination and variance heterogeneity in a complete random design using a small wave filter",Iraqi Journal of Statistical Science,No.18,Issue.10.

[2]. Bruce, A. G., & Gao, H. Y. (1995, September). WaveShrink: Shrinkage functions and thresholds. In wavelet applications in signal and image processing III (Vol. 2569, pp. 270-281). SPIE.‏

[3]. Dehda, B., & Melkemi, K. (2017). Image denoising using new wavelet thresholding , Function. Journal of Applied Mathematics and Computational Mechanics, 16(2), 55-65.‏

[4]. Gao, H. Y. (1997). Wavelet Shrinkage Denoising Using the Non-Negative Garrote," Mathsoft. Inc. Tech. Rep.‏

[5]. Gençay, R., Selçuk, F., & Whitcher, B. J. (2001). An introduction to wavelets and other Filtering methods in finance and economics. Elsevier.‏

[6]. Han, G., & Xu, Z. (2016). Electrocardiogram signal denoising based on a new improved, Wavelet thresholding. Review of Scientific Instruments, 87(8), 084303.‏

[7]. He, C., Xing, J. C., & Yang, Q. L. (2014). Optimal wavelet basis selection for wavelet denoising of structural vibration signal. In Applied Mechanics and Materials (Vol. 578, pp. 1059-1063). Trans Tech Publications Ltd.‏

[8]. In, F., & Kim, S. (2013). An introduction to wavelet theory in finance: a wavelet multiscale Approach. World Scientific.‏

[9]. Montgomery, D. C. (2020). Design and analysis of experiments. John Wiley & sons.‏

[10]. Nason, G. P. (Ed.). (2008). Wavelet methods in statistics with R. New York, NY: Springer New York.‏

[11]. Taher, M. M., & Ridha, S. M. (2022). The suggested threshold to reduce data noise for A factorial experiment. International Journal of Nonlinear Analysis and Applications, 13(1), 3861-3872.

[12]. Tammireddy, P. R., & Tammu, R. (2014). Image reconstruction using wavelet Transform with extended fractional Fourier transform.‏

[13]. Tang, H., Liu, Z. L., Chen, L., & Chen, Z. Y. (2013). Wavelet image denoising based the new threshold function. In Applied Mechanics and Materials (Vol. 347, pp. 2231-2235). Trans Tech Publications Ltd.‏

[14]. Zaeni, A., Kasnalestari, T., & Khayam, U. (2018, October). Application of wavelet Transformation symlet type and coiflet type for partial discharge signals denoising. In 2018 5th International Conference on Electric Vehicular Technology (ICEVT) (pp. 78- 82). IEEE.‏

[15]. Zhang, Y., Zhou, H., Dong, Y., & Wang, L. (2021). Restraining EMI of Displacement Sensors Based on Wavelet Fuzzy Threshold Denoising. In Signal and Information Processing, Networking, and Computers (pp. 543-551). Springer, Singapore.‏

1. Introduction

The wavelet transform has been discussed in addressing pollution and heterogeneity (Ali & Mawlood,2010)for the complete randomized design using wavelet filters and some types of threshold rules. In addition, a threshold was Suggested to reduce the observations noise for a factorial experiment by (Taher&Sabah,2022) compared with the Universal Threshold. In this research, the application of different levels of analysis, through the use of the Daubechies transform of wavelengths from 2 to 3 and the transformation of Coiflets of wavelengths from 1 to 2 with hard, soft, and non-negative threshold rules, and comparison of the results.

2. Discrete Wavelet transform (DWT)

Whereas

Represents a high-pass filter (In & Kim, 2013)

Whereas Represents a low-pass filter, (In & Kim, 2013).and it is related to through

The following figure shows the division of data into two components

Low passfilter

Approximation coefficients

X[n]

High pass filter

Detail coefficients

Data filter coefficients

Figure (1): High-pass filter and low-pass filter for x vector

In general, the discrete wavelet transform is used with data that contain discrete variables and have discrete outputs.

3. The Transfer levels Multi Resolution Analysis

Through the observations vector, the levels of analysis for this experiment will be (Nason, 2008).

The next step is calculating the elaboration and measurement coefficients for the other levels (Nason,2008).

This process is called the multiscale transform algorithm, Through the simple experiment applied, we will have the following levels { (j) , (j-1) , (j-2) , (j-3) , (j-4) , (j-5)}.

4. Stages of Discrete Wavelet Transform Using Orthogonality

Whereas

observation value (j) from treatment (i), The general arithmetic mean, The effect of the i-treatment for this observation, The amount of random error , number of treatments number of blocks

The vector x contains all the observations of the experiment. One of the critical conditions in the wavelet transformation is The size of the observations fulfills the following condition.

4.2.Applying Coiflets transform and the wavelet Daubechies on observation  the Experiment, we will have wavelet coefficients that can be represented Using an orthogonal matrix , and multiplying it by the observations vector, the Wavelet Coiflets and Daubechies coefficients vector of the following form is Obtained.

whereas

The wavelet coefficients vector

Orthogonal matrix

Vector observation

From the formula (9), we get a vector coefficient of a discrete wavelet which can be represented in the following form.

Approximation coefficients

4.3. We note that the formula (6) can be obtained through which the values of the original Data Depending on the orthogonality condition of the Discrete wavelet coefficients.

To clarify what was mentioned above, we take the following figure, which shows the discrete wavelet transformation coefficients for the data for four levels .

DWT₁

Approximation

CAA

CAD

Level 1

Level 2

DWT₂

CAAA

CAAD

CAD

Level 3

DWT₃

CAAAA

CAAAD

CAAD

CAD

Level 4

DWT₄

Details

Figure(2):The wavelet transform levels for size N = 32

Several types of Wavelets exist through the above offer About Discrete Wavelet transformation and orthogonality. In this research, we will address: Daubechies Wavelet, Coiflets Wavelet

4.4.Threshold rules (hard, soft, non-negative) are applied with threshold limit Universal Threshold and suggested threshold.

Equation (7) represents the Universal Threshold (UT) value (Gençay et al.,2001)(Zhang et al., 2021), while equation (8) represents the Suggested Threshold (ST) (Taher & Ridha, 2022).

5. Daubechies Wavelet

where's

Db: Abbreviation for the name of the researcher Daubechies

N: Wavelet rank

father wavelet Daubechies

Where Daubechies scaling function coefficients

Mother wavelet Daubechies

Where Daubechies wavelet function coefficients

Figure Next represents a Daubechies Wavelet of several lengths

Figure(3): the wave function and The scale function of a wave Daubechies

6. Coiflets Wavelet

They are discrete wavelets designed by Ingrid Daubechies, at the request of Ronald Coifman (Tammireddy & Tammu, 2014).

where's

: Abbreviation for the name of the researcher Ronald Coifman

: Wavelet rank

This family of wavelets is characterized by the presence of a relationship that relates the length of the filter with its rank

It is also built from the father function (father wavelet) and mother function (mother wavelet), in the following formula (Tammireddy & Tammu, 2014)

father wavelet Coiflets

Where Coiflets scaling function coefficients

Mother wavelet Coiflets

Where Coiflets wavelet function coefficients

Figure(4): the wave function and The scale function of a wave Coiflets

7. Threshold Rules

Three types of threshold rules will be applied in this research, namely:

7.1. Hard Threshold Rule

One of the types of threshold rules, it is applied in discrete wavelet transform and takes the following form (Dehda & Melkemi, 2017),( Zaeni et al.,2018).

7.2. Soft Threshold Rule

It is another type of threshold rule and can be written as (Bruce & Gao , 1995).

And It is considered one of the standard techniques for processing observational noise (Tang et al., 2013) (Han & Xu, 2016).

7.3 Non-Negative Garrote Threshold Rule

It is also the base of wavelet Shrink introduced by Gao

8. Evaluation Criteria

For the purpose of comparing the results between the Coiflets Wavelet Transformation and the Daubechies Wavelet Transformation, several criteria were applied,

8.1.The mean square error for design is defined by the following formula(Montgomery, 2020):

Where total sum of squares , treatments Sum of squares , blocks Sum of squares

8.2. The coefficient of variation(cv)is defined by the following formula: (Montgomery, 2020):

8.3. The mean of squares for original and transform observation is defined by the following
Formula (He et al.,2014):

, observation value of experiment , Experiment data after transformation

8.4.The Signal-to-noise ratio (SNR) is defined by the following formula(He et al.,2014):

Where’s: ,

and The following diagram shows the steps of the wavelet transform, with the comparison criteria

start

Input Data

Coiflets Wavelet

Daubchies Wavelet

Calculation of the Universal Threshold value

Calculation of the Suggested Threshold value

Evaluation Criteria's

MSe ,CV. MSe_(w),SNR

Best Vector

Application Hard rule with Universal

Application Hard rule with Suggested

Application Non N rule with Universal

Application Non N rule with Suggested

Application Soft rulewith Universal

Application Soft rulewith Suggested

*The Scheme was prepared by the researcher

Figure (5): Stages of The Wavelet Transform in this research

9. The practical side

By applying equations (10-19) on experiment observations, The results are shown in table (1), Which represents a summary of wavelet transformation when the level of analysis is (J-4)

Table (1): The best results of wavelet transformation when the

transformation level is (j-4) for The first experiment

wavelets	Data	MSe	CV	MSe(w)	SNR	Test Normal data	Test Homogeneity of variance
	X	16769	0.30			normal	Pv>0.05
DbN	XD2h3st	16439	0.30	4599.2	10.21	P.v>0.15		P.v=0.93
	XD2h4st	13560	0.27	7310.5	8.19	P.v=0.14		P.v=0.89
	XD2s2st	10373	0.27	7595.1	8.03	P.v=0.09		P.v=0.91
	XD2s3st	8634	0.27	15902	4.82	P.v>0.15		P.v=0.91
	XD2n2st	13908	0.28	2767.4	12.41	P.v>0.15		P.v=0.92
	XD2n3st	11832	0.27	7496.4	8.08	P.v>0.15		P.v=0.91
	XD3h3st	16313	0.30	5589.8	9.36	P.v=0.09		P.v=0.82
	XD3s2st	10953	0.28	11498.9	6.23	P.v=0.02		P.v=0.17
	XD3n2st	14023	0.28	4695.9	10.12	P.v=0.13		P.v=0.96
	XD3n3st	12544	0.28	10820.9	6.49	P.v=0.02		P.v=0.97
Coif N	XC1h2st	16207	0.29	2099.5	13.61	P.v>0.15		P.v=0.94
	XC1h3st	16762	0.30	5108.8	9.75	P.v>0.15		P.v=0.99
	XC1s2st	9374	0.26	10088.8	6.79	P.v>0.15		P.v=0.96
	XC1n2st	12623	0.27	3812.5	11.02	P.v>0.15		P.v=0.97
	XC1n3st	9621	0.24	10004.1	6.83	P.v=0.08		P.v=0.96
	XC2h3st	15856	0.29	3893.4	10.93	P.v>0.15		P.v=0.98
	XC2s2st	9109	0.26	9520.9	7.05	P.v=0.1		P.v=0.95
	XC2n2st	12508	0.26	3457.2	11.45	P.v>0.15		P.v=0.95
	XC2n3st	9154	0.23	9257.2	7.17	P.v=0.02		P.v=0.98

By applying equations (10-19) on experiment observations, The results are shown in table (2), Which represents a summary of wavelet transformation when the level of analysis is (J-5)

Table (2): The best results of wavelet transformation when the

transformation level is (j-5) for The first experiment

wavelets	Data	MSe	CV	MSe(w)	SNR	Test Normal data	Test Homogeneity of variance
	X	16769	0.30			normal	Pv>0.05
DbN	XD2h2th	17011	0.30	1587.1	14.83	P.v>0.15	P.v=0.89
	XD2h3th	17315	0.31	5925.5	9.11	P.v>0.15	P.v=0.98
	XD2s2th	10950	0.31	13444.8	5.55	P.v>0.15	P.v=0.96
	XD2n2th	14574	0.30	4127.1	10.68	P.v>0.15	P.v=0.95
	XD2n3th	13487	0.32	10955.3	6.44	P.v>0.15	P.v=0.95
	XD3h2th	18447	0.32	1448.4	15.22	P.v>0.15	P.v=0.92
	XD3h3th	14503	0.28	4554.9	10.25	P.v>0.15	P.v=0.98
	XD3s2th	12318	0.33	12287.2	5.94	P.v>0.15	P.v=0.95
	XD3n2th	15620	0.31	3806.6	11.03	P.v>0.15	P.v=0.96
	XD3n3th	15619	0.34	9804.4	6.92	P.v>0.15	P.v=0.97
Coif N	XC1h2th	16246	0.29	1522.9	15.01	P.v>0.15	P.v=0.85
	XC1s2th	13728	0.35	13167.1	5.64	P.v>0.15	P.v=0.88
	XC1n2th	16864	0.32	4029.3	10.78	P.v>0.15	P.v=0.92
	XC1n3th	16400	0.35	11277.3	6.31	P.v>0.15	P.v=0.91
	XC2h2th	15950	0.29	1486.1	15.11	P.v>0.15	P.v=0.86
	XC2h3th	18164	0.32	6210.8	8.90	P.v>0.15	P.v=0.95
	XC2s2th	12904	0.34	13800.4	5.43	P.v>0.15	P.v=0.86
	XC2n2th	16326	0.32	4230	10.57	P.v>0.15	P.v=0.89
	XC2n3th	15105	0.34	11486.1	6.23	P.v>0.15	P.v=0.91

9.1 Results

Where XD2h3st :Represents the second-order Daubechies wavelet transformation filter with the hard rule and suggested threshold,And so on for the rest of the vectors

10. Conclusion

1-When applying The Discrete wavelet transform at the level (j-4), it led to a decline in the value of MSe and CV for all cases, but at the level (j-5), there was an increase and decrease in the values of the Criteria MSe and CV,in this experiment.

2-Obtaining the best results when applying the hard threshold rule with the Universal and suggested threshold according to standards.

Statistics

Article View: 341

PDF Download: 367