CN101923169A - Atmospheric model for high-precision orbit confirmation of phased array radar - Google Patents
Atmospheric model for high-precision orbit confirmation of phased array radar Download PDFInfo
- Publication number
- CN101923169A CN101923169A CN2010102378952A CN201010237895A CN101923169A CN 101923169 A CN101923169 A CN 101923169A CN 2010102378952 A CN2010102378952 A CN 2010102378952A CN 201010237895 A CN201010237895 A CN 201010237895A CN 101923169 A CN101923169 A CN 101923169A
- Authority
- CN
- China
- Prior art keywords
- theta
- partiald
- rho
- cos
- prime
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to an atmospheric model for high-precision orbit confirmation of a phased array radar, which aims at providing an atmospheric model applicable to the high-precision orbit confirmation of a phased array radar. Mathematic singularities in DTM 94 atmospheric model density calculation are removed. The invention has the technical scheme that the expression form of the atmospheric density Rho in the atmospheric model is shown in the specifications, wherein subscript i equals to 1, 2......5 which respectively correspond to five gas atoms or molecules of H, He, O, N2 and O2.
Description
Technical field
The present invention relates to the Atmospheric models field, be specifically related to the Atmospheric models in the phased-array radar high-precision orbit.
Background technology
In Satellite Orbit Determination research, near-earth satellite particularly, atmospherical drag are a kind of perturbation sources that must consider.We relate to atmospheric density ρ and the atmospheric density partial derivative with respect to satellite position to the consideration of atmospherical drag in the phased-array radar high-precision orbit
Calculating.Atmospheric density ρ is provided by Atmospheric models, and therefore the research of Atmospheric models receive much concern.The people such as C.Berger of France are in the article " Improvement of theempirical thermospheric model DTM:DTM94-a comparative review of various temporalvariations and prospects in space geodesy app1ications " of Journalof Geodesy (1998) 72:161-178, the Atmospheric models of the still comparatively general level that is in a leading position, the i.e. hot Atmospheric models of DTM94 have in the world been proposed now.Domestic Jiang Hu, yellow alkali also has article " DTM94 Atmospheric models and application thereof " that these Atmospheric models are introduced in Yunnan Observatory platform 1999 the 4th phases of periodical, and its effect is affirmed.The hot Atmospheric models of DTM94 have characteristics such as precision height, computing velocity be fast, easy to use.
Yet we find in the process of exploitation phased-array radar high-precision orbit software, and in the use, the hot Atmospheric models of DTM94 are owing to adopt spherical coordinates as basic variable, so its atmospheric density is calculated ball the two poles of the earth and had incalculable singular point.Simultaneously, if directly the partial derivative of DTM94 atmospheric density with respect to the solid square coordinate in ground calculated in differentiate, then can on denominator, occur
The factor, the molecule of overhead partial derivative and denominator can appear as zero situation simultaneously at the two poles of the earth, cause near the partial derivative in the sky, the two poles of the earth to calculate and may collapse.Because the existence of this singular point causes correlation computations to lose efficacy near the two poles of the earth, has had a strong impact on the versatility and the robustness of software.Therefore we need seek a kind of improving one's methods, and can make phased-array radar high-precision orbit software be unlikely to singular point problem partial failure because of Atmospheric models.
Consider that from physical angle it is continuous that atmospheric density distributes, the singular point that do not have exists.And the hot Atmospheric models of DTM94 are when calculating atmospheric density and atmospheric density with respect to the partial derivative of the solid square coordinate in ground and since adopt be spherical coordinates as basic variable, so, have incalculable singular point in the terrestrial pole sky.This is mainly reflected in following several:
Sunday item:
D={a
21P
11+a
22P
31+a
23P
51+(a
24P
11+a
25P
21)cos[Ω(d-a
18)]}cosωt
+{a
26P
11+a
27P
31+a
28P
51+(a
29P
11+a
30P
21)cos[Ω(d-a
18)]}sinωt
Half cycle day item:
SD={a
31P
22+a
32P
32cos[Ω(d-a
18)]}cos2ωt
+{a
33P
22+a
34P
32cos[Ω(d-a
18)]}sin2ωt
1/3 Sunday item:
TD=a
35P
33cos3ωt+a
36P
33sin3ωt;
The geomagnetic activity item:
λ in the top formula,
Be the sub-satellite point longitude and latitude, t is the local true solar time of sub-satellite point, a
NmBe normal value coefficient, can from the hot Atmospheric models coefficient table of DTM94, look into and get, ω=2 pi/2s, 4 (hour
-1), Ω=2 π/365 (day
-1) also be constant, P
NmBe about
Association rein in allow polynomial expression.As seen, in the sky, the two poles of the earth, the local true solar time t of substar longitude λ and substar all loses definition, becomes above several mathematics singular point, causes the DTM94 mathematical model can't calculate at the two poles of the earth.But it is not an essential singularity, and the appropriate mathematical method just should be able to be overcome and solve by reference.
Summary of the invention
The object of the present invention is to provide a kind of Atmospheric models that are applicable in the phased-array radar high-precision orbit, eliminate the mathematics singular point in the density calculation of DTM94 Atmospheric models.The present invention is based on the hot Atmospheric models of DTM94,, the atmospheric density of any can be calculated, realize that the no singular point of atmospheric density calculates improving of atmospheric density calculating method in the DTM94 Atmospheric models.Further goal of the invention of the present invention is to provide the no singular point computing method of atmospheric density with respect to ground solid square coordinate partial derivative.
The present invention uses for reference (Pines, 1973) and eliminates the used method of the terrestrial gravitation humorous expression formula singular point of gesture ball, with satellite the earth's core square r and direction cosine θ
1, θ
2, θ
3The spherical coordinates r that replaces it, λ,
As the basic variable of the hot Atmospheric models of DTM94, make λ and t no longer to occur in the formula.Make square coordinate and direction cosine that definition all be arranged in the arbitrfary point, therefore calculate meaning point in office and can realize, do not have mathematics singular point problem through the atmospheric density of transforming.Make that simultaneously atmospheric density can not be zero the collapse because of the molecule denominator with respect to the calculating of ground solid square coordinate partial derivative simultaneously.
The technical scheme that realizes the object of the invention is:
Atmospheric models in the phased-array radar high-precision orbit, the expression-form of atmospheric density ρ is:
Subscript i=1,2 ... 5, the corresponding H of difference, H
e, O, N
2, O
2This five kinds of gas atoms or molecule.
In the following formula,
A
VR=1/ (6.022 * 10
23) be the inverse of avogadros constant;
M
iBe the molecular weight of neutral gas composition i, M
1=1, M
2=4, M
3=16, M
4=28, M
5=32;
n
i(h)=A
iExp (G
i(L)) f
i(h), h is the floor level of satellite;
Following formula the right symbol description is as follows:
This wherein,
T
120=380K represents 120 kilometers highly temperature at place,
T
∞=A
T(1+G
T(L)), the expression Outer Atmospheric Temperature,
α
iBe the thermal diffusion coefficient of neutral gas composition i, α
1=-0.40, α
2=-0.38, α
3=α
4=α
5=0.0,
γ
i=M
ig
120/(σRT
∞),
R=831.4 is a gas law constant,
g
120=944.6626 (cm/s
2), the gravitational acceleration at 120 kilometers places of expression floor level,
σ=T '
120/ (T
∞-T
120), represent relative vertical temperature gradient,
T '
120=14.348, the thermograde at 120 kilometers places of expression,
R
ρ=6356.77km, expression earth polar radius,
e
2=2 ε-ε
2, expression earth meridian circle excentricity,
ε=1/298.25, expression earth meridian circle ellipticity.
And other coefficients:
A
5=4.775×10
10(cm
-3),
A
1, A
2, A
3, A
4And A
TValue is looked in can the coefficient table by DTM94 model (Berger et al, 1998) and is got, and promptly is respectively the α of each neutral gas composition and temperature T correspondence in the table
1Value;
For G
i(L), be described as follows:
G
5(L)=0,
G
1(L), G
2(L), G
3(L), G
4(L) and G
T(L) unified below with G (L) expression, because their difference is with top A
1, A
2, A
3, A
4And A
TSimilar, only be the value of when looking into DTM94 model coefficient table, looking into each neutral gas composition and temperature T correspondence respectively.
G(L)=LAT+F+M+β(AN1+SAN1+AN2+SAN2+D+SD+TD);
Wherein,
The latitude item:
LAT=a
2H
20+a
3H
40+a
37H
10;
The solar activity item:
The geomagnetic activity item:
Symmetry anniversary item:
AN1=(a
9+a
10H
20)cos[Ω(d-a
11)];
Symmetry item half a year:
SAN1=(a
12+a
13H
20)cos[Ω(d-a
14)];
Asymmetric anniversary item:
AN2=(a
15H
10+a
16H
30+a
17H
50)cos[Ω(d-a
18)];
Asymmetric half a year item:
SAN2=a
19H
10cos[2Ω(d-a
20)];
Sunday item:
D=(-B
1cosλ
s+B
2sinλ
s)θ
1-(B
1sinλ
s+B
2cosλ
s)θ
2;
Half cycle day item:
1/3 Sunday item:
Coefficient is described as follows in the top formula:
B
1=a
21H
11+a
22H
31+a
23H
51+(a
24H
11+a
25H
21)cos[Ω(d-a
18)],
B
2=a
26H
11+a
27H
31+a
28H
51+(a
29H
11+a
30H
21)cos[Ω(d-a
18)],
B
3=a
31H
22+a
32H
32cos[Ω(d-a
18)],
B
4=a
33H
22+a
34H
32cos[Ω(d-a
18)],
B
5=a
35H
33,
B
6=a
36H
33;
a
i(i=2,3 ... 39) value can look into DTM94 model coefficient table (Berger et, al.1998);
The serve as reasons day of year number of year first calculation of d;
F
10.7Be sun 10.7cm radiant flux, with d-1 days values;
Be sun 10.7cm radiant flux mean value, average in three solar rotation cycles (81 days) before the d;
K
ρBe 3 hours geomagnetic indicess, get calculate constantly before
Hour value, wherein
Being the geocentric latitude of sub-satellite point, is unit with the degree;
Ω=2π/365(day
-1);
H
10=θ
3,
H
11=1
H
21=3θ
3,
H
22=3,
H
32=15θ
3,
H
33=15
About being described as follows of coordinate amount:
Satellite the earth's core square r and direction cosine θ
1, θ
2, θ
3With the pass of its spherical coordinates be:
X, y, z are the square coordinate of satellite in body-fixed coordinate system, and they can calculate in the sky, the two poles of the earth.
With
The expression square coordinate of the sun in body-fixed coordinate system, the pass of it and spherical coordinates is:
As a further improvement on the present invention, in the described Atmospheric models, further comprise the computing method of atmospheric density with respect to ground solid square coordinate partial derivative, atmospheric density with respect to the computing method of ground solid square coordinate partial derivative is:
If X
1=x, X
2=y, X
3=z represents the square coordinate of satellite in body-fixed coordinate system, then has:
For j=1,2:
For j=3:
In the top formula,
B′
1=a
22H′
31+a
23H′
51+a
25H′
21cos[Ω(d-a
18)]
B′
2=a
27H′
31+a
28H′
51+a
30H′
21cos[Ω(d-a
18)]
B′
3=a
32H′
32cos[Ω(d-a
18)]
B′
4=a
34H′
32cos[Ω(d-a
18)]
H?′
10=1
H′
20=3θ
3
H′
21=3
H′
31=15θ
3
H′
32=15
The coefficient that does not directly illustrate in the above formula all has explanation in atmospheric density is calculated.
Beneficial effect of the present invention is the hot Atmospheric models in conjunction with DTM94, provided its atmospheric density and atmospheric density no singular point computing formula with respect to ground solid square coordinate partial derivative, for versatility and the robustness that improves low orbit satellite orbit determination software provides a practical tool, computing method of the present invention are expressed the no singular point of other Atmospheric models has reference function too.
Embodiment
Be described further below in conjunction with embodiment.
The present invention is successful Application in phased-array radar high-precision orbit software, and the method for the no singular point transformation of process makes total system to move continuously incessantly.
Atmospheric models in the phased-array radar high-precision orbit, the expression-form of atmospheric density ρ is:
Subscript i=1,2 ... 5, the corresponding H of difference, H
e, O, N
2, O
2This five kinds of gas atoms or molecule.
In the following formula,
A
VR=1/ (6.022 * 10
23) be the inverse of avogadros constant;
M
iBe the molecular weight of neutral gas composition i, M
1=1, M
2=4, M
3=16, M
4=28, M
5=32;
n
i(h)=A
iExp (G
i(L)) f
i(h), h is the floor level of satellite;
Following formula the right symbol description is as follows:
This wherein,
T
120=380K represents 120 kilometers highly temperature at place,
T
∞=A
T(1+G
T(L)), the expression Outer Atmospheric Temperature,
α
iBe the thermal diffusion coefficient of neutral gas composition i, α
1=-0.40, α
2=-0.38, α
3=α
4=α
5=0.0,
γ
i=M
ig
120/(σRT
∞),
R=831.4 is a gas law constant,
g
120=944.6626 (cm/s
2), the gravitational acceleration at 120 kilometers places of expression floor level,
σ=T '
120/ (T
∞-T
120), represent relative vertical temperature gradient,
T '
120=14.348, the thermograde at 120 kilometers places of expression,
R
ρ=6356.77km, expression earth polar radius,
Expression satellite ground height,
e
2=2 ε-ε
2, expression earth meridian circle excentricity,
ε=1/298.25, expression earth meridian circle ellipticity.
And other coefficients:
A
5=4.775×10
10(cm
-3),
A
1, A
2, A
3, A
4And A
TValue is looked in can the coefficient table by DTM94 model (Berger et al, 1998) and is got, and promptly is respectively the α of each neutral gas composition and temperature T correspondence in the table
1Value;
For G
i(L), be described as follows:
G
5(L)=0,
G
1(L), G
2(L), G
3(L), G
4(L) and G
T(L) unified below with G (T) expression, because their difference is with top A
1, A
2, A
3, A
4And A
TSimilar, only be the value of when looking into DTM94 model coefficient table, looking into each neutral gas composition and temperature T correspondence respectively.
G(L)=LAT+F+M+β(AN1+SAN1+AN2+SAN2+D+SD+TD);
Wherein,
The latitude item:
LAT=a
2H
20+a
3H
40+a
37H
10;
The solar activity item:
The geomagnetic activity item:
Symmetry anniversary item:
AN1=(a
9+a
10H
20)cos[Ω(d-a
11)];
Symmetry item half a year:
SAN1=(a
12+a
13H
20)cos[Ω(d-a
14)];
Asymmetric anniversary item:
AN2=(a
15H
10+a
16H
30+a
17H
50)cos[Ω(d-a
18)];
Asymmetric half a year item:
SAN2=a
19H
10cos[2Ω(d-a
20)];
Sunday item:
D=(-B
1cosλ
s+B
2sinλ
s)θ
1-(B
1sinλ
s+B
2cosλ
s)θ
2;
Half cycle day item:
1/3 Sunday item:
Coefficient is described as follows in the top formula:
B
1=a
21H
11+a
22H
31+a
23H
51+(a
24H
11+a
25H
21)cos[Ω(d-a
18)],
B
2=a
26H
11+a
27H
31+a
28H
51+(a
29H
11+a
30H
21)cos[Ω(d-a
18)],
B
3=a
31H
22+a
32H
32cos[Ω(d-a
18)],
B
4=a
33H
22+a
34H
32cos[Ω(d-a
18)],
B
5=a
35H
33,
B
6=a
36H
33;
a
i(i=2,3 ... 39) value can look into DTM94 model coefficient table (Berger et, al.1998);
The serve as reasons day of year number of year first calculation of d;
F
10.7Be sun 10.7cm radiant flux, with d-1 days values;
Be sun 10.7cm radiant flux mean value, average in three solar rotation cycles (81 days) before the d;
K
ρBe 3 hours geomagnetic indicess, get calculate constantly before
Hour value, wherein
Being the geocentric latitude of sub-satellite point, is unit with the degree;
Ω=2π/365(day
-1);
H
10=θ
3,
H
11=1
H
21=3θ
3,
H
22=3,
H
32=15θ
3,
H
33=15
About being described as follows of coordinate amount:
Satellite the earth's core square r and direction cosine θ
1, θ
2, θ
3With the pass of its spherical coordinates be:
X, y, z are the square coordinate of satellite in body-fixed coordinate system, and they can calculate in the sky, the two poles of the earth.
With
The expression square coordinate of the sun in body-fixed coordinate system, the pass of it and spherical coordinates is:
As a further improvement on the present invention, in the described Atmospheric models, further comprise the computing method of atmospheric density with respect to ground solid square coordinate partial derivative, atmospheric density with respect to the computing method of ground solid square coordinate partial derivative is:
If X
1=x, X
2=y, X
3=z represents the square coordinate of satellite in body-fixed coordinate system, then has:
For j=1,2:
For j=3:
In the top formula,
B′
1=a
22H′
31+a
23H′
51+a
25H′
21cos[Ω(d-a
18)]
B′
2=a
27H′
31+a
28H′
51+a
30H′
21cos[Ω(d-a
18)]
B′
3=a
32H′
32cos[Ω(d-a
18)]
B′
4=a
34H′
32cos[Ω(d-a
18)]
H′
10=1
H′
20=3θ
3
H′
21=3
H′
31=15θ
3
H′
32=15
The coefficient that does not directly illustrate in the above formula all has explanation in atmospheric density is calculated.
Claims (2)
1. the Atmospheric models in the phased-array radar high-precision orbit is characterized in that, the expression-form of atmospheric density ρ is:
Subscript i=1,2 ... 5, the corresponding H of difference, H
e, O, N
2, O
2This five kinds of gas atoms or molecule;
In the following formula,
A
VR=1/ (6.022 * 10
23) be the inverse of avogadros constant;
M
iBe the molecular weight of neutral gas composition i, M
1=1, M
2=4, M
3=16, M
4=28, M
5=32;
n
i(h)=A
iExp (G
i(L)) f
i(h), h is the floor level of satellite;
Following formula the right symbol description is as follows:
This wherein,
T
120=380K represents 120 kilometers highly temperature at place,
T
∞=A
T(1+G
T(L)), the expression Outer Atmospheric Temperature,
α
iBe the thermal diffusion coefficient of neutral gas composition i, α
1=-0.40, α
2=-0.38, α
3=α
4=α
5=0.0,
γ
i=M
ig
120/(σRT
∞),
R=831.4 is a gas law constant,
g
120=944.6626 (cm/s
2), the gravitational acceleration at 120 kilometers places of expression floor level,
σ=T
120/ (T-T
120), represent relative vertical temperature gradient,
T '
120=14.348, the thermograde at 120 kilometers places of expression,
R
ρ=6356.77km, expression earth polar radius,
e
2=2 ε-ε
2, expression earth meridian circle excentricity,
ε=1/298.25, expression earth meridian circle ellipticity;
And other coefficient:
A
5=4.775×10
10(cm
-3),
A
1, A
2, A
3, A
4And A
TValue is looked in can the coefficient table by the DTM94 model and is got, and promptly is respectively the α of each neutral gas composition and temperature T correspondence in the table
1Value;
For G
i(L), be described as follows:
G
5(L)=0,
G
1(L), G
2(L), G
3(L), G
4(L) and G
T(L) unified below with G (L) expression, because their difference is with top A
1, A
2, A
3, A
4And A
TSimilar, only be the value of when looking into DTM94 model coefficient table, looking into each neutral gas composition and temperature T correspondence respectively;
G(L)=LAT+F+M+β(AN1+SAN1+AN2+SAN2+D+SD+TD);
Wherein,
The latitude item:
LAT=a
2H
20+a
3H
40+a
37H
10;
The solar activity item:
The geomagnetic activity item:
Symmetry anniversary item:
AN1=(a
9+a
10H
20)cos[Ω(d-a
11)];
Symmetry item half a year:
SAN1=(a
12+a
13H
20)cos[Ω(d-a
14)];
Asymmetric anniversary item:
AN2=(a
15H
10+a
16H
30+a
17H
50)cos[Ω(d-a
18)];
Asymmetric half a year item:
SAN2=a
19H
10cos[2Ω(d-a
20)];
Sunday item:
D=(-B
1cosλ
s+B
2sinλ
s)θ
1-(B
1sinλ
s+B
2cosλ
s)θ
2;
Half cycle day item:
1/3 Sunday item:
Coefficient is described as follows in the top formula:
B
1=a
21H
11+a
22H
31+a
23H
51+(a
24H
11+a
25H
21)cos[Ω(d-a
18)],
B
2=a
26H
11+a
27H
31+a
28H
51+(a
29H
11+a
30H
21)cos[Ω(d-a
18)],
B
3=a
31H
22+a
32H
32cos[Ω(d-a
18)],
B
4=a
33H
22+a
34H
32cos[Ω(d-a
18)],
B
5=a
35H
33,
B
6=a
36H
33;
a
i(i=2,3 ... 39) value can be looked into DTM94 model coefficient table;
The serve as reasons day of year number of year first calculation of d;
F
10.7Be sun 10.7cm radiant flux, with d-1 days values;
K
ρBe 3 hours geomagnetic indicess, get calculate constantly before
Hour value, wherein
Being the geocentric latitude of sub-satellite point, is unit with the degree;
Ω=2π/365(day
-1);
H
10=θ
3,
H
11=1
H
21=3θ
3,
H
22=3,
H
32=15θ
3,
H
33=15
About being described as follows of coordinate amount:
Satellite the earth's core square r and direction cosine θ
1, θ
2, θ
3With the pass of its spherical coordinates be:
X, y, z are the square coordinate of satellite in body-fixed coordinate system;
With
The expression square coordinate of the sun in body-fixed coordinate system, the pass of it and spherical coordinates is:
2. the Atmospheric models in the phased-array radar high-precision orbit according to claim 1, it is characterized in that, described Atmospheric models further comprise the expression-form of atmospheric density with respect to ground solid square coordinate partial derivative, and atmospheric density with respect to the expression-form of ground solid square coordinate partial derivative is:
If X
1=x, X
2=y, X
3=z represents the square coordinate of satellite in body-fixed coordinate system, then has:
For j=1,2:
For j=3:
In the top formula,
B′
1=a
22H′
31+a
23H′
51+a
25H′
21cos[Ω(d-a
18)]
B′
2=a
27H′
31+a
28H′
51+a
30H′
21cos[Ω(d-a
18)]
B′
3=a
32H′
32cos[Ω(d-a
18)]
B′
4=a
34H′
32cos[Ω(d-a
18)]
H′
10=1
H′
20=3θ
3
H′
21=3
H′
31=15θ
3
H′
32=15
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102378952A CN101923169A (en) | 2010-07-27 | 2010-07-27 | Atmospheric model for high-precision orbit confirmation of phased array radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102378952A CN101923169A (en) | 2010-07-27 | 2010-07-27 | Atmospheric model for high-precision orbit confirmation of phased array radar |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101923169A true CN101923169A (en) | 2010-12-22 |
Family
ID=43338208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102378952A Pending CN101923169A (en) | 2010-07-27 | 2010-07-27 | Atmospheric model for high-precision orbit confirmation of phased array radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101923169A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540180A (en) * | 2012-01-02 | 2012-07-04 | 西安电子科技大学 | Space-based phased-array radar space multi-target orbit determination method |
CN109323698A (en) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | Space target meteorology multi-model tracking and guiding technology |
-
2010
- 2010-07-27 CN CN2010102378952A patent/CN101923169A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540180A (en) * | 2012-01-02 | 2012-07-04 | 西安电子科技大学 | Space-based phased-array radar space multi-target orbit determination method |
CN109323698A (en) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | Space target meteorology multi-model tracking and guiding technology |
CN109323698B (en) * | 2018-12-03 | 2021-05-11 | 中科星图(西安)测控技术有限公司 | Space target merle multi-model tracking and guiding method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104792340B (en) | A kind of star sensor installation error matrix and navigation system star ground combined calibrating and the method for correction | |
Edwards et al. | Maximizing net power in circular turns for solar and autonomous soaring aircraft | |
CN102322862B (en) | Method for determining absolute orbit and relative orbit of formation flight satellite | |
CN102175260B (en) | Error correction method of autonomous navigation system | |
CN104165640A (en) | Near-space missile-borne strap-down inertial navigation system transfer alignment method based on star sensor | |
CN107609267B (en) | A kind of moon Finite Thrust repeatedly captures track implementation method | |
CN108181916A (en) | The control method and device of moonlet relative attitude | |
CN104484493B (en) | A kind of airship returns to predetermined drop point regression orbit method for designing | |
Zhao et al. | Minimum-fuel station-change for geostationary satellites using low-thrust considering perturbations | |
CN102923323A (en) | Design method for low-energy transit among interplanetary fixed orbits based on invariant manifold | |
CN104932546B (en) | Ship borne type self-adapting solar energy cell panel loading device and adaptive approach | |
Horsley | An investigation into using differential drag for controlling a formation of CubeSats | |
CN102116630A (en) | Mars probe on-board quick and high-precision determination method | |
CN101923169A (en) | Atmospheric model for high-precision orbit confirmation of phased array radar | |
He et al. | Solution set calculation of the Sun-perturbed optimal two-impulse trans-lunar orbits using continuation theory | |
CN103438892B (en) | A kind of astronomical autonomous orbit determination algorithm based on EKF of improvement | |
Pacheco et al. | Superrotation of the ionosphere and quiet time zonal ion drifts at low and middle latitudes observed by Republic of China Satellite‐1 (ROCSAT‐1) | |
CN103274066B (en) | Design method of escape orbit starting from Halo track and used for detecting deep space target | |
CN106326576A (en) | Yaw estimating method for whole spacecraft offset angle momentum under any reference system | |
Liu et al. | Mission-based optimization of insulation layer for the solar array on the stratospheric airship | |
Wang et al. | Cloud‐tracked winds for the first Mars Global Surveyor mapping year | |
Gong et al. | Solar sail transfer trajectory from L1 point to sub-L1 point | |
CN103149936B (en) | Combined attitude determination method for UPF (user port function) algorithm optimized by DNA (deoxyribonucleic acid) algorithm | |
Chang et al. | Trajectory Optimization for Dynamic Soaring Remotely Piloted Aircraft with Under-Wing Solar Panels | |
CN106595657B (en) | Axial symmetry aircraft attitude measurement device and its measurement method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20101222 |