CN112596056A - Passive synthetic aperture radiation source positioning method based on multi-view processing - Google Patents

Passive synthetic aperture radiation source positioning method based on multi-view processing Download PDF

Info

Publication number
CN112596056A
CN112596056A CN202110213207.7A CN202110213207A CN112596056A CN 112596056 A CN112596056 A CN 112596056A CN 202110213207 A CN202110213207 A CN 202110213207A CN 112596056 A CN112596056 A CN 112596056A
Authority
CN
China
Prior art keywords
radiation source
positioning
signal
distance
target
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
CN202110213207.7A
Other languages
Chinese (zh)
Other versions
CN112596056B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202110213207.7A priority Critical patent/CN112596056B/en
Publication of CN112596056A publication Critical patent/CN112596056A/en
Application granted granted Critical
Publication of CN112596056B publication Critical patent/CN112596056B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9041—Squint mode

Landscapes

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

Abstract

The invention provides a passive synthetic aperture radiation source positioning method based on multi-view processing, which comprises the steps of reading/simulating a radiation source target signal received by a radar generated by an instrument, and carrying out down-conversion and demodulation on the received signal to obtain a Doppler signal of a radiation source target; aiming at the determined positioning area, according to the relation between the Doppler center frequency, the modulation frequency and the squint angle, matching filters under different squint angles are designed to obtain radiation source positioning images under different squint angles; and performing incoherent accumulation on the radiation source positioning images at different oblique angles to obtain a radiation source positioning image of the positioning area, and splicing the radiation source positioning images in different areas to obtain a high-coverage positioning image. According to the invention, the data volume is reduced through downsampling, the positioning efficiency is improved, the signal-to-noise ratio of the positioning image is improved through incoherent accumulation of the positioning images under different squint angles, and by adopting the technical scheme of the invention, the high-precision, high-resolution and high-sensitivity rapid positioning of the radiation source target is realized.

Description

Passive synthetic aperture radiation source positioning method based on multi-view processing
Technical Field
The invention relates to a radiation source positioning technology, in particular to a passive synthetic aperture radiation source positioning method based on multi-view processing.
Background
The position of the radiation source is important characteristic information of the electromagnetic radiation source, has a relatively stable state, is closely related to battlefield situation, mission planning and combat action, and is an important basis for distinguishing different radiation sources from one another. The traditional radiation source positioning method mainly comprises the following types: frequency measurement positioning, lateral positioning and time difference positioning. However, the method has the problems that high coverage area, high precision and high sensitivity cannot be simultaneously achieved. To address this problem, researchers have proposed passive synthetic aperture location methods. The passive synthetic aperture positioning method realizes the high-precision and high-sensitivity positioning of the target of the radiation source through the matching filtering in the azimuth direction. However, the problem that the linear frequency modulation matched filter is not applicable any more in a longer synthetic aperture time exists, and the long synthetic aperture time is long, the azimuth bandwidth is large, the signal processing time is long, and the positioning efficiency is low.
Disclosure of Invention
The purpose of the invention is: the problem that a linear matching filter is not applicable any more and the positioning efficiency is low within a long synthetic aperture time is solved.
The invention adopts the following scheme for solving the technical problems:
a passive synthetic aperture radiation source positioning method based on multi-view processing comprises the following steps:
step 101: the method comprises the steps that an instrument reads/simulates a radiation source target signal received by a radar, and the received signal is subjected to down-conversion and demodulation to obtain a Doppler signal of a radiation source target;
step 102: aiming at the determined positioning area, according to the relation between the Doppler center frequency, the modulation frequency and the squint angle, matching filters under different squint angles are designed to obtain radiation source positioning images under different squint angles;
step 103: and performing incoherent accumulation on the radiation source positioning images at different oblique angles to obtain a radiation source positioning image of the positioning area, and splicing the radiation source positioning images in different areas to obtain a high-coverage positioning image.
The specific process of the step 101 is as follows:
setting relevant parameters of a radiation source, specifically comprising: the modulation mode and the carrier frequency of each radiation source signal are the same, and the modulation mode of each radiation source signal is binary phase shift keying BPSK, and the carrier frequency is
Figure 929455DEST_PATH_IMAGE001
Of 1 at
Figure 293266DEST_PATH_IMAGE002
A radiation source signal is
Figure 745107DEST_PATH_IMAGE003
，
Figure 874737DEST_PATH_IMAGE004
，
Figure 610481DEST_PATH_IMAGE005
The time is represented by the time of day,
Figure 636206DEST_PATH_IMAGE006
is as follows
Figure 90190DEST_PATH_IMAGE002
A baseband symbol signal of the respective radiation source signal,
Figure 944882DEST_PATH_IMAGE007
the total number of the radiation source targets;
setting up scene and airborne radar parameter specifically includes: the target of the radiation source radiates electromagnetic signals to the periphery on the ground surface and has the height of
Figure 285865DEST_PATH_IMAGE009
Speed of airborne radar
Figure 13649DEST_PATH_IMAGE010
Fly at a constant speed, and
Figure 954929DEST_PATH_IMAGE011
the target distance of each radiation source is
Figure 629624DEST_PATH_IMAGE012
Of 1 at
Figure 339960DEST_PATH_IMAGE002
The target azimuth distance of each radiation source is
Figure 441908DEST_PATH_IMAGE013
；
Thereby obtaining signals received by the airborne radar
Figure 667222DEST_PATH_IMAGE014
：
Figure 614450DEST_PATH_IMAGE015
Wherein
Figure 444871DEST_PATH_IMAGE016
Is shown as
Figure 248879DEST_PATH_IMAGE011
The instantaneous distance of each source target to the radar,
Figure 915484DEST_PATH_IMAGE017
is shown as
Figure 452688DEST_PATH_IMAGE011
The received signal strength of the individual radiation source targets,
Figure 622769DEST_PATH_IMAGE018
is zero mean and variance of
Figure 378104DEST_PATH_IMAGE019
The white noise of the gaussian is generated,
Figure 846DEST_PATH_IMAGE020
represents the speed of light;
receiving signal according to signal carrier frequency
Figure 617773DEST_PATH_IMAGE021
Down conversion to obtain signal
Figure 688366DEST_PATH_IMAGE022
：
Figure 303018DEST_PATH_IMAGE023
Wherein
Figure 209794DEST_PATH_IMAGE024
；
For down-converted signal
Figure 614099DEST_PATH_IMAGE025
Square demodulation to obtain Doppler signal of radiation source target
Figure 555510DEST_PATH_IMAGE026
：
Figure 137802DEST_PATH_IMAGE027
Wherein
Figure 531874DEST_PATH_IMAGE028
Is a normal plural number;
for Doppler signal
Figure 677553DEST_PATH_IMAGE029
Cut off and sample to obtain discrete signal
Figure 207892DEST_PATH_IMAGE030
：
Figure 226663DEST_PATH_IMAGE031
Wherein
Figure 357299DEST_PATH_IMAGE032
，
Figure 57402DEST_PATH_IMAGE033
Which is indicative of the sampling frequency, is,
Figure 707826DEST_PATH_IMAGE034
intercepting the signal in a time range of
Figure 615608DEST_PATH_IMAGE035
，
Figure 984273DEST_PATH_IMAGE036
Meaning that the rounding is done down,
Figure 19225DEST_PATH_IMAGE037
，
Figure 55314DEST_PATH_IMAGE038
to represent
Figure 133997DEST_PATH_IMAGE039
A matrix of complex numbers is maintained.
The specific process of step 102 is as follows:
determining a location area based on beam pointing directions and spatial coordinate positions of receiving antennas
Figure 193220DEST_PATH_IMAGE040
：
Distance to distance range
Figure 297442DEST_PATH_IMAGE041
Azimuth distance range
Figure 923725DEST_PATH_IMAGE042
Wherein, in the step (A),
Figure 392884DEST_PATH_IMAGE043
is expressed by
Figure 470561DEST_PATH_IMAGE044
As a starting point, the method comprises the following steps of,
Figure 627742DEST_PATH_IMAGE045
is an end point, interval
Figure 310527DEST_PATH_IMAGE046
The sampling is carried out and the sampling is carried out,
Figure 747325DEST_PATH_IMAGE047
to represent
Figure 827145DEST_PATH_IMAGE048
The matrix of real numbers is then maintained,
Figure 210853DEST_PATH_IMAGE049
the total distance is counted to the distance unit;
beam pointing and location area based on receiving antenna
Figure 528571DEST_PATH_IMAGE040
Determining squint angle range for multi-view sub-images
Figure 136269DEST_PATH_IMAGE050
Wherein, in the step (A),
Figure 109911DEST_PATH_IMAGE051
total number of squint angles;
for positioning area
Figure 94047DEST_PATH_IMAGE040
Sum matched filter duration
Figure 548162DEST_PATH_IMAGE052
Obtaining an oblique view angle
Figure 326762DEST_PATH_IMAGE053
The effective time range of the signals to be processed is as follows:
Figure 115596DEST_PATH_IMAGE054
wherein
Figure 841106DEST_PATH_IMAGE055
Is the speed of the radar, and,
Figure 149728DEST_PATH_IMAGE056
，
Figure 614076DEST_PATH_IMAGE057
；
order to
Figure 578621DEST_PATH_IMAGE058
，
Figure 170139DEST_PATH_IMAGE059
If, if
Figure 582535DEST_PATH_IMAGE060
And is
Figure 906200DEST_PATH_IMAGE061
Then, then
Figure 420358DEST_PATH_IMAGE062
，
Figure 805115DEST_PATH_IMAGE063
，
Figure 822750DEST_PATH_IMAGE064
，
Figure 379633DEST_PATH_IMAGE065
(ii) a If it is not
Figure 568038DEST_PATH_IMAGE066
Or
Figure 235779DEST_PATH_IMAGE067
The truncated Doppler discrete signal obtained in step 101
Figure 107920DEST_PATH_IMAGE068
Oblique angle cannot be obtained
Figure 101284DEST_PATH_IMAGE069
Lower positioning area
Figure 573723DEST_PATH_IMAGE070
The positioning image of (1) needs to skip the oblique view
Figure 779576DEST_PATH_IMAGE071
The next oblique view angle is performed in the next positioning process
Figure 975065DEST_PATH_IMAGE072
The positioning process of (2);
for positioning area
Figure 123019DEST_PATH_IMAGE070
Considering the influence of square demodulation, based on the demodulated Doppler center frequency, instantaneous modulation frequency and squint angle
Figure 833486DEST_PATH_IMAGE073
In relation to (2)
Figure 46292DEST_PATH_IMAGE074
，
Figure 407872DEST_PATH_IMAGE075
Obtaining an oblique view
Figure 477460DEST_PATH_IMAGE076
Center frequency of lower Doppler
Figure 675223DEST_PATH_IMAGE077
And instantaneous frequency modulation rate
Figure 675409DEST_PATH_IMAGE078
；
For the signal
Figure 642228DEST_PATH_IMAGE079
Are mixed to obtain
Figure 882716DEST_PATH_IMAGE080
；
For the signal
Figure 20305DEST_PATH_IMAGE081
Low-pass filtering to obtain
Figure 371652DEST_PATH_IMAGE082
Denotes a cut-off frequency of
Figure 192978DEST_PATH_IMAGE083
The low-pass filter of (1);
for the signal
Figure 869947DEST_PATH_IMAGE084
Down-sampling to obtain
Figure 25990DEST_PATH_IMAGE085
，
Figure 915449DEST_PATH_IMAGE086
，
Figure 60122DEST_PATH_IMAGE087
Is the multiple of the down-sampling,
Figure 620242DEST_PATH_IMAGE088
，
Figure 14315DEST_PATH_IMAGE089
，
Figure 973043DEST_PATH_IMAGE090
，
Figure 768961DEST_PATH_IMAGE091
represents rounding up;
according to miningSample frequency
Figure 974683DEST_PATH_IMAGE092
And downsampling multiple
Figure 856052DEST_PATH_IMAGE094
Determining a location area
Figure 87313DEST_PATH_IMAGE095
The division of the azimuth distance is:
Figure 252584DEST_PATH_IMAGE096
wherein
Figure 176678DEST_PATH_IMAGE097
The total number of azimuth distance units;
design squint angle
Figure 748604DEST_PATH_IMAGE098
The following matched filter:
Figure 32824DEST_PATH_IMAGE099
wherein
Figure 803334DEST_PATH_IMAGE100
；
For distance in the radial direction
Figure 898329DEST_PATH_IMAGE101
Using oblique angles of view
Figure 754289DEST_PATH_IMAGE102
Lower corresponding matched filter
Figure 779883DEST_PATH_IMAGE103
For the Doppler signals obtained by processing
Figure 404900DEST_PATH_IMAGE104
Performing matched filtering to obtain:
Figure 670796DEST_PATH_IMAGE105
wherein
Figure 263320DEST_PATH_IMAGE106
For the convolution calling function in Matlab software,
Figure 905654DEST_PATH_IMAGE107
；
will be provided with
Figure 385177DEST_PATH_IMAGE108
As a result of the positioning of the azimuth distance range:
Figure 274504DEST_PATH_IMAGE109
wherein
Figure 839478DEST_PATH_IMAGE110
，
Figure 754344DEST_PATH_IMAGE111
The expression is taken to be the minimum value,
Figure 619532DEST_PATH_IMAGE112
representing an oblique angle of view
Figure 679761DEST_PATH_IMAGE113
Distance in the direction of distance
Figure 732030DEST_PATH_IMAGE114
The lower azimuth positioning result;
according to different distance directions
Figure 716167DEST_PATH_IMAGE114
Downward azimuth distance positioning result
Figure 425409DEST_PATH_IMAGE115
Obtain an oblique view angle
Figure 204009DEST_PATH_IMAGE113
Lower positioning area
Figure 946837DEST_PATH_IMAGE095
Positioning image of
Figure 718353DEST_PATH_IMAGE116
。
The specific process of step 103 is as follows:
the different squint angles obtained in step 102
Figure 26975DEST_PATH_IMAGE117
The radiation source positioning image of the positioning area is obtained by non-coherent accumulation of the lower radiation source positioning image
Figure 242055DEST_PATH_IMAGE118
：
Figure 3338DEST_PATH_IMAGE119
；
Wherein
Figure 578545DEST_PATH_IMAGE120
；
According to the location area
Figure 944935DEST_PATH_IMAGE095
Range of (1)
Figure 65338DEST_PATH_IMAGE121
，
Figure 828763DEST_PATH_IMAGE122
And positioning the image
Figure 958393DEST_PATH_IMAGE118
Drawing a positioning area
Figure 976028DEST_PATH_IMAGE095
And judging whether a radiation source target exists or not according to the positioning image, and if the target exists, finding the target and determining the position of the radiation source target.
The invention has the beneficial effects that: the invention provides a passive synthetic aperture radiation source positioning method based on multi-view processing, because the signal duration of matched filtering under different squint angles is shorter, the signal bandwidth is narrower, the positioning efficiency can be improved by reducing the data volume through downsampling, and the signal-to-noise ratio of the positioning image can be improved by incoherent accumulation of the positioning images under different squint angles, thereby realizing the rapid positioning of the radiation source target with high precision, high resolution and high sensitivity.
Drawings
FIG. 1 is a schematic illustration of the positioning principle of a radiation source in an embodiment of the present invention;
FIG. 2 is a schematic flow chart illustrating the positioning of a radiation source in an embodiment of the present invention;
FIG. 3 is a schematic view of a geometric model of the positioning of a radiation source in an embodiment of the present invention;
FIG. 4 is a schematic diagram illustrating effective processing time of signals at different squint angles according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of an azimuth distance matching process according to an embodiment of the present invention;
FIG. 6 is a schematic diagram illustrating a geometric relationship between positioning images at different oblique angles according to an embodiment of the present disclosure;
FIG. 7a is a simulation result of an embodiment of the present invention;
FIG. 7b is a top view of a positioning image according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
Fig. 1 is a schematic diagram of a positioning principle of a radiation source in an embodiment of the present invention, as shown in fig. 1, in a long synthetic aperture time, a doppler history of a radar relative to a target exhibits a nonlinear characteristic, and a chirp matched filter is no longer applicable.
Fig. 2 is a schematic view of a positioning process of a radiation source in an embodiment of the present invention, and as shown in fig. 2, the technical scheme of the present invention includes the following steps:
step 101: the instrument reads/simulates a radiation source target signal received by the radar, and the Doppler signal of the radiation source target is obtained by down-conversion and demodulation of the received signal.
Specifically, FIG. 3 is a schematic diagram of a geometric model of the positioning of the radiation source according to an embodiment of the present invention, as shown in FIG. 3, the target of the radiation source is in the detection region
Figure 532911DEST_PATH_IMAGE123
Electromagnetic signals are radiated inwards and all around, the modulation modes and carrier frequencies of all radiation source signals are the same, the modulation types of all radiation source signals are Binary Phase Shift Keying (BPSK), and the carrier frequencies are all Binary Phase Shift Keying (BPSK)
Figure 721316DEST_PATH_IMAGE124
Of 1 at
Figure 123478DEST_PATH_IMAGE125
A radiation source signal is
Figure 261199DEST_PATH_IMAGE126
，
Figure 238251DEST_PATH_IMAGE127
，
Figure 461422DEST_PATH_IMAGE128
The time is represented by the time of day,
Figure 932854DEST_PATH_IMAGE129
is as follows
Figure 925081DEST_PATH_IMAGE125
A baseband symbol signal of the respective radiation source signal,
Figure 10718DEST_PATH_IMAGE130
the total number of the radiation source targets. Has a height of
Figure 721185DEST_PATH_IMAGE009
Speed of airborne radar
Figure 730729DEST_PATH_IMAGE131
Flying at constant speed, the signal received by the target radiation of the radiation source is
Figure 92309DEST_PATH_IMAGE132
，
Figure 161896DEST_PATH_IMAGE133
Is shown as
Figure 94080DEST_PATH_IMAGE125
The received signal strength of the individual radiation source targets,
Figure 377424DEST_PATH_IMAGE134
is shown as
Figure 609822DEST_PATH_IMAGE125
The instantaneous distance of each source target to the radar,
Figure 584731DEST_PATH_IMAGE135
is as follows
Figure 269790DEST_PATH_IMAGE136
The distance of the radiation sources to the distance,
Figure 870405DEST_PATH_IMAGE137
is as follows
Figure 894993DEST_PATH_IMAGE125
The azimuthal distance of each of the radiation sources,
Figure 40803DEST_PATH_IMAGE138
is zero mean and variance of
Figure 196847DEST_PATH_IMAGE139
The white noise of the gaussian is generated,
Figure 617464DEST_PATH_IMAGE140
indicating the speed of light.
Receiving signal according to signal carrier frequency
Figure 293296DEST_PATH_IMAGE141
Down conversion to obtain signal
Figure 875587DEST_PATH_IMAGE142
：
Figure 722189DEST_PATH_IMAGE143
Wherein
Figure 415339DEST_PATH_IMAGE144
；
For down-converted signal
Figure 211256DEST_PATH_IMAGE145
Square demodulation to obtain Doppler signal of radiation source target
Figure 479295DEST_PATH_IMAGE146
：
Figure 95084DEST_PATH_IMAGE147
Wherein
Figure 591925DEST_PATH_IMAGE148
Is a normal plural number;
for Doppler signal
Figure 694879DEST_PATH_IMAGE149
Cut off and sample to obtain discrete signal
Figure 618973DEST_PATH_IMAGE150
：
Figure 987637DEST_PATH_IMAGE151
Wherein
Figure 22589DEST_PATH_IMAGE152
，
Figure 42367DEST_PATH_IMAGE153
Which is indicative of the sampling frequency, is,
Figure 137362DEST_PATH_IMAGE154
intercepting the signal in a time range of
Figure 727743DEST_PATH_IMAGE155
，
Figure 18916DEST_PATH_IMAGE156
Meaning that the rounding is done down,
Figure 643932DEST_PATH_IMAGE157
，
Figure 909829DEST_PATH_IMAGE158
to represent
Figure 987506DEST_PATH_IMAGE159
A matrix of complex numbers is maintained.
Step 102: and aiming at the determined positioning area, designing matched filters under different oblique angles according to the relation between the Doppler center frequency, the modulation frequency and the oblique angles, and obtaining radiation source positioning images under different oblique angles.
Specifically, as shown in fig. 3, the positioning area is determined according to the beam direction and the spatial coordinate position of the receiving antenna
Figure 884967DEST_PATH_IMAGE161
：
Distance to distance range
Figure 364490DEST_PATH_IMAGE162
Azimuth distance range
Figure 801287DEST_PATH_IMAGE163
Wherein, in the step (A),
Figure 818791DEST_PATH_IMAGE164
is expressed by
Figure 999236DEST_PATH_IMAGE165
As a starting point, the method comprises the following steps of,
Figure 598845DEST_PATH_IMAGE166
is an end point, interval
Figure 472123DEST_PATH_IMAGE167
The sampling is carried out and the sampling is carried out,
Figure 508081DEST_PATH_IMAGE168
to represent
Figure 492218DEST_PATH_IMAGE169
The matrix of real numbers is then maintained,
Figure 149595DEST_PATH_IMAGE170
the total distance is counted to the distance unit;
beam pointing and location area based on receiving antenna
Figure 177463DEST_PATH_IMAGE171
Determining squint angle range for multi-view sub-images
Figure 451449DEST_PATH_IMAGE172
Wherein, in the step (A),
Figure 973697DEST_PATH_IMAGE173
total number of squint angles;
FIG. 4 is a schematic diagram of effective signal processing time at different squint angles according to an embodiment of the present invention, as shown in FIG. 4, for a positioning region
Figure 547898DEST_PATH_IMAGE171
Sum matched filter duration
Figure 746667DEST_PATH_IMAGE174
Obtaining an oblique view angle
Figure 976791DEST_PATH_IMAGE175
The effective time range of the signals to be processed is as follows:
Figure 37151DEST_PATH_IMAGE176
wherein
Figure 715126DEST_PATH_IMAGE177
Is the speed of the radar, and,
Figure 835529DEST_PATH_IMAGE178
，
Figure 349687DEST_PATH_IMAGE179
；
order to
Figure 213738DEST_PATH_IMAGE180
，
Figure 746219DEST_PATH_IMAGE181
If, if
Figure 240785DEST_PATH_IMAGE182
And is
Figure 976660DEST_PATH_IMAGE183
Then, then
Figure 893670DEST_PATH_IMAGE184
，
Figure 296969DEST_PATH_IMAGE185
，
Figure 759174DEST_PATH_IMAGE186
，
Figure 429016DEST_PATH_IMAGE187
(ii) a If it is not
Figure 900449DEST_PATH_IMAGE188
Or
Figure 892675DEST_PATH_IMAGE189
The truncated Doppler discrete signal obtained in step 101
Figure 791361DEST_PATH_IMAGE190
Oblique angle cannot be obtained
Figure 751096DEST_PATH_IMAGE191
Lower positioning area
Figure 963902DEST_PATH_IMAGE192
The positioning image of (1) needs to skip the oblique view
Figure 810636DEST_PATH_IMAGE193
The next oblique view angle is performed in the next positioning process
Figure 129490DEST_PATH_IMAGE194
The positioning process of (2);
for positioning area
Figure 61674DEST_PATH_IMAGE192
Considering the influence of square demodulation, based on the demodulated Doppler center frequency, instantaneous modulation frequency and squint angle
Figure 140489DEST_PATH_IMAGE195
In relation to (2)
Figure 841729DEST_PATH_IMAGE196
，
Figure 534747DEST_PATH_IMAGE197
Obtaining an oblique view
Figure 954227DEST_PATH_IMAGE198
Center frequency of lower Doppler
Figure 305574DEST_PATH_IMAGE199
And instantaneous frequency modulation rate
Figure 641746DEST_PATH_IMAGE200
；
For the signal
Figure 53136DEST_PATH_IMAGE201
Are mixed to obtain
Figure 959912DEST_PATH_IMAGE202
；
For the signal
Figure 114950DEST_PATH_IMAGE203
Low-pass filtering to obtain
Figure 243312DEST_PATH_IMAGE204
Denotes a cut-off frequency of
Figure 825603DEST_PATH_IMAGE205
The low-pass filter of (1);
for the signal
Figure 219675DEST_PATH_IMAGE206
Down-sampling to obtain
Figure 427671DEST_PATH_IMAGE207
，
Figure 223589DEST_PATH_IMAGE208
，
Figure 711202DEST_PATH_IMAGE209
Is the multiple of the down-sampling,
Figure 45101DEST_PATH_IMAGE210
，
Figure 541941DEST_PATH_IMAGE211
，
Figure 192365DEST_PATH_IMAGE212
，
Figure 382038DEST_PATH_IMAGE213
represents rounding up;
according to the sampling frequency
Figure 740250DEST_PATH_IMAGE214
And downsampling multiple
Figure 775202DEST_PATH_IMAGE215
Determining a location area
Figure 748975DEST_PATH_IMAGE216
The division of the azimuth distance is:
Figure 93237DEST_PATH_IMAGE217
wherein
Figure 683618DEST_PATH_IMAGE218
Is the total number of azimuth distance elements.
Design squint angle
Figure 787841DEST_PATH_IMAGE219
The following matched filter:
Figure 412857DEST_PATH_IMAGE220
wherein
Figure 928021DEST_PATH_IMAGE221
；
FIG. 5 is a schematic diagram of an azimuth distance matching process in the embodiment of the present invention, as shown in FIG. 5, for distance distances
Figure 208961DEST_PATH_IMAGE222
Using oblique angles of view
Figure 585715DEST_PATH_IMAGE223
Lower corresponding matched filter
Figure 783347DEST_PATH_IMAGE224
For the Doppler signals obtained by processing
Figure 220145DEST_PATH_IMAGE225
Performing matched filtering to obtain:
Figure 50698DEST_PATH_IMAGE226
wherein
Figure 480411DEST_PATH_IMAGE227
For the convolution calling function in Matlab software,
Figure 814440DEST_PATH_IMAGE228
；
will be provided with
Figure 625401DEST_PATH_IMAGE229
As a result of the positioning of the azimuth distance range:
Figure 926939DEST_PATH_IMAGE230
wherein
Figure 911075DEST_PATH_IMAGE231
，
Figure 365190DEST_PATH_IMAGE232
The expression is taken to be the minimum value,
Figure 143790DEST_PATH_IMAGE233
representing an oblique angle of view
Figure 932624DEST_PATH_IMAGE234
Distance in the direction of distance
Figure 454872DEST_PATH_IMAGE235
The lower azimuth positioning result;
according to different distance directions
Figure 966756DEST_PATH_IMAGE235
Downward azimuth distance positioning result
Figure 165525DEST_PATH_IMAGE236
Obtain an oblique view angle
Figure 192387DEST_PATH_IMAGE234
Lower positioning area
Figure 518326DEST_PATH_IMAGE216
Positioning image of
Figure 681454DEST_PATH_IMAGE237
。
Step 103: and performing incoherent accumulation on the radiation source positioning images at different oblique angles to obtain a radiation source positioning image of the positioning area, and splicing the radiation source positioning images in different areas to obtain a high-coverage positioning image.
Specifically, fig. 6 is a schematic view of a geometric relationship between positioning images at different oblique angles in the embodiment of the present invention, and as shown in fig. 6, the positioning images at different oblique angles obtained in step 102 are processed
Figure 537544DEST_PATH_IMAGE238
The radiation source positioning image of the positioning area is obtained by non-coherent accumulation of the lower radiation source positioning image
Figure 786123DEST_PATH_IMAGE239
：
Figure 915753DEST_PATH_IMAGE240
；
Wherein
Figure 651496DEST_PATH_IMAGE241
；
According to the location area
Figure 942800DEST_PATH_IMAGE216
Range of (1)
Figure 678675DEST_PATH_IMAGE242
，
Figure 861264DEST_PATH_IMAGE243
And positioning the image
Figure 998984DEST_PATH_IMAGE239
Drawing a positioning area
Figure 664452DEST_PATH_IMAGE216
And judging whether a radiation source target exists or not according to the positioning image, and if the target exists, finding the target and determining the position of the radiation source target.
The treatment effect of the present invention is further shown below with reference to specific examples.
The effect display is displayed by a simulation experiment, and the simulation parameters set by the embodiment comprise:
signal carrier frequency: 7.2445 GHz; code rate: 1000 bound/s, airborne radar speed 250 m/s; sampling time range: 0-80 s; sampling frequency: 50 kHz; down-sampling multiple: 100, respectively; an imaging area: distance to distance range: 9900 m: 1 m: 10100m, azimuth distance range: 9900 m: 0.5 m: 10100 m; radiation source target 1: distance 10000m, azimuth distance 10000m, and received signal strength 1; radiation source target 2: distance 10002m, azimuth distance 10002m, received signal strength 1; noise variance: 10; the signal-to-noise ratio for each radiation source is: -10 dB; squint angle range for positioning: -40, place: 0.1 part: 401; duration of the matched filter: 1 s; the low pass filter has a cut-off frequency of 300Hz, an order of 30 and is constructed by using a hamming window.
The simulation generated signal data length is 4000001, and the positioning process takes about 23 s.
Fig. 7a and 7b are simulation results of the embodiment of the present invention, fig. 7a is a positioning result, fig. 7b is a top view of a positioning image, and it can be seen from the positioning result that there are two radiation source targets: target 1: a distance 10000m in a direction and a distance 10000m in an azimuth direction; target 2: the range distance is 10002m, and the azimuth distance is 10002m, which are consistent with the real target. Therefore, the invention can realize the rapid positioning of the target of the radiation source with high precision, high resolution and high sensitivity without considering the influence of fixed frequency deviation brought by down-conversion.
The above detailed description further illustrates the object, technical solutions and advantages of the present invention, and it should be understood that the above embodiments are only examples of the present invention and are not intended to limit the scope of the present invention, and modifications and equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims (4)

1. A method for positioning a passive synthetic aperture radiation source based on multi-view processing, the method comprising the steps of:
step 101: the method comprises the steps that an instrument reads/simulates a radiation source target signal received by a radar, and the received signal is subjected to down-conversion and demodulation to obtain a Doppler signal of a radiation source target;
step 102: aiming at the determined positioning area, according to the relation between the Doppler center frequency, the modulation frequency and the squint angle, matching filters under different squint angles are designed to obtain radiation source positioning images under different squint angles;
step 103: and performing incoherent accumulation on the radiation source positioning images at different oblique angles to obtain a radiation source positioning image of the positioning area, and splicing the radiation source positioning images in different areas to obtain a high-coverage positioning image.
2. The method according to claim 1, wherein the specific process of step 101 is as follows:
setting relevant parameters of a radiation source, specifically comprising: the modulation mode and the carrier frequency of each radiation source signal are the same, and the modulation mode of each radiation source signal is binary phase shift keying BPSK, and the carrier frequency is
Figure 626555DEST_PATH_IMAGE001
Of 1 at
Figure 789552DEST_PATH_IMAGE002
A radiation source signal is
Figure 533517DEST_PATH_IMAGE003
，
Figure 911409DEST_PATH_IMAGE004
，
Figure 188454DEST_PATH_IMAGE005
The time is represented by the time of day,
Figure 651797DEST_PATH_IMAGE006
is as follows
Figure 199453DEST_PATH_IMAGE007
A baseband symbol signal of the respective radiation source signal,
Figure 618802DEST_PATH_IMAGE008
the total number of the radiation source targets;
setting up scene and airborne radar parameter specifically includes: the target of the radiation source radiates electromagnetic signals to the periphery on the ground surface and has the height of
Figure 124869DEST_PATH_IMAGE009
Speed of airborne radar
Figure 544349DEST_PATH_IMAGE010
Fly at a constant speed, and
Figure 348226DEST_PATH_IMAGE011
the target distance of each radiation source is
Figure 169552DEST_PATH_IMAGE012
Of 1 at
Figure 846521DEST_PATH_IMAGE011
The target azimuth distance of each radiation source is
Figure 2564DEST_PATH_IMAGE013
；
Thereby obtaining signals received by the airborne radar
Figure 95285DEST_PATH_IMAGE014
：
Figure 302276DEST_PATH_IMAGE015
Wherein
Figure 340026DEST_PATH_IMAGE016
Is shown as
Figure 734098DEST_PATH_IMAGE011
The instantaneous distance of each source target to the radar,
Figure 427248DEST_PATH_IMAGE017
is shown as
Figure 472433DEST_PATH_IMAGE011
The received signal strength of the individual radiation source targets,
Figure 428888DEST_PATH_IMAGE018
is zero mean and variance of
Figure 310256DEST_PATH_IMAGE019
The white noise of the gaussian is generated,
Figure 56364DEST_PATH_IMAGE020
represents the speed of light;
receiving signal according to signal carrier frequency
Figure 972368DEST_PATH_IMAGE021
Down conversion to obtain signal
Figure 630882DEST_PATH_IMAGE022
：
Figure 452076DEST_PATH_IMAGE023
Wherein
Figure 752608DEST_PATH_IMAGE024
；
For down-converted signal
Figure 523118DEST_PATH_IMAGE025
Square demodulation to obtain Doppler signal of radiation source target
Figure 352533DEST_PATH_IMAGE026
：
Figure 658094DEST_PATH_IMAGE027
Wherein
Figure 496737DEST_PATH_IMAGE028
Is a normal plural number;
for Doppler signal
Figure 856174DEST_PATH_IMAGE029
Cut off and sample to obtain discrete signal
Figure 636917DEST_PATH_IMAGE030
：
Figure 714595DEST_PATH_IMAGE031
Wherein
Figure 356929DEST_PATH_IMAGE032
，
Figure 288982DEST_PATH_IMAGE033
Which is indicative of the sampling frequency, is,
Figure 725779DEST_PATH_IMAGE034
intercepting the signal in a time range of
Figure 290753DEST_PATH_IMAGE035
，
Figure 736778DEST_PATH_IMAGE036
Meaning that the rounding is done down,
Figure 320074DEST_PATH_IMAGE037
，
Figure 927773DEST_PATH_IMAGE038
to represent
Figure 183305DEST_PATH_IMAGE039
A matrix of complex numbers is maintained.
3. The method according to claim 1, wherein the specific process of step 102 is as follows:
determining a location area based on beam pointing directions and spatial coordinate positions of receiving antennas
Figure 154060DEST_PATH_IMAGE040
：
Distance to distance range
Figure 873754DEST_PATH_IMAGE041
Azimuth distance range
Figure 652354DEST_PATH_IMAGE042
Wherein, in the step (A),
Figure 191920DEST_PATH_IMAGE043
is expressed by
Figure 963436DEST_PATH_IMAGE044
As a starting point, the method comprises the following steps of,
Figure 475320DEST_PATH_IMAGE045
is an end point, interval
Figure 424821DEST_PATH_IMAGE046
The sampling is carried out and the sampling is carried out,
Figure 700951DEST_PATH_IMAGE047
to represent
Figure 26890DEST_PATH_IMAGE048
The matrix of real numbers is then maintained,
Figure 393280DEST_PATH_IMAGE049
the total distance is counted to the distance unit;
beam pointing and location area based on receiving antenna
Figure 762951DEST_PATH_IMAGE040
Determining squint angle range for multi-view sub-images
Figure 277108DEST_PATH_IMAGE050
Wherein, in the step (A),
Figure 406739DEST_PATH_IMAGE051
total number of squint angles;
for positioning area
Figure 424373DEST_PATH_IMAGE040
Sum matched filter duration
Figure 430857DEST_PATH_IMAGE052
Obtaining an oblique view angle
Figure 166731DEST_PATH_IMAGE053
The effective time range of the signals to be processed is as follows:
Figure 834473DEST_PATH_IMAGE054
wherein
Figure 221461DEST_PATH_IMAGE055
Is the speed of the radar, and,
Figure 949246DEST_PATH_IMAGE056
，
Figure 172417DEST_PATH_IMAGE057
；
order to
Figure 830800DEST_PATH_IMAGE058
，
Figure 823027DEST_PATH_IMAGE059
If, if
Figure 987292DEST_PATH_IMAGE060
And is
Figure 432180DEST_PATH_IMAGE061
Then, then
Figure 690991DEST_PATH_IMAGE062
，
Figure 6566DEST_PATH_IMAGE063
，
Figure 810574DEST_PATH_IMAGE064
，
Figure 526114DEST_PATH_IMAGE065
(ii) a If it is not
Figure 73770DEST_PATH_IMAGE066
Or
Figure 40589DEST_PATH_IMAGE067
The truncated Doppler discrete signal obtained in step 101
Figure 15498DEST_PATH_IMAGE068
Oblique angle cannot be obtained
Figure 153087DEST_PATH_IMAGE069
Lower positioning area
Figure 504434DEST_PATH_IMAGE070
The positioning image of (1) needs to skip the oblique view
Figure 591339DEST_PATH_IMAGE071
The next oblique view angle is performed in the next positioning process
Figure 251996DEST_PATH_IMAGE072
The positioning process of (2);
for positioning area
Figure 158772DEST_PATH_IMAGE070
Considering the influence of square demodulation, based on the demodulated Doppler center frequency, instantaneous modulation frequency and squint angle
Figure 251493DEST_PATH_IMAGE073
In relation to (2)
Figure 442172DEST_PATH_IMAGE074
，
Figure 24463DEST_PATH_IMAGE075
Obtaining an oblique view
Figure 418535DEST_PATH_IMAGE076
Center frequency of lower Doppler
Figure 111685DEST_PATH_IMAGE077
And instantaneous frequency modulation rate
Figure 357203DEST_PATH_IMAGE078
；
For the signal
Figure 110395DEST_PATH_IMAGE079
Are mixed to obtain
Figure 991763DEST_PATH_IMAGE080
；
For the signal
Figure 737871DEST_PATH_IMAGE081
Low-pass filtering to obtain
Figure 388296DEST_PATH_IMAGE082
Denotes a cut-off frequency of
Figure 312389DEST_PATH_IMAGE083
The low-pass filter of (1);
for the signal
Figure 133584DEST_PATH_IMAGE084
Down-sampling to obtain
Figure 168536DEST_PATH_IMAGE085
，
Figure 673466DEST_PATH_IMAGE086
，
Figure 34041DEST_PATH_IMAGE087
Is the multiple of the down-sampling,
Figure 139269DEST_PATH_IMAGE088
，
Figure 915595DEST_PATH_IMAGE089
，
Figure 540611DEST_PATH_IMAGE090
，
Figure 324284DEST_PATH_IMAGE091
represents rounding up;
according to the sampling frequency
Figure 401961DEST_PATH_IMAGE092
And down-sampling multiple to determine location area
Figure 44295DEST_PATH_IMAGE093
The division of the azimuth distance is:
Figure 523818DEST_PATH_IMAGE094
wherein
Figure 413146DEST_PATH_IMAGE095
The total number of azimuth distance units;
design squint angle
Figure 978119DEST_PATH_IMAGE096
The following matched filter:
Figure 158565DEST_PATH_IMAGE097
wherein
Figure 7441DEST_PATH_IMAGE098
；
For distance in the radial direction
Figure 615140DEST_PATH_IMAGE099
Using oblique angles of view
Figure 870672DEST_PATH_IMAGE100
Lower corresponding matched filter
Figure 838497DEST_PATH_IMAGE101
For the Doppler signals obtained by processing
Figure 292612DEST_PATH_IMAGE102
Performing matched filtering to obtain:
Figure 336791DEST_PATH_IMAGE103
wherein
Figure 876357DEST_PATH_IMAGE104
For the convolution calling function in Matlab software,
Figure 871643DEST_PATH_IMAGE105
；
will be provided with
Figure 180264DEST_PATH_IMAGE106
As a result of the positioning of the azimuth distance range:
Figure 129766DEST_PATH_IMAGE107
wherein
Figure 405895DEST_PATH_IMAGE108
，
Figure 731834DEST_PATH_IMAGE109
The expression is taken to be the minimum value,
Figure 894962DEST_PATH_IMAGE110
representing an oblique angle of view
Figure 467895DEST_PATH_IMAGE111
Distance in the direction of distance
Figure 982053DEST_PATH_IMAGE112
The lower azimuth positioning result;
according to different distance directions
Figure 846104DEST_PATH_IMAGE112
Downward azimuth distance positioning result
Figure 129318DEST_PATH_IMAGE113
Obtain an oblique view angle
Figure 935469DEST_PATH_IMAGE111
Lower positioning area
Figure 671343DEST_PATH_IMAGE093
Positioning image of
Figure 276768DEST_PATH_IMAGE114
。
4. The method according to claim 1, wherein the specific process of step 103 is as follows:
the different squint angles obtained in step 102
Figure 666686DEST_PATH_IMAGE116
The radiation source positioning image of the positioning area is obtained by non-coherent accumulation of the lower radiation source positioning image
Figure 394470DEST_PATH_IMAGE118
：
Figure 617641DEST_PATH_IMAGE120
；
Wherein
Figure 89074DEST_PATH_IMAGE122
；
According to the location area
Figure 330568DEST_PATH_IMAGE124
Range of (1)
Figure 166937DEST_PATH_IMAGE126
，
Figure 877404DEST_PATH_IMAGE128
And positioning the image
Figure 136216DEST_PATH_IMAGE118
Drawing a positioning area
Figure 248529DEST_PATH_IMAGE124
And judging whether a radiation source target exists or not according to the positioning image, and if the target exists, finding the target and determining the position of the radiation source target.
CN202110213207.7A 2021-02-26 2021-02-26 Passive synthetic aperture radiation source positioning method based on multi-view processing Active CN112596056B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110213207.7A CN112596056B (en) 2021-02-26 2021-02-26 Passive synthetic aperture radiation source positioning method based on multi-view processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110213207.7A CN112596056B (en) 2021-02-26 2021-02-26 Passive synthetic aperture radiation source positioning method based on multi-view processing

Publications (2)

Publication Number Publication Date
CN112596056A true CN112596056A (en) 2021-04-02
CN112596056B CN112596056B (en) 2021-10-08

Family

ID=75207639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110213207.7A Active CN112596056B (en) 2021-02-26 2021-02-26 Passive synthetic aperture radiation source positioning method based on multi-view processing

Country Status (1)

Country Link
CN (1) CN112596056B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113625273A (en) * 2021-06-29 2021-11-09 西安电子科技大学 Aliasing digital signal synthetic aperture positioning method
CN113848549A (en) * 2021-09-15 2021-12-28 电子科技大学 Radiation source positioning method based on synthetic aperture technology
CN114624649A (en) * 2022-03-11 2022-06-14 北京理工大学 Method for positioning airborne passive synthetic aperture radiation source insensitive to residual frequency offset
CN115856888A (en) * 2022-12-07 2023-03-28 北京理工大学 Radiation source positioning method based on back projection
CN115877318A (en) * 2022-12-07 2023-03-31 北京理工大学 Radiation source positioning method based on multi-aperture cross positioning

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111352067A (en) * 2020-04-01 2020-06-30 北京理工大学 Passive positioning method for multi-target passive synthetic aperture
CN111736187A (en) * 2020-06-22 2020-10-02 北京理工大学 High-precision high-sensitivity single-satellite GNSS positioning method based on passive synthetic aperture

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111352067A (en) * 2020-04-01 2020-06-30 北京理工大学 Passive positioning method for multi-target passive synthetic aperture
CN111736187A (en) * 2020-06-22 2020-10-02 北京理工大学 High-precision high-sensitivity single-satellite GNSS positioning method based on passive synthetic aperture

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ANG LI ET AL.: ""Passive Synthetic Aperture High-Precision Radiation Source Location by Single Satellite"", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
HOSSEIN HEIDARI SHARIF ABADI ET AL.: ""A Two-Dimensional Spectrum Estimation of Target Point in Passive SAR Based on Pulsed Sources"", 《THE 22ND IRANIAN CONFERENCE ON ELECTRICAL ENGINEERING (ICEE 2014)》 *
张莉婷 等: ""基于被动合成孔径的单星无源高精度定位方法"", 《航天电子对抗》 *
朱建丰 等: ""基于合成孔径阵列的雷达辐射源被动定位技术研究"", 《电子学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113625273A (en) * 2021-06-29 2021-11-09 西安电子科技大学 Aliasing digital signal synthetic aperture positioning method
CN113625273B (en) * 2021-06-29 2023-12-22 西安电子科技大学 Synthetic aperture positioning method for aliasing digital signals
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
CN114624649A (en) * 2022-03-11 2022-06-14 北京理工大学 Method for positioning airborne passive synthetic aperture radiation source insensitive to residual frequency offset
CN114624649B (en) * 2022-03-11 2024-03-19 北京理工大学 Method for positioning airborne passive synthetic aperture radiation source insensitive to residual frequency offset
CN115856888A (en) * 2022-12-07 2023-03-28 北京理工大学 Radiation source positioning method based on back projection
CN115877318A (en) * 2022-12-07 2023-03-31 北京理工大学 Radiation source positioning method based on multi-aperture cross positioning
CN115877318B (en) * 2022-12-07 2024-04-16 北京理工大学 Radiation source positioning method based on multi-aperture cross positioning
CN115856888B (en) * 2022-12-07 2024-04-19 北京理工大学 Radiation source positioning method based on back projection

Also Published As

Publication number Publication date
CN112596056B (en) 2021-10-08

Similar Documents

Publication Publication Date Title
CN112596056A (en) Passive synthetic aperture radiation source positioning method based on multi-view processing
Zhang et al. Extending reliability of mmwave radar tracking and detection via fusion with camera
CN113009442B (en) Method and device for identifying multipath target of radar static reflecting surface
CN213581338U (en) Radar-guided photoelectric target detection system
CN111487597B (en) Universal electronic reconnaissance interception simulation method based on time-space-frequency energy data
CN109143239A (en) A kind of imaging method of the circumferential synthetic aperture radar based on one-dimensional range profile
CN105403886A (en) Automatic extraction method for airborne SAR scaler image position
CN113625273A (en) Aliasing digital signal synthetic aperture positioning method
CN115267777A (en) Single-star strabismus positioning method based on passive synthetic aperture
CN115856888B (en) Radiation source positioning method based on back projection
JP6923799B2 (en) Radar device
US6816109B1 (en) Method for automatic association of moving target indications from entities traveling along known route
CN116125488A (en) Target tracking method, signal fusion method, device, terminal and storage medium
CN111965637B (en) Method for analyzing and positioning projectile body based on double/multi-base radar
CN116068503A (en) Combined calibration method and device for millimeter wave radar and laser radar and terminal equipment
US6356227B1 (en) Smearing compensation apparatus for a radar system
US9176218B1 (en) Gain projection geo-location systems and methods
CN116827458B (en) Multi-angle-based frequency offset estimation and interference source positioning method
CN112835025A (en) Target speed ambiguity resolution method and device based on radar detection and terminal equipment
US20210063572A1 (en) Spatial recognition device, spatial recognition method, and program
CN105487053B (en) A kind of land photomask manufacturing process based on surveillance radar
CN109917334A (en) The method that aircraft platforms intercept and capture sensibility is quantitatively evaluated
CN115326093B (en) Method and device for calculating track of vehicle coordinate system, storage medium and electronic equipment
CN113296062B (en) SAR deception template generation method, system and application based on retrieval optimal value
CN113625260B (en) Three-dimensional speed measurement system and method

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