A reverberation time measurement method, system, device, and medium

Abstract

The application provides a reverberation time measurement method, system, device and medium. The method comprises the following steps: calculating the energy value of each frequency band in a sound wave signal to obtain a frequency band energy value; creating a useful sound wave energy model, obtaining an ambient noise energy value of a current environment, inputting the frequency band energy value and the ambient noise energy value into the useful sound wave energy model to obtain a useful sound wave energy value; taking the maximum value of the useful sound wave energy value as a starting point, selecting the useful sound wave energy value in a preset range as a fitting range, fitting the useful sound wave energy value in the fitting range to obtain a fitting straight line; and performing reverberation calculation according to the slope of the fitting straight line to obtain a reverberation time. The application can more accurately estimate the noise power in the test stage, and can also shorten the working time of measurement.

Description

Technical Field

[0001] This invention relates to the field of loudspeaker technology, and in particular to a method, system, device and medium for measuring reverberation time. Background Technology

[0002] The generally accepted definition of reverberation time is: when the sound energy density decreases to 1 / 10 of its original value. 6 The reverberation time is the time required for the sound pressure level to decay by 60 dB. The reverberation time of a certain frequency band is the time required for the residual sound in a room to reach a stable state after the sound source stops emitting sound, and for the average sound energy density to decay from its original value to one part per million after repeated absorption by sound-absorbing materials. It is expressed as T60 or RT. Too short a reverberation time results in a dry, dull, and unnatural sound; too long a reverberation time makes the sound muddy; a suitable reverberation time results in a smooth and pleasant sound. Reverberation time is an important evaluation indicator for the quantitative estimation of sound energy in acoustic design.

[0003] The steady-state noise cutoff method is the most common method for measuring reverberation time. This method first establishes a stable sound field in the room using a sound source, then suddenly stops the sound source, and uses a microphone to measure the attenuation of the sound pressure level in the room, recording the attenuation curve. Finally, the reverberation time is calculated by taking 60 dB of sound pressure level to drop from the attenuation curve. However, this method has a drawback: it is difficult to achieve an absolutely quiet environment in a real room during the reverberation time measurement process; there is often ambient noise, and this ambient noise does not remain stable for a long time, but can only be considered stable for a short period.

[0004] In some existing methods for measuring reverberation time, after the sound pressure level has decayed to a minimum and stabilized, the microphone continues to collect a signal, which is then used as ambient noise to subtract the ambient noise power when estimating the reverberation time. However, in practice, this method sometimes fails to guarantee that the ambient noise remains stable for a considerable period after the sound pressure level has decayed to a minimum. This constraint affects the accuracy of the reverberation time measurement. Summary of the Invention

[0005] The purpose of this invention is to provide a method, system, device, and medium for measuring reverberation time, so as to reduce the impact of environmental noise estimation errors on the accuracy of reverberation time measurement.

[0006] To achieve the above objectives, in a first aspect, embodiments of the present invention provide a method for measuring reverberation time, the method comprising:

[0007] Calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value;

[0008] Create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model to obtain the useful sound wave energy value;

[0009] Starting from the maximum value of the useful sound wave energy, select a range of useful sound wave energy values ​​within a preset range as the fitting range, and fit the useful sound wave energy values ​​within the fitting range to obtain a fitting straight line;

[0010] The reverberation time is obtained by calculating the reverberation based on the slope of the fitted straight line.

[0011] Further, the step of selecting a range of useful sound wave energy values ​​within a preset range as the fitting range, starting from the maximum value of the useful sound wave energy value, includes:

[0012] The first difference is obtained by subtracting the maximum value of the useful sound wave energy value from the first energy correction value. In the time after the time corresponding to the maximum value of the useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range from the first difference is taken as the first useful sound wave energy value.

[0013] The difference between the first useful sound wave energy value and the second energy correction value is obtained to obtain the second difference value. In the time after the time corresponding to the first useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the second difference value is taken as the second useful sound wave energy value.

[0014] The time range between the first moment corresponding to the first useful sound wave energy value and the second moment corresponding to the second useful sound wave energy value is used as the fitting range.

[0015] Furthermore, after obtaining the fitted straight line, the process also includes:

[0016] Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error;

[0017] Determine whether the current mean square error is less than the preset initial mean square error;

[0018] If so, the current mean square error is multiplied by a preset adjustment step size factor to obtain the adjusted noise energy, and the environmental noise energy is updated to the sum of the adjusted noise energies of each iteration.

[0019] Further, after obtaining the fitted straight line, the process includes:

[0020] Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error;

[0021] Determine whether the current mean square error is less than the preset minimum mean square error;

[0022] If so, the minimum mean square error is updated to the current mean square error, and the slope of the fitted line is updated to the slope of the fitted line at the time of fitting.

[0023] Further, fitting the useful sound wave energy values ​​within the fitting range includes:

[0024] The energy value corresponding to each time point within the fitting range is used as the ordinate, and the corresponding time point is used as the abscissa to obtain the fitting point;

[0025] The least squares method is used to fit the fitted points.

[0026] Further, the step of inputting the frequency band energy value and the environmental noise energy value into the useful sound wave energy model to obtain the useful sound wave energy value includes:

[0027] The useful sound wave energy value is obtained by performing difference and logarithmic operations on the frequency band energy value and the environmental noise energy in the useful sound wave energy model.

[0028] Furthermore, the reverberation time is calculated based on the slope of the fitted straight line, specifically using the following formula:

[0029] t RT =-60 / (k·f) s )

[0030] Among them, t RT Let f be the reverberation time, k be the slope of the fitted line, and f be the reverberation time. s This is the preset sampling rate.

[0031] Secondly, embodiments of the present invention provide a reverberation time measurement system, the system comprising:

[0032] The energy calculation module is used to calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value.

[0033] The noise removal module is used to create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model, and obtain the useful sound wave energy value.

[0034] The linear fitting module is used to select a range of useful sound wave energy values ​​within a preset range as the fitting range, starting from the maximum value of the useful sound wave energy value, and to fit the useful sound wave energy values ​​within the fitting range to obtain a fitted linear line.

[0035] The reverberation time calculation module is used to calculate the reverberation time based on the slope of the fitted straight line.

[0036] Thirdly, embodiments of the present invention also provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the above-described method.

[0037] Fourthly, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the above-described method.

[0038] This invention provides a method, system, device, and medium for measuring reverberation time. The method includes: calculating the energy values ​​of each frequency band in a sound wave signal to obtain frequency band energy values; creating a useful sound wave energy model, obtaining the environmental noise energy value of the current environment, and inputting the frequency band energy values ​​and the environmental noise energy value into the useful sound wave energy model to obtain useful sound wave energy values; selecting a range of useful sound wave energy values ​​within a preset range as a fitting range, using the maximum value of the useful sound wave energy values ​​as a starting point, and fitting the useful sound wave energy values ​​within the fitting range to obtain a fitting straight line; and calculating reverberation based on the slope of the fitting straight line to obtain the reverberation time. This invention can more accurately estimate the noise power during the testing phase and can also shorten the measurement working time. Attached Figure Description

[0039] Figure 1 This is a flowchart illustrating a reverberation time measurement method according to an embodiment of the present invention;

[0040] Figure 2 This is a schematic diagram of the fast Fourier transform operation in a reverberation time measurement method according to an embodiment of the present invention;

[0041] Figure 3 This is a system block diagram of a reverberation time measurement system according to an embodiment of the present invention;

[0042] Figure 4 This is an internal structural diagram of the computer device in an embodiment of the present invention. Detailed Implementation

[0043] To make the objectives, technical solutions, and beneficial effects of this application clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Obviously, the embodiments described below are only part of the embodiments of the present invention and are used to illustrate the present invention, but are not intended to limit the scope of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0044] In one embodiment, such as Figure 1 As shown, the present invention provides a method for measuring reverberation time. The method includes:

[0045] S11. Calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value;

[0046] The microphone transmits sound waves through a sound source. Once the sound field reaches a steady state, the sound source stops emitting sound. The sound receiving device then transmits the sound at a sampling rate f. s The received signal is digitized and stored in memory. If 10 seconds of data are stored, if f... s =48kHz, the data length stored in memory over 10 seconds is 480KB. An N-point Fast Fourier Transform (FFT) is performed sequentially on the data in memory from the first data point to the last data point, as follows: Figure 2 As shown. During an FFT operation, N data points need to be fetched from memory each time. In the first FFT operation, N data points are fetched from the first data point in memory; subsequent fetches begin from the Mth data point after the previous starting point. That is, in the first FFT operation, data is fetched from the 1st data point to the Nth data point; in the second FFT operation, data is fetched from the (M+1)th data point to the Mth data point; and in the third FFT operation, data is fetched from the (2M+1)th data point to the (2M+N+1)th data point.

[0047] The number of FFTs that can be performed on the data in memory is:

[0048]

[0049] Where K is the number of FFTs that can be performed, L is the length of the data in memory, N is the length of the data truncated each time, and M is the distance between the starting points.

[0050] If we take N = 8192, then K = 983. That is, performing an FFT every 0.01s yields a total of 983 FFT results. Therefore, the reverberation time of each frequency band can be calculated individually with an accuracy of 0.01s.

[0051] The FFT result W of the k-th segment of the acoustic signal in memory k (i), i = 1, 2, ..., N, k = 1, 2, ..., K. A specific frequency band W is selected. k (i), i = N0, N0+1, ..., N1, calculate the energy of this specific frequency band:

[0052]

[0053] Where i represents a frequency point in the frequency band, and N0, N0+1, ..., N1 are the various frequency points in the frequency band.

[0054] The above-described process of energy calculation using Fast Fourier Transform can make fuller use of the characteristics of each frequency band in the sound wave signal, thus improving the accuracy of reverberation time calculation.

[0055] S12. Create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model, and obtain the useful sound wave energy value.

[0056] The method for calculating the useful sound wave energy value is to perform difference and logarithmic operations on the frequency band energy value and the environmental noise energy in the useful sound wave energy model to obtain the useful sound wave energy value.

[0057] Specifically, set the initial energy W of the ambient noise. noise =0, the noise energy W noise This will be updated in future updates. The energy of the frequency band minus the environmental noise energy W. noise Taking the logarithm gives the useful sound energy value E(k).

[0058]

[0059] Among them, W noise This refers to environmental noise energy.

[0060] Removing environmental noise here can reduce its impact on reverberation time measurement.

[0061] S13. Starting from the maximum value of the useful sound wave energy value, select the useful sound wave energy value within a preset range as the fitting range, and fit the useful sound wave energy value within the fitting range to obtain a fitting straight line.

[0062] The specific method for selecting the fitting range is as follows:

[0063] The first difference is obtained by subtracting the maximum value of the useful sound wave energy value from the first energy correction value. In the time after the time corresponding to the maximum value of the useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range from the first difference is taken as the first useful sound wave energy value.

[0064] The difference between the first useful sound wave energy value and the second energy correction value is obtained to obtain the second difference value. In the time after the time corresponding to the first useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the second difference value is taken as the second useful sound wave energy value.

[0065] The time range between the first moment corresponding to the first useful sound wave energy value and the second moment corresponding to the second useful sound wave energy value is used as the fitting range.

[0066] For example, take the maximum value max{E(k)} of the useful sound wave energy value E(k) and its corresponding time t. max From time t max Initially, the first energy correction value is 5dB, then the first difference is max{E(k)}-5dB. In E(k), the energy E(k0) that is closest to max{E(k)}-5dB is selected, and its corresponding time t0 is obtained.

[0067] Starting from time t0, select the energy E(k1) from E(k) that is closest to E(k0) - 20dB, corresponding to time t. 20 Alternatively, take the energy E(k2) that is closest to E(k0) - 30dB, corresponding to time t. 30 ;

[0068] The energy value corresponding to each time point within the fitting range is used as the ordinate, and the corresponding time point is used as the abscissa to obtain the fitting point;

[0069] The least squares method is used to fit the fitted points.

[0070] From the energy E(k0) at time t0 to the energy at time t 20 The energy E(k1) or the energy E(k0) from time t0 to time t 30 The energy E(k2) is denoted as y at each point. i , i = 1, 2, ..., L;

[0071] Fitted straight line The coefficients of the straight line for the j-th fitting are obtained by the least squares method, where y is the energy value at each time step and x is the time step number.

[0072]

[0073]

[0074] After obtaining the fitted line, the following steps are also included:

[0075] Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error;

[0076] Determine whether the current mean square error is less than the preset initial mean square error;

[0077] If so, the current mean square error is multiplied by a preset adjustment step size factor to obtain the adjusted noise energy, and the environmental noise energy is updated to the sum of the adjusted noise energies of each iteration.

[0078] For example, to calculate the mean square error of the j-th iteration,

[0079] Save the mean square error e1 calculated in the first time, that is, the initial mean square error; initialize the minimum mean square error e min = e1;

[0080] If the mean square error e in the j-th time j < e1, then adjust the noise energy where μ is the adjustment step factor, set by experiment, such as μ = 0.01; if not, then calculate the reverberation time;

[0081] After obtaining the fitting straight line, including:

[0082] Calculate the mean square error of multiple useful sound wave energy values within the fitting range to obtain the current mean square error;

[0083] Judge whether the current mean square error is less than the preset minimum mean square error;

[0084] If so, update the minimum mean square error to the current mean square error, and update the slope of the fitting straight line to the slope of the fitting straight line during this fitting.

[0085] For example, accumulate the adjusted noise energy for M times and update the noise energy Subtract the updated noise energy before taking the logarithm next time.

[0086] If the mean square error e in the j-th time j < e min , then update the minimum mean square error: e min = e j , and update the slope k = k j .

[0087] The above process of comparing and updating the mean square error can exclude the accumulated noise energy during the calculation process to better exclude the influence of noise energy on the calculation of reverberation time.

[0088] S14. Perform reverberation calculation according to the slope of the fitting straight line to obtain the reverberation time.

[0089] Specifically, it is calculated using the following formula:

[0090] t RT = -60 / (k · f s )

[0091] where t RT is the reverberation time, k is the slope of the fitting straight line, and f s is the preset sampling rate.

[0092] For the calculation of the reverberation time, other calculation formulas can also be used, which will not be elaborated here. It should be noted that the reverberation time here is the reverberation time relative to a certain frequency band.

[0093] This invention provides a method for measuring reverberation time, which can more accurately estimate noise power during the testing phase. Specifically, it includes: calculating the energy values ​​of each frequency band in a sound wave signal to obtain frequency band energy values; creating a useful sound wave energy model, obtaining the environmental noise energy value of the current environment, and inputting the frequency band energy value and the environmental noise energy value into the useful sound wave energy model to obtain useful sound wave energy values; using the maximum value of the useful sound wave energy value as a starting point, selecting a range of useful sound wave energy values ​​within a preset range as a fitting range, fitting the useful sound wave energy values ​​within the fitting range to obtain a fitted straight line; and calculating reverberation based on the slope of the fitted straight line to obtain the reverberation time. Compared with existing technologies, this invention also avoids the need for additional time to estimate environmental noise power, shortening the working time of reverberation measurement.

[0094] Based on the above-described reverberation time measurement method, this invention also provides a reverberation time measurement system, such as... Figure 3 As shown, the system includes:

[0095] Energy calculation module 1 is used to calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value.

[0096] The noise removal module 2 is used to create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model, and obtain the useful sound wave energy value.

[0097] The noise removal module 2 is also used to perform difference and logarithmic operations on the frequency band energy value and the environmental noise energy in the useful sound wave energy model to obtain the useful sound wave energy value.

[0098] The linear fitting module 3 is used to select a range of useful sound wave energy values ​​within a preset range as the fitting range, starting from the maximum value of the useful sound wave energy value, and to fit the useful sound wave energy values ​​within the fitting range to obtain a fitted linear line.

[0099] The linear fitting module 3 is also used to subtract the maximum value of the useful sound wave energy value from the first energy correction value to obtain the first difference value, and to take the useful sound wave energy value that is within the preset deviation range from the first difference value as the first useful sound wave energy value in the time after the time corresponding to the maximum value of the useful sound wave energy value.

[0100] The difference between the first useful sound wave energy value and the second energy correction value is obtained to obtain the second difference value. In the time after the time corresponding to the first useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the second difference value is taken as the second useful sound wave energy value.

[0101] The time range between the first moment corresponding to the first useful sound wave energy value and the second moment corresponding to the second useful sound wave energy value is used as the fitting range.

[0102] The linear fitting module 3 is also used to take the energy value corresponding to each time point within the fitting range as the vertical axis and the corresponding time point as the horizontal axis to obtain the fitting point;

[0103] The least squares method is used to fit the fitted points.

[0104] The reverberation time calculation module 4 is used to calculate the reverberation time based on the slope of the fitted straight line.

[0105] Reverberation time calculation module 4 also calculates the reverberation time using the following formula:

[0106] t RT =-60 / (k·f) s )

[0107] Among them, t RT Let f be the reverberation time, k be the slope of the fitted line, and f be the reverberation time. s This is the preset sampling rate.

[0108] For specific limitations regarding a reverberation time measurement system, please refer to the limitations regarding a reverberation time measurement method described above, which will not be repeated here. Each module in the above system can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in hardware or independently of the processor in a computer device, or stored in software in the memory of a computer device, so that the processor can call and execute the corresponding operations of each module.

[0109] Figure 4 An internal structural diagram of a computer device is shown in one embodiment. This computer device may specifically be a terminal or a server. Figure 4 As shown, the computer device includes a processor, memory, network interface, display, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The network interface is used to communicate with external terminals via a network connection. The display screen can be an LCD screen or an e-ink display screen. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad mounted on the computer device casing, or an external keyboard, touchpad, or mouse.

[0110] Those skilled in the art will understand that Figure 4 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computing devices may include more or fewer components than shown in the diagram, or combine certain components, or have the same component arrangement.

[0111] In summary, this invention provides a method, system, device, and medium for measuring reverberation time. The method includes: calculating the energy values ​​of each frequency band in a sound wave signal to obtain frequency band energy values; creating a useful sound wave energy model, obtaining the environmental noise energy value of the current environment, and inputting the frequency band energy values ​​and the environmental noise energy value into the useful sound wave energy model to obtain useful sound wave energy values; using the maximum value of the useful sound wave energy value as a starting point, selecting a range of useful sound wave energy values ​​within a preset range as a fitting range, fitting the useful sound wave energy values ​​within the fitting range to obtain a fitted straight line; and calculating reverberation based on the slope of the fitted straight line to obtain the reverberation time. This invention can more accurately estimate the noise power during the testing phase and can also shorten the measurement working time.

[0112] The various embodiments in this specification are described in a progressive manner. For directly identical or similar parts of the embodiments, refer to each other. Each embodiment focuses on its differences from other embodiments. In particular, the system embodiments are basically similar to the method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments. It should be noted that the technical features of the above embodiments can be combined arbitrarily. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as the combination of these technical features does not contradict each other, it should be considered within the scope of this specification.

[0113] The embodiments described above are merely preferred embodiments of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various improvements and substitutions without departing from the technical principles of this invention, and these improvements and substitutions should also be considered within the scope of protection of this application. Therefore, the scope of protection of this patent application should be determined by the scope of the claims.

Claims

1. A method for measuring reverberation time, characterized in that, The method includes: Calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value; Create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model to obtain the useful sound wave energy value; The first difference is obtained by subtracting the maximum value of the useful sound wave energy value from the first energy correction value. In the time after the time corresponding to the maximum value of the useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range from the first difference is taken as the first useful sound wave energy value. The difference between the first useful sound wave energy value and the second energy correction value is obtained to obtain the second difference value. In the time after the time corresponding to the first useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the second difference value is taken as the second useful sound wave energy value. The time range between the first moment corresponding to the first useful sound wave energy value and the second moment corresponding to the second useful sound wave energy value is used as the fitting range. The useful sound wave energy values ​​within the fitting range are fitted to obtain a fitting straight line. Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error; Determine whether the current mean square error is less than the preset minimum mean square error; If so, then update the minimum mean square error to the current mean square error, and update the slope of the fitted line to the slope of the fitted line at this time. Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error; Determine whether the current mean square error is less than the preset initial mean square error; If so, the current mean square error is multiplied by a preset adjustment step size factor to obtain the adjusted noise energy, and the environmental noise energy is updated to the sum of the adjusted noise energies of each iteration. The reverberation time is obtained by calculating the reverberation based on the slope of the fitted straight line.

2. The reverberation time measurement method according to claim 1, characterized in that, The fitting of the useful acoustic energy values ​​within the fitting range includes: The energy value corresponding to each time point within the fitting range is used as the ordinate, and the corresponding time point is used as the abscissa to obtain the fitting point; The least squares method is used to fit the fitted points.

3. The method for measuring reverberation time according to claim 1, characterized in that, The step of inputting the frequency band energy value and the environmental noise energy value into the useful sound wave energy model to obtain the useful sound wave energy value includes: The useful sound wave energy value is obtained by performing difference and logarithmic operations on the frequency band energy value and the environmental noise energy in the useful sound wave energy model.

4. The method for measuring reverberation time according to claim 1, characterized in that, The reverberation time is calculated based on the slope of the fitted straight line, using the following formula: ； in, For reverberation time, The slope of the fitted line, This is the preset sampling rate.

5. A reverberation time measurement system, characterized in that, For performing a reverberation time measurement method as described in any one of claims 1-4, the system comprises: The energy calculation module is used to calculate the energy value of each frequency band in a sound wave signal to obtain the frequency band energy value. The noise removal module is used to create a useful sound wave energy model, obtain the environmental noise energy value of the current environment, input the frequency band energy value and the environmental noise energy value into the useful sound wave energy model, and obtain the useful sound wave energy value. The linear fitting module is used to subtract the maximum value of the useful sound wave energy value from the first energy correction value to obtain the first difference value. In the time after the time corresponding to the maximum value of the useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the first difference value is taken as the first useful sound wave energy value. The difference between the first useful sound wave energy value and the second energy correction value is obtained to obtain the second difference value. In the time after the time corresponding to the first useful sound wave energy value, the useful sound wave energy value that is within the preset deviation range of the second difference value is taken as the second useful sound wave energy value. The time range between the first moment corresponding to the first useful sound wave energy value and the second moment corresponding to the second useful sound wave energy value is used as the fitting range. The useful sound wave energy values ​​within the fitting range are fitted to obtain a fitting straight line. Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error; Determine whether the current mean square error is less than the preset minimum mean square error; If so, then update the minimum mean square error to the current mean square error, and update the slope of the fitted line to the slope of the fitted line at this time. Calculate the mean square error of multiple useful acoustic wave energy values ​​within the fitting range to obtain the current mean square error; Determine whether the current mean square error is less than the preset initial mean square error; If so, the current mean square error is multiplied by a preset adjustment step size factor to obtain the adjusted noise energy, and the environmental noise energy is updated to the sum of the adjusted noise energies of each iteration. The reverberation time calculation module is used to calculate the reverberation time based on the slope of the fitted straight line.

6. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 4.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 4.

Images