Fig 11,top: FR normalized to band average (0dB), in 25Hz wide band centered at 125Hz; Bottom: Level distribution in band; SD=4dB.

Stochastic nature of room acoustics 1

The Schroeder frequency:

 fc  ~2000·(T/V)0.5

This is an akuTEK article by M Skålevik           www.akutek.info main page   article site index page

(1)

Average spacing vs Quantization of level

N = 4p / 3·( f / c)3·V

Text Box: (7)

Number of modes below frequency f (Hz):

Text Box: (8)

Fig 8: FR in a 90m3 room, fc ~170Hz

Fig 9: FR level distribution

Fig 10: Distribution above and below fc

Text Box: (10)

The data in Table 11 does not indicate that the level distribution in the 70-170Hz region is of a different nature than in the 200Hz-2.8kHz region. We conclude that in this case, from the study of peak spacing and level distribution, the Schroeder Frequency does not represent an observable cross-over frequency between two significantly different frequency regions.

We shall need more data to continue this investigation, searching for an observable cross-over frequency.

 

70-170Hz

170-270Hz

270-370 Hz

200-2.8k Hz

Standard deviation

5.4dB

5.5dB

4.8dB

5.0dB

5% percentile

-8

-10

-6

-9

95% percentile

8

8

9

7

Deviation from reference

17%

8%

18%

0% (Ref)

First we shall study the level distribution (dark blue curve, Fig 9) of the FR in the 90 cubic meter room, above the Schroeder Frequency (200Hz to 2.8kHz), along with a normal distribution having the same average and standard deviation as the FR. The bell-shape reflects the Gaussian properties, confirming Schroeder’s theory. The FR in this region can be considered stochastic and we shall use its distribution as reference in the following low frequency analysis. Note the asymmetry (skewness= –1.1dB) of the level distribution as must be expected since the dB-scale goes from - to 6dB as sound pressure goes from 0 to 2 times the reference pressure at 0dB. The level distribution curve deviates 9% from the normal distribution curve (RMS deviation). 90% of the FR levels are between –9 and +7dB.

Now turning to the region below the Schroeder Frequency fc (Fig 10), the distribution curves (red curve 70-170Hz) seem to deviate considerably from the curve of the reference region 200-2.8k Hz (dotted curve). The deviation (RMS) is 17% in the 70-170Hz region, which is more than in the 170-270Hz region. However, the deviation in the 270-370Hz region is even greater (18%), meaning that the deviation below the Schroeder Frequency is not significant.

Standard deviation is 5.4Hz in the 70-170Hz region. Neither this or the 5% and 95% percentiles show significant differences from the values above the Schroeder Frequency. Table 11 below summarize the key properties of the different frequency regions.

Case study—Level distribution

Case study—Observed cross-over

In general, we want to support theory by observations. In particular, we want to support the theoretical cross-over frequency suggested by Schroeder by results with proper statistical significance. Procedure in steps as follows:

Search after an observable cross-over frequency

Step 6 in the procedure above has resulted in an observed cross-over frequency, detected by a best fit line method, at 107Hz, which is equal to 0.63 times the Schroeder Frequency of the room.

In further work, the result will be compared with measurements from other rooms, and discussed in relation to possible theoretical explanations.

For a start, figure 15 presents a case with no detectable crossover, for a single point in a different room than the case above. The volume of the room is 50m3 and the reverberation time in the possible crossover region is 0.39s.

Text Box: (11)

Fig 12: SD(dB) in 25Hz bands from one FR. Red dot indicates the 125Hz band.

Text Box: (12)
Text Box: (13)

Fig 13, top: SD (dB) in 25Hz bands, room average and 5% confidence intervals from 23 FRs in the 125Hz band (red dot). Bottom: Spatial distribution of SD’s in the 125Hz band

Text Box: (14)

Fig 14: Best fit lines for two frequency regions cross at f=107Hz

Fig 15: Similar to Fig 14, but a different room (V=50m3, T=0.39s), single measurement, no detectable cross over.

Text Box: (15)

Case study, observed Cross-Over-Frequency, conclusion

Discussion

In this theoretical discussion of significant regions of the frequency response, the frequency region above the Schroeder Frequency, for which Schroeder’s stochastic theory is valid (Gaussian level distribution, 4/T average spacing between maxima, etc.), will be referred to as the High Frequency (or HF) region. Below the HF region there is a transition region, below which the Low Frequency (LF) region is of a nature that is significantly different from the HF region.

Historically, the Schroeder Frequency has evolved as a limit above which the frequency response is likely to be dominated by stochastic processes. In 1954, Schroeder suggested a criteria of a 10-fold modal overlap on average, resulting in a critical frequency of 4000·(T/V)0.5. In his 1962 paper, he refers to several measurement reports indicating that the critical frequency may be as low as the current one (1), corresponding to an average 3-fold modal overlap, fig 16. However, this does not in it self exclude the possibility of stochastic nature at even lower frequencies. While the Schroeder Frequency serves well as a lower limit for the HF region, it is less suitable as a higher limit for the LF  region. The latter is expected to be dominated by average spacing between maxima that is wider than the modal bandwidth.

A possible two-sided cross-over frequency could be the frequency having an average mode spacing equal to the average modal bandwidth ~2.2/T, fig 17, which could be the center of a more or less wide and symmetrical transition region. However, below this cross-over frequency the average mode spacing is Df > 2.2/T, and in particular there will occur the mode spacing Df ~4/T, fig 18, which according to Schroeder is a characteristic feature of the HF region, contradicting this cross-over frequency candidate. In some cases it will in fact be impossible to distinguish an observation of mode spacing equal to 4/T and a maxima spacing equal to 4/T. This implies that in the Low Frequency region, the mode spacing must be on average wider than 4/T. Further, when average modal spacing is more than two times the modal bandwidth, Df >4.4/T, fig 19, more than 50% of the frequency response will depend on direct energy from the source, referred to below as <50% modal density. Corresponding wavelengths will be larger than the square root of the room absorption area, l>A0.5.

We may conclude that the transition region can be defined by several mode spacing criteria, with corresponding frequencies computed from (6). Moreover, the upper frequency limit of the region in which the near-field of a sound source is significantly influenced by reverberant sound pressure, i.e. the near-field extends into the reverberant field (krc < 1), falls in the very middle of this transition region. A summary is given in Table 20 below. Corresponding average mode spacing and modal overlap is illustrated in fig’s 16-19.

Continued in Part 2

First published 19.03.2009, latest change 12.10.2012

Related:

Tonal response in rooms (Microstructure of room acoustic FR), a presentation by M Skålevik
Schroeder Frequency Revisited  FA 2011 presentation  29.06.2011
Schroeder Frequency Revisited  (Forum Acusticum 2011 paper) 17.04.2011

Go to www.akutek.info main page     Go to article site index page

Mode spacing Df (Hz)~

f (Hz)~

f (Hz)~ re Schroeder fc

Comment

0.73/T

2000·(V/T)0.5

fc

3-fold overlap

2.2/T

1200·(V/T)0.5

0.6·fc

Single overlap

3.5/T

1000·(V/T)0.5

0.5·fc

Reverberant near-field

4/T

900·(V/T)0.5

0.45·fc

HF feature

>4.4/T

800·(V/T)0.5

< 0.4·fc

<50% modal density

Text Box: (17)
Text Box: (18)
Text Box: (19)

Fig’s 16-19: Illustration of two adjacent modes with mode spacing normalized to unity. Degree  of modal overlap and density decreasing from top to bottom, with 3-fold overlap in (16) to 50% modal density in (19).

Text Box: (16)
Text Box: (9)

Share this:

Step

Problem

Result

1

Define a proper statistic measure that is sensitive to Gaussian level distribution features, e.g. standard deviation, 5% percentile, skewness

Choice: The standard deviation (SD) within a proper band of the FR.

2

Choose a proper frequency bandwidth for the statistic measure, which is wide enough to contain the expected mode spacing and peak spacing (4/T) in the FR-curve, and still narrow enough to provide resolution in the region of interest, e.g. below 200Hz

Choice: 25Hz bandwidth, wider than mode spacing for f>37Hz and peak spacing for T>0.17s, resolution of 7 bands below 200Hz

3

With the chosen statistic measure, assign the computed value to the center frequency of each frequency band between 20Hz and 2.8 kHz of a frequency response measured in a point of a room

SD computed for each 25Hz band (like in fig 11), for each measured FR (plot in fig 12). In the following referred to as the FR band  deviation.

4

Since we are approaching the cross-over frequency as a room specific property (related to V and T), it is natural to continue the search by studying the room average of each frequency band.  

Room average from 23 points are plotted for each frequency band in fig 13, with error bars indicating the confidence interval (0.05 significance and spatial standard deviation over the 23 points).

5

Given the room average of the FR band deviation as plotted in (13), is it possible to find a cross-over frequency below which the FR band deviation shows a group of values that significantly differ from the group of values at higher frequencies? This may not be straightforward, given the gradually rising values towards lower frequencies in (13).

Based on the FR deviations in 25Hz bands (12) over 23 points in the room, with confidence intervals of 5% significance, it is not possible to detect a critical frequency below 800Hz.

6

Given the result of step 5, above, and noting the rising curve toward low frequencies, we turn to a different approach. We separate the data plotted in (13) in a low frequency (LF) part and a high frequency (HF) part, studying the best fit line of the LF part with an arbitrary upper frequency limit? For the HF region, data up to 2.8kHz is included.

By trial and error it was found that the best fit line of the LF region that  showed the optimum correlation (Pearson) of R2=0.85, meets the best fit line of the HF region at f=107Hz and level 4.1dB (fig 14).