International Journal of Digital Multimedia Broadcasting

VolumeÂ 2008Â (2008), Article IDÂ 523649, 9 pages

http://dx.doi.org/10.1155/2008/523649

## On the Efficient Broadcasting of Heterogeneous Services over Band-Limited Channels: Unequal Power Allocation for Wavelet Packet Division Multiplexing

CNIT-MCLab, Department of Electrical and Electronic Engineering, University of Cagliari, Piazza d'Armi, Cagliari 09124, Italy

Received 5 September 2007; Revised 2 November 2007; Accepted 4 December 2007

Academic Editor: MassimilianoÂ Laddomada

Copyright Â© 2008 Maurizio Murroni. 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

Multiple transmission of heterogeneous services is a central aspect of broadcasting technology. Often, in this framework, the design of efficient communication systems is complicated by stringent bandwidth constraint. In wavelet packet division multiplexing (WPDM), the message signals are waveform coded onto wavelet packet basis functions. The overlapping nature of such waveforms in both time and frequency allows improving the performance over the commonly used FDM and TDM schemes, while their orthogonality properties permit to extract the message signals by a simple correlator receiver. Furthermore, the scalable structure of WPDM makes it suitable for broadcasting heterogeneous services. This work investigates unequal error protection (UEP) of data which exhibit different sensitivities to channel errors to improve the performance of WPDM for transmission over band-limited channels. To cope with bandwidth constraint, an appropriate distribution of power among waveforms is proposed which is driven by the channel error sensitivities of the carried message signals in case of Gaussian noise. We address this problem by means of the genetic algorithms (GAs), which allow flexible suboptimal solution with reduced complexity. The mean square error (MSE) between the original and the decoded message, which has a strong correlation with subjective perception, is used as an optimization criterion.

#### 1. Introduction

Unequal error protection (UEP) is a channel coding technique used to increase the robustness of data that exhibit different sensitivities to transmission errors. This is often the case of digital multimedia compressed streams such as JPEG2000 [1] or MPEG [2]. Due to the extensive use of predictive and variable length codes, a compressed stream is in general more vulnerable to data losses and transmission errors, which can desynchronize the decoder causing spatial and temporal error propagation [3]. In broadcasting, feedback channel is not available, thus UEP relies on differentiated forward error correction (FEC) coding [4]: depending on their sensitivities to channel errors, data are protected with codes with higher or lower error correcting capabilities. Reed-Solomon (RS) or Turbo Codes (TC) are frequently used [5, 6], but also more performing techniques, based on Rate-Compatible (RC) codes [7], have been proposed by the research community. Unequal power allocation (UPA) is an alternative UEP technique which is deployed when, for several reasons, FEC coding is not efficient [8]. For broadcasting multiplexed communications (e.g., DVB, DAB), for instance, the available channel bandwidth per service is a key constraint and the use of FEC-based UEP schemes is barely suitable. In fact, FEC is a discrete nature coding scheme. It is subjected to some constraint which restricts the protection level (i.e., the code rate) only to a set of fixed values. Therefore, the overhead introduced by FEC codes can be a significant limitation for the efficient use of the bandwidth. On the other hand, UPA aims at distributing the available budget power over the parts of the stream, according to their sensitivities to channel error, to achieve improved final quality on transmitted data without any increase of the transmission bandwidth. Basically, UPA is performed by assigning different power weights to the data according to their â€śimportanceâ€ť (i.e., channel error sensitivities) within the stream: higher transmission power is assigned to more sensible data. As to this, UPA is a â€ścontinuousâ€ť process in the sense that weights are chosen in a real set with an accuracy which can be a priori selected and in theory infinite. Therefore, against FEC, UPA allows more flexibility in the protection of sensible data.

Wavelet packet modulation for orthogonally multiplexed communication was introduced as a promising technique to improve performance of conventional FDM and TDM schemes in both Gaussian and impulsive noises [9â€“11]. The properties of wavelet packets are exploited to embed data into waveforms which are mutually orthogonal both in time and frequency. Several studies conducted on this technology have shown that opportune design allows minimizing the energy of timing error interferences, which impair conventional TDM systems [10]. The overlapping bandpass nature of the transmission pulses (i.e., wavelets) allows better exploitation of the bandwidth respect to classical FDM [10], and it also intrinsically mitigates fading effects [12]. Moreover, due to the scalability of its structure, wavelet packets permit to multiplex data with different format (e.g., JPEG2000 and MPEG-2), therefore being a desirable choice for broadcasting heterogeneous services.

In this work, a UPA scheme for wavelet packet division multiplexing (WPDM) is proposed. UPA applied to WPDM consists on assigning different power to wavelet packets according to the importance of the message signals carried on. In other words, considering a generic bit pattern, individual bits are weighted differently taking the channel conditions (i.e., the signal-to-noise ratio (SNR)) into account and transmitted on separate wavelet packets. As to the optimization, we use the mean square error in the parameter domain

with
and being the transmitted and decoded parameter, respectively. The nontrivial
complexity of the problem does not allow closed-form analytical solution,
which, thus, has to be sought by numerical approach. In literature, solutions
based on the gradient algorithm have been proposed [8]. The complexity of such
optimization methods increases with the size (i.e., number of bits) of the
frame to be transmitted. In this work, we address UPA by exploiting the
potentialities of the *Genetic Algorithm* (GA) to reduce the computation complexity. The use of GA as to the
weights optimization is one of the novel aspects of this work. A genetic
algorithm [13] is a search technique used in computing to find true or
approximate solutions to optimization and search problem. GAs are
extensively used in literature in different application fields of communication
engineering such as network design, unicast, and multicast routing [14â€“16]. They allow finding
iterated numerical solution to complex problems with accuracy dependent on the
number of iterations selected. The major advantage of genetic algorithms is
their flexibility and robustness as a global search method. They can deal with
highly nonlinear problems and nondifferentiable functions as well as functions
with multiple local optima. They are also readily amenable to parallel
implementation, which renders them appropriate in real-time adaptive
communications, extensively used for reconfigurable broadcasting services.

Results show that the proposed UPA-WPDM scheme allows increasing resilience of data which exhibit different sensitivities to channel errors during their transmission over AWGN channel. The performance improvement in terms of quality achieved in the parameter domain (i.e., ) has been proved against an equally distributed WPDM- and FEC-based UEP systems, in the presence of similar bandwidth constraint. Moreover, the bandwidth gain for target quality (i.e., fixed ) at a fixed bit error rate has been evaluated beside UEP FEC-based techniques.

In the following section, an overview on the WPDM technology is given. Section 3 formally defines UPA for WPDM by describing in detail the weighting optimization procedure and the GA-based proposed solution. The performance of the proposed UPA-WPDM scheme on Gaussian channel is analyzed and compared to equally power distributed equivalent schemes and to channel coding UEP systems in Section 4. Conclusions follow in Section 5.

#### 2. Wavelet Packet Division Multiplexing

WPDM is a multiple signal transmission technique in which the message signals are waveform-coded onto wavelet packet basis functions for transmission. To define the wavelet packet basis functions, we refer to wavelet multiresolution analysis (MRA), the details of which can be found in a number of textbooks [17â€“23] and tutorial articles [24â€“31].

Let be a unit-energy real causal FIR filter of
length which
is orthogonal to its even translates; that is, ,
where is the Kronecker delta, and let be the (conjugate) quadrature mirror filter
(QMF), .
If satisfies some mild technical conditions [17, 31], we can use an iterative algorithm to find the function for an arbitrary interval .
Subsequently, we can define the family of functions , , in the following (binary) tree-structured
manner:
where .
For any given tree structure, the function at the *leafs* of the tree forms a *wavelet packet*. They have a finite
duration, ,
and are self- and mutually-orthogonal at integer multiples of dyadic intervals,
and hence they are a natural choice for scalable multiplexing applications [9, 10]. In Figure 1, the wavelet packet functions (a) and the relevant power
spectrum (b) for three-level (i.e., eight size wavelet packet) standard 12-tap
Daubechies filters decomposition [23].

In WPDM, binary messages have polar representation (i.e., ), waveform-coded by pulse amplitude modulation (PAM) of and then added together to form the composite signal . WPDM can be implemented using a transmultiplexer and a single modulator [10] as Figure 2 illustrates for a two-level decomposition. In this case, where , with being the set of terminal index pairs and the equivalent sequence filter from the th terminal to the root of the tree, which can be found recursively from (2). The original message can be recovered from using

An example of WPDM tree for a system that can be used for broadcasting heterogeneous services is shown in Figure 3(a). In this case, the transmission system uses two wavelet packets composed by two and four waveforms (i.e., wavelets), respectively. In Figure 3(b), the relevant subband structure is displayed: the total bandwidth is equally shared between the two packets, but a different partitioning (two against four) is implemented within each packet. Differently formatted streams can be transmitted by associating them to the appropriate wavelet packets.

#### 3. Unequal Power Allocation For WPDM

Without loss of generality to model, a generic bitstream exhibits different error sensitivities to channel conditions, we consider a discrete periodic (period ) memoryless source : and an analog to digital process AD: with , , being the LSB. Each is then multiplied with the specific weight of the diagonal matrix . The weighted bit pattern is then transmitted by a th order WPDM over a channel affected by additive white Gaussian noise (AWGN) with zero mean and variance . The signal at the receiver front end is with as in (1) and .

After demodulation, the distributed vector is , where , represents the demodulated noise along the signal message components (i.e., relevant noise). Following decision based on Maximum Likelihood (ML) criterion, the estimate is produced by inverse digital to analog (DA) process. A sketch of the system is depicted in Figure 4.

#### 3.1. Weight Optimization

Considering bipolar binary representation , if bits in are inverteddue to AWGN, a wrong decision is made at the receiver, thus producing a distortion . Aim of the optimization process is to calculate optimal weights in the sense of a minimized expected value . Assuming ergodicity, it is possible to calculate as follow: where are the different possible parameter values, the occurrence of the reproduction levels , and the transition probabilities between transmitted and received bit patterns. Due to the orthogonal properties of WPDM waveforms and to the independence of the noise samples, the transition probabilities are [4, 32]: By imposing , we can write [8] and impose the following constraint on the weights WPDM is based on binary amplitude modulation, thus, the bit error probabilities in (6) are [33]

Mathematically, the optimization problem is to minimize (5) under the constraint (7). In other words, UPA raises () the immunity to noise channel for more significant bits, paying as a counterpart lower robustness () on less significant one, to achieve average improved performance on the transmission of parameter in the sense of minimum expected distortion .

The complexity of the above optimization
problem, which increases with the size of frames ,
does not allow closed form solutions. Therefore, to identify the
solution, we use a numerical approach based on *Genetics Algorithms* (GAs).

#### 3.2. Genetics Alghoritms (GAs)

GAs are implemented as a computer simulation in which a population of
abstract representations (*chromosomes*)
of candidate solutions (*genes*) to an
optimization problem evolves toward better solutions. The evolution usually
starts from a population of randomly generated chromosomes and happens in
generations. In each generation, the fitness of every chromosome in the
population is evaluated, multiple chromosomes are stochastically selected from
the current population (based on their fitness), and modified (mutated or
recombined) to form a new population. The new population is then used in the
next iteration of the algorithm.

In
the proposed system, the chromosomesare defined as arrays of genes .
The range of possible values of is constrained by (7). An initial population of chromosomes is randomly selected. The
fitness function is as defined as in (5). Two operations are allowed to
determine the evolution of the initial population: *crossover* (with probability ) used
to interchange the elements of two chromosomes and *mutation* (with probability ) which modify the value of one or more genes within a chromosome with
the aim of leading the search out of local optima. In particular, the most
fitting part of the population is selected and directly inserted
in the new generation, while the rest of the population is
discarded and replaced by a subpopulation created by means of the crossover and
mutation operators. In the case of two identical chromosomes resulted after the
crossover and mutation operations, two individuals are randomly generated. The
termination condition is satisfied once either the algorithmreaches a selected number of iterations (IT) or the fitness function maintains the same value for iterations. At the end of the process, the
chromosome with low score in the fitness function (i.e., lower distortion on
the reconstructed frame) will be selected for the transmission.

Figure 5 gives an example of the crossover and mutation operations.

In this particular case, chromosomes are composed by four genes; at iteration the crossover operator swaps the first two genes of the chromosomes and as they were at iteration , whereas the mutation varies the chromosome by multiplying the second and fourth genes for the quantity with , respectively. The flowchart of the proposed GA is shown in Figure 6.

The accuracy of such approach is
strictly dependent on the values of IT and IT_{
MAX
}, whereas the complexity of the algorithm depends also on the
definition of chromosomes, on the size of the initial population and on
the and probabilities. Chromosomes are arrays of genes
which are real values. The higher the precision on the representation of the
genes (i.e., the number of decimal digits used to approximate real values), the higher the accuracy achieved by
the UPA, but also, the higher the complexity of the algorithm. Similarly, big-size
populations guarantee higher performance, but also lead to time consuming
processing. A critical matter is the selection of and probabilities: high values can determine
instability of the GA which could diverge, whereas, on the other side, low
values likely lead to slow convergence.

#### 4. Results

A WPDM system which deploys two packets of size is used to multiplex two streams having same rate, but different format (see Figure 7). Standard Daubechies minimum-phase scaling filters of length [31], which guarantee short delay and substantial capacity advantage over conventional FDM systems [10], are deployed. Without loss of generality, to model the parameters and delivered at time we use zero-mean () Gaussian sources and with unitary variance (). AD/DA processes deploy natural binary mapping based on uniform quantizers.

At first, we
have run some preliminary tests to analyze the importance of the GA parameters.
The crossover operator was allowed to interchange genes whereas
the mutation occurred on genes,
being the operator which produces the
integer part of the argument. In other words, at each iteration, a maximum of
40% of the chromosome parents could appear on the next generation of
chromosomes and only 10% of a chromosome could vary. According to this, was varied in the range , and in the range 0.3*Ă·*0.7 and 0.01*Ă·*0.3, respectively. Finally for mutation varied within the range . The maximum difference in terms of
fitness function value among all the solutions was observed to be less than 5%.
Therefore, the following considerations can be made: huge-size populations bring
to better solutions at the expense of a higher-processing time; the probability is suggested to be set equal to or higher than 0.1, whereas above to avoid an excessive number of iterations; the probability does not sort significant effects in the range used. As to
the outcome from the preliminary tests on GA behaviour applied to the UPA
problem, in the following experiments the population was composed by chromosomes, eight decimal digits were used to
represent genes (i.e., ), the probability , and , whereas and .
For the sake of clearness, Table 1 summarizes the parameter setting for the
experiments.

Achieved quality in the parameters domain is expressed in terms of the signal-to-noise ratio () measured in decibel with as in (1). is evaluated at varying average bit error probabilities with as in (8).

We have compared the proposed UPA with
a benchmark equal power allocation (EPA) WPDM system and an UEP scheme based on
FEC coding. In the latter system, we have deployed Reed-Solomon (RS) codes [33]. RS codes are *nonbinary
cyclic* codes with symbols made up of -bit sequences, where is
any positive integer having a value greater than 2. codes on -bit symbols exist for all and for which where is the number of data symbols
being encoded, and is the total number of code symbols in the encoded
block. The error-correcting capability of the generic code is .
UEP is implemented by protecting
data with codes with higher- or lower-code rate . At varying the channel error rate, for every
WPDM channel, an appropriate code is
selected for data protection according to the sensitivity to channel errors of
the data carried on. More significant data (e.g., MSB) are protected by codes
with higher error-correcting capabilities (i.e., higher-code rates). In
particular, for any average error rate the
optimization procedure aims at selecting the codes so
that the is minimized under the bound of constant average code rate .

For our experiments, we have selected and which corresponds to an increase of the total bandwidth of about 16%. To reduce the complexity of the coding process, we have fixed the number of code symbols in the encoded block . The average error correcting capability of the system is therefore symbols per codeword. In other words, on the average, such a scheme is able to correct up to 3 symbols that contain errors in a codeword. Tables 2 and 3 report the details (i.e., actual code rate and error correcting capability ) of the codes used at for the transmission of and , respectively.

In Figures 8 and 9, we refer to UEP RS-based coding as RS (38,32). The analysis of the graphics reveals that UPA outperforms EPA along all the variation ranges of the average bit error probability within the transmitted frame with a peak gain of 6.84?dB at in case of . Same behaviour is noticeable with respect to RS coding for , with 3.57?dB the peak gain for and for . For all the systems perform similarly with slight prevalence of the RS coding which is more evident for . Superior performance in case of transmission can be justified by the higher precision obtained by a finer power distribution performed with eight weights with respect to a coarser allocation based on only four weights as for . More generally, the UPA prevalence is due to the capability of the optimization procedure to obtain high accuracy by selecting weights in a range of real values.

In Figures 10 and 11 show how for severe channel conditions the weights relevant to higher significant bits (i.e., , and ) are emphasized with respect to all the others. For approaching a decrease of the above weights corresponds to an increase of and which become also higher than 1. For all the weights converge to equal unitary value, but still remaining slightly different for .

Figure 12 shows the percentage bandwidth gain achieved by UPA with respect to UEP based on RS coding for target quality (i.e., fixed ) on the transmitted parameters and , at fixed , for the WPDM system used for experiments as represented in Figure 5. A minimum bandwidth gain above 20% is noticeable whereas similar high variations are observed in both cases. This is due to the discrete nature of RS codes, which are constrained to only a definite set of possible code rates. On the other hand, UPA is a continuous process which guarantees more flexibility in the protection of sensitive data.

In order to assess the suitability of the proposed scheme for real applications, such as audio and video broadcasting, as a further test, we have considered the specific multiplexed transmission of a standard image and a stereo-audio sequence. Referring to the system proposed in Figure 7, we have used the well-known image â€śLenaâ€ť of size in RGB format coded at 8?bpp per color component (see Figure 13), as a transmission source . We have measured the quality on the reconstructed image by standard PSNR metric expressed in decibel. On the other hand, we have ripped a 5 seconds from a stereo-audio CD signal sampled at 44.1?KHz coded at 16?bps and used as a source . For the evaluation of the quality on the received audio signal, we have used the perceptual evaluation of audio quality (PEAQ) strategy [34]. PEAQ is technique recommended by the ITU, which evaluates the quality of an audio signal by a single number, called objective difference grade (ODG), which varies within a range , with 0 the highest quality score. PEAQ has proven to be more performing than conventional metrics based on mean square error on the evaluation of the performance of the conventional audio codecs [34].

Table 4 shows the results achieved in case of . The quality on the reconstructed image is slightly below 30?dB, whereas the PEAQ measured on the received audio sequence is just up . This result is in line with the typical performance of low-bit-rate audio and video codecs. For the transmission of audiovideo at a rate of 64?Kbit/s, MP3/MPEG-4 codecs achieve PSNR approaching 30?dB for the reconstructed frames [2] and PEAQ of around for the audio sequence [35]. Since conventional DAB and DVB broadcasting systems work at the proposed system could be an alternative solution for the broadcasting of multimedia heterogeneous contents in case of extremely hard transmission condition, when only little quality requirements are set.

#### 5. Conclusion

In this work, we have presented an orthogonal multiple transmission system based on wavelet packet modulation suitable for the resilient broadcasting of data which demonstrate different sensitivities to transmission errors. A novel unequal error protection technique based on differentiated allocation of the transmission power over the modulated waveforms allows improving the final quality of the received parameters in case of AWGN channel, without any increase of the transmission bandwidth. The optimization of the weights has relied on Genetic Algorithms which allowed to achieve reduced complexity. Due to its scalability properties, the proposed scheme is able to provide for multiple transmissions of heterogeneous services which can be independently protected according to their specific format. Therefore, unequal power allocation applied to wavelet packet division multiplexing offers improved flexibility to broadcaster. Nevertheless, it is worthy to point out that particular attention has to be given to the design of the wavelet filters which are real-valued and under the approach of UPA could impair the performance of the transmission in case of wireless systems. In fact, the proposed UPA scheme may increase the dynamic range of the input signals to the WPDM modulator in Figure 4. Since is real-causal FIR filter, the bigger input amplitude range may increase the complexity of these filters. This may be a disadvantage of UPA for implementation.

Future work on this subject will investigate the capability of the proposed scheme to deal with real-time varying transmission conditions including the presence of fading effects and the broadcasting of reconfigurable heterogeneous services.

#### Acknowledgment

The author would like to thank the three anonymous referees for their constructive comments and suggestions.

#### References

- D. S. Taubman and M. W. Marcellin,
*JPEG2000: Image Compression Fundamentals, Standards and Practice*, Kluwer Academic, Norwell, Mass, USA, 2002. - F. Pereira and T. Ebrahimi,
*The MPEG-4 Book*, IMSC Press Multimedia Series, Prentice-Hall, Englewood Cliffs, NJ, USA, 2002. - K. R. Rao and J. J. Hwang,
*Techniques & Standard for Image Video & Audio Coding*, Prentice-Hall, Upper Saddle River, NJ, USA, 1996. - J. G. Proakis,
*Digital Communications*, McGraw-Hill, New York, NY, USA, 4th edition, 2001. - A. Natu and D. Taubman, â€śUnequal protection of JPEG2000 code-streams in wireless channels,â€ť in
*Proceedings of IEEE Global Telecommunications Conference (GLOBECOM '02)*, vol. 1, pp. 534â€“538, Taipei, Taiwan, November 2002. - B. A. Banister, B. Belzer, and T. R. Fischer, â€śRobust image transmission using JPEG2000 and turbo-codes,â€ť
*IEEE Signal Processing Letters*, vol. 9, no. 4, pp. 117â€“119, 2002. View at Publisher Â· View at Google Scholar - X. Pan, A. Cuhadar, and A. H. Banihashemi, â€śCombined source and channel coding with JPEG2000 and rate-compatible low-density parity-check codes,â€ť
*IEEE Transactions on Signal Processing*, vol. 54, no. 3, pp. 1160â€“1164, 2006. View at Publisher Â· View at Google Scholar - T. Brüggen and P. Vary, â€śUnequal error protection by modulation with unequal power allocation,â€ť
*IEEE Communications Letters*, vol. 9, no. 6, pp. 484â€“486, 2005. View at Publisher Â· View at Google Scholar - A. R. Lindsey, â€śWavelet packet modulation for orthogonally multiplexed communication,â€ť
*IEEE Transactions on Signal Processing*, vol. 45, no. 5, pp. 1336â€“1339, 1997. View at Publisher Â· View at Google Scholar - K. M. Wong, J. Wu, T. N. Davidson, and Q. Jin, â€śWavelet packet division multiplexing and wavelet packet design under timing error effects,â€ť
*IEEE Transactions on Signal Processing*, vol. 45, no. 12, pp. 2877â€“2890, 1997. View at Publisher Â· View at Google Scholar - K. M. Wong, J. Wu, T. N. Davidson, Q. Jin, and P.-C. Ching, â€śPerformance of wavelet packet-division multiplexing in impulsive and Gaussian noise,â€ť
*IEEE Transactions on Communications*, vol. 48, no. 7, pp. 1083â€“1086, 2000. View at Publisher Â· View at Google Scholar - T. N. Davidson, A.-J. Schott, and K. M. Wong, â€śBranch-hopped wavelet packet division multiplexing,â€ť in
*Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '98)*, vol. 6, pp. 3233â€“3236, Seattle, Wash, USA, May 1998. View at Publisher Â· View at Google Scholar - D. Whitley, â€śA genetic algorithm tutorial,â€ť
*Statistics and Computing*, vol. 4, no. 2, pp. 65â€“85, 1994. View at Publisher Â· View at Google Scholar - L. Atzori and A. Raccis, â€śNetwork capacity assignment for multicast services using genetic algorithms,â€ť
*IEEE Communications Letters*, vol. 8, no. 6, pp. 403â€“405, 2004. View at Publisher Â· View at Google Scholar - C. W. Ahn and R. S. Ramakrishna, â€śA genetic algorithm for shortest path routing problem and the sizing of populations,â€ť
*IEEE Transactions on Evolutionary Computation*, vol. 6, no. 6, pp. 566â€“579, 2002. View at Publisher Â· View at Google Scholar - R. Bhattacharya, P. Venkateswaran, S. K. Sanyal, and R. Nandi, â€śGenetic algorithm based efficient routing scheme for multicast networks,â€ť in
*Proceedings of IEEE International Conference on Personal Wireless Communications (ICPWC '05)*, pp. 500â€“504, New Delhi, India, January 2005. View at Publisher Â· View at Google Scholar - I. Daubechies,
*Ten Lectures on Wavelets*, SIAM, Philadelphia, Pa, USA, 1992. - A. N. Akansu and R. A. Haddad,
*Multiresolution Signal Decomposition: Transforms, Subbands, Wavelets*, Academic Press, Boston, Mass, USA, 1992. - C. K. Chui,
*An Introduction to Wavelets*, Academic Press, San Diego, Calif, USA, 1992. - P. P. Vaidyanathan,
*Multirate Systems and Filter Banks*, Prentice-Hall, Upper Saddle River, NJ, USA, 1993. - G. G. Walter,
*Wavelets and Other Orthogonal Systems with Applications*, CRC Press, Boca Raton, Fla, USA, 1994. - M. Vetterli and J. Kovacevic,
*Wavelets and Subband Coding*, Prentice-Hall, Englewood Cliffs, NJ, USA, 1995. - G. Strang and T. Nguyen,
*Wavelets and Filter Banks*, Wellesley-Cambridge Press, Wellesley, Mass, USA, 1996. - G. Strang, â€śWavelets and dilation equations: a brief introduction,â€ť
*SIAM Review*, vol. 31, no. 4, pp. 614â€“627, 1989. View at Publisher Â· View at Google Scholar Â· View at MathSciNet - O. Rioul and M. Vetterli, â€śWavelets and signal processing,â€ť
*IEEE Signal Processing Magazine*, vol. 8, no. 4, pp. 14â€“38, 1991. View at Publisher Â· View at Google Scholar - B. Jawerth and W. Sweldens, â€śAn overview of wavelet based multiresolution analyses,â€ť
*SIAM Review*, vol. 36, no. 3, pp. 377â€“412, 1994. View at Publisher Â· View at Google Scholar Â· View at MathSciNet - A. Cohen and J. Kovacevic, â€śWavelets: the mathematical background,â€ť
*Proceedings of the IEEE*, vol. 84, no. 4, pp. 514â€“522, 1996. View at Publisher Â· View at Google Scholar - N. Hess-Nielsen and M. V. Wickerhauser, â€śWavelets and time-frequency analysis,â€ť
*Proceedings of the IEEE*, vol. 84, no. 4, pp. 523â€“540, 1996. View at Publisher Â· View at Google Scholar - K. Ramchandran, M. Vetterli, and C. Herley, â€śWavelets, subband coding, and best bases,â€ť
*Proceedings of the IEEE*, vol. 84, no. 4, pp. 541â€“560, 1996. View at Publisher Â· View at Google Scholar - S. G. Mallat, â€śA theory for multiresolution signal decomposition: the wavelet representation,â€ť
*IEEE Transactions on Pattern Analysis and Machine Intelligence*, vol. 11, no. 7, pp. 674â€“693, 1989. View at Publisher Â· View at Google Scholar Â· View at Zentralblatt MATH - I. Daubechies, â€śOrthonormal bases of compactly supported wavelets,â€ť
*Communications on Pure and Applied Mathematics*, vol. 41, no. 7, pp. 909â€“996, 1988. View at Publisher Â· View at Google Scholar Â· View at Zentralblatt MATH Â· View at MathSciNet - A. Papoulis and S. U. Pillai,
*Probability, Random Variables, and Stochastic Processes*, McGraw-Hill, New York, NY, USA, 2002. - S. Benedetto and E. Biglieri,
*Principles of Digital Transmission: With Wireless Applications*, Kluwer Academic, Norwell, Mass, USA, 1999. - Recommendation ITU-R BS.1387-1, 11/01, â€śMethod for objective measurements of perceived audio quality (PEAQ)â€ť.
- M. Salovarda, I. Bolkovac, and H. Domitrovic, â€śEstimating perceptual audio system quality using PEAQ algorithm,â€ť in
*Proceedings of the 18th International Conference on Applied Electromagnetics and Communications (ICECom '05)*, vol. 1, pp. 1â€“4, Dubrovnik, Croatia, October 2005.