Sound Source Omnidirectional Positioning Calibration Method Based on Microphone Observation Angle

To get source azimuth from microphone observation angle of view in a complex real environment, this article, on the basis of the analysis of geometric positioning method, established a seven-element microphone array model and proposed a sound source omnidirectional positioning calibration method based on microphone observation angle. By using a seven-element array to invert the position and angle of a sound source, the relative time delay value of a pair of microphones on the vertical axis of the coordinate system is used to determine the elevation angle polarity and realize the omnidirectional sound source positioning. The array parameters, sound velocity, array size, horizontal deflection angle, elevation angle, and sound source are analyzed, and the error method is proposed. The sound source data was measured using the microphone array perspective, and a new Cartesian coordinate system was established based on the observation angle of view for omnidirectional positioning calibration of the sound source. The simulation results show that the position error of the method is about 0.01% and the angle error is about 0.005%, with high calibration accuracy. The actual measurement results show that this method can effectively calibrate the sound source azimuth, the error rate of the source coordinates is around 10%, the horizontal declination angle error is less than 5%, and the elevation angle error is less than 8%. Appropriately increasing the spacing of the array will have a better calibration effect in an actual complex experimental environment.


Introduction
Sound is everywhere, along with people's daily life.In recent years, with industrial, civil, and military field, a sharp increase in demand for the positioning system, voice recognition systems [1,2], and sound source detection system [3] has become a new research hotspot.A sound localization system is actually the first to receive through a device the sound source signal, then the received signal, and sound source signal processing technology is adopted to improve the relevant technology to identify the target source location [4,5].At present, the sound source positioning technology is mainly divided into two types [6]: active sound source positioning and passive sound source positioning.Compared with active sound source positioning technology, passive sound source positioning technology has strong concealment and anti-interference.This paper is based on a passive sound source positioning algorithm.It is difficult to locate the sound source rapidly and precisely because of the current sound source localization method, which directly affects the azimuth calibration of the sound source.In order to solve the above problems, we should improve the sound source localization accuracy in a complex environment and introduce a sound source azimuth calibration method, which has great practical significance.
Traditional sound source localization technology is mainly used in speech signal processing and military field [7,8].Due to the foreign investment in the technology, the starting point is earlier than our country, so the foreign application in this technology is relatively mature.Equipment such as the UK Ferranti company Picker helicopter ISC company equipment alarm system, the United States' PALS ERAM remote passive acoustic positioning system (GPS), and the United States anti-vehicle intelligent mines [9], and so on, all the success of different applications in the field of military equipped with our R&D and design of passive acoustic positioning system, has the higher performance of the system and has more ability of target detection and recognition.In the field of speech signal processing, foreign researchers have conducted a lot of in-depth research.Oppenheim and Lim [10] and Oppenheim et al. [11,12] also demonstrated the contribution of phase-only signal reconstruction on human perception with regard to the quality of the reconstructed signal, exemplified for audio, image, crystallographic, and speech processing applications.Liu et al. studied the perceptual importance of the phase spectrum specifically on the intervocalic stop consonants, and they also reported the strong dependency of the perception of intervocalic stops on phase information [13].In 1985, Flanagan, in a large conference speech enhancement in the application of the speech signal processing technology, obtained the good actual effect.In 1992, Silverman applied the technology to the speaker's speech recognition system, and it also achieved good results.In addition, the research on the acoustic source localization algorithm also has a great breakthrough, in which the beamforming method based on the microphone array [14,15] appears early and can be applied in practice.In view of the limited resolution of the problems of the algorithm, after in-depth analysis and research, scholars put forward the algorithm of high resolution [16,17], such as the minimum variance method, the characteristics of the structure method, and the signal subspace method [18].
Further research by researchers Paliwal and Alsteris [19] elaborates and demonstrates the importance of phase information in human hearing.Alsteris and Paliwal [20] show that the short-term phase spectrum can significantly improve the speech articulation in a small window of 20 to 40 milliseconds, which is very useful.Alsteris and Paliwal [21] based on automatic speech recognition (ASR) proved that in all cases, the resolution is superior to short weight of stimulus provided by the clarity, short time of amplitude, and phase spectrum information knowledge which lead to superior human speech recognition performance.Paliwal et al. [22] reported based on the analysis of the match or not match window for AMS during the processing of the amplitude and phase spectrum estimation and phase spectrum compensation (PSC) method that the minimum mean square error (MMSE) short-term spectrum amplitude (STSA) estimation to replace the noise amplitude spectrum and improve the voice quality is estimated based on the accurate phase spectrum.The overview of microphone array speech processing is proposed by Brandstein and Ward [23] in 2001, and it is considered that the phase difference between channels is the main source of information for delay estimation and source location application.Kudriashov [24] proposed the basic method for the three-dimensional positioning of sound source in distance-span-elevation coordinates by combining the time delay estimation results.
Although the research in the field of sound source localization technology in our country lags behind that of advanced countries but with the rapid development of modern technology in our country, the researchers are also actively involved in the study of the technology.Under the support of the National Happiness Technology 863, some institutions in our country have carried out the research on the smart bomb [25].There is no lack of some technical strength of the high-tech enterprises and study of the sound source localization technology, represented by the company, establishing a sound source localization based on microphone array product research and development projects.Ju et al. [26] combine the characteristics of the subarray algorithm and the microphone array signal processing in the far-field narrowband signal and propose the sound source 3D localization algorithm based on the microphone array.Ju et al. [27] extended the classical MUSIC algorithm to the speech array processing to realize the three-dimensional localization of the sound source.Ding et al. [28] studied the discharge characteristic of the air gap to the level of the high altitude air gap discharge voltage calibration which provides basis, but with the increase of altitude, rod-plate and barbar clearance operation impact discharge voltage is reduced, and bar-bar discharge gap voltage margin lower than the error bar-plate gap is enlarged.In Xing et al.'s paper [29], in the data of the atmospheric electric field calibration method, in order to improve the consistency of the electric network data, the scope of the thunderstorm monitoring is expanded, to effectively solve the problem of the observatory data unification and avoid a severe distortion in electric network; the calibration method is limited to an altitude of the vertical direction but is not from the perspectives of atmospheric electric field instrument observation of thunder cloud monitoring.Wang et al. [30] put forward a kind of array by using auxiliary array element each dependent on an amplitude-phase error since the new method of correction is applicable to arbitrary array geometry, and its computational complexity is small, but only limited to the two-dimensional space.
Based on the acoustic source geometry localization algorithm, this paper studies the sound source omnidirectional positioning calibration method based on microphone observation angle.Based on the array, the relative time delay of a pair of microphones on the vertical axis is used to realize the omnidirectional sound source location.The array parameters are analyzed, and the improved method is put forward, using the data from the angle which measured sound source array based on observation angle of view, on the basis of the forward modeling data, carried out omnidirectional positioning calibration of the sound source, and analyzed the simulation data and actual measurement experiment in a complex real environment.

Seven-Element Microphone Array Model
The sound source geometric localization algorithm [31][32][33] and position calibration methods [29,30,34] are based on a microphone array model [23,[35][36][37][38]; the establishment of a reasonable set microphone array model and array element spacing can better match algorithm, reduce the complexity 2 Complexity of the calculation, and improve the efficiency of the operation and sound comprehensive orientation calibration precision.A seven-element microphone array model was established, and seven microphones were used to compose the array as shown in Figures 1-4: Figure 1 shows the array model where the sound source S is not on the x-0-z plane, 0-y-z plane, and z-axis; Figure 2 shows the array model of the sound source S on the x-0-z plane; Figure 3 shows the array model of the sound source S on the 0-y-z plane; Figure 4 shows the array model of the sound source S on the z-axis.The sound source propagation velocity is c, and the space position of sound source S is x, y, z , and r is the distance from the sound source to the array center M0.The time for S to propagate to microphones M0(0,0,0), M1(a,0,0), M2(0,a,0), M3(−a,0,0), M4(0,−a,0), M5(0,0,a), and M6(0,0,−a) is t0, t1, t2, t3, t4, t5, and t6, respectively.Based on the model, the relative delay values of 5 groups are set: T 1 = t1 − t0, T 2 = t2 − t0, T 3 = t3 − t0, T 4 = t4 − t0, and T 5 = t6 − t5.H 1 is the horizontal deflection of the sound source S relative to the x-positive half axis, and the range is between 0 and 360 degrees.H 2 is the elevation angle formed by the sound source S and the x-y-0 plane (range between −90 and 90 degrees).If the value of H 2 is positive, that is, T 5 is greater than zero, it means that the sound source is above the x-y-0 plane.If it is negative, that is, T 5 is less than zero, it means that the sound source is below the x-y-0 plane; this is the key to realizing the omnidirectional sound source positioning calibration.

Sound Source Geometric Localization Algorithm
A suitable microphone array model can be established in order to obtain a smaller error and further prepare the omnidirectional sound source positioning calibration.The localization algorithm of sound source azimuth needed to undertake choosing according to the array structure and adopt different localization algorithms will have different positioning precisions and computing speeds, accompanied by the algorithm of the computation size and position precision of algorithm.
Based on the location method of arrival delay, there is a maximum likelihood method based on search [39][40][41] and geometric localization algorithm [31][32][33]: the search method has a large amount of calculations and does not apply to real-time systems; if the sound source geometric localization algorithm is adopted, the algorithm has high positioning accuracy and the arrival delay itself has strong anti-interference, which can reduce the influence of the attenuation and interference of the actual sound propagation process on the accuracy of the sound source omnidirectional positioning calibration.3 Complexity 3.1.Sound Source S Is Not on the x-0-z Plane, 0-y-z Plane, and z-Axis.According to Figure 1, the seven-element microphone array model whose sound source is not on the x-0-z plane, 0-yz plane, and z-axis represents the distance from the sound source S to M0, M1, M2, M3, and M4, and its expression is as follows: x + a 2 + y 2 + z 2 = r + cT 3 2 , x, y, and z are the sound source coordinate parameters of the sound source geometric localization algorithm based on the array of seven-element microphones, and formula (1) is the linear analysis basis of the omnidirectional positioning calibration method of sound source.According to (1), it can be concluded that x′ is the simulation value of the sound source coordinate parameter x, and parameter y′ is the simulation value of the source coordinate parameter y.
According to the actual situation, since the magnitudes of sum e = T 1 + T 3 and m − e = T 2 + T 4 and sum m are both about 10 −6 , then the order of magnitude c 2 me is 10 −7 , and the order of magnitude 4a 2 − c 2 n is only 10 1 ; relatively speaking, c 2 me can be ignored, and the same reason c 2 m m − e can be ignored, so (2) can be simplified for Using (1) again, we can get the expression of the distance from the sound source to the center of the array M0, that is, In the formula, a is the distance between the array elements, c is the sound source propagation speed, and r′ is the simulation parameter of the sound source coordinate parameter r. a and c are theoretical constants.Since (3) and ( 4) are only related to relative delay value T 1 to T 4 , they can be used in simulation experiments for the sound source omnidirectional positioning calibration method and error analysis for the localization accuracy of sound source.
According to Figure 1, the seven-element microphone array model where the sound source S is not on the x-0-z plane, 0-y-z plane, and z-axis can be used to derive the correlation equation: 5) based on the polar coordinate relationship of the three-dimensional geometry is the theoretical basis for solving the horizontal declination angle and the elevation angle.Among them, H 1 ′ is the simulation value of the horizontal declination angle H 1 in the omnidirectional positioning of the sound source, H 2 ′ is the simulation value of the elevation angle H 2 , r ′ is the simulation value of the acoustic source coordinate parameter r, and z ′ is the simulation value of the sound source coordinate parameter z.
Using ( 5), the value of horizontal declination H 1 ′ is According to the positive and negative of T 5 , the elevation value H 2 ′ can be obtained by polarity judgment.Using (5) again, the following expression can be obtained: In the formula, T 5 is the time delay between the sound source S arriving at the microphones M5 and M6.After getting H 2 ′ according to the polarity judgment of T 5 , ( 5) is used to determine the sound source coordinate parameter z′, and then the parameters x′ and y′ are combined to realize the omnidirectional positioning of the sound source and further be used for omnidirectional positioning calibration of the sound source.
Based on the seven-element microphone array model shown in Figure 1, the theoretical angle value can be solved, and its expression is as follows: 4 Complexity Sound Source S Is on the x-0-z Plane.According to the seven-element microphone array model of the sound source S on the x-0-z plane in Figure 2, because only the value of the sound source coordinate parameter y′ is changed, when y′ is equal to 0 and x′ and z′ are not 0, which is equivalent to only changing the time delay value T 1 to T 5 , using ( 3) and ( 4), the expression of the sound source coordinate parameters x′, y′, and z′ is as follows: x Based on Figure 2, the value of horizontal declination H 1 ′ can be obtained, and its expression is as follows: 10), it can be seen that the value of H 1 ′ is completely determined by the positive and negative values of the sound source coordinate parameter x ′ , which conforms to the solid geometry shown in Figure 2.
According to the time delay value T 5 and the positive and negative of the sound source coordinate parameter x ′ , the elevation value H 2 ′ can be judged by the polarity and the value H 2 ′ can be obtained by using formula (5).Its expression is as follows: Similarly, after determining H 2 ′ according to the polarities of T 5 and x ′ , use (5) to determine the sound source coordinate parameter z ′ and combine parameters x ′ and y ′ to achieve omnidirectional localization of the sound source and further use it for the sound source omnidirectional positioning calibration.Using (8), the theoretical values of horizontal declination H 1 and elevation H 2 can be solved.
3.3.Sound Source S Is on the 0-y-z Plane.According to the seven-element microphone array model of the sound source S on the 0-y-z plane in Figure 3, because only the value of the sound source coordinate parameter x′ is changed, when x′ is equal to 0 and y ′ and z ′ are not 0, which is equivalent to only changing the time delay value T 1 to T 5 , again using ( 3) and ( 4), the expression of the sound source coordinate parameters x′, y′, and z′ is as follows: Based on Figure 3, the value of horizontal declination H 1 ′ can be obtained, and its expression is as follows: From (13), it can be seen that the value of H 1 ′ is completely determined by the positive and negative values of the sound source coordinate parameter y ′ , which conforms to the solid geometry shown in Figure 3.
According to the time delay value T 5 and the positive and negative of the sound source coordinate parameter y ′ , the elevation value H 2 ′ can be judged by the polarity and the value H 2 ′ can be obtained by using formula (5).Its expression is as follows: Similarly, after determining H 2 ′ according to the polarities of T 5 and y′, use (5) to determine the sound source coordinate parameter z ′ and combine parameters x ′ and y ′ to achieve omnidirectional localization of the sound source and further use it for the sound source omnidirectional 5 Complexity positioning calibration.Again using (8), the theoretical values of horizontal declination H 1 and elevation H 2 can be solved.
3.4.Sound Source S Is on the z-Axis.According to the sevenelement microphone array model of the sound source S on the z-axis in Figure 4, because only the values of the sound source coordinate parameters x′ and y′ are changed and x′ and y′ are both 0 and z′ are not 0 at this time, which is equivalent to changing only the time delay value T 1 to T 5 .By using (3) and ( 4), the expressions for parameters x ′ , y ′ , and z ′ are as follows: Based on Figure 4, the value of horizontal declination H 1 ′ can be obtained, and its expression is as follows: According to the positive and negative values of the time delay value T 5 , the values of the elevation angle value H 2 ′ and the sound source coordinate parameter z ′ can be determined and obtained, and the related expressions are as follows: Again using (8), the theoretical values of horizontal declination H 1 and elevation H 2 can be solved.

Error Analysis of Sound Source Positioning Calibration and the Way to Reduce the Error
Through the analysis of azimuth positioning formulas (3) and ( 4), it can be known that the azimuth positioning performance of the target sound source is related to the time delay estimation error, the size of the microphone array, and the effective speed of sound.
The statistical characteristics of time delay are the same, assuming that the standard deviation of the time delay estimation error is ξ T and independent of each other.

Horizontal Declination Angle Error Analysis. The estimation error of the horizontal declination angle H 1 caused by the delay estimation error is
Using H 1 = arctan T 4 − T 2 / T 3 − T 1 , one can obtain the partial derivative of the horizontal declination angle H 1 relative to the delay value T 1 to T 4 , and its expression is as follows: Equation ( 19) can be used to derive the equation for the estimation error ξ H 1 T of horizontal declination angle H 1 caused by the delay estimation error, which is Equation ( 20) is simulated.When the estimation error at the time delay is 1 μs, the error curve shown in Figure 5 can be obtained.
As can be seen from Figure 5, the accuracy of the horizontal declination angle estimation has nothing to do with the horizontal declination angle value, but it is affected by the time delay estimation error, the microphone array size, and the elevation angle value.When the size of the microphone array a increases, the estimation error of the horizontal declination angle H 1 decreases; between the estimation errors of the horizontal declination angle H 1 and the elevation angle H 2 , the estimation error of the horizontal 6 Complexity declination angle decreases as the elevation angle decreases.Therefore, this is a more accurate measurement of lowaltitude and ground sound sources than aerial sound source measurements.
When the array spacing a is equal to 2 meters and the elevation angles are 15 degrees, 30 degrees, 45 degrees, and 60 degrees, the relationship between the estimation error of horizontal declination angle H 1 and the time delay estimation error is shown in Figure 6.
From Figure 6, we can see that when the value of the elevation H 2 is certain, the estimation error of horizontal declination angle H 1 increases linearly with the increase of the time delay estimation error; as the elevation angle value H 2 increases, the linear variation between the two increases more sharply, indicating that the estimation error of the delay has a greater influence on the measurement accuracy of the horizontal declination angle H 1 , and its measurement accuracy may drastically decrease.

Elevation Error Analysis. The estimation error of the elevation angle H 2 caused by the delay estimation error is
Using one can find the partial derivative of the elevation angle H 2 with respect to the delay value T 1 to T 4 , which has the following expression: Equation ( 22) can be used to derive the equation for the estimation error ξ H 2 T of the elevation angle H 2 caused by the delay estimation error, which is Equation ( 23) is simulated.When the estimation error at the time delay is 1 μs, the error curve shown in Figure 7 can be obtained: As can be seen from Figure 7, the estimation accuracy of the elevation angle has nothing to do with the horizontal declination value, but it is also affected by the time delay estimation error, the microphone array size, and the elevation angle value.When the size of the microphone array a increases, the estimation error of the elevation angle H 2 decreases; the elevation angle H 2 itself changes the angle of elevation, and as the elevation angle increases, the estimation error of the elevation angle decreases, that is, the accuracy of elevation angle measurement will increase.
When the array spacing a is equal to 2 meters and the elevation angles are 15 degrees, 30 degrees, 45 degrees, and 60 degrees, the relationship between the estimation error of elevation angle H 2 and the time delay estimation error is shown in Figure 8.

Complexity
Figure 8 shows that when the value of elevation angle H 2 is certain, its own estimation error will increase linearly with the increase of time delay estimation error; as the value of elevation angle H 2 decreases, the linear change between the two will become more severe, indicating that the time delay estimation error has a greater influence on the accuracy of the elevation angle H 2 , and its measurement accuracy may drastically decrease.

Error Analysis of the Distance from the Sound Source to
the Microphone Array Center.The standard deviation of the estimation error of the distance r from the sound source to the center of the array is Once again using (1), the expression of any relative delay can be obtained as follows: Among them, i is equal to 1 to 4, and T i is the relative delay of the sound source S to reach the two microphones.From (4), it can be seen that the expression of the distance r from the sound source to the array center M0 is related to T i , and therefore, T i is solved first.
Let a/r 2 − 2a/r cos H 2 cos H 1 − π/2 i − 1 = x, a → 2a; equation ( 25) can be expressed as Taylor's formula is used to expand ( 26) and ignore highorder infinitesimal parts, which is Let cos H 2 cos H 1 − π/2 i − 1 = A, where i is equal to 1 to 4 and, at the same time, summed the two sides of the above equation.Because of ∑ 4 i=1 A = 0 and ∑ 4 i=1 A 2 = cos 2 H 2 , according to (27), the expression is as follows: According to (4) about the distance r from the sound source to the center of the array M0, using (28), r has the following expression: Using (29) to find the relevant partial derivative and substituting it into (24), we can conclude that In (30), the relative ranging error ξ rT of the sound source is related to the accuracy of time delay estimation, the size of the microphone array, the distance r from the sound source to the center of the array, and the elevation angle.However, irrelevant to the horizontal deviation angle, the variation of the horizontal deviation angle does not cause the change of relative ranging error.8 Complexity Equation ( 30) is simulated, the elevation angle H 2 is 45 degrees, the delay estimation error is 1 microsecond, and the error curve shown in Figure 9 can be obtained.
As can be seen from Figure 9, when the time delay estimation accuracy and the microphone array size are certain, the relative ranging error increases as the distance r from the sound source to the array center increases, and the relative ranging performance based on the array degrades.This situation is suitable for near-field sound source location; increasing the array size a, the relative range error is reduced, and the relative ranging accuracy will be improved, but the array size should not be less than 2 meters; otherwise, the relative ranging error will rise sharply.
When the elevation angle H 2 is 30 degrees and 60 degrees, respectively, the relative ranging error is 0.05 radians and 0.1 radians, respectively, and the distance between the sound source and the center of the array is 100 meters; the relationship between the time delay estimation error and the microphone array size is shown in Figure 10.
As can be seen from Figure 10, when the relative ranging error, elevation angle H 2 , and microphone array size a are combined and the microphone array size is constant, the relative ranging error has a greater influence on the time delay estimation error than the elevation angle size; when the size is large, the increase of the array size a will cause an exponential increase in the error of the time delay estimation; especially when the size of the array is larger than 8 meters, the error increases more significantly.Therefore, it is recommended that the array size should not be too large.
When the elevation angles H 2 are 30 degrees and 30 degrees, respectively, and the relative ranging error is 0.05 radians and 0.1 radians, respectively, and the array size is 2 meters, the relationship between the time delay estimation error and the distance from the sound source to the center of the array is shown in Figure 11.
It can be seen from Figure 11 that when the distance r from the sound source to the center of the array is small, the relationship between it and the time delay estimation error changes drastically.In particular, when the sound source is very close to the center of the array, the time delay estimation error tends to be infinite; this agrees with the actual situation.This indicates that it is not appropriate to locate the sound source that is too close to the center of the array.Before integrating the error analysis, when the distance between the sound source and the array center reaches 8 to 10 meters, it will be better for the sound source to be fully calibrated.

Error Analysis of the Effect of Sound Velocity on
Calibration Accuracy of Sound Source Location.The speed at which sonars travel in different media is different.At a temperature of 273.16K and an atmospheric pressure of 105 Pa, the speed of sound c in the air is equal to 331.45 plus or minus 0.05 meters per second.In the actual environment, the influence of the ambient temperature on the speed of sound can not be ignored.Under the condition that the atmospheric pressure and other conditions remain relatively stable, the formula for calculating the sound velocity in the air at different temperatures is Among them, c 0 = 331 45 meters per second and T 0 = 273 16 K.
According to (31), relevant conclusions can be drawn: the speed of sound propagation in air is different at different temperatures, and the speed of sound increases with increasing temperature.Therefore, it is necessary to perform temperature measurement according to different environments during actual application to improve calibration accuracy.

The Way to
Reduce the Error.Based on the error analysis of parameters such as the horizontal declination angle, elevation angle, the distance between the sound source and the center of the array, and sound velocity based on the sevenelement microphone array model, the following five-point error reduction method is proposed: (1) From the analysis of the estimation error of the horizontal declination angle H 1 and the elevation angle H 2 , the microphone array size a can be appropriately increased to reduce the estimation error of both, thereby improving the measurement accuracy of both (2) From the analysis of the estimation error of the distance from the sound source to the center of the array, the size of the array can be appropriately increased to reduce the relative ranging error, thereby improving the relative ranging accuracy, and the general array size should not be less than 2 meters (3) From the curve of the relation between array size a and time delay estimation error, it can be seen that the general array size should not be larger than 10 meters (4) It can be seen from the plot of the distance r from the sound source to the center of the array element and the time delay estimation error that when r reaches 8 to 10 meters, the sound source has a better omnidirectional positioning calibration effect

Sound Source Omnidirectional Positioning Calibration
Consider that in the actual environment, locating the sound source based on the microphone array cannot be directly transmitted to the spatial location information of the actual sound source of the observation party.Using the omnidirectional data of the sound source measured from the perspective of the array of seven-element microphones, a new coordinate system was rebuilt based on the observation angle based on the microphone, and the azimuth position of the sound source was calibrated.
5.1.Sound Source Positioning Calibration Model.Let the observation angle N be located in the coordinate system of the original microphone array, and denote it as X, Y, Z .
The sound source positioning calibration model based on the observation angle of the microphone is shown in Figure 12.
According to Figure 12, we convert X, Y, Z to (0,0,0) and use this as a starting point to create a new coordinate system, converting x′, y′, z′ to x′ − X, y′ − Y, z′ − Z , which can be equivalent to X ′ , Y ′ , Z ′ .Let the angle between the sound source S and the positive half of x1 in the new coordinate system be the horizontal declination angle h 1 in the range of 0 to 360 degrees, and the angle between the source S and the plane of the new coordinate system x1, y1, 0 is the elevation angle h 2 (range from −90 to 90 degrees).If the value of h 2 is positive, it means that the sound source is above the plane of N; if it is negative, it is below it; the distance from X ′ , Y ′ , Z ′ to observation angle N is R ′ .

Sound Source Positioning Calibration
Method.According to Figure 12, based on the sound source positioning calibration model, using the measured omnidirectional positioning coordinates of the sound source, no matter where the 10 Complexity observation angle N is located, the sound source positioning calibration can be performed.Its expression is as follows: After using (32) to determine the omnidirectional positioning coordinates X′, Y′, Z′ of the sound source based on the observation angle, according to the actual situation, using the analytical method of the three-dimensional geometry, determine the positive and negative coordinates of the sound source parameters X ′ and Y ′ , and obtain the horizontal declination angle h 1 and the elevation angle h 2 .The value of the expression is as follows: From (33), it can be seen that for any spatial position of the sound source, the horizontal declination angle and the elevation angle value can be obtained.In combination with (32), the omnidirectional positioning calibration of the sound source can be realized.

The Simulation Results
Based on the array of seven-element microphones, the sound source geometric localization algorithm was introduced into the model, the sound source azimuth was simulated on the basis of the forward data, the sound source data measured from the array observation angle was used for sound source localization calibration, and two experiments were performed: Experiment 1 is a simulative experiment of a omnidirectional localization algorithm based on a sevenelement microphone array, and Experiment 2 is a simulation experiment of an omnidirectional positioning calibration method for a sound source.
In Experiment 1, in order to verify the positioning accuracy of the sound source geometric positioning algorithm, a seven-element microphone array model was used to simulate it.According to Figures 1-4, the parameters are set: c is equal to 340 meters per second, and a is equal to 2 meters.The four forward datasets are as follows: x equals 80 meters, y equals −60 meters, and z equals 120 meters; x equals −50 meters, y equals 0 meters, and z equals −130 meters; x equals 0 meters, y equals 70 m, and z equals 150 meters; and x equals 0 meters, y equals 0 meters, and z equals −80 meters.The omnidirectional localization simulation results of a sound source based on a seven-element microphone array are shown in Table 1.
From Table 1, it can be seen that although the simulation location of the sound source is different from the theoretical location, the error rate of the sound source coordinate parameters x ′ , y ′ , and z ′ is less than 0.01%, which is still within the acceptable range; compared with the theoretical angle, the deviation of simulation is also smaller, and the angle error rate is about 0.005%, which indicates that the sound source has a high geometric localization accuracy and has a very good sound source omnidirectional positioning effect.
In Experiment 2, using ( 32) and ( 33), the sound source was omnidirectionally positioned and calibrated to verify the accuracy of the calibration method.The basic parameters are set as follows: c is equal to 340 meters per second, a is equal to 2 meters, and two sets of forward data are selected (unit: meters): 60, −40, 80 and −50, 0, −60 .The set observation angle N is (unit: m): 20, 30, 40 , and the results of the omnidirectional positioning calibration of the sound source are shown in Table 2.
As can be seen from Table 2, the calibrated sound source position is close to the theoretical value.The sound source localization calibration is based on the sound source geometric localization algorithm.The calibrated sound source localization error rate is lower than 0.01%, and the error rate of horizontal declination angle h 1 and elevation angle h 2 is less

Actual Measurement Results
In the actual measurement experiment, this article uses the acoustic sensor to construct the seven-element microphone array passive acoustic data acquisition system, write the program on the Keil4 software platform, use Flymcu to receive the serial port transmission data, and record the five relative time delay values displayed in the text interface, namely, T 1 , T 2 , T 3 , T 4 , and T 5 .Use the measured relative delay values to calculate the calibrated sound source azimuth X′, Y′, Z′ , horizontal declination angle h 1 , and elevation angle h 2 .At different time points in April 2018, the indoor test sites were selected at the Nanjing University of Information Science & Technology, Pukou District, Nanjing, using a portable Bluetooth stereo to simulate sound sources.Based on the coordinate system where the seven-element microphone array is located, with the microphone M0 as the origin, the position where the three-dimensional coordinates have been measured: (−0.5 m,2 m,1.6 m), (1.8 m,3.7 m,−0.9 m), and (1.4 m,2.3 m,2.6 m), emits sound source for actual test experiment.The actual sound source omnidirectional positioning calibration experiment diagram is shown in Figure 13.The observation angle N is located at (1.6 m,2.7 m,1 m), corresponding to the three groups of sound source coordinates, and the actual array spacing a is 0.4 m, 0.6 m, and 0.8 m.The relative delay values, sound source localization, and angle calibration data are shown in Tables 3-6.
From Tables 3-6, it can be seen that when there is less noise and reverberation in the indoor environment, the measured calibration data based on the seven-element microphone array model and the error rate are in the controllable range although they are different from the theoretical data.The actual measured data is more reliable.
As can be seen from Tables 3-5, the actual measured horizontal declination angle h 1 and elevation angle h 2 are significantly different from the theoretical data.In combination with (7), the horizontal declination angle H 1 ′ is related to the sound source coordinate parameters x ′ and y ′ ; however, the value of elevation H 2 ′ is not only related to the sound source 12 Complexity coordinate parameter x ′ and the distance r ′ from the sound source to the center of the array but also related to H 1 ′ .When there is an error in the measured sound source coordinate parameters, the horizontal declination angle h 1 and the elevation angle h 2 obtained after calibration must also have errors.
Although there are some errors, compensation can be made by increasing the accuracy of the relative delay measurement.
From Table 6, it can be seen that with the increase of the actual spacing of the array elements, due to the existence of large sites and model volume factors, only the maximum array element spacing of 0.8 meters is taken, but the error rate of sound source coordinates, horizontal declination angle, and elevation angle has a tendency of getting lower and lower, the error rate of sound source coordinates is around 10%, the horizontal declination angle error is less than 5%, and the elevation angle error is less than 8%, which is close to the theoretical data.Combined with the error analysis part of the text, if we further increase the array element spacing between 2 and 10 meters, we can more fully demonstrate the superiority of the sound source omnidirectional positioning calibration method based on microphone observation angle.
In summary, in the actual measurement, the sound source omnidirectional positioning calibration method based on microphone observation angle is close to the simulation effect.Through actual measurement experiments, it is proved that the actual calibration accuracy can be more ensured after appropriately increasing the array element spacing.

Conclusions
Based on the seven-element microphone array model, the sound source geometric localization algorithm was introduced into the array model, and the sound source data are measured by array angle.This paper proposes a sound source omnidirectional positioning calibration method based on microphone observation angle.By determining the polarity of the elevation angle, omnidirectional sound source positioning can be achieved.The errors caused by the array parameters are analyzed, and some improved methods are proposed.A new coordinate system is established.From the perspective of the observation point, then the location of the sound source is calibrated.Simulation and actual measurement experiments show that the method has high positioning accuracy and can obtain accurate calibration results.
The seven-element microphone array is arranged in a cross array on a three-dimensional coordinate axis, and the array element spacing between the microphones can be adjusted to reduce the array element density in the actual arrangement area.Appropriately increasing the array element spacing can also improve the accuracy of the actual sound source omnidirectional positioning calibration.Four kinds of array models were established to solve the problem of calculating blind spots in the sound source position in the actual environment.Based on the model, the derivation of the relevant formulas for the omnidirectional positioning calibration of the sound source is carried 13 Complexity out, and the methods were applied to actual measurement experiments to make the sound source positioning calibration more comprehensive.
Of course, there is still much room for improvement in the accuracy of the calibration of sound source azimuths based on the array angle.In a real complex environment, microphones should be placed reasonably.The relative time delay value will still be affected by attenuation and outside interference during sound source propagation; there are certain errors that affect sound source positioning calibration.There are many problems involved.Assuming that the model is also ideal, how to transfer the actual indoor measurement experiment to the field for performance testing in a complex real environment requires further exploration and research in the following work.

Figure 4 :
Figure 4: Sound source S on the z-axis.

Figure 1 :
Figure1: Sound source S is not on the x-0-z plane, 0-y-z plane, and z-axis.

Figure 2 :
Figure 2: Sound source S on the x-0-z plane.

Figure 3 :
Figure 3: Sound source S on the 0-y-z plane.

Figure 5 :
Figure 5: Estimation error of the horizontal declination angle H 1 vs. array size a and elevation angle H 2.

Figure 6 :Figure 7 :
Figure 6: The relationship between the estimation error of horizontal declination angle H 1 and the time delay estimation error.

Figure 8 :
Figure 8: The relationship between the estimation error of elevation angle H 2 and the time delay estimation error.

Figure 9 :
Figure 9: The relation curve between the relative ranging error and the array size; the distance from the sound source to the center of the array.

Figure 10 :Figure 11 :
Figure 10: The relationship between microphone array size and time delay estimation error.

Figure 12 :
Figure 12: Sound source positioning calibration model based on microphone observation angle.

Table 2 :
The simulation results of the sound source omnidirectional positioning calibration method.

Table 1 :
Sound source omnidirectional localization simulation results based on seven-element microphone array.

Table 3 :
The actual calibration data on April 2nd.

Table 4 :
The actual calibration data on April 9th.

Table 5 :
The actual calibration data on April 12th.

Table 6 :
The error rate of the actual measured sound source omnidirectional positioning calibration data.