CN110231616A - Sea Moving objects detection and location method based on big-dipper satellite radiation source - Google Patents

Sea Moving objects detection and location method based on big-dipper satellite radiation source Download PDF

Info

Publication number
CN110231616A
CN110231616A CN201910278789.XA CN201910278789A CN110231616A CN 110231616 A CN110231616 A CN 110231616A CN 201910278789 A CN201910278789 A CN 201910278789A CN 110231616 A CN110231616 A CN 110231616A
Authority
CN
China
Prior art keywords
signal
big
result
echo
radiation source
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
CN201910278789.XA
Other languages
Chinese (zh)
Other versions
CN110231616B (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.)
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 CN201910278789.XA priority Critical patent/CN110231616B/en
Publication of CN110231616A publication Critical patent/CN110231616A/en
Application granted granted Critical
Publication of CN110231616B publication Critical patent/CN110231616B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The present invention discloses a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source, is applied to Radar Moving Target detection field, for the problem low as external sort algorithm target echo signal-to-noise ratio using big-dipper satellite of the existing technology;The present invention is all made of symmetrical Keystone to the target echo from different big-dipper satellites received first and converts, and obtains range migration correction result;Then using sub-frame processing and oblique-FFT is gone to the echo-signal after each range migration correction, obtains the accumulation result of each echo-signal;Secondly the accumulation result of each echo is biradical apart from back projection to local coordinate according to difference, obtain the accumulation of multiple echoes as a result, when realizing more big-dipper satellites as radiation source, the effective accumulation for the target echo signal energy that different big-dipper satellites obtain;Finally according to the accumulation of multiple echoes of step S3 as a result, being detected and being positioned to target.

Description

Sea Moving objects detection and location method based on big-dipper satellite radiation source
Technical field
The invention belongs to Radar Moving Target detection field, in particular to a kind of passive radar fortune based on satellite radiation source Moving-target detection and location technology.
Background technique
Passive radar is a kind of special bis- (more) bases Passive Radar Systems.It is different from traditional radar, passive radar Not transmission of electromagnetic signals itself, but target acquisition is carried out using the electromagnetic wave of target emanation.No due to passive radar itself Emit electromagnetic wave, therefore it is with strong concealment;Meanwhile can reduce volume, reduce cost.In recent years, as passive radar shines The increase of the signal kinds and quantity in source is penetrated, the research of passive radar also increasingly increases.
Existing passive radar be largely using ground signal as signal source, but since it can only cover land area, And sea area cannot be completely covered, it is not particularly suited for sea and monitors on a large scale.In response to this problem, it can choose big-dipper satellite letter Number be used as radiation source.Big Dipper satellite signal has the advantages that other spatial electromagnetic signals are incomparable: (1) satellite as radiation source Signal can cover sea;(2) strong interference immunity;(3) satellite-signal feature is it is known that be convenient for target detection;(4) big-dipper satellite is China Independent research, use Big Dipper satellite signal high safety.Therefore being monitored using Big Dipper satellite signal to sea is one Feasible developing direction.
Since the power of big-dipper satellite transmitting signal is lower and extremely remote apart from ground, so that Big Dipper satellite signal reaches in ground Power density is low when ball surface, to keep the sea Moving Target Return Signal-to-Noise received extremely low, moves mesh for sea Target detection brings difficulty.
In order to improve the signal-to-noise ratio of sea Moving Target Return signal, document " Multi-frame fractional Fourier transform technique for moving target detection with space-based Passive radar.Iet Radar Sonar&Navigation, 2017,11 (5): 822-828. " proposes a kind of long observation The method for carrying out energy accumulation to target echo signal under time, first adjust the distance compressed echo-signal of this method carry out single order Keystone transformation, realizes range migration correction;Then the data after range migration correction are divided into multiframe, to each frame number According to FrFT is used, the coherent accumulation to moving target energy is realized;Frequency alignment and incoherent product are carried out between last multiframe data It is tired, make signal energy High accumulation to the domain Doppler frequency center-doppler frequency rate (DC-DFR).But the single order that the method uses Keystone transformation can only eliminate single order range walk, and there are still can not accurately realize range migration school to second order range curvature Just.Meanwhile there is also computationally intensive, the cumbersome disadvantages of process for the method.In addition, the method only has studied single satellite conduct The case where radiation source, cannot achieve the positioning of sea moving target.
Summary of the invention
In order to solve the above technical problems, the invention proposes a kind of, the passive radar sea based on big-dipper satellite radiation source is transported Moving-target detection and localization method, can effectively realize using more big-dipper satellites as the detection of the sea moving target of radiation source and Positioning.
A kind of the technical solution adopted by the present invention are as follows: passive radar sea moving target inspection based on big-dipper satellite radiation source Survey and localization method, comprising:
S1, be all made of symmetrical Keystone to the target echo from different big-dipper satellites received and convert, obtain away from Result is corrected from migration;
S2, using sub-frame processing and oblique-FFT is gone to the echo-signal after each range migration correction, obtains each echo-signal Accumulation result;
S3, by the accumulation result of each echo according to difference it is biradical apart from back projection to local coordinate, and in local seat No-coherence cumulating is carried out to the accumulation result of each echo under mark system;
S4, according to the no-coherence cumulating of multiple echoes of step S3 as a result, being detected and being positioned to target.
Further, step S1 includes:
S11, system parameter initialization, the parameter of the initialization include: distance to sample frequency, the target observation time, Equivalent signal pulse-recurrence time;
S12, using more big-dipper satellites as radiation source, direct wave and target echo are enrolled respectively with two antennas And it demodulates;
S13, relevant treatment is carried out to obtaining each big-dipper satellite as the corresponding direct wave of radiation source and target echo, obtain Each big-dipper satellite as the corresponding target echo of radiation source distance to the compressed result of pulse;
S14, symmetrical Keystone transformation is carried out to the compressed result of pulse to the distance of each target echo, obtains each mesh Mark the corresponding range migration correction result of echo.
Further, step S14 include it is following step by step:
A1, distance is carried out to the corresponding range migration correction result of each target echo to Fast Fourier Transform (FFT), obtain each Target echo is corresponding apart from frequency domain-orientation time-domain signal;
A2, corresponding to each target echo positive second order Keystone transformation is carried out apart from frequency domain-orientation time-domain signal;
A3, corresponding to each target echo negative second order Keystone transformation is carried out apart from frequency domain-orientation time-domain signal;
The result that A4, the step A2 corresponding to each target echo are obtained is multiplied with the obtained result of step A3;
Multiplied result that A5, the step A4 corresponding to each target echo are obtained carries out distance to inverse fast Fourier transform, To obtain the corresponding range migration correction result of each target echo.
Further, step S2 includes:
S21, the corresponding range migration correction result of some target echo is taken same apart from gate signal, it is corresponding obtains its Bearing signal;
S22, it goes tiltedly to handle and is divided into multiframe to what the bearing signal was unified;
S23, the coherent accumulation of signal and the no-coherence cumulating of interframe signal in frame are carried out to the bearing signal after going tiltedly;
S24, step S23- is executed to the Distance Remaining gate signal of the corresponding range migration correction result of the target echo S24 obtains single echo-signal in the accumulation result of the biradical distance-Doppler tune frequency domain of difference.
Further, step S22 include it is following step by step:
B1, the range that doppler frequency rate is set and tune separation values;
B2, one group of linear FM signal is constituted according to tune separation values;
B3, the bearing signal of taking-up is multiplied from the linear FM signal of different frequency modulation rates respectively, obtains 2D signal;
B4, with the time span of a determining frame, the obtained 2D signal of step B3 is divided into several frames.
Further, step S23 include it is following step by step:
C1, corresponding value of every frame signal each orientation time is directly added, is obtained every under the conditions of different tune frequency values The zero-frequency value of frame signal;
C2, the zero-frequency value of every frame signal under the conditions of different tune frequency values is subjected to no-coherence cumulating, obtained in different frequency modulation The zero-frequency value of signal under the conditions of rate value.
Further, step S3 includes:
S31, building local coordinate, reference axis is respectively position X, position Y and doppler frequency rate;
S32, by Ocean Scenes grid division, obtain the space coordinate of each grid;
S33, it is corresponding that different big-dipper satellites are obtained according to big-dipper satellite and the spatial position of receiving station for each grid The biradical distance of difference, and then the corresponding distance of different satellites is obtained to the time;
S34, the corresponding distance of different satellites that step S33 is obtained is brought into incoherent product that step S23 is obtained to the time Tire out as a result, obtaining the corresponding one-dimensional doppler frequency rate signal of each grid;
S35, the corresponding one-dimensional doppler frequency rate signal of big-dipper satellite each under current grid is directly subjected to incoherent product It is tired;
The no-coherence cumulating of each big-dipper satellite is as a result, obtain in local under S36, all grids obtained according to step S35 Three dimensional signal beneficial effects of the present invention in coordinate system: method of the invention, firstly, being defended to what is received from different Beidous The target echo of star is all made of symmetrical Keystone transformation, obtains range migration correction result;Secondly, to each range migration correction Echo-signal afterwards is using framing and goes oblique-FFT processing, obtains the accumulation result of each echo-signal;Then, by the product of each echo Tired result is biradical apart from back projection to local coordinate according to difference, obtains the accumulation result of multiple echoes;Finally, according to more The accumulation of a echo is as a result, detecting target and being positioned;Big-dipper satellite is efficiently solved as external sort algorithm target echo The low problem of signal-to-noise ratio;The method of the present invention has following advantages:
1, the detection and positioning present invention firstly provides multiple big-dipper satellites as the sea moving target of external sort algorithm;
2, using symmetrical Keystone convert, can in the case where unknown Doppler parameter simultaneously remove range walk and Range migration correction is better achieved in range curvature, while making Doppler frequency center zero;
3, when carrying out the accumulation processing of single satellite backward energy, directly corresponding value of each orientation time is added, greatly Calculation amount is reduced greatly;
4, according to the biradical distance of difference, the accumulation result for the sea Moving Target Return that multiple big-dipper satellites are obtained is reversed Local coordinate is projected to, and direct under local coordinate to the accumulation result of the sea Moving Target Return of each big-dipper satellite No-coherence cumulating is carried out, effective accumulation to different echo-signal energy is better achieved.
Detailed description of the invention
Fig. 1 is flow diagram provided in an embodiment of the present invention;
Fig. 2 is the passive radar geometric configuration schematic diagram provided in an embodiment of the present invention based on big-dipper satellite radiation source;
Fig. 3 is energy accumulation result schematic diagram of the target echo provided in an embodiment of the present invention in local coordinate;
Fig. 4 is the X-Y diagrammatic cross-section at accumulation result peak value provided in an embodiment of the present invention.
Specific embodiment
The present invention is extremely low using big-dipper satellite as the Moving Target Return Signal-to-Noise of radiation source received in order to solve The problem of, increase the observation time of target.Under long observation time, range walk and range curvature cannot be ignored, this hair It is bright to be converted using symmetrical Keystone, in the case where unknown Doppler parameter, range walk and range curvature can be eliminated simultaneously, Realize accurate range migration correction.Due to the variation of the target scattering characteristics under long observation time, can not directly be concerned with Accumulation, the present invention use sub-frame processing method, signal in frame is considered as it is relevant, in frame signal carry out coherent accumulation, frame Between signal carry out no-coherence cumulating.The present invention goes tiltedly to handle to the signal after symmetrical Keystone conversion process using unified, For removing oblique signal completely, its energy can all be run up at Doppler frequency center by Fourier transformation, and it is symmetrical Keystone is converted so that the Doppler frequency center of signal becomes zero.It therefore, can be different direction to the time for signal in frame Signal be directly added, to realize energy accumulating at zero-frequency;Simultaneously as different frame goes tiltedly to handle using unified, it is more General Le mass center is constant, therefore the incoherent processing of interframe does not need to carry out Doppler frequency center compensation deals.Consideration is defended with single Beidou Star is bad as the target echo energy accumulation effect of radiation source, and target can not be detected in noise background.The present invention use with More big-dipper satellites carry out target echo energy accumulation as radiation source.When using more big-dipper satellites as radiation source, by upper The processing of the step of face is available using different big-dipper satellites as the accumulation of the sea Moving Target Return of radiation source as a result, again that it is anti- To local coordinate is projected to, no-coherence cumulating is carried out.The present invention has target echo energy accumulation effect good, can carry out simultaneously The outstanding advantages of target detection and positioning.
For convenient for those skilled in the art understand that technology contents of the invention, with reference to the accompanying drawing to the content of present invention into one Step is illustrated.
As shown in Fig. 2, receiving station is fixed, the big-dipper satellite as radiation source selects the geometric configuration schematic diagram of the present embodiment For middle round Earth's orbit (MEO) satellite, the B3I signal for selecting big-dipper satellite to broadcast is as external sort algorithm signal, carrier frequency fcFor 1268.520MHz, signal bandwidth 20.46MHz.Select 5 big-dipper satellites as radiation source, sea fortune in the present embodiment Moving-target is expressed as (1000,0,0) m in reference moment position vector.
The technical solution adopted by the present invention process is as shown in Figure 1, include six parts, successively are as follows: 1, system parameter initialization; 2, direct wave and target echo are enrolled, line-spacing descriscent pulse of going forward side by side compression;3, symmetrical Keystone transformation;4, that unifies removes Xie Chu Reason and framing;5, single satellite accumulates;6, the 2D signal back projection in step 5 is sat to local according to difference biradical distance Mark system, and no-coherence cumulating is carried out to the accumulation result for the sea Moving Target Return that different big-dipper satellites obtain.
The method that the present invention mainly uses emulation experiment is verified, and all steps, conclusion are all verified just on Matlab Really.
It is of the invention the specific implementation process is as follows:
1, system parameter initializes:
The parameter specifically initialized includes: distance to sample rate fs, equivalent signal pulse recurrence frequency PRF, target observation Time Ta;F in the present embodimentsValue is 50MHz, and PRF value is 1000Hz, T in the present embodiment in the present embodimentaValue is 60s。
2, direct wave and target echo are enrolled, line-spacing descriscent pulse of going forward side by side compression:
1. the big-dipper satellite number as radiation source is M=5, the case where the m big-dipper satellite is radiation source is considered, to straight It is enrolled up to wave, demodulation obtains base band direct-path signal sd,m(τ, η), enrolls target echo, and demodulation obtains base band mesh Mark echo-signal sr,m(τ, η), wherein τ is distance to the time, and η is the orientation time, and η value range is [- 30,30] s, m= 1,2 ..., M, subscript d indicate that direct wave, subscript r indicate target scattering echo.
2. carrying out relevant treatment to direct-path signal and target echo signal, target echo signal distance is obtained to pulse pressure Result after contracting are as follows:
Wherein, subscript pc indicates pulse compression, σηFor the complex scattering coefficients of target, for value with time change, c is light Speed, ρ () are the cross-correlation function of echo-signal and through reference signal wave, and ρ () is triangular wave in the present embodiment, and R (η) is Difference is biradical apart from history, may be expressed as:
R (η)=Rt(η)+Rr(η)-Rb(η)
Wherein, Rt(η) and Rr(η) respectively indicates big-dipper satellite and receiving station's range-to-go, Rb(η) indicates that Beidou is defended The distance between star and receiving station, R (η) can be indicated with the second Taylor series are as follows:
Wherein, R0Indicate R (η) value at η=0, A, B respectively indicate R (η) about the single order of η, second dervative in η=0 The value at place.
3, symmetrical Keystone transformation:
1. being obtained distance to Fast Fourier Transform (FFT) apart from frequency domain-orientation to the compressed result progress distance of pulse The signal in domain:
Wherein, P () is the Fourier transformation of ρ (), and has:
In above formula, fτCoupling with η causes range walk, fτWith η2Coupling cause range curvature.
2. the result in 1. is carried out positive second order Keystone transformationSimultaneously for transformed knot Fruit replaces η therein with ηm, it obtains:
Wherein,
From the above equation, we can see that positive second order Keystone transformation eliminates fτWith η2Coupling, that is, eliminate range curvature.
3. the result in 1. is carried out negative second order Keystone transformationSimultaneously for transformed knot Fruit replaces η therein with ηm, it obtains:
Wherein,
From the above equation, we can see that negative second order Keystone transformation has also been removed range curvature.
4. by 2. result with 3. in result be multiplied, obtain:
Wherein,
From the above equation, we can see that symmetrical Keystone transformation may be implemented to eliminate range walk and range curvature simultaneously.
5. carrying out distance to inverse fast Fourier transform to the result 4. obtained:
sSKT,m(τ, η)=IFFTrg{SSKT,m(fτ,η)}
Wherein, IFFTrgIndicate distance to inverse Fourier transform.
4, that unifies goes tiltedly processing and framing:
1. taking out sSKT,mSame distance obtains its corresponding bearing signal, is denoted as s to the time in (τ, η)m(η), by step 3 can obtain smThe Doppler frequency center of (η) is zero, therefore only needs to consider the true value of doppler frequency rateAndValue not Know.Setting includesDoppler frequency rate range, doppler frequency rate range is [- 1.25,1.25] Hz/ in the present embodiment S, interval delta fdrFor 0.005Hz/s;
2. according to the doppler frequency rate value f in 1.dr, one group of LFM of composition (linear frequency modulation, Linear frequency modulation) signal, i.e. sLFM(fdr, η) and=exp (- j π fdrη2), by this group of LFM signal sLFM(fdr, η) and sm(η) is multiplied, and obtains To 2D signal sdechirp,m(fdr,η);
3. choosing the time span T that 2s is every frame in the present embodimentfra, by sdechirp,m(fdr, η) it is divided into along the orientation time 30 frames consider the i-th frame signal sm,i(fdr, η), the range for corresponding to η is [- 32+2i, -30+2i] s.Since signal in frame can be with Ignore complex scattering coefficients σηBring influences, can be by the i-th frame signal in the case where not considering amplitude and ignoring constant phase It indicates are as follows:
5, coherent accumulation is carried out to signal in frame, interframe signal carries out no-coherence cumulating:
1. by every frame signal sm,i(fdr, η) and in different fdrUnder each orientation time corresponding value be added, obtain different fdrUnder The zero-frequency value of every frame signal, i.e.,As selected doppler frequency rate fdrWith actual Doppler Frequency modulation rateIt is equal, i.e., when going oblique completely, only considers in the value range of η, ignore constant phase, it can be by every frame signal It indicates are as follows:
From the above equation, we can see that can by Fourier transformation by the energy accumulation of the signal at zero-frequency, directly will not here Zero-frequency value required for corresponding value addition obtains with the orientation time accumulates result to get to signal energy in frame.
2. will 1. obtained in result carry out interframe no-coherence cumulating, obtain different fdrUnder the conditions of 60s signal corresponding zero Frequency is worth, i.e.,Since when going oblique completely, each frame signal has run up to energy at zero-frequency, because This interframe signal does not need to carry out Doppler frequency center compensation, directly carries out the zero-frequency value that no-coherence cumulating obtains and believes to get to 60s Number energy accumulation result.
It is unknown to the time since there are the distances of target, it needs from sSKT,mLetter of the different distance to the time is taken out in (τ, η) Number obtain corresponding sm(η), and 4~5 process flow is repeated, obtain the two dimension accumulation result of single echo.By different echoes It is accumulated, obtains the two dimension accumulation result s of multiple echoesm(τ,fdr), m=1,2 ..., M;
6, local coordinate is constructed, according to the biradical distance of difference, by the 2D signal back projection in 5 to local coordinate system System, and no-coherence cumulating is carried out to the accumulation result for the sea Moving Target Return that different big-dipper satellites obtain:
1. constructing local coordinate, coordinate is respectively position X, Y and doppler frequency rate;
2. by Ocean Scenes grid division, if the space coordinate of each grid is (x, y, 0).
3. taking a grid from 2., is fixed due to satellite position it is known that receiving station location, it is corresponding that different satellites can be obtained Difference distance, and then obtain the corresponding fast time τ of different satellites1,...,τm,...,τM
4. the result in 3. is substituted into the 2D signal s in step 5m(τ,fdr) in, obtain the corresponding one-dimensional signal of grid sm(τ=τm,fdr), no-coherence cumulating is directly carried out to the corresponding one-dimensional signal of different big-dipper satellites, i.e.,
5. obtaining three dimensional signal s by 3. and 4. handling to 2. middle total-gridtotal(X,Y,fdr), according to accumulation as a result, The detection of target can be achieved.Target positioning can be achieved at the same time in ground X and Y value according to corresponding to peak value.
From the figure 3, it may be seen that different big-dipper satellites obtain when the present invention is realized well using more big-dipper satellites as radiation source Target echo signal energy effective accumulation.The X-Y diagrammatic cross-section at accumulation result peak value as shown in Figure 4 is it is found that originally It is (1000,0) that target, which can be obtained, in the position of Sea background in invention, is consistent with actual position, realizes positioning well.Cause This, the present invention can be effectively realized using more big-dipper satellites as the detection and positioning of the sea moving target of radiation source.
Those of ordinary skill in the art will understand that the embodiments described herein, which is to help reader, understands this hair Bright principle, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.For ability For the technical staff in domain, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made Any modification, equivalent substitution, improvement and etc. should be included within scope of the presently claimed invention.

Claims (8)

1. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source characterized by comprising
S1, symmetrical Keystone transformation is all made of to the target echo from different big-dipper satellites received, obtains distance and moves Dynamic correction result;
S2, using sub-frame processing and oblique-FFT is gone to the echo-signal after each range migration correction, obtains the product of each echo-signal Tired result;
S3, by the accumulation result of each echo according to difference it is biradical apart from back projection to local coordinate, and in local coordinate Under no-coherence cumulating is carried out to the accumulation result of each echo;
S4, according to the no-coherence cumulating of multiple echoes of step S3 as a result, being detected and being positioned to target.
2. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 1, It is characterized in that, step S1 includes:
S11, system parameter initialization, the parameter of the initialization includes: distance to sample frequency, the target observation time, equivalent The signal pulse repetition time;
S12, using more big-dipper satellites as radiation source, direct wave and target echo are enrolled and are solved respectively with two antennas It adjusts;
S13, relevant treatment is carried out to obtaining each big-dipper satellite as the corresponding direct wave of radiation source and target echo, obtain each north Struggle against satellite as the corresponding target echo of radiation source distance to the compressed result of pulse;
S14, symmetrical Keystone transformation is carried out to the compressed result of pulse to the distance of each target echo, obtains each target and returns The corresponding range migration correction result of wave.
3. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 2, It is characterized in that, step S14 include it is following step by step:
A1, distance is carried out to the corresponding range migration correction result of each target echo to Fast Fourier Transform (FFT), obtain each target Echo is corresponding apart from frequency domain-orientation time-domain signal;
A2, corresponding to each target echo positive second order Keystone transformation is carried out apart from frequency domain-orientation time-domain signal;
A3, corresponding to each target echo negative second order Keystone transformation is carried out apart from frequency domain-orientation time-domain signal;
The result that A4, the step A2 corresponding to each target echo are obtained is multiplied with the obtained result of step A3;
The multiplied result that A5, the step A4 corresponding to each target echo are obtained carries out distance to inverse fast Fourier transform, thus Obtain the corresponding range migration correction result of each target echo.
4. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 3, It is characterized in that, step S2 includes:
S21, the corresponding range migration correction result of some target echo is taken same apart from gate signal, obtains its corresponding orientation Signal;
S22, it goes tiltedly to handle and is divided into multiframe to what the bearing signal was unified;
S23, the coherent accumulation of signal and the no-coherence cumulating of interframe signal in frame are carried out to the bearing signal after going tiltedly;
S24, step S23-S24 is executed to the Distance Remaining gate signal of the corresponding range migration correction result of the target echo, obtained To single echo-signal the biradical distance-Doppler tune frequency domain of difference accumulation result.
5. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 4, It is characterized in that, step S22 include it is following step by step:
B1, the range that doppler frequency rate is set and tune separation values;
B2, one group of linear FM signal is constituted according to tune separation values;
B3, the bearing signal of taking-up is multiplied from the linear FM signal of different frequency modulation rates respectively, obtains 2D signal;
B4, with the time span of a determining frame, the obtained 2D signal of step B3 is divided into several frames.
6. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 5, It is characterized in that, the range of doppler frequency rate described in step B1 includes the true value of doppler frequency rate.
7. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 6, It is characterized in that, step S23 include it is following step by step:
C1, corresponding value of every frame signal each orientation time is directly added, obtains every frame letter under the conditions of different tune frequency values Number zero-frequency value;
C2, the zero-frequency value of every frame signal under the conditions of different tune frequency values is subjected to no-coherence cumulating, obtained in different tune frequency values Under the conditions of signal zero-frequency value.
8. a kind of sea Moving objects detection and location method based on big-dipper satellite radiation source according to claim 7, It is characterized in that, step S3 includes:
S31, building local coordinate, reference axis is respectively position X, position Y and doppler frequency rate;
S32, by Ocean Scenes grid division, obtain the space coordinate of each grid;
S33, for each grid, according to big-dipper satellite and the spatial position of receiving station, obtain the corresponding difference of different big-dipper satellites Biradical distance, and then the corresponding distance of different satellites is obtained to the time;
S34, the corresponding distance of different satellites that step S33 is obtained is brought into no-coherence cumulating knot that step S23 is obtained to the time Fruit obtains the corresponding one-dimensional doppler frequency rate signal of each grid;
S35, the corresponding one-dimensional doppler frequency rate signal of big-dipper satellite each under current grid is directly subjected to no-coherence cumulating;
The no-coherence cumulating of each big-dipper satellite is as a result, obtain in local coordinate system under S36, all grids obtained according to step S35 Three dimensional signal in system.
CN201910278789.XA 2019-04-09 2019-04-09 Sea surface moving target detection and positioning method based on Beidou satellite radiation source Active CN110231616B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910278789.XA CN110231616B (en) 2019-04-09 2019-04-09 Sea surface moving target detection and positioning method based on Beidou satellite radiation source

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910278789.XA CN110231616B (en) 2019-04-09 2019-04-09 Sea surface moving target detection and positioning method based on Beidou satellite radiation source

Publications (2)

Publication Number Publication Date
CN110231616A true CN110231616A (en) 2019-09-13
CN110231616B CN110231616B (en) 2021-01-12

Family

ID=67860767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910278789.XA Active CN110231616B (en) 2019-04-09 2019-04-09 Sea surface moving target detection and positioning method based on Beidou satellite radiation source

Country Status (1)

Country Link
CN (1) CN110231616B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848549A (en) * 2021-09-15 2021-12-28 电子科技大学 Radiation source positioning method based on synthetic aperture technology
CN114114235A (en) * 2021-11-26 2022-03-01 哈尔滨工业大学 Target detection method and system based on single-channel passive radar
CN114488132A (en) * 2022-02-22 2022-05-13 北京卫星信息工程研究所 Target detection method of satellite-borne GNSS-S radar
CN115201803A (en) * 2022-07-13 2022-10-18 北京理工大学 Passive positioning method and system based on satellite over-the-top time measurement

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011174875A (en) * 2010-02-25 2011-09-08 Mitsubishi Electric Corp Passive radar device
CN104569968A (en) * 2014-12-10 2015-04-29 合肥工业大学 Passive imaging system based on GPS multi-star irradiation and ground single station receiving
CN104678417A (en) * 2015-02-11 2015-06-03 中国地质大学(北京) Target positioning method based on reflection signals of navigation satellite
CN105137420A (en) * 2015-09-09 2015-12-09 电子科技大学 Multiframe accumulated noncoherent MIMO radar detection and positioning method
CN109100762A (en) * 2018-07-10 2018-12-28 电子科技大学 Bistatic radar sea moving target detecting method based on big-dipper satellite radiation source

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011174875A (en) * 2010-02-25 2011-09-08 Mitsubishi Electric Corp Passive radar device
CN104569968A (en) * 2014-12-10 2015-04-29 合肥工业大学 Passive imaging system based on GPS multi-star irradiation and ground single station receiving
CN104678417A (en) * 2015-02-11 2015-06-03 中国地质大学(北京) Target positioning method based on reflection signals of navigation satellite
CN105137420A (en) * 2015-09-09 2015-12-09 电子科技大学 Multiframe accumulated noncoherent MIMO radar detection and positioning method
CN109100762A (en) * 2018-07-10 2018-12-28 电子科技大学 Bistatic radar sea moving target detecting method based on big-dipper satellite radiation source

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PENGHUI HUANG: "A Fast SAR Imaging Method for Ground Moving Target Using a Second-Order WVD Transform", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
ZHONGYU LI: "Satellite-based PCL for moving target detection using keystone transform and FrFT", 《2016 CIE INTERNATIONAL CONFERENCE ON RADAR (RADAR)》 *
王俊林: "外辐射源雷达***中***处理增益分析", 《数据采集与处理》 *
范梅梅: "基于北斗卫星信号的无源雷达可行性研究", 《信号处理》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848549A (en) * 2021-09-15 2021-12-28 电子科技大学 Radiation source positioning method based on synthetic aperture technology
CN113848549B (en) * 2021-09-15 2023-06-23 电子科技大学 Radiation source positioning method based on synthetic aperture technology
CN114114235A (en) * 2021-11-26 2022-03-01 哈尔滨工业大学 Target detection method and system based on single-channel passive radar
CN114488132A (en) * 2022-02-22 2022-05-13 北京卫星信息工程研究所 Target detection method of satellite-borne GNSS-S radar
CN114488132B (en) * 2022-02-22 2023-01-20 北京卫星信息工程研究所 Target detection method of satellite-borne GNSS-S radar
CN115201803A (en) * 2022-07-13 2022-10-18 北京理工大学 Passive positioning method and system based on satellite over-the-top time measurement
CN115201803B (en) * 2022-07-13 2024-04-16 北京理工大学 Passive positioning method and system based on satellite overhead moment measurement

Also Published As

Publication number Publication date
CN110231616B (en) 2021-01-12

Similar Documents

Publication Publication Date Title
CN110231616A (en) Sea Moving objects detection and location method based on big-dipper satellite radiation source
CN105445730B (en) A kind of Sea Current inverting Spaceborne SAR System and its method based on angle diversity
CN109100762B (en) Bistatic radar sea surface moving target detection method based on Beidou satellite radiation source
CN105785327B (en) Frequency diversity array synthetic aperture radar high-resolution wide swath imaging method
US9442188B2 (en) Negative pseudo-range processing with multi-static FMCW radars
US4172255A (en) HF coastal current mapping radar system
CN105866750B (en) A kind of method of multiple GPS satellite weak echo signal detections
CN108398676A (en) A kind of external illuminators-based radar Dim moving target detection method
CN107966688A (en) Broadband radar target velocity solution blur method based on phase interference techniques
CN109507668B (en) Double-base imaging method based on navigation satellite signals
Huang et al. Radar high speed small target detection based on keystone transform and linear canonical transform
Li et al. Passive radar array with low-power satellite illuminators based on fractional Fourier transform
CN103018719B (en) Generation method of OTH (Over-The-Horizon) radar transmitting waveform
EP3771921B1 (en) Method for processing non-interrupted phase synchronization signal of bistatic sar based on coded signal
CN104678386B (en) Method for detecting target by utilizing relevant power of GNSS(global navigation satellite system) sea surface reflection signal
CN109814100A (en) SAR Ground moving target imaging method based on sub-aperture parameter Estimation
CN103592647A (en) Array three-dimensional SAR data acquisition method
CN101813772A (en) Array beamforming method by quickly expanding and dragging broadband frequency domain
Fang et al. Improved passive SAR imaging with DVB-T transmissions
CN110133646B (en) NLCS imaging-based multi-channel two-pulse clutter cancellation method for bistatic forward-looking SAR
Walterscheid et al. Challenges and first results of an airborne passive SAR experiment using a DVB-T transmitter
CN108549064A (en) External sort algorithm moving-target detection method based on Doppler frequency fuzzy compensation in arteries and veins
CN108919206B (en) External radiation source radar polarization filtering method based on subcarrier processing
Gurgel et al. Remote sensing of surface currents and waves by the HF radar WERA
CN113945901A (en) Engineering implementation method for suppressing radio frequency interference through compressed sensing

Legal Events

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