CN102519454A - Selenocentric direction correction method for sun-earth-moon navigation - Google Patents

Abstract

A moon center direction correction method for sun-earth-moon navigation. Firstly, the moon center is directly calculated by using the fan-shaped field of view of the integrated sun-earth-moon navigation sensor to scan the moon's pulse time without considering the periodic change of the moon's bright area seen by the satellite. Then, considering the bright area of the moon actually visible to the satellite, calculate the measurement compensation amount corresponding to the change of the bright area according to the measurement principle, correct the pulse time when the fan-shaped field of view of the integrated sun-earth-moon navigation sensor scans the moon, and recalculate the corrected and The direction of the center of the moon is compensated, and the invention greatly improves the navigation accuracy of the autonomous navigation system based on the measurement information of the sun, the earth and the moon.

Description

The moon heart adjustment in direction method of ground month navigation in a kind of day

Technical field

The invention belongs to spacecraft independent navigation field, relate to the moon heart adjustment in direction method of ground month navigation in a kind of day.

Background technology

Navigation information is carried out aerial mission for satellite and is had great importance.Up to the present, satellite navigation both domestic and external mainly relies on uphole equipment to accomplish, and this navigate mode exists many shortcomings such as expensive, the autonomous viability difference of operation.Just because of recognizing with the ground observing and controlling is the defective that main traditional mode exists in operation and management, the autonomous navigation of satellite Study on Technology has been caused the extensive concern of Chinese scholars.The U.S., Russia, European Space Agency etc. just begin studying autonomous navigation of satellite technology from the eighties of last century sixties; Successively having proposed a plurality of is the autonomous navigation system scheme of main background with the applied satellite; And developed corresponding spaceborne surveying instrument, comprise Horizon scanner, space sextant or the like.China is also carrying out the research and the exploration of spacecraft autonomous navigation technology always, but starts late, and most studies also only rests on the demonstration of method stage.

From considerations such as quality, cost, power consumption and batch processes, miniaturization, integrated, the low-cost development trend that has become following autonomous navigation sensor.The MANS autonomous navigation system of U.S. Microcosm company development, utilize special-purpose integrated sensor according to day,, month confirm the track and the three-axis attitude of spacecraft in real time in the rail measurement data, be the autonomous navigation system on the complete meaning.Because the sun, the earth, the moon are the nature celestial bodies, the clear and identification easily of its characteristics of motion, the navigation resource is unrestricted.Therefore, study this autonomous navigation technology and have important engineering application value." space test platform-D " spacecraft that USAF is launched as far back as in March, 1994 just carries this system to have carried out in the rail flight test.

Traditional independent navigation algorithm based on day ground month measurement exists the low problem of navigation accuracy, influences and restricted the development and the application of the autonomous navigation system of measuring based on ground moon day greatly.Along with the development of aerospace industry, increasingly high requirement has been proposed for the precision of satellite navigation.For example be Navsat, geodetic satellite (GEOS) and the reconnaissance satellite etc. of the service of national economy and military technology; All need survey the rail system its high precision track and space tracking at space motion is provided; Thereby for the user accomplishes the precision navigation location; Measure the accurate position mark of ground point position, perhaps effectively monitoring, tracking target are also confirmed exact position of target or the like.In a word, realize high-precision autonomous navigation of satellite, the navigation accuracy that how to improve the autonomous navigation system of measuring based on ground moon day more becomes the challenge problem in international aerospace navigation, guidance and control field.

Summary of the invention

Technology of the present invention is dealt with problems and is: the deficiency that overcomes prior art; The moon heart adjustment in direction method of ground month navigation in a kind of day is provided; The integrated day ground month fan-shaped visual field scanning of navigation sensor revised to the pulse of the moon constantly; Thereby improved the precision of definite month heart direction, finally improved the precision of autonomous navigation system.

Technical solution of the present invention is: the moon heart adjustment in direction method performing step of ground month navigation in a kind of day is following:

(1) calculates the moon heart direction position angle φ in integrated day ground month navigation sensor measurement coordinate system revise _MWith angle of pitch δ _M, formula is following

${\mathrm{\φ}}_M={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}+t_{L2-\mathrm{MOON}}}{2}-t_{M1-M-\mathrm{REF}})---\left(22\right)$

δ _M＝arctg(sinμ _Mctgβ _L) (23)

In the formula, ω _ROTThe sweep speed of integrated day ground month navigation sensor; t _L1-MOONMonth fan-shaped visual field 1 of navigation sensor, integrated day ground scans the pulse moment of the moon; t _L2-MOONMonth fan-shaped visual field 2 of navigation sensor, integrated day ground scans the pulse moment of the moon; t _M1-M2-REFThe integrated day ground month navigation sensor plane of symmetry is through the pulse moment of reference point on integrated day ground month navigation sensor equatorial plane;

D _L1-L2Month fan-shaped visual field 1 of navigation sensor, integrated day ground was ahead of the angle of month fan-shaped visual field 2 of navigation sensor, integrated day ground on the ground month navigation sensor equatorial plane on integrated; β _LFan-shaped visual field 1 of month navigation sensor, integrated day ground and fan-shaped visual field 2 are with respect to the angle of integrated day ground month navigation sensor scanning axes of rotation skew.

(2) the coordinate E of moon heart direction in integrated day ground month navigation sensor measurement coordinate system before calculate revising _Sm, formula is following

$E_{\mathrm{sm}}=\left[\begin{array}{c}{\sin \mathrm{\δ}}_M\\ \cos {\mathrm{\δ}}_M\sin {\mathrm{\φ}}_M\\ \cos {\mathrm{\δ}}_M\cos {\mathrm{\φ}}_M\end{array}\right]---\left(24\right)$

(3) calculate to revise vector that satellite before month heart direction points to the moon three component x in integrated day ground month navigation sensor measurement coordinate system _m, y _m, z _m, formula is following

$\left[\begin{array}{c}{x}_m\\ y_m\\ z_m\end{array}\right]=L_{\mathrm{ms}}{E}_{\mathrm{sm}}---\left(25\right)$

In the formula, L _MsFor the moon arrives satellite distance.

(4) calculate to revise before month heart direction angle μ between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ＝arccos(z _m/L _ms) (26)

(5) calculate poor between maximum visual angle, moon clear zone and the actual angle of visibility.

γ＝arccos(R _m/L _ms)-li (27)

In the formula, R _mBe the moon radius of a ball; Li is actual angle of visibility, is to obtain according to a day ground month ephemeris computation of storing on satellite orbit that forecasts on the star and the star.

(6) calculate offset angle λ, formula is following

λ＝arccos[(L _ms ²-R _m ²+R _m ²cosγ)/L _ms ²] (28)

(7) calculate intermediate variable η, formula is following

η＝arccos{[L _ms ²-R _m ²+R _m ²cos(li)]/L _ms ²}?(29)

(8) calculate to revise back month heart direction angle μ ' between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ′＝μ+[arccos(R _m/L _ms)-η-λ] (30)

(9) calculate satellite and point to the projection r of moon vector at the integrated day ground month navigation sensor equatorial plane _k, formula is following

$r_k=\sqrt{{L_{\mathrm{ms}}}^2-{x_m}^2}---\left(31\right)$

(10) calculate intermediate variable Q, formula is following

$Q=\sqrt{{x_m}^2+{(L_{\mathrm{ms}}-r_k)}^2}---\left(32\right)$

(11) calculate intermediate variable T, formula is following

$T=\frac{\sqrt{Q^2{r_k}^2{\left(\cos \left({\mathrm{\β}}_L\right)\right)}^2-(L_{\mathrm{ms}}-r_k)({L_{\mathrm{ms}}}^2{r}_k-L_{\mathrm{ms}}{x_m}^2-{L_{\mathrm{ms}}}^3-r_k{Q}^2)}-{\mathrm{Qr}}_k\cos \left({\mathrm{\β}}_L\right)}{L_{\mathrm{ms}}-r_k}---\left(33\right)$

(12) calculate scanning corner ε, formula is following

$\mathrm{\ϵ}=\arccos \left(\frac{{L_{\mathrm{ms}}}^2+{r_k}^2+{x_m}^2-T^2}{{2L}_{\mathrm{ms}}{r}_k}\right)---\left(34\right)$

(13) calculate the angle ε _ L1 and the ε _ L2 of the projection on integrated day ground month navigation sensor equatorial plane of the fan-shaped visual field of integrated day ground month navigation sensor 1, fan-shaped visual field 2 and month heart direction, be divided into following 4 kinds of situation

1. if y _m＜0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(35\right)$

2. if y _m＜0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\ϵ}-\mathrm{\μ})\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(36\right)$

3. if y _m>=0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=2\mathrm{\π}-\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(37\right)$

4. if y _m>=0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\μ}+\mathrm{\ϵ})\\ \mathrm{\ϵ}-L2=\mathrm{\ϵ}-\mathrm{\μ}\end{array}\right.---\left(38\right)$

(14) the integrated day ground moon fan-shaped visual field 1 of navigation sensor and 2 after the calculating correction month heart direction scans the pulse moment t ' of the moon _L1-MOONAnd t ' _L2-MOON, formula is following

$\left\{\begin{array}{c}{t}_{L1-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L1/{\mathrm{\ω}}_{\mathrm{ROT}}\\ t_{L2-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L2/{\mathrm{\ω}}_{\mathrm{ROT}}\end{array}\right.---\left(39\right)$

(15) calculate revised month heart direction position angle φ ' in integrated day ground month navigation sensor measurement coordinate system _MWith angle of pitch δ ' _M, formula is following

${\mathrm{\φ}}_M^{\′}={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}^{\′}+t_{L2-\mathrm{MOON}}^{\′}}{2}-t_{M1-M2-\mathrm{REF}})---\left(40\right)$

δ′ _M＝arctg(sinμ′ _Mctgβ _L) (41)

In the formula, ${\mathrm{\μ}}_M^{\′}=\frac{1}{2}{\mathrm{\ω}}_{\mathrm{ROT}}(t_{L2-\mathrm{MOON}}^{\′}-t_{L1-\mathrm{MOON}}^{\′})-\frac{D_{L1-L2}}{2}.$

(16) calculate revised month heart direction coordinate E ' in integrated day ground month navigation sensor measurement coordinate system _Sm, formula is following

$E_{\mathrm{sm}}^{\′}=\left[\begin{array}{c}\sin {\mathrm{\δ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\sin {\mathrm{\φ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\cos {\mathrm{\φ}}_M^{\′}\end{array}\right]---\left(42\right)$

Thereby accomplish moon heart adjustment in direction to day ground month navigation.

The present invention's advantage compared with prior art: the shape that prior art is ignored the moon clear zone that can observe on the star is the actual conditions that constantly change, and the moon heart direction that obtains through integrated day ground month navigation sensor scanning survey information has certain angular deviation with an actual moon heart direction.The present invention introduces a month heart direction with the visible moon clear zone of the satellite of predicting on star information and confirms algorithm; The integrated day ground month fan-shaped visual field scanning of navigation sensor revised to the pulse of the moon constantly; Thereby improved the precision of definite month heart direction, finally improved the precision of autonomous navigation system.

Description of drawings

The synoptic diagram of the visual field of integrated day ground month navigation sensor of Fig. 1;

The synoptic diagram of the measurement coordinate system of integrated day ground month navigation sensor of Fig. 2;

Fig. 3 is the synoptic diagram to the geometric relationship of measuring by the moon;

Fig. 4 is a month heart directional compensation synoptic diagram;

Fig. 5 is preceding month heart position angle of compensation and elevation angle error;

Fig. 6 is for adopting position angle and the elevation angle error after the inventive method compensates.

Embodiment

A ground month navigation sensor comprised two double cone scanning type earth sensors and two sector display formula life sensors in integrated day.Integrated day ground month navigation sensor can be responsive to earth heat radiation, again can be responsive to day, month visible light.

In Fig. 1, the visual field that infrared visual field 1 and infrared visual field 2 are double cone scanning type earth sensors, the fan-shaped visual field of visible light is the visual field of sector display formula life sensor.

What the present invention mainly adopted is the integrated day sector display formula life sensor in the ground month sensor; They utilize the si-photodiode detector can be responsive to the sun and the moon, and the moment that in fan-shaped visual field, occurs according to the sun and the moon can calculate the orientation of its direction vector with respect to satellite.The visual field of each visible light sensor is that 72 degree are long, the wide fan-shaped shape of slit of 2.5 degree, and it is with respect to scanning axes of rotation skew 16 degree angles.In the optical head scanning process, fan-shaped visual field scanning cross the space for and scanning rotating shaft angle spend to the spherical zone zone between 87 degree 23.

As shown in Figure 2, the measurement coordinate system of integrated day ground month navigation sensor.Initial point O _sThe intersection point of the scanning rotating shaft of integrated day ground month navigation sensor and integrated day ground month sensor equatorial plane, X _sAlong scanning rotor shaft direction, Z _sIn the sensor equatorial plane, and make the reference point that is fixed on sensor be positioned at O _s-Z _sX _sIn the plane, Y _sMake O _s-X _sY _sZ _sConstitute right hand orthonormal system.

The fan-shaped visual field of two sector display formula life sensors shown in Fig. 3 is designated as fan-shaped visual field 1 and fan-shaped visual field 2 respectively, the symmetrical plane of fan-shaped visual field 1 and 2 is designated as the M1-M2 face.Fan-shaped visual field 1 was ahead of fan-shaped visual field 2 on the ground month navigation sensor equatorial plane on integrated angle is designated as D _L1-L2, see Fig. 3.Fan-shaped visual field 1 and 2 angles with respect to the scanning axes of rotation skew are designated as β _L, see Fig. 3.Here D _L1-L2=4 degree, β _L=16 degree.

Integrated day ground month navigation sensor mainly contains three times to the basic observed quantity of the moon: the pulse that the moon is swept in fan-shaped visual field 1 is t constantly _L1-MOON, the fan-shaped visual field 2 pulse t constantly that sweeps to the moon _L2-MOON, the pulse of plane of symmetry M1-M2 through reference point t constantly _M1-M2-REFCan obtain position angle and the elevation angle of moon heart direction according to these three observed quantities with respect to integrated day ground month navigation sensor measurement coordinate system, and the coordinate of the satellite vector that points to the moon in integrated day ground month navigation sensor measurement coordinate system.

The sweep speed of integrated day ground month navigation sensor is designated as ω _ROT, ω _ROT=240 rev/mins=8 π radian per seconds, getting sensor can only be rotated counterclockwise.

In Fig. 4, provided the principle of the correction moon of the present invention heart direction.Wherein, solid line representes that partly dotted portion is represented the edge, clear zone under the situation that the moon of simulating after the over-compensation be not blocked from the edge, the actual clear zone of the moon that satellite is seen.

Known integrated day ground month navigation sensor scanning survey information t _L1-MOON, t _L2-MOON, t _M1-M2-REF, then the present invention calculates the moon heart direction E ' that contains compensation rate _SmMethod step following:

(1) calculates the moon heart direction position angle φ in integrated day ground month navigation sensor measurement coordinate system revise _MWith angle of pitch δ _M, formula is following

${\mathrm{\φ}}_M={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}+t_{L2-\mathrm{MOON}}}{2}-t_{M1-M2-\mathrm{REF}})---\left(43\right)$

δ _M＝arctg(sinμ _Mctgβ _L) (44)

In the formula, ${\mathrm{\μ}}_M=\frac{1}{2}{\mathrm{\ω}}_{\mathrm{ROT}}(t_{L2-\mathrm{MOON}}-t_{L1-\mathrm{MOON}})-\frac{D_{L1-L2}}{2}.$

(2) the coordinate E of moon heart direction in integrated day ground month navigation sensor measurement coordinate system before calculate revising _Sm, formula is following

$E_{\mathrm{sm}}=\left[\begin{array}{c}{\sin \mathrm{\δ}}_M\\ \cos {\mathrm{\δ}}_M\sin {\mathrm{\φ}}_M\\ \cos {\mathrm{\δ}}_M\cos {\mathrm{\φ}}_M\end{array}\right]---\left(45\right)$

(3) calculate to revise vector that satellite before month heart direction points to the moon three component x in integrated day ground month navigation sensor measurement coordinate system _m, y _m, z _m, formula is following

$\left[\begin{array}{c}{x}_m\\ y_m\\ z_m\end{array}\right]=L_{\mathrm{ms}}{E}_{\mathrm{sm}}---\left(46\right)$

In the formula, L _MsFor the moon arrives satellite distance.

(4) calculate to revise before month heart direction angle μ between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ＝arccos(z _m/L _ms) (47)

(5) calculate poor between maximum visual angle, moon clear zone and the actual angle of visibility.

γ＝arccos(R _m/L _ms)-li (48)

In the formula, R _mBe the moon radius of a ball; Li is actual angle of visibility, is to obtain according to a day ground month ephemeris computation of storing on satellite orbit that forecasts on the star and the star.

(6) calculate offset angle λ, formula is following

λ＝arccos[(L _ms ²-R _m ²+R _m ²cosγ)/L _ms ²] (49)

(7) calculate intermediate variable η, formula is following

η＝arccos{[L _ms ²-R _m ²+R _m ²cos(li)]/L _ms ²} (50)

(8) calculate to revise back month heart direction angle μ ' between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ′＝μ+[arccos(R _m/L _ms)-η-λ] (51)

(9) calculate satellite and point to the projection r of moon vector at the integrated day ground month navigation sensor equatorial plane _k, formula is following

$r_k=\sqrt{{L_{\mathrm{ms}}}^2-{x_m}^2}---\left(52\right)$

(10) calculate intermediate variable Q, formula is following

$Q=\sqrt{{x_m}^2+{(L_{\mathrm{ms}}-r_k)}^2}---\left(53\right)$

(11) calculate intermediate variable T, formula is following

$T=\frac{\sqrt{Q^2{r_k}^2{\left(\cos \left({\mathrm{\β}}_L\right)\right)}^2-(L_{\mathrm{ms}}-r_k)({L_{\mathrm{ms}}}^2{r}_k-L_{\mathrm{ms}}{x_m}^2-{L_{\mathrm{ms}}}^3-r_k{Q}^2)}-{\mathrm{Qr}}_k\cos \left({\mathrm{\β}}_L\right)}{L_{\mathrm{ms}}-r_k}---\left(54\right)$

(12) calculate scanning corner ε, formula is following

$\mathrm{\ϵ}=\arccos \left(\frac{{L_{\mathrm{ms}}}^2+{r_k}^2+{x_m}^2-T^2}{{2L}_{\mathrm{ms}}{r}_k}\right)---\left(55\right)$

(13) calculate the angle ε _ L1 and the ε _ L2 of the projection on integrated day ground month navigation sensor equatorial plane of the fan-shaped visual field of integrated day ground month navigation sensor 1, fan-shaped visual field 2 and month heart direction, be divided into following 4 kinds of situation

1. if y _m＜0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(56\right)$

2. if y _m＜0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\ϵ}-\mathrm{\μ})\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(57\right)$

3. if y _m>=0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=2\mathrm{\π}-\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(58\right)$

4. if y _m>=0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\μ}+\mathrm{\ϵ})\\ \mathrm{\ϵ}-L2=\mathrm{\ϵ}-\mathrm{\μ}\end{array}\right.---\left(59\right)$

(14) the integrated day ground moon fan-shaped visual field 1 of navigation sensor and 2 after the calculating correction month heart direction scans the pulse moment t ' of the moon _L1-MOONAnd t ' _L2-MOON, formula is following

$\left\{\begin{array}{c}{t}_{L1-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L1/{\mathrm{\ω}}_{\mathrm{ROT}}\\ t_{L2-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L2/{\mathrm{\ω}}_{\mathrm{ROT}}\end{array}\right.---\left(60\right)$

(15) calculate revised month heart direction position angle φ ' in integrated day ground month navigation sensor measurement coordinate system _MWith angle of pitch δ ' _M, formula is following

${\mathrm{\φ}}_M^{\′}={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}^{\′}+t_{L2-\mathrm{MOON}}^{\′}}{2}-t_{M1-M2-\mathrm{REF}})---\left(61\right)$

δ′ _M＝arctg(sinμ′ _Mctgβ _L) (62)

In the formula, ${\mathrm{\μ}}_M^{\′}=\frac{1}{2}{\mathrm{\ω}}_{\mathrm{ROT}}(t_{L2-\mathrm{MOON}}^{\′}-t_{L1-\mathrm{MOON}}^{\′})-\frac{D_{L1-L2}}{2}.$

(16) calculate revised month heart direction coordinate E ' in integrated day ground month navigation sensor measurement coordinate system _Sm, formula is following

$E_{\mathrm{sm}}^{\′}=\left[\begin{array}{c}\sin {\mathrm{\δ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\sin {\mathrm{\φ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\cos {\mathrm{\φ}}_M^{\′}\end{array}\right]---\left(63\right)$

So far, calculated revised month heart direction E ' _Sm, compare the moon heart direction E before revising _Sm, its precision is higher.

Simulation example: the flying height of satellite is the near-circular orbit of 500km, and three attitude angle are zero, and the nominal orbit parameter is a=6878.14km, e=0.001, i=97.4 °, Ω=157.5 °, ω=90 °, M=0 °.

Attitude misalignment, the random measurement noise mean square deviation of the moon is taken as 0.05 ° (3 σ).Track condition is estimated initial value: there is the 10km error in semi-major axis a, the same nominal value of other state initial value.

The emulation initial time is got 0 o'clock on the 1st January in 2010, simulation time 3600s, step-length 60s.Simulation result shown in Fig. 5～7, wherein position angle and height angular unit degree of being.

Fig. 5 has provided the error of position angle and the elevation angle in integrated day ground month navigation sensor measurement coordinate system of month heart direction vector before the compensation; Fig. 6 has provided position angle and the elevation angle error after the compensation.

Can find out by Fig. 5 and Fig. 6; Compensation forward angle error is about 0.06 °～0.4 °; The elevation angle error is between-0.08 °～0.08 °, and compensation back bearing error is about-0.09 °～0.13 °, and the elevation angle error is between-0.045 °～0.045 °.It is thus clear that the moon clear zone shape through satellite is seen compensates, and can improve definite precision of moon heart direction greatly, especially azimuthal definite precision is more than doubled.

The content of not doing to describe in detail in the instructions of the present invention belongs to those skilled in the art's known technology.

Claims (1)

1. the moon heart adjustment in direction method of navigating by one kind of day ground moon is characterized in that step is following:

(1) calculates the moon heart direction position angle φ in integrated day ground month navigation sensor measurement coordinate system revise _MWith angle of pitch δ _M, formula is following

${\mathrm{\φ}}_M={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}+t_{L2-\mathrm{MOON}}}{2}-t_{M1-M-\mathrm{REF}})---\left(1\right)$

δ _M＝arctg(sinμ _Mctgβ _L) (2)

In the formula, ω _ROTThe sweep speed of integrated day ground month navigation sensor; t _L1-MOONMonth fan-shaped visual field 1 of navigation sensor, integrated day ground scans the pulse moment of the moon; t _L2-MOONMonth fan-shaped visual field 2 of navigation sensor, integrated day ground scans the pulse moment of the moon; t _M1-M2-REFThe integrated day ground month navigation sensor plane of symmetry is through the pulse moment of reference point on integrated day ground month navigation sensor equatorial plane;

D _L1-L2Month fan-shaped visual field 1 of navigation sensor, integrated day ground was ahead of the angle of month fan-shaped visual field 2 of navigation sensor, integrated day ground on the ground month navigation sensor equatorial plane on integrated; β _LFan-shaped visual field 1 of month navigation sensor, integrated day ground and fan-shaped visual field 2 are with respect to the angle of integrated day ground month navigation sensor scanning axes of rotation skew;

(2) the coordinate E of moon heart direction in integrated day ground month navigation sensor measurement coordinate system before calculate revising _Sm, formula is following

$E_{\mathrm{sm}}=\left[\begin{array}{c}{\sin \mathrm{\δ}}_M\\ \cos {\mathrm{\δ}}_M\sin {\mathrm{\φ}}_M\\ \cos {\mathrm{\δ}}_M\cos {\mathrm{\φ}}_M\end{array}\right]---\left(3\right)$

(3) calculate to revise vector that satellite before month heart direction points to the moon three component x in integrated day ground month navigation sensor measurement coordinate system _m, y _m, z _m, formula is following

$\left[\begin{array}{c}{x}_m\\ y_m\\ z_m\end{array}\right]=L_{\mathrm{ms}}{E}_{\mathrm{sm}}---\left(4\right)$

In the formula, L _MsFor the moon arrives satellite distance;

(4) calculate to revise before month heart direction angle μ between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ＝arccos(z _m/L _ms) (5)

(5) calculate poor between maximum visual angle, moon clear zone and the actual angle of visibility.

γ＝arccos(R _m/L _ms)-li (6)

In the formula, R _mBe the moon radius of a ball; Li is actual angle of visibility, is to obtain according to a day ground month ephemeris computation of storing on satellite orbit that forecasts on the star and the star;

(6) calculate offset angle λ, formula is following

λ＝arccos[(L _ms ²-R _m ²+R _m ²cosγ)/L _ms ²] (7)

(7) calculate intermediate variable η, formula is following

η＝arccos{[L _ms ²-R _m ²+R _m ²cos(li)]/L _ms ²}(8)

(8) calculate to revise back month heart direction angle μ ' between the reference point on the projection of the integrated day ground month navigation sensor equatorial plane and the integrated day ground moon navigation sensor equatorial plane, formula is following

μ′＝μ+[arccos(R _m/L _ms)-η-λ] (9)

(9) calculate satellite and point to the projection r of moon vector at the integrated day ground month navigation sensor equatorial plane _k, formula is following

$r_k=\sqrt{{L_{\mathrm{ms}}}^2-{x_m}^2}---\left(10\right)$

(10) calculate intermediate variable Q, formula is following

$Q=\sqrt{{x_m}^2+{(L_{\mathrm{ms}}-r_k)}^2}---\left(11\right)$

(11) calculate intermediate variable T, formula is following

$T=\frac{\sqrt{Q^2{r_k}^2{\left(\cos \left({\mathrm{\β}}_L\right)\right)}^2-(L_{\mathrm{ms}}-r_k)({L_{\mathrm{ms}}}^2{r}_k-L_{\mathrm{ms}}{x_m}^2-{L_{\mathrm{ms}}}^3-r_k{Q}^2)}-{\mathrm{Qr}}_k\cos \left({\mathrm{\β}}_L\right)}{L_{\mathrm{ms}}-r_k}---\left(12\right)$

(12) calculate scanning corner ε, formula is following

$\mathrm{\ϵ}=\arccos \left(\frac{{L_{\mathrm{ms}}}^2+{r_k}^2+{x_m}^2-T^2}{{2L}_{\mathrm{ms}}{r}_k}\right)---\left(13\right)$

(13) calculate the angle ε _ L1 and the ε _ L2 of the projection on integrated day ground month navigation sensor equatorial plane of the fan-shaped visual field of integrated day ground month navigation sensor 1, fan-shaped visual field 2 and month heart direction, be divided into following 4 kinds of situation

1. if y _m＜0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(14\right)$

2. if y _m＜0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\ϵ}-\mathrm{\μ})\\ \mathrm{\ϵ}-L2=\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(15\right)$

3. if y _m>=0 and μ>=ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-\mathrm{\μ}-\mathrm{\ϵ}\\ \mathrm{\ϵ}-L2=2\mathrm{\π}-\mathrm{\μ}+\mathrm{\ϵ}\end{array}\right.---\left(16\right)$

4. if y _m>=0 and μ＜ε, then

$\left\{\begin{array}{c}\mathrm{\ϵ}\_L1=2\mathrm{\π}-(\mathrm{\μ}+\mathrm{\ϵ})\\ \mathrm{\ϵ}-L2=\mathrm{\ϵ}-\mathrm{\μ}\end{array}\right.---\left(17\right)$

(14) calculate to revise pulse that integrated day fan-shaped visual field 1 of ground month navigation sensor and fan-shaped visual field 2 after month heart direction scan moon t ' constantly _L1-MOONAnd t ' _L2-MOON, formula is following

$\left\{\begin{array}{c}{t}_{L1-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L1/{\mathrm{\ω}}_{\mathrm{ROT}}\\ t_{L2-\mathrm{MOON}}^{\′}=\mathrm{\ϵ}\_L2/{\mathrm{\ω}}_{\mathrm{ROT}}\end{array}\right.---\left(18\right)$

(15) calculate revised month heart direction position angle φ ' in integrated day ground month navigation sensor measurement coordinate system _MWith angle of pitch δ ' _M, formula is following

${\mathrm{\φ}}_M^{\′}={\mathrm{\ω}}_{\mathrm{ROT}}(\frac{t_{L1-\mathrm{MOON}}^{\′}+t_{L2-\mathrm{MOON}}^{\′}}{2}-t_{M1-M2-\mathrm{REF}})---\left(19\right)$

δ′ _M＝arctg(sinμ′ _Mctgβ _L) (20)

In the formula, ${\mathrm{\μ}}_M^{\′}=\frac{1}{2}{\mathrm{\ω}}_{\mathrm{ROT}}(t_{L2-\mathrm{MOON}}^{\′}-t_{L1-\mathrm{MOON}}^{\′})-\frac{D_{L1-L2}}{2}.$

(16) calculate revised month heart direction coordinate E ' in integrated day ground month navigation sensor measurement coordinate system _Sm, formula is following

$E_{\mathrm{sm}}^{\′}=\left[\begin{array}{c}\sin {\mathrm{\δ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\sin {\mathrm{\φ}}_M^{\′}\\ \cos {\mathrm{\δ}}_M^{\′}\cos {\mathrm{\φ}}_M^{\′}\end{array}\right]---\left(21\right).$

Images