Journal of Computer Networks and Communications

Volume 2011, Article ID 235259, 10 pages

http://dx.doi.org/10.1155/2011/235259

## A Novel Technique for Transmission of M-Ary Signal through Wireless Fading Channel Using Wavelet Denoising

^{1}EEE Department, ADUST, Dhaka, Bangladesh^{2}WebSatMedia Pte. Ltd., Technopark, Singapore 469003^{3}School of Engineering and Computer Science, IUB, Dhaka, Bangladesh

Received 30 March 2011; Accepted 24 June 2011

Academic Editor: Youyun Xu

Copyright © 2011 Md. Zahangir Alam et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The paper proposes a novel technique for reducing noise in M-ary signal transmission through wireless fading channel using wavelet denoising that play the key role. The paper also explains that the conventional threshold-based technique is not capable of denoising M-ary quadrature amplitude modulated (M-QAM) signals having multilevel wavelet coefficients through wireless fading channels. A detailed step by step wavelet decomposition and reconstruction processes are discussed here to transform a signal function into wavelet coefficients using simulation software like MATLAB. A 16-QAM modulated symbol through a Rician fading channel is weighted by a control variable of complex form to force the mean of each detail coefficient except low frequency component to zero to enhance noiseless property. The bit error rate (BER) performance of the simulation results are furnished to show the effectiveness of the proposed technique. The root mean square of the deviation of the reconstruct signal from the original signal is used to express the effectiveness of the proposed technique. The traditional denoising provides very high value (above 90%) of the percentage root mean square difference (PDR) and the proposed technique provides only 10% PDR value for the symbol through a noisy channel. The result of the simulation study reveals that the BER performance can be increased using an appropriate control variable to force the mean of each detail coefficient to zero.

#### 1. Introduction and Literature Review

The signal transmitted through a wireless channel may arrive at the receiver through a number of paths with different amplitude, phase, and time delay due to multipath fading. The relative movement of either transmitter or receiver or both causes time varying multipath fading. Equalizer is used to minimize the fading effect, but the design complexity of the equalizer increases with the increase of data rate. Wavelet decomposition is in use in the area of image processing, image compression, segmentation, and denoising [1–5]. The wavelet coefficients are estimated from the noisy signals using maximum a posteriori (MAP) estimator [6]. In [6], four new bivariate distributions are proposed to wavelet coefficients of natural images to calculate the dependency between a coefficient and its parent. The wavelet denoising algorithm is extremely useful to reconstruct the original image signal from its noisy form. The main idea of denoising is to subtract the threshold value from all coefficients greater than the threshold one and set other coefficients to zero [7]. The wavelet transform of a signal is passed through a threshold that removes the coefficients below a certain value. Most of the energy of a signal is associated in the low-frequency region, and consequently the approximate coefficients are also contained in this region. Thus the high-frequency noise terms [8] are removed through this process. Wavelet denoising based on threshold scheme can be found in [9–11]. Wavelet soft-thresholding algorithm has been proposed for image denoising to reconstruct a noisy image in [12]. Fractal image denoising method is used to generate an original image from the noisy one discussed in [13]. A robust wavelet denoising-based estimator using a robust loss function is used to denoise a wireless signal through a Gaussian noisy channel [14]. Traditional wavelet denoising uses a threshold value [7] to reconstruct noisy wavelet coefficients by allowing actual noise coefficient if greater than the threshold and zero otherwise.

In this paper, the authors attempt to apply wavelet decomposition technique to generate wavelet coefficients of a noisy signal passing through a fading channel. The multipath fading channel has been modeled recently as reported in [15]. The reconstruction of the original symbol from the noisy transmitted symbol through the fading channel using wavelet denoising technique is discussed in detail in this paper. Wavelet decomposition and reconstruction methods are also derived in [16–18]. The wavelet decomposition and reconstruction techniques are discussed here to generate wavelet coefficients for a noisy symbol and reconstruct the signal from the coefficients by a computer simulation using a programming language like C, Java, or software like MATLAB. The threshold-based denoising technique is used to reconstruct the modulated symbol. The problem of the traditional denoising is that the wavelet coefficients of 16-QAM or 64-QAM modulated symbols consist of multiple-bits and the thresholding technique is not able to estimate the noisy coefficients either greater or less than the threshold value. The authors also propose, in this paper, a denoising technique for multiple-bits symbol using a complex weighting control variable to force the mean of each high-frequency detail coefficient to zero. The transmitted 16-QAM symbol through a Rician fading channel is received in MATLAB simulation environment. The received noisy symbol is transformed to generate wavelet coefficients based on the wavelet decomposition technique. The proposed M-ary denoising technique predicts the actual detail coefficients from the received noisy coefficients. The transmitted signal is reconstructed from the predicted wavelet coefficients according to the decomposition technique. Bit error rate (BER) performances of M-ary signals with additive white Gaussian noise (AWGN) through the Rician channel are calculated and graphically shown. The M-ary signal is also used in antenna diversity environment to obtain better BER performance. The simulation study reveals that the proposed technique provides better BER performance of M-ary signal for a multipath fading channel, depending on the weighting control variables. The PDR of the 16-WAM symbol through an AWGN channel for traditional denoising technique is approximately 82.5% at dB. The proposed technique provides PDR value of 17.3% of the same symbol through the same channel at the same value of . Thus the proposed technique provides excellent BER performance along with the antenna diversity technique that can be used in land mobile telephone system in 3G or 4G or in WiMAX physical layer to transmit better quality noise free audio or video signal at very high data rate.

The organization of the paper is as follows. In Section 2, the modeling of a multipath fading channel is discussed to characterize the wireless channel. Discrete wavelet transform technique is discussed in Section 3 to generate the wavelet coefficients of a 16-QAM modulated symbol through the fading channel, and wavelet reconstruction is also discussed in this section. Traditional denoising and its complexity for M-ay are discussed in Section 4. The M-ary symbols and its application in wireless signal processing are proposed in Sections 5 and 6. Results of the simulation study are shown in Section 7. Finally, the conclusion of the work is drawn in Section 8.

#### 2. Modeling of Wireless Channel

A QAM bandpass signal corresponding to the transmitted base-band signals and with angular frequency can be expressed as [19] The result of the signal in (1) transmitted through direct and reflected paths between the transmitter and receiver appears as shown in Figure 1. The reflected paths are the delayed versions of the signal at the receiver, and each path undergoes scattering and the locally scattered and reflected components combine at the receiver, the phenomenon being known as multipath fading. If the direct path of Figure 1 is absent, then the nonline-of-sight (NLOS) paths are characterized by Rayleigh distribution, and the line-of-sight (LOS) propagation paths are characterized by Rician distribution.

The carrier frequency is shifted from the actual one due to the relative movement of either transmitter or receiver or of both. The resultant received signal frequency is given by where is the Doppler frequency shift and is defined by where is the wavelength, is the vehicle speed, and is the angle between the direction of wave arrival and the line of the vehicle movement. Due to the effect of Doppler shift, the modulating signal will be [20] The received signal due to the superposition of say paths can be written as where, , , and are the attenuation factor and carrier phases. The fading path having a delay of , where is the th propagation path length and is the velocity of light. The received signal can be written in the form The received signal can be rewritten as where Hence, is the th set of complex path gains, is the excess delay of th path, and () is the delta function. For the Rician channel, the fading process can be written in the form [21] where is the Rician -factor of the th path, is the average envelope power for th path. The received complex signal consists of a large number of plane waves as in (7) and the received complex envelope can be written as Here, and are Gaussian random process with nonzero means and , respectively. If the random variables and have the same variance , then the received complex envelope has a Rician distribution of the form [22] here, is the direct signal amplitude and is the zero order modified Bessel function of the first kind and is defined as The Rician distribution of the received signal is shown in Figure 2.

#### 3. Wavelet Transform Algorithm

The wavelet transform provides the time-frequency representation of the signal through decomposition by passing the time domain signal through various high-pass and low-pass filters. Each filter gives output signal of various frequency bands. The wavelet coefficients after the decomposition process can be used for further processing to obtain the desired level of signal.

The coefficients are the function of scale and position and are the sum over the signal time multiplied by scale and shifted version of the wavelet function (). Filter technique is used to represent time-scale representation and when the signal passes through filters then its feature is changed by upsampling and downsampling operations. In downsampling operation some of the samples of the signal are removed and in upsampling some of the new samples are added. The wavelet transform is to compute a set of coefficients and each coefficient is the inner product of input function and basic wavelet [23]:
The continuous wavelet transform is defined [16] as:
where and represents the scale and position, _{ a,b} is defined as
The transform becomes
The discrete wavelet transform can be used to reconstruct the original signal from the wavelet coefficients [17] as:
where, .

The family of function is an orthonormal basis of the space of finite energy signals, the multiresolution analysis (MRA) of is a closed sequence [18] as Detail spaces are associated with the approximation spaces, and is orthogonal complement of in as As scaling function is associated with , the integer translate of the family of generates . Similarly, the integer translate of the family of generates . By replacing by , we can decompose into its components as where and can be written as [17] with Let us consider two discrete filters and with their impulse response, and .

Now let us define the following equations: Now, replace by , then we have The discarding of odd coefficient is called downsampling, and the decomposition algorithm is shown in Figure 3.

Similarly, the reconstruction of can be written as Now, considering the same discrete filters and with impulse responses , and , we have Setting odd entries to zero, we have We have, Now, writing , and , the sequence and are called upsample and the reconstruction formula is as:

The reconstruction step is shown in Figure 4.

Multiresolution analysis analyze the signal in different frequency band and enables to observe signal in both time and frequency domain. The three level multiresolution analysis of is shown in Figure 5.

The analysis gives the three detail coefficients , , and , and the approximate coefficient . Here is the half resolution of , is the half resolution of . These give the detail information of the signal, and the three levels multiresolution step can be written mathematically as where in all cases. The multiresolution wavelet reconstruction algorithm is shown in Figure 6.

The multiresolution reconstruction produces final three details (, , and ), and approximation . Mathematically , , , and can be written as where is the upsample of defined as The reconstructed signal can be written as

#### 4. Problem Statement

Wavelet denoising is used to predict approximate detail coefficients by allowing minimum error between the detailed coefficients of threshold noisy signal and the original. A noisy random signal is represented as where is the original signal and is the noise with mean zero. Now, decompose the noisy signal using three levels wavelet decomposition to generate three detail coefficients and one approximate coefficient. A threshold value is applied to each detail coefficient to allow the actual value of the noisy detail coefficient if it greater than the threshold and force to zero otherwise, mathematically can be written as [7] where is the threshold value of level. In first stage the threshold can be chosen by finding the minimum error between the detail of noisy and original signal. Later, a loop is applied for each level to find the corresponding threshold value that gives the minimum error. The noisy signal is denoised by using the threshold level to each detail coefficients after applying three level wavelet transform to generate three detail coefficients. Finally the signal is regenerated using the wavelet reconstruction technique. The result of the denoising of the noisy signal is shown in Figure 7.

In the denoising technique, let the detail coefficients have two levels a and b as in Figure 8(a). When noise affects the signal then the maximum and minimum level may be shifted due to addition of the noise. The threshold level c is used to allow the coefficient equal to b when the noisy coefficient is equal or greater than c, and equal to a if the noisy coefficient less than c. If the coefficients have multilevel for M-ary modulated signal, that is, another two level d and e beyond a and b then the denoising technique is unable to reconstruct the original signal because if the noisy coefficient is less than c and a then the technique forces the coefficient equal to a instead of d as shown in Figure 8(b).

#### 5. Proposed Multilevel Denoising Technique

A modulated signal (where ) is affected by noise and interference (where , and are any real value), and the resultant signal can be written as The mean of the white noise is zero, and the mean of the each detail wavelet coefficient is zero [8]. The standard deviation of detail coefficients () of at level (, ; is decomposition level) is as where is the length of the detail coefficients, and mean is defined as The standard deviation of detail coefficients of the noisy signal can be written as where is the components added with the original detail coefficient due the effect of noise. The mean of noisy detail coefficient can be written as To estimate the original detail coefficients, the mean of noisy detail coefficient should be zero, and the standard deviation can be written as where, .

A complex weighting vector , where and are any real number, is used to minimize random noise by a factor of in [15]. In this paper, we used the complex weighting vector to minimize the noise in the detail coefficients by the same amount to generate noise-free wavelet coefficient by approximating the mean of detail coefficient to zero. The modulated signal is multiplied by a complex variable , , , then the resultant signal can be written as . The noisy received signal can be written as If , then Here, the noise and interference term in (46) of the received signal through a fading channel is reduced by a factor of . The mean of the detail coefficients of less noisy signal in (46) is as The noise with each of the detail coefficients is due to the noise with the original signal . The complex number reduced the noise term associated with by a factor of . If and then the mean , and the complex number may be defined as mean control variable. The noise and interference of detail coefficient is reduced by a factor of , and the reconstructed signal from the noise free detail coefficients may be considered as free from noise and interference as the mean of each detail coefficient tends to zero.

#### 6. Proposed Denoising Method in Signal Processing

The noise and interference associated with the wavelet detail coefficient can be discarded using the proposed denoising technique as in (42) to (47). The transmitted signal is multiplied by the complex control variable . The mean of the detail coefficient is kept at a level approximately equal to zero by reducing the noise term by a factor of as according to (46). The value of the mean control variable is generated according to the required mean value of each detail coefficient. The transmitted signal can be obtained by using wavelet reconstruction as shown in Figure 9. In this figure the modulated symbol is transmitted through a noisy channel and the noisy signal at the receiver is processed by using three level wavelet decomposition according to (30). The transmitted modulated symbol is reconstructed by using wavelet reconstruction method of the denoised wavelet coefficients according to (31).

Wavelet denoising can be used in antenna diversity technique to obtain better quality signal through a noisy channel. In antenna diversity, signal is transmitted through more than one antenna and received by a single antenna, on the other way one antenna transmits signal and more than one antenna receive the signal. The first process is called transmit diversity and second is called receive diversity. A combiner is used at the receiver to add the entire incoming signal in phase or search the path having highest signal strength [24]. The wavelet denoising method is used to further discarding noise from combiner output signal, the entire process is shown in Figure 10. In transmit diversity, let the signal be sent from the transmitter to receiver through two diversity path and , the maximum ratio combiner (MRC) combines the two received fading signal , , according to [24], as In two branches transmit diversity, symbol from antenna zero can be denoted by and from antenna one by and in the next symbol period signal the transmitted signal can be represented as , and . The received signal for channel response and can be written as here , and are the noise signal for 1st and 2nd path. Now, we select the control variable of form , then the resultant transmitted symbol is The received signal can be expressed as Dividing (35) by , and considering , we have The received signal from other diversity path can be written with applying the same operation as first antenna, we have the received signal as The MRC combiner combined these two received as in (48) The received signal in (54) is then passed through the wavelet decomposition process to generate the detail coefficients. Here, the control variable () is adjusted to approximate the mean of each detail coefficient equal to zero to obtain better channel capacity. The entire process is shown in Figure 10. The PDR is used to examine the percentage of the improvement in the reconstructed signal from the original as: original 16-QAM symbol and : reconstructed 16-QAM symbol.

#### 7. Results

The simulation results in MATLAB environment of the BER performance of the AWGN channel for the 16-QAM symbol is shown in Figure 11. Here, the three-level wavelet decomposition technique is applied on the noisy received signal to generate the appropriate detail coefficients. The wavelet reconstruction technique is applied using (31) to generate the detail coefficients which ultimately generate the transmitted signal. The generated transmitted signal is found to be the same as the one received directly from the noisy channel. The multilevel denoising technique shown in Figure 9 is used to increase the performance of the channel. The BER of AWGN channel using different values of the control variable to approximate each detail coefficient equal to zero is shown in Figure 11 for different values of . The BER through the direct path is 0.2620 at dB and the wavelet decomposition provides the same BER performance at the same value. Now, the controlling value of is used for which, the BER value is 0.085 for the same parameter. The control variables , , and of values of , , and , respectively, to provide BER values of 0.082, 0.0630, and 0.0290 at of 3 dB. The mean of three details coefficients for the above condition is . The control variable allows the mean of the detail coefficients to tend to zero and this phenomenon reduces the noise and interference effects. The BER performance through a Rician fading channel is shown in Figure 12 using the same control variable as of Figure 11 with sampling frequency of 10 MHz, Doppler frequency shift of 100 Hz, and the Rician factor of 3.

The control variable in antenna diversity is used to reduce the noise level by a factor more than the previous case because diversity provides some noise reduction performance. Here, the BER performance in antenna diversity for the control values of values and is shown in Figure 13. The BER is 0.027 for two branch transmit diversity through the AWGN channel with dB. Now, the BER for the value of control variable and with wavelet denoising in the same branch transmit diversity through the same channel is 0.008 and 0.001, respectively, with considering the same value of . Here, the simulation is performed by sending a 16-QAM symbol through an AWGN channel with mean zero and variance . The phase offsets for the two paths are 30 and 35, respectively. The PDR of the AWGN channel is calculated using (55) by considering as transmitted 16-QAM symbol. The denoising technique in [7] is used to reconstruct , and we get the value of PDR which is above 90% for any value of . Similarly, the proposed denoising technique is used to reconstruct for the same symbol through the same channel. The PDR of the same channel for simply transmitting the symbol through the channel without any signal processing technique, and using proposed denoising technique at different values of control variables is shown in Figure 14 at dB. It is seen from the result that the proposed denoising technique provides the lowest deviation of reconstruct signal from the original one, because traditional denoising provides PDR of about 92%, direct transmission provides PDR of 82.5%, and the proposed multilevel denoising provides a PDR of 17.3% for a signal having multipoint in the constellation diagram with dB for all cases.

#### 8. Conclusion and Future Work

Three main issues regarding wireless signal denoising through a multipath fading channel are addressed in this paper. Firstly, the wavelet decomposition and reconstruction of a 16-QAM symbol through a noisy channel using MATLAB simulation software. The BER performance of a modulated signal through AWGN and Rician channel and the performance using wavelet denoising technique for the same signal through the same channel is the same as shown from Figures 11 and 12. It is seen from the result that the wavelet decomposition and reconstruction techniques discussed in this paper for the simulation study reconstruct the signal using the wavelet decompose coefficient with 100% accuracy. Secondly, the traditional threshold-based denoising technique is used to denoise 16-QAM modulated noisy wavelet detail coefficients, the traditional denoising method is not capable to generate noise free coefficients due to the fact that the one threshold level cannot predict a noisy coefficient from multilevel values as discussed in Section 4. Thirdly, a control variable based multilevel wavelet denoising technique is proposed to denoise a complex modulated symbol through a fading channel, and the BER performance is shown using appropriate control variable to force each detail coefficient to zero by reconstructing noise free transmitted symbol.

The simulation result shows that the BER performance is increased with increasing tendency of the mean of the detail coefficients close to zero by selecting appropriate control variable as shown in Figures 11 and 12. The proposed multilevel denoising technique is also used in antenna diversity to obtain a far better BER performance, and the simulation study shows that the proposed denoising technique along with the antenna diversity to provide far better BER performance than the one shown in Figures 11 and 12. The PDR of the proposed multilevel denoising technique shows that it provides accurate reconstruction of a complex modulated symbol through a multipath fading channel because it deviates the estimate signal from the parent one by only 10%.

In the future, the analysis can be done to search for far better a control variable depending on the condition of the wireless channel characteristics to obtain far more accurate results. The proposed technique can be used in multiple input multiple output (MIMO) case to send image signal at very high date rate with far better BER results through WiMAX physical layer.

#### References

- J. M. Shapiro, “Embedded image coding using zerotrees of wavelet coefficients,”
*IEEE Transactions on Signal Processing*, vol. 41, no. 12, pp. 3445–3462, 1993. View at Publisher · View at Google Scholar · View at Scopus - H. Choi and R. Baraniuk, “Multiscale texture segmentation using wavelet-domain hidden Markov models,” in
*Proceedings of the 32nd Asilomar Conference on Signals, Systems & Computers*, vol. 2, pp. 1692–1697, November 1998. View at Scopus - J. Liu and P. Moulin, “Image denoising based on scale-space mixture modeling of wavelet coefficients,” in
*Proceedings of the International Conference on Image Processing (ICIP '99)*, pp. 386–390, Kobe, Japan, October 1999. View at Scopus - M. K. Mihçak, I. Kozintsev, K. Ramchandran, and P. Moulin, “Low-complexity image denoising based on statistical modeling of wavelet coefficients,”
*IEEE Signal Processing Letters*, vol. 6, no. 12, pp. 300–303, 1999. View at Publisher · View at Google Scholar · View at Scopus - J. Portilla, V. Strela, M. J. Wainwright, and E. P. Simoncelli, “Adaptive Wiener denoising using a Gaussian scale mixture model in the wavelet domain,” in
*Proceedings of the IEEE International Conference on Image Processing (ICIP '01)*, pp. 37–40, October 2001. View at Scopus - L. Sendur and I. W. Selesnick, “Bivariate shrinkage functions for wavelet-based denoising exploiting interscale dependency,”
*IEEE Transactions on Signal Processing*, vol. 50, no. 11, pp. 2744–2756, 2002. View at Publisher · View at Google Scholar · View at Scopus - M. Alfaouri and K. Daqrouq, “ECG signal denoising by wavelet transform thresholding,”
*American Journal of Applied Sciences*, vol. 5, no. 3, pp. 276–281, 2008. View at Google Scholar · View at Scopus - D. Lee Fugal,
*Conceptual Wavelets in Digital Signal Processing*, Space & Signals Technologies LLC, 2006. - F. Abramovich, T. Sapatinas, and B. W. Silverman, “Wavelet thresholding via a Bayesian approach,”
*Journal of the Royal Statistical Society Series B*, vol. 60, no. 4, pp. 725–749, 1998. View at Google Scholar · View at Scopus - A. Chambolle, R. A. DeVore, N. Y. Lee, and B. J. Lucier, “Nonlinear wavelet image processing: variational problems, compression, and noise removal through wavelet shrinkage,”
*IEEE Transactions on Image Processing*, vol. 7, no. 3, pp. 319–335, 1998. View at Google Scholar · View at Scopus - M. S. Crouse, R. D. Nowak, and R. G. Baraniuk, “Wavelet-based statistical signal processing using hidden Markov models,”
*IEEE Transactions on Signal Processing*, vol. 46, no. 4, pp. 886–902, 1998. View at Google Scholar · View at Scopus - S. G. Chang, B. Yu, and M. Vetterli, “Adaptive wavelet thresholding for image denoising and compression,”
*IEEE Transactions on Image Processing*, vol. 9, no. 9, pp. 1532–1546, 2000. View at Google Scholar · View at Scopus - M. Ghazel, G. H. Freeman, and E. R. Vrscay, “Fractal-wavelet image denoising revisited,”
*IEEE Transactions on Image Processing*, vol. 15, no. 9, pp. 2669–2675, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Sardy, P. Tseng, and A. G. Bruce, “Robust wavelet denoising,”
*IEEE Transactions on Signal Processing*, vol. 49, no. 6, pp. 1146–1152, 2001. View at Publisher · View at Google Scholar · View at Scopus - Md. Z. Alam, C. Chittaranjon Patra, C. Patra, and M. Abdus Sobhan, “Modeling and performance analysis of a wireless fading channel,” in
*Proceedings of the 2nd ISECS International Coloquium on Computing, Communication, Control, and Management (CCCM '09)*, Sanya, China, August 2009. - L. Debnath and D. Bhatta,
*Integral Transform and Their Application*, Chapman & Hall/CRC, 2nd edition, 2007. - A. Boggess and F. J. Narcowich,
*A First Course in Wavelet with Fourier Analysis*, Prentice Hall, Upper Saddle River, NJ, USA, 2001. - G. Oppenheim,
*Wavelets and Their Applications*, ISTE Ltd., London, UK, 1st edition, 2007. - J. G. Proakis,
*Digital Communication*, McGraw Hill, New York, NY, USA, 3rd edition, 1995. - H. Schulze and C. Luders,
*Theory and Applications of OFDM and CDMA-Wideband Wireless Comunications*, John Wiley & Sons, London, UK, 2005. - M. C. Jeruchim, P. Balaban, and K. S. Shanmugan,
*Simulation of Communication Systems*, Kluwer Academic/Plenum, New York, NY, USA, 2nd edition, 2000. - G. K. Stuber,
*Principles of Mobile Communication*, Kluwer Academic Publishers, 2001. - S.-T. Bow,
*Pattern Recognition and Image Processing*, Marcel Dekker, New York, NY, USA, 2nd edition, 2002. - S. M. Alamouti, “A simple transmit diversity technique for wireless communications,”
*IEEE Journal on Selected Areas in Communications*, vol. 16, no. 8, pp. 1451–1458, 1998. View at Google Scholar · View at Scopus