CN106707273A - Method for detecting multi-station radar signal fusion based on Neyman-Pearson rule digitalizing - Google Patents

Method for detecting multi-station radar signal fusion based on Neyman-Pearson rule digitalizing Download PDF

Info

Publication number
CN106707273A
CN106707273A CN201710050308.0A CN201710050308A CN106707273A CN 106707273 A CN106707273 A CN 106707273A CN 201710050308 A CN201710050308 A CN 201710050308A CN 106707273 A CN106707273 A CN 106707273A
Authority
CN
China
Prior art keywords
signal
radar
radar station
quantization
label
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
Application number
CN201710050308.0A
Other languages
Chinese (zh)
Other versions
CN106707273B (en
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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201710050308.0A priority Critical patent/CN106707273B/en
Publication of CN106707273A publication Critical patent/CN106707273A/en
Application granted granted Critical
Publication of CN106707273B publication Critical patent/CN106707273B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/87—Combinations of radar systems, e.g. primary radar and secondary radar
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/04—Systems determining presence of a target
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/35—Details of non-pulse systems
    • G01S7/352—Receivers
    • G01S7/354—Extracting wanted echo-signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention belongs to the technical field of radar and discloses a method for detecting multi-station radar signal fusion based on Neyman-Pearson rule digitalizing. The method comprises the following steps: arranging a multi-station radar system, including a signal fusion detection center and N radar stations; confirming a digitalizing bit b for a received echo signal by each radar station, thereby acquiring the quantity of digitalizing zones: M=2b; giving a false alarm probability of the signal fusion detection center and confirming the digital thresholds of the N radar stations; digitalizing the echo signal by the N radar stations according to the respective digital thresholds and transmitting a corresponding mark after the digitalizing to the signal fusion detection center; detecting a target by the signal fusion detection center according to the false alarm probability and the corresponding mark after the digitalizing for the echo signal by the N radar stations; and reducing the communication bandwidth required when the radar station echo signal is transmitted to the signal fusion detection center.

Description

Based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection method
Technical field
The present invention relates to Radar Technology field, more particularly to it is a kind of based on how graceful Pearson came (Neyman-Pearson) criterion The multistation Radar Signal Fusion detection method of quantization, can be used to reduce radar and signal fused in the detection of multistation Radar Signal Fusion Communication bandwidth between center.
Background technology
Multistation radar fence is made up of a signal fused center and multiple radar stations.Compared with single radar, multistation radar Net has the advantages that bigger investigative range, positioning precision higher, stronger survival ability and stronger antijamming capability.Cause This, the research of multistation radar fence becomes a new study hotspot.In the research of multistation radar fence, the signal of multistation radar Fusion detection method is an important research direction.
For current existing research, the signal fused detection method of multistation radar is broadly divided into three kinds：Multistation radar The signal level fusion detection method of flight path level fusing method, the judgement level fusing method of multistation radar and multistation radar.Multistation thunder The flight path level fusing method for reaching refers to each radar station transmits to signal fused center the flight path of oneself, and signal fused center is by institute There is the Track Fusion of radar station an into flight path.The judgement level fusing method of multistation radar refers to each radar station by respective judgement Result is transmitted to signal fused center, and signal fused center finally judges having for target according to the court verdict of all radars Nothing.The signal level fusion detection method of multistation radar refers to that each radar station transmits to fusion center, letter the echo-signal of oneself Number fusion center judges the presence or absence of target according to all radar station echo-signals.
In detection performance, the fusion method of the judgement level fusing method better than flight path level of multistation radar, and multistation radar Signal level fusion detection method again better than adjudicate level fusing method.At present, the flight path level fusing method of multistation radar into It is ripe, and it has been applied to actual engineering system.The goal in research of next step is that signal level fusing method is applied into new one The netted radar system in generation.For the signal level fusion detection method of multistation radar, due to each radar station, need will be all of Echo-signal is transmitted to signal fused center, it is therefore desirable to big communication bandwidth.But, between radar station and signal fused center Communication bandwidth it is limited, so new quantization method must be studied to reduce required communication band when radar echo signal is transmitted It is wide.And while communication bandwidth is reduced, it is ensured that the detection performance that signal fused center has had.Returned for different radars Ripple signal model, it has been proposed that various quantization criterions and corresponding Quantitative fusion detection algorithm.Such as, in gaussian signal model Under, to maximize Fisher information as criterion is quantified；Under Rayleigh signal model, to maximize Pasteur's distance as criterion is entered Row quantization etc..But, existing multistation radar Quantitative fusion detection algorithm detection performance loss is larger, it is therefore desirable to further grind Study carefully.
The content of the invention
Deficiency it is an object of the invention to be directed to above-mentioned prior art, proposes a kind of based on how graceful Pearson criterion quantifies Multistation Radar Signal Fusion detection method, to ensure under conditions of given false-alarm probability, to maximize signal fused center Detection probability for criterion quantify radar station echo-signal, reduce radar station echo-signal transmit to signal fused center when institute The communication bandwidth for needing.
To reach above-mentioned purpose, the present invention is adopted the following technical scheme that and is achieved.
It is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection method, it is characterised in that it is described Method comprises the following steps：
Step 1, sets netted radar system, and the netted radar system includes a signal fused inspection center and N number of Radar station；
Step 2, determines the quantization digit b of echo-signal of each radar station to receiving, and is obtained according to quantization digit b Quantized interval number M=2b, the corresponding label of each quantized interval is expressed as m, m=0,1,2 ..., M-1, the quantization digit b Meet：1≤b<Ch/fs, wherein, ChIt is the maximum communication bandwidth between each radar station and signal fused inspection center, fsIt is every The highest sample frequency of individual radar station；
Step 3, the false-alarm probability at Setting signal fusion detection center, according to the false-alarm probability and the quantized interval Number determines the quantization threshold of N number of radar station；
Step 4, N number of radar station quantifies according to each self-corresponding quantization threshold to echo-signal, and will be right after quantization The label answered is transmitted to signal fused inspection center；
Step 5, the signal fused inspection center is according to the false-alarm probability, N number of radar station to the echo-signal amount of carrying out Corresponding label, detects to target after change.
Advantages of the present invention is as follows：The present invention due under conditions of given false-alarm probability directly maximizing signal fused The detection probability at center is that criterion quantifies radar echo signal, therefore, compared with the conventional method, in the condition of same communication bandwidth Under, signal fused center can be made to obtain preferably detection property.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is provided in an embodiment of the present invention a kind of based on how the multistation Radar Signal Fusion that graceful Pearson criterion quantifies is examined Survey method realizes schematic flow sheet；
Fig. 2 is to realize schematic flow sheet when the present invention solves quantization threshold using majorized function；
Fig. 3 is that the inventive method is shown with detection probability contrast of the Pasteur apart from quantization method is maximized when quantization digit is 1 It is intended to；
Fig. 4 is that the inventive method is shown with detection probability contrast of the Pasteur apart from quantization method is maximized when quantization digit is 2 It is intended to；
Fig. 5 is that the inventive method is shown with detection probability contrast of the Pasteur apart from quantization method is maximized when quantization digit is 3 It is intended to；
Fig. 6 Fig. 3 be quantization digit be 4 when the inventive method with maximize Pasteur apart from quantization method and non-quantized letter The detection probability contrast schematic diagram of number fusion detection method.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made Embodiment, belongs to the scope of protection of the invention.
The embodiment of the present invention provide it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection method, As shown in figure 1, methods described comprises the following steps：
Step 1, sets netted radar system, and the netted radar system includes a signal fused inspection center and N number of Radar station.
Specifically, giving a netted radar system, netted radar system includes a signal fused center and N number of thunder Up to station；N number of radar station is designated as 1 to n-th radar station, the maximum allowed between each radar station and signal fused center Communication bandwidth is all identical；The maximum communication bandwidth allowed between radar station and signal fused center is designated as Ch, maximum communication band C widehUnit for bps；Highest sample frequency in N number of radar station is designated as fs, highest sample frequency fsUnit be conspicuous Hereby.
Step 2, determines the quantization digit b (quantizations of each radar station of echo-signal of each radar station to receiving Number is identical), and quantized interval number M=2 is obtained according to quantization digit bb, the corresponding label of each quantized interval is expressed as m, m =0,1,2 ..., M-1, the quantization digit b meet：1≤b<Ch/fs, wherein, ChIt is that each radar station and signal fused are detected Maximum communication bandwidth between center, fsIt is the highest sample frequency of each radar station.
On the basis of above formula is met, the value of quantization digit b is bigger, and the detection performance at signal fused center is better；But It is, when quantization digit b increases to a certain extent, further to increase quantization digit b and the detection performance at signal fused center is changed Kind very little, therefore, the experience value of quantization digit b is 3~5.
Step 3, the false-alarm probability at Setting signal fusion detection center, according to the false-alarm probability and the quantized interval Number determines the quantization threshold of N number of radar station.
Step 3 specifically includes following sub-step：
The echo-signal that i-th radar station is received is designated as x for N number of radar station by (3a)i, the noise of echo-signal Than being designated as λi, the test statistics of echo-signal is designated as yi=xi/μi, wherein, μiIt is i-th noise power of radar station, i=1, 2,…,N；
In practice, if signal to noise ratio λiWith noise power μiActual value it is unknown, can by estimate method obtain it Value.
M-1 quantization threshold of i-th radar station is designated as T by (3b) successivelyi,1,Ti,2,…,Ti,M-1, then i-th radar station The corresponding prior probability expression formula of m-th quantized interval be：
Wherein, viRepresent the test statistics y of the echo-signal that i-th radar station is receivediCorresponding label after quantization, H1Represent and assume that target is present, H0Represent and assume no target, P () represents probable value, and exp represents exponential function, P (vi=m | H1) represent and assume that target is present under conditions of, the test statistics y of the echo-signal that i-th radar station is receivediIt is right after quantization The label v for answeringiThe probability of=m, P (vi=m | H0) represent assume it is aimless under the conditions of, what i-th radar station was received returns The test statistics y of ripple signaliCorresponding label v after quantizationiThe probability of=m；
(3c) makes quantized interval corresponding label m=0,1,2 ..., M-1, so as to obtain i-th M quantization of radar station Interval corresponding prior probability expression formula；
I=1,2 ..., N are made, so as to obtain the prior probability expression formula that N number of radar station distinguishes corresponding M quantized interval；
(3d) distinguishes the prior probability expression formula of corresponding M quantized interval according to N number of radar station, obtains signal fused inspection Expression formula L (the v of the likelihood ratio of measured center1,v2,…,vN) be:
Wherein, m=0,1,2 ..., M-1, for any one radar station, the corresponding prior probability expression of its quantized interval Formula has M kind possibilities, so that the likelihood ratio of signal fused inspection center has MNPossibility is planted, ∏ is represented and even multiply symbol；
(3e) is according to the likelihood ratio expression formula of the signal fused inspection center, the false-alarm of the signal fused inspection center Probability, is calculated the likelihood ratio threshold value t and corresponding probable value α of likelihood ratio threshold value：
Wherein,Represent that the likelihood ratio to meeting signal fused inspection center exceedes likelihood ratio threshold value t All combinations sued for peace,Represent the likelihood to meeting signal fused inspection center Than all combinations more than likelihood ratio threshold value t assuming that the probability under without goal condition is sued for peace.
(3f) is according to the likelihood ratio of the signal fused inspection center, the likelihood ratio threshold value and the likelihood ratio thresholding It is worth corresponding probable value, obtains the detection probability P of signal fused inspection centerd：
(3g) is criterion to maximize the detection probability of the signal fused inspection center, builds and solves N number of radar station point Not corresponding M-1 quantization threshold Ti,1,Ti,2,…,Ti,M-1Optimized model：
s.t.0<Ti,1<…<Ti,M-1, i=1 ..., N
Wherein, f (Pd) it is detection probability PdFunction ensureing to minimize function f (Pd) be equivalent to maximize detection probability Pd；
(3h) solves the Optimized model, obtains N number of radar station and distinguishes corresponding M-1 quantization threshold Ti,1,Ti,2,…, Ti,M-1, i=1,2 ..., N.
Specifically, the Optimized model in step (3g) is the nonlinear programming problem of belt restraining, can use existing non- Linear optimization Algorithm for Solving；Because the fminimax functions in MATLAB softwares can fast and effectively solve the problem, therefore, This example solves the quantization threshold T for obtaining N number of radar station using the fminimax functions in MATLAB softwaresi,1,Ti,2,…, Ti,M-1, i=1,2 ..., N.
As shown in Fig. 2 in sub-step (3g), according to quantization threshold Ti,1,Ti,2,…,Ti,M-1, the optimization of i=1,2 ..., N Model solution obtains the quantization threshold T of N number of radar stationi,1,Ti,2,…,Ti,M-1, the detailed process of i=1,2 ..., N is：
(3g1) sets cycle-index Nc, an interim storage matrix T for N rows M-1 row is set, minimum target function is set Value fmin=1, cycle-index mark k=1 is set；
(3g2) sets random quantization threshold initial valueI=1,2 ..., N；
(3g3) is according to quantization threshold initial valueI=1,2 ..., N, using in MATLAB softwares Fminimax functions, the quantization threshold value for obtaining this circulation is solved according to the Optimized model in step (3g)I=1,2 ..., the target function value f after N and minimumk；
(3g4) if minimize after target function value fk<fmin, then interim storage matrix is made
Minimum target functional value fmin=fk；
(3g5) is if cycle-index mark k<Nc, then the value of cycle-index mark k adds 1, and goes to step (3g2)；It is no Then, step (3g6) is gone to；
(3g6) is according to interim storage matrix T, the quantization threshold T after being solvedi,1=T (i, 1), Ti,2=T (i, 2) ..., Ti,M-1=T (i, M-1), i=1,2 ..., N；T (i, j) represents the element of the i-th row jth row of interim storage matrix T.
Exemplary, in sub-step (3g), build and solve the corresponding M-1 quantization threshold T of N number of radar station differencei,1, Ti,2,…,Ti,M-1Optimized model：
s.t.0<Ti,1<…<Ti,M-1, i=1 ..., N
Wherein, f (Pd) it is detection probability PdFunction ensureing to minimize function f (Pd) be equivalent to maximize detection probability Pd, function f (Pd) expression formula be f (Pd)=1-PdOr f (Pd)=1/Pd。
Step 4, N number of radar station quantifies according to each self-corresponding quantization threshold to echo-signal, and will be right after quantization The label answered is transmitted to signal fused inspection center.
Step 4 specifically includes following sub-step：
(4a) i-th radar station is according to its corresponding M-1 amountChange doorLimit Ti,1,Ti,2,…,Ti,M-1, i-th radar station is connect The test statistics y of the echo-signal for receivingiQuantified；If test statistics yiMeet Ti,m≤yi<Ti,m+1, then system is checked Metering yiLabel v after quantizationi=m；
(4b) i-th radar station is by test statistics yiLabel v after quantizationiTransmit to signal fused inspection center；
(4c) makes i=1,2 ..., N, so that N number of radar station is by the label v after quantization1,v2,…,vNTransmit to signal fused Inspection center.
Step 5, the signal fused inspection center is according to the false-alarm probability, N number of radar station to the echo-signal amount of carrying out Corresponding label, detects to target after change.
Step 5 specifically includes following sub-step：
(5a) obtains signal fused according to quantized interval number M and quantized interval corresponding label m=0,1,2 ..., M-1 What the N number of label of inspection center was sued for peace is possible to value vsum={ 0,1,2 ..., N × (M-1) }, wherein, × represent multiplication sign；
Assuming that signal fused inspection center label summation value is v during no targetsumProbability be P (vsum|H0)；
(5b) is according to the false-alarm probability Pfa, solved by following formula and obtain detection threshold g and probable value γ：
Wherein,Represent to meeting label summation value vsumAll labels more than detection threshold g Summation value vsumAssuming that the probability under without goal condition is sued for peace,Expression is asked meeting label With value vsumAll labels summation value v equal to detection threshold gsumAssuming that the probability under without goal condition is asked With.
The label v of N number of radar station that (5c) signal fused inspection center will receive1,v2,…,vNSummation, obtains label SumIf label sum vf>G, then judgement has target；If label sum vf=g, then adjudicated with probable value γ There is target；If label sum vf<G, then adjudicate no target.
Effect of the invention is tested by following simulation comparison and further illustrated：
1. simulation parameter is set：Netted radar system includes 1 signal fused inspection center and 5 radar stations, this 5 The echo-signal of radar station has identical signal to noise ratio, and the span of signal to noise ratio is 0dB to 20dB, the value interval of signal to noise ratio It is 0.5dB；Quantization digit b=4, i.e. quantized interval number M=16；The false-alarm probability of signal fused inspection center is set to 10-6； When the quantization threshold of each radar station is solved, cycle-index N is setc=20.This example quantifies with maximization Pasteur's distance Signal fused detection method is contrasted.
2. emulation content：
Quantization digit b=1, i.e. quantized interval number M=2 are set, are set according to simulation parameter, calculate signal to noise ratio value model Enclose the corresponding detection probability curve of interior the inventive method as shown in the square marks line in Fig. 3；Give maximization Pasteur's distance Detection probability curve during quantization is as shown in the addition marks line in Fig. 3.
Quantization digit b=2, i.e. quantized interval number M=4 are set, are set according to simulation parameter, calculate signal to noise ratio value model Enclose the corresponding detection probability curve of interior the inventive method as shown in the square marks line in Fig. 4；Give maximization Pasteur's distance Detection probability curve during quantization is as shown in the addition marks line in Fig. 4.
Quantization digit b=3, i.e. quantized interval number M=8 are set, are set according to simulation parameter, calculate signal to noise ratio value model Enclose the corresponding detection probability curve of interior the inventive method as shown in the square marks line in Fig. 5；Give maximization Pasteur's distance Detection probability curve during quantization is as shown in the addition marks line in Fig. 5.
Quantization digit b=4, i.e. quantized interval number M=16 are set, are set according to simulation parameter, calculate signal to noise ratio value In the range of the corresponding detection probability curve of the inventive method as shown in the square marks line in Fig. 6；Give maximization Pasteur away from From detection probability curve when quantifying as shown in the addition marks line in Fig. 6；In order to contrast, the detection given when not quantifying is general Rate curve is as shown in the circular mark line in Fig. 6.
3. analysis of simulation result：
From Fig. 3 to Fig. 6, identical quantization digit is given, detected with the signal fused for maximizing Pasteur's distance quantization Method is compared, and the inventive method can obtain preferably detection performance.In the case where detection probability is 0.5, with maximization bar The signal fused detection method that family name's distance quantifies is compared, and when quantization digit is 1, the detection performance of this method is carried as shown in Figure 3 Rise 2.081dB；When quantization digit is 2, the detection performance boost 0.723dB of this method as shown in Figure 4；When quantization digit is 3 When, the detection performance boost 0.306dB of this method as shown in Figure 5；When quantization digit is 4, the detection of this method as shown in Figure 6 Performance boost 0.106dB.
In addition, it will be appreciated from fig. 6 that when quantization digit be 4 when, the detection probability curve of the inventive method with do not quantify when inspection Survey probability curve to be sufficiently close to, in the case where detection probability is 0.5, the detection performance loss of the inventive method only has 0.051dB。
Summary emulation experiment can be seen that the present invention according to how graceful Pearson criterion quantifies to signal, and existing What is had is compared based on the quantization method for maximizing Pasteur's distance, under identical quantization digit, can obtain preferably detection property Energy.
The above, specific embodiment only of the invention, but protection scope of the present invention is not limited thereto, and it is any Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all contain Cover within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.

Claims (5)

1. it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection method, it is characterised in that the side Method comprises the following steps：
Step 1, sets netted radar system, and the netted radar system includes a signal fused inspection center and N number of radar Stand；
Step 2, determines the quantization digit respectively b of echo-signal of each radar station to receiving, and is obtained according to quantization digit b To quantized interval number M, M=2b, the corresponding label of each quantized interval is expressed as m, m=0,1,2 ..., M-1；Institute State quantization digit b satisfactions：1≤b ＜ Ch/fs；Wherein, ChFor the maximum between each radar station and signal fused inspection center is led to Letter bandwidth, fsIt is the highest sample frequency of each radar station；
Step 3, the false-alarm probability at Setting signal fusion detection center is true according to the false-alarm probability and the quantized interval number Fixed N number of each self-corresponding quantization threshold of radar station；
Step 4, N number of radar station quantifies respectively according to each self-corresponding quantization threshold to the echo-signal being respectively received, And corresponding label after quantization is transmitted separately to signal fused inspection center；
Step 5, after the signal fused inspection center quantifies according to the false-alarm probability, N number of radar station to echo-signal Corresponding label, detects to target, and then knows in the echo-signal that N number of radar station is received with the presence or absence of target.
2. it is according to claim 1 it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection side Method, it is characterised in that step 3 specifically includes following sub-step：
The echo-signal that i-th radar station is received is designated as x for N number of radar station by (3a) respectivelyi, i-th radar station is connect The signal to noise ratio of the echo-signal for receiving is designated as λi, the test statistics of the echo-signal that i-th radar station is received is designated as yi= xi/μi；Wherein, μiIt is i-th noise power of radar station, i=1,2 ..., N；
M-1 quantization threshold of i-th radar station is designated as T by (3b) successivelyI, 1, TI, 2..., TI, m..., TI, M-1, then i-th thunder It is up to the corresponding prior probability expression formula of m-th quantized interval stood：
P ( v i = m | H 1 ) = exp &lsqb; - T i , m / ( 1 + &lambda; i ) &rsqb; - exp &lsqb; - T i , m + 1 / ( 1 + &lambda; i ) &rsqb; P ( v i = m | H 0 ) = exp ( - T i , m ) - exp ( - T i , m + 1 )
Wherein, viRepresent the test statistics y of the echo-signal that i-th radar station is receivediCorresponding label, H after quantization1Table Show that hypothesis target is present, H0Represent and assume no target, P () represents probable value, and exp represents exponential function；P(vi=m | H1) Represent under conditions of assuming that target is present, the test statistics y of the echo-signal that i-th radar station is receivediCorrespondence after quantization Label viThe probability of=m；P(vi=m | H0) represent assume it is aimless under the conditions of, the echo that i-th radar station is received The test statistics y of signaliCorresponding label v after quantizationiThe probability of=m；M=0,1,2 ..., M-1；
(3c) makes quantized interval corresponding label m=0,1,2 ..., M-1, so as to obtain i-th M quantized interval of radar station Corresponding prior probability expression formula；
I=1,2 ..., N are made, so as to obtain the prior probability expression formula that N number of radar station distinguishes corresponding M quantized interval；
(3d) distinguishes the prior probability expression formula of corresponding M quantized interval according to N number of radar station, obtains during signal fused detects Expression formula L (the v of the likelihood ratio of the heart1, v2..., vN) be：
L ( v 1 , v 2 , ... , v N ) = &Pi; i = 1 N P ( v i = m | H 1 ) P ( v i = m | H 0 )
Wherein, m=0,1,2 ..., M-1, for any one radar station, the corresponding prior probability expression formula of its quantized interval has M Possibility is planted, so that the likelihood ratio of signal fused inspection center has MNPossibility is planted, Π is represented and even multiply symbol；
(3e) is general according to the likelihood ratio expression formula of the signal fused inspection center, the false-alarm of the signal fused inspection center Rate, is calculated the likelihood ratio threshold value t and corresponding probable value α of likelihood ratio threshold value respectively：
P f a = &Sigma; L ( v 1 , v 2 , ... , v N ) > t &Pi; i = 1 N P ( v i = m | H 0 ) + &alpha; &times; &Sigma; L ( v 1 , v 2 , ... , v N ) = t &Pi; i = 1 N P ( v i = m | H 0 )
Wherein,Represent that the likelihood ratio to meeting signal fused inspection center exceedes all of likelihood ratio threshold value t Combination is sued for peace；
(3f) is according to the likelihood ratio of the signal fused inspection center, the likelihood ratio threshold value and the likelihood ratio threshold value pair The probable value answered, obtains the detection probability P of signal fused inspection centerd：
P d = &Sigma; L ( v 1 , v 2 , ... , v N ) > t &Pi; i = 1 N P ( v i = m | H 1 ) + &alpha; &times; &Sigma; L ( v 1 , v 2 , ... , v N ) = t &Pi; i = 1 N P ( v i = m | H 1 )
(3g) is criterion to maximize the detection probability of the signal fused inspection center, builds the N number of radar station of solution right respectively The M-1 quantization threshold T for answeringI, 1, TI, 2..., TI, M-1Optimized model：
m i n T i , 1 , T i , 2 , ... , T i , M - 1 i = 1 , ... , N f ( P d )
S.t.0 ＜ TI, 1＜ ... ＜ TI, M-1, i=1 ..., N
Wherein, f (Pd) it is detection probability PdFunction ensureing to minimize function f (Pd) be equivalent to maximize detection probability Pd；
(3h) solves the Optimized model, obtains N number of radar station and distinguishes corresponding M-1 quantization threshold TI, 1, TI, 2..., TI, M-1, I=1,2 ..., N.
3. it is according to claim 1 it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection side Method, it is characterised in that in sub-step (3g), builds and solves the corresponding M-1 quantization threshold T of N number of radar station differenceI, 1, TI, 2..., TI, M-1Optimized model：
m i n T i , 1 , T i , 2 , ... , T i , M - 1 i = 1 , ... , N f ( P d )
S.t.0 ＜ TI, 1＜ ... ＜ TI, M-1, i=1 ..., N
Wherein, f (Pd) it is detection probability PdFunction, and minimize function f (Pd) be equivalent to maximize detection probability Pd, function f (Pd) expression formula be f (Pd)=1-PdOr f (Pd)=1/Pd。
4. it is according to claim 1 it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection side Method, it is characterised in that step 4 specifically includes following sub-step：
(4a) i-th radar station is according to its corresponding M-1 quantization threshold TI, 1, TI, 2..., TI, M-1, i-th radar station is received The test statistics y of the echo-signal for arrivingiQuantified；If test statistics yiMeet TI, m≤yi＜ TI, m+1, then system is checked Metering yiLabel v after quantizationi=m；
(4b) i-th radar station is by test statistics yiLabel v after quantizationiTransmit to signal fused inspection center；
(4c) makes i=1,2 ..., N, so that N number of radar station is by the label v after quantization1, v2..., vNTransmit to signal fused detection Center.
5. it is according to claim 1 it is a kind of based on how graceful Pearson criterion quantify multistation Radar Signal Fusion detection side Method, it is characterised in that step 5 specifically includes following sub-step：
(5a) obtains signal fused detection according to quantized interval number M and quantized interval corresponding label m=0,1,2 ..., M-1 What the N number of label in center was sued for peace is possible to value vsum={ 0,1,2 ..., N × (M-1) }；Wherein, × represent multiplication sign；
Assuming that signal fused inspection center label summation value is v during no targetsumProbability be P (vsum|H0)；
(5b) is according to the false-alarm probability Pfa, solved by following formula and obtain detection threshold g and probable value γ：
P f a = &Sigma; v s u m > g P ( v s u m | H 0 ) + &gamma; &times; &Sigma; v s u m = g P ( v s u m | H 0 )
Wherein,Represent to meeting label summation value vsumAll labels summation more than detection threshold g takes Value vsumAssuming that the probability under without goal condition is sued for peace,Represent to meeting label summation value vsumAll labels summation value v equal to detection threshold gsumAssuming that the probability under without goal condition is sued for peace；
The label v of N number of radar station that (5c) signal fused inspection center will receive1, v2..., vNSummation, obtains label sumIf label sum vf＞ g, then adjudicating in the echo-signal that N number of radar station is received has target；If label it And vf=g, then have target in adjudicating the echo-signal that N number of radar station is received with probable value γ；If label sum vf＜ g, then Adjudicating in the echo-signal that N number of radar station is received does not have target.
CN201710050308.0A 2017-01-23 2017-01-23 Based on how graceful Pearson criterion quantization multistation Radar Signal Fusion detection method Active CN106707273B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710050308.0A CN106707273B (en) 2017-01-23 2017-01-23 Based on how graceful Pearson criterion quantization multistation Radar Signal Fusion detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710050308.0A CN106707273B (en) 2017-01-23 2017-01-23 Based on how graceful Pearson criterion quantization multistation Radar Signal Fusion detection method

Publications (2)

Publication Number Publication Date
CN106707273A true CN106707273A (en) 2017-05-24
CN106707273B CN106707273B (en) 2019-05-21

Family

ID=58908860

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710050308.0A Active CN106707273B (en) 2017-01-23 2017-01-23 Based on how graceful Pearson criterion quantization multistation Radar Signal Fusion detection method

Country Status (1)

Country Link
CN (1) CN106707273B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107607926A (en) * 2017-10-31 2018-01-19 西安电子科技大学 A kind of distributed radar low traffic calibration signal merges target acquisition processing method
CN108039940A (en) * 2017-11-24 2018-05-15 天津大学 A kind of invalid CCE elimination methods of LTE system PDCCH
CN110956191A (en) * 2018-09-27 2020-04-03 财团法人工业技术研究院 Classifier, classification method and classification system based on probability fusion
CN111812591A (en) * 2020-06-16 2020-10-23 南京航空航天大学 Target detection method based on Bayesian estimation
CN113406583A (en) * 2021-06-22 2021-09-17 电子科技大学长三角研究院(衢州) Approximate calculation method for cloud MIMO radar target detection probability
CN113434816A (en) * 2020-09-21 2021-09-24 重庆工商大学 Method for detecting signal under noise enhancement neman-pearson criterion
CN113534096A (en) * 2021-07-19 2021-10-22 东莞理工学院 LiDAR signal geometric feature extraction method and system based on spline function
US11984964B2 (en) 2021-05-10 2024-05-14 Raytheon Company Decentralized control via adaptive importance encoding

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7079072B1 (en) * 1987-01-23 2006-07-18 Raytheon Company Helicopter recognition radar processor
CN102298140A (en) * 2011-05-16 2011-12-28 中国人民解放军海军航空工程学院 Radar fence anti-interference usefulness assessment method based on detection probability geometric distribution curve
KR101546421B1 (en) * 2015-02-16 2015-08-24 에스티엑스엔진 주식회사 Adaptive constant false alarm rate processing method
CN105425223A (en) * 2015-11-11 2016-03-23 西安电子科技大学 Detection method of sparse distance extension radar target in generalized Pareto clutter
CN105954739A (en) * 2016-04-20 2016-09-21 电子科技大学 Knowledge-aided nonparametric constant false alarm detection method
CN106199588A (en) * 2016-06-24 2016-12-07 西安电子科技大学 The multistation Radar Signal Fusion detection method quantified based on Pasteur's distance

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7079072B1 (en) * 1987-01-23 2006-07-18 Raytheon Company Helicopter recognition radar processor
CN102298140A (en) * 2011-05-16 2011-12-28 中国人民解放军海军航空工程学院 Radar fence anti-interference usefulness assessment method based on detection probability geometric distribution curve
KR101546421B1 (en) * 2015-02-16 2015-08-24 에스티엑스엔진 주식회사 Adaptive constant false alarm rate processing method
CN105425223A (en) * 2015-11-11 2016-03-23 西安电子科技大学 Detection method of sparse distance extension radar target in generalized Pareto clutter
CN105954739A (en) * 2016-04-20 2016-09-21 电子科技大学 Knowledge-aided nonparametric constant false alarm detection method
CN106199588A (en) * 2016-06-24 2016-12-07 西安电子科技大学 The multistation Radar Signal Fusion detection method quantified based on Pasteur's distance

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
彭志刚等: "分布式ＭＩＭＯ 雷达信号级量化融合检测方法研究", 《现代雷达》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107607926B (en) * 2017-10-31 2020-07-03 西安电子科技大学 Method for detecting and processing low-traffic quasi-signal fusion target of distributed radar
CN107607926A (en) * 2017-10-31 2018-01-19 西安电子科技大学 A kind of distributed radar low traffic calibration signal merges target acquisition processing method
CN108039940A (en) * 2017-11-24 2018-05-15 天津大学 A kind of invalid CCE elimination methods of LTE system PDCCH
CN108039940B (en) * 2017-11-24 2020-12-11 天津大学 Invalid CCE (control channel element) removing method for PDCCH (physical downlink control channel) of LTE (long term evolution) system
CN110956191B (en) * 2018-09-27 2023-08-18 财团法人工业技术研究院 Classifier, classification method and classification system based on probability fusion
CN110956191A (en) * 2018-09-27 2020-04-03 财团法人工业技术研究院 Classifier, classification method and classification system based on probability fusion
CN111812591A (en) * 2020-06-16 2020-10-23 南京航空航天大学 Target detection method based on Bayesian estimation
CN113434816A (en) * 2020-09-21 2021-09-24 重庆工商大学 Method for detecting signal under noise enhancement neman-pearson criterion
US11984964B2 (en) 2021-05-10 2024-05-14 Raytheon Company Decentralized control via adaptive importance encoding
CN113406583A (en) * 2021-06-22 2021-09-17 电子科技大学长三角研究院(衢州) Approximate calculation method for cloud MIMO radar target detection probability
CN113406583B (en) * 2021-06-22 2022-08-02 电子科技大学长三角研究院(衢州) Approximate calculation method for cloud MIMO radar target detection probability
CN113534096B (en) * 2021-07-19 2023-09-15 东莞理工学院 LiDAR signal geometric feature extraction method and system based on spline function
CN113534096A (en) * 2021-07-19 2021-10-22 东莞理工学院 LiDAR signal geometric feature extraction method and system based on spline function

Also Published As

Publication number Publication date
CN106707273B (en) 2019-05-21

Similar Documents

Publication Publication Date Title
CN106707273A (en) Method for detecting multi-station radar signal fusion based on Neyman-Pearson rule digitalizing
Morss et al. Examining the use of weather forecasts in decision scenarios: Results from a US survey with implications for uncertainty communication
CN108664632A (en) A kind of text emotion sorting algorithm based on convolutional neural networks and attention mechanism
Feijóo et al. Simulation of correlated wind speeds: A review
CN102413474B (en) Self-adaption trust management system and method of cognitive radio network
Zheng et al. Channel non-line-of-sight identification based on convolutional neural networks
CN107180260B (en) Short wave communication frequency selecting method based on Evolutionary Neural Network
CN107967487A (en) A kind of colliding data fusion method based on evidence distance and uncertainty
CN106199588B (en) Multistation Radar Signal Fusion detection method based on Pasteur&#39;s distance quantization
Hu et al. Two‐stage constant false alarm rate detection for distributed multiple‐input–multiple‐output radar
CN110222513A (en) A kind of method for monitoring abnormality of Above-the-line, device and storage medium
Huang et al. An information diffusion technique to assess integrated hazard risks
Chagas et al. Genetic algorithms and simulated annealing optimization methods in wireless sensor networks localization using artificial neural networks
Rhif et al. Deep learning models performance for NDVI time series prediction: a case study on north west Tunisia
CN116050504A (en) Wind power short-term prediction model based on deep learning
Bengio et al. A connectionist system for medium-term horizon time series prediction
CN107371175A (en) A kind of self-organizing network fault detection method using cooperation prediction
Lee et al. Evaluating the effectiveness of sequential aviation security screening policies
CN109686429A (en) Physician visits period recommended method and device
Tokuyama et al. The effect of using attribute information in network traffic prediction with deep learning
Deng et al. Multi‐period probabilistic‐scenario risk assessment of power system in wind power uncertain environment
CN106304108B (en) MIMO radar mobile platform fast reserve dispositions method based on variation monitoring demand
Bayo et al. The percolating cluster is invisible to image recognition with deep learning
Huang et al. Ship trajectory anomaly detection based on multi-feature fusion
Arif et al. Machine Learning and Deep Learning Based Network Slicing Models for 5G Network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant