The Journal of the Acoustical Society of Korea. 31 May 2016. 216-222
https://doi.org/10.7776/ASK.2016.35.3.216

ABSTRACT


MAIN

  • I. Introduction

  • II. Multichannel delay- compensated FXLMS algorithm

  • III. Real-time ANC implementation

  • IV. Experimental Results

  • V. Conclusions

I. 서  론

ANC (Active Noise Control) cancels the primary noise by generating and combining an anti-noise (with equal amplitude but opposite phase) at the location of the error microphone, and it efficiently attenuates low frequency noise with benefits in size and cost.[1,2] Over the past few decades a great progress has been made in ANC, of which main purpose is generating a quiet zone, where the noise is cancelled. Generally, the quiet zone is generated around the position of error microphone, however, the size of this quiet zone is usually limited especially in a single channel (which means one reference microphone, one cancelling speaker and one error microphone) ANC system.[3] Thus generating a quiet zone with sufficient size and good performance is of considerable significance.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F1.jpg

Fig. 1. Block diagram of multichannel delay-compensated FXLMS algorithm.

This paper applies a multichannel ANC structure to create a quiet zone with satisfactory size. The multichannel FXLMS algorithm[2] has been widely used for ANC application.[4-6] To improve the performance of quiet zone, a multichannel delay-compensated FXLMS (Filtered-X Least Mean Square) algorithm[7,8] with faster convergence rate and better performance is selected for our multichannel ANC. The proposed multi-channel ANC employes one reference microphone, two cancelling speakers and four error microphones, details on it are discussed in later section. The real-time system is designed and implemented in TMS320C6713 floating-point DSP board, other devices, including microphones and speakers, are connected to DSP board by the self-designed circuit.

The system is tested to cancel 100 to 500 Hz tonal frequency noises generated by a function generator and primary loudspeaker. Experiment results prove that our system generates a quiet zone with sufficient size and maximum 24 dB noise attenuation.

The organization of the paper is as follows: Section II presents the multichannel delay-compensated FXLMS algorithm. Section III describes the real-time system setup and implementation. Section IV illustrates the results the quiet zone. The conclusions are summarized in Section V.

II. Multichannel delay- compensated FXLMS algorithm

FXLMS algorithm has been widely used in practical ANC due to its robustness and simplicity, numerous algorithms are researched to improve its performance, but most of them achieve better performance with a significant computational complexity increase, which cannot be applied for real-time implementation. In the proposed system, a multichannel delay-compensated FXLMS algorithm is applied. It achieves faster convergence rate and better performance with only double computation load, and guarantees the possibility of real-time implementation. The multichannel delay- compensated FXLMS is illustrated in Fig. 1. It employs http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC193E.gif reference sensors to form the reference signal vector http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC197E.gif, generates http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC198E.gif cancelling signals http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC19CE.gif to drive the corresponding secondary sources, and distributes http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC19FE.gif error sensors over desired locations to measure the residual noise components http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1A2D.gif.

The wide arrows represent an array of signals that are symbolically expressed as vectors. The matrix http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1BB5.gif represents http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1BF5.gif primary path transfer functions, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1C44.gif, from the reference signal http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1C74.gif to each error sensor output http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1CC3.gif. The matrix http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1D02.gif represents http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1D42.gif secondary-path transfer functions, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1D72.gif, from http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1D92.gif secondary sources to http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1DB2.gif error sensors. There are http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1DF2.gif possible feedforward channels, each demanding a separate adaptive filter, and these adaptive filters are represented by the matrix http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1E31.gif. The delay-compensated FXLMS algorithm uses the estimated plant models http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1EAF.gif to subtract the contribution of the secondary path from the error signals, so that the estimated primary field signals http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1F2D.gif are obtained, indeed, the adaptive filters http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1FAB.gif then try to predict the estimated signals http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC1FEB.gif instead of the original http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC202A.gif signals.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F2.jpg

Fig. 2. SNR curves of FXLMS and delay-compensated FXLMS.

Using the following additional notations, we then obtain the multichannel delay-compensated FXLMS algorithm described by Eqs.(1) ~ (5), where http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC204A.gif is the length of finite impulse response (FIR) adaptive filters in http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC208A.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC20AA.gif is the length of FIR adaptive filters in http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC20DA.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2129.gif is the value at time http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2149.gif of the http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC216A.gifth cancelling noise signal, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC21D8.gif is the value at time http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC21F8.gif of the http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2218.gifth coefficient in http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2239.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2259.gif is the value of the http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2279.gifth coefficient of http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC22B9.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC22F8.gif is the value at time http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2318.gif of the filtered reference signal, http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2348.gif is the estimate of the value at time http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2369.gif of the primary sound field http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC23E7.gif at the http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2407.gifth error sensor, and http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2437.gif is the alternative error signal.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC24D4.gif

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC25A0.gif (1)

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC25FF.gif

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC26BB.gif (2)

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC272A.gif

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC27E6.gif (3)

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2893.gif

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC296F.gif (4)

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2A3B.gif (5)

Table 1 presents the number of multiples required for an iteration of both algorithms. Comparing with multichannel FXLMS algorithm, the delay-compensated FXLMS algorithm has almost twice computational load, but much faster convergence speed and better performance (see Fig. 2). 10 fan noises, 10 factory noises and 10 car noises are tested to the performances of two algorithms. The SNR is defined as the ratio between the power of residual noise and power of primary noise. The average SNR curves of 30 simulations are presented in Fig. 2.

Table 1. Computational complexity (number of multiplies) of multichannel algorithms.

FXLMS

Delay-compensated FXLMS

Analytical computational load

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PICCA5D.gif

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PICCB39.gif

Computational load of http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PICCBA7.gif

141

231

Computational load of http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PICCBE6.gif

774

1544

III. Real-time ANC implementation

A single channel real-time ANC system[3] is presented in our previous research. It applies one reference microphone, one cancelling speaker, one error microphone and FXLMS algorithm. Its structure is shown in Fig. 3.

The experiment result in the next section proves a quiet zone is generated around the location of error microphone. But it must be noticed that the area of quiet zone is limited, and it is inefficient and meaningless for practical usage. To increase the size of quiet zone, multichannel ANC system with one reference microphone, two cancelling speakers and four error microphones is designed.

The proposed multichannel system with multichannel delay-compensated FXLMS algorithm is shown in Fig. 4.

The primary noise is generated by function generator and output from a powered studio monitor speaker- YAMAHA HS50, which has a frequency response from 100 to 20 kHz with a maximum power output of 70 W.

The reference microphone and error microphones are unidirectional condenser microphones with 2.2 kΩ impedance, -46±4 dB sensitivity and 70 to 20 kHz frequency response range. The cancelling loudspeakers have the same specifications as the primary loudspeaker.

The system is implemented in TMS320C6713 DSP board, the board employs 12-bit 8-channel A/D convertor AD7888 and 12-bit 4-channel D/A convertor AD7564, both analog input voltage of A/D convertor and analog output voltage of D/A convertor are 0 to 2.5 V. That is the reason we designed +1.25 V bias voltage before the ADC signal input and -1.25 V bias voltage after the DAC signal output. The anti-aliasing filter and reconstruction filter are 4-order low-pass Butterworth filters; their passband frequency is 750 Hz and stopband frequency is set to 1000 Hz with -9 dB stopband attenuation. For the operational amplifier, a balance input microphone amplifier circuit with adjustable gain is designed using high fidelity audio operational amplifier IC LME 49740.

To better understand our system, the equipment setup of multichannel ANC is illustrated in Fig. 5. To simulate the 3 dimensional environment, our system is implemented in a 15×10×3 m room, the heights of microphones and speakers are 1.5 m with their stands.

IV. Experimental Results

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F3.jpg

Fig. 3. The structure of single channel real-time ANC system.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F4.jpg

Fig. 4. The structure of proposed multichannel real-time ANC system.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F5.jpg

Fig. 5. System setup of multichannel real-time ANC.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F6.jpg

Fig. 6. The quiet zone generated by single channel ANC system.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F7.jpg

Fig. 7. The quiet zone generated by multichannel ANC system with FXLMS.

http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/Figure_ASK_35_03_07_F8.jpg

Fig. 8. The quiet zone generated by proposed multi-channel ANC system.

The real-time ANC system is tested to cancel tonal noises, with frequency range from 100 to 500 Hz, generated by function generator. The system uses 8000 Hz sampling rate, the length of adaptive filter http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2AB9.gif and http://static.apub.kr/journalsite/sites/ask/2016-035-03/N0660350307/images/PIC2AF8.gif are selected as 64 and 128 respectively. A SPL meter is used to measure the quiet zone in terms of noise attenuation (dB). The results are shown in Figs. 6 ~ 8, which present the quiet zone in a 60×60 cm area generated by single channel ANC, multichannel ANC with FXLMS and proposed multichannel ANC with delay- compensated FXLMS.

A coordinate system, as shown in Fig. 5, is to better present the location of equipment. In single channel ANC implementation. The primary noise speaker is placed in (200, 0) (in centimeter), the reference microphone is placed in (150, 0), the cancelling speaker is placed in (40, 40) and error microphone is set at (0, 0). From Fig. 6, a quiet zone with maximum 20 dB attenuation is generated around the location of error microphone. But the area of quiet zone is limited, the quiet zone with 12dB noise attenuation is smaller than a 10×10 cm area, that is inefficient and meaningless for practical usage.

For multichannel ANC systems, the primary noise speaker is placed in (350, 0), the reference microphone is placed in (300, 0), the cancelling speakers are symmetrically placed in (0, -120) and (0, 120), the error microphone are distributed in (20, 20), (-20, 20), (-20, -20) and (20, -20). Comparing the results in Figs. 6 and 7, a multichannel structure successfully enlarges the quiet zone to a satisfactory size. From Figs. 7 and 8, the noise attenuation is evidently improved with multichannel delay-compensated FXLMS algorithm. From Fig. 8, the quiet zone with 12 dB noise attenuation is up to a 40×40 cm area. This proves the quiet zone performance of proposed ANC system is successfully enhanced in both in terms of zone size and noise attenuation.

V. Conclusions

Generating a quiet zone, in which the primary noise is cancelled, is the main objective of ANC. For practical application, this quiet zone must be large enough to satisfy the needs. To enhance the noise cancellation in quiet zone, including enlargement of zone size and improvement of noise cancellation gain, the proposed multichannel real-time ANC system is designed. The proposed system employs multichannel structure, which includes one reference micro-phone, two cancelling speakers and four error microphones, to enlarge the quiet zone; a multichannel delay-compensated FXLMS is then applied to improve the noise attenuation.

The quiet zones generated by ANC system are measured by SPL meter in terms of noise attenuation in dB. The results measured in a 60×60 cm area prove that, comparing with the single channel ANC system and multichannel ANC with FXLMS algorithm, the proposed ANC system successfully enhances the quiet zone performance in terms of zone size and noise attenuation.

References

1
1.S. J. Elliot, Signal Processing for Active Control (Academic Press, London, UK, 2001),  pp. 6-13.
2
2.S. M. Kuo and D. R. Morgan, Active Noise Control Systems: Algorithms and DSP Implementations (Wiley, New York, NY, 1996), pp. 5-16.
3
3.X. Mu, J. Ko, and J. Rheem, “Modified FxLMS algorithm for active noise control and its real-time implementation,” IEE. Kr. 50, 172-176 (2013).
4
4.A. Esperance, M. Bouchard, B. Paillard, C. Guigou, and A. Boudreau, “Active noise control in large circular duct using an error sensors plane,” Applied Acoust. 57, 357-374, (1999).
5
5.K. Beemanpally, K. Pottim, and S. Kuo, “Multi-channel hybrid active noise control system for infant incubators,” in Proc. IEEE Int. Conf. 69-76 (2010).
6
6.T. Jung, J. Kim, K. Kim, and S. Nam, “Active noise reduction system using multi-channel ANC,” 11th Int. Conf. on Control, Automation and Systems, (2011).
7
7.E. Bjarnason, “Active noise cancellation using a modified form of the filtered-x LMS algorithm,” in Proc. 6th Eur. Signal Processing Conf. 2, 1053-1056 (1992).
8
8.M. Bouchard and S. Norcross, “Computational load reduction of fast convergence algorithms for multichannel active noise control,” Signal Proce. 83, 121-134 (2003).
페이지 상단으로 이동하기