CN103235294A - Method for estimating weak signal separation on basis of positioning for external radiation sources - Google Patents

Method for estimating weak signal separation on basis of positioning for external radiation sources Download PDF

Info

Publication number
CN103235294A
CN103235294A CN2013101068558A CN201310106855A CN103235294A CN 103235294 A CN103235294 A CN 103235294A CN 2013101068558 A CN2013101068558 A CN 2013101068558A CN 201310106855 A CN201310106855 A CN 201310106855A CN 103235294 A CN103235294 A CN 103235294A
Authority
CN
China
Prior art keywords
signal
target
nlms
time delay
signals
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
Application number
CN2013101068558A
Other languages
Chinese (zh)
Inventor
李万春
张和发
李莽
李立萍
魏平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN2013101068558A priority Critical patent/CN103235294A/en
Publication of CN103235294A publication Critical patent/CN103235294A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a method for estimating weak signal separation on the basis of positioning for external radiation sources. The method includes steps of performing direct-wave and multi-path cancellation for a sample sequence x (t) of a main channel by means of NLMS (normalized least mean square) to obtain output signals y<NLMS>(t); estimating the number of source signals of the output signals y<NLMS>(t) by an AIC (Akaike information criterion) or an MDL (minimum description length) criterion, preliminarily whitening the output signals y<NLMS>(t) to obtain whitened signals z(t), and updating a separation matrix W by an ICA (independent component analysis) algorithm to obtain M channels of estimated signals y(t); respectively performing time-frequency two-dimensional correlation for the M channels of estimated signals y<1>(t), ..., y<M>(t) and reference signals s<ref>(t); and searching a peak point in absolute correlation values phi<i>(tau, f) corresponding to the M channels of estimated signals, and determining time delay and Doppler frequency shift which correspond to the peak point as time delay and Doppler frequency shift of a target when a signal-to-noise ratio of a certain channel of signals is larger than or equal to a threshold value.

Description

A kind of feeble signal based on location, external radiation source is separated method of estimation
Technical field
The invention belongs to signal processing technology, be specifically related to the feeble signal isolation technics.
Background technology
The passive location technology, because transmitting-receiving separates, itself does not initiatively transmit, and only is that information by passive reception external radiation source realizes the technology of locating.Compare with active localization method, it is far away that the passive location method has operating distance, and interference performance is strong, has the ability of lobe-on-receive, location and tracking to target.The passive location technology can be divided into radiation source location and location, external radiation source.The radiation source location mainly is to utilize the signal of target self institute radiation to position, location, external radiation source then is to utilize third party's radiation or external radiation, general civilian TV signal, broadcast singal, signal of communication, the satellite that adopts known geographical locations, and civilian radiation source such as cellular base station etc.Be characterized in, if utilize a plurality of external radiation source signals, will greatly strengthen the power range of receiving station and to the reconnaissance capability of low detectable target.
Outside in the research of radiation source positioning system, be that the passive location systematic study in external radiation source is paid close attention to by lot of domestic and foreign research institution and scholar with the broadcast singal, become a research focus.The principle signal of such positioning system as shown in Figure 1, radar is to receiving direct wave, the echo of moving target, the direct wave of other common-frequency broadcasting launching tower transmission and the echo of other object reflection that including of signal comes from broadcast emission tower itself, such as the still life echo of various buildings reflections.Therefore, the target echo signal that radar receives in the signal is extremely faint, and is submerged in the various interference.In actual applications, be the enough power ranges of assurance system, receiver need have high sensitivity to receive faint moving-target reflected signal.On the other hand, because the region of civilian broadcast system has broadcast emission tower at the emission broadcast singal at target detection zone periphery, and, when the employed broadcast singal of these Signals ﹠ Systems frequency together, will form co-channel interference.Therefore, the signal that the broadcast singal of launching except the launching tower of known location in the signal of receiving in receiving station's antenna main lobe and this signal interested moving-target of process and still life reflect, also comprise the common-frequency broadcasting signal of close region broadcast emission tower, this signal becomes the undesired signal of positioning system.
Traditional co-channel interference signal inhibition method is by airspace filter the direction of arrival of co-channel interference to be carried out trap, thereby its filtering from receive signal is gone out, and again the signal after the filtering is done time delay, Doppler shift, direction of arrival DOA estimation to target.Yet, when the aspect ratio of co-channel interference signal and target near the time, traditional method can also filter out target simultaneously, can not make estimation to the parameters of target effectively again.
Summary of the invention
Technical matters to be solved by this invention is, proposes under a kind of strong jamming and the low signal-to-noise ratio condition can to realize extracting the weak target echoed signal based on the feeble signal separation method of location, broadcast singal external radiation source.
The present invention is that the technical scheme that solves the problems of the technologies described above the lock employing is, a kind of feeble signal based on location, external radiation source is separated method of estimation, may further comprise the steps:
Step 1: the take over party carries out bandpass sampling respectively to the modulation signal that the pending main channel array modulation signal that receives and reference antenna receive, and obtains main channel sample sequence x (t) and reference sample sequence s Ref(t), wherein t represents sampling instant, and sampling length is L;
Step 2: utilize normalization minimum mean-square two to take advantage of NLMS implementation direct wave and multipath to offset to main channel sample sequence x (t) and obtain output signal output signal y NLMS(t);
Step 3: utilize red pond information criterion AIC or minimum description length MDL criterion to output signal y NLMS(t) implement the source signal number and estimate, estimate source signal number M, then to output signal y NLMS(t) carry out signal z (t) after prewhitening obtains albefaction, use independent component analysis ICA algorithm to upgrade separation matrix W, the M road estimated signal y (t) that obtains, y (t)=Wz (t)=y 1(t) ..., y M(t);
Step 4: with M road estimated signal y 1(t) ..., y M(t) respectively with reference signal s Ref(t) do that time-frequency two-dimensional is relevant to obtain corresponding correlation | Ψ i(τ, f) |, | &Psi; i ( &tau; , f ) | = | &Sigma; t = 0 L - 1 y i ( t ) s ref * ( t - &tau; ) e j 2 &pi;f / L | , I=3 ..., M, * represents conjugation, and L is bandpass sampling length, and τ is time delay, f Doppler shift;
Step 5: respectively at the correlation of M road estimated signal correspondence | Ψ i(τ, f) | middle search peak point, with the time delay of each road signal peak value point correspondence and Doppler shift as a useful signal, calculate each road signal to noise ratio (S/N ratio), as a certain road signal to noise ratio (S/N ratio) during more than or equal to threshold value, time delay and the Doppler shift that time delay and the Doppler of its peak point correspondence is defined as a target frequently then, otherwise judge and do not have target.
The present invention has adopted based on the whitening approach that is applicable to the blind extraction of feeble signal, can be implemented in and to realize under strong jamming and the low signal-to-noise ratio condition extracting interested weak target echoed signal, help back two dimension coherent matched filter to extract time delay and the Doppler shift of target.
Further, the present invention also carries out following steps by echoed signal match estimating target direction of arrival DOA after step 5:
Step 6: utilize the time delay of target and Doppler shift that reference signal is carried out match, obtain target echo signal s (t),
Figure BDA00002986801900022
τ wherein dFor estimating the time delay of target, f dFor estimating the Doppler shift of target;
Step 7: the DOA that utilizes target echo signal s (t) and main channel sample sequence x (t) to carry out known waveform estimates, estimates the DOA of target.
Use echoed signal match estimating target direction of arrival DOA can realize the measurement parameter pairing, avoid the data association, reduced calculated amount, more accurately estimating target DOA.
The invention has the beneficial effects as follows, can be implemented in and accurately extract interested weak target echoed signal under strong jamming and the low signal-to-noise ratio condition, to realize location and the tracking to target.
Description of drawings
Fig. 1: based on the system principle synoptic diagram of location, broadcasting external radiation source.
Fig. 2: signal handling principle figure of the present invention.
Fig. 3: the time-frequency two-dimensional correlogram of estimating target one during initial time.
Fig. 4: the time-frequency two-dimensional correlogram of estimating target two during initial time.
Embodiment
Signal processing flow is as shown in Figure 2:
Step 1: the take over party carries out bandpass sampling respectively to the modulation signal that the pending main channel array modulation signal that receives and reference antenna receive, and obtains main channel sample sequence x (t) and reference sample sequence s Ref(t), wherein t represents sampling instant, and sampling length is L;
Step 2: utilize normalization minimum mean-square two to take advantage of NLMS implementation direct wave and multipath to offset to main channel sample sequence x (t) and obtain output signal output signal y NLMS(t); Direct wave and multipath offset to the known technology in the existing co-channel interference signal inhibition method, do not give unnecessary details at this;
Step 3: utilize red pond information criterion AIC or minimum description length MDL criterion to output signal y NLMS(t) implement the source signal number and estimate, estimate source signal number M, then to output signal y NLMS(t) carry out signal z (t) after prewhitening obtains albefaction, use the quick independent component analysis nc-FastICA of non-round complex signal algorithm to upgrade separation matrix W, the M road estimated signal y (t) that obtains, y (t)=Wz (t)=y 1(t) ..., y M(t); Other existing ICA algorithms are equally applicable to the calculating of separation matrix W, and are quicker than other ICA algorithm at this nc-FastICA algorithm of selecting for use;
Step 4: with M road estimated signal y 1(t) ..., y M(t) respectively with reference signal s Ref(t) do that time-frequency two-dimensional is relevant to obtain corresponding correlation | Ψ i(τ, f) |, | &Psi; i ( &tau; , f ) | = | &Sigma; t = 0 L - 1 y i ( t ) s ref * ( t - &tau; ) e j 2 &pi;f / L | , I=1 ..., M, * represents conjugation, and L is bandpass sampling length, and τ is time delay, f Doppler shift;
Step 5: respectively at the correlation of M road estimated signal correspondence | Ψ i(τ, f) | middle search peak point (τ d, f d), with each road signal peak value point (τ d, f d) corresponding time delay and Doppler shift be as a useful signal, calculates each road signal to noise ratio (S/N ratio)
Figure BDA00002986801900041
Figure BDA00002986801900042
Be illustrated in the estimated signal of i road and do not comprise | Ψ id, f d) | all | Ψ i(τ, f) | average, as a certain road signal to noise ratio (S/N ratio) during more than or equal to threshold value 13dB, time delay and the Doppler shift that time delay and the Doppler of its peak point correspondence is defined as a target frequently then, otherwise judge and do not have target;
Step 6: utilize the time delay of target and Doppler shift that reference signal is carried out match, obtain target echo signal s (t),
Figure BDA00002986801900043
τ wherein dFor estimating the time delay of target, f dFor estimating the Doppler shift of target;
Step 7: the DOA that utilizes target echo signal s (t) and main channel sample sequence x (t) to carry out known waveform estimates, estimates the DOA of target.1 to M road signal carries out relevantly among the target echo signal s (t) that at first calculates match and the main channel sample sequence x (t), obtains correlation matrix
Figure BDA00002986801900044
r ^ i = ( 1 L &Sigma; t = 1 L x i ( t ) s * ( t ) ) ( &Sigma; t = 1 L s ( t ) s H ( t ) ) - 1 I=1 ..., M, * represents conjugation, H represents conjugate transpose; Or
Figure BDA00002986801900046
Figure BDA00002986801900047
Be matrix matched filtering function, i:i+m kThe i to i+m of the 1st row of matrix is got in-1 expression k-1 row, H ^ = ( &Sigma; t = 1 L x ( t ) s H ( t ) ) ( &Sigma; t = 1 L s ( t ) s H ( t ) ) - 1 , m kFor carrying out consequent level and smooth submatrix array number earlier; It is level and smooth to make forward direction, constructor array compound matrice W,
Figure BDA00002986801900049
R=M-m k, before and after remaking to smoothly
Figure BDA000029868019000410
Wherein,
Figure BDA000029868019000411
Be m k* m kOppose the angular moment battle array; Afterwards, ask W FbAutocorrelation matrix, and carry out eigenwert and decompose EVD, obtain the eigenwert decomposition result
Figure BDA000029868019000412
Calculate the order R of Q, get in the matrix V the 1st to R row and obtain matrix V k=V (:, 1:R); The pseudo-spectrum of structure is carried out spectrum peak search at last
Figure BDA000029868019000413
a 1(θ) be array flow pattern vector, the maximum of points θ of search MaxBe exactly the DOA of target.
In order to verify effect of the present invention, be the stereophonic FM broadcasting signal with the independent random signal before mixing, the sound signal passband that is used for modulation is 30Hz~15KHz, maximum frequency deviation ± 75Hz, the pre-emphasis time constant is 50 μ s, the sample frequency of simulate signal is 882KHz, and modulation back signal passband is 200KHz.The sound signal that is used for modulation in first frequency modulation broadcasting simulate signal is one section symphony, the centre frequency of modulation back signal is 200KHz, second fm broadcast signal is one section soprano's solo for the sound signal of modulating, the centre frequency of modulation back signal is 200.3KHz, and the power of two fm broadcast signals all remains 1.Two signals highly overlap on frequency domain, can't be separated with the bandpass filter of routine.Emulation platform operating system is WINXP, adopts MATLAB to compile.
The Track In Track experiment that the signal processing flow of embodiment carries out under the binocular target situation, simulate signal parameter and simulate signal target component are as follows:
Table 1: simulate signal parameter
Signal Time delay (e-4) Amplitude fading (dB) DOA(°) Doppler shift (Hz)
Main direct wave s 1 0 0 -50 ?
Target echo s 2 3.064 -50 -40 17.81
Target echo s 3 3.2285 -50 -30 30.85
Interference due to the adjacent station s 4 0 -20 49 ?
Main platform one multipath s 5 2.26e-2 -20 20 ?
Main platform two multipath s 6 4.52e-2 -25 10 ?
Table 2: simulate signal target component
Parameter Distance (km) The orientation (°) Speed (m/s) Heading (°) Decay (dB)
Target one 92 -40 90 165 -50
Target two 82 -30 150 95 -50
In the experiment, the distance of establishing between frequency modulation cell site and the radar station is 60km, and establishing initial time is 0, and every the 10s once sampling, the sampling duration is 1s, 9 groups of sampled datas.
Table 3 and table 4 are listed the contrast between target one and target two parameters estimated values and the actual value respectively.
Table 3: the parameter comparison of target one
Figure BDA00002986801900051
Table 4: the parameter comparison of target two
Figure BDA00002986801900052
The time-frequency two-dimensional correlogram of Fig. 3 and Fig. 4 provide respectively when being initial time estimating target one, target two, more accurate and peak value is clearly to the parameter estimation of target one as can be seen, and from table 3 and table 4 as can be seen present embodiment can be good at estimating the parameters of target, the truer flight path of simulated target.

Claims (5)

1. the feeble signal based on location, external radiation source is separated method of estimation, it is characterized in that, may further comprise the steps:
Step 1: the take over party carries out bandpass sampling respectively to the modulation signal that the pending main channel array modulation signal that receives and reference antenna receive, and obtains main channel sample sequence x (t) and reference sample sequence s Ref(t), wherein t represents sampling instant, and sampling length is L;
Step 2: utilize normalization minimum mean-square two to take advantage of NLMS implementation direct wave and multipath to offset to main channel sample sequence x (t) and obtain output signal output signal y NLMS(t);
Step 3: utilize red pond information criterion AIC or minimum description length MDL criterion to output signal y NLMS(t) implement the source signal number and estimate, estimate source signal number M, then to output signal y NLMS(t) carry out signal z (t) after prewhitening obtains albefaction, use independent component analysis ICA algorithm to upgrade separation matrix W, the M road estimated signal y (t) that obtains, y (t)=Wz (t)=y 1(t) ..., y M(t);
Step 4: with M road estimated signal y 1(t) ..., y M(t) respectively with reference signal s Ref(t) do that time-frequency two-dimensional is relevant to obtain corresponding correlation | Ψ i(τ, f) |, | &Psi; i ( &tau; , f ) | = | &Sigma; t = 0 L - 1 y i ( t ) s ref * ( t - &tau; ) e j 2 &pi;f / L | , I=3 ..., M, * represents conjugation, and L is bandpass sampling length, and τ is that time delay, f are Doppler shift;
Step 5: respectively at the correlation of M road estimated signal correspondence | Ψ i(τ, f) | middle search peak point, with the time delay of each road signal peak value point correspondence and Doppler shift as a useful signal, calculate each road signal to noise ratio (S/N ratio), as a certain road signal to noise ratio (S/N ratio) during more than or equal to threshold value, time delay and the Doppler shift that time delay and the Doppler of its peak point correspondence is defined as a target frequently then, otherwise judge and do not have target.
2. a kind of feeble signal based on location, external radiation source is separated method of estimation according to claim 1, it is characterized in that, also comprises after the step 5:
Step 6: utilize the time delay of target and Doppler shift that reference signal is carried out match, obtain target echo signal s (t),
Figure FDA00002986801800012
τ wherein dFor estimating the time delay of target, f dFor estimating the Doppler shift of target;
Step 7: the DOA that utilizes target echo signal s (t) and main channel sample sequence x (t) to carry out known waveform estimates, estimates the DOA of target.
3. a kind of feeble signal based on location, external radiation source is separated method of estimation according to claim 1, it is characterized in that described independent component analysis ICA algorithm is specially the quick independent component analysis nc-FastICA of non-round complex signal algorithm.
4. a kind of feeble signal based on location, external radiation source is separated method of estimation according to claim 1, it is characterized in that each road signal to noise ratio (S/N ratio) is
Figure FDA00002986801800021
Wherein, Be illustrated in the estimated signal of i road and do not comprise | Ψ id, f d) | all | Ψ i(τ, f) | average.
5. separate method of estimation as a kind of feeble signal based on location, external radiation source as described in the claim 4, it is characterized in that described threshold value is 13dB.
CN2013101068558A 2013-03-29 2013-03-29 Method for estimating weak signal separation on basis of positioning for external radiation sources Pending CN103235294A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013101068558A CN103235294A (en) 2013-03-29 2013-03-29 Method for estimating weak signal separation on basis of positioning for external radiation sources

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013101068558A CN103235294A (en) 2013-03-29 2013-03-29 Method for estimating weak signal separation on basis of positioning for external radiation sources

Publications (1)

Publication Number Publication Date
CN103235294A true CN103235294A (en) 2013-08-07

Family

ID=48883345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013101068558A Pending CN103235294A (en) 2013-03-29 2013-03-29 Method for estimating weak signal separation on basis of positioning for external radiation sources

Country Status (1)

Country Link
CN (1) CN103235294A (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105787081A (en) * 2016-03-03 2016-07-20 电子科技大学 Radiation platform correlation method based on radiation source space positions
CN106054143A (en) * 2016-07-12 2016-10-26 西安电子科技大学 Outer radiation source radar cofrequency interference elimination method
CN107144811A (en) * 2017-05-12 2017-09-08 电子科技大学 A kind of cyclic subspace direction-finding method of single channel receiving array signal
CN107515382A (en) * 2017-08-21 2017-12-26 电子科技大学 A kind of moving-target DOA trackings based on phase compensation
CN107850655A (en) * 2015-05-22 2018-03-27 诺基亚技术有限公司 Packet prepares
CN108919239A (en) * 2018-04-10 2018-11-30 华南理工大学 It is a kind of to offset method based on double ultrasonic direct waves for combining adaptive-filterings
CN109471065A (en) * 2018-09-28 2019-03-15 中国电子科技集团公司第三十六研究所 A kind of direction-finding method of coherent signal
CN109521457A (en) * 2018-09-25 2019-03-26 中国辐射防护研究院 A kind of gamma emitter item division methods and system based on information criterion
CN109642949A (en) * 2016-07-13 2019-04-16 德克萨斯仪器股份有限公司 Method and apparatus for narrowband ranging system
CN111030762A (en) * 2019-12-23 2020-04-17 北京华力创通科技股份有限公司 Method and device for measuring chirp signal power
CN111198374A (en) * 2020-01-13 2020-05-26 哈尔滨工程大学 Doppler sensitive signal moving target underwater sound detection method based on space-time-frequency joint interference suppression
CN111323744A (en) * 2020-03-19 2020-06-23 哈尔滨工程大学 Target number and target angle estimation method based on MDL (minimization drive language) criterion
CN111381208A (en) * 2018-12-29 2020-07-07 西安电子科技大学 Method for estimating direction of arrival of known waveform information source
CN114280593A (en) * 2019-07-22 2022-04-05 华为技术有限公司 Radar system and vehicle

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张和发: "《盲信号分离技术及其在无源定位中的应用》", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107850655A (en) * 2015-05-22 2018-03-27 诺基亚技术有限公司 Packet prepares
CN105787081A (en) * 2016-03-03 2016-07-20 电子科技大学 Radiation platform correlation method based on radiation source space positions
CN105787081B (en) * 2016-03-03 2019-02-15 电子科技大学 A kind of radiation platform correlating method based on radiation source spatial position
CN106054143A (en) * 2016-07-12 2016-10-26 西安电子科技大学 Outer radiation source radar cofrequency interference elimination method
CN109642949A (en) * 2016-07-13 2019-04-16 德克萨斯仪器股份有限公司 Method and apparatus for narrowband ranging system
CN109642949B (en) * 2016-07-13 2024-04-19 德克萨斯仪器股份有限公司 Method and apparatus for a narrowband ranging system
CN107144811A (en) * 2017-05-12 2017-09-08 电子科技大学 A kind of cyclic subspace direction-finding method of single channel receiving array signal
CN107515382A (en) * 2017-08-21 2017-12-26 电子科技大学 A kind of moving-target DOA trackings based on phase compensation
CN107515382B (en) * 2017-08-21 2019-12-06 电子科技大学 Moving target DOA tracking method based on phase compensation
CN108919239A (en) * 2018-04-10 2018-11-30 华南理工大学 It is a kind of to offset method based on double ultrasonic direct waves for combining adaptive-filterings
CN108919239B (en) * 2018-04-10 2022-07-26 华南理工大学 Ultrasonic direct wave cancellation method based on double-combination adaptive filtering
CN109521457A (en) * 2018-09-25 2019-03-26 中国辐射防护研究院 A kind of gamma emitter item division methods and system based on information criterion
CN109471065B (en) * 2018-09-28 2021-05-25 中国电子科技集团公司第三十六研究所 Direction finding method for coherent signals
CN109471065A (en) * 2018-09-28 2019-03-15 中国电子科技集团公司第三十六研究所 A kind of direction-finding method of coherent signal
CN111381208A (en) * 2018-12-29 2020-07-07 西安电子科技大学 Method for estimating direction of arrival of known waveform information source
CN111381208B (en) * 2018-12-29 2021-11-09 西安电子科技大学 Method for estimating direction of arrival of known waveform information source
CN114280593A (en) * 2019-07-22 2022-04-05 华为技术有限公司 Radar system and vehicle
CN114280593B (en) * 2019-07-22 2024-05-14 华为技术有限公司 Radar system and vehicle
CN111030762A (en) * 2019-12-23 2020-04-17 北京华力创通科技股份有限公司 Method and device for measuring chirp signal power
CN111198374A (en) * 2020-01-13 2020-05-26 哈尔滨工程大学 Doppler sensitive signal moving target underwater sound detection method based on space-time-frequency joint interference suppression
CN111323744A (en) * 2020-03-19 2020-06-23 哈尔滨工程大学 Target number and target angle estimation method based on MDL (minimization drive language) criterion
CN111323744B (en) * 2020-03-19 2022-12-13 哈尔滨工程大学 Target number and target angle estimation method based on MDL (minimization drive language) criterion

Similar Documents

Publication Publication Date Title
CN103235294A (en) Method for estimating weak signal separation on basis of positioning for external radiation sources
CN110412559B (en) Non-coherent fusion target detection method for MIMO radar of distributed unmanned aerial vehicle
CN109188470B (en) GNSS deception jamming detection method based on convolutional neural network
CN103713277B (en) A kind of radiation source localization method of position-based information field
CN103217670B (en) Outer radiation source weak signal detection method based on PCA (principal component analysis)
CN106125053A (en) Pulse Doppler radar polarization anti jamming method
CN104007421B (en) Rowland C passive radar TOA method of estimation based on total variation and compressed sensing
CN104515971A (en) Airborne single-station passive positioning method for multiple broadband targets
CN104142496B (en) Based on the statistics MIMO radar multi-target orientation method that connected domain divides
CN105891815A (en) Combined estimation algorithm based on broadcast signal passive positioning
CN109444896B (en) Underwater sound array positioning system and positioning method thereof
CN104155662B (en) The mutual disturbance restraining method of self adaptation based on GNSS correlation peak detector
Chabriel et al. Adaptive target detection techniques for OFDM-based passive radar exploiting spatial diversity
CN107272030B (en) Target detection method based on information fusion of multiple GPS satellites
CN103293518A (en) Positioning and detection method of radiation source outside broadcast signals
CN105866750A (en) Method for detecting multiple GPS (global positioning system) satellite weak echo signals
CN102087354A (en) Passive radar grouping LS-CLEAN weak target detection method
CN102353947A (en) Method for estimating target echo signal subspaces of passive radars based on CSA-MWF (correlation subtraction algorithm-multistage wiener filter)
CN107561507A (en) A kind of clutter cancellation method of external illuminators-based radar
CN105891817A (en) Distributed passive radar target detection method under direct wave-free condition
CN105607096A (en) Dual-satellite time difference and frequency difference positioning method and positioning device
CN114415210B (en) Multi-dimensional information combined satellite navigation deception jamming source direct positioning method
CN106054165A (en) Distributed passive radar target detection realization method
CN104597434A (en) Multi-frame coherent TBD method for improving envelope shifting compensation and fractional order Fourier transformation
Lee et al. Neural network-based ranging with LTE channel impulse response for localization in indoor environments

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130807