CN105183959B - A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability - Google Patents
A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability Download PDFInfo
- Publication number
- CN105183959B CN105183959B CN201510527488.8A CN201510527488A CN105183959B CN 105183959 B CN105183959 B CN 105183959B CN 201510527488 A CN201510527488 A CN 201510527488A CN 105183959 B CN105183959 B CN 105183959B
- Authority
- CN
- China
- Prior art keywords
- satellite
- antenna
- angle
- blocked
- passes
- 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.)
- Active
Links
Abstract
The method that a kind of number obtained when satellite is in-orbit to carry out real pass of the present invention passes antenna occlusion probability, first, the sum of the subtended angle over the ground biography antenna of satellite, which is calculated, will not be blocked maximum lateral swinging angle；Then, STK simulating scenes are established, number is drawn and passes the region that antenna may be blocked；Secondly, probability is blocked using area ratio and angle ratio, calculating；Finally, according to result of calculation, the formulation that number passes antenna occlusion probability level is provided.The present invention proposes a kind of new thought analyzed number and pass antenna occlusion probability, area is created than the concept with angle ratio, and give using area than calculating the method that number passes antenna occlusion probability with angle ratio, with reference to the powerful two-dimentional three dimensional display capabilities of STK, realize that logarithm passes the analysis of antenna occlusion probability using very small amount of calculation, satellite user and satellite development side is contributed to assess and examine satellite operation on orbit efficiency, formulating index of correlation for satellite user and satellite development side provides important evidence.
Description
Technical field
The invention belongs to spacecraft mission planning field, and in particular to a kind of number biography obtained when the in-orbit progress of satellite is real to be passed
The method of antenna occlusion probability.
Background technology
Remote sensing satellite, particularly high-resolution earth observation Optical remote satellite, substantial amounts of figure will be produced during its work
As data.But the restriction of satellite volume, quality, power consumption etc. is limited by, it is used for preserving the solid-state memory of data on star
The amount of capacity of (referred to as depositing admittedly) is severely limited.Therefore, using camera while being imaged, number passes antenna (hereinafter referred to as over the ground
Number passes antenna) while the real imaging pattern that passes of earth station's passback is taken seriously.In order to ensure enough passback speed and transmission
Gain, it is desirable to which number passes antenna during passing in fact, is realized and rotated by the driving of itself servo control mechanism, allows antenna beam to refer to all the time
Ground station.Enter because the imaging of high-resolution earth observation remote sensing satellite is substantially in the case where posture side-sway biases offline mode
OK, therefore the real mode of operation that passes also generally all is that side-sway passes in fact.
In order to ensure the enough energy, satellite needs to install solar energy in large area windsurfing, and windsurfing is in satellite transit process
It is middle to be rotated according to the incident angle of sunshine, to improve charge efficiency.In order to prevent influenceing antenna, over the ground optics over the ground
The working field of view over the ground of the equipment such as camera, windsurfing are often installed on the left and right sides of satellite, and its installation axle is parallel to track system Y
Axle, it is vertical with satellite transit direction with plane-parallel, it is vertical with satellite orbit face.But when satellite is inclined in wide-angle side-sway
When putting flight, windsurfing will carry out wide-angle deflection together with celestial body.Because windsurfing installation axle is parallel to track system Y-axis, then side
Side windsurfing can be upturned away from ground direction during pendulum, and opposite side windsurfing can bow downwards close to ground direction.Closely
This side windsurfing in face, it is possible to which the Wave beam forming that meeting logarithm passes antenna blocks, and makes it can not be during passing in fact with " seeing "
Face station, so as to which real biography task can not be completed.
Because satellite sailboard area is huge, and side-sway in fact pass be high-resolution earth observation Optical remote satellite can not or
Scarce operation on orbit mode, therefore eclipse phenomena is very difficult to avoid completely.And from the point of view of conventional type design experience, pass through satellite sheet
The optimization of body configuration design can not settle the matter once and for all.
Satellite user can receive the presence of eclipse phenomena, it may be desirable to it is as low as possible to block probability.Therefore need by imitative
The means such as true analysis, calculate the in-orbit number carried out when side-sway passes in fact of satellite and pass antenna occlusion probability, and as satellite system
System index, period is examined by the method for long-time statistical in orbit after satellite launch.
The simulating analysis that conventional logarithm passes antenna occlusion probability was more inflexible, i.e., complete mechanical according to satellite structure
Type design result, the largest enveloping condition that an antenna will not be blocked is provided, without considering any other factor.It is such
As a result, cause to analyze come to block probability often very high, be not substantially inconsistent with real in-orbit actual conditions.
The content of the invention
The technical problems to be solved by the invention are：Overcome the deficiencies in the prior art, there is provided one kind obtain satellite it is in-orbit enter
The method for the number biography antenna occlusion probability gone when passing in fact, can rapidly be calculated and block generally according to the configuration design point of satellite
Rate.
The technical scheme is that：A kind of number obtained when the in-orbit progress of satellite is real to be passed passes the side of antenna occlusion probability
Method, step are as follows：
1) satellite maximum lateral swinging angle is obtained, and number passes antenna not by the maximum lateral swinging angle θ under circumstance of occlusion；
2) calculate and obtain the maximum lateral swinging angle that the sum of the subtended angle over the ground biography antenna of satellite will not be blocked, judge maximum side
Whether pivot angle, which is more than number, passes the maximum lateral swinging angle θ that antenna will not be blocked, if the satellite maximum side-sway obtained in step 1)
Angle is more than counted number in step 2) and passes the maximum lateral swinging angle that antenna will not be blocked, then into step 3)；Otherwise, show to count
Passing antenna will not be blocked certainly, and it is 0% to block probability, now obtain number when satellite side-sway passes in fact and pass antenna occlusion probability
Number passes antenna and will not be blocked when index, i.e. side-sway pass in fact；
3) viewing area of satellite and earth station is drawn in STK；
4) three-dimensional feature needed for being set in STK
41) angle that number passes antenna zero vector and satellite-earth station's vector is drawn in STK；
42) satellite body system YZ planes are set in STK；
5) number passes the region that may be blocked of antennas and is divided into two pieces, respectively positioned at Satellite offline east side and west side,
It is respectively defined as eastern section and west side that number passes the region that antenna may be blocked；The lateral swinging angle that satellite is set in STK is number
The maximum positive lateral swinging angle that antenna will not be blocked is passed, draws the eastern section that number passes the region that antenna may be blocked；
6) the maximum minus side pivot angle for setting the attitude of satellite not to be blocked for number biography antenna in STK, draws number and passes day
The west section in the region that line may be blocked；
7) calculate number when satellite side-sway passes in fact and pass antenna occlusion probability
71) calculate and obtain area ratio
Determine number that step 6) obtains pass antennas may be blocked region eastern section and west section area, calculate eastern section and west section
The area in region accounts for the ratio of whole earth station viewing area, obtain area than calculation formula it is as follows：
Area ratio≤(eastern section area+west section area)/ground gross area
72) calculate and obtain angle ratio
The maximum side that will not be blocked using the number biography antenna in the satellite maximum lateral swinging angle and step 2) in step 1)
Pivot angle, calculating acquisition angle ratio≤(satellite maximum lateral swinging angle-number passes the maximum lateral swinging angle that antenna will not be blocked)/satellite is most
Big lateral swinging angle=1- numbers pass maximum lateral swinging angle/satellite maximum lateral swinging angle that antenna will not be blocked；
73) obtain and block Probability p=area than * angles ratio≤(eastern section area+west section area)/ground gross area * (1- numbers
Pass maximum lateral swinging angle/satellite maximum lateral swinging angle that antenna will not be blocked)；
74) now obtain number when satellite side-sway passes in fact and pass antenna occlusion probability level, i.e., when satellite side-sway angle absolute value
When being spent no more than θ, number passes antenna and will not be blocked when side-sway passes in fact, and when when satellite side-sway angle, absolute value is spent more than θ, side-sway is real
The probability that number biography antenna is blocked during biography is no more than p.
The sum of subtended angle over the ground of step 2) Satellite passes the maximum lateral swinging angle circular that antenna will not be blocked
For：
Subtended angle=asin (earth radius/(satellite orbital altitude+earth radius))+surplus angle over the ground；Number passes antenna not
The maximum lateral swinging angle that can be blocked=number passes antenna not by the subtended angle over the ground of hard-over-satellite under circumstance of occlusion.
The specific method for drawing the viewing area of satellite and earth station in step 3) in STK is as follows：
31) simulating scenes are established；
32) satellite orbit is established；
33) earth station is established；
34) earth station antenna is established；
35) viewing area of satellite and earth station is shown.
The present invention compared with prior art the advantages of be：
(1) biggest advantage of the present invention, it is to pass antenna occlusion probability analysis method in traditional number to only account for satellite structure
On the basis of type design, satellite and earth station's relative position relation have further been considered, systematically analyzed when blocking generation
Prerequisite condition element so that the analysis result significantly in-orbit real work state of close satellite；
(2) the maximum bright spot of the present invention is to create area than the concept with angle ratio., can be more by area ratio
The probability for first precondition --- occlusion area --- that antenna occlusion occurs is given exactly；, can be with by angle ratio
Relatively accurately give the probability for another precondition --- side-sway angle --- that antenna occlusion occurs.Area ratio and angle
The synthesis of ratio, also can, which provides to count, passes the probability that antenna blocks；
(3) present invention is simple to operate, and for computing function powerful full utilization STK as instrument, amount of calculation is very small, is not required to
Program, it is only necessary to which the input parameter for carrying out simulation calculation is provided.During analysis, by constantly micro- manually in STK
Adjust scenario parameters and observe the satellitosis situation of change in two-dimentional Three-dimensional Display frame, obtain analysis result；
(4) The present invention gives a kind of formulation of new antenna occlusion probability level, the formulation is relative to conventional index
For formulation, description is more accurate；
(5) present invention has preferable general applicability.Function is passed in fact with side-sway present invention may apply to various
Satellite carries out related simulation analysis, and has been verified in the development process for having launched satellite model.
Brief description of the drawings
Fig. 1 is that number of the present invention passes the FB(flow block) that antenna rotates planing method.
Embodiment
As shown in figure 1, it is the FB(flow block) of the inventive method, details are as follows for specific steps：
1) satellite maximum lateral swinging angle is obtained, and number passes antenna not by the hard-over under circumstance of occlusion
Satellite maximum lateral swinging angle is proposed by satellite user as one of most important system index of satellite, is design of satellites
Important input parameter.
Number passes antenna not by the hard-over under circumstance of occlusion, refers to when number passes antenna rotation and is no more than the angle, its
Rotate anyway, all without being blocked.The angle is bigger, and it is lower that number passes the probability that antenna is blocked.The angle and satellite
Configuration is related in itself, is provided by satellite configuration design engineer after satellite configuration design is completed.
Therefore, this step two angles to be obtained, all it is as the invention method input to be used ginseng
Number.
2) sum of subtended angle over the ground for calculating satellite passes the maximum lateral swinging angle that antenna will not be blocked
21) subtended angle over the ground of satellite is calculated
The subtended angle over the ground of satellite, refer to the angle between satellite-the earth's core vector satellite-earth's surface point of contact vector, be one
Acute angle.The angle is only relevant with satellite orbital altitude, and orbit altitude is higher, and subtended angle is smaller over the ground.
Calculation formula is as follows：
Subtended angle over the ground
=asin (earth radius/(satellite orbital altitude+earth radius))
But consider earth surface characteristic, with reference to conventional engineering experience, typically require on the basis of above-mentioned result of calculation
The certain surplus angle of upper increase, generally takes 2~3 degree.Then above-mentioned formula is changed into：
Subtended angle over the ground
=asin (earth radius/(satellite orbital altitude+earth radius))+surplus angle
22) calculate number and pass the maximum lateral swinging angle that antenna will not be blocked
Number passes the maximum lateral swinging angle that antenna will not be blocked, and refers to when satellite side-sway is no more than the angle, and number passes day
Line rotates anyway, and its wave beam is all without being blocked.The angle is less than or equal to satellite maximum lateral swinging angle, and is defended when it is closer
During star maximum lateral swinging angle, number passes the probability that antenna is blocked will be lower.
Calculation formula is as follows：
Number passes the maximum lateral swinging angle that antenna will not be blocked
=number passes antenna not by the subtended angle over the ground of hard-over-satellite under circumstance of occlusion
23) judge whether maximum lateral swinging angle is more than number and passes the maximum lateral swinging angle that antenna will not be blocked
If " the satellite maximum lateral swinging angle " that obtains in step 1) is more than, counted in step 22) " number passes antenna will not be by
To the maximum lateral swinging angle blocked ", then into step 3)；Otherwise, then mean that number passes antenna and will not be blocked certainly, block general
Rate is 0%, is directly entered step 8).
Block probability for 0% when, number when satellite side-sway passes in fact can be provided and pass antenna occlusion probability level, as defending
The system index of star design.Its index formulation is as follows：
● number passes antenna and will not be blocked when side-sway passes in fact；
3) viewing area of satellite and earth station is drawn in STK
According to satellite orbit parameter and the geographical location information of earth station, pass through STK softwares (Satellite ToolKit)
Draw the viewing area of satellite and earth station.STK is that a aerospace Simulation analyzes business software, has two-dimentional three-dimensional visible circle
Face, powerful, easy to operate, in the world and China's space industry widely uses.
It is that STK versions are 8.1 used in patent of the present invention.The STK of different editions operationally can be distinguished slightly.Due to
The main contents of this and non-invention patent, will not be repeated here.With the process of STK acquisition satellites and earth station viewing area such as
Under：
31) simulating scenes are established
After opening STK softwares, the File-New on click tools column, a simulating scenes are generated.
32) satellite orbit is established
Satellite orbit is established in this step, its detailed process is as follows：
A. the Insert-New on click tools column, selects Satellite in an open window, and generation one is newly defended
Star.Cancel is selected in an open window；
B. the satellite just generated is double-clicked in the Object Browser columns in left side, satellite attribute is opened and window is set；
C. Basic-Orbit is clicked in attribute sets window, opens track and interface is set.In Semimajor Axis
The design load of satellite orbit semi-major axis is inserted, it is other to use default value；
D. 2D Graphics-Pass are clicked in attribute sets window, opens two-dimentional window and interface is set.In Visible
Descending is selected in Sides, i.e., only rail section drops in display.Click on OK and close window；
E. now in 2D Graphics interfaces i.e. it can be seen that the satellite orbit just established；
F. in 3D Graphics interfaces, the View From/To buttons of left side eye shape are clicked on, open View
From/To for Scenario1 windows, the satellite just established is selected in View From windows.Selection lower section
Unconstrained Rotation options, Show " Track Ball " options are not selected.Click on OK and close window；
G. now in 3D Graphics i.e. it can be seen that the dummy satellite just established.It can be put using mouse right and left key
It is big to reduce visual field or adjustment visual field orientation；
33) earth station is established
Earth station is established in this step, its detailed process is as follows：
A. the Insert-New on click tools column, Facility is selected in an open window, generate a new ground
Stand；
B. the earth station just generated is double-clicked in the Object Browser columns in left side, earth station's attribute is opened and window is set
Mouthful；
C. Basic-Position is clicked in attribute sets window, opens earth station location and interface is set.
Latitude, the longitude and altitude of earth station are inserted in Latitude, Longitude and Altitude, it is other to use default value.Click on
OK。
D. now in 2D Graphics i.e. it can be seen that the earth station just established；
34) earth station antenna is established
Earth station antenna is established in this step, its detailed process is as follows：
A. the earth station just generated is clicked
B. the Insert-New on click tools column, Sensor is selected in an open window, generate a new antenna；
C. the antenna just generated is double-clicked in the Object Browser columns in left side, antenna property is opened and window is set；
D. Basic-Definition is clicked in attribute sets window, opens antenna basic setup window.In Sensor
Simple Conic are selected in Type, the design load of earth station antenna visual field is set in Cone Angle.Click on OK；
35) viewing area of satellite and earth station is shown
Satellite and earth station (being in fact earth station antenna, acute pyogenic infection of finger tip is hereafter simplified with earth station) are shown in this step
Viewing area.Its detailed process is as follows：
A. the antenna just generated is double-clicked in the Object Browser columns in left side, antenna property is opened and window is set；
B. 2D Graphics-Position are clicked in attribute sets window, opens day line projection and window is set.
Object Altitude are selected in Project to columns, select to establish just now in Projection Altitude Object
Satellite.Click on OK；
C. the viewing area of satellite and earth station can be now shown in 2D Graphics interfaces, the region is one
It is circular.
According to above step, you can draw the viewing area of satellite and earth station.Typically, required according to crew station design,
When satellite is in the viewing area, earth station antenna is certainly it can be seen that satellite；Meanwhile required according to design of satellites, when
When satellite is in the viewing area, if not considering, windsurfing blocks, and the number on satellite passes antenna also certainly it can be seen that ground
Stand.
4) three-dimensional features such as required vector, angle and plane are set in STK
STK softwares have powerful three dimensional display capabilities.This step in STK three-dimensional scenics by by adding Computer Aided Design
The three-dimensional features such as required vector, angle and plane, it is ready the step of come for behind the inventive method.Detailed process is such as
Under：
41) angle that number passes antenna zero vector and satellite-earth station's vector is drawn in STK
The number to be drawn in this step passes the angle of antenna zero vector and satellite-earth station's vector, and its detailed process is such as
Under：
A. in STK, the satellite to be analyzed is clicked in the Object Browser columns of left side, then up toolbar
Upper selection Satellite-Vector Geometry, open " Vector Geometry Tool for satellites name " window；
B. in an open window, right side Create New Vector buttons are clicked on.In an open window, in Name
Column gives a vector name；In Type columns selection Fixed in Axes；The Body of satellite is selected in Reference Axes
(namely body coordinate system)；Cartesian is selected in Type；Number is filled out in X, Y, Z and passes antenna zero vector in whole star
Component under body coordinate system, the component are design of satellites value, typically can all be designed to (0,0,1), i.e., number passes antenna and is in zero
During position substar is pointed to along satellite+Z axis.OK is clicked on, establishes vector completion.Here what is established is that number passes antenna zero vector；
C. still in an open window, right side Create New Vector buttons are clicked on.In the given arrow in Name columns
Measure name；Displacement is selected on Type columns；In the Center of Origin Point selection satellites；In Destination
Point selects the Center of earth station.OK is clicked on, establishes vector completion.Here what is established is satellite-earth station's vector；
D. still in an open window, right side Create New Angle buttons are clicked on.In an open window, exist
Name columns give an angle name；In Type columns selection Between Vectors；Step b is selected to build in From Vector
Vertical number passes antenna zero vector；Satellite-earth station's the vector for selecting step c to establish in To Vector.OK is clicked on, is established
Angle is completed.Here what is established is the angle of above-mentioned two vector；
E. Close is clicked on, closes " Vector Geometry Tool for satellites name " window；
F. in STK, the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and sets
Window；
G. 3D Graphics-Vector are clicked in attribute sets window, Add is clicked in the interface, what is newly opened
Two vector angles that satellite-earth station vector step d of step c foundation is established are double-clicked in window respectively, OK is clicked on and closes window
Mouthful.Here it is that two vector angles that the satellite-earth station vector step d for establishing step c is established are added to 3D Graphics
In；
H. satellite-earth station of step g additions is found in the Name columns in the 3D Graphics-Vector windows upper left corner
Vector, its corresponding square frame in the show columns of right side is clicked, occurred in square frame one small to hooking.Here it is to add step g
Satellite-earth station's vector shown in 3D Graphics；
I. the angle of step g additions is found in the Name columns in the 3D Graphics-Vector windows upper left corner, clicks it
The corresponding square frame in the show columns of right side, occur in square frame one small to hooking.Here it is in 3D fields by the angle that step g is added
Shown in scape；
J. the square frame on the left of Show Angle Value columns is clicked on right side, is occurred in square frame one small to hooking.Here
It is to show the angle value of the step i angles shown in 3D scenes；
K. OK is clicked on.It can now show that satellite-earth station's vector, and number pass antenna zero in 3D Graphics
The angle and its angle value of vector and satellite-earth station's vector.
42) satellite body system YZ planes are set in STK
The satellite body system YZ planes to be drawn in this step, its detailed process are as follows：
A. in STK, the satellite to be analyzed is clicked in the Object Browser columns of left side, then up toolbar
Upper selection Satellite-Vector Geometry, open " Vector Geometry Tool for satellites name " window；
B. in an open window, right side Create New Plane buttons are clicked on.In an open window, on Name columns
Give a plane name；Quadrant is selected on Type columns；YZ is selected in Quadrant columns；In Reference System
The Body (namely body coordinate system) of middle selection satellite.OK is clicked on, establishes plane completion.Here what is established is satellite body system
YZ planes；
C. Close is clicked on, closes " Vector Geometry Tool for satellites name " window；
D. in STK, the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and sets
Window；
E. 3D Graphics-Vector are clicked in attribute sets window, Add is clicked in the interface, what is newly opened
The satellite body system YZ planes that step b is established are double-clicked in window, OK is clicked on and closes window.Here it is the satellite for establishing step b
Body series YZ planes are added in 3D Graphics；
F. the satellite body system YZ of step e additions is found in the Name columns in the 3D Graphics-Vector windows upper left corner
Plane, its corresponding square frame in the show columns of right side is clicked, occurred in square frame one small to hooking.Here it is to add step e
Plane shown in 3D Graphics；
G. the square frame on the left of Translucent Plane columns is clicked on right side, is occurred in square frame one small to hooking.Here
It is that the plane that step f is shown is shown as transparent in 3D scenes.Transparency can be carried out by the Translucency on right side
Set；
H. OK is clicked on.Satellite body system YZ planes can be now shown in 3D Graphics.
5) eastern section that number passes the region that antenna may be blocked is drawn in STK
Utilize the three-dimensional of the input parameter of step 1), the result of calculation of step 2), the viewing area of step 3) and step 4)
Feature, number can be drawn on STK and passes the region that antenna may be blocked.Region herein, refers in particular to that " satellite and earth station can
In viewed area, number passes the region that antenna may be blocked ".
Operation content in this step is that drop rail Duan Shichuan or rail lift Duan Shichuan is relevant with satellite orbit.Due to current state
The track of interior most Optical remote satellite is all designed according to drop rail Duan Shichuan, therefore this explanation is also to drop rail
Illustrated exemplified by section.The method of rail lift section and drop rail section are basically identical, it is only necessary to the sign of lateral swinging angle are adjusted, due to using very
It is few, repeated no more in this explanation.
Because the region that number biography antenna may be blocked is divided into eastern section and west section, therefore described respectively.It is east first
Section.Its detailed process is as follows：
51) the maximum positive lateral swinging angle for setting the attitude of satellite not to be blocked for number biography antenna in STK
The maximum positive lateral swinging angle that set the attitude of satellite not to be blocked for number biography antenna in this step, its detailed process
It is as follows：
A. in STK, the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and sets
Window；
B. Basic-Attitude is clicked in attribute sets window, in the interface, is selected on first Type column
Fixed in Axes, on Reference Axes columns, selection " satellite name VVLH ", YPR is selected on second Type column
Angles, inserted on Roll columns " number passes the maximum lateral swinging angle that will not be blocked of antennas " counted in step 2) on the occasion of；
C. OK is clicked on, now it can be seen that the deflection situation of the attitude of satellite in 3D Graphics.
52) show that number passes antenna and may be blocked the eastern section in region
Detailed process is as follows：
A. in STK, the Start buttons on click tools column, satellite is made to orbit.In 2D Graphics
Observation, until satellite is flown in earth station viewing area.Decrease Time Step buttons on click tools column make its fortune
Scanning frequency degree slows down.Observed in 3D Graphics, until " satellite-earth station's vector " in step 41)-h is located at step 42)
In satellite body system YZ planes in；
B. the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set.
Attribute sets in window and clicks on Basic-Orbit.Constantly adjustment RAAN parameters, the Apply in interface is clicked on after having changed every time
Button, modification is set to come into force；While modification in 3D Graphics observation of steps 41) in-j " number pass antenna zero vectors with
The angle angle value of satellite-earth station's vector ", until its be equal to obtained in step 1) " number passes antenna not under by circumstance of occlusion
Hard-over ".It should be noted that the satellite orbit after modification RAAN must assure that satellite is in the east of earth station；
C. previous step 52 is constantly repeated)-a and step 52)-b, the position in orbit to satellite and RAAN parameters are carried out
Fine setting, while 3D Grahpics are observed, until two above-mentioned conditions are all met, i.e.,：In step 41)-h " satellite-
Earth station's vector " is located in the satellite body system YZ planes in step 42)；In step 41)-j " number pass antenna zero vectors with
The angle angle value of satellite-earth station's vector " " counts equal to what is obtained in step 1) and passes antenna not by the maximum under circumstance of occlusion turn
Angle "；
D. in 2D Graphics, it can be seen that the earth station viewing area of circle is divided into by line tracking under the star of satellite
Two pieces.That less piece of east side is exactly " number pass antennas may be blocked the eastern section in region ".
6) west section that number passes the region that antenna may be blocked is drawn in STK
It is similar with step 5) method, using the input parameter of step 1), the result of calculation of step 2), step 3) visible area
Domain and the three-dimensional feature of step 4), the west section that number passes the region that antenna may be blocked can be drawn on STK.Its specific mistake
Journey is as follows：
61) the maximum minus side pivot angle for setting the attitude of satellite not to be blocked for number biography antenna in STK
The maximum minus side pivot angle that set the attitude of satellite not to be blocked for number biography antenna in this step, its detailed process
It is as follows：
A. in STK, the satellite to be analyzed is clicked in the Object Browser columns of left side, is replicated, is pasted,
The satellite of a duplication is will appear from Object Browser columns；
B. the duplication satellite to be analyzed is double-clicked, satellite attribute is opened and window is set；
C. Basic-Attitude is clicked in attribute sets window, in the interface, is selected on first Type column
Fixed in Axes, on Reference Axes columns, selection " satellite name VVLH ", YPR is selected on second Type column
Angles, the negative value of " number passes the maximum lateral swinging angle that antenna will not be blocked " counted in step 2) is inserted on Roll columns；
D. OK is clicked on, now it can be seen that replicating the deflection situation of the attitude of satellite in 3D Graphics.
62) show that number passes antenna and may be blocked the west section in region
Detailed process is as follows：
A. in STK, the Start buttons on click tools column, satellite is made to orbit, in 2D Graphics
Observation, until it is flown in earth station viewing area.Decrease Time Step on click tools column make its speed of service
Slow down.Observed in 3D Graphics, until " satellite-earth station's vector " in step 41)-h is located at defending in step 42)
In star body series YZ planes；
B. the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set.
Attribute sets in window and clicks on Basic-Orbit.Constantly adjustment RAAN parameters, the Apply in interface is clicked on after having changed every time
Button, modification is set to come into force；While modification in 3D Graphics observation of steps 41) in-j " number pass antenna zero vectors with
The angle angle value of satellite-earth station's vector ", until its be equal to obtained in step 1) " number passes antenna not under by circumstance of occlusion
Hard-over ".It should be noted that the track after modification RAAN must assure that satellite is in the west of earth station；
C. previous step 62 is constantly repeated)-a and step 62)-b, the position in orbit to satellite and RAAN parameters are carried out
Fine setting, until two above-mentioned conditions are all met, i.e.,：" satellite-earth station's vector " in step 41)-h is located at step
42) in the satellite body system YZ planes in；" the folder of number biography antenna zero vector and satellite-earth station's vector in step 41)-j
Angle angle value " is equal to " number passes antenna not by the hard-over under circumstance of occlusion " obtained in step 1)；
D. 2D Graphics are entered, it can be seen that the earth station viewing area of circle is divided into by line tracking under the star of satellite
Three pieces.That less piece of east side is exactly " number pass antennas may be blocked the eastern section in region " obtained in step 52)-d.West
That less piece of side is exactly " number pass antennas may be blocked the west section in region ".Two pieces of area sums, it is that " number passes antenna can
Can be blocked region ".
7) calculate number when satellite side-sway passes in fact and pass antenna occlusion probability
After step 5) and step 6) terminate, have been obtained for " number pass antennas may be blocked region ", be divided into eastern section and
West section.Below with " the satellite maximum lateral swinging angle " and step 2) in the result and step 1) " number passes antennas will not be by
The maximum lateral swinging angle blocked ", calculate number when satellite side-sway passes in fact and pass antenna occlusion probability.It is comprised the following steps that：
71) reference area ratio
In 2D Graphics, observation of steps 62)-the d obtained eastern section and west of " number pass antennas may be blocked region "
Section.Visually the area in estimation eastern section and west section region accounts for the ratio of whole earth station viewing area, obtains area ratio.
Satellite only when " number pass antennas may be blocked region " obtained in step 62)-d is interior, passes in fact in side-sway
When be possible to number occurs to pass antenna situation about being blocked；Meanwhile " satellite and the earth station that satellite only calculates in step 3)
Side-sway could be carried out in viewing area " to pass in fact.Therefore area ratio is defined as the ratio between above-mentioned two region area.
Calculation formula is as follows：
Area ratio
≤ (eastern section area+west section the area)/gross area
72) angle ratio is calculated
Utilize " the maximum side that number biography antenna will not be blocked in the satellite maximum lateral swinging angle and step 2) in step 1)
Pivot angle ", obtain angle ratio.
Satellite only in " number passes the maximum lateral swinging angle that will not be blocked of antennas " during lateral swinging angle exceedes step 2),
Side-sway is possible to that the situation that number biography antenna is blocked occurs when passing in fact；Meanwhile lateral swinging angle during satellite side-sway reality biography will not surpass
" the satellite maximum lateral swinging angle " crossed in step 1).Therefore by angle ratio be defined as both differences and " satellite maximum lateral swinging angle " it
Than.
Calculation formula is as follows：
Angle ratio
≤ (maximum lateral swinging angle-maximum lateral swinging angle that will not be blocked)/maximum lateral swinging angle
The maximum lateral swinging angle that=1- will not be blocked/maximum lateral swinging angle
73) obtain and block probability
Using area ratio and angle ratio, number when satellite side-sway passes in fact can be obtained and pass antenna occlusion probability.
Calculation formula is as follows：
Block probability
=area is than * angles ratio
≤ (eastern section area+west section the area)/gross area
* (the maximum lateral swinging angle that 1- will not be blocked/maximum lateral swinging angle)
What above-mentioned formula calculated, be exactly number biography antenna occlusion probability when satellite side-sway passes in fact, its result of calculation is not one
The exact value of individual determination, and it is less than a value.That is, the result being calculated is to provide a maximum for blocking probability
Value.
74) provide number when satellite side-sway passes in fact and pass antenna occlusion probability level
At this point it is possible to the number provided when satellite side-sway passes in fact passes antenna occlusion probability level, the system as design of satellites
Index.Its index formulation is as follows：
● when satellite side-sway angle absolute value is spent no more than " XX ", number passes antenna and will not be blocked when side-sway passes in fact；
● when satellite side-sway angle absolute value is spent more than " XX ", number passes the probability that antenna is blocked and is no more than when side-sway passes in fact
“YY”。
XX and YY herein is calculated by step 2) in invention method and step 7) -3 respectively.
8) final result is exported
According to analysis before, the final result of step 8) can directly be reached by step 2), also can be by step 2)~step
7) gradually obtain, whether the Rule of judgment depending on step 2) meets.
Case study on implementation
Assuming that simulated conditions are as follows：
● satellite transit is in 600km orbit altitudes
It is as follows according to the inventive method, calculating process：
1) satellite maximum lateral swinging angle is obtained, and number passes antenna not by the hard-over under circumstance of occlusion
In the implementation case, satellite maximum lateral swinging angle sum is passed into antenna and not distinguished by the hard-over under circumstance of occlusion
It is set to：
Satellite maximum lateral swinging angle：35°
Number passes antenna not by the hard-over under circumstance of occlusion：90°
Above-mentioned parameter is the assumption value provided according to the experience in design of conventional model.In actual design, two parameter is equal
The real data for the strict demonstration and design of should learning from else's experience.
2) sum of subtended angle over the ground for calculating satellite passes the maximum lateral swinging angle that antenna will not be blocked
21) subtended angle over the ground of satellite is calculated
Earth radius takes 6378km, orbit altitude 600km.Then subtended angle over the ground can be calculated according to formula.Surplus angle takes
For 3 °.As a result it is：
Subtended angle over the ground
=asin (earth radius/(satellite orbital altitude+earth radius))-surplus angle
- 3 ° of=asin (6378km/ (6378km+600km))
≈66°-3°
=63 °
22) calculate number and pass the maximum lateral swinging angle that antenna will not be blocked
It is 90 ° to obtain " number passes antenna not by the hard-over under circumstance of occlusion " by step 1), then can be calculated according to formula
Number passes antenna：
Number passes the maximum lateral swinging angle that antenna will not be blocked
=number passes antenna not by the subtended angle over the ground of hard-over-satellite under circumstance of occlusion
=90 ° -63 °
=27 °
23) judge whether maximum lateral swinging angle is more than number and passes the maximum lateral swinging angle that antenna will not be blocked
From step 1), satellite maximum lateral swinging angle is 35 °.
From step 22), the maximum lateral swinging angle that number passes antenna and will not be blocked is 27 °.
Therefore, whether maximum lateral swinging angle is more than the maximum lateral swinging angle that number biography antenna will not be blocked, then can not directly enter
Enter step 8), step 3) should be entered.
3) viewing area of satellite and earth station is drawn in STK
Satellite is obtained with STK and the process of earth station viewing area is as follows：
31) simulating scenes are established
After opening STK softwares, the File-New on click tools column, a simulating scenes are generated.The scene of acquiescence is entitled
Scenario1；
32) satellite orbit is established
According to above-mentioned steps, satellite orbit is established, wherein Semimajor Axis insert the design load of semi-major axis：6378km
+ 600km=6978km.Satellite name is the Satellite1 that software default assigns.
33) earth station is established
According to above-mentioned steps, earth station is established, wherein Latitude, Longitude and Altitude is inserted successively
(0deg, 0deg, 0km).Earth station's name is the Facility1 that software default assigns.Generation one is positioned at 0 degree of east longitude, north latitude 0
Degree, height above sea level 0km earth station.
34) earth station antenna is established
According to above-mentioned steps, earth station antenna is established, wherein Cone Angle insert 85deg.Earth station antenna name is
The Sensor1 that software default assigns.One field range of generation is the earth station antenna of 85 degree of semi-cone angle, that is to say, that the ground
It is 5 degree to stand to the ground minimum elevations of satellite.
35) viewing area of satellite and earth station is shown
According to above-mentioned steps, satellite Satellite1 and earth station antenna Sensor1 viewing area are drawn.
4) three-dimensional features such as required vector, angle and plane are set in STK
41) angle that number passes antenna zero vector and satellite-earth station's vector is drawn in STK
The number to be drawn in this step passes the angle of antenna zero vector and satellite-earth station's vector, and its detailed process is such as
Under：
A. in STK, Satellite1 is clicked in the Object Browser columns of left side, is then up selected on toolbar
Satellite-Vector Geometry are selected, open " Vector Geometry Tool for Satellite1 " windows；
B. in an open window, right side Create New Vector buttons are clicked on.In an open window, in Name
In column LingWeiShiLiang is named to vector；In Type columns selection Fixed in Axes；Selected in Reference Axes
Select satellite Satellite1 Body (namely body coordinate system)；Cartesian is selected in Type；Filled out in X, Y, Z
(0,0,1), i.e., point to substar along satellite+Z axis.OK is clicked on, establishes vector completion.Here what is established is that number passes antenna zero arrow
Measure LingWeiShiLiang；
C. still in an open window, right side Create New Vector buttons are clicked on.It is named on Name columns to vector
DiMianZhanShiLiang；Displacement is selected on Type columns；In Origin Point selection satellites Satellite1
Center；In Destination Point selections earth station Facility1 Center.OK is clicked on, establishes vector completion.
That establish here is satellite-earth station vector DiMianZhanShiLiang；
D. still in an open window, right side Create New Angle buttons are clicked on.In an open window, exist
JiaJiao is named to angle in Name columns；In Type columns selection Between Vectors；Select to build just now in From Vector
Vertical number passes antenna zero vector LingWeiShiLiang；Satellite-earth station's the arrow established just now is selected in To Vector
Measure DiMianZhanShiLiang.OK is clicked on, establishes angle completion.That establish here is the angle JiaJiao of above-mentioned two vector；
E. Close is clicked on, closes " Vector Geometry Tool for Satellite1 " windows；
F. in STK, Satellite1 is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set
Mouthful；
G. 3D Graphics-Vector are clicked in attribute sets window, Add is clicked in the interface, what is newly opened
Satellite-earth station vector the DiMianZhanShiLiang and two vector angle JiaJiao established just now is double-clicked in window, is clicked on
OK closes window.Here it is that satellite-earth station the vector DiMianZhanShiLiang and step d for establishing step c is established
Two vector angle JiaJiao are added in 3D Graphics；
H. satellite-earth station's the arrow added just now is found in the Name columns in the 3D Graphics-Vector windows upper left corner
DiMianZhanShiLiang is measured, clicks its corresponding square frame in the show columns of right side, is occurred in square frame one small to hooking.This
In be that the satellite-earth station vector DiMianZhanShiLiang for adding step g is shown in 3D Graphics；
I. the angle JiaJiao added just now is found in the Name columns in the 3D Graphics-Vector windows upper left corner,
Its corresponding square frame in the show columns of right side is clicked, is occurred in square frame one small to hooking.Here it is the angle for adding step g
JiaJiao is shown in 3D scenes；
J. the square frame on the left of Show Angle Value columns is clicked on right side, is occurred in square frame one small to hooking.Here
It is to show the angle value for the angle JiaJiao that step i is shown in 3D scenes；
K. OK is clicked on, satellite-earth station's vector can be now shown in 3D Graphics
DiMianZhanShiLiang, and number pass antenna zero vector LingWeiShiLiang and satellite-earth station's vector
DiMianZhanShiLiang angle JiaJiao and its angle value.
42) satellite body system YZ planes are set in STK
The satellite body system YZ planes to be drawn in this step, its detailed process are as follows：
A. in STK, the satellite to be analyzed is clicked in the Object Browser columns of left side, then up toolbar
Upper selection Satellite-Vector Geometry, open " Vector Geometry Tool for Satellite1 " windows
Mouthful；
B. in an open window, right side Create New Plane buttons are clicked on.In an open window, on Name columns
YZ-PingMian is named to plane；Quadrant is selected on Type columns；YZ is selected in Quadrant columns；In Reference
Satellite Satellite1 Body (namely body coordinate system) is selected in System.OK is clicked on, establishes plane completion.Here build
That vertical is satellite body system YZ planes YZ-PingMian；
C. Close is clicked on, closes " Vector Geometry Tool for Satellite1 " windows；
D. in STK, the satellite to be analyzed is double-clicked in the Object Browser columns of left side, satellite attribute is opened and sets
Window；
E. 3D Graphics-Vector are clicked in attribute sets window, Add is clicked in the interface, what is newly opened
The satellite body system YZ plane YZ-PingMian established just now are double-clicked in window, OK is clicked on and closes window.Here it is by step b
The satellite body system YZ planes YZ-PingMian of foundation is added in 3D Graphics；
F. the satellite body system YZ added just now is found in the Name columns in the 3D Graphics-Vector windows upper left corner
Plane YZ-PingMian, its corresponding square frame in the show columns of right side is clicked, occurred in square frame one small to hooking.Here it is
The plane that step e is added is shown in 3D Graphics；
G. the square frame on the left of Translucent Plane columns is clicked on right side, is occurred in square frame one small to hooking.Here
It is that the plane that step f is shown is shown as transparent in 3D scenes.Transparency is arranged to by the Translucency on right side
80%；
H. OK is clicked on.Satellite body system YZ planes YZ-PingMian can be now shown in 3D Graphics.
5) eastern section that number passes the region that antenna may be blocked is drawn in STK
Utilize the three-dimensional of the input parameter of step 1), the result of calculation of step 2), the viewing area of step 3) and step 4)
Feature, number can be drawn on STK and passes the region that antenna may be blocked.Region herein, refers in particular to that " satellite and earth station can
In viewed area, number passes the region that antenna may be blocked ".
It is eastern section first.Its detailed process is as follows：
51) the maximum positive lateral swinging angle for setting the attitude of satellite not to be blocked for number biography antenna in STK
The maximum positive lateral swinging angle that set the attitude of satellite not to be blocked for number biography antenna in this step, its detailed process
It is as follows：
A. in STK, Satellite1 is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set
Mouthful；
B. Basic-Attitude is clicked in attribute sets window, in the interface, is selected on first Type column
Fixed in Axes, on Reference Axes columns, selection " satellite name VVLH ", YPR is selected on second Type column
Angles, inserted on Roll columns counted number in step 2) pass the maximum lateral swinging angles that will not be blocked of antennas on the occasion of 27deg；
C. OK is clicked on, now it can be seen that the deflection situation of the attitude of satellite in 3D Graphics.
52) show that number passes antenna and may be blocked the eastern section in region
Detailed process is as follows：
A. in STK, the Start buttons on click tools column, satellite Satellite1 is made to orbit.In 2D
Observed in Graphics, until satellite Satellite1 is flown in earth station Facility1 viewing areas.On click tools column
Decrease Time Step slow down its speed of service.Observed in 3D Graphics, until satellite-earth station's vector
DiMianZhanShiLiang is located in satellite body system YZ planes YZ-PingMian.
B. Satellite1 is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set.In attribute
Set in window and click on Basic-Orbit.Constantly adjustment RAAN parameters, the Apply buttons in interface are clicked on after having changed every time,
Modification is set to come into force；Observed number passes antenna zero vector LingWeiShiLiang with defending in 3D Graphics while modification
Star-earth station vector DimianZhanShiLiang angle JiaJiao angle value, until it is equal to what is obtained in step 1)
" number passes antenna not by the hard-over under circumstance of occlusion " --- 90 °.It should be noted that the satellite orbit after modification RAAN must
It must ensure that satellite Satellite1 is in earth station Facility1 east；
C. previous step 52 is repeated)-a and step 52)-b, the position in orbit to satellite and RAAN parameters carry out micro-
Adjust, while observe 3D Grahpics, until two above-mentioned conditions are all met, i.e.,：Satellite-earth station's vector
DiMianZhanShiLiang is located in satellite body system YZ planes YZ-PingMian；Number passes antenna zero vector
LingWeiShiLiang and satellite-earth station vector DimianZhanShiLiang angle JiaJiao angle values are equal to step
1) " number passes antenna not by the hard-over under circumstance of occlusion " obtained in --- 90 °.
D. in 2D Graphics, it can be seen that the earth station viewing area of circle is divided into by line tracking under the star of satellite
Two pieces.That less piece of east side is exactly " number pass antennas may be blocked the eastern section in region ".
6) west section that number passes the region that antenna may be blocked is drawn in STK
It is similar with step 5) method, using the input parameter of step 1), the result of calculation of step 2), step 3) visible area
Domain and the three-dimensional feature of step 4), the west section that number passes the region that antenna may be blocked can be drawn on STK.Its specific mistake
Journey is as follows：
61) the maximum minus side pivot angle for setting the attitude of satellite not to be blocked for number biography antenna in STK
The maximum minus side pivot angle that set the attitude of satellite not to be blocked for number biography antenna in this step, its detailed process
It is as follows：
A. in STK, Satellite1 is clicked in the Object Browser columns of left side, is replicated, is pasted, in Object
The satellite Satellite2 of a duplication is will appear from Browser columns；
B. Satellite2 is double-clicked, satellite attribute is opened and window is set；
C. Basic-Attitude is clicked in attribute sets window, in the interface, is selected on first Type column
Fixed in Axes, " Satellite2VVLH " is selected on Reference Axes columns, and YPR is selected on second Type column
Angles, inserted on Roll columns the negative value of " number passes the maximum lateral swinging angle that will not be blocked of antennas " counted in step 2)-
27deg；
D. OK is clicked on, now it can be seen that the deflection situation of Satellite2 postures in 3D Graphics.
62) show that number passes antenna and may be blocked the west section in region
Detailed process is as follows：
A. in STK, the Start buttons on click tools column, satellite Satellite2 is made to orbit.In 2D
Observed in Graphics, until satellite Satellite2 is flown in earth station Facility1 viewing areas.On click tools column
Decrease Time Step slow down its speed of service.Observed in 3D Graphics, until satellite-earth station's vector
DiMianZhanShiLiang is located in satellite body system YZ planes YZ-PingMian.
B. Satellite2 is double-clicked in the Object Browser columns of left side, satellite attribute is opened and window is set.In attribute
Set in window and click on Basic-Orbit.Constantly adjustment RAAN parameters, the Apply buttons in interface are clicked on after having changed every time,
Modification is set to come into force；Observed number passes antenna zero vector LingWeiShiLiang with defending in 3D Graphics while modification
Star-earth station vector DiMianZhanShiLiang angle JiaJiao angle values, until it is equal to what is obtained in step 1)
" number passes antenna not by the hard-over under circumstance of occlusion " --- 90 °.It should be noted that the track after modification RAAN must be protected
Card satellite is in the west of earth station；
C. previous step 62 is repeated)-a and step 62)-b, the position in orbit to satellite and RAAN parameters carry out micro-
Adjust, until two above-mentioned conditions are all met, i.e.,：Satellite-earth station vector DiMianZhanShiLiang is located at satellite
In body series YZ planes YZ-PingMian；Number passes antenna zero vector LingWeiShiLiang and satellite-earth station's vector
DiMianZhanShiLiang angle JiaJiao angle values are equal to " counting and pass antenna not by circumstance of occlusion for the middle acquisition of step 1)
Under hard-over " --- 90 °.
D. 2D Graphics are entered, it can be seen that the earth station viewing area of circle is divided into by line tracking under the star of satellite
Three pieces.That less piece of east side is exactly " number pass antennas may be blocked the eastern section in region " obtained in step 52)-d.West
That less piece of side is exactly " number pass antennas may be blocked the west section in region ".Two pieces of area sums, it is that " number passes antenna can
Can be blocked region ".
7) calculate number when satellite side-sway passes in fact and pass antenna occlusion probability
After step 5) and step 6) terminate, have been obtained for " number pass antennas may be blocked region ", be divided into eastern section and
West section.Below with " the satellite maximum lateral swinging angle " and step 2) in the result and step 1) " number passes antennas will not be by
The maximum lateral swinging angle blocked ", calculate number when satellite side-sway passes in fact and pass antenna occlusion probability.It is comprised the following steps that：
71) reference area ratio
In 2D Graphics, observation of steps 62)-the d obtained eastern section and west of " number pass antennas may be blocked region "
Section.Visually the area in estimation eastern section and west section region accounts for the ratio of whole earth station viewing area, obtains area ratio.It is estimated,
Shaded area accounts for the 10% of whole earth station viewing area area
Calculation formula is as follows：
Area ratio
≤ (eastern section area+west section the area)/gross area
≈ 10% (visually estimation)
72) angle ratio is calculated
Utilize " the maximum side that number biography antenna will not be blocked in the satellite maximum lateral swinging angle and step 2) in step 1)
Pivot angle ", obtain angle ratio.
Calculation formula is as follows：
Angle ratio
≤ (maximum lateral swinging angle-maximum lateral swinging angle that will not be blocked)/maximum lateral swinging angle
The maximum lateral swinging angle that=1- will not be blocked/maximum lateral swinging angle
=1-27 °/35 °
≈ 22.86%
73) obtain and block probability
Using area ratio and angle ratio, number when satellite side-sway passes in fact can be obtained and pass antenna occlusion probability.
Calculation formula is as follows：
Block probability
=area is than * angles ratio
≤ 10%*22.86%
≈ 3%
What above-mentioned formula calculated, be exactly number biography antenna occlusion probability when satellite side-sway passes in fact, its result of calculation is not one
The exact value of individual determination, and it is less than a value.That is, the result being calculated is to provide a maximum for blocking probability
Value.
74) provide number when satellite side-sway passes in fact and pass antenna occlusion probability level
At this point it is possible to the number provided when satellite side-sway passes in fact passes antenna occlusion probability level, the system as design of satellites
Index.Its index formulation is as follows：
● when satellite side-sway angle absolute value is no more than 27 °, number passes antenna and will not be blocked when side-sway passes in fact；
● when satellite side-sway angle absolute value is spent more than 27 °, number passes the probability that antenna is blocked and is no more than when side-sway passes in fact
3%.
8) final result is exported
According to analysis before, the final result of step 8) can directly be reached by step 2), also can be by step 2)~step
7) gradually obtain, whether the Rule of judgment depending on step 2) meets.In the implementation case, because the condition of step 8) obtains
Meet, therefore final result is gradually obtained by step 2)~step 7).
The content not being described in detail in description of the invention, especially with respect to STK operating method, belong to this area skill
The known technology of art personnel.
Claims (3)
1. a kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability, it is characterised in that step is as follows：
1) satellite maximum lateral swinging angle is obtained, and number passes antenna not by the hard-over under circumstance of occlusion；
2) calculate and obtain the maximum lateral swinging angle that the sum of the subtended angle over the ground biography antenna of satellite will not be blocked, judge satellite maximum side
Whether pivot angle, which is more than number, passes the maximum lateral swinging angle θ that antenna will not be blocked, if the satellite maximum side-sway obtained in step 1)
Angle is more than counted number in step 2) and passes the maximum lateral swinging angle θ that antenna will not be blocked, then into step 3)；Otherwise, show
Number passes antenna and will not be blocked certainly, and it is 0% to block probability, and it is general that number when now obtaining satellite side-sway reality biography passes antenna occlusion
Number passes antenna and will not be blocked when rate index, i.e. side-sway pass in fact；
3) viewing area of satellite and earth station is drawn in STK；
4) three-dimensional feature needed for being set in STK
41) angle that number passes antenna zero vector and satellite-earth station's vector is drawn in STK；
42) satellite body system YZ planes are set in STK；
5) region that number biography antenna may be blocked is divided into two pieces, respectively positioned at Satellite offline east side and west side, difference
It is defined as eastern section and west section that number passes the region that antenna may be blocked；The lateral swinging angle of satellite is set to pass day for number in STK
The maximum positive lateral swinging angle that line will not be blocked, draw the eastern section that number passes the region that antenna may be blocked；
6) the maximum minus side pivot angle that the setting attitude of satellite will not be blocked for number biography antenna in STK, drawing number biography antenna can
The west section in the region that can be blocked；
7) calculate number when satellite side-sway passes in fact and pass antenna occlusion probability
71) calculate and obtain area ratio
Determine number that step 6) obtains pass antennas may be blocked region eastern section and west section area, calculate eastern section and west section region
Area account for the ratio of whole earth station viewing area, obtain area than calculation formula it is as follows：
Area ratio≤(eastern section area+west section area)/ground gross area
72) calculate and obtain angle ratio
The maximum lateral swinging angle that will not be blocked using the number biography antenna in the satellite maximum lateral swinging angle and step 2) in step 1),
Calculate and obtain angle ratio≤(satellite maximum lateral swinging angle-number passes the maximum lateral swinging angle that antenna will not be blocked)/satellite maximum side
Pivot angle=1- numbers pass maximum lateral swinging angle/satellite maximum lateral swinging angle that antenna will not be blocked；
73) obtain and block Probability p=area than * angles ratio≤(eastern section area+west section area)/ground gross area * (1- numbers biography days
The maximum lateral swinging angle that line will not be blocked/satellite maximum lateral swinging angle)；
74) now obtain number when satellite side-sway passes in fact and pass antenna occlusion probability level, i.e., absolute value does not surpass when satellite side-sway angle
When crossing number biography antenna not by the maximum lateral swinging angle θ under circumstance of occlusion, number passes antenna and will not be blocked when side-sway passes in fact, works as satellite
Lateral swinging angle absolute value exceed number pass antenna not by the maximum lateral swinging angle θ under circumstance of occlusion when, side-sway in fact pass when number pass antenna by
The probability blocked is no more than p.
2. the method that a kind of number obtained when the in-orbit progress of satellite is real to be passed according to claim 1 passes antenna occlusion probability,
It is characterized in that：The sum of subtended angle over the ground of step 2) Satellite passes maximum lateral swinging angle that antenna will not the be blocked specifically side of calculating
Method is：
Subtended angle=asin (earth radius/(satellite orbital altitude+earth radius))+surplus angle over the ground；Number passes antennas will not be by
Antenna is passed not by the subtended angle over the ground of hard-over-satellite under circumstance of occlusion to the maximum lateral swinging angle blocked=number.
3. the method that a kind of number obtained when the in-orbit progress of satellite is real to be passed according to claim 1 passes antenna occlusion probability,
It is characterized in that：The specific method for drawing the viewing area of satellite and earth station in step 3) in STK is as follows：
31) simulating scenes are established；
32) satellite orbit is established；
33) earth station is established；
34) earth station antenna is established；
35) viewing area of satellite and earth station is shown.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510527488.8A CN105183959B (en) | 2015-08-25 | 2015-08-25 | A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510527488.8A CN105183959B (en) | 2015-08-25 | 2015-08-25 | A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105183959A CN105183959A (en) | 2015-12-23 |
| CN105183959B true CN105183959B (en) | 2018-03-09 |
Family
ID=54906038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510527488.8A Active CN105183959B (en) | 2015-08-25 | 2015-08-25 | A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105183959B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106095731B (en) * | 2016-08-18 | 2018-04-20 | 中国电子科技集团公司第五十四研究所 | A kind of satellite health blocks the computational methods that number over the ground passes the antenna time |
| CN106647704B (en) * | 2016-12-19 | 2019-11-12 | 北京空间飞行器总体设计部 | A kind of GEO track number passes antenna tracking accuracy and determines method |
| FR3069523A1 (en) * | 2017-07-27 | 2019-02-01 | Prodose | METHOD OF MAKING A NETWORK FOR THE SUPPLY OF THE INTERNET IN PARTICULAR ON THE SURFACE OF THE GLOBE TERRESTRIAN, A PLANE FOR CARRYING OUT IT |
| CN113691301B (en) * | 2021-08-23 | 2022-08-02 | 深圳海卫通网络科技有限公司 | Method and device for measuring shielding of ship satellite communication antenna |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0618541A1 (en) * | 1993-04-01 | 1994-10-05 | International Business Machines Corporation | Interactive graphics computer system for planning star-sensor-based satellite attitude maneuvers |
| CN102700728A (en) * | 2012-06-18 | 2012-10-03 | 上海卫星工程研究所 | Method for determining sheltering of solar cell array during in-orbit running of satellite and application of method |
| CN103970991A (en) * | 2014-04-29 | 2014-08-06 | 上海卫星工程研究所 | Method for analyzing on-orbit sheltered condition of solar battery arrays of stationary orbiting satellite |
| CN104410361A (en) * | 2014-10-10 | 2015-03-11 | 中国空间技术研究院 | A satellite solar wing occlusion testing method |
-
2015
- 2015-08-25 CN CN201510527488.8A patent/CN105183959B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0618541A1 (en) * | 1993-04-01 | 1994-10-05 | International Business Machines Corporation | Interactive graphics computer system for planning star-sensor-based satellite attitude maneuvers |
| CN102700728A (en) * | 2012-06-18 | 2012-10-03 | 上海卫星工程研究所 | Method for determining sheltering of solar cell array during in-orbit running of satellite and application of method |
| CN103970991A (en) * | 2014-04-29 | 2014-08-06 | 上海卫星工程研究所 | Method for analyzing on-orbit sheltered condition of solar battery arrays of stationary orbiting satellite |
| CN104410361A (en) * | 2014-10-10 | 2015-03-11 | 中国空间技术研究院 | A satellite solar wing occlusion testing method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105183959A (en) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105183959B (en) | A kind of method that number obtained when the in-orbit progress of satellite is real to be passed passes antenna occlusion probability | |
| CN102607526B (en) | Target posture measuring method based on binocular vision under double mediums | |
| CN109592079A (en) | A kind of spacecraft coplanar encounter of limiting time becomes rail strategy and determines method | |
| CN104573251A (en) | Method for determining full-field-of-view apparent spectral radiance of satellite-borne optical remote sensor | |
| CN104154919A (en) | Method for autonomous measurement of pose of tripod structure of solar panel on non-cooperative spacecraft | |
| CN105628055B (en) | A kind of deep space probe landing autonomous optical navigation target imaging simulation system | |
| CN106352566A (en) | Tower type solar thermal power generation sun-tracking analysis system | |
| CN109612438B (en) | Method for determining initial orbit of space target under constraint of virtual coplanar condition | |
| CN108680165A (en) | Target aircraft attitude determination method based on optical imagery and device | |
| CN110146093A (en) | Binary asteroid detection independently cooperates with optical navigation method | |
| CN103150414A (en) | STL (standard template library)-based analysis method on sheltering to solar wing by communications satellite large antenna | |
| Souza et al. | A 3D-gis extensionf for sky view factors assessment in urban environment | |
| CN103679780A (en) | A space target real time simulation method | |
| CN105844700A (en) | System for acquiring three-dimensional point clouds in outdoor scene | |
| Basmadji et al. | Development of ground station for a terrain observer-hardware in the loop simulations | |
| CN110220492A (en) | A kind of sub-satellite point solar elevation simulating analysis and system | |
| Roganov et al. | The task of increasing the number of training situations as one of the areas for improving the ergatic software and hardware systems “Aviation simulator” | |
| CN106095731B (en) | A kind of satellite health blocks the computational methods that number over the ground passes the antenna time | |
| CN106372036B (en) | A kind of satellite sun wing blocks the computational methods of number biography antenna time over the ground | |
| CN103968835A (en) | Simulating method of refraction star | |
| Liang et al. | Accuracy analysis of oblique photogrammetry measurement in 3D modeling of power line selection design | |
| CN114485620A (en) | Orbital dynamics fused asteroid detector autonomous visual positioning system and method | |
| CN103544677A (en) | Space-air-ground integration situational expression engine and shaking elimination method | |
| Cuzzocrea | Analysis and validation of spaceborne synthetic imagery using a vision-based pose initialization algorithm for non-cooperative spacecrafts | |
| CN108897029A (en) | Noncooperative target short distance Relative Navigation vision measurement system index evaluating 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |