CN103135092A - Micro aperture acoustic array moving objective orientation method - Google Patents
Micro aperture acoustic array moving objective orientation method Download PDFInfo
- Publication number
- CN103135092A CN103135092A CN2013100467438A CN201310046743A CN103135092A CN 103135092 A CN103135092 A CN 103135092A CN 2013100467438 A CN2013100467438 A CN 2013100467438A CN 201310046743 A CN201310046743 A CN 201310046743A CN 103135092 A CN103135092 A CN 103135092A
- Authority
- CN
- China
- Prior art keywords
- array
- moving target
- signal
- orientation method
- pattern matrix
- 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.)
- Pending
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention relates to a micro aperture acoustic array moving objective orientation method. The micro aperture acoustic array moving objective orientation method comprises the following steps: selecting a mike array, ascertaining a signal broadband B= fH-fL of a moving objective, wherein, fH is a highest value of the signal broadband, and fL is the lowest value of the signal broadband, confirming a stream pattern matrix, enabling the frequency coefficient of the stream pattern matrix to use a mathematical relation between the stream pattern matrix and the mike array, through an angle searching method, and acquiring an azimuth angle anticipation of the moving object. The micro aperture acoustic array moving objective orientation method utilizes an orthogonality relation between the stream pattern matrix and a mike signal, starting to search the azimuth angle, and so that effective azimuth angle anticipation of the moving objective is acquired. At the premise of guaranteeing not to change orientation precision, the micro aperture acoustic array moving objective orientation method greatly reduces required computation complexity of anticipating the azimuth angle of the moving objective.
Description
Technical field
The present invention relates to the array signal process technique field, particularly relate to a kind of micropore diameter acoustic array moving target orientation method.
Background technology
The goal orientation method is in field extensive application such as radar, sonar, ground Unattended Ground Sensor networks.The acoustical signal that the main movement targets such as wheeled vehicle, creeper truck, low flyer produce is the exemplary wideband random signal, and the orientation method for wideband random signal mainly comprises at present: 1. based on the orientation method 2. of irrelevant signal based on the orientation method 3. of the coherent signal maximum likelihood method for broadband sound source.Wherein broadband signal need to be decomposed narrow band data on nonoverlapping bands based on the orientation method of irrelevant signal, each arrowband is processed respectively, estimate thereby obtain the position angle, the method thinking is simple, but precision is relatively poor, and computational complexity is large; Orientation method based on coherent signal need to focus on the signal space on nonoverlapping Frequency point in frequency band the reference frequency point, obtain the data covariance of single-frequency point after focusing, the method that the recycling narrow band signal is processed is carried out the position angle and is estimated, the method or need to pre-estimate, perhaps need complicated focus method, method complexity and computational complexity are large; Maximum likelihood method operational precision for broadband sound source is high, applied widely, but computational complexity is huge, is widely used thereby limited to a great extent it.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of micropore diameter acoustic array moving target orientation method, can guarantee the constant prerequisite decline low computational complexity of orientation accuracy.
The technical solution adopted for the present invention to solve the technical problems is: a kind of micropore diameter acoustic array moving target orientation method is provided, comprises the following steps:
(1) select microphone array, determine the signal bandwidth B=f of moving target
H-f
L, wherein, f
HBe the mxm. of signal bandwidth, f
LMinimum for signal bandwidth;
(2) determine the flow pattern matrix, make the frequency parameter of flow pattern matrix
(3) utilize mathematical relation between flow pattern matrix and microphone array column signal, the position angle that obtains moving target by the angle searching mode is estimated.
Satisfy between the signal bandwidth B of the maximum diameter of hole D of microphone array and moving target in described step (1)
Wherein, V
SoundBe the velocity of sound.
The normalization frequency spectrum that the frequency parameter f of the flow pattern matrix in described step (2) satisfies the signal broadband makes its probability density function standard deviation minimum.
Adopt multi-signal sorting algorithm in described step (3), beamforming algorithm, maximum likelihood algorithm or maximum entropy algorithm realize that the angle searching mode obtains the position angle estimation of moving target.
Microphone array in described step (1) is classified linear array, circle battle array, planar array or triangle battle array as.
Beneficial effect
Owing to having adopted above-mentioned technical scheme, the present invention compared with prior art, have following advantage and good effect: the present invention utilizes the orthogonality relation of flow pattern matrix and microphone signal, carries out the position angle search, and effective position angle that can obtain moving target is estimated.Orientation method of the present invention both need not that each Frequency point was carried out separately the position angle to be estimated, need not again frequency focusing is carried out the position angle estimation again after reference frequency, thereby greatly reduces computational complexity.The present invention is guaranteeing under the constant prerequisite of orientation accuracy, has reduced greatly the needed computational complexity of moving target position angle estimation.
Embodiment
Below in conjunction with specific embodiment, further set forth the present invention.Should be understood that these embodiment only to be used for explanation the present invention and be not used in and limit the scope of the invention.Should be understood that in addition those skilled in the art can make various changes or modifications the present invention after the content of having read the present invention's instruction, these equivalent form of values fall within the application's appended claims limited range equally.
Embodiments of the present invention relate to a kind of micropore diameter acoustic array moving target orientation method, comprise the following steps:
(1) select microphone array, determine the signal bandwidth B=f of moving target
H-f
L, wherein, f
HBe the mxm. of signal bandwidth, f
LMinimum for signal bandwidth.Wherein, satisfy between the signal bandwidth B of the maximum diameter of hole D of microphone array and moving target
Wherein, V
SoundBe the velocity of sound.Guarantee
Can be so that the unchangeability of the mathematical relation of described flow pattern matrix A and microphone signal, and
The position angle estimated result of the less gained of value more accurate.Microphone array can be but be not limited to as Types Below: linear array, circle battle array, planar array, triangle battle array.
(2) determine the flow pattern matrix, make the frequency parameter of flow pattern matrix
Wherein, the normalization frequency spectrum that the frequency parameter f of flow pattern matrix satisfies the signal broadband makes its probability density function standard deviation hour, and the position angle estimated accuracy is the highest, and the method directional effect is best.Flow pattern matrix A=[a
1(ω
0) a
2(ω
0) ... a
n(ω
0)], wherein,
i=1,2,...,n
。In formula, ω
0=2 * π * f is reference point: τ take the center of circle for uniform circular array
ki=r/V
Sound(cos (2 (k-1)-θ
i) cos φ
i), k=1,2 ..., M, wherein, θ
iAnd φ
iRepresent respectively position angle and the angle of pitch, r is radius of a circle.
(3) utilize mathematical relation between flow pattern matrix and microphone array column signal, the position angle that obtains moving target by the angle searching mode is estimated.Wherein, when selecting the Multiple Signal Classification method, the mathematical relation between flow pattern matrix and microphone array column signal is the steering vector a in the flow pattern matrix
i(ω
0) with the noise subspace quadrature of sensor array signal.But this step can adopt and be not limited to the position angle that following method realizes that the angle searching mode obtains moving target and estimate, multi-signal sorting algorithm, beamforming algorithm, maximum likelihood algorithm and maximum entropy algorithm.
Need to prove, for different algorithms, above-mentioned mathematical relation is not identical yet, when adopting the multi-signal sorting algorithm in following specific embodiment, is orthogonality relation between flow pattern matrix and microphone array column signal.
The below further illustrates the present invention with a specific embodiment.
1. select quaternary circle battle array microphone array, circle battle array aperture D is 0.04m, THE VELOCITY OF SOUND IN AIR V
SoundBe 340m/s.
2. the signal bandwidth B of main target wheeled vehicle, creeper truck and low flyer can be limited to 0～400Hz, at this moment,
Satisfy the condition less than 1/20.This moment, the frequency parameter f of flow pattern matrix A can be chosen as 200Hz.
3. after definite flow pattern matrix A, the MUSIC spectrum of the above-mentioned echo signal of search can realize the effective orientation to target.It is worth mentioning that, in above-described embodiment, microphone array can also be linear array, planar array and triangle battle array, and realization can also be beamforming algorithm, maximum likelihood algorithm and maximum entropy algorithm to the method for effective orientation of target.
Use Matlab software to carry out emulation to said method, contrast the working time of other broadband goal orientation methods, result is as shown in the table:
Annotate: TCT is a kind of typical coherent signal processing method
Be not difficult to find, the present invention utilizes the orthogonality relation of flow pattern matrix and microphone signal, carries out the position angle search, and effective position angle that can obtain moving target is estimated.Orientation method of the present invention both need not that each Frequency point was carried out separately the position angle to be estimated, need not again frequency focusing is carried out the position angle estimation again after reference frequency, thereby greatly reduces computational complexity.The present invention is guaranteeing under the constant prerequisite of orientation accuracy, has reduced greatly the needed computational complexity of moving target position angle estimation.
Claims (5)
1. a micropore diameter acoustic array moving target orientation method, is characterized in that, comprises the following steps:
(1) select microphone array, determine the signal bandwidth B=f of moving target
H-f
L, wherein, f
HBe the mxm. of signal bandwidth, f
LMinimum for signal bandwidth;
(2) determine the flow pattern matrix, make the frequency parameter of flow pattern matrix
(3) utilize mathematical relation between flow pattern matrix and microphone array column signal, the position angle that obtains moving target by the angle searching mode is estimated.
2. micropore diameter acoustic array moving target orientation method according to claim 1, is characterized in that, satisfies between the signal bandwidth B of the maximum diameter of hole D of microphone array and moving target in described step (1)
Wherein, V
SoundBe the velocity of sound.
3. micropore diameter acoustic array moving target orientation method according to claim 1, is characterized in that, the normalization frequency spectrum that the frequency parameter f of the flow pattern matrix in described step (2) satisfies the signal broadband makes its probability density function standard deviation minimum.
4. micropore diameter acoustic array moving target orientation method according to claim 1, it is characterized in that, adopt multi-signal sorting algorithm in described step (3), beamforming algorithm, maximum likelihood algorithm or maximum entropy algorithm realize that the angle searching mode obtains the position angle estimation of moving target.
5. micropore diameter acoustic array moving target orientation method according to claim 1, is characterized in that, the microphone array in described step (1) is classified linear array, circle battle array, planar array or triangle battle array as.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013100467438A CN103135092A (en) | 2013-02-05 | 2013-02-05 | Micro aperture acoustic array moving objective orientation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013100467438A CN103135092A (en) | 2013-02-05 | 2013-02-05 | Micro aperture acoustic array moving objective orientation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103135092A true CN103135092A (en) | 2013-06-05 |
Family
ID=48495175
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2013100467438A Pending CN103135092A (en) | 2013-02-05 | 2013-02-05 | Micro aperture acoustic array moving objective orientation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103135092A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103994820A (en) * | 2014-04-21 | 2014-08-20 | 中国科学院上海微系统与信息技术研究所 | Moving target identification method based on micro-aperture microphone array |
| CN104035065A (en) * | 2014-06-23 | 2014-09-10 | 河北工业大学 | Sound source orienting device on basis of active rotation and method for applying sound source orienting device |
| CN106526533A (en) * | 2016-11-14 | 2017-03-22 | 中国科学院上海微系统与信息技术研究所 | Microporous MEMS acoustic array sensor and use method thereof |
| CN108683972A (en) * | 2014-10-20 | 2018-10-19 | 索尼公司 | sound processing system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101184328A (en) * | 2007-04-18 | 2008-05-21 | 中国科学院上海微系统与信息技术研究所 | Wireless sensing network node locating method |
| CN101349742A (en) * | 2008-08-29 | 2009-01-21 | 西安电子科技大学 | Method for optimizing space between broad band phased array elements and measuring frequency and direction of frequency domain multiple targets |
| CN102435299A (en) * | 2011-08-12 | 2012-05-02 | 中国科学院上海微系统与信息技术研究所 | Mini-type vibration sensor array and application thereof to ground target tracking |
-
2013
- 2013-02-05 CN CN2013100467438A patent/CN103135092A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101184328A (en) * | 2007-04-18 | 2008-05-21 | 中国科学院上海微系统与信息技术研究所 | Wireless sensing network node locating method |
| CN101349742A (en) * | 2008-08-29 | 2009-01-21 | 西安电子科技大学 | Method for optimizing space between broad band phased array elements and measuring frequency and direction of frequency domain multiple targets |
| CN102435299A (en) * | 2011-08-12 | 2012-05-02 | 中国科学院上海微系统与信息技术研究所 | Mini-type vibration sensor array and application thereof to ground target tracking |
Non-Patent Citations (2)
| Title |
|---|
| 刘云: "目标定向及多波束实现", 《中国优秀博硕士学位论文全文数据库（硕士） 信息科技辑》 * |
| 蒋敏 等: "多输入多输出阵列目标方位估计方法", 《系统工程与电子技术》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103994820A (en) * | 2014-04-21 | 2014-08-20 | 中国科学院上海微系统与信息技术研究所 | Moving target identification method based on micro-aperture microphone array |
| CN103994820B (en) * | 2014-04-21 | 2016-05-11 | 中国科学院上海微系统与信息技术研究所 | A kind of moving target recognition methods based on micropore diameter microphone array |
| CN104035065A (en) * | 2014-06-23 | 2014-09-10 | 河北工业大学 | Sound source orienting device on basis of active rotation and method for applying sound source orienting device |
| CN108683972A (en) * | 2014-10-20 | 2018-10-19 | 索尼公司 | sound processing system |
| US11172292B2 (en) | 2014-10-20 | 2021-11-09 | Sony Corporation | Voice processing system |
| CN108683972B (en) * | 2014-10-20 | 2022-08-09 | 索尼公司 | Sound processing system |
| CN106526533A (en) * | 2016-11-14 | 2017-03-22 | 中国科学院上海微系统与信息技术研究所 | Microporous MEMS acoustic array sensor and use method thereof |
| CN106526533B (en) * | 2016-11-14 | 2019-03-05 | 中国科学院上海微系统与信息技术研究所 | A kind of micropore diameter MEMS acoustic matrix sensor and its application method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109997375B (en) | Concentric differential microphone array and associated beamforming | |
| CN102411136B (en) | Phase interferometer direction finding method for ambiguity resolution by extension baselines | |
| CN111123192B (en) | Two-dimensional DOA positioning method based on circular array and virtual extension | |
| CN105912791B (en) | DOA estimation method based on local search in virtual relatively prime array | |
| CN101995574B (en) | Near field focusing beam forming positioning method | |
| CN104360310B (en) | A kind of multiple target near-field sources localization method and apparatus | |
| CN104777450B (en) | A kind of two-stage MUSIC microphone array direction-finding method | |
| CN102147458B (en) | Method and device for estimating direction of arrival (DOA) of broadband sound source | |
| CN103308889A (en) | Passive sound source two-dimensional DOA (direction of arrival) estimation method under complex environment | |
| CN104166120B (en) | A kind of acoustic vector justifies battle array robust wideband MVDR direction estimation methods | |
| CN104730513A (en) | Multistage sub-array focusing MVDR wave beam forming method | |
| CN103018730A (en) | Distributed sub-array wave arrival direction estimation method | |
| CN103135092A (en) | Micro aperture acoustic array moving objective orientation method | |
| US10670694B1 (en) | System and method for hybrid-weighted, multi-gridded acoustic source location | |
| CN103353588B (en) | Two-dimensional DOA (direction of arrival) angle estimation method based on antenna uniform planar array | |
| CN103969630B (en) | A kind of robust wideband Beamforming Method constant based on frequency response | |
| CN106093866A (en) | A kind of sound localization method being applicable to hollow ball array | |
| CN106249244A (en) | A kind of acoustic vector circle battle array wide-band coherent signal source direction estimation method | |
| CN104898086A (en) | Sound intensity estimation sound source orientation method applicable for minitype microphone array | |
| CN101252382A (en) | Wide frequency range signal polarizing and DOA estimating method and apparatus | |
| CN105572629B (en) | A kind of two-dimentional direction-finding method of low computational complexity suitable for General Cell structure | |
| CN103901421B (en) | Underwater sound array SMI-MVDR Estimation of Spatial Spectrum method based on diagonal angle off-load | |
| Zhang et al. | Acoustic source localization via subspace based method using small aperture MEMS arrays | |
| CN105572631A (en) | Maximum likelihood target DOA estimation method based on multi-wave potential combined treatment | |
| CN112305497B (en) | DOA estimation direction finding fuzzy elimination method for near field microphone array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130605 |