1、 Rec. ITU-R S.1559 1 RECOMMENDATION ITU-R S.1559 Methodology for computing the geographical distribution of maximum downlink equivalent power flux-density levels generated by non- geostationary fixed-satellite service systems using circular orbits (Question ITU-R 236/4) (2002) The ITU Radiocommunica
2、tion Assembly, considering a) that interference from non-geostationary (GSO) systems into GSO networks will likely have significant variations both geographically and temporally; b) that geographic distributions of interference levels caused by non-GSO systems vary with traffic loading, cell layouts
3、, and constellation configurations; c) that interference from non-GSO systems into GSO networks can vary with beam scheduling algorithms of the non-GSO system; d) that non-GSO beam scheduling algorithms may likely be adjusted to meet traffic loading variations and other factors; e) that non-GSO inte
4、rference levels are dependent on the geographic location of the GSO ground station and GSO satellite location; f) that information on the geographical distribution of maximum downlink equivalent power flux-density (epfd) levels generated by non-GSO systems may be useful to GSO systems designers in o
5、rder to determine the expected level of non-GSO interference, recommends 1 that the method described in Annex 1 could be used to construct geographic distributions of maximum epfdlevels generated within a GSO satellite coverage area by a non-GSO system that uses circular orbits; 2 that Annex 1 could
6、 be used by administrations, including GSO/fixed-satellite service (FSS) network operators, to provide them with guidance for designing links for a specific geographical location given the following Notes: NOTE 1 Validation software has been specified to check compliance of a non-GSO system with epf
7、dlimits in Article 22 of the Radio Regulations (see Recommendation ITU-R S.1503). This software generates epfd levels which are meant to be an absolute envelope and are not representative of what would be generated over time by the non-GSO FSS system in operation. NOTE 2 If the methodology in Annex
8、1 is applied using the assumptions and satellite pfd mask approach as defined in the Recommendation ITU-R S.1503 validation approach, geographical distribution will be obtained of the absolute maximum epfdenvelope. In the case where the method in Annex 1 is applied to simulation models more represen
9、tative of the non-GSO operation, the maps generated are likely to vary during the lifetime of the non-GSO system as operating parameters will change. 2 Rec. ITU-R S.1559 ANNEX 1 Algorithm for calculating geographic distribution of maximum epfdinterference levels caused by non-GSO systems 1 Introduct
10、ion Since non-GSO interference will vary geographically it could be useful for administrations to be able to quantify this interference using software that can estimate the geographic and temporal distribution of non-GSO FSS interference. The analysis will represent the maximum interference expected
11、 over the life of the non-GSO system. Analyses of typical scenarios provide little value since the interference levels can have extremely large variations. The methodology presented in this Annex could be used to generate maps showing representative maximum non-GSO interference power levels that cou
12、ld be received on any land area. The maps will provide guidance for GSO operators for designing links. These maps would provide GSO operators with additional knowledge of maximum epfd locations. This Annex gives a method for calculating the geographic distribution of maximum epfd levels for non-GSO
13、systems with circular orbits that use an exclusion angle for GSO interference mitigation. For these non-GSO systems the maximum epfdinterference levels (for GSO ground stations with 3 m or greater antennas) occur on or near where the non-GSO spacecraft is inline between the GSO ground station and th
14、e GSO spacecraft. For non-GSO systems where this is not the case this approach does not apply. The maximum epfdis found from a computer simulation of the non-GSO system. However, since the value is near an inline situation, only a small number of time points need to be simulated to determine the max
15、imum epfd. Only points that are within the 10 dB beamwidth of the GSO ground station need to be considered since in these types of non-GSO systems locations outside of the 10 dB beamwidth will almost certainly be at least 10 dB below the maximum epfdvalue found on the Earths surface. Since most non-
16、GSO systems have unique characteristics, the proposed method does not describe how to simulate the specific non-GSO system. Instead it computes the time periods needed to be simulated. A resolution of 1 longitude by 1 latitude is sufficient resolution to show the general variations in interference l
17、evels, however finer resolution may be necessary for more detailed studies, especially in the case of larger antennas. It is noted that non-GSO systems with continuously repeating ground tracks will have maximum epfdvariations over a much smaller area. 2 Outline of algorithm The method given can be
18、used to calculate time periods during a non-GSO simulation when a spacecraft passes though a given beamwidth of a GSO ground antenna. First a piece-wise linear approximation to the intersection of the GSO ground antenna beam with the non-GSO sphere is Rec. ITU-R S.1559 3 computed. For each point on
19、the intersection, exact orbital parameters of a non-GSO spacecraft are computed such that the spacecraft will intersect the point on the next orbit. By comparing a given spacecrafts orbital parameters with the orbital parameters at each point, the exact times when the satellite passes through the an
20、tenna beam can be determined. 3 Symbols used ER: Earth rotational coordinates. An orthogonal three-dimensional coordinate system centred on Earths centre. ERs z-axis passes through the North Pole and the x-axis passes through the longitude 0 latitude 0 point GT: three dimensional location of GSO gro
21、und station in ER coordinates G: three dimensional location of GSO satellite in ER coordinates Gx: x component of vector G Gy: y component of vector G Gz: z component of vector G : normalizes a vector (i.e. V is equivalent to V / |V| ) : dot product operator : cross product operator Min(Qi) determin
22、es the minimum of Qifor all i Max(Qi) Determines the maximum of Qifor all i Ax: x axis of coordinate system pointing from GSO ground station to the GSO Ay:y axis of coordinate system pointing from GSO ground station to the GSO Az:z axis of coordinate system pointing from GSO ground station to the GS
23、O : half of the 3 dB beamwidth of GSO ground station antenna Vi: vectors forming a cone with half angle and centre pointing in z direction Wi: vectors forming a cone with half angle and centre pointing in GSO direction Qi:locations on the non-GSO sphere that intersect vectors Wi0: mean anomaly = 0 0
24、: perigee argument at an initial moment r: perigee argument precession rate : initial longitude r: ascending node longitude precession rate e:Earths rotation rate 4 Rec. ITU-R S.1559 k:k-th longitude crossing for a non-GSO spacecraft S: semi-major axis Tia: number of times during an ascending orbit
25、when the test spacecraft passes through QiTid: number of times during a descending orbit when the test spacecraft passes through Qiia: initial longitude of orbits that cross point Qiwhile ascending in latitude id: initial longitude of orbits that cross point Qiwhile descending in latitude tk: number
26、 of times when a non-GSO satellite ascends across the equatorial plane n: number of vectors in the cone. 4 Calculation of intersection of GSO ground terminal antenna beam with non-GSO sphere Vishown in Fig. 1, defines a set of vectors that sweep out a conic shape centred at the origin with half of t
27、he 3 dB beamwidth of . The centre of the cone points toward the ER coordinate systems z axis. Vi= (cos(i) sin(), sin(i) sin(), cos() (1) where: i = 0, ., n i= 2i/(n) = angle between any two adjacent vectors in Fig. 1 (2) : half of the 3 dB beamwidth. 1559-01zyxV0V5V1V2V4V3FIGURE 1Visual representati
28、on of vector Vi1/2 angle beamwidthRec. ITU-R S.1559 5 Next an orthogonal coordinate system Ax, Ay, Azis computed such that the Azaxis points from the GSO ground terminal toward the GSO. Az= (G GT) (3) T = (0, 0, 1) (4) Ax= (T Az) (5) Ay= Az Ax(6) The cone is pointed toward the GSO yielding the vecto
29、rs Wishown in Fig. 2. This is done using a simple coordinate transform. Wi= Ax(Vix) + Ay(Viy) + Az(Viz) (7) 1559-02GzyxW5W1W4W3W2W0AzAyAxGSOGSO ground stationGTFIGURE 2Visual representation of vector WiThe vectors Wiare extended from GT toward G, using the methods defined in equations (8) to (14), u
30、ntil they intersect the non-GSO sphere at points Qias shown in Fig. 3. Liin equation (8), represents a line beginning at GT and extending in the direction of Wi. Liwill intersect the non-GSO sphere when | Li| is equal to the non-GSO semi-major axis. Equations (10) to (13) solve the quadratic equatio
31、n shown in equation (9) to find the intersection point Qi. 6 Rec. ITU-R S.1559 Li= Wi+ GT (8) Li Li= S2(9) ai= Wi Wi(10) bi= 2 Wi GT (11) ci= GT GT S2(12) i= (bi+ (bi2 2 aici)0.5) / (2 ai) (13) Qi= iWi+ G (14) 1559-03Q0Q5Q1Q4Q2Q3Sphere with radius equal to thenon-GSO orbit radiusGSOGSO ground statio
32、nGTEarths centreFIGURE 3Intersection of vector Wiwith the non-GSO sphereRec. ITU-R S.1559 7 5 Calculation of non-GSO orbital parameters and times when satellite and beam intersection points (Qi) are co-located A test spacecraft with the orbital parameters of a non-GSO satellite and a mean anomaly of
33、 zero is chosen such that it will intersect Qion the next orbit. Equation (15) gives an expression for the orbit location of the test spacecraft in ER coordinates. =)(sin)(sin)cos()cos()(sin)cos()(cos()cos()(sin)(sin)cos()(cos(iSiSiSzyx(15) where: = 0+ 0+ (r+ 2/T ) t = + (r e) t 0: mean anomaly = 0
34、0: perigee argument at an initial moment r: perigee argument precession rate : initial longitude r: ascending node longitude precession rate e: Earths rotation rate S : semi-major axis i : inclination angle. For a Qipoint that is lower in latitude then the inclination angle, there will be two unique
35、 orbits of the test spacecraft that intersect Qi. One occurs when the spacecraft is ascending in latitude and the other when the spacecraft is descending. For each of these orbits, both the initial longitude (time = 0) and intersection times with Qiare computed. Since the z coordinate of the test sp
36、acecraft is S sin() sin(i), can have two possible solutions. One when the spacecraft is ascending in latitude and the other when the spacecraft is descending in latitude. i1= sin1(iz/(S sin(i) (16) i2= i1Tij, the number of times during the orbit when the test spacecraft passes through Qiis given by:
37、 Tij= (ij 0) / (r+ 2/P) for j = 1, 2 (17) where P is orbit period. From Tijthe initial longitude is computed by solving the simultaneous equations for x and y in equation (15). ij= cos1(Qix+ Qiy) / (R(cos2(ij) + sin2(ij) rTi eTi(18) where R is orbit radius. 8 Rec. ITU-R S.1559 6 Calculating longitud
38、e crossings for a given simulated satellite The times at which a satellite crosses the equatorial plane is calculated as follows: tk= (2k 0 0) / (r+ 2/P) (19) where: k = 1, , M M : number of orbits before satellite repeats 0: mean anomaly = 0 0: perigee argument at an initial moment r: perigee argum
39、ent precession rate P : orbit period. The corresponding longitude ascending crossings are given by: k= + (r e)tk(20) 7 Determining if an orbit crosses through the beam of the GSO ground station antenna If a satellites ascending longitude crossing is in-between the min(ij) and max(ij) for any j, then
40、 the satellite will cross the beam of the GSO antenna at some point during the next orbit. It is important to note that the modulo nature of longitude is taken into account when making this calculation. 8 Determining the simulation time of a spacecraft crossing through the beam of the GSO ground sta
41、tion antenna The simulated time of an orbit crossing an antenna beam is calculated by intersecting the orbit with the piecewise line segments of Qias shown in Fig. 4. If an ascending longitude crossing is in-between the min(ij) and max(ij) then the orbit will cross the beam of the GSO antenna. The o
42、rbit will intersect the Qiat two locations. By interpolating the corresponding Tijthe intersection times are approximated. Let: Ni= (i + 1) mod n (21) For any non-zero duration crossing there will be one value of i = i1 such that ij k Nij.Also there will be one value of i = i2 such that ij k Nij. NO
43、TE 1 Points near the equator might have only one intersection since there may be a crossing that begins on one orbit and ends on the next. The corresponding intersection time after the equatorial crossing can be approximated from a linear interpolation of the known intersection times Ti. Rec. ITU-R
44、S.1559 9 Let: 1= (k i1j) / (Ni1j i1j) (22) int_t1= 1TNi1j+ (1 1)Ti1j(23) 2= (k i2j) / (Ni2j i2j) (24) int_t2= 2TNi2j+ (1 2)Ti2j(25) Since int_t1and int_t2represent the times since a longitude crossing the time in the simulation is computed as follows: Start_timek= min(int_t1,int_ t2) + tk(26) End_ti
45、mek= max(int_ t1,int_ t2) + tk(27) for all k such that min(ij) k max(ij) (28) These Start and End times are then computed for all of the non-GSO spacecraft. 1559-04Sphere with radius equal to thenon-GSO orbit radiuststarttendInitiallongitudenon-GSO orbitEquatorialplaneFIGURE 4Intersection of non-GSO
46、 orbit with antenna beam10 Rec. ITU-R S.1559 9 Example computation Non-GSO systems will naturally have a variation in maximum epfdlevels based on latitude, distance from the nearest non-GSO gateway and elevation angle from the GSO ground station to the supporting GSO spacecraft. An example of such a
47、 map is shown in Fig. 5. This map assumes the use of hypothetical pfd masks according to Recommendation ITU-R S.1503 (a fully loaded system, an envelope of all possible scheduling algorithms), and non-repeating ground tracks. It should be noted that the methodology is also valid for repeating ground
48、-tracks and other operating assumptions. These maps may be used for sensitive GSO link planning purposes, as they represent an upper bound on the interference distribution, even though the maps generated by the system in operation may well show generally lower levels of maximum epfdfor any geographi
49、cal location. 1559-05epfddB(W/(m2 40 kHz)160 to 162162 to 164164 to 166166 to 168168 to 170FIGURE 5Example of geographic distribution of maximum epfd levels10 Future work This Recommendation only treats the case of non-GSO systems that use circular orbits. It would be useful to further develop the methodology so that it could be used for the c
