CN104316945A - Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering - Google Patents
Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering Download PDFInfo
- Publication number
- CN104316945A CN104316945A CN201410657732.8A CN201410657732A CN104316945A CN 104316945 A CN104316945 A CN 104316945A CN 201410657732 A CN201410657732 A CN 201410657732A CN 104316945 A CN104316945 A CN 104316945A
- Authority
- CN
- China
- Prior art keywords
- satellite
- kalman filtering
- estimation
- unscented kalman
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/428—Determining position using multipath or indirect path propagation signals in position determination
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 invention relates to a satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering. The effects of accurate estimation of time delay difference in a time-varying and non-Gaussian-noise environment and the stability and the real-time performance of a time difference positioning algorithm are mainly achieved. The positioning method includes the steps that (1), preprocessing such as zero mean and normalization is carried out on received satellite interference signals, and the high-order cumulants of the received satellite interference signals are calculated; (2), an improved time delay estimation method based on four-order cumulants minimum mean square errors is used for estimating the time delay difference generated by retransmitting the interference signals through three satellites; (3), according to estimation of the time delay difference, initial estimation of the position of the interference signals is obtained by using an analytical method, and then the positioning result of the previous moment is modified and updated by using an unscented Kalman filtering method. Compared with the prior art and method, a satellite interference source can be positioned in the time-varying and non-Gaussian-noise environment, and the positioning method is high in convergence rate, high in stability, easy to implement and small in calculation quantity.
Description
One technical field
The present invention relates to a kind of Satellite tool kit method, particularly a kind of satellite interference source three star problem method based on Higher Order Cumulants and Unscented kalman filtering.
Two background technologies
Satellite communication has the features such as broad covered area, communication distance is long, communication line is stable, communication band is wide, capacity is large, has played indispensable effect in national defence and commercial communication field.But because satellite transit is on disclosed space orbit, even destroy so satellite system is easily interfered.Wherein, that have a mind to from ground or Radio frequency interference (RFI) is day by day serious to the threat of telstar unintentionally.No matter the elimination of which kind of interference, first all accurately must locate interference source, utilize interference source localization method, when not changing telstar duty, test, estimate the approximate location of ground interference source, coordinate other means more on this basis, effectively can determine the accurate location of ground interference source.So Satellite tool kit method is to ensureing the effective of satellite communication system and reliability service, play a very important role.
Double-Star Positioning System technology is the main method of the not clear interference source in location at present, and it is by realizing interference source location to the estimation of time of arrival poor (TDOA) and arrival rate difference (FDOA).But there is the Accurate Measurement difficulty of frequency shift parameters in Double-Star Positioning System technology, and affects seriously by factors such as attitude of satellite changes, generally needs to set up multiple reference station to carry out error correction, the positioning precision that just can be comparatively satisfied with.
Samsung time difference positioning method utilizes the satellite group of space three satellite compositions in a distance, the signal of measurement ground interference source propagates into the mistiming between two between satellite, calculate two hyperboloids formed thus, and utilize the intersection point of itself and interference source place earth sphere, draw interference source position.Different with Double-Star Positioning System algorithm, this is that one utilizes two with reference to satellites, the method for being located by twice calculation delay, only make use of delay parameter, relative to frequency shift parameters, the mensuration of delay parameter is comparatively easily with accurate, finally higher to the positioning precision of interference source.
The performance of Samsung time difference positioning method depends on the estimated accuracy to difference time of arrival.Mainly contain about delay time estimation method: generalized correlation method, maximum likelihood estimate, least-mean-square error algorithm and high-order statistic method etc.Wherein, generalized correlation method is early stage delay time estimation method, and its principle is simple, calculated amount is also very little, but the precision estimated is not high, so more in the applications that background condition meets and accuracy requirement is not too high.Maximum likelihood estimate generally needs the probability density understanding signal in advance, and in practice, and the probability distribution of signal is difficult to obtain often, and maximum likelihood estimate has very large limitation when applying.Least mean-square error Time Delay Estimation Method is widely applicable, and calculated amount is little, and its real-time is better, practical, but they can not suppress the impact of correlation noise at all.The Higher Order Cumulants perseverance of Gaussian process is zero, therefore high-order statistic method is very effective when estimating the signal time delay difference with Gaussian Correlative Noise, but in practical application, the noise being mingled in signal may be non-gaussian Complex Noise, and signal to noise ratio (S/N ratio) also may change in time, during high-order statistic method pair, when change, non-Gaussian noise, satisfactory performance can not be obtained.Therefore, time delay estimation problem under change when existing delay parameter method of estimation needs to solve, non-Gaussian noise situation.
On the other hand, Localization Estimate Algorithm of TDOA can be summed up as the Solve problems of Nonlinear System of Equations.It solves positioning equation according to one group of time difference of measuring, and solves position of interference source exactly from geometrically saying by two crossing time difference hyperboloids and ground spheres intersect, its essence is to combine three ternary quadratic equations to solve.At present, that commonly uses in Localization Estimate Algorithm of TDOA has analytical method, Newton iteration method, least-squares iteration method, particle filter etc.Analytical method a kind ofly unknown number is directly calculated to the method solved, and generally has two solutions, can produce location ambiguity.Newton iterative relates to the calculating of trigonometric function, and calculated amount is larger.Least-squares iteration legal position stability is poor, and not easily restrains.Particle filter needs to adopt a large amount of particle samplers to go the error of approaching to reality to distribute, and calculated amount is large, and therefore real-time is poor.Therefore, need to improve existing Localization Estimate Algorithm of TDOA to solve the problem such as stability, real-time.
Three summary of the invention
The object of the invention is the problem such as stability, real-time in order to solve time-vary delay system tracking and Localization Estimate Algorithm of TDOA in Satellite tool kit delay parameter method of estimation, estimate to combine with Unscented kalman filtering by based on the fourth order cumulant least mean-square error time delay improved, propose a kind of satellite interference source three star problem method based on Higher Order Cumulants and Unscented kalman filtering.
Realizing step of the present invention is:
The first, the pre-service such as zero-mean, normalization is carried out to the satellite interference signal received in ground receiving station, and calculates its Higher Order Cumulants;
The second, utilize improve based on fourth order cumulant least mean-square error delay time estimation method, estimate undesired signal and forward the delay inequality produced through three satellites;
3rd, according to the estimation of delay inequality, utilize analytical method to obtain the initial estimation of target location, then utilized the positioning result of Unscented kalman filtering method to a upper moment to carry out correction and upgrade.
Beneficial effect of the present invention is: 1, have employed improvement based on fourth order cumulant least mean-square error delay time estimation method, the method effectively can suppress Gaussian noise and can accurately estimate time-vary delay system parameter; 2, utilize Unscented kalman filtering method constantly to revise renewal to positioning result, the method fast convergence rate, stability are high, and realize simple, and calculated amount is little.
Four accompanying drawing explanations
Fig. 1 is satellite interference source Samsung positioning using TDOA schematic diagram;
Fig. 2 be improve based on fourth order cumulant least mean-square error delay time estimation method theory diagram;
Fig. 3 is the positioning result correction renewal process flow diagram based on Unscented kalman filtering.
Five embodiments
Below in conjunction with accompanying drawing, the invention will be further described.
1, Fig. 1 gives satellite interference source Samsung positioning using TDOA system schematic.The undesired signal that ground interference source is launched, through the forwarding of different satellite, receives by processing enter is unified.First processing enter carries out the pre-service such as zero-mean, normalization to signal, then the mistiming parameter that three tunnel undesired signals arrive is estimated, recycling analytical method obtains the initial estimation of target location, finally utilizes the positioning result of Unscented kalman filtering method to a upper moment to carry out correction and upgrades.
2, Fig. 2 give improvement based on fourth order cumulant least mean-square error delay time estimation method theory diagram.
The model that undesired signal time delay is estimated can be expressed as follows:
Wherein k is discrete-time variable, and D is time delay, and A is decay factor, and x (k), y (k) are respectively to be had different delay and has superposed the undesired signal of noise, and s (k) is undesired signal, w
1(k) and w
2k () is respectively ground unrest, and suppose that it is uncorrelated, zero-mean independent with undesired signal mutually.
According to convolution theory, inhibit signal s (k-D) can be represented as
Wherein
sinc function is the even function being axis of symmetry with it to the maximum, and meets
Namely, when time variable n is tending towards infinite, sinc function decays to 0.In actual applications, usually n is limited in a suitable scope, namely can is inhibit signal approximate representation
Replace infinite with an abundant large integer p, and make p > | D|.
So inhibit signal model can be changed into
Wherein,
b
i=Asinc(i-D) (6)
For zero-mean stationary stochastic process, the as follows from fourth order cumulant of stochastic variable x (k) can be defined:
C
xxxx(τ,0,0)=cum[x(k),x(k+τ),x(k),x(k)] (7)
Definition stochastic variable x (k), the Cross fourth order cumulant of y (k) is as follows:
C
xyxx(τ,0,0)=cum[x(k),y(k+τ),x(k),x(k)] (8)
Formula (5) is substituted into (8), can obtain:
And the fourth order cumulant of Gaussian noise is zero, then
Weight coefficient vector is made to be:
W=[b
-p,b
-p+1,…,b
p]
T (11)
Cross fourth order cumulant vector is:
C
xy=[C
xyxx(-p,0,0),C
xyxx(-p+1,0,0),…,C
xyxx(p,0,0)]
T (12)
From fourth order cumulant vector be:
Then formula (10) can be expressed as with vector form:
C
xy=C
4x·W (14)
Wherein C
4xfor (2p+1) × (2p+1) dimension, W is that (2p+1) × 1 is tieed up, C
xyfor (2p+1) × 1 is tieed up.
Input signal x (k) is through calculating C from fourth order cumulant
xxxx(τ, 0,0), then obtain outputing signal C by auto-adaptive fir filter
xzxx(τ, 0,0), and the C calculated with the Cross fourth order cumulant of x (k) and y (k)
xyxx(τ, 0,0) subtracts each other, and obtains error e (k).
Definition adaptive error function
for
Wherein, weight coefficient vector is
The gradient of error function is:
Adopt steepest descent method to realize the iterative process of optimized parameter vector in adaptive process.Make initial weight vector
each numerical value changing current vector with the increment being proportional to negative gradient vector, namely
Wherein,
μ is converging factor, is used for regulating adaptive stability and speed.
Formula (17) represents with scalar, then
Wherein
In adaptive process, make FIR filter coefficient close to the time delays parameter of channel, namely
make error function
minimum, in last hunting time delay parameter estimated value
maximum value, thus indirect estimation time delay D.
3, Fig. 3 gives and carries out correction renewal process flow diagram based on Unscented kalman filtering method to positioning result
Unscented kalman filtering is a kind of filtering based on deterministic sampling.It avoids by the method for deterministic sampling carries out linearization by nonlinear equation, and can obtain and be accurate to the average suitable with three rank taylor series expansions and variance, performance exceedes EKF, and calculated amount increases few.
Its step is as follows:
1. initialization
In formula,
be respectively interference position with X (0) and estimate initial value and actual position initial value, P
0for error initial value.
2. Sigma value and weighting coefficient is calculated
In formula, χ is Sigma value, W weighting coefficient,
for the average of location variable X (k), L is the quantity of Sigma point; λ=γ
2(L+ κ)-L is a compositely proportional factor; Constant γ determines that Sigma point is at mean value
spreading range, it often get one very little on the occasion of, as 0≤γ≤10
-4; Constant κ is also a scale factor, and normal value is 3-L; β with
distribution relevant, when for Gaussian distribution, β=2;
for the i-th row of On Square-Rooting Matrices.
3. the time upgrades
In formula, Φ is state-transition matrix,
for error matrix, Q is noise covariance matrix.
4. measurement updaue
In formula, h (.) is for measuring transforming function transformation function.
In formula, R is state-noise covariance matrix.
Claims (3)
1., based on a satellite interference source three star problem method for Higher Order Cumulants and Unscented kalman filtering, the steps include: the first, the pre-service such as zero-mean, normalization is carried out to the satellite interference signal received, and calculate its Higher Order Cumulants; The second, utilize improve based on fourth order cumulant least mean-square error delay time estimation method, estimate undesired signal and forward the delay inequality produced through three satellites; 3rd, according to the estimation of delay inequality, utilize analytical method to obtain the initial estimation of undesired signal position, then utilized the positioning result of Unscented kalman filtering method to a upper moment to carry out correction and upgrade.
2. method according to claim 1, it is characterized in that: propose improvement based on fourth order cumulant least mean-square error delay time estimation method, utilize two-way undesired signal from fourth order cumulant and Cross fourth order cumulant to build error function, utilize steepest descent method to find the minimum value of error function, make filter coefficient close to the delay parameter of channel, finally find the maximum value in delay parameter estimated value.
3. method according to claim 1, it is characterized in that: propose the positioning result correction update method based on Unscented kalman filtering device, the undesired signal position utilizing analytical method to calculate is as the initial value of Unscented kalman filtering, the time adopting the method for deterministic sampling to complete filtering upgrades and measurement updaue, thus constantly carries out correction renewal to the positioning result in a upper moment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410657732.8A CN104316945A (en) | 2014-11-13 | 2014-11-13 | Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410657732.8A CN104316945A (en) | 2014-11-13 | 2014-11-13 | Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104316945A true CN104316945A (en) | 2015-01-28 |
Family
ID=52372198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410657732.8A Pending CN104316945A (en) | 2014-11-13 | 2014-11-13 | Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316945A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105609112A (en) * | 2016-01-15 | 2016-05-25 | 苏州宾果智能科技有限公司 | Sound source positioning method and apparatus and time delay estimation method and apparatus |
CN105939184A (en) * | 2016-03-04 | 2016-09-14 | 哈尔滨工业大学深圳研究生院 | UKF based aerospace DTN network bundle transmission delay estimation algorithm |
CN108919628A (en) * | 2018-05-15 | 2018-11-30 | 原时(荆门)电子科技有限公司 | A kind of Kalman filtering and Fir filtering control method and system to combining for atomic clock |
CN109946725A (en) * | 2017-12-20 | 2019-06-28 | 慧众行知科技(北京)有限公司 | A kind of satellite positioning method and system |
CN110879406A (en) * | 2019-11-29 | 2020-03-13 | 慧众行知科技(北京)有限公司 | Reference time difference determining method and system |
CN112731480A (en) * | 2020-11-10 | 2021-04-30 | 北京航空航天大学 | Ground signal source direct positioning method based on low-earth orbit satellite |
CN114665991A (en) * | 2022-05-23 | 2022-06-24 | 中国海洋大学 | Short wave time delay estimation method, system, computer equipment and readable storage medium |
-
2014
- 2014-11-13 CN CN201410657732.8A patent/CN104316945A/en active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105609112A (en) * | 2016-01-15 | 2016-05-25 | 苏州宾果智能科技有限公司 | Sound source positioning method and apparatus and time delay estimation method and apparatus |
CN105939184A (en) * | 2016-03-04 | 2016-09-14 | 哈尔滨工业大学深圳研究生院 | UKF based aerospace DTN network bundle transmission delay estimation algorithm |
CN105939184B (en) * | 2016-03-04 | 2019-10-22 | 哈尔滨工业大学深圳研究生院 | The method of empty day DTN network bundle propagation delay time estimation based on Unscented kalman filtering |
CN109946725A (en) * | 2017-12-20 | 2019-06-28 | 慧众行知科技(北京)有限公司 | A kind of satellite positioning method and system |
CN108919628A (en) * | 2018-05-15 | 2018-11-30 | 原时(荆门)电子科技有限公司 | A kind of Kalman filtering and Fir filtering control method and system to combining for atomic clock |
CN110879406A (en) * | 2019-11-29 | 2020-03-13 | 慧众行知科技(北京)有限公司 | Reference time difference determining method and system |
CN112731480A (en) * | 2020-11-10 | 2021-04-30 | 北京航空航天大学 | Ground signal source direct positioning method based on low-earth orbit satellite |
CN112731480B (en) * | 2020-11-10 | 2023-09-29 | 北京航空航天大学 | Ground signal source direct positioning method based on low-orbit satellite |
CN114665991A (en) * | 2022-05-23 | 2022-06-24 | 中国海洋大学 | Short wave time delay estimation method, system, computer equipment and readable storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104316945A (en) | Satellite interference source three-satellite positioning method based on high-order cumulants and unscented Kalman filtering | |
CN109298388B (en) | Direct over-the-horizon target geographic coordinate estimation method based on azimuth information | |
CN107677272B (en) | AUV (autonomous Underwater vehicle) collaborative navigation method based on nonlinear information filtering | |
CN105954712B (en) | The direct localization method of the multiple target of associated wireless electric signal complex envelope and carrier phase information | |
CN104714244A (en) | Multi-system dynamic PPP resolving method based on robust self-adaption Kalman smoothing | |
CN105353345B (en) | A kind of single channel radar signal passive location method based on high-order observed quantity | |
CN106599368A (en) | FastSLAM method based on particle proposal distribution improvement and adaptive particle resampling | |
CN105301601A (en) | Global navigation satellite system (GNSS) ionosphere delayed three-dimensional modeling method suitable for global area | |
CN103278813A (en) | State estimation method based on high-order unscented Kalman filtering | |
CN103955600A (en) | Target tracking method and truncated integral Kalman filtering method and device | |
CN111551897B (en) | TDOA (time difference of arrival) positioning method based on weighted multidimensional scaling and polynomial root finding under sensor position error | |
CN110657806A (en) | Position resolving method based on CKF, chan resolving and Savitzky-Golay smooth filtering | |
Chhetri et al. | Scheduling multiple sensors using particle filters in target tracking | |
CN103675880B (en) | Lasting air navigation aid under a kind of satellite-signal congestion situations | |
Rahman et al. | Earth-centered Earth-fixed (ECEF) vehicle state estimation performance | |
CN103296995A (en) | Unscented transformation and unscented Kalman filtering method in any-dimension high order (>/=4) | |
CN109212466A (en) | A kind of broadband direction-finding method based on quantum dragonfly mechanism of Evolution | |
Zhou et al. | A window-recursive approach for GNSS kinematic navigation using pseudorange and Doppler measurements | |
CN113219406B (en) | Direct tracking method based on extended Kalman filtering | |
CN110456335A (en) | A kind of MIMO radar object localization method neural network based and system | |
KR101887877B1 (en) | Multistatic passive coherent location system and method for estimating target location using fm network | |
CA2682438A1 (en) | Cognitive tracking radar | |
Chen et al. | Study on GPS/INS loose and tight coupling | |
Song et al. | Square-root cubature FastSLAM algorithm for mobile robot simultaneous localization and mapping | |
Dey et al. | Joint estimation of parameters and states of nonlinear systems using adaptive divided difference filter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 210007 Qinhuai District, Jiangsu Province, after the standard Camp No. 18, No. Applicant after: The 36th Institute of Central Military Commission Equipment Development Department Address before: 210007 Qinhuai District, Jiangsu Province, after the standard Camp No. 18, No. Applicant before: No. 63 Inst. of the Headquarters of the Genearal Staff of C.P. L. A. |
|
COR | Change of bibliographic data | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150128 |