CN115201803B - Passive positioning method and system based on satellite overhead moment measurement - Google Patents
Passive positioning method and system based on satellite overhead moment measurement Download PDFInfo
- Publication number
- CN115201803B CN115201803B CN202210825979.0A CN202210825979A CN115201803B CN 115201803 B CN115201803 B CN 115201803B CN 202210825979 A CN202210825979 A CN 202210825979A CN 115201803 B CN115201803 B CN 115201803B
- Authority
- CN
- China
- Prior art keywords
- distance
- azimuth
- radar
- frequency modulation
- satellite
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005259 measurement Methods 0.000 title claims abstract description 22
- 238000013508 migration Methods 0.000 claims abstract description 27
- 230000005012 migration Effects 0.000 claims abstract description 27
- 230000005855 radiation Effects 0.000 claims abstract description 16
- 230000006835 compression Effects 0.000 claims abstract description 14
- 238000007906 compression Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims description 14
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a passive positioning method and a passive positioning system based on satellite overhead time measurement, which realize the positioning method of a radiation source by measuring the satellite overhead time and solve the problem of low positioning precision of the traditional passive positioning method. The range migration curve is obtained through pulse compression of the range direction of the radar signal, and the over-top moment of the satellite is obtained through secondary fitting of the range migration curve.
Description
Technical Field
The invention belongs to the field of radiation source positioning, and particularly relates to a passive positioning method and system based on satellite overhead moment measurement.
Background
In radar countermeasure, radar signal reconnaissance aims at detecting, locating and identifying a radar signal source of an opposite party. The invention provides a method for realizing radiation source by measuring satellite overhead time in radar countermeasure. The current measurement information for passive positioning comprises direction finding positioning, time difference positioning and frequency finding positioning, wherein the direction finding positioning needs to obtain more accurate angle finding precision by increasing the length of an antenna base line, and has high requirements on a load antenna; the time difference positioning requires at least three observation stations, and has high requirements on common-view conditions and time synchronization; the frequency measurement positioning is carried out, the error is influenced by the frequency measurement precision, and the error is in kilometer level. The accurate measurement of the satellite over-radiation source top moment is realized, and is a key link in the radiation source positioning.
Disclosure of Invention
The invention aims to provide a passive positioning method and a passive positioning system based on satellite overhead moment measurement, in particular to a pulse Doppler radar radiation source, belongs to the field of radiation source positioning, and solves the problem of low positioning precision of a traditional passive positioning method. The range migration curve is obtained through pulse compression of the range direction of the radar signal, and the over-top moment of the satellite is obtained through secondary fitting of the range migration curve.
In order to solve the technical problems, the technical scheme of the invention is as follows:
A passive positioning method based on satellite overhead time measurement, the method comprising:
performing distance-direction matched filtering on the received radar baseband signal to obtain radar data after pulse compression;
obtaining a range migration curve of the radar by extracting peak information of the radar data after pulse compression;
Performing quadratic term fitting on the range migration curve to obtain the over-top moment and the azimuth distance;
Obtaining azimuth signals by extracting phases of peak information of the radar data after pulse compression;
estimating the frequency modulation rate of the azimuth signal, determining the search step of the frequency modulation according to the positioning precision, generating a plurality of groups of matched filters, carrying out matched filtering on the azimuth signal to obtain a matched filtering result under different frequency modulation, wherein the frequency modulation rate corresponding to the peak value is the frequency modulation rate of the azimuth signal, and obtaining the distance-to-distance of the radar through the relation between the frequency modulation rate and the distance-to-distance;
And positioning the radar radiation source through the azimuth distance and the distance.
Further, the radar baseband signal S 0 (τ, η) is specifically:
Where τ is distance to time, η is a complex constant at a near azimuth moment, coefficient a 0, ω r (τ) is a rectangular pulse signal, T r is pulse duration, where/> is received signal strength, θ is an angle to line of sight measured in an oblique plane, β bw is azimuth beam width, R (η) is a distance of a target from an observation platform, c is a speed of light, f 0 is a carrier frequency of the target radar, K r is an intra-pulse frequency of the target radar, η c is a beam center deviation moment.
Further, the radar data S rc (τ, η) after pulse compression is specifically:
Src(τ,η)＝IFFT{S0(f0,η)H(f0)}＝A0sinc(τ-R(η)/c)ωa(η-ηc)×exp(-j2πf0R(η)/c);
Wherein S 0(f0, η) is the fourier transform of S 0 (τ, η);
Further, performing quadratic term fitting on the range migration curve to obtain the over-top moment and the azimuth distance;
Wherein the fitting formula is:
f(x)＝P0x2+P1x+P2；
And f (x) is the range migration value obtained in the second step, x is the azimuth moment corresponding to the range migration value in the second step, and finally the fitted range migration curve is obtained.
Further, the over-top moment is
Further, the azimuth distance R a＝tp. V.
Further, the expression of the matched filter is:
Wherein K is the frequency modulation of the matched filter, T s is the synthetic aperture time, the generated matched filter and the obtained azimuth signal are subjected to matched filtering to obtain a matched filtering result under different frequency modulation, the frequency modulation corresponding to the peak value is the frequency modulation of the azimuth signal, and the distance-to-distance of the radar is obtained by the relation between the frequency modulation and the distance-to-distance R 0, and the relation is that:
where v is the flight speed of the satellite platform and λ is the wavelength of the radiation source.
A passive positioning system based on satellite overhead time measurement, the system comprising: one or more processors;
A memory for storing one or more programs;
The one or more programs, when executed by the one or more processors, cause the one or more processors to perform a passive positioning method based on satellite overhead time measurements as described above.
Compared with the prior art, the invention has the advantages that:
The invention relates to a passive positioning method and a passive positioning system based on satellite overhead moment measurement, which are used for positioning a radar radiation source by using a passive synthetic aperture positioning technology.
According to the method, the accurate estimation of the over-top moment is realized through the secondary fitting of the range migration curve, and the positioning accuracy is improved.
Drawings
FIG. 1, schematic diagram of radar signal localization;
FIG. 2, radar signal after range pulse compression;
FIG. 3, range migration curve;
FIG. 4, fitted range migration curve;
Fig. 5, azimuth signal.
Detailed Description
The following describes specific embodiments of the present invention with reference to examples:
It should be noted that the structures, proportions, sizes and the like illustrated in the present specification are used for being understood and read by those skilled in the art in combination with the disclosure of the present invention, and are not intended to limit the applicable limitations of the present invention, and any structural modifications, proportional changes or size adjustments should still fall within the scope of the disclosure of the present invention without affecting the efficacy and achievement of the present invention.
Also, the terms such as "upper," "lower," "left," "right," "middle," and "a" and the like recited in the present specification are merely for descriptive purposes and are not intended to limit the scope of the invention, but are intended to provide relative positional changes or modifications without materially altering the technical context in which the invention may be practiced.
Example 1:
as shown in fig. 1, the satellite flies over the radar target as shown in fig. 1, and the received radar baseband signal is
S0(τ,η)＝A0ωr(τ-R(η)/c)ωa(η-ηc)×exp(-j2πf0R(η)/c)exp(jπKr(τ-R(η)/c)2)
Where τ is distance to time, η is a complex constant at a near azimuth moment, coefficient a 0, ω r (τ) is a rectangular pulse signal, T r is pulse duration, where/> is received signal strength, θ is an angle to line of sight measured in an oblique plane, β bw is azimuth beam width, R (η) is a distance of a target from an observation platform, c is a speed of light, f 0 is a carrier frequency of the target radar, K r is an intra-pulse frequency of the target radar, η c is a beam center deviation moment.
1) The received radar baseband signal is subjected to distance-wise matched filtering to obtain radar data S rc (τ, η) after pulse compression, as shown in fig. 2:
Src(τ,η)＝IFFT{S0(f0,η)H(f0)}＝A0sinc(τ-R(η)/c)ωa(η-ηc)×exp(-j2πf0R(η)/c);
Wherein S 0(f0, eta) is the Fourier transform of S 0 (tau, eta)
2) And extracting peak information of the radar data after pulse compression to obtain a range migration curve of the radar, wherein the range migration curve is shown in figure 3.
3) And performing quadratic term fitting on the range migration curve according to the following formula to obtain the over-top moment of the range migration curve.
f(x)＝P0x2+P1x+P2；
Wherein f (x) is represented by the range migration value obtained in 2), x represents the corresponding azimuth moment in 2), and finally a fitted range migration curve is obtained, as shown in fig. 4.
The over-top moment is
Azimuth distance R a＝tp ·v.
4) The azimuth signal is obtained by extracting the phase at the peak information of the pulse-compressed radar data, as shown in fig. 5.
5) Estimating the frequency modulation rate of the azimuth signal, determining the search step of the frequency modulation rate according to the positioning accuracy, thereby generating a plurality of groups of filters, wherein the expression of the matched filters is as follows
Wherein K is the frequency modulation of the matched filter, T s is the synthetic aperture time, the generated matched filter and the obtained azimuth signal are subjected to matched filtering to obtain a matched filtering result under different frequency modulation, the frequency modulation corresponding to the peak value is the frequency modulation of the azimuth signal, and the distance-to-distance of the radar is obtained by the relation between the frequency modulation and the distance-to-distance R 0, wherein the relation is that
Where v is the flight speed of the satellite platform and λ is the wavelength of the radiation source.
6) And positioning the radar radiation source is realized through the azimuth distance and the distance.
Example 2:
the invention is described in detail below in connection with specific signal examples:
In the simulation experiment, the carrier frequency f 0 of the radar radiation source transmitting signal is 2GHz, the radar is pulse Doppler radar, the orbit height h of the satellite-borne receiver is 500km, the distance R 0 is 707km, the flying speed v of the satellite platform is 7000m/s, the synthetic aperture time T s is 0.69s, the pulse duration T r is 40us, the intra-pulse modulation frequency K r of the target radar is 5e11Hz/s, the pulse repetition frequency PRF is 1KHz, the baseband signal of the radar is sampled at the sampling rate of 100MHz, the distance pulse compression is carried out on the radar signal of the baseband according to the step 1), the amplitude of the output radar data is taken, as shown in fig. 2, and the range migration curve of the radar is obtained according to the step 2) as shown in fig. 3. Obtaining a fitted range migration curve according to the step 3), as shown in fig. 4, and obtaining the overhead time t p as 1.6432s and the azimuth distance R a as 11503m. The result of the matched filtering at different tuning frequencies according to step 4) and step 5) is shown in fig. 5, and tuning frequency K is 461.8Hz/s. Finally, the distance to the radar radiation source is 707377m. Finally, obtaining the position of the radar radiation source as a distance direction distance 707377m and an error 270m; the azimuthal distance was 11503m and the error was 3m.
While the preferred embodiments of the present invention have been described in detail, the present invention is not limited to the above embodiments, and various changes may be made without departing from the spirit of the present invention within the knowledge of those skilled in the art.
Many other changes and modifications may be made without departing from the spirit and scope of the invention. It is to be understood that the invention is not to be limited to the specific embodiments, but only by the scope of the appended claims.
Claims (8)
1. A passive positioning method based on satellite overhead time measurement, the method comprising:
performing distance-direction matched filtering on the received radar baseband signal to obtain radar data after pulse compression;
obtaining a range migration curve of the radar by extracting peak information of the radar data after pulse compression;
Performing quadratic term fitting on the range migration curve to obtain the over-top moment and the azimuth distance;
Obtaining azimuth signals by extracting phases of peak information of the radar data after pulse compression;
estimating the frequency modulation rate of the azimuth signal, determining the search step of the frequency modulation according to the positioning precision, generating a plurality of groups of matched filters, carrying out matched filtering on the azimuth signal to obtain a matched filtering result under different frequency modulation, wherein the frequency modulation rate corresponding to the peak value is the frequency modulation rate of the azimuth signal, and obtaining the distance-to-distance of the radar through the relation between the frequency modulation rate and the distance-to-distance;
And positioning the radar radiation source through the azimuth distance and the distance.
2. A passive positioning method based on satellite overhead time measurement according to claim 1, wherein the radar baseband signal S 0 (τ, η) is specifically:
Where τ is distance to time, η is near azimuth time, coefficient A 0 is a complex constant, ω r (τ) is a rectangular pulse signal, T r is pulse duration, where/> is received signal strength, θ is the angle to line of sight measured in the oblique plane, β bw is azimuth beam width, R (η) is the distance of the target from the observation platform, c is the speed of light, f 0 is the carrier frequency of the target radar, K r is the intra-pulse frequency of the target radar, η c is the beam center offset time.
3. The passive positioning method based on satellite overhead time measurement according to claim 1, wherein the radar data S rc (τ, η) after pulse compression is specifically:
Src(τ,,η)＝IFFT{S0(f0,η)H(f0)}＝A0sinc(τ-R(η)/c)ωa(η-ηc)×exp(-j2πf0R2(η)/c);
wherein S 0(f0, η) is the fourier transform of S 0 (τ, η);
4. The passive positioning method based on satellite overhead moment measurement according to claim 1, wherein a distance migration curve is subjected to quadratic term fitting to obtain overhead moment and azimuth distance; wherein the fitting formula is:
f(x)＝P0x2+P1x+P2；
And f (x) is the range migration value obtained in the second step, x is the azimuth moment corresponding to the range migration value in the second step, and finally the fitted range migration curve is obtained.
5. The passive positioning method based on satellite overhead time measurement according to claim 1, wherein the overhead time is
6. A passive positioning method based on satellite overhead time measurement according to claim 1, characterized in that the azimuth distance R a＝tp ·v.
7. The passive positioning method and system based on satellite overhead time measurement according to claim 1, wherein the expression of the matched filter is:
Wherein K is the frequency modulation of the matched filter, T s is the synthetic aperture time, the generated matched filter and the obtained azimuth signal are subjected to matched filtering to obtain a matched filtering result under different frequency modulation, the frequency modulation corresponding to the peak value is the frequency modulation of the azimuth signal, and the distance-to-distance of the radar is obtained by the relation between the frequency modulation and the distance-to-distance R 0, and the relation is that:
where v is the flight speed of the satellite platform and λ is the wavelength of the radiation source.
8. A passive positioning system based on satellite overhead time measurement, the system comprising: one or more processors;
A memory for storing one or more programs;
The one or more programs, when executed by the one or more processors, cause the one or more processors to perform a passive positioning method based on satellite overhead time measurements as recited in any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210825979.0A CN115201803B (en) | 2022-07-13 | 2022-07-13 | Passive positioning method and system based on satellite overhead moment measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210825979.0A CN115201803B (en) | 2022-07-13 | 2022-07-13 | Passive positioning method and system based on satellite overhead moment measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115201803A CN115201803A (en) | 2022-10-18 |
| CN115201803B true CN115201803B (en) | 2024-04-16 |
Family
ID=83579375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210825979.0A Active CN115201803B (en) | 2022-07-13 | 2022-07-13 | Passive positioning method and system based on satellite overhead moment measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115201803B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012233824A (en) * | 2011-05-06 | 2012-11-29 | Toshiba Corp | Passive radar device, guiding device and radio wave detecting method |
| JP2013029420A (en) * | 2011-07-28 | 2013-02-07 | Mitsubishi Electric Corp | Passive radar system |
| CN108693545A (en) * | 2018-03-29 | 2018-10-23 | 中国人民解放军国防科技大学 | Abnormal target positioning method based on satellite-borne ADS-B |
| CN110231616A (en) * | 2019-04-09 | 2019-09-13 | 电子科技大学 | Sea Moving objects detection and location method based on big-dipper satellite radiation source |
| CN111551922A (en) * | 2020-05-18 | 2020-08-18 | 西安电子科技大学 | Three-dimensional space double/multi-base radar high-speed target detection method |
| CN111736187A (en) * | 2020-06-22 | 2020-10-02 | 北京理工大学 | High-precision high-sensitivity single-satellite GNSS positioning method based on passive synthetic aperture |
| CN112560295A (en) * | 2021-02-20 | 2021-03-26 | 北京理工大学 | Satellite equivalent velocity calculation method for passive synthetic aperture positioning |
| CN112904326A (en) * | 2021-01-29 | 2021-06-04 | 哈尔滨工业大学 | Satellite-borne passive positioning method based on virtual aperture |
-
2022
- 2022-07-13 CN CN202210825979.0A patent/CN115201803B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012233824A (en) * | 2011-05-06 | 2012-11-29 | Toshiba Corp | Passive radar device, guiding device and radio wave detecting method |
| JP2013029420A (en) * | 2011-07-28 | 2013-02-07 | Mitsubishi Electric Corp | Passive radar system |
| CN108693545A (en) * | 2018-03-29 | 2018-10-23 | 中国人民解放军国防科技大学 | Abnormal target positioning method based on satellite-borne ADS-B |
| CN110231616A (en) * | 2019-04-09 | 2019-09-13 | 电子科技大学 | Sea Moving objects detection and location method based on big-dipper satellite radiation source |
| CN111551922A (en) * | 2020-05-18 | 2020-08-18 | 西安电子科技大学 | Three-dimensional space double/multi-base radar high-speed target detection method |
| CN111736187A (en) * | 2020-06-22 | 2020-10-02 | 北京理工大学 | High-precision high-sensitivity single-satellite GNSS positioning method based on passive synthetic aperture |
| CN112904326A (en) * | 2021-01-29 | 2021-06-04 | 哈尔滨工业大学 | Satellite-borne passive positioning method based on virtual aperture |
| CN112560295A (en) * | 2021-02-20 | 2021-03-26 | 北京理工大学 | Satellite equivalent velocity calculation method for passive synthetic aperture positioning |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115201803A (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7969345B2 (en) | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters | |
| CN105572635B (en) | The passive method for rapidly positioning in list station based on least square method | |
| CN102565784B (en) | Method of moving-target relocation and velocity ambiguity resolution based on velocity synthetic aperture radar (VSAR) system | |
| CN109856635B (en) | CSAR ground moving target refocusing imaging method | |
| EP1505408A1 (en) | A method for SAR processing without INS data | |
| Ruegg et al. | Capabilities of dual-frequency millimeter wave SAR with monopulse processing for ground moving target indication | |
| CN104749570B (en) | It is a kind of to move constant airborne biradical synthetic aperture radar target localization method | |
| CN105717480A (en) | Single-observer passive rapid location method based on phase difference | |
| CN114545411B (en) | Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization | |
| CN105717479A (en) | Single-station passive quick positioning method based on fractional order Fourier transformation | |
| CN102121990B (en) | Space-time analysis-based target rotation speed estimating method for inverse synthetic aperture radar | |
| CN112904326B (en) | Satellite-borne passive positioning method based on virtual aperture | |
| CN105954751A (en) | Ka FMCW SAR moving object detection method and apparatus | |
| CN107271997B (en) | Airborne multi-channel CSSAR ground moving target motion parameter estimation method | |
| CN105044667A (en) | Double-satellite tracking method, device and system for moving target | |
| CN105676174A (en) | Single-station passive rapid positioning method based on phase difference and least square method | |
| CN111751822B (en) | Time and phase synchronization method, system, equipment and medium | |
| CN115201803B (en) | Passive positioning method and system based on satellite overhead moment measurement | |
| Brooker et al. | W-band airborne interrupted frequency modulated CW imaging radar | |
| CN104459651B (en) | Method for estimating length of equivalent baseline of airborne SAR-GMTI system | |
| Nye et al. | Measuring the change in thickness of the Antarctic ice sheet | |
| CN105717481B (en) | Based on the passive method for rapidly positioning in list station for going tiltedly to handle | |
| Munoz-Ferreras et al. | Monitoring road traffic with a high resolution LFMCW radar | |
| CN102466790B (en) | Centrifugal acceleration measurement method for airborne electronic support measures (ESM) system | |
| Luomei et al. | Motion compensation for multirotors minisar system |
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 |