CN101403791B - Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof - Google Patents
Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof Download PDFInfo
- Publication number
- CN101403791B CN101403791B CN2008101374641A CN200810137464A CN101403791B CN 101403791 B CN101403791 B CN 101403791B CN 2008101374641 A CN2008101374641 A CN 2008101374641A CN 200810137464 A CN200810137464 A CN 200810137464A CN 101403791 B CN101403791 B CN 101403791B
- Authority
- CN
- China
- Prior art keywords
- dsp
- data
- spectrum
- bus
- finding
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a super-resolution direction-finding device used for rapidly and real time estimating a spatial spectrum and a direction-finding method thereof. The super-resolution direction-finding device mainly consists of a principal-subordinate DSP parallel signal processing module formed by a group of DSP, a logic control module formed by FPGA and a communication interface module, wherein, the communication interface module is used for receiving sampling data of all channels from a digital receiver and reporting calculated angles to a host computer; a principal DSP and a subordinate DSP in the parallel signal processing module are connected in a bus-sharing close coupling mode, and a data bus is utilized to transmit the middle data among DSP; a DSP bus and the logic control module are mutually connected by the data bus, an address bus and a control bus. The invention can be used for rapidly and precisely finding direction for a plurality of radiant sources and further precisely positioning the radiant sources, and also can be used for tracking and positioning moving radiant sources according to requirements, or can be arranged on a moving device to track and position the radiant sources, and simultaneously, the invention also can be used for the boat positioning, the maritime search and rescue, the radar exploration in a civil aviation airport, the direction finding of undesired signals, and the like.
Description
(1) technical field
What the present invention relates to is a kind of direction-finding device, the invention still further relates to a kind of direction-finding method.The apparatus and method that specifically a kind of wave electromagnetic radiation source is surveyed.
(2) background technology
The existing device that radiation source is carried out direction finding mainly adopts the phase-interferometer principle, uses than phase system and carries out direction finding.Shortcoming than phase system is to carry out direction finding to individual signals, and direction finding precision is not high enough, also occurs many-valued fuzziness easily at high band.Recently the Estimation of Spatial Spectrum method that proposes can break through the Rayleigh restriction, realizes high-precision super-resolution direction finding.But the calculated amount of this method is very big, and causes it to be difficult to be applied to the occasion that real-time requires.The speed that how to improve Estimation of Spatial Spectrum direction finding algorithm is the gordian technique that this technology is applied to engineering reality to realize the real-time of direction finding.
(3) summary of the invention
The object of the present invention is to provide a kind of fast in real time, fast real-time space spectrum estimation ultra-resolution direction-finding device that direction finding precision is high.The present invention also aims to provide a kind of fast real-time space spectrum estimation ultra-resolution direction-finding method.
The object of the present invention is achieved like this:
The principal and subordinate DSP parallel signal processing module that fast real-time space spectrum estimation ultra-resolution direction-finding device of the present invention mainly is made of one group of DSP, the Logic control module and communication interface modules three parts that are made of FPGA are formed; Communication interface modules is reported and submitted to host computer from a sampled data and an angle that calculates that digital receiver receives each passage; The close coupled system that principal and subordinate DSP in the parallel signal processing module adopts bus to share connects, and utilizes data bus transmission intermediate data between the DSP; Interconnected between dsp bus and the Logic control module by data, address and control bus.
Direction-finding method based on fast real-time space spectrum estimation ultra-resolution direction-finding device of the present invention is: at first receive the data of digital receiver sampling, data are carried out pre-service, obtain correlation matrix then; Correlation matrix is carried out feature decomposition, draw the noise characteristic subspace, and utilize eigenwert to estimate the number of radiation source; The spatial spectrum value that calculates the given angle zone with the direction steering vector and the noise characteristic subspace of different incidence angles degree correspondence then; Again the spectrum value is carried out spectrum peak search, the position correspondence at the spectrum peak that obtains the true incident direction of radiation source.
In order to improve arithmetic speed, guarantee real-time, direction-finding device of the present invention and method have adopted following technology: the one, the close coupled system that hardware configuration adopts the multi-disc dsp bus to share, multi-disc DSP parallel processing.Utilize data bus transmission intermediate data between the DSP, realize the data transmission zero-waiting, guarantee the high-speed real-time computing.The 2nd, spectrum peak search adopts the mode of coarse search and smart search combination.Coarse search is found out peak point with big step-length in region-wide, smart search with small step progress line search, is found out accurate peak point near the peak point that coarse search is found out.When having guaranteed precision, improved real-time like this.The 3rd, for the emitter Signals of given antenna array structure and certain frequency, the required steering vector of computer memory spectrum value adopts lookup table mode.DSP computer memory spectrum need be used the steering vector of noise subspace and each angle.After each Data Update, noise subspace can change, but steering vector is constant.Because it is bigger to calculate the operand of steering vector, the required time is very long, so behind the system power-on reset, DSP at first calculate all big step length searchings the steering vector of all angles of process, make the steering vector form of a two dimension according to the angle of pitch and position angle, leave on the sheet of DSP in the internal memory.Like this, when the pairing spatial spectrum value in a certain angle of pitch of each calculating and position angle, parse the relative address of luffing angle dimension in the steering vector form earlier according to the angle of pitch, parse the relative address of orientation angles dimension in the steering vector form then according to the position angle, again with two relative address additions, just obtain the specific address of the pairing steering vector of this position angle and the angle of pitch in the steering vector form, call steering vector by this specific address table look-at and come computer memory spectrum value, and need not calculate steering vector in real time, save a lot of times, greatly improved real-time.
The present invention can carry out accurately direction finding fast to a plurality of radiation sources simultaneously, and then it is accurately located, and can also follow the tracks of, locate mobile radiation source as required, perhaps is contained on the mobile device radiation source is followed the tracks of, located.Also can be used for the location and the maritime search and rescue of ship simultaneously, the detection of Civil Aviation Airport radar and the direction finding of undesired signal etc.
(4) description of drawings
Fig. 1 is the theory diagram of fast real-time space spectrum estimation ultra-resolution direction-finding device of the present invention;
Fig. 2 is the process flow diagram of fast real-time space spectrum estimation ultra-resolution direction-finding method of the present invention.
(5) embodiment
For example the present invention is done description in more detail below in conjunction with accompanying drawing:
In conjunction with Fig. 1, the Logic control module and communication interface modules three parts that say principal and subordinate DSP parallel signal processing module that fast real-time space spectrum estimation ultra-resolution direction-finding device of the present invention mainly is made of four DSP, are made of FPGA are formed.Communication interface modules is finished from digital receiver and is received the sampled data of each passage and the angle that calculates is reported and submitted to host computer.Principal and subordinate DSP parallel signal processing module is the nucleus module of signal processor, finishes whole signal processing algorithms high-speed real-time.Algorithm mainly is an orthogonality of utilizing signal subspace and noise subspace in the sampled data, forms sharp-pointed needle-like spatial spectrum in given angular regions, and the pairing angle in spectrum peak is exactly the angle of radiation source incident.Find out the position at the spectrum peak in the zone then by spectrum peak search, also just can accurately estimate the direction of radiation source.Logic control module is finished look-at-me generation, marking signal generation, data buffering and functions such as conversion, communication interface logic control.
In conjunction with Fig. 2, fast real-time space spectrum estimation ultra-resolution direction-finding method of the present invention is: after system powered on, processor at first carried out initialization.At once calculate the steering vector form after initialization finishes, the result is left in internal memory on the sheet of DSP.Calculated after the steering vector, processor is waited for the initiation command that digital receiver sends after sampling finishes.In case initiation command is arranged, processor master DSP receives the data of digital receiver sampling and begins to enter handling procedure.These data are by 5 passages, 100 matrixes of 5 * 100 that the sampling snap constitutes.Processor at first carries out relevant treatment to the data matrix and obtains 5 * 5 correlation matrix.Correlation matrix is carried out feature decomposition, obtain eigenwert and proper vector.Utilize eigenwert to estimate the number of radiation source by the information theory criterion.Utilize information source number that proper vector is divided, obtain signal characteristic vector and noise feature vector respectively.Noise feature vector is also with regard to corresponding noise characteristic subspace, and main DSP gives each sheet from DSP by data bus with noise feature vector.Principal and subordinate DSP begins computer memory spectrum value in angular regions separately then.Calculating the required steering vector of using of spectrum value directly calls in internal memory by lookup table mode.After the spatial spectrum value in whole zone was all calculated and finished, each DSP composed the peak coarse search with big step-length in angular regions separately.After search finishes, send spectrum peak position and spectrum peak that coarse search obtains back to main DSP again from DSP.Main DSP compares all coarse search spectrum peaks, determines the pairing spectrum of emitter Signals incident peak, and then uses the smart search of small step progress row around these spectrum peaks.Smart search finishes and just obtains the accurate incident angle of emitter Signals.And then angle reported host computer.Report the next initiation command of waiting for digital receiver after finishing, enter next processing procedure again.
The entire process process is owing to adopted the hardware configuration and the efficient software algorithm fast of the multi-disc DSP parallel processing of shared bus, can in very short time, calculate the incident angle of radiation source, therefore when having very high direction finding precision, also has very strong real-time.
Claims (1)
1. direction-finding method based on the fast real-time space spectrum estimation ultra-resolution direction-finding device, the principal and subordinate DSP parallel signal processing module that the fast real-time space spectrum estimation ultra-resolution direction-finding device mainly is made of one group of DSP, form by Logic control module and communication interface modules three parts that FPGA constitutes, communication interface modules is reported and submitted to host computer from a sampled data and an angle that calculates that digital receiver receives each passage, the close coupled system that principal and subordinate DSP in the parallel signal processing module adopts bus to share connects, utilize data bus transmission intermediate data between the DSP, between dsp bus and the Logic control module by data, address and control bus are interconnected; At first receive the data of digital receiver sampling, data are carried out pre-service, obtain correlation matrix then; Correlation matrix is carried out feature decomposition, draw the noise characteristic subspace, and utilize eigenwert to estimate the number of radiation source; The spatial spectrum value that calculates the given angle zone with the direction steering vector and the noise characteristic subspace of different incidence angles degree correspondence then; Again the spectrum value is carried out spectrum peak search, the position correspondence at the spectrum peak that obtains the true incident direction of radiation source; It is characterized in that its concrete steps are: at first carry out initialization, at once calculate the steering vector form after initialization finishes, the result is left in internal memory on the sheet of DSP; Calculated after the steering vector, wait for the initiation command that digital receiver sends after sampling finishes, in case initiation command is arranged, the data that main DSP receives the digital receiver sampling also begin to enter handling procedure, and these data are 5 * 100 the matrixes that are made of 5 passages, 100 sampling snaps; Processor at first carries out relevant treatment to the data matrix and obtains 5 * 5 correlation matrix; Correlation matrix is carried out feature decomposition, obtain eigenwert and proper vector; Utilize eigenwert to estimate the number of radiation source by the information theory criterion; Utilize information source number that proper vector is divided, obtain signal characteristic vector and noise feature vector respectively; Noise feature vector is also with regard to corresponding noise characteristic subspace, and main DSP gives each sheet from DSP by data bus with noise feature vector; Principal and subordinate DSP begins computer memory spectrum value in angular regions separately then; Calculating the required steering vector of using of spectrum value directly calls in internal memory by lookup table mode; After the spatial spectrum value in whole zone was all calculated and finished, each DSP composed the peak coarse search with big step-length in angular regions separately; After search finishes, send spectrum peak position and spectrum peak that coarse search obtains back to main DSP again from DSP; Main DSP compares all coarse search spectrum peaks, determines the pairing spectrum of emitter Signals incident peak, and then uses the smart search of small step progress row around these spectrum peaks; Smart search finishes and just obtains the accurate incident angle of emitter Signals; And then angle reported host computer; Report the next initiation command of waiting for digital receiver after finishing, enter next processing procedure again.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008101374641A CN101403791B (en) | 2008-11-06 | 2008-11-06 | Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008101374641A CN101403791B (en) | 2008-11-06 | 2008-11-06 | Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101403791A CN101403791A (en) | 2009-04-08 |
| CN101403791B true CN101403791B (en) | 2011-01-26 |
Family
ID=40537862
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2008101374641A Expired - Fee Related CN101403791B (en) | 2008-11-06 | 2008-11-06 | Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101403791B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103217681B (en) * | 2013-03-19 | 2015-05-06 | 中国科学院声学研究所 | Tree-shaped topological mechanism multiprocessor sonar signal processing method |
| CN103278805B (en) * | 2013-05-21 | 2016-01-13 | 江苏科技大学 | Based on the search lighting Poewr control method that radio frequency is stealthy when guiding without priori |
| CN103795478B (en) * | 2014-01-16 | 2015-08-05 | 北京科技大学 | A kind of many primary users number detection method based on typical association analysis |
| CN105954707B (en) * | 2016-04-20 | 2018-06-12 | 东南大学 | A kind of method for the reduction direction finding multivalue problem that battle array is received based on isosceles trapezoid |
| CN109633525B (en) * | 2018-11-29 | 2021-01-19 | 上海无线电设备研究所 | Quick uniform circular array space spectrum direction finding realization method |
| CN111948598B (en) * | 2020-07-30 | 2022-11-29 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Method and device for detecting space domain interference signal |
| CN112821905B (en) * | 2020-12-30 | 2022-01-14 | 北京航空航天大学杭州创新研究院 | ADS-B receiver digital gain compensation system based on EVS and compensation method thereof |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101105525A (en) * | 2007-07-06 | 2008-01-16 | 哈尔滨工程大学 | Pure phase type broad frequency band microwave radiation source direction finding system and method |
| CN101109799A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Sense finding device for solid base line broad band microwave radiation source and sense finding method thereof |
| CN101109798A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Accurate sense finding device for P/L waveband radiation source and sense finding method thereof |
| CN101109797A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Five-time frequency interval ultra-wide frequency band sense finding method based on virtual array extension |
-
2008
- 2008-11-06 CN CN2008101374641A patent/CN101403791B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101105525A (en) * | 2007-07-06 | 2008-01-16 | 哈尔滨工程大学 | Pure phase type broad frequency band microwave radiation source direction finding system and method |
| CN101109799A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Sense finding device for solid base line broad band microwave radiation source and sense finding method thereof |
| CN101109798A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Accurate sense finding device for P/L waveband radiation source and sense finding method thereof |
| CN101109797A (en) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | Five-time frequency interval ultra-wide frequency band sense finding method based on virtual array extension |
Non-Patent Citations (4)
| Title |
|---|
| 司锡才等.舰载通信空间谱估计超分辨测向技术研究.《哈尔滨工程大学学报》.2002,第23卷(第5期),全文. * |
| 徐家雄.二维DOA估计的高速并行实现.《电子科技大学硕士论文》.2005,第30、38、41-46、48页. * |
| 白蕾，司锡才.空间圆阵列天线超分辨测向技术研究.《系统工程与电子技术》.2001,第23卷(第5期),全文. * |
| 谢纪岭,司锡才,唐建红.基于多级维纳滤波器的二维测向算法及DSP实现.《宇航学报》.2008,第29卷(第1期),第318页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101403791A (en) | 2009-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101403791B (en) | Fast real-time space spectrum estimation ultra-resolution direction-finding device and method thereof | |
| CN109212471B (en) | Positioning base station, system and method | |
| CN102411136B (en) | Phase interferometer direction finding method for ambiguity resolution by extension baselines | |
| CN110058202B (en) | High-performance hybrid positioning method based on angle and signal arrival time difference estimation and implementation device thereof | |
| CN107015191B (en) | One kind single dipole polarization sensitization array dimensionality reduction DOA estimation method under multi-path jamming environment | |
| CN102435978B (en) | Direction finding device for phase interferometer and phase spectrum interferometer based on multiple baselines | |
| CN111693978B (en) | Scatter detection method based on MIMO millimeter wave radar | |
| CN109633525A (en) | A kind of quick uniform circular array spatial spectrum direction finding implementation method | |
| CN107918122A (en) | A kind of portable scene surveillance radar system | |
| CN105676171A (en) | Single channel dual base station ultra-short wave signal spatial positioning method | |
| CN114488064B (en) | Distance-speed joint estimation method | |
| CN106896360A (en) | A kind of FPGA implementation method of SAR signal processing algorithms | |
| CN106707255B (en) | phased array radar simulation system and method | |
| CN109471068B (en) | Unmanned aerial vehicle positioning system and method based on radio frequency array signal DOA estimation | |
| CN111812630B (en) | System and method for detecting target and estimating DOA (direction of arrival) of external radiation source radar when interference remains | |
| CN103336279A (en) | Missile-borne SAR (synthetic aperture radar) imaging real-time signal processing system | |
| CN109444901A (en) | A kind of multiple submatrixes SAS sub-aperture image method under isomerous environment | |
| Huang et al. | An efficient FPGA implementation for 2-D MUSIC algorithm | |
| Pei et al. | QTT Ultra-wideband Signal Acquisition and Baseband Data Recording System Design Based on the RFSoC Platform | |
| CN113820701B (en) | High-frame-frequency rapid target detection method based on two-dimensional frequency-phase scanning array | |
| Gao et al. | Design of direction finding system of rotating multi-baseline phase interferometer based on GPU and DDS | |
| Wu et al. | Indoor single station 3D localization based on L-shaped sparse array | |
| CN202750234U (en) | Three-dimensional positioning system under synergistic action of wireless sensor networks and ultrasonic waves | |
| Xiaoguang | Signal Processing System of Back-Projection Algorithm with Multi GPU s | |
| CN216485477U (en) | Radar sparse array antenna |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110126 Termination date: 20171106 |