background image

GOLDBERG AND MÄKIVIRTA

 

AUTOMATED IN-SITU EQUALISATION

 

 

AES 23RD CONFERENCE, May 23-25, 2003 

rithm discussed in Section 3 and statistical data for 
each measurement before and after equalisation was 
recorded. The statistical data is analysed to study how 
the objective quality of the system magnitude re-
sponse has been improved by using the proposed algo-
rithm for setting the room response controls. 

Table 10. Acoustic measurement system parameters. 

Parameter 

Equipment / Setting 

Measurement System 

WinMLS2000 [39] 

Microphone 

Neutrik 3382 [40] 

Sample rate, 

f

s

 48 

kHz 

MLS sequence order 

14 (16) 

Averages 1 
Impulse response length 

0.341 s (1.36 s) 

Time window 

Half-cosine 

FFT size 

16384 (65536) 

Frequency resolution 

2.93 Hz (0.733 Hz) 

 

4.1.  Statistical Data Analysis 

A further statistical analysis was conducted on all of 
the loudspeakers in the study. The bandwidths of the 
frequency bands used are shown in Table 11. The 
bandwidths ‘

LF

’, ‘

MF

’ and ‘

HF

’ are later referred to 

collectively as the ‘

subbands’

 and correspond roughly 

to the bandwidths for each driver in the three-way sys-
tems. 

Table 11. Frequency band definitions the statistical 
data analysis; 

f

LF

 is the frequency of the lower –3 dB 

limit of the frequency range. 

 

Frequency Range Limit 

Bandwidth Name 

Low 

High 

Broadband 

f

LF

 15 

kHz 

LF 

f

LF

 400 

Hz 

MF 

400 Hz 

3.5 kHz 

HF  

3.5 kHz 

15 kHz 

 
For each loudspeaker, the broadband (Table 11) mag-
nitude response data median value is standardised to 
0 dB. 
The statistical descriptors recorded before and after 
equalisation for each loudspeaker and in each fre-
quency band defined in Table 11 are the minimum, 
maximum and range of the magnitude dB values. Also 
for the magnitude pressure values in each bandwidth 
(Table 11), the median, 5% & 95% percentiles and 
quartiles are recorded. In addition, the root-mean-
square (RMS) deviation of the pressure from the me-
dian in each bandwidth is calculated: the value is ex-
pressed in dB. 

These statistical descriptors are compared for each 
subband to study the in-band flatness improvement 
due to equalisation. The median values for each sub-
band are compared to study the broadband tonal bal-
ance improvement. This is indicated by a reduction of 
the median value differences. 
 

4.2.  Example of Statistical Data Analysis 

Figure 7 in Appendix C shows a case example where 
room response control settings are calculated accord-
ing to the optimisation algorithm. The equalisation 
target is a flat magnitude response (straight line at 
0 dB level). The in-situ frequency response of the 
loudspeaker was recorded before equalisation, i.e. 
when all the room response controls were set to their 
default position, which has no effect to the response. 
The appropriate room response control settings were 
calculated using the optimisation algorithm, applied to 
the loudspeaker and the corrected in-situ frequency 
response plotted. The loudspeaker’s passband (trian-
gles) and the frequency band of equalisation (crosses) 
are indicated on the graphical output. The proposed 
room response control settings are shown and the ef-
fect of these settings is visualised in the response plot. 
The treble tilt, midrange level and bass tilt controls 
have been set. The equalisation corrects the low fre-
quency alignment and improves the linearity across 
the whole passband. 
Figure 8 in Appendix C shows a statistical analysis of 
the same loudspeaker presented in graphical form. 
The upper three plots were calculated before equalisa-
tion and the lower three plots after equalisation. The 
plots display the values of percentiles in the magni-
tude value distribution (box plot), the histogram of 
values and the fit of the magnitude values to normal 
distribution before and after equalisation. These plots 
clearly show that the distribution in magnitude data 
has been reduced. This is illustrated by the reduced 
range in the box plot and the value histogram, as well 
as a better fit to a normal distribution in the normal 
probability plot. 
 

4.3. Results 

A total of 63 loudspeakers were measured before and 
after equalisation. Of these, 12 were small two-way 
systems, 22 were two-way systems, 30 were three-
way systems and three were large systems. 
Depending on the product type, not all of the room 
response controls are available (Tables 1–4). Table 12 
shows the number times the controls were used when 
available on the loudspeaker. The midrange level con-
trol is used most frequently and the bass roll-off the 
least. 

Содержание Frequency Response Optimisatio

Страница 1: ...ce of calibrating active loud speakers Even with experienced system calibrators significant variance between calibrations can be seen Furthermore with a number of different people cali brating loudspeaker systems additional variance in results will occur For these reasons an automated calibration method was developed to ensure consis tency of calibrations Presented first in this paper is the discr...

Страница 2: ...al to significantly improve perceived sound quality The practical challenge is the selection of the best settings for the low order in situ equaliser Despite advances in psychoacoustics it is difficult to quantify what the listener actually perceives the sound quality to be or to optimise equalisation based on that evaluation 13 15 Because of this in situ equalisa tion typically attempts to obtain...

Страница 3: ...re determined by the crossover filters The bass tilt control compensates for a bass boost seen when the loudspeaker is loaded by large nearby boundaries 33 36 This typically happens when a loudspeaker is placed next to or mounted into an acoustically hard wall This filter is a first order shelv ing filter The bass roll off control compensates for a bass boost often seen at the very lowest frequenc...

Страница 4: ...her than using a least squares type objective function is that the bass roll off tends to assume maximum attenuation to minimise the RMS deviation This type of objective function does not yield the best setting as subjectively a loss of bass extension is perceived This stage of the optimiser algorithm takes six filtering steps three for small two way models 3 2 2 Midrange Level to Treble Level Rat...

Страница 5: ...tics of experienced calibration engineers The resulting num ber of filtering steps has been dramatically reduced for the larger systems Table 9 and even the relatively simple two way systems show a substantial improve ment when compared to the number of filtering steps needed by direct search method as summarised in Table 5 There are two main reasons for the improve ment the constraint of not allo...

Страница 6: ... there is also a 3 dB roll off with 50 Hz being down by 1 dB and 40 Hz by 2 dB Tolerance lines are set to 3 dB with additional leeway at low and high fre quencies 1 An example of the room equaliser settings output for the large system optimised in Figure 1 is shown in Figure 2 The optimised result is displayed in green and dark grey boxes The green boxes are room re sponse controls that should be ...

Страница 7: ...d tonal bal ance improvement This is indicated by a reduction of the median value differences 4 2 Example of Statistical Data Analysis Figure 7 in Appendix C shows a case example where room response control settings are calculated accord ing to the optimisation algorithm The equalisation target is a flat magnitude response straight line at 0 dB level The in situ frequency response of the loudspeak...

Страница 8: ...subbands show no changes or a slight increase of the RMS deviation Three way systems show a clear reduction in most cases of both the quartile difference Figure 13 and RMS deviation Figure 14 for the broadband and LF subband Slight and equal numbers of increases and reductions are seen for MF and HF subbands A similar trend is seen for the three large systems in cluded in this study Figure 15 16 M...

Страница 9: ...n the median level for the LF subband A similar outcome is noted sepa rately for each loudspeaker type However only in the three way systems is an improvement seen also in the MF and HF subband variance 25 to 75 Percentile Difference Change due to Equalisation All models 3 2 1 0 1 Broadband LF MF HF Level dB RMS Deviation Change due to Equalisation All models 5 4 3 2 1 0 1 Broadband LF MF HF Level...

Страница 10: ...d room response control settings The settings achieve im proved equalisation in the form of a smaller RMS de viation from the target response The improvement is not limited by the optimisation method but by the room response controls which are not intended to cor rect for narrow band deviations in the frequency re sponse Examples of these are response variations re sulting from acoustic issues suc...

Страница 11: ...Thesis at the Helsinki University of Technol ogy 41 8 REFERENCES 1 Genelec Oy http www genelec com 2003 Feb 2 Walker R Equalisation of Room Acoustics and Adaptive Systems in the Equalisation of Small Rooms Acoustics Proc 15th Int Conf paper 15 005 1998 Oct 3 Cox T J and D Antonio P Determining Op timum Room Dimensions for Critical Listening Envi ronments A New Methodology presented in 110th Conv A...

Страница 12: ...2001 Sep 29 Mäkivirta A Antsalo P Karjalainen M and Välimäki V Low Frequency Modal Equalisation of Loudspeaker Room Responses presented in 111th Conv Audio Eng Soc preprint 5480 2001 Sept 30 Karjalainen M Esquef P A A Antsalo P Mäkivirta A and Välimäki V Frequency Zooming ARMA Modelling of Resonant and Reverberant Sys tems J Audio Eng Soc vol 50 pp 1012 1029 2002 Dec 31 Moore B C J Glasberg B R Pl...

Страница 13: ...User Inputs Model Database Stored Measurement Microphone Compensation CTRL M Measurement Dump Reset Graph and Outputs Get Model Number Apply Mic Compensation Remove DC Window FFT and Smooth Load Impulse Response Set DIPtimisation Range Display Original Freq Response Display Target Response Calculate Target Resp Stored Measurement CLOSE DIPtimiser 1 2 Figure 5 Software flow chart part 1 CLOSE Set F...

Страница 14: ...5 2003 14 Is Large System Is Small System Load Filters Model Filters Preset BRO Find ML TL Ratio Set BL BT wrt ML TL Reset BRO Set TT Display Final Tone Control Settings Display Final Frequency Response Set BT Is 3 way System 1 2 Figure 5 continued Software flow chart part 2 Y N N Y ...

Страница 15: ...GOLDBERG AND MÄKIVIRTA AUTOMATED IN SITU EQUALISATION AES 23RD CONFERENCE May 23 25 2003 15 APPENDIX B SOFTWARE GRAPHICAL USER INTERFACE Figure 6 Software graphical user interface at start up ...

Страница 16: ... AND MÄKIVIRTA AUTOMATED IN SITU EQUALISATION AES 23RD CONFERENCE May 23 25 2003 16 APPENDIX C CASE EXAMPLE STATISTICAL GRAPHS Figure 7 Case example optimisation results Figure 8 Case example statistical output ...

Страница 17: ...5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB Low Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB Low Frequency 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1029A 1029A 1029A 1029A 1029A 1029A...

Страница 18: ...2 0 1 5 1 0 0 5 0 0 0 5 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB High Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB High Frequency 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 102...

Страница 19: ...to Equalisation 7 6 5 4 3 2 1 0 1 2 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB Low Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB Low Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB ...

Страница 20: ...Equalisation 7 6 5 4 3 2 1 0 1 2 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB High Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB High Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A 1029A Level dB H...

Страница 21: ...ence Before Equalisation 0 2 4 6 8 10 12 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB Low Frequency 25 to 75 Percentile Difference Change due to Equalisation 4 5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB Broadband 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1030A 1030A 1031A...

Страница 22: ... 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB High Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB High Frequency 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A ...

Страница 23: ...ange due to Equalisation 7 6 5 4 3 2 1 0 1 2 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB Low Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB Low Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB Low Frequ...

Страница 24: ...e due to Equalisation 7 6 5 4 3 2 1 0 1 2 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB High Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB High Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1030A 1030A 1031A 1031A 1031A 1031A 1031A 1031A 1031A 1032A 1032A Level dB High Frequ...

Страница 25: ...tion 0 2 4 6 8 10 12 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB Low Frequency 25 to 75 Percentile Difference Change due to Equalisation 4 5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB Broadband 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 ...

Страница 26: ... 0 5 0 0 0 5 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB High Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB High Frequency 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 S30D S30D S30D S30D 1037B 1037B 10...

Страница 27: ...tion 7 6 5 4 3 2 1 0 1 2 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB Low Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB Low Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 10...

Страница 28: ...n 7 6 5 4 3 2 1 0 1 2 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB High Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 1038A 1038A 1038A 1039A Level dB High Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 S30D S30D S30D S30D 1037B 1037B 1037B 1037B 1037B 1038A 1038A 103...

Страница 29: ...evel dB Low Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 1036A 1036A 1036A Level dB Low Frequency 25 to 75 Percentile Difference Change due to Equalisation 4 5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 1036A 1036A 1036A Level dB Broadband 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1036A 1036A 1036A Level dB Low Frequency 25 to 75 Percentile Diff...

Страница 30: ...ifference Change due to Equalisation 4 5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 0 5 1036A 1036A 1036A Level dB High Frequency 25 to 75 Percentile Difference Before Equalisation 0 2 4 6 8 10 12 1036A 1036A 1036A Level dB High Frequency 25 to 75 Percentile Difference After Equalisation 0 2 4 6 8 10 12 1036A 1036A 1036A Level dB High Frequency 25 to 75 Percentile Difference Change due to Equalisation 4 ...

Страница 31: ... 1036A 1036A Level dB Broadband RMS Deviation Change due to Equalisation 7 6 5 4 3 2 1 0 1 2 1036A 1036A 1036A Level dB Low Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1036A 1036A 1036A Level dB Low Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1036A 1036A 1036A Level dB Low Frequency RMS Deviation Change due to Equalisation 7 6 5 4 3 2 1 0 1 2 1036A 1036A 1036A Level dB Fig...

Страница 32: ...036A 1036A Level dB Midrange RMS Deviation Change due to Equalisation 7 6 5 4 3 2 1 0 1 2 1036A 1036A 1036A Level dB High Frequency RMS Deviation Before Equalisation 0 2 4 6 8 10 1036A 1036A 1036A Level dB High Frequency RMS Deviation After Equalisation 0 2 4 6 8 10 1036A 1036A 1036A Level dB High Frequency RMS Deviation Change due to Equalisation 7 6 5 4 3 2 1 0 1 2 1036A 1036A 1036A Level dB Fig...

Отзывы: