Research Article  Open Access
Sequentially Adapted Parallel Feedforward Active Noise Control of Noisy Sinusoidal Signals
Abstract
A large class of acoustic noise sources has an underlying periodic process that generates a periodic noise component, and thus their acoustic noise can in general be modeled as the sum of a periodic signal and a randomly fluctuating signal (usually a broadband background noise). Active control of periodic noise (i.e., for a mixture of sinusoids) is more effective than that of random noise. For mixtures of sinusoids in a background broadband random noise, conventional FXLMSbased single filter method does not reach the maximum achievable Noise Attenuation Level (). In this paper, an alternative approach is taken and the idea of a parallel active noise control (ANC) architecture for cancelling mixtures of periodic and random signals is presented. The proposed ANC system separates the noise into periodic and random components and generates corresponding antinoises via separate noise cancelling filters, and tends to reach consistently. The derivation of is presented. Both the separation and noise cancellation are based on adaptive filtering. Experimental results verify the analytical development by showing superior performance of the proposed method, over the singlefilter approach, for several cases of sinusoids in white noise.
1. Introduction
Active noise control (ANC) is a technique of cancelling acoustic noise by generating an appropriate antinoise signal using loudspeakers, and directing it towards the region where noise cancellation is required. Rapid progress in digital signal processors (DSPs) and sensor technologies as well as new realtime adaptive control algorithm designs [1â€“3] has opened the door to new approaches in ANC. Active control of sinusoidal signals has been an important research and development topic [4] and can lead to different system architectures [5] and design simplifications [6]. The standard approach to realize an ANC system is to employ the FXLMS algorithm in feedforward configuration [4, 7] as portrayed in Figure 1. denotes the primary path and denotes the secondary path. From Figure 1, it follows that the exact solution to the noise cancelling filter (the control filter) is given by Since the stable transfer function is in general nonminimum phase, should in general be a stable but noncausal system, and hence unrealizable unless is perfectly predictable. The noncausality restriction no longer applies to a periodic signal since its past and future values can be determined from the knowledge of the period or alternatively from the frequencies, amplitudes, and phases of the constituent sinusoids. Hence a periodic noise signal can ideally be cancelled perfectly by designing an appropriate (predictive) cancelling filter in the ANC system. For linearly modeled ANC system with stationary signals, the causality restriction places an upper limit on the achievable noise attenuation level. The noise attenuation level (NAL) is defined as the ratio of noise power before cancellation to the noise power after cancellation. Computation of the maximum achievable NAL for some given stationary noise process is derived in this paper. The maximum achievable NAL is denoted by . For an additive mixture of periodic and random signals, the can be calculated by assuming that the periodic component is cancelled completely. However experience with the FXLMS algorithm, using single cancelling filter, for additive mixtures of periodic and random signals shows that the NAL achieved is usually less than the . In this paper an ANC structure for such mixtures, that can achieve performances much closer to , is introduced. The notion of separating the periodic and random noise components and generating corresponding antinoise components is developed. Due to the simultaneous and additive way in which the two antinoise components are generated and combined, the proposed method is termed as parallel feedforward ANC. It should be noted that in this paper, the adaptation scheme is sequential, that is, the two noise cancelling filters are computed by a twostep adaptation procedure.
This idea, which was initially introduced in [8] in an empirical fashion, is given a rigorous analytical flavor with supporting experimental results in this paper. The proposed separationbased parallel ANC system has three stages: the first stage estimates the sinusoidal parameters and separates the periodic part from the random part of the noise signal using adaptive noise cancelling technique. The second and third stages operate simultaneously, generating the corresponding antinoise signals for periodic part and random part, respectively. Experimental results using different mixtures of noise signals show improved performance of the proposed ANC method.
Active noise control of a mixture of sinusoids with or without harmonic relation is a wellstudied topic. An indepth study of the convergence properties of the FXLMS algorithm when the noise to be controlled is a harmonically related mixture of sinusoids is presented in [9]. The idea of using a parallel bank of control filters, each handling a subset of the sinusoids is presented and analyzed. The more general case of active noise control of nonharmonically related sinusiods is analysed in a novel manner in [10] by employing a characteristic equation derived from an equivalent statespace model. A faster converging version of the FXLMS algorithm is also derived. A rootlocusbased convergence analysis of an FXLMSbased ANC system for nonharmonically related sinusoids is done in [11]. A phase tuning approach based on an integrative internal model principle is developed in [12] for active noise control of sinusoidal disturbances. It should be noted that the methods developed in [9â€“12] are applicable to pure sinusoidal mixtures without any broadband random noise.
The purpose of this paper is to introduce a new active noise control (ANC) architecture for effective active noise control of sinusoidal mixtures in a random background noise. The key idea is that for sinusoidal signals in random noise, it is possible to separate the sinusoidal and random parts and generate separate antinoise signals in a parallel fashion. Since periodic signals can be expressed as a mixture of sinusoids, the proposed method is also applicable to periodic signals with an additive random noise. The effectiveness of noise attenuation is expressed interms of the Noise Attenuation Level (NAL), which measures the relative level of attenuated noise with respect to the actual noise. The notion of, and an exact formula for, maximum achievable Noise Attenuation Level () for sinusoidal signals in white noise is presented. The proposed ANC system design is shown to achieve performances close to . The key ideas communicated are as follows:
(1)Given a feedforward ANC setup (one reference microphone and one error microphone setup) and a stationary acoustic noise input, there exists a maximum achievable noise attenuation level () if linear estimation with an MMSE criterion is employed. The method to calculate this upper limit is developed for sinusoidal mixtures in random noise.(2)In the active noise control of a mixture of sinusoids and a random noise with a single noise cancelling filter (or control filter), the adaptive algorithm finds an optimum filter that is neither equal to the optimum filter when sinusoids alone are present nor equal to the optimum filter when random noise alone is present. Alternatively, assuming that the sinusoidal and random parts can be separated, better cancellation can be achieved by formulating two mean square errors (MSEs) criteriaone for the sinusoidal component and the other for the random component.(3)For the periodic signal part, perfect cancellation can be achieved by equalizing the nonminimum phase secondary path with a noncausal equalizer in a causal manner since for periodic signals, the future values can be determined from the knowledge of the period or alternatively from the frequencies, amplitudes, and phases of the constituent sinusoids. For the random part the best cancellation achieveable is given by (4)For a mixture of sinusoids and random noise, it is possible to separate the periodic and random components and have a parallel ANC architecture. This architecture results in two separate control filters, each of which is optimum for the respective component and hence the overall system can reach better than the conventional singlefilter configuration. It seems that such an approach has not been considered in the literature so far.The paper is organized as follows. The notational framework and various definitions are introduced in â€śNotations and Definitionsâ€ť section, towards the end of the paper. Section 2 introduces the proposed parallel feed forward ANC architecture. Section 3 develops the methodology to calculate . Section 4 presents and discusses the experimental results. Concluding remarks are made in Section 5. The reader is advised to peruse the â€śNotations and Definitionsâ€ť section before proceeding to read Sections 2 and 3.
2. The Parallel Feedforward ANC Method (PFANC)
2.1. Motivation
The key idea of this paper is that for sinusoidal plus random noise mixtures, separating the sinusoids and random noise and generating antinoise signals separately improves the performance of the ANC system. The motivation for this approach lies in the observation that the optimum control (cancelling) filters for the periodic part and the random part are not the same. Assuming perfect separation of into and , and separate generation of antinoise signals for the periodic and random parts, the expression for the composite antinoise signal generated by the ANC system is where and are the noise cancelling filters of and , respectively, and denotes the convolution operation.
The best noise cancellation is achieved when and are optimum in the MMSE sense. The superscript Â° denotes the optimality of any solution in the rest of the manuscript. From (1), it is obvious that for the tones at , the optimum filter is any such that
For the random part , the optimum solution is the causal Wiener filter (whose exact form is shown later in (51), approximating the noncausal function .
Let denote the residual noise for the proposed parallel antinoise generation approach. Then assuming perfect separation of into and , is given by
Since can be written as the sum of and , (4) can be written as Thus consists of two components: from the sinusoidal component and from the random component. It should be noted that when , there is perfect cancellation of , that is, . The best cancellation of occurs when with the corresponding minimum error . Since and are statistically independent, the MMSE is
where is the expectation operator.
Let denote the residual noise for the singlefilter approach. The antinoise is given by The expression for is then
The best cancellation (in the MMSE sense) is achieved by minimising with respect to . Let the MMSE solution be with the corresponding error In general, the MMSE solution, will be different from both and . Since only will result in an error sequence with no tonal component, will always result in a residual tonal component in . The residual tonal component is denoted by and the residual random noise component is denoted by . Thus (9) can be written as where and .
The expression for MMSE owing to statistical independence is then
It should be recalled that the best attenuation of the random component is attained when equals . Since it is not the case in general, it follows that Since and are positive quantities, it follows that Since from (6), the MMSE for the proposed method , it follows from (13) that Thus the MMSE of the proposed method is always lesser than the MMSE of singlefilter method.
It should also be noted that when a single adaptive filter is used, the adaptive algorithm tends to cancel the dominating signal component, that is, the tones more effectively. As a result, the random component may be attenuated to a lesser degree, unaffected, or even reinforced. In fact, the adaptive filter in the singlefilter case shows a nonWiener behavior [13], the analysis of which is nontrivial. In contrast, the proposed ANC structure offers a simple means of enhancing performance.
2.2. Parallel Feedforward ANC (PFANC): Problem Formulation
The parallel feedforward ANC is formulated, as two subproblems, as follows.
(1)Given an acoustic noise signal , devise a method of separating into and , where is periodic, and is random.(2)Design an ANC system that has two parallel adaptive noise cancelling filters and for generating antinoise signals separately. It should be noted that, for analytical purposes, we consider a very short convergence time and hence drop the time indices that should accompany and , owing to their timevarying (adaptive) nature. These two cancelling filters are parts of and , respectively, (see Figure 2). The overall system will monitor the error microphone output to adjust the weights of these cancelling filters.In the proposed approach, the first subproblem is addressed in three steps: a frequency analyzer block which estimates the sinusoidal frequencies, a synthesizer block which generates a reference signal based on the estimated frequencies, and an adaptive noise canceller block which separates into and . The second subproblem is solved by setting up two FXLMSbased ANC systems. The method is elaborated in Figure 2. A more detailed explanation of the proposed approach is presented in subsequent sections.
2.3. PFANC Method
The proposed PFANC method can be divided into three major parts.
(1)Estimation of sinusoidal parameters(2)Separation of sinusoidal and random components of the noise signal.(3)Generation of antinoise signals corresponding to the sinusoidal and random signal components.2.3.1. Estimation of Sinusoidal Parameters
This part estimates the sinusoidal parameters from the measured input noise signal . Let samples of be observed. The estimation problem can then be stated as
Given that can be modeled as
and given observed samples of estimate , , . The determination of is implicit in the problem.
The above estimation problem is an important and a wellstudied one [14]. Since there are many methods available, and since this is an important step in the proposed scheme, the tools and techniques that were found to be particularly effective are described. The first step is detecting the number of sinusoids, which is essentially an order estimation problem. The next step is estimation of , , . An informationtheoretic approach [15] for the first step and a subspacebased high resolution parametric approach (ESPRIT) [14] for the second step are prescribed.
Detection of L
The first step is the detection of , the number of sinusoids, which amounts to order estimation. The gist of all the order estimation methods is defining an information measure as a function of the number of sinusoids. The number of sinusoids that best explains the observations, minimizes the information measure. Out of the host of order estimation techniques available [16], the method described in [15] is particularly effective for the sinusoids in noise scenario. The information measure is based on the Minimum Description Length (MDL) principle introduced in [17]. The reasons for adopting this method are twofold: (1) the method has been shown to be consistent, that is, it asymptotically approaches the optimum MAP rule [15, 17] and (2) there is no need to choose a penalty function [15].
The following quantities need to be defined, before introducing the technique.
In general for any quantity , denotes the estimated value of the quantity whereas the uncapped denotes the true value. For example, if is the vector of frequencies. then is the vector of estimated frequencies. Let us define the vectors , , , , and matrix is the , that is, order covariance matrix of the observations . Superscripts H and denote Hermitian and transpose operations, respectively.
However in practice, the covariance matrix of order is calculated from observations () of , that is, from . The observations data length should be greater than the covariance matrix order so that an accurate estimate of can be computed. The expression for the sample covariance matrix is [14] where, The eigenstructure of has all the information about the angular frequencies . Additionally let be the eigenvalues of sorted in decreasing magnitude. The number of sinusoids is then determined from the (for complex exponentials) that minimizes [15] The number of sinusoids can be found from , since for real sinusoids there are complex exponentials.
The minimum is found by computing the for values of ranging from . The eigenvalues are computed from an eigenvalue decomposition of . An example plot of for is shown in Figure 4.
Estimation of
The estimation of can be done using various ways like MUSIC, ESPRIT, FFT, and so forth, [14]. The method which has been shown to have very good statistical performance [14] is the nonlinear least squares (NLS) method in which the optimal values of minimizes the least squares measure,
In fact, for Gaussian noise scenarios the NLS is the best (in terms of minimum variance) estimator that achieves the CramerRao bound. The estimation accuracy of NLS is given by
For ANC problems, could be allowed to be reasonably large and as will be shown later, the PFANC method demonstrates superior performance for reasonable SNR values (â€‰dB). For example, for â€‰dB and , will be of the order of . Thus, very precise estimates can be expected from the NLS method. The drawback is that the minimization in (20) is a complex multidimensional search problem and could be computationally intractable. This drawback can be addressed by using the subspace methods, whose estimation performance is close to the best performance shown in (21) and is computationally efficient [14]. Out of the multitude of subspace methods (like MUSIC, ESPRIT, and MinNorm), ESPRIT has been shown to have superior statistical performance [14, 18]. Moreover, ESPRIT does not give spurious frequencies once the correct order is determined using (19).
ESPRIT Calculations
For completeness of discussion, the steps in frequency estimation using the ESPRIT method are as follows. (1)Calculate the sample covariance matrix using (17).(2)Calculate the eigenvalue decomposition of
where, (3)Let be the matrix composed of the eigenvectors corresponding to the eigenvalues ,
Define
(4)The frequencies are estimated by solving the system of linear equations
whose least squares solution is given by
The argument of the complex eigenvalues of gives the frequencies . Let be the eigenvalue of . Then the frequency estimates are
The amplitudes and phases can be estimated as
where , . The matrix is calculated from the estimated frequencies â€™s according to
indicates that the first samples of the data frame are used for ESPRIT calculations.
The estimated frequencies are then used to generate a reference periodic signal shown in Figure 2 using The estimates and can be further improved, by estimating from , by using the adaptive cancellation technique as explained in the next section.
2.3.2. Signal Separation
This block essentially estimates from in (30) by trying to cancel from (i.e., by estimating from ). This step adaptively finetunes the parameters from (30) to their actual values and also separates into and . It can be easily seen from Figure 2 that after convergence, the error signal is and the cancelling signal becomes , thus achieving separation. The update equation for the weight vector is given by An initial portion of the reference signal is used to estimate the tones and achieve component separation. The separated signal components will be denoted by and to indicate that they are estimates. For high , the estimates are highly accurate. By incorporating the fact that the effect of an arbitrary transfer function on a sinusoid can be modeled by two taps [4, 9], computationally efficient implementations of the above separation scheme can be derived.
2.3.3. Generation of Antinoise Signal
This block generates the antinoise (i.e., the cancelling) signal in two steps via two parallel ANC systems, and . uses as the reference and generates an antinoise which attempts to cancel . , after passing through becomes . The error signal is given by where denotes the desired noise to be cancelled and denotes the antinoise generated by as shown in Figure 2. The adaptive filter and linearly operate on the periodic signal to produce which lies in the same subspace as that of . Since is uncorrelated with , convergence in this case essentially means that is minimised in the mean square sense. Thus after convergence, from (32) The update equation for the weight vector is given by where is the signal vector obtained after filtering with . After convergence is attained in , the adaptation is stopped and is fixed to the steadystate value. The system continues to generate antinoise which cancels . The error microphone's output is . The next step is to switch on whose reference input is the separated random part and error input is which it tries to minimize by generating the antinoise signal .
The update equation for the weight vector is given by where is the signal obtained after filtering with . After convergence the error is and should not contain any component of the periodic noise as long as the acoustic system is stationary. Since we have access to only one error sequence, the procedure for finding the optimum and is sequential. However, this sequential adaptive tuning procedure can be repeated when the error sequence energy goes above a prespecified threshold, since the changes in the acoustic system will result in an increase in the energy of the error sequence.
It should be noted that the adaptation methods presented here for the implementation of and are based on simple FXLMS for the sake of simplicity and quick experimentation. For the tonal active noise control part , the analysis and methods described in [9â€“12] can be incorporated. The proposed method can also be extended to accommodate multiple reference microphones in which case the reduced reference signal set procedure developed in [19] can lessen the controller dimension and computational complexity.
3. Calculation of
The goal of this section is to sketch the method for calculating . is obtained, assuming that the adaptive ANC system has converged and that the adaptive ANC system is linear timeinvariant (stationary), as shown in Figure 5(a). First the for wide sense stationary random noise is derived and then expanded to accommodate sinusoidal mixtures in a wide sense stationary random noise. This derivation is important since it gives an insight into how well the proposed method performs when compared to conventional ANC structure with a single cancelling filter. calculation is also very useful during actual realtime implementation of an ANC system. Based on a few basic measurements on an actual realtime ANC setup, the computation of gives an indication of the expected performance. This comes in handy when experimenting with realtime ANC systems like the acoustic laboratory shown in Figure 9. Once the maximum limit is known, the sources of error in an underperforming system can be analyzed.
(a)
(b)
In the three subsections that follow, the expressions for are derived for three cases: when the acoustic noise is (i) widesense stationary, (ii) a sum of sinusoidal signals, and (iii) a mixture of sinusoids and a widesense stationary random noise, that is, combination of (i) and (ii).
3.1. for WideSense Stationary Signals
In linear estimation problems involving stationary signals, Wiener filtering provides the best estimate of a linear timeinvariant (LTI) system based on (minimum mean square error) MMSE criterion [20, 21]. In a linear stationary ANC system as shown in Figure 5(a), the transfer function often has poles outside the unit circle, that is, when is nonminimumphase. In such cases, a stable noncausal Wiener filter is the best estimate of . An expression for MMSE using noncausal Wiener filter is given in [4, 19]. However if is not perfectly predictable, a noncausal Wiener filter is not realizable. The realizeable MMSE estimate is then the causal and stable Wiener filter that approximates .
In order to derive a measurable expression for , let be the causal stable IIR Wiener filter that best approximates : with the superscript and subscript denoting the optimality and causality, respectively, of the solution.
The expression for is obtained as follows:
Let be the crosspower spectral density of two jointly stationary processes and . The expression for is where denotes the crosscovariance sequence. denotes the power spectral density of a stationary process and can be expressed using (37) after appropriate substitution.
Assuming the signals in the system have real positive (power spectrum density) PSD satisfying PaleyWiener condition [21], the stable causal Wiener solution is given in [20, 22]. Referring to Figure 5(b) the causal Wiener solution is where denotes the causal part of the quantity inside. The subscript c stands for causality. is the minimumphase and causal spectral factor of , that is, The next few steps are devoted to the computation of and in terms of , and the input PSD . The transfer function that maps to is and the transfer function that maps to is . Thus using the theory of linear filtering of random processes [21, 23] is related to as follows: It should be noted that
where, the spectral factorization of is is minimumphase and causal. can be decomposed into and as follows: is a polynomial that contains only the zeros outside the unit circle and is responsible for the noncausality of . Hence it follows that is a minimum phase, stable and rational function. Thus (41) becomes
is a polynomial of positive powers of and hence is noncausal. The causal can be expressed as where is the order of .
It should be noted that any stable transfer function can also be written as a product of a stable minimumphase and a stable allpass system. In this case
By definition, is a minimumphase system and will be denoted by .
is an allpass transfer function and will be denoted as . It should be noted that the inverse of is . The above calculations lead to the familiar minimumphase, allpass decomposition [22]
Equation (44) can now be written as The minimumphase and causal is computed as Substituting (49) and (40) in (38) which can be written as and the corresponding output error of the ANC system is computed by where denotes the convolution operation.
When using numerical software, it should be noted that even though the form of in (51) is IIR, it should be considered as a very long FIR filter. The calculation of (51) can be done entirely in the frequency domain by using numerical software like MATLAB. The minimumphase, maximumphase decomposition can be done employing cepstral techniques [22]. The causal part can be found by taking the (Inverse Fast Fourier Transform) IFFT of the quantity inside the bracket and discarding the negative index terms. The maximum possible NAL, that is, the upper bound, can then be computed using for where is the frame length and the frame position and denotes the 2norm operation.
It should be noted that given in (53) can be computed as long as the spectral factorization of (42) is valid.
3.2. for a Mixture of Sinusoids with No Random Noise
When contains sinusoidal components with no random noise, the power spectral density of becomes singular. Hence the technique developed in Section 3.1 cannot be applied. The way out of this difficulty is by recognising that for a pure mixture of sinusoids or a purely periodic signal without any background random noise, the effect of nonminimum phase secondary path can be perfectly equalised. Hence perfect cancellation can be acheived. The nonminimum phase nature essentially translates to a noncausal equaliser which can be realized for purely periodic signals because the future values of the signal are known by indexing the past, that is, by knowing one past period of the signal. For example placing a (noncausal) equalizer as shown in Figure 6 can perfectly equalize a nonminimum phase secondary path. As a result, the effective secondary path becomes the minimumphase , which can be exactly equalized since it is perfectly invertible. Factorizing as , the equalizer is obtained as where is the order of (after sufficient truncation) when expressed in FIR form. Thus . The equalization operations are summarized as where is the period and . Here denotes the integer floor. For cases for which , Thus in theory, a pure mixture of sinusoids can be perfectly cancelled. Thus arbitrarily high noise attenuation level can be reached, that is, .
3.3. for Mixture of Sinusoids in WideSense Stationary Random Noise
When contains sinusoidal components in addition to the widesense stationary random noise, the upper bound becomes a function of and the associated with the random noise part. The derivation below is fundamentally based on the fact that every stationary sequence can be expressed as a sum of a perfectly predictable sequence and a stochastic sequence which are orthogonal to each other (Woldâ€™s Decomposition) [24]. In this case the sinusoidal components constitute the perfectly predictable part and the background broadband noise component constitutes the stochastic part.
The noise attenuation level by definition is given by where is the power of , is the power of , is the power of the tonal component (after convergence) of , and is the power of the random component (after convergence) of .
Due to linearity, . Since the deterministic component and the random component are independent, The performance gain due to the proposed separationbased parallel feedforward method can be understood from (60). By using a dedicated cancelling filter, the periodic part can be almost perfectly canceled, that is, . Perfect cancellation of the periodic part is not possible with conventional single filter ANC. Hence the denominator of (60) is always higher for the conventional method. On the assumption that the periodic part can be perfectly cancelled, (60) becomes in (61) defines an upper bound for the ANC system since it is the maximum NAL possible under the assumptions of perfect separation and perfect cancellation of the periodic part . Thus for a mixture of periodic tones in random noise, the maximum achievable NAL, that is, the is given by . Denoting the NAL of an ANC system to a random noise input by , the expression in (61) can be written as The above expression generalizes the form of maximum achievable NAL derived for stationary random noise in Section 3.1 to stationary random noise along with sinusoids. It should be noted that of (62) is calculated from (51)â€“(53) after setting , that is, . To see the effect of the primary path it should be noticed that where â€™s are the frequencies of the sinusoids in .
In case of a single tone (), and for The expression becomes For the parallel ANC system, (66) can be written as
That is, the improvement in NAL of a mixture of periodic and random noise signals is greater than that of the performance of a purely random noise by a factor of . If we assume , then . It seems that this way of characterizing maximum achievable performances of ANC system is new. Furthermore, is low when there are deep notches in at the frequency , in which case the relative position of reference microphone should be changed.
As a summary, the procedure described in the paper is as follows.
Step 1. Given a feedforward ANC set up as shown in Figure 9, record and estimate and . The same ANC setup can be used for both measurements. The primary path can be measured switching off the antinoise generation (i.e., ) and by recording the measurements from the reference microphone and from error microphone, and estimating the transfer function that maps to . For measurement of , methods given in [4] can be employed.
Step 2. If does not have tones use (51) through (53) to calculate . The conventional singlefilter approach is sufficient for this case. Additionally the optimum Wiener filter can also be calculated from (51) and can be used to initialize the adaptive filter .
Step 3. If has tones along with random noise, use (67) to calculate the . This step needs calculation of the power spectral density of which can be generated from using the adaptive noise cancellationbased separation that was developed in Sections 2.3.1 and 2.3.2.
Step 4. Implement the PFANC scheme. requires the calculation of the period from the autocorrelation sequence of and the calculation of from (55). Additionally the optimum Wiener filter can also be calculated from (51) and can be used to initialize the adaptive filter .
4. Experimental Results
Simulation of the proposed method, using FXLMS for adaptation, was implemented by synthesizing for different SNR values and different linear combinations of multiple sinusoids. Single and multitone signals of various frequencies were synthesized and of varying power was added. The and used in these experiments were obtained from [4] and characterize real life acoustic transfer functions (Figure 7).
Each experiment was done on a 100000 sample long generated at a sampling rate of 16â€‰kHz. The first 10000 samples were used for the frequency estimation part, that is, determination of , , , . The next 20000 samples were used to train the adaptive noise canceller that is, to separate and as discussed in Section 2.3.2. The rest of the samples were used to generate antinoises via and . The performance of the proposed method is compared with the conventional FXLMSbased ANC [4]. The results are shown in Table 1. The last column shows the maximum achievable noise attenuation level using (62). was found to be 19.58â€‰dB for the and used. It can be seen that the NAL obtained by the proposed method is closer to than that obtained by the conventional singlefilter method. For low SNR values (like 10â€‰dB) the performance improvement is less noticeable when compared to higher SNR values. For higher SNR values excellent performance of the proposed method can be observed and thus the proposed method targets high SNR applications. It should also be noticed that perfect separation and cancellation of the periodic part may not be achieved in practice resulting in maximum NAL being less than that given by (62).

Experiments 21 through 24 were carried out using a signal derived from functional magnetic resonance imaging (fMRI) acoustic noise [25, 26]. fMRI acoustic noise has SNRâ€™s 20â€‰dB and hence amenable to the proposed approach. Seven significant tones from a recorded fMRI acoustic noise were estimated and the sinusoidal component was created from the seven tones. White noise was then added with different SNR values. Figures 8(a) and 8(b) present a comparison of noise cancellation performance between the conventional singlefilter method and the proposed PFANC for an 8 tones case and a 10 tones case, respectively. As it can be seen, PFANC considerably reduces the background noise in addition to sinusoidal noise.
(a)
(b)
5. Conclusion
In this paper the focus has been on active noise control of acoustic signals that have a periodic/harmonic structure along with a random noise. The key idea stems from the observation that the periodic component and the random component may need different noise cancelling filters for the best performance. This leads to the idea of separation of the periodic and random components and the proposed parallel form of ANC, in which the reference input is separated into periodic and random components and separate antinoise signals are parallely generated. A novel adaptive filteringbased implementation of signal separation and antinoise generation is developed and simulated. A simple but an insightful formula is derived to quantify performance. The method is shown to work very well with single and multitone signals in white noise.
Some important and recurring signals, blocks, and notation that are used in the paper are listed in the Notations and Definitions section. Other definitions are made as they are introduced.
Notations and Definitions
It gives the signal energy. For a signal , the energy is given by , where is the expectation operator.  
f:  The subscript â€śâ€ť for any quantity denotes the component of , corresponding to the subscripted variable. 
:  The version of acoustic noise recorded by a reference microphone placed somewhere in the noisy area. This measurement is known as the reference signal in ANC literature. 
:  The periodic component of noise which can be written as a sum of sinusoids, that is, . 
:  The random noise component. It is assumed to be zeromean and stationary, where . 
:  The period of the periodic component. can be measured from the autocorrelation sequence of . 
:  The (signaltonoise ratio) SNR of periodic component to the random component, that is, . 
:  The individual SNR of sinusoidal components, that is, . It follows that . 
:  The measurement from error microphone. The cancellation is achieved around the error microphone. The error signal corresponding to a particular component will have a subscript to denote that. For example, the component of error signal corresponding to will be denoted as and its power will be . 
:  The stable linear transfer function representing the primary path , which is the path through the acoustic cavity between reference and error microphones. The corresponding impulse response is . 
:  The stable linear transfer function representing the secondary path . models system functions of microphones, digital and analog circuits, the cancelling loudspeaker, and the acoustic gaps between the loudspeaker and error microphone. is stable and assumed to have no zeros on the unit circle. 
:  An estimate of , which is placed in the reference signal path. In all the ensuing discussions, is considered as an accurate model, that is, =. 
:  The noise signal after passing through representing the noise signal to be canceled, that is, the desired signal . It should be noted that . 
:  The main adaptive noise cancelling (control) filter which generates the appropriate antinoise signal using the error and reference input signals. 
:  The control signal generated by . gets shaped into when it reaches the error microphone due to the effect of the secondary path . 
:  The maximum achieveable Noise Attenuation Level possible for a given ANC. algorithm/strategy. The for an input reference signal will be denoted as 
:  Noise Attenuation Level (NAL) of an ANC system which has the reference input , where . Depending on the context, can also indicate the maximum achievable NAL for the given reference signal . 
Acknowledgments
This study was supported by the VA IDIQ contract number VA549P0027 awarded and administered by the Dallas, TX VA Medical Center. The content of this paper does not necessarily reflect the position or the policy of the U.S. government, and no official endorsement should be inferred.
References
 M. Pawelczyk, â€śAdaptive noise control algorithms for active headrest system,â€ť Control Engineering Practice, vol. 12, no. 9, pp. 1101â€“1112, 2004. View at: Publisher Site  Google Scholar
 M. O. Tokhi, M. A. Hossain, and M. H. Shahed, Parallel Computing for RealTime Signal Processing and Control, Springer, London, UK, 2003.
 M. A. Hossain and M. O. Tokhi, â€śRealtime design constraints in implementing active vibration control algorithms,â€ť International Journal of Automation and Computing, vol. 3, no. 3, pp. 252â€“262, 2006. View at: Google Scholar
 S. M. Kuo and D. R. Morgan, Active Noise Control Systems: Algorithms and DSP Implementation, John Willey & Sons, New York, NY, USA, 1996.
 S. Johansson, S. Nordebo, and I. Claesson, â€śConvergence analysis of a twinreference complex leastmeansquares algorithm,â€ť IEEE Transactions on Speech and Audio Processing, vol. 10, no. 4, pp. 213â€“221, 2002. View at: Publisher Site  Google Scholar
 M. Pawelczyk, â€śFeedforward algorithms with simplified plant model for active noise control,â€ť Journal of Sound and Vibration, vol. 255, no. 1, pp. 77â€“95, 2002. View at: Publisher Site  Google Scholar
 M. O. Tokhi and S. M. Veres, Active Sound and Vibration Control: Theory and Applications, Institution of Electrical Engineers, London, UK, 2002.
 G. Kannan, A. A. Milani, and I. M. S. Panahi, â€śActive noise control of noisy periodic signals using signal separation,â€ť in Proceedings of IEEE International Conference on Acoustic, Speech, and Signal Processing (ICASSP '08), pp. 1617â€“1620, Las Vegas, Nev, USA, MarchApril 2008. View at: Publisher Site  Google Scholar
 S. M. Kuo and A. B. Puvvala, â€śEffects of frequency separation in periodic active noise control systems,â€ť IEEE Transactions on Audio, Speech and Language Processing, vol. 14, no. 5, pp. 1857â€“1866, 2006. View at: Publisher Site  Google Scholar
 Y. Hinamoto and H. Sakai, â€śAnalysis of the filteredX LMS algorithm and a related new algorithm for active control of multitonal noise,â€ť IEEE Transactions on Audio, Speech and Language Processing, vol. 14, no. 1, pp. 123â€“130, 2006. View at: Publisher Site  Google Scholar
 L. Vicente and E. Masgrau, â€śNovel FxLMS convergence condition with deterministic reference,â€ť IEEE Transactions on Signal Processing, vol. 54, no. 10, pp. 3768â€“3774, 2006. View at: Publisher Site  Google Scholar
 Y. Zhang, P. G. Mehta, R. R. Bitmead, and C. R. Johnson Jr. , â€śDirect adaptive control for tonal disturbance rejection,â€ť in Proceedings of the IEEE American Control Conference, vol. 3, pp. 1480â€“1482, Philadelphia, Pa, USA, June 1998. View at: Publisher Site  Google Scholar
 N. J. Bershad and J. C. M. Bermudez, â€śSinusoidal interference rejection analysis of an LMS adaptive feedforward controller with a noisy periodic reference,â€ť IEEE Transactions on Signal Processing, vol. 46, no. 5, pp. 1298â€“1313, 1998. View at: Publisher Site  Google Scholar
 P. Stoica and R. Moses, Introduction to Spectral Analysis, PrenticeHall, Upper Saddle River, NJ, USA, 1997.
 M. Wax and T. Kailath, â€śDetection of signals by information theoretic criteria,â€ť IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 33, no. 2, pp. 387â€“392, 1985. View at: Publisher Site  Google Scholar
 P. Stoica and Y. SelĂ©n, â€śModelorder selection: a review of information criterion rules,â€ť IEEE Signal Processing Magazine, vol. 21, no. 4, pp. 36â€“47, 2004. View at: Publisher Site  Google Scholar
 J. Rissanen, â€śModeling by shortest data description,â€ť Automatica, vol. 14, no. 5, pp. 465â€“471, 1978. View at: Publisher Site  Google Scholar
 P. Stoica and T. SĂ¶derstrĂ¶m, â€śStatistical analysis of MUSIC and subspace rotation estimates of sinusoidal frequencies,â€ť IEEE Transactions on Signal Processing, vol. 39, no. 8, pp. 1836â€“1847, 1991. View at: Publisher Site  Google Scholar
 W. Dehandschutter, J. van Herbruggen, J. Swevers, and P. Sas, â€śRealtime enhancement of reference signals for feedforward control of random noise due to multiple uncorrelated sources,â€ť IEEE Transactions on Signal Processing, vol. 46, no. 1, pp. 59â€“69, 1998. View at: Publisher Site  Google Scholar
 M. H. Hayes III, Statistical Digital Signal Processing and Modeling, John Willey & Sons, New York, NY, USA, 1996.
 A. Papoulis and S. U. Pillai, Probability, Random Variables, and Stochastic Processes, McGrawHill, New York, NY, USA, 4th edition, 2001.
 J. G. Proakis and D. G. Manolakis, Digital Signal Processing, PrenticeHall, Upper Saddle River, NJ, USA, 4th edition, 2006.
 S. Elliott, Signal Processing for Active Control, Academic Press, New York, NY, USA, 2001.
 A. Papoulis, â€śLevinson's algorithm, Wold's decomposition, and spectral estimation,â€ť SIAM Review, vol. 27, no. 3, pp. 405â€“441, 1985. View at: Publisher Site  Google Scholar
 A. A. Milani, G. Kannan, I. M. S. Panahi, R. Briggs, and K. Gopinath, â€śWeight stacking analysis of delayless subband adaptive algorithms for fMRI active noise cancellation,â€ť in Proceedings of the IEEE Dallas Engineering in Medicine and Biology Workshop (DEMBS '07), pp. 134â€“137, Dallas, Tex, USA, November 2007. View at: Publisher Site  Google Scholar
 V. R. Ramachandran, G. Kannan, A. A. Milani, and I. M. S. Panahi, â€śSpeech enhancement in functional MRI environment using adaptive subband algorithms,â€ť in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07), vol. 2, pp. 341â€“344, Honolulu, Hawaii, USA, April 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2009 Govind Kannan 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.