CN100405012C - Satellite direction finder and direction finding method - Google Patents
Satellite direction finder and direction finding method Download PDFInfo
- Publication number
- CN100405012C CN100405012C CNB2003101102396A CN200310110239A CN100405012C CN 100405012 C CN100405012 C CN 100405012C CN B2003101102396 A CNB2003101102396 A CN B2003101102396A CN 200310110239 A CN200310110239 A CN 200310110239A CN 100405012 C CN100405012 C CN 100405012C
- Authority
- CN
- China
- Prior art keywords
- partiald
- satellite
- phi
- platform
- rho
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The present invention relates to a satellite direction detector, in which a rotation platform is installed on a supporting frame, with satellite carrier phase measurement and receiver antennae installed on both ends of an arm of the rotation platform. The rotation platform is provided with a telescope, a level rotation angle electronic recorder, a vertical angel electronic recorder and a wave carrier phase measuring receiver. The synchronous sampling outputs of the devices installed on the rotation platform are all connected with a data processing unit. The direction detecting method is driving the rotation platform to rotate firstly. Then the two antennae measure wave carrier phase and give out second impulses at the rate of once a second, and meanwhile the rotation platform records rotation angles according the second impulses. Then the data processing unit processes data according to a drafted math model capable of weakening system errors, and works out the azimuth angle of a dial plate zero position. Thus, when the telescope points to any direction, the azimuth angle of the direction can be read out from the dial plate. The present invention can be used for any satellite navigation system, including a geosynchronous fixed-point satellite navigation system used for azimuth determination. The present invention can stably determine directions using only one satellite, and can provide verification and enhance direction determination accuracy using two or more satellites.
Description
Technical field
The present invention relates to a kind of bearing transmission unit and assay method, particularly relate to a kind of device and direction-finding method that utilizes satellite to carry out direction-finding.
Background technology
Existing satellite navigation system as GPS of USA, Muscovite GLONASS and geostationary satellite navigational system, has location (providing the on-site coordinate of user) function, but do not have orientation (fixing the position) function in system design.
Except the needs location, also need orientation a large amount of in dual-use.Usually the user need provide position coordinate and orientation simultaneously.Utilize the Navsat location and fix the position and to satisfy this requirement, and have round-the-clock, quick and high-precision characteristics.
Bring into use gps satellite to carry out direction-finding (I at first propose) from China in 1987; Its requirement is the GPS receiver to be set two stations, the carrier phase of the observation satellite that (contains 4) more than 4, and observation certain hour (for example 0.5~1 hour).The ultimate principle of this method is to utilize 4 above satellites to carry out two stations three-dimensional relative positioning, the geodetic azimuth of reduction two station lines (abbreviation limit) again.Because there is the phase place complete cycle number (blur level) of observation for the first time in the carrier phase observation to satellite, must be no less than 4 satellites resolves and requires to observe certain hour to change the apparent direction of satellite to wait for satellite motion simultaneously, thereby the change equation coefficient can be separated system of equations.This method has obtained widespread use.
Using any other satellite navigation system (comprising the geostationary satellite navigational system) and necessary moonscope decreased number to is carried out orientation, is the important technology approach of breaking away from the directed restriction of present applying GPS.And synchronous satellite navigation system has only two fixed point geostationary satellites, and it does not satisfy the basic demand of 4 above satellites of observation; The apparent direction of satellite is almost constant in observation time in addition, can not resolve blur level (listed system of equations rank defect).That is to say and continue to use existing theory and method can not obtain oriented solution.
Summary of the invention
The purpose of this invention is to provide and a kind ofly utilize single satellite (comprising earth synchronous and fixed communications satellite) can measure the satellite goniometer of target azimuth and utilize this goniometer to carry out the method for direction finding.
For achieving the above object, the technical solution used in the present invention is: a kind of satellite goniometer, being equipped with on support can be by the rotation platform of drive motor driven in rotation, rotation platform comprises the support arm that stretch out at two ends, satellite carrier phase measurement receiver antenna is equipped with at the support arm two ends, telescope, level angle electronic logger, vertical angle electronic logger and carrier phase measurement receiver are housed on the rotation platform, and the synchronized sampling output of carrier phase measurement receiver, horizontal angle electronic logger and vertical angle electronic logger all is connected with data processing unit.
Platform brachium 0.75-2m, platform vertical stability 5 ", platform rotating speed 1 cycle per minute, horizontal angle registration accuracy 2 " and, vertical angle registration accuracy 3 ", three quadrature bias 5 " and, collimation telescope precision 3 ".
The port number of two antennas is 4, adopts common local oscillator, phase measurement accuracy≤0.01 week, fundamental frequency degree of stability 10
-8, fundamental frequency accuracy 10
-6, each channel sample postpones≤100 μ s, and two antenna channels sampling delay are poor≤10 μ s.
The method that the present invention utilizes described satellite goniometer to carry out direction finding is: the rotation of driven in rotation platform, make two antennas carry out carrier phase measurement and provide pulse per second (PPS) by 1 time/second, rotation platform carries out the corner record synchronously by pulse per second (PPS) simultaneously, data processing unit is pressed the many groups phase place behind the cycle slip editor of being gathered, platform corner data and in conjunction with known satellite and survey station coordinate data, carry out data processing according to the mathematical model of drafting that can slacken systematic error, calculate the position angle of scale zero-bit, after this, any direction of pointing of the telescope all can be read the position angle of this direction from scale.
The described azimuthal mathematical model of scale zero-bit of resolving is:
F is a carrier frequency in the formula, and c is the light velocity, D (t
i) be the platform scale corresponding readings that is connected with receiver, D
0For pointing to the circle reading in positive north, dd is the deformation of antenna distance d, a
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch.
The scale zero-bit is azimuthal resolves plumb line deviation simultaneously resolving, and iteration carries out, and resolves the mathematical model that plumb line deviation separates to be:
F is a carrier frequency in the formula, and c is the light velocity, D (t
i) be the platform scale corresponding readings that is connected with receiver, ξ, η are two quadrature component a of plumb line deviation
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch.
The removing method of antenna phase center deviation is to make the rotation in revolution of two antennas, and rotation platform rotates two all antennas from circling.
After adopting such scheme, the present invention has following good effect:
1, the present invention can comprise that earth synchronous and fixed communications satellite navigational system carries out direction-finding with any satellite navigation system.It only uses a satellite just can be stably directed, uses two or two above satellites to provide to check and improves orientation accuracy.
2, satellite goniometer of the present invention does not need any orientation initial value in the direction finding process, only survey station coordinate and the co-ordinates of satellite that needs navigational system to provide.
3, the present invention utilizes the controlled motion of antenna to solve the key issue that the geostationary satellite phase measurement is difficult to resolve blur level, utilizes a large amount of redundancy observation to weaken the influence to orientation of multinomial systematic error simultaneously, improves direction finding precision.
4, the relation of application antenna rotation and revolution can slacken the influence of the tangential component of antenna phase deviation to orientation.
5, the duty of this instrument belongs to low kinetic measurement, because its motion is controlled, laws of motion are strong, can solve the cycle slip problem that may occur easily and reliably.
6, the automaticity height of instrument does not need accurate leveling, easy operating.
7, this instrument has the ability in fast measuring orientation, simulation calculation show when observation time is 15 minutes, (satisfy primary demand) precision be about 5 ".Observation time can further reduce, and just precision reduces, for example 3 minutes 20 seconds the time, and about 17 "; Can satisfy part has relatively high expectations but the not high demand of precision to observation time.Can strengthen arm head and observation time when precision prescribed is higher solves.
8, have telescope, can carry out for example demarcation of inertial navigation system for point instrumentation.
Although 9 satellite goniometers are for using earth synchronous and fixed communications satellite Navigation System Design, its method and apparatus can be used for other equally as gps satellite orientation (using GPS receiver and data processing software instead).Because GPS visible satellite number is many, precision can also further improve.
10, the satellite goniometer also can be used for attitude determination.
The present invention can use the synchronous pinpoint navigation of earth system and carry out orientation, has solved navigational system and can only locate non-orientable problem, satisfies the problem in dual-use middle needs fast measuring orientation.
Description of drawings
Fig. 1 is the structural representation of satellite goniometer of the present invention;
Fig. 2 is a direction finding principle schematic of the present invention;
Fig. 3 is the process flow diagram that resolves of direction-finding method of the present invention.
Fig. 4 is the theory diagram of simulation calculation.
Embodiment
The satellite goniometer is to adopt carrier phase measurement to resolve the specialized equipment in orientation separately.Its emphasis solve to adopt single satellite to resolve the Carrier Phase Ambiguity Resolution problem of orientation and geostationary satellite initial epoch, carries out direction finding to use two earth synchronous and fixed communications satellites.As shown in Figure 1, being equipped with on support 10 can be by the rotation platform 12 of drive unit 9 driven in rotation, satellite carrier phase measurement receiver antenna 5 is equipped with at support arm 1 two ends at rotation platform 12, level angle electronic logger 8, vertical angle electronic logger 4 and carrier phase measurement receiver 6 are housed on the rotation platform 12, and level angle electronic logger 8 and vertical angle electronic logger 4 are used to write down the corner reading.Also be provided with telescope holder 2 on rotation platform 12, telescope 3 is housed on the telescope holder 2, telescope is used for determining the required position angle of sensing, back, positive north.The output of carrier phase measurement receiver 6, horizontal angle electronic logger 8 and vertical angle electronic logger 4 all is connected with data processing unit 7.Accumulator 11 also is housed on rotation platform, is used to each several part that working power is provided.Data processing unit 7 adopts single-chip microcomputer or other similar computing units, the program of cooperation direction-finding method establishment according to the present invention gets final product, the composition of data processing unit and program composition belong to those skilled in the art's conventional knowledge, so the present invention no longer discloses flow chart.
The goniometer the key technical indexes is as follows:
The goniometer the key technical indexes is to draw according to the impact analysis to resultnat accuracy of desired precision and sub-indicator.
2~14 minutes direction finding time (not containing preliminary work)
Direction finding precision 30~10 " (time is long, the precision height)
The rotation platform technical requirement
Platform brachium 0.75-2m (optional) (the long precision height of brachium)
Platform vertical stability 5 "
Platform rotating speed 1 cycle per minute
Horizontal angle registration accuracy 2 "
Vertical angle registration accuracy 3 "
Three quadrature bias 5 "
Telescope (or autocollimation) pointing accuracy 3 "
The technical requirement of phase measurement receiver
Antenna 2
Port number 4
Phase measurement accuracy≤0.01 week
Fundamental frequency degree of stability 10
-8
Each channel sample postpones≤100 μ s
Channel sample postpones poor≤10 μ s in twos
Its direction finding principle as shown in Figure 2, S is a satellite among the figure, P is the survey station of known coordinate, P
1, P
2For placing two carrier phase receiver antennas on the rotation platform, N be the north to, A is the position angle that two receiver lines point to, ρ is the distance of satellite to survey station, d is a brachium.Adopt the horizontal system of coordinates, the X-axis energized north, the Z axle points to local reference ellipsoid normal, and Y-axis is by the right-handed coordinate system definition.
Two receiver antenna carrier phase observed quantities poor that places rotation platform to adopt common local oscillator, promptly
F is a carrier frequency in the formula, and c is the light velocity, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch, A
iBe i P epoch
2P
1The position angle of direction, the platform scale corresponding readings that is connected with receiver is D (t
i), then
D in the formula
0For pointing to the circle reading in positive north.
The precision in scale summary orientation can remain within 1~2 °, and (2) formula can change and turns to
Such equation can have tens of to hundreds of according to rotating speed, sample frequency and the Measuring Time of rotation platform.In the equation, known satellite and survey station coordinate time
Be given value, D (t
i) be given value, unknown parameter to be separated only has N (t
1) and D
0Wherein the coefficient of N is always 1, D
0Coefficient different epoch under the situation that rotation platform constantly rotates bigger relative variation can be arranged, equation can be separated.
In principle, Dan Xingke carries out the orientation and resolves, and two satellites can improve reliability.In addition, although should theory and method be mainly to be to solve the peculiar problem of earth synchronous and fixed communications satellite navigational system, it also is applicable to other navigational system.
Error analysis and model
(1) basic observed quantity: the difference of the phase place that two antennas are surveyed and platform corner reading.
(2) main error source analysis and processing:
(a) antenna phase center deviation
Adopt the poor of two antenna phase values in the data processing, the antenna phase center deviation refers to the poor of the two.
Horizontal component: can represent by plane vector argument da
a, mould da
r
Vertical component: scalar da
h
Two antennas may have different values, and its technical indicator is 5mm.May be different for different satellite (direction).Da wherein
rAnd da
hIn model, consider da as system deviation
tAdopt the antenna rotation to eliminate (seeing 4).
(b) antenna distance error: (with nominal value poor) dd
r, the deformation that causes comprising process and assemble error, Yin Wendu; Consider in model as system deviation.
(c) phase measurement error: the receiver technical indicator was 0.01 week, was 0.014 week for the difference of two antennas.In model, consider as stochastic error.
(d) horizontal angle measuring error: index is 3 "; Consider as stochastic error.
(e) phase measurement and angle measurement time synchronization error: index is 10 μ s, presses the antenna angular velocity of rotation, is 0.2 to the influence of orientation "; Can ignore.
(f) two antenna samples are asynchronous: index is 10 μ s, has different value (because of signal intensity may be different) for two satellites.Consider in model as system deviation.
(g) the asynchronous stability of two antenna samples: index was 0.01 week.Consider in model as system deviation, be assumed to the time and be linear change, δ t
i=a
iT has different value (because of signal intensity may be different) for two satellites.
(h) complete cycle counting: DNN has different value for two satellites.Consider in model as system deviation.
(i) orientation initial value error: required orientation initial value in resolving, solve from observation data, the angle reading of getting the phase differential maximal value correspondence of single star is the reading of satellite direction, calculates satellite aximuth according to this, their difference is the orientation initial value.
(j) attitude and plumb line deviation: the deviation of platform turning axle and coordinate system z axle, wherein two of plumb line deviation quadrature component ξ, η can reach 10 usually " or bigger (domestic about 30 " in), the leveling error that comprises instrument, it is the attitude of platform, allow to reach 300 ", both sums are still represented with ξ, η.The influence that it has identical magnitude to the direction-finding precision is considered to be resolved in model as system deviation; Or say, can calculate the attitude angle of relative normal simultaneously.
The parameter that need resolve
Scale zero point: the circle reading of direct north.
Plumb line deviation or attitude: ξ, η.
Phase place complete cycle counting: N is to different satellite differences.But it has comprised other constant error, as the vertical component of antenna deviation.
Brachium error: actual brachium and nominal value poor.
Antenna phase center deviation (three-component)
The asynchronous stability error of antenna samples: time t once or quadratic polynomial.
Mathematical model
From (3) formula, adopt the difference of two receiver carrier phase observed quantities of common local oscillator to be
N in the formula
12(t
1) be complete cycle number scale number poor of first epoch of two receiving antennas (receiver) constantly, A
iBe i P epoch
2P
1The position angle of direction, D
0The circle reading of direct north, the platform scale corresponding readings that is connected with receiver is D (t
i), then
The deformation dd that considers antenna distance d is then:
Can solve the circle reading D that points to positive north
0With antenna distance deformation dd
Can separate two quadrature component ξ, η of attitude misalignment on this basis
The orientation is separated and is separated iteration with plumb line deviation and carry out in the real work.
For other Systematic error sources that may exist (being determined by test data analyzer), part can be paid attention in data processing, and its mathematical model form is
A in the formula
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative.The elimination of antenna phase center deviation
The antenna phase center deviation belongs to systematic bias, and wherein the component (abbreviation normal direction) on the plane that constitutes perpendicular to antenna connection and normal (vertical line) can't solve or eliminate by data processing; Its influence very big again (the 2mm deviation will cause about 200 when brachium 1m " azimuthal error) must take specific process to be eliminated.The antenna phase center deviation can be decomposed into aforesaid da
r, da
hAnd da
aDa
hComponent is considered in model as system deviation, da
rAnd da
aAdopt following method to eliminate.
Two antennas take the mode of rotation to work, and platform rotation (revolution) 2 all antennas are from circling (rotation period is 2 times of the period of revolution).Rotate in two weeks at platform like this, the antenna condition of each position all has and is 180 ° of reverse states with it, since the time very short (1 minute at interval), da in the antenna phase center deviation
rAnd da
aComponent to the contribution of separating equate, opposite, will automatically eliminate this error in separating this moment.
Because da
hComponent is the component along the antenna turning axle, and with the antenna irrelevant to rotation, each position that the antenna public affairs circle all has and is 180 ° of reverse correspondence positions with it, they to the contribution of separating equate, opposite, will eliminate this error automatically in separating this moment.
Cycle slip editor and scale orientation initial value D
0Calculating
When antenna at the uniform velocity revolved round the sun, two antenna phase differences became varies with sinusoidal function, were not difficult to find cycle slip through data fitting, and corrected the complete cycle numeration deviation that causes because of possible cycle slip, promptly finished the cycle slip editor.This moment, phase differential was varies with sinusoidal function.Calculate the orientation (co-ordinates of satellite and survey station coordinate are known) of the satellite of surveying; Getting the circle reading of maximal value correspondence of the phase differential behind the cycle slip editor and the difference of satellite position is scale orientation initial value.Press simulation result, initial value can deviation 9 °, asks fixed method to satisfy accuracy requirement like this.
The method that the present invention utilizes the satellite goniometer to carry out direction finding is: the rotation of driven in rotation platform, make two antennas carry out carrier phase measurement and provide pulse per second (PPS) by 1 time/second, rotation platform carries out the corner record synchronously by pulse per second (PPS) simultaneously, data processing unit is pressed the many groups phase place behind the cycle slip editor of being gathered, platform corner data and in conjunction with known satellite and survey station coordinate data, carry out data processing according to the mathematical model of drafting that can slacken systematic error, calculate the position angle of scale zero-bit, after this, any direction of pointing of the telescope all can be read the position angle of this direction from scale.
It resolves flow process as shown in Figure 3, and wherein satellite position, survey station position provide for satellite navigation system, after coordinate transform to orientation and attitude algorithm; Phase observations value and platform corner value are respectively phase receivers and platform automatic corner recording-member provides, and the phase observations value is carried out initial value behind the editor of cycle slip definite, and transports to orientation and attitude algorithm in the lump with the platform corner; After solving orientation and attitude, display orientation and attitude, and with the output of universal data interface agreement, in order to being connected with miscellaneous equipment.
Simulation calculation
For the stability of the correctness of verifying mathematical model, system of equations and additional parameter to the effect of impair system sexual deviation, carried out comprising that constringency performance, elementary solution, individual event and the multinomial systematic error of iteration comprehensively resolved.What below provide is the Computer Simulation that comprises that antenna phase deviation, plumb line deviation (attitude), brachium deformation, phase observations amount fuzzy parameter comprehensively resolve.
Emulated data generation and orientation are resolved by the software of two independent operatings and are finished respectively.The emulated data that produces provides with the form of file, and the orientation is resolved software and read in file data (simulated measurement) and resolve.
Raw data file comprises:
Co-ordinates of satellite, survey station coordinate, phase observations value and corresponding circle reading thereof.
Comprise in the emulated data:
Fuzzy parameter (hundred all magnitudes), parallax range error (5mm magnitude), antenna phase center deviation (5mm magnitude), 0.01 week of phase observations amount stochastic error, instrument attitude misalignment (300 " magnitude).Simulated conditions: observation time: 14 minutes, systematic bias was set to:
Brachium roll-call latitude longitude error orientation initial value brachium error ξ η
1.00 C30015D 30.0 100.0 .014 2.000 .005 300. -300.
Blur level N
1Blur level N
2Postpone
1Postpone
2Da
R1Da
R2Da
A1Da
A2Da
H1Da
H2
10.0 20.0 1.000 2.000 .005 .003 .110 1.440 .004 .002
Wherein subscript 1,2 is represented two antennas respectively.
(1) emulation of no errors of measurements
At first do not have observation (at random) error emulation (having whole systematic bias), with the correctness of proof theory and mathematical model, it the results are shown in Table 1
Table 1 no errors of measurements simulation result
Scale orientation at zero point (°) | Azimuthal error (") | Attitude ξ (") | ξ error (") | Attitude η (") | η error (") |
-2.0000 | 0.06 | 300.41 | 0.09 | -300.94 | 0.06 |
The scale orientation at zero point that solves and the difference of attitude value and setting value are 0.1 " in, illustration method, model and simulation software are correct.
(2) emulation of observational error is arranged
Observation (at random) error is arranged, carry out 30 (different random error) emulation altogether, the results are shown in Table 2
The precision statistics of 30 orientation of table 2 and attitude algorithm
Scale orientation at zero point (°) | Azimuthal error (") | Attitude ξ (") | ξ error (") | Attitude η (") | η error (") |
-1.999366 | 2.281560 | 296.345940 | -3.654060 | -301.685572 | -1.685572 |
-2.001437 | -5.174682 | 304.715953 | 4.715953 | -287.139541 | 12.860459 |
-2.000141 | -.507404 | 293.490322 | -6.509678 | -296.257756 | 3.742244 |
-2.001318 | -4.746476 | 303.728986 | 3.728986 | -284.360341 | 15.639659 |
-1.999158 | 3.032143 | 292.483685 | -7.516315 | -307.256073 | -7.256073 |
-1.999999 | .003292 | 303.712022 | 3.712022 | -306.419107 | -6.419107 |
-2.003391 | -12.207758 | 288.877216 | -11.122784 | -301.412243 | -1.412243 |
-1.998536 | 5.271786 | 303.287024 | 3.287024 | -323.258591 | -23.258591 |
-2.001264 | -4.550459 | 289.618218 | -10.381782 | -311.257296 | -11.257296 |
-2.000279 | -1.005353 | 290.265214 | -9.734786 | -283.636863 | 16.363137 |
-1.999839 | .578325 | 297.750507 | -2.249493 | -276.865259 | 23.134741 |
-2.000505 | -1.818004 | 304.108587 | 4.108587 | -282.325367 | 17.674633 |
-2.002048 | -7.371408 | 299.211752 | -.788248 | -305.702332 | -5.702332 |
-1.998548 | 5.227831 | 295.368320 | -4.631680 | -311.227274 | -11.227274 |
-2.001186 | -4.270437 | 304.366264 | 4.366264 | -297.740244 | 2.259756 |
-1.999739 | .938479 | 312.068257 | 12.068257 | -287.317450 | 12.682550 |
-2.001772 | -6.378321 | 292.862821 | -7.137179 | -290.199219 | 9.800781 |
-2.000079 | -.284429 | 290.839434 | -9.160566 | -320.024362 | -20.024362 |
-1.997920 | 7.487009 | 312.167344 | 12.167344 | -313.311746 | -13.311746 |
-2.000619 | -2.229697 | 295.494287 | -4.505713 | -298.460754 | 1.539246 |
-2.001039 | -3.740401 | 303.098407 | 3.098407 | -304.903176 | -4.903176 |
-2.000512 | -1.841796 | 287.368578 | -12.631422 | -304.193416 | -4.193416 |
-2.000588 | -2.116840 | 296.032151 | -3.967849 | -293.461157 | 6.538843 |
-2.002165 | -7.792380 | 293.519429 | -6.480571 | -308.503401 | -8.503401 |
-1.998806 | 4.299325 | 288.416108 | -11.583892 | -312.863066 | -12.863066 |
-1.998250 | 6.298911 | 309.628992 | 9.628992 | -307.434863 | -7.434863 |
-2.000854 | -3.073617 | 302.636745 | 2.636745 | -307.221923 | -7.221923 |
-2.001890 | -6.803619 | 302.187916 | 2.187916 | -312.343793 | -12.343793 |
-1.999119 | 3.171209 | 310.322943 | 10.322943 | -303.069537 | -3.069537 |
-1.999972 | .100164 | 301.962328 | 1.962328 | -272.435744 | 27.564256 |
-2.000345 | 4.744350 | 298.864525 | 7.273446 | -300.409582 | 12.460633 |
(mean value) | (standard deviation) | (mean value) | (standard deviation) | (mean value) | (standard deviation) |
Simulation calculation shows, has solved the fuzzy parameter in the equation, the root-mean-square value of orientation error: 4.74 ".In the simulation calculation direction-finding precision be about 5 seconds illustrative system errors to slacken effect obvious.Wherein differ 1.2 between the mean value that resolves of attitude misalignment and the setting value " and 0.4 ", illustrate that attitude algorithm does not have systematic bias.
Theoretical further emulation explanation: precision and brachium are inversely proportional to, and promptly the long more precision of brachium is high more; Precision just becomes to become with observation time, and promptly the lengthening of observation time helps the raising (being approximately subduplicate relation) of precision.
(3) emulation of different regions
Because the height difference of different regions satellite, mensuration precision to the orientation has certain influence, direction-finding precision for different regions in the observer nation, selection is distributed in China, east, south, west, five areas, north, the systematic error setting is identical with single star direction-finding with observation time, has carried out the emulation of direction-finding.Two satellites (80 °, 140 °) are adopted in this emulation, totally 30 emulation, and its result's statistics sees Table 3.
The precision of table 3 different regions direction-finding
The area | Latitude | Longitude | The satellite number | Simulation times | Bearing accuracy |
The middle part | 30.0° | 100.0° | 2 | 30 | 4.74 |
The east | 30.0° | 120.0° | 2 | 30 | 4.53 |
Western | 30.0° | 80.0° | 2 | 30 | 5.33 |
The south | 20.0° | 100.0° | 2 | 30 | 4.61 |
Northern | 40.0° | 100.0° | 2 | 30 | 4.86 |
On average | 4.81 |
As can be seen, the direction-finding result of different regions has certain difference, its root-mean-square value about 4 "~5 " range changing.
(4) emulation of different observation times
3 minutes 20 seconds emulation of observation time
The area | Latitude | Longitude | The satellite number | Simulation times | Bearing accuracy |
The middle part | 30.0° | 100.0° | 2 | 30 | 17.09 |
The east | 30.0° | 120.0° | 2 | 30 | 18.72 |
Western | 30.0° | 80.0° | 2 | 30 | 14.14 |
The south | 20.0° | 100.0° | 2 | 30 | 17.00 |
Northern | 40.0° | 100.0° | 2 | 30 | 17.20 |
On average | 16.83 |
1 minute 40 seconds emulation of observation time
The area | Latitude | Longitude | The satellite number | Simulation times | Bearing accuracy |
The middle part | 30.0° | 100.0° | 2 | 30 | 28.04 |
The east | 30.0° | 120.0° | 2 | 30 | 29.23 |
Western | 30.0° | 80.0° | 2 | 30 | 25.19 |
The south | 20.0° | 100.0° | 2 | 30 | 27.83 |
Northern | 40.0° | 100.0° | 2 | 30 | 28.28 |
On average | 27.71 |
Above data show, can shorten observation time, and just precision decreases, to adapt to the specific (special) requirements of different user.
8. expanded application
(1) dynamic application
Above-mentioned is orientation and attitude determination in stationary state, and it can be expanded and be used for dynamically, for example the naval vessel.
Different with static application, the orientation of carrier, attitude and customer location are constantly to change in the dynamic application.For the continuous variation of customer location, can in the listed equation of each sampling instant, adopt different navigation position by known survey station coordinate.
When the continuous variation of orientation and attitude, can select one with reference to the moment (for example selecting the last sampling instant constantly) as reference, use inertia device (the light pricker gyro lower) and obtain orientation and the attitude variation constantly of sampling instant relative reference, observed quantity is revised as price.Obtained separating this moment is corresponding to constantly.Carrier transient bearing and attitude obtain with the relative reference variable quantity correction constantly that adds inertia device output with reference to value constantly.
Here only use short time (for example a few minutes) the output difference of inertia device, its zero point drift (main error source of inertia device) can not cause substantial effect to modified value.Greatly reduce the technical requirement of inertial navigation device.Because orientation and attitude are resolved according to Navsat, the error of inertia device can not accumulated.
(2) use other satellite navigation systems
Though above-mentioned at the discussion of earth synchronous and fixed communications satellite navigational system, it is equally applicable to other satellite navigation system.In difference that adopt to adopt other satellite system only is that the centre halfback's championship that establishes an equation is put and directly obtained in not being by the navigational system text, but the user utilizes the navigational parameter calculating satellite position of navigation message.In fact, said method is almost constant this singularity problem of apparent place that has solved geostationary satellite, does not hinder the versatility to other satellite.
Claims (6)
1. satellite goniometer, it is characterized in that, be that be equipped with on support can be by the rotation platform of drive motor driven in rotation, rotation platform comprises the support arm that stretch out at two ends, satellite carrier phase measurement receiver antenna is equipped with at the support arm two ends, telescope is housed on the rotation platform, the level angle electronic logger, vertical angle electronic logger and carrier phase measurement receiver, the carrier phase measurement receiver, the synchronized sampling output of horizontal angle electronic logger and vertical angle electronic logger all is connected with data processing unit, this data processing unit is pressed the many groups phase place behind the cycle slip editor of being gathered, platform corner data and in conjunction with known satellite and survey station coordinate data, carry out data processing according to the mathematical model of drafting that can slacken systematic error, calculate the position angle of scale zero-bit, after this, the any direction of pointing of the telescope all can be read the position angle of this direction from scale, and the wherein said azimuthal mathematical model of scale zero-bit of resolving is:
F is a carrier frequency in the formula, and c is the light velocity, and ρ is the distance of satellite to survey station, D (t
i) be the platform scale corresponding readings that is connected with receiver, D
0For pointing to the circle reading in positive north, dd is the deformation of antenna distance d, a
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch.
2. satellite goniometer as claimed in claim 1 is characterized in that, platform brachium 0.75-2m, platform vertical stability 5 ", platform rotating speed 1 cycle per minute, horizontal angle registration accuracy 2 ", vertical angle registration accuracy 3 ", three quadrature bias 5 ", collimation telescope precision 3 ".
3. satellite goniometer as claimed in claim 1 is characterized in that the port number of two antennas is 4, adopts common local oscillator, phase measurement accuracy≤0.01 week, fundamental frequency degree of stability 10
-8, fundamental frequency accuracy 10
-6, each channel sample postpones≤100 μ s, and two antenna channels sampling delay are poor≤10 μ s.
4. the method for utilizing the described satellite goniometer of claim 1 to carry out direction finding, it is characterized in that, the rotation of driven in rotation platform, make two antennas carry out carrier phase measurement and provide pulse per second (PPS) by 1 time/second, rotation platform carries out the corner record synchronously by pulse per second (PPS) simultaneously, data processing unit is pressed the many groups phase place behind the cycle slip editor of being gathered, platform corner data and in conjunction with known satellite and survey station coordinate data, carry out data processing according to the mathematical model of drafting that can slacken systematic error, calculate the position angle of scale zero-bit, after this, the any direction of pointing of the telescope all can be read the position angle of this direction from scale, and the wherein said azimuthal mathematical model of scale zero-bit of resolving is:
F is a carrier frequency in the formula, and c is the light velocity, and ρ is the distance of satellite to survey station, D (t
i) be the platform scale corresponding readings that is connected with receiver, D
0For pointing to the circle reading in positive north, dd is the deformation of antenna distance d, a
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch.
5. method as claimed in claim 4 is characterized in that, the scale zero-bit is azimuthal resolves plumb line deviation simultaneously resolving, and iteration carries out, and resolves the mathematical model that plumb line deviation separates to be:
F is a carrier frequency in the formula, and c is the light velocity, D (t
i) be the platform scale corresponding readings that is connected with receiver, ξ, η are two quadrature components of plumb line deviation, a
iBe system deviation δ φ
iTo φ
1 s(t
i)-φ
2 s(t
i) partial derivative, N
12(t
1) be complete cycle number scale number poor of the phase observations value that two receiving antennas are constantly obtained first epoch.
6. as claim 4 or 5 described methods, it is characterized in that the removing method of antenna phase center deviation is to make the rotation in revolution of two antennas, rotation platform rotates two all antennas from circling.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101102396A CN100405012C (en) | 2003-12-23 | 2003-12-23 | Satellite direction finder and direction finding method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101102396A CN100405012C (en) | 2003-12-23 | 2003-12-23 | Satellite direction finder and direction finding method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1632465A CN1632465A (en) | 2005-06-29 |
CN100405012C true CN100405012C (en) | 2008-07-23 |
Family
ID=34843165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101102396A Expired - Fee Related CN100405012C (en) | 2003-12-23 | 2003-12-23 | Satellite direction finder and direction finding method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100405012C (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101051872B (en) * | 2007-05-16 | 2011-05-04 | 西安理工大学 | Method for obtaining ASF(TOA), combined ASF(TD) and ASF(TOA) measurement |
CN102062852A (en) * | 2010-11-15 | 2011-05-18 | 西安邮电学院 | Method for realizing rapid orientation by utilizing multiple antennae |
CN102322844B (en) * | 2011-06-13 | 2013-05-08 | 南京中科天文仪器有限公司 | Pointing correction method of telescope based on planar interpolating model |
CN102411137B (en) * | 2011-07-28 | 2015-04-22 | 周必友 | Method for positioning signal source and device thereof |
CN106153021B (en) * | 2016-06-30 | 2019-01-25 | 湖南中森通信科技有限公司 | A kind of north finding method and equipment based on network RTK |
CN108181630B (en) * | 2017-11-30 | 2020-11-03 | 东南大学 | Beidou double-antenna rotation rapid orientation method |
CN108459297A (en) * | 2017-12-29 | 2018-08-28 | 北京日月九天科技有限公司 | A kind of northern instrument of precisely finger that can quickly find direct north |
CN108226960A (en) * | 2017-12-29 | 2018-06-29 | 北京日月九天科技有限公司 | A kind of astronomical telescope that Global Navigation Satellite System direction-measuring function is configured |
-
2003
- 2003-12-23 CN CNB2003101102396A patent/CN100405012C/en not_active Expired - Fee Related
Non-Patent Citations (6)
Title |
---|
利用几何约束快速求解整周模糊度. 逯亮清,胡小平.国防科技大学学报,第23卷第2期. 2001 |
利用几何约束快速求解整周模糊度. 逯亮清,胡小平.国防科技大学学报,第23卷第2期. 2001 * |
双星定姿技术试验研究. 郑伟,吴杰,胡小平,赵健康.空间科学学报,第20卷第3期. 2000 |
双星定姿技术试验研究. 郑伟,吴杰,胡小平,赵健康.空间科学学报,第20卷第3期. 2000 * |
双星定姿的建模与仿真. 吴杰,任萱,胡小平,郑伟,郑冲.空间科学学报,第19卷第2期. 1999 |
双星定姿的建模与仿真. 吴杰,任萱,胡小平,郑伟,郑冲.空间科学学报,第19卷第2期. 1999 * |
Also Published As
Publication number | Publication date |
---|---|
CN1632465A (en) | 2005-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Petrovski et al. | Digital satellite navigation and geophysics: A practical guide with GNSS signal simulator and receiver laboratory | |
CN101246024B (en) | Method for external field fast calibrating miniature multi-sensor combined navigation system | |
CN113311436B (en) | Method for correcting wind measurement of motion attitude of laser wind measuring radar on mobile platform | |
CN110058236A (en) | It is a kind of towards three-dimensional Ground Deformation estimation InSAR and GNSS determine Quan Fangfa | |
CN101975872B (en) | Method for calibrating zero offset of quartz flexible accelerometer component | |
CN103822629B (en) | Positioning system based on multi-directional polarized light navigation sensor and positioning method of positioning system | |
CN101881619B (en) | Ship's inertial navigation and astronomical positioning method based on attitude measurement | |
CN101187568A (en) | Multi-position strapping north-seeking system direction effect calibration method | |
CN101833080A (en) | Method for measuring attitude of carrier by using additional constraint condition of GPS system | |
CN102207386A (en) | North-finding method based on orientation effect error compensation | |
CN106405592B (en) | Vehicle-mounted Beidou carrier phase cycle slips detection and restorative procedure and system | |
CN102426017A (en) | Star-sensor-based method for determining attitude of carrier relative to geographical coordinate system | |
CN110727003A (en) | Pseudo-range simulation method of Beidou satellite navigation system | |
CN100405012C (en) | Satellite direction finder and direction finding method | |
CN103529482A (en) | Method for determining dynamic accelerated speed of carrier precisely | |
CN102998713A (en) | Satellite gravity gradient inversion method based on power spectrum half analysis | |
CN104567802B (en) | The survey line formula land-sea height transfer method of integrated boat-carrying gravity and GNSS | |
CN106959457B (en) | GLONASS almanac parameter estimation method for satellite navigation | |
Wright et al. | The effectiveness of global positioning system electronic navigation | |
CN106767934B (en) | The obliquity sensor installation parameter scaling method of astronomical positioning system horizontal measurement | |
Zhang et al. | Measuring the deflection of the vertical via local reference point surveying and pointing calibration of a VLBI telescope: A case study at the Urumqi station | |
Kumar et al. | The global positioning system: Popular accuracy measures | |
CN104655095B (en) | A kind of method that utilization optical fibre gyro determines geographic latitude | |
JPH07294623A (en) | Position-measuring device and moving body equipped with position-measuring device | |
Karaush et al. | Calibration of GNSS receiver with GLONASS signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080723 Termination date: 20131223 |