CN105510935A - Device used for navigation - Google Patents
Device used for navigation Download PDFInfo
- Publication number
- CN105510935A CN105510935A CN201510937376.XA CN201510937376A CN105510935A CN 105510935 A CN105510935 A CN 105510935A CN 201510937376 A CN201510937376 A CN 201510937376A CN 105510935 A CN105510935 A CN 105510935A
- Authority
- CN
- China
- Prior art keywords
- satellite
- module
- navigation
- calculate
- observed quantity
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
-
- 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
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 discloses a device used for navigation. The device comprises a capture and tracking processing module, a message extraction module, an observation quantity extraction module, a positioning calculation module and an output module. Satellite signals are continuously tracked via the capture and tracking processing module so that navigation telegraph message bit data of a satellite are obtained; and satellite navigation telegraph messages and observation quantity are obtained via the message extraction module and the observation quantity extraction module, then positioning calculation is performed via the positioning calculation module, and outputting is performed via the output module according to a standard format NMEA format formulated by marine electronic equipment so that the navigation telegraph message bit data of the satellite can be automatically and continuously obtained for positioning and the situation that navigation is not timely can be avoided, and thus the device used for navigation has the technical effects of being timely and accurate.
Description
Technical field
The application relates to technical field of satellite communication, particularly relates to a kind of device for navigating.
Background technology
Along with the development of science and technology, communication technology of satellite have also been obtained development at full speed, and the kind of electronic product also gets more and more, and people have also enjoyed the various facilities that development in science and technology brings.Present people can pass through various types of electronic equipment, enjoy the comfortable life along with development in science and technology brings.
Such as navigating instrument is the product of modern science and technology, navigating instrument is used on automobile, for location, navigation and amusement, along with the universal of automobile and the construction of road, the economic interaction of intercity is more frequent, automatic navigator seems very important, and the navigation that accurately location, navigation, amusement function roll into one more can meet the demand of car owner, becomes the basic equipment on car.
Such as, but the navigation of current navigating instrument out of true, often there will be navigation way mistake or navigation situation not in time, occur that navigation postpones, and can not obtain follow-up course in time, make car owner be at a loss in the process of navigation; Or by car owner's misguidance in the sea, increase car owner dangerous.
Summary of the invention
For solving the problems of the technologies described above, the invention provides a kind of device for navigating, comprising:
Acquisition and tracking processing module, continues to follow the tracks of for carrying out satellite-signal, and obtains the navigation message bit data of satellite;
Message extraction module, for receiving described navigation message bit data, and goes out satellite navigation message from described navigation message bit extracting data;
Observed quantity extraction module, for going out observed quantity from described navigation message bit extracting data;
Positioning calculation module, for carrying out co-located according to described observed quantity and described satellite navigation message, obtains navigation results;
Output module, for navigation results described in the standard format NMEA formatted output formulated with electronic equipment according to sea.
Preferably, described satellite navigation message comprises ephemeris, satellite health, time system, the satellite almanac of satellite.
Preferably, described observed quantity extraction module, specifically for obtaining the pseudo range observed quantity of tracked satellite or Doppler measurements or carrier phase observed quantity according to Base-Band Processing data.
Preferably, described pseudo range observed quantity is utilize the transmission time of satellite-signal to be multiplied by transmission range that the light velocity obtains signal; The transmission time of described satellite-signal refers to satellite-signal transmitting time t
svdeduct receiver local zone time t
u.
Preferably, described Doppler measurements is after the frequency values obtaining carrier wave NCO, by the frequency values of described carrier wave NCO and the nominal intermediate frequency value f of radio-frequency front-end
csubtract each other acquisition.
Preferably, described positioning calculation module specifically comprises:
Computing module, for going out the position of satellite according to the ephemeris computation of described satellite;
Positioning calculation module, to position customer location for utilizing the position of described satellite and described observed quantity and resolves.
Preferably, described computing module specifically for:
The first step: the mean angular velocity n calculating satellite motion;
Second step: the mean anomaly M of calculating observation moment satellite
s;
M
s=M
0+ n (t-t
oe); In formula, M
0for reference moment t
oetime mean anomaly, provided by broadcast ephemeris;
3rd step: calculate eccentric anomaly: E
s=M
s+ e
ssinE
s;
4th step: calculate anomaly f
s:
5th step: calculate to rise and hand over elongation
in formula, ω is the argument of perigee, is provided by broadcast ephemeris;
6th step: calculate perturbation correction member;
7th step: to a liter friendship elongation
satellite radius vector r, orbit inclination i carry out perturbation and correct;
8th step: calculate the position of satellite in track areal coordinate system;
9th step: the longitude L of calculating observation moment ascending node;
Tenth step: calculate the position of satellite in spherical coordinate system instantaneously;
11 step: calculate the position of satellite in conventional terrestrial coordinate system.
Preferably, described mean angular velocity n is obtained by following steps:
First according to the parameter provided in broadcast ephemeris
calculate with reference to moment t
oemean angular velocity n
0:
in formula, GM is the product of gravitational constant G and earth gross mass M, GM=3.9860047 × 1014m
3/ s
2, then according to the mean angular velocity n:n=n of perturbation parameter △ n calculating observation moment satellite given in broadcast ephemeris
0+ △ n.
Preferably, described calculating perturbation correction member comprises the following steps:
By the C provided in broadcast ephemeris
uc, C
us, C
rc, C
rs, C
ic, C
is6 perturbation parameters, obtain the perturbation correction member δ that the liter caused due to the second order zonal harmonic coefficient item of earth gravitational field bit function hands over elongation μ
μ, satellite radius vector r perturbation correction member δ
r, and the perturbation correction member δ of inclination of satellite orbit i
i, computing formula is as follows:
Preferably, described to a liter friendship elongation
the step that satellite radius vector r, orbit inclination i carry out perturbation correction is:
r=a
s(1-e
scosE
s)+δ
r
In formula: a
sfor the major radius of satellite orbit, i
0for t
oethe orbit inclination in moment, is provided by the kepler parameters in broadcast ephemeris;
for the rate of change of i, provided by the perturbation parameter in broadcast ephemeris.
By one or more technical scheme of the present invention, the present invention has following beneficial effect or advantage:
The invention provides a kind of device for navigating, comprising: acquisition and tracking processing module, message extraction module, observed quantity extraction module, positioning calculation module, output module.Carry out continuing to follow the tracks of the navigation message bit data obtaining satellite to satellite-signal by acquisition and tracking processing module; By message extraction module, observed quantity extraction module, obtain satellite navigation message and observed quantity, then positioned by positioning calculation module and resolve, and exported by the standard format NMEA form that output module is formulated according to sea electronic equipment, the navigation message bit data of acquisition satellite that can be automatic continuous position, and then avoid occurring that navigation situation not in time occurs, and have technique effect promptly and accurately.
Accompanying drawing explanation
Fig. 1 is the structural representation of device for navigating in the embodiment of the present invention.
Embodiment
In order to make the application the technical staff in the technical field more clearly understand the application, below in conjunction with accompanying drawing, by specific embodiment, technical scheme is described in detail.
A kind of device for navigating that the present invention describes, is mainly divided into four parts principle: acquisition and tracking processing section, message and observed quantity Extraction parts, positioning calculation and system export.
Based on above-mentioned principle, for the structure of device of navigating as shown in Figure 1.Device for navigating mainly comprises: acquisition and tracking processing module 11, message extraction module 12, observed quantity extraction module 13, positioning calculation module 14, output module 15.
Acquisition and tracking processing module 11, continues to follow the tracks of for carrying out satellite-signal, and obtains the navigation message bit data of satellite;
Message extraction module 12, for receiving described navigation message bit data, and goes out satellite navigation message from described navigation message bit extracting data;
Observed quantity extraction module 13, for going out observed quantity from described navigation message bit extracting data.
Positioning calculation module 14, for carrying out co-located according to described observed quantity and described satellite navigation message, obtains navigation results;
Output module 15, for navigation results described in the standard format NMEA formatted output formulated with electronic equipment according to sea.
Introduce the concrete function of modules below respectively.
In acquisition and tracking processing module 11, acquisition and tracking data processing comprises ALT-CH alternate channel state, adjustment is followed the tracks of carrier loop and code loop circuit state, maintains correct tracking mode, drawn the information such as correct Doppler, code phase and carrier phase, also can carry out the navigation message bit data of bit synchronous and satellite.
In the measurement moment, acquisition and tracking data processing module can encapsulate the observation value obtained, and passes to next unit (message extraction module 12, observed quantity extraction module 13) by system interface, carries out message and resolves and extract with observed quantity.
Message extraction module 12 for receiving described navigation message bit data, and goes out the navigation messages such as satellite ephemeris, almanac, UTC parameter, Ionospheric Parameters from described navigation message bit extracting data.And observed quantity extraction module 13, then mainly try to achieve the observed quantities such as the pseudorange of satellite, Doppler, carrier phase.
Lower mask body introduces message extraction module 12 and observed quantity extraction module 13.
Message extraction module 12.
Message extraction module 12, mainly decodes to bit stream information according to GPS and BD2 (a kind of positioning system) interface control document, obtains satellite navigation message.Satellite navigation message provides the data such as ephemeris, satellite health, time system, satellite almanac of concerned satellites for user, is that user positions and navigate necessary basic data.For the structure of different global position system navigation message, need to adopt different decoding processes, correctly to obtain navigation message efficiently.
On the L1 wave band of GPS, the transfer rate of people's code navigation message is 50bit/s, with " frame " for unit outwards sends.The length of every frame is 1500bit, has broadcast a frame and has needed 30s.Every frame comprises 5 subframes, and each subframe all comprises 300bit, and the time of broadcasting is 6s.Each subframe can be divided into again 10 words, and each word is made up of 30bit.First, second and third subframe of navigation message is the ephemeris parameter of satellite, and every 30s repeats once, and its content upgrades once every 2h.Fourth, fifth subframe is the almanac data of satellite, and every 30s overturns one page, and its content includes 25 pages, so the navigation message sending complete set needs 12.5min, and then repeats.
The Z counting of leading character, subframe is comprised in gps navigation message bit stream information, and the content such as the ephemeris of satellite, almanac data.Bit stream is synchronously searched for, until find leading character " 10001011b ".Then according to checking algorithm, information is verified, if not by verification, abandon this data; If by verification, then current data is effective, according to navigation message form, can parse Ionospheric Parameters, ephemeris and almanac data etc.
Extraction and the GPS of BD2 navigation message are similar, also comprise and separate frame structure and information analysis.
Observed quantity extraction module 13.
Observed quantity extraction module 13, obtains the observed quantity such as pseudorange, Doppler, carrier phase of tracked satellite according to Base-Band Processing data.
Pseudo range observed quantity is utilize the transmission time of satellite-signal to be multiplied by transmission range that the light velocity obtains signal.Here transmission time refers to satellite-signal transmitting time t
svdeduct receiver local zone time t
u.Local zone time t
udirectly provided by receiver local clock, and satellite-signal transmitting time t
svneed the counting by subframe or word string, and could accurately calculate after providing local pseudo-code phase by code tracking loop.
The extraction of Doppler measurements is after the frequency values obtaining carrier wave NCO, by the nominal intermediate frequency value f of itself and radio-frequency front-end
ccompare, its difference is exactly Doppler shift observed quantity.
Positioning calculation module 14.
Positioning calculation module 14, carries out co-located according to observed quantity information and the navigation message obtained, and utilizes filtering algorithm to make positioning result more reliable and more stable, meet high sensitivity demand.
Positioning calculation module 14 specifically comprises:
Computing module, for going out the position of satellite according to the ephemeris computation of described satellite;
Positioning calculation module 14, to position customer location for utilizing the position of described satellite and described observed quantity and resolves.
Lower mask body introduces the concrete enforcement principle of computing module, and its distance resolving satellite is carried out step by step, in addition, and can the distance of resolving GPS satellite or the distance of BD2 (a kind of positioning system) satellite.
1) calculating of gps satellite locus
The first step: the mean angular velocity n calculating satellite motion.
First, according to the parameter provided in broadcast ephemeris
calculate with reference to moment t
oemean angular velocity n
0:
in formula, GM is the product of gravitational constant G and earth gross mass M, and its value is GM=3.9860047 × 1014m
3/ s
2.Then according to the mean angular velocity n:n=n of perturbation parameter △ n calculating observation moment satellite given in broadcast ephemeris
0+ △ n.
Second step: the mean anomaly M of calculating observation moment satellite
s.
The cycle of operation due to satellite is 12 hours, adopts time of satellite perigee passing t
0when calculating mean anomaly M, extrapolating, interval is maximum reaches 6 hours.And broadcast ephemeris upgrades once for every 2 hours, with reference to when being constantly located at the central moment, extrapolation interval is less than or equal to 1 hour.So use t
oereplace time of satellite perigee passing t
0after, extrapolation interval will reduce greatly, also can obtain the higher result of precision with better simply model.
M
s=M
0+n(t-t
oe)
In formula, M
0for reference moment t
oetime mean anomaly, provided by broadcast ephemeris.
3rd step: calculate eccentric anomaly.
The Kepler's equation represented with radian is: E
s=M
s+ e
ssinE
s.
4th step: calculate true anomaly.
According to Keplerian orbit equation, anomaly f can be obtained
swith eccentric anomaly E
sbetween relation:
in formula, e
sfor the excentricity of satellite orbit, provided by broadcast ephemeris.True anomaly can be obtained thus and calculate common-used formula:
5th step: calculate to rise and hand over elongation
Utilize formula
calculate to rise and hand over elongation, in formula, ω is the argument of perigee, is provided by broadcast ephemeris.
6th step: calculate perturbation correction member.
C is given in broadcast ephemeris
uc, C
us, C
rc, C
rs, C
ic, C
is6 perturbation parameters, can obtain the liter caused due to the second order zonal harmonic coefficient item of earth gravitational field bit function accordingly and hand over elongation
perturbation correction member δ
μ, satellite radius vector r perturbation correction member δ
r, and the perturbation correction member δ of inclination of satellite orbit i
i.Computing formula is as follows:
7th step: to a liter friendship elongation
satellite radius vector r, orbit inclination i carry out perturbation and correct.
r=a
s(1-e
scosE
s)+δ
r,
In formula: a
sfor the major radius of satellite orbit, i
0for t
oethe orbit inclination in moment, is provided by the kepler parameters in broadcast ephemeris.
for the rate of change of i, provided by the perturbation parameter in broadcast ephemeris.
8th step: calculate the position of satellite in track areal coordinate system.
In track plane right-angle coordinate, the planimetric rectangular coordinates of (true origin is positioned at the earth's core, and X-axis points to ascending node) satellite is:
x
p=rcosμ
y
p=rsinμ
9th step: the longitude L of calculating observation moment ascending node.
If with reference to moment t
oetime ascending node right ascension be
ascending node is △ Ω to the rate of change of time, and so the right ascension of ascending node Ω of observation moment t should be:
In formula, △ Ω can provide from the perturbation parameter of broadcast ephemeris.If start time this week (when Sun. 0) Greenwich sidereal time is GASTweek, then observe moment Greenwich sidereal time be:
GAST=GAST
week+ω
et,
ω in formula
efor rotational-angular velocity of the earth, its value is ω
e=7.292115 × l0-5rad/s, t are the time (s) in this week.The longitude so just can trying to achieve observation moment ascending node is:
Order
Then have:
L=Ω
0+△Ω(t-t
oe)-ω
et=Ω
0+(△Ω-ω
e)t-△Ωt
oe。
Tenth step: calculate the position of satellite in spherical coordinate system instantaneously.
After the geodetic longitude L of known ascending node and the inclination angle i of orbit plane, just try to achieve the position of satellite in body-fixed coordinate system by twice rotation:
11 step: calculate the position of satellite in conventional terrestrial coordinate system.
The observation moment position of satellite in conventional terrestrial coordinate system
The form of BD2 satellite broadcasting ephemeris and the broadcast ephemeris of gps satellite substantially similar, therefore the calculating of its satellite spatial position is also similar with upper joint, repeats no more here.
Positioning calculation module 14.
Positioning calculation module 14, uses satellite position and pseudo range observed quantity, is calculated the position of receiver current time by solving equation.If the position of known satellite i is r
i(x
i, y
i, z
i), the pseudorange of this satellite is p
iif receiver location is r
u(x
u, y
u, z
u), receiver clock-offsets is b
u.For bimodulus positioning system, system time difference △ T also as a unknown quantity, then can obtain measuring equation:
For bimodulus positioning system, the position of known at least five satellites and pseudorange, substitute into above formula and do linearization, can combine and solve, obtain r
u(x
u, y
u, z
u) and the value of receiver clock-offsets.
Based on above-mentioned location ultimate principle, in software program, least square and kalman filter method is used to carry out solving of receiver location and speed.Then, positioning result is detected.If by this detection, then export this positioning result, witness marker puts 1, otherwise does not export this positioning result, and witness marker sets to 0.
Least square method is the basic skills of positioning calculation, and it can try to achieve optimum solution containing between error and each measured value of noise, makes the sum of square of deviations of all measured values minimum.
Meet at least four stars under single system station-keeping mode, under dual system station-keeping mode when at least five stars, according to the method that GPS and BD2 interface control document provides, utilize ephemeris parameter to calculate position and the speed of participation position location satellite.Owing to comprising satellite clock correction and ionosphere, tropospheric delay error in pseudorange, ephemeris parameter before positioning calculation, is first utilized to calculate satellite clock correction; For GPS and BD2 satellite, all ionosphere delay can be calculated according to ionosphere delay parameter; And estimate according to correlation model tropospheric delay.These correcting values are utilized to correct pseudorange afterwards.After being rejected by pseudorange larger for measuring error, Newton iteration method is used to solve, until iteration terminates.
Kalman filtering algorithm, it is a kind of Linear Minimum Variance Filter method, it adopts state equation and measurement equation to describe stochastic linear system, according to the criterion that the variance of estimated state error is minimum, the measuring value from measured noise pollution, estimate each state in each moment of system in real time.It comprises prediction and revises two processes.Kalman filter adopts recursive form in computing method, is easy to realize on computers.
If the system state equation after discretize and measurement equation are respectively:
In formula:
X
krepresent k moment (n × 1) dimension system state vector; X
k-1it is then the system state vector in k-1 moment; Φ
k|k-1for the transition matrix of etching system during k-1; W
k-1the k-1 moment (n × 1) maintain system noise vector; Γ
k-1that system noise drives matrix; Z
krepresent that (m × 1) dimension in k moment measures vector; H
kfor (m × n) in k moment ties up observing matrix; V
kit is k moment (m × 1) dimension measurement noise vector.
Meanwhile, W
k-1and V
kmutual incoherent zero-mean white noise sequence, that is:
E[W
k-1]=0,
E[V
k]=0,
Wherein, Q
kfor system noise variance matrix, be nonnegative definite, R
kfor measurement noises variance matrix, it is positively definite matrix; δ
kjkroneker function:
If discrete system meets above-mentioned requirements, the state optimization of discrete system is estimated to adopt following steps to try to achieve:
State one-step prediction:
State estimation:
Filter gain:
One-step prediction square error:
Estimate square error:
Output module 15, for navigation results described in the standard format NMEA formatted output formulated with electronic equipment according to sea.
NMEA form:
1)GGA-PsitionData。
Functional description: output statement.Locator data is described.This statement comprise locate to receiver, data relevant when surveying.If only BD, gps satellite are used for location compute, transmitting identifier is BD, GP, if employ dual system satellite to obtain location compute, transmits identifier GN.
Form:
$--GGA,hhmmss.ss,llll.ll,a,yyyyy.yy,a,x,xx,x.x,x.x,U,x.x,U,xxxx,x.x,x.x*hh<CR><LF>。
Table 1
Numbering | Implication | Span | Unit | Remarks |
1 | Location moment (UTC time) | |||
2 | Latitude | |||
3 | Latitude direction | N/S | N-north latitude, S-south latitude | |
4 | Longitude | |||
5 | Longitudinal | E/W | E-east longitude, W-west longitude | |
6 | State indicates | 0-8 | Note 1 | |
7 | Participate in the satellite number of location | |||
8 | HDOP value | |||
9 | Antenna geodetic height | |||
10 | Antenna geodetic height unit | Rice | ||
11 | Height anomaly (CGS2000) | |||
12 | Height anomaly unit | Rice | ||
13 | Differential data length of time | |||
14 | No. ID, difference platform | |||
15 | VDOP value |
Note 1: state instruction (this data field can not be empty)
(1) when this statement identifier is GP, state indicates: 0-station-keeping mode is unavailable or invalid; 1-GPSSPS pattern, location is effectively; 2-differential GPS SPS pattern, location is effectively; 3-GPSPPS pattern, location is effectively; In real time dynamically (RTK), system is in RTK pattern 4-, has fixing complete cycle number; The RTK that 5-floats, system is in RTK pattern, and complete cycle number is floating; 6-estimation mode (dead reckoning); 7-manual input mode; 8-simulator pattern.
(2) when this statement identifier is BD, state indicates: 0-location is unavailable or invalid; 1-is without Differential positioning, and location effectively; 2-Differential positioning, location is effectively; 3-double frequency is located, and location effectively.
(3) when this statement identifier is GN, state indicates: 0-location is unavailable or invalid; 1-compatible positioning, location effectively.
(4), time without positioning result, locating information field is empty.
2)RMC-MinimumSpecificData。
Functional description: the simplest navigation transmission data.Form:
$--RMC,hhmmss.ss,A,llll.ll,a,yyyyy.yy,a,x.x,x.x,ddmmyy,x.x,a,a*hh<CR><LF>。
Table 2
3)GSV-VisibleSatellitesData。
Functional description: output statement.This statement comprises visual satellite number, satellite identification number, the elevation angle, position angle and signal to noise ratio (S/N ratio) (C/N0) value.Each transmission mostly is 4 satellites most, and the statement sum of transmission and the statement number of transmission show in the first two field.If many BD, gps satellite can be seen, use GSV statement respectively, identify the BD satellite seen with identifier BD, represent with GP and see gps satellite.GN identifier does not use together with this statement.
Form: $--GSV, x, x, xx, xx, xx, xxx, x.x ... * hh<CR><LFGreatT. GreaT.GT.
Table 3
Numbering | Implication | Span | Unit | Remarks |
1 | GSV statement sum | |||
2 | Current GSV statement sequence number | |||
3 | Satellite number in the visual field | |||
4 | Satellite number | |||
5 | Satellite elevation angle | Degree | ||
6 | Satellite aximuth | Degree | ||
7 | Signal to noise ratio (S/N ratio) | dB-Hz | ||
…… | Repeat 4 ~ 7 fields | Other satellite information |
4)VTG-ForwardDirectionandGroundSpeedData。
Functional description: this statement comprises actual heading relative to ground and speed.If only BD, gps satellite are used for location compute, transmitting identifier is BD, GP, if employ dual system satellite to obtain location compute, transmits identifier GN.
Form: $ GPVTG,
xxx.x,
t,
xxx.x,
m, xxx.x,
n,
xxx.x,
k* hh<CR><LFGreatT. GreaT.GT.
Table 4
Numbering | Implication | Span | Unit | Remarks |
1,2 | Course over ground | Degree | Geographical north | |
3,4 | Course over ground | Degree | Magnetic north | |
5,6 | Ground speed | Joint | 10 --> | |
7,8 | Ground speed | Thousand ms/h | ||
5 | Satellite elevation angle | Degree | ||
6 | Satellite aximuth | Degree | ||
7 | Signal to noise ratio (S/N ratio) | dB-Hz | ||
…… | Repeat 4 ~ 7 fields | Other satellite information |
5)GLL-GeographicTimeandPositionData。
Functional description: output statement.Terrestrial coordinate locating information, the latitude of carrier, longitude, positioning time and state.If only BD, gps satellite are used for location compute, transmitting identifier is BD, GP, if employ dual system satellite to obtain location compute, transmits identifier GN.
Form: $--GLL, llll.ll, a, yyyyy.yy, a, hhmmss.ss, A, x*hh<CR><LFGreat T.GreaT.GT.
Table 5
Note 1: pattern indicates: 0-automatic mode; 1-difference modes; 2-estimates (dead reckoning) pattern; 3-manual input mode; 4-simulator pattern.
6)GSA-DOPandVisibleSatellitesData。
Functional description: output statement.This statement comprises the satellite and dilution of precision (DOP) value used in the navigation calculation of subscriber equipment mode of operation, the report of GGA statement.When only resolving position by BD, gps satellite system, representing with identifier BD respectively and transmitting BD satellite dilution of precision and effective satellite number, representing with GP and transmit gps satellite dilution of precision and effective satellite number.When integrated use BD, GPS are to obtain location compute, many GSA statements can be produced, each GSA statement uses BD, GP as identifier respectively, to distinguish in data combination the satellite of the different satellite systems used, and every bar statement has PDOP, HDOP, VDOP and the TDOP of the combination satellite system for location compute.GN identifier does not use together with this statement.
Form: $--GSA, a, x, xx ..., xx, x.x, x.x, x.x*hh<CR><LFGre atT.GreaT.GT.
Table 6
Note 1: pattern indicates:
M-is manual, forces to be used for 2D or 3D pattern;
A-is automatic, allows 2D/3D auto-changing.
Note 2: select pattern:
1-location is unavailable or invalid;
2-2D;
3-3D。
7)ZDA-DataandTime
Functional description: output statement.Terrestrial coordinate locating information, the latitude of carrier, longitude, positioning time and state.If only BD, gps satellite are used for location compute, transmitting identifier is BD, GP, if employ dual system satellite to obtain location compute, transmits identifier GN.
Form: $--ZDA,
hhmmss.ss,
dd, mm, yyyy, xx, yy*CC.
Table 7
Numbering | Implication | Span | Unit | Remarks |
1 | Time (UTC) | Hour Minute Second | Note 1 | |
2 | Time (UTC) | Day month year | Note 2 | |
3 | Local time zone hour | -13~13 | Hour | |
4 | Local time zone minute | 0~59 | Minute |
By one or more embodiment of the present invention, the present invention has following beneficial effect or advantage:
The invention provides a kind of device for navigating, comprising: acquisition and tracking processing module, message extraction module, observed quantity extraction module, positioning calculation module, output module.Carry out continuing to follow the tracks of the navigation message bit data obtaining satellite to satellite-signal by acquisition and tracking processing module; By message extraction module, observed quantity extraction module, obtain satellite navigation message and observed quantity, then positioned by positioning calculation module and resolve, and exported by the standard format NMEA form that output module is formulated according to sea electronic equipment, the navigation message bit data of acquisition satellite that can be automatic continuous position, and then avoid occurring that navigation situation not in time occurs, and have technique effect promptly and accurately.
Although described the preferred embodiment of the application, one of ordinary skilled in the art once obtain the basic creative concept of cicada, then can make other change and amendment to these embodiments.So claims are intended to be interpreted as comprising preferred embodiment and falling into all changes and the amendment of the application's scope.
Obviously, those skilled in the art can carry out various change and modification to the application and not depart from the spirit and scope of the application.Like this, if these amendments of the application and modification belong within the scope of the application's claim and equivalent technologies thereof, then the application is also intended to comprise these change and modification.
Claims (10)
1. the device for navigating, is characterized in that, comprising:
Acquisition and tracking processing module, continues to follow the tracks of for carrying out satellite-signal, and obtains the navigation message bit data of satellite;
Message extraction module, for receiving described navigation message bit data, and goes out satellite navigation message from described navigation message bit extracting data;
Observed quantity extraction module, for going out observed quantity from described navigation message bit extracting data;
Positioning calculation module, for carrying out co-located according to described observed quantity and described satellite navigation message, obtains navigation results;
Output module, for navigation results described in the standard format NMEA formatted output formulated with electronic equipment according to sea.
2. device as claimed in claim 1, it is characterized in that, described satellite navigation message comprises ephemeris, satellite health, time device, the satellite almanac of satellite.
3. device as claimed in claim 1, is characterized in that, described observed quantity extraction module, specifically for obtaining the pseudo range observed quantity of tracked satellite or Doppler measurements or carrier phase observed quantity according to Base-Band Processing data.
4. device as claimed in claim 3, is characterized in that, described pseudo range observed quantity is utilize the transmission time of satellite-signal to be multiplied by transmission range that the light velocity obtains signal; The transmission time of described satellite-signal refers to satellite-signal transmitting time t
svdeduct receiver local zone time t
u.
5. device as claimed in claim 3, it is characterized in that, described Doppler measurements is after the frequency values obtaining carrier wave NCO, by the frequency values of described carrier wave NCO and the nominal intermediate frequency value f of radio-frequency front-end
csubtract each other acquisition.
6. device as claimed in claim 2, it is characterized in that, described positioning calculation module specifically comprises:
Computing module, for going out the position of satellite according to the ephemeris computation of described satellite;
Positioning calculation module, to position customer location for utilizing the position of described satellite and described observed quantity and resolves.
7. device as claimed in claim 6, is characterized in that, described computing module specifically for:
The first step: the mean angular velocity n calculating satellite motion;
Second step: the mean anomaly M of calculating observation moment satellite
s;
M
s=M
0+ n (t-t
oe); In formula, M
0for reference moment t
oetime mean anomaly, provided by broadcast ephemeris;
3rd step: calculate eccentric anomaly: E
s=M
s+ e
ssinE
s;
4th step: calculate anomaly f
s:
5th step: calculate to rise and hand over elongation
in formula, ω is the argument of perigee, is provided by broadcast ephemeris;
6th step: calculate perturbation correction member;
7th step: to a liter friendship elongation
satellite radius vector r, orbit inclination i carry out perturbation and correct;
8th step: calculate the position of satellite in track areal coordinate system;
9th step: the longitude L of calculating observation moment ascending node;
Tenth step: calculate the position of satellite in spherical coordinate system instantaneously;
11 step: calculate the position of satellite in conventional terrestrial coordinate system.
8. device as claimed in claim 7, it is characterized in that, described mean angular velocity n is obtained by following steps:
First according to the parameter provided in broadcast ephemeris
calculate with reference to moment t
oemean angular velocity n
0:
in formula, GM is the product of gravitational constant G and earth gross mass M, GM=3.9860047 × 1014m
3/ s
2, then according to the mean angular velocity n:n=n of perturbation parameter △ n calculating observation moment satellite given in broadcast ephemeris
0+ △ n.
9. device as claimed in claim 7, is characterized in that, described calculating perturbation correction member comprises the following steps:
By the C provided in broadcast ephemeris
uc, C
us, C
rc, C
rs, C
ic, C
is6 perturbation parameters, obtain the perturbation correction member δ that the liter caused due to the second order zonal harmonic coefficient item of earth gravitational field bit function hands over elongation μ
μ, satellite radius vector r perturbation correction member δ
r, and the perturbation correction member δ of inclination of satellite orbit i
i, computing formula is as follows:
10. device as claimed in claim 7, is characterized in that, described to a liter friendship elongation
the step that satellite radius vector r, orbit inclination i carry out perturbation correction is:
R=a
s(1-e
scosE
s)+δ
r, in formula: a
sfor the major radius of satellite orbit, i
0for t
oethe track in moment inclines
Angle, is provided by the kepler parameters in broadcast ephemeris;
for the rate of change of i, provided by the perturbation parameter in broadcast ephemeris.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937376.XA CN105510935A (en) | 2015-12-15 | 2015-12-15 | Device used for navigation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937376.XA CN105510935A (en) | 2015-12-15 | 2015-12-15 | Device used for navigation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105510935A true CN105510935A (en) | 2016-04-20 |
Family
ID=55719044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510937376.XA Pending CN105510935A (en) | 2015-12-15 | 2015-12-15 | Device used for navigation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510935A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107643527A (en) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | GPS regards signal imitation generation method altogether and GPS regards signal simulator altogether |
CN108267765A (en) * | 2018-03-13 | 2018-07-10 | 北京沙谷科技有限责任公司 | Use the MEMS integrated navigation modeling methods of imperfect GNSS information |
CN108415044A (en) * | 2018-02-09 | 2018-08-17 | 广州吉欧电子科技有限公司 | A kind of processing system, processing method and GNSS receiving devices |
CN109085627A (en) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | A method of improving Global Navigation Satellite System single frequency receiving positioning accuracy |
CN112255648A (en) * | 2020-10-20 | 2021-01-22 | 腾讯科技(深圳)有限公司 | Motion state detection method and device, electronic equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101132191A (en) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | Baseband signal processing method for GNSS receiver |
CN101592723A (en) * | 2008-05-30 | 2009-12-02 | 凹凸电子(武汉)有限公司 | GPS receiver and localization method thereof |
CN101609140A (en) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | A kind of compatible navigation receiver positioning system and localization method thereof |
-
2015
- 2015-12-15 CN CN201510937376.XA patent/CN105510935A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101132191A (en) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | Baseband signal processing method for GNSS receiver |
CN101592723A (en) * | 2008-05-30 | 2009-12-02 | 凹凸电子(武汉)有限公司 | GPS receiver and localization method thereof |
CN101609140A (en) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | A kind of compatible navigation receiver positioning system and localization method thereof |
Non-Patent Citations (2)
Title |
---|
工测顽童: "GPS卫星的坐标计算", 《百度文库》 * |
谢钢: "《GPS原理与接收机设计》", 31 July 2009 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107643527A (en) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | GPS regards signal imitation generation method altogether and GPS regards signal simulator altogether |
CN108415044A (en) * | 2018-02-09 | 2018-08-17 | 广州吉欧电子科技有限公司 | A kind of processing system, processing method and GNSS receiving devices |
CN108415044B (en) * | 2018-02-09 | 2020-09-15 | 广州吉欧电子科技有限公司 | Processing system, processing method and GNSS receiving equipment |
CN108267765A (en) * | 2018-03-13 | 2018-07-10 | 北京沙谷科技有限责任公司 | Use the MEMS integrated navigation modeling methods of imperfect GNSS information |
CN108267765B (en) * | 2018-03-13 | 2021-08-03 | 北京沙谷科技有限责任公司 | MEMS combined navigation modeling method using incomplete GNSS information |
CN109085627A (en) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | A method of improving Global Navigation Satellite System single frequency receiving positioning accuracy |
CN109085627B (en) * | 2018-06-11 | 2019-07-12 | 北京讯腾智慧科技股份有限公司 | A method of improving Global Navigation Satellite System single frequency receiving positioning accuracy |
CN112255648A (en) * | 2020-10-20 | 2021-01-22 | 腾讯科技(深圳)有限公司 | Motion state detection method and device, electronic equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Skog et al. | In-car positioning and navigation technologies—A survey | |
US8788200B2 (en) | Method and system for a data interface for aiding a satellite positioning system receiver | |
US7940210B2 (en) | Integrity of differential GPS corrections in navigation devices using military type GPS receivers | |
CN105510935A (en) | Device used for navigation | |
CN100437142C (en) | Error separation method based on foundation strength system and foundation strength system | |
CN101950025B (en) | Data quality monitoring method for local enhancing system | |
CN105044747B (en) | Time synchronization device and method based on multi-satellite common view and filtering | |
CN105467415B (en) | A kind of SUAV RTK relative positioning methods constrained based on difference pressure altitude | |
CN101395443B (en) | Hybrid positioning method and device | |
CN102253399B (en) | Doppler differential compensation velocity measurement method utilizing carrier phase central value | |
CN101833080A (en) | Method for measuring attitude of carrier by using additional constraint condition of GPS system | |
Johnson et al. | Where in the world are my field plots? Using GPS effectively in environmental field studies | |
CN103837879A (en) | Method for realizing high-precision location based on Big Dipper system civil carrier phase combination | |
Dovis et al. | Recent advancement on the use of global navigation satellite system-based positioning for intelligent transport systems [guest editorial] | |
CN106707322A (en) | RTK/SINS-based high dynamic positioning attitude-determining system and method | |
Bhardwaj et al. | Satellite navigation and sources of errors in positioning: a review | |
CN105510946A (en) | Quick resolving method for integer ambiguity of BDS satellite carrier phase | |
CN105572705A (en) | Navigation device | |
Wells et al. | NAVSTAR performance analysis | |
CN105510936A (en) | Satellite-borne GNSS combined orbit determination method | |
CN105572706A (en) | Automatic navigation driving system | |
CN103645488B (en) | System and method for determining GPS receiver position | |
Yu | Improved positioning of land vehicle in ITS using digital map and other accessory information | |
CN108351423A (en) | Railway localization method and system | |
CN106199655A (en) | A kind of vector tracking method based on federated 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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160420 |