CN111679244A - Direct sound time-frequency point selection method based on plane wave relative density - Google Patents
Direct sound time-frequency point selection method based on plane wave relative density Download PDFInfo
- Publication number
- CN111679244A CN111679244A CN202010400193.5A CN202010400193A CN111679244A CN 111679244 A CN111679244 A CN 111679244A CN 202010400193 A CN202010400193 A CN 202010400193A CN 111679244 A CN111679244 A CN 111679244A
- Authority
- CN
- China
- Prior art keywords
- plane wave
- time
- density
- frequency point
- spherical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000010187 selection method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 15
- 230000001131 transforming effect Effects 0.000 claims abstract description 3
- 230000009466 transformation Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 2
- 230000006870 function Effects 0.000 description 15
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a direct sound time-frequency point selection method based on plane wave relative density, which comprises the following steps: firstly, establishing a ball array output model in a reverberation environment, and transforming the ball array output model into a time-frequency domain through short-time Fourier transform in order to apply the sparsity of a voice signal; obtaining the spherical harmonic coefficient of the plane wave density of the signal through spherical harmonic domain plane wave decomposition; utilizing inverse spherical Fourier transform to construct the plane wave density of each time-frequency point at any certain angle, wherein the higher the density is, the more likely the time-frequency point is dominated by single direct sound; after modulus is taken, traversing the angle, and calculating the ratio of the maximum plane wave density to the second maximum plane wave density as the relative density of the plane waves; and selecting the time frequency point with high relative density of the plane wave so as to obtain the final DOA estimation result of the direct sound. According to the method, a plane wave density function in a spherical harmonic domain is considered, the complex calculated amount caused by covariance matrix eigenvalue decomposition is reduced, time frequency points with high relative density are selected, and the DOA estimation precision is improved.
Description
Technical Field
The invention relates to a direct sound time-frequency point selection method based on plane wave relative density, which is applied to the technical fields of voice enhancement, video conference, robot hearing and the like.
Background
As one of important research directions of array signal processing, Direction of arrival (DOA) estimation of signals is extremely important. Many researchers are working on speech enhancement, sound scene analysis. Of these studies, more and more focus is on the study of DOA estimation in real environments. In a reverberant environment, reflected sound confuses direct sound due to multipath effects of sound propagation, so that it is difficult to accurately obtain a direct sound DOA estimation result.
In recent years, a common method for improving the robustness of reverberation is to process only signal segments containing direct sound, while rejecting those signal segments that are disturbed by reflections. For this purpose, Short-time fourier Transform (STFT) is commonly used for array signals, the signals are transformed into the time-frequency domain, and then time-frequency points dominated by direct sound are selected. With this method, the robustness to reverberation depends to a large extent on the ability to correctly select the time-frequency points dominated by the direct sound.
At present, spherical microphone arrays are widely applied to sound field analysis. The spherical array has three-dimensional symmetry, so that the sound field can be analyzed more comprehensively and the DOA of the sound source can be estimated conveniently. The spherical array acoustic field analysis is based on spherical fourier transform, i.e. the acoustic field is decomposed into orthogonal basis functions of spherical harmonics. Through spherical harmonic decomposition, the steering matrix decouples the independent portions of angle and frequency. Mohan et al propose a coherent detection method by selecting time-frequency points with an effective rank of 1, i.e. selecting time-frequency points with only a single active sound source within the time-frequency points. Since the covariance matrix is time-smoothed, when the effective rank is 1, a sound source with little coherence can be detected, and thus, the performance is good even in a low reverberation environment. However, under the condition of high reverberation, the selected time-frequency points contain direct sound and coherent reverberation signals, and the DOA estimation positioning result is influenced. Rafaely et al propose to use direct-path significance (DPD) test to select time-frequency points dominated by direct sound. The DPD test adopts frequency smoothing based on a spherical array, and the influence of coherent signals in a reverberation environment is obviously reduced. And selecting the time-frequency point corresponding to the high ratio by calculating the ratio of the maximum eigenvalue of the local covariance matrix to the second maximum eigenvalue. This test is for an array of balls, which eliminates the need to compute a focus matrix during frequency smoothing. However, the key of the frequency point when the eigenvalue ratio is selected is to select a time-frequency point when the energy of the maximum eigenvalue is relatively large, and at this time, the dominant sound source in the time-frequency point may be a reverberant sound with high energy, thereby reducing the estimation performance of the following Multiple Signal Classification (MUSIC). Madmoni et al propose a method for selecting time-frequency points based on plane wave similarity based on the MUSIC method. The method mainly comprises the step of calculating the similarity between a first eigenvector of a covariance matrix and each angle plane wave component in each time frequency point. The higher the similarity is, the higher the accuracy of the time frequency point selected by the time frequency point is, the DOA estimation performance of the final direct sound is improved, but the calculation amount is increased. In practical scene application, due to limited computing resources and strict real-time requirements, a high-complexity method for selecting the time frequency point is inevitably limited.
Disclosure of Invention
The invention aims to: aiming at the defects of the prior art, the method for selecting the direct sound time-frequency points based on the relative density of the plane waves is provided, the calculated amount of time-frequency point selection is reduced, and the estimation of the DOA of the direct sound is facilitated.
In order to achieve the above object, the idea of the present invention is:
firstly, converting a signal received based on a spherical array in a reverberation environment into a time-frequency domain; then calculating the spherical harmonic coefficient of the plane wave density function; then calculating the ratio of the maximum value of the plane wave density to the second maximum value; and finally, selecting a time frequency point corresponding to the high ratio, and calculating to obtain a DOA estimation result of the direct sound.
Firstly, establishing a signal model based on a spherical array in a reverberation environment, and converting a signal into a time-frequency domain through short-time Fourier transform in order to utilize the sparsity of a voice signal; then, performing spherical Fourier transform on the time-frequency domain signal, and performing plane wave decomposition to obtain spherical harmonic coefficients of a plane wave density function; calculating the plane wave density of each angle through inverse spherical Fourier transform, and calculating the ratio of the maximum value to the second maximum value after modulus taking, namely the plane wave relative density; and finally, selecting a time-frequency point corresponding to the high ratio in order to select a time-frequency point with dominant direct sound, and calculating to obtain a direct sound DOA estimation result.
According to the inventive concept, the technical scheme adopted by the invention is as follows:
a direct sound time-frequency point selection method based on plane wave relative density comprises the following steps:
1) establishing a spherical microphone array output model in a reverberation environment, and converting an obtained signal into a time-frequency domain through short-time Fourier transform in order to utilize the sparsity of a voice signal;
2) performing spherical Fourier transform and plane wave decomposition on the time-frequency domain signals obtained in the step 1) to obtain spherical harmonic coefficients of a plane wave density function;
3) calculating the spherical harmonic coefficient of the plane wave density function obtained in the step 2) to obtain the plane wave density of any certain angle corresponding to each time-frequency point by applying the inverse transformation of the spherical Fourier transform;
4) traversing the angle after modulus extraction of the plane wave density obtained in the step 3), and calculating the ratio of the maximum value to the second maximum value to obtain the plane wave relative density;
5) selecting the time frequency point corresponding to the relatively high density of the plane wave according to the ratio obtained in the step 4), and thus calculating to obtain a final direct sound DOA estimation result.
Compared with the prior art, the method has the following obvious and prominent substantive characteristics and remarkable advantages:
according to the method, the covariance matrix does not need to be calculated, the characteristic value decomposition is not involved, the calculated amount is reduced, the fact that the direct sound is not only high-density plane waves is considered, the maximum value of the density of the plane waves accounts for the absolute advantage is also considered, the time-frequency point corresponding to the high-density plane waves is selected in a ratio mode, and the DOA estimation precision of the final direct sound is improved.
Drawings
Fig. 1 is a flowchart of a direct sound time-frequency point selection method based on plane wave relative density according to the present invention.
Fig. 2 is a schematic diagram of a coordinate system of the spherical microphone array of the present invention.
FIG. 3 is a schematic diagram of selecting time-frequency points corresponding to a relatively high density of plane waves according to the present invention.
Detailed Description
For a better understanding of the technical solution of the present invention, the following detailed description of the preferred embodiments of the present invention is provided in conjunction with the accompanying drawings:
example one
Referring to fig. 1-3, a direct sound time-frequency point selection method based on plane wave relative density includes the following steps:
1) establishing a spherical microphone array output model in a reverberation environment, and converting an obtained signal into a time-frequency domain through short-time Fourier transform in order to utilize the sparsity of a voice signal;
2) performing spherical Fourier transform and plane wave decomposition on the time-frequency domain signals obtained in the step 1) to obtain spherical harmonic coefficients of a plane wave density function;
3) calculating the spherical harmonic coefficient of the plane wave density function obtained in the step 2) to obtain the plane wave density of any certain angle corresponding to each time-frequency point by applying the inverse transformation of the spherical Fourier transform;
4) traversing the angle after modulus extraction of the plane wave density obtained in the step 3), and calculating the ratio of the maximum value to the second maximum value to obtain the plane wave relative density;
5) selecting the time frequency point corresponding to the relatively high density of the plane wave according to the ratio obtained in the step 4), and thus calculating to obtain a final direct sound DOA estimation result.
Example two
In this embodiment, referring to fig. 1, a process of a direct sound time-frequency point selection method based on a plane wave relative density is that, in a reverberation environment, a signal received by a spherical array is transformed to a time-frequency domain through short-time fourier transform, and then spherical harmonic coefficients corresponding to a plane wave density function of the received signal are obtained through spherical fourier transform and plane wave decomposition; at the moment, the steering matrix only contains sound source angle information, a plane wave density function corresponding to any certain angle is calculated in each time frequency point by solving inverse Fourier transform, the relative density of a plane wave is calculated by considering that the direct sound is not only the high density of the plane wave, but also the plane wave with the maximum plane wave density occupying the absolute dominant plane wave, a time frequency point higher than a threshold value is selected, and the DOA estimation of the direct sound is obtained, wherein the specific implementation steps are as follows:
1) as shown in fig. 2, a spherical microphone array output model in a reverberation environment is established, and in order to utilize the sparsity of a voice signal, an obtained signal is converted into a time-frequency domain through short-time fourier transform, so that direct sound estimation is converted into a time-frequency point for selecting a direct sound leading part, which is specifically as follows:
assuming that D signal sources are incident to a uniform spherical array with an array element number of L, signals received by the spherical array are expressed in a matrix form:
p(t)=V(k,Φ)s(t)+n(t) (1)
wherein,representing the sound pressure signal received by the ball array, s (t) ═ s1(t),...,sD(t)]TIs a vector of the sound source signal,indicating a steering matrix containing DOA and frequency information, n (t) ═ n1(t),...,nL(t)]TIs a mean of 0 and a variance ofWhite additive gaussian noise.
Transform equation (2) to the time-frequency domain using short-time fourier transform (STFT), expressed as:
p(τ,ω)=V(k,Φ)s(τ,ω)+n(τ,ω) (2)
where τ is the time index and ω is the frequency index. Since the steering matrix is only related to DOA and frequency information, the steering matrix is unchanged after a short time fourier transform.
2) Performing spherical Fourier transform and plane wave decomposition on the time-frequency domain signals obtained in the step 1) to obtain spherical harmonic coefficients of a plane wave density function, wherein the specific process is as follows:
for spherical arrays, the steering matrix is decomposed into DOA and frequency components. Equation (2) can be expressed as:
p(τ,ω)=Y(Ω)B(k)YH(Φ)s(τ,ω)+n(τ,ω) (3)
wherein,
Y(Φ)=[yT(Φ1),...,yT(ΦD)]T(4)
is D × (N +1)2Phi of the sound source of dimensiond=(θd,φd) The information on the angle of the light source,each elementRepresenting spherical harmonics of order N and m, N being the maximum spherical harmonic order, and similarly, Y (Ω) is L × (N +1)2Array element of dimensional spherical arrayAngle information of (2). (N +1)2×(N+1)2The diagonal matrix of dimensions b (k) contains the radial function of the plane wave scattering from a rigid sphere.
And (4) calculating to obtain a plane wave density function through spherical harmonic domain plane wave decomposition according to the formula:
wherein, anm(τ, ω) is a length of (N +1)2The spherical harmonic coefficient of the plane wave density function of (1).
3) Calculating the spherical harmonic coefficient of the plane wave density function obtained in the step 2) by using the inverse transformation of the spherical Fourier transform to obtain each angle theta corresponding to each time-frequency pointj(J1.. J.) the corresponding plane wave density is as follows:
a(τ,ω,Θj)=yT(Θ)a(τ,ω) (6)
4) traversing the angle after modulus of the plane wave density obtained in the step 3), and calculating the ratio of the maximum value to the second maximum value to obtain the plane wave relative density, wherein the specific process is as follows:
DEN(τ,ω,Θj)=|a(τ,ω,Θj)| (7)
where DEN (τ, ω) represents the mode of the calculated plane wave density at (τ, ω). The larger the value of DEN (τ, ω), the more likely it is that the direct sound is dominant. However, if (τ, ω) is dominant in reverberant sound, the DEN (τ, ω) also has a maximum value, which is not sufficiently prominent.
In order to select the time frequency point with the relatively prominent maximum value, a plane wave density ratio is set, and the expression is as follows:
wherein RA (τ, ω) represents the plane wave density ratio at this time-frequency point (τ, ω), Θj1And Θj2Corresponding to the angles at which the DEN gets the maximum and second largest value, respectively.
5) Selecting the time frequency point corresponding to the relatively high density of the plane wave from the ratio obtained in the step 4), thereby obtaining a final direct sound DOA estimation result by calculation, wherein the specific process is as follows:
D*={(τ,ω):RA(τ,ω)>TH} (9)
wherein D is*The selected time frequency point set is shown, and TH is a threshold value. Calculating the final DOA estimated value in the selected time frequency pointsExpressed as:
in summary, the method for selecting the direct sound time-frequency point based on the relative density of the plane wave comprises the following steps: firstly, establishing a ball array output model in a reverberation environment, and transforming the ball array output model into a time-frequency domain through short-time Fourier transform in order to apply the sparsity of a voice signal; obtaining the spherical harmonic coefficient of the plane wave density of the signal through spherical harmonic domain plane wave decomposition; utilizing inverse spherical Fourier transform to construct the plane wave density of each time-frequency point at any certain angle, wherein the higher the density is, the more likely the time-frequency point is dominated by single direct sound; after modulus is taken, traversing the angle, and calculating the ratio of the maximum plane wave density to the second maximum plane wave density as the relative density of the plane waves; and selecting the time frequency point with high relative density of the plane wave so as to obtain the final DOA estimation result of the direct sound. The method considers the plane wave density function in the spherical harmonic domain, reduces the complex calculated amount caused by covariance matrix eigenvalue decomposition, considers that direct sound is not only the plane wave high density but also the maximum value of the plane wave density occupies the absolute advantage, provides a method for the plane wave relative density, selects the time frequency point with high relative density, and improves the DOA estimation precision.
The embodiments of the present invention have been described above with reference to the accompanying drawings, but the present invention is not limited to the above embodiments, and various changes and modifications can be made according to the purpose of the invention, and all changes, modifications, substitutions, combinations or simplifications made according to the spirit and principle of the technical solution of the present invention shall be equivalent substitution ways, so long as the invention meets the purpose of the present invention, and the technical principle and inventive concept of the method for selecting the time-frequency point of direct sound based on the relative density of plane waves shall fall within the protection scope of the present invention.
Claims (6)
1. A direct sound time-frequency point selection method based on plane wave relative density is characterized by comprising the following steps:
1) establishing a spherical microphone array output model in a reverberation environment, and converting an obtained signal into a time-frequency domain through short-time Fourier transform in order to utilize the sparsity of a voice signal;
2) performing spherical Fourier transform and plane wave decomposition on the time-frequency domain signals obtained in the step 1) to obtain spherical harmonic coefficients of a plane wave density function;
3) calculating the spherical harmonic coefficient of the plane wave density function obtained in the step 2) to obtain the plane wave density of any certain angle corresponding to each time-frequency point by applying the inverse transformation of the spherical Fourier transform;
4) traversing the angle after modulus extraction of the plane wave density obtained in the step 3), and calculating the ratio of the maximum value to the second maximum value to obtain the plane wave relative density;
5) selecting the time frequency point corresponding to the relatively high density of the plane wave according to the ratio obtained in the step 4), and thus calculating to obtain a final direct sound DOA estimation result.
2. The method for selecting the direct sound time-frequency point based on the relative density of the plane wave according to claim 1, wherein the spherical microphone array output model in the reverberation environment is established in the step 1), and specifically comprises the following steps:
assuming that D signal sources are incident to a uniform spherical array with an array element number of L, signals received by the spherical array are expressed in a matrix form:
p(t)=V(k,Φ)s(t)+n(t) (1)
wherein,representing the sound pressure signal received by the ball array, s (t) ═ s1(t),...,sD(t)]TIs a vector of the sound source signal,indicating a steering matrix containing DOA and frequency information, n (t) ═ n1(t),...,nL(t)]TIs a mean of 0 and a variance ofAdditive white gaussian noise of (1);
transforming equation (2) to the time-frequency domain using a short-time fourier transform, expressed as:
p(τ,ω)=V(k,Φ)s(τ,ω)+n(τ,ω) (2)
wherein τ is a time index and ω is a frequency index; since the steering matrix is only related to DOA and frequency information, the steering matrix is unchanged after a short time fourier transform.
3. The method for selecting the direct sound time-frequency point based on the relative density of the plane wave as claimed in claim 1, wherein the spherical harmonic coefficient of the plane wave density function is obtained by performing spherical Fourier transform and plane wave decomposition on the time-frequency domain signal obtained in step 1), and the specific operation steps are as follows:
for a spherical array, the steering matrix is decomposed into DOA and frequency components; equation (2) can be expressed as:
p(τ,ω)=Y(Ω)B(k)YH(Φ)s(τ,ω)+n(τ,ω) (3)
wherein,
Y(Φ)=[yT(Φ1),...,yT(ΦD)]T(4)
is D × (N +1)2Phi of the sound source of dimensiond=(θd,φd) The information on the angle of the light source,each elementRepresenting spherical harmonics of order N and m, N being the maximum spherical harmonic order, and similarly, Y (Ω) being L × (N +1)2Array element of dimensional spherical arrayAngle information of (a); (N +1)2×(N+1)2The diagonal matrix of dimensions B (k) contains the radial function of the plane wave scattering from a rigid sphere;
and (4) calculating to obtain a plane wave density function through spherical harmonic domain plane wave decomposition according to the formula:
wherein, anm(τ, ω) is a length of (N +1)2The spherical harmonic coefficient of the plane wave density function of (1).
4. The method for selecting direct sound time-frequency points based on plane wave relative density according to claim 1, wherein the spherical harmonic coefficients of the plane wave density function obtained in the step 2) are subjected to inverse spherical Fourier transform to calculate and obtain each angle theta corresponding to each time-frequency pointj(J1.. J.) the specific operation procedure for the plane wave density is as follows:
a(τ,ω,Θj)=yT(Θ)a(τ,ω) (6)。
5. The method for selecting the direct sound time-frequency point based on the relative density of the plane wave as claimed in claim 1, wherein the relative density of the plane wave is obtained by taking the modulus of the density of the plane wave obtained in step 3), traversing the angle, and calculating the ratio of the maximum value to the second maximum value, and the operation steps are as follows:
DEN(τ,ω,Θj)=|a(τ,ω,Θj)| (7)
where DEN (τ, ω) represents the mode of the calculated plane wave density at (τ, ω); the larger the value of DEN (τ, ω), the more likely (τ, ω) is that the direct sound is dominant; if (τ, ω) is dominant in reverberant sound, there is also a maximum in DEN (τ, ω), which is not sufficiently prominent;
in order to select the time frequency point with the relatively prominent maximum value, a plane wave density ratio is set, and the expression is as follows:
6. The method for selecting the direct sound time-frequency point based on the relative density of the plane wave according to claim 1, wherein the ratio obtained in the step 4) is used for selecting the time-frequency point corresponding to the relatively high density of the plane wave, so as to calculate and obtain the final direct sound DOA estimation result, and the method comprises the following operation steps:
D*={(τ,ω):RA(τ,ω)>TH} (9)
wherein D is*Representing the selected time frequency point set, TH being a threshold value; calculating the final DOA estimated value in the selected time frequency pointsExpressed as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010400193.5A CN111679244B (en) | 2020-05-13 | 2020-05-13 | Direct sound time-frequency point selection method based on plane wave relative density |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010400193.5A CN111679244B (en) | 2020-05-13 | 2020-05-13 | Direct sound time-frequency point selection method based on plane wave relative density |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111679244A true CN111679244A (en) | 2020-09-18 |
CN111679244B CN111679244B (en) | 2023-11-17 |
Family
ID=72434188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010400193.5A Active CN111679244B (en) | 2020-05-13 | 2020-05-13 | Direct sound time-frequency point selection method based on plane wave relative density |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111679244B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687304A (en) * | 2021-07-07 | 2021-11-23 | 浙江大华技术股份有限公司 | Direct sound detection method, system and computer readable storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592628A (en) * | 2013-11-12 | 2014-02-19 | 上海大学 | Multi-sound-source positioning method based on formation of real value weight beam in spherical harmonic domain |
CN106066490A (en) * | 2016-05-25 | 2016-11-02 | 中国石油大学(北京) | Prestack inversion of Density method and device based on spherical wave |
JP2017130899A (en) * | 2016-01-22 | 2017-07-27 | 日本電信電話株式会社 | Sound field estimation device, method therefor and program |
-
2020
- 2020-05-13 CN CN202010400193.5A patent/CN111679244B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592628A (en) * | 2013-11-12 | 2014-02-19 | 上海大学 | Multi-sound-source positioning method based on formation of real value weight beam in spherical harmonic domain |
JP2017130899A (en) * | 2016-01-22 | 2017-07-27 | 日本電信電話株式会社 | Sound field estimation device, method therefor and program |
CN106066490A (en) * | 2016-05-25 | 2016-11-02 | 中国石油大学(北京) | Prestack inversion of Density method and device based on spherical wave |
Non-Patent Citations (3)
Title |
---|
周晓凤 等: "基于密度聚类的直达声波达方向估计", 《工业控制计算机》 * |
孙长伟 等: "混响环境下改进的球谐波域L1-SVD声源定位算法", 《武汉大学学报( 理学版)》 * |
黄青华 等: "基于球谐域MVDR波束形成技术的远场多声源定位", 《新型工业化》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687304A (en) * | 2021-07-07 | 2021-11-23 | 浙江大华技术股份有限公司 | Direct sound detection method, system and computer readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN111679244B (en) | 2023-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109597022B (en) | Method, device and equipment for calculating azimuth angle of sound source and positioning target audio | |
CN109087664B (en) | Speech enhancement method | |
Herzog et al. | Eigenbeam-ESPRIT for DOA-vector estimation | |
Moore et al. | Direction of arrival estimation using pseudo-intensity vectors with direct-path dominance test | |
CN111044973A (en) | MVDR target sound source directional pickup method for microphone matrix | |
Ren et al. | A novel multiple sparse source localization using triangular pyramid microphone array | |
Xia et al. | Noise reduction method for acoustic sensor arrays in underwater noise | |
CN111679244B (en) | Direct sound time-frequency point selection method based on plane wave relative density | |
CN113593596B (en) | Robust self-adaptive beam forming directional pickup method based on subarray division | |
Dwivedi et al. | Joint doa estimation in spherical harmonics domain using low complexity cnn | |
Dai et al. | Multiple speech sources localization in room reverberant environment using spherical harmonic sparse Bayesian learning | |
Gao et al. | An order-aware scheme for robust direction of arrival estimation in the spherical harmonic domain | |
CN105138743B (en) | A kind of coherent signal source direction estimation method based on single vector hydrophone | |
Çöteli et al. | Multiple sound source localization with rigid spherical microphone arrays via residual energy test | |
SongGong et al. | Robust indoor speaker localization in the circular harmonic domain | |
CN108269581B (en) | Double-microphone time delay difference estimation method based on frequency domain coherent function | |
Dwivedi et al. | Hybrid sh-cnn-mp approach for super resolution doa estimation | |
CN110824484B (en) | Array element position estimation method based on constant modulus algorithm | |
Pan et al. | Multiple spherical arrays design for acoustic source localization | |
Chen et al. | Multi-channel end-to-end neural network for speech enhancement, source localization, and voice activity detection | |
Huang et al. | Two‐dimensional DOA estimation using one‐dimensional search for spherical arrays | |
Huang et al. | Source localization with minimum variance distortionless response for spherical microphone arrays | |
Liu et al. | Sound source localization and speech enhancement algorithm based on fixed beamforming | |
Bin et al. | The combination of spectrum subtraction and cross-power spectrum phase method for time delay estimation | |
Xiao et al. | Acoustic Beamforming via Interference-Plus-Noise Covariance Matrix Construction for Interferences and Noise Attenuation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |